Next Article in Journal
Study on the Algae Lysis Method of White Rot Fungi Algae Control System
Next Article in Special Issue
Advances in Transboundary Aquifer Assessment
Previous Article in Journal
Maximizing Sustainability in Reservoir Operation under Climate Change Using a Novel Adaptive Accelerated Gravitational Search Algorithm
Previous Article in Special Issue
Investigation of the Origin of Hueco Bolson and Mesilla Basin Aquifers (US and Mexico) with Isotopic Data Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydrogeomorphologic Mapping of the Transboundary San Pedro Aquifer: A Tool for Groundwater Characterization

by
José Ismael Minjárez Sosa
,
Grisel Alejandra Gutiérrez Anguamea
,
Rogelio Monreal
,
Francisco Javier Grijalva Noriega
and
Elia M. Tapia-Villaseñor
*
Departamento de Geología, Universidad de Sonora, Hermosillo 83000, Mexico
*
Author to whom correspondence should be addressed.
Water 2022, 14(6), 906; https://doi.org/10.3390/w14060906
Submission received: 1 October 2021 / Revised: 11 February 2022 / Accepted: 14 February 2022 / Published: 14 March 2022
(This article belongs to the Special Issue Advances in Transboundary Aquifer Assessment)

Abstract

:
Hydrogeomorphology is an emerging discipline that studies the relationship between landforms and hydrology, focusing on groundwater and surface water interactions. This study presents the methodology for the elaboration of a hydro-geomorphological map oriented to illustrate the relationships between the aquifer components and geomorphological characteristics in the United States-Mexico Transboundary San Pedro Aquifer (TSPA). This information contributes to a further understanding of the TSPA, facilitates the location of groundwater recharge and discharge zones, is useful for the development of sustainable groundwater management strategies, and could be useful in developing conceptual and numerical groundwater models for the region.

1. Introduction

Granular and fractured aquifers represent an important source of fresh water in arid and semi-arid regions that are highly dependent on groundwater resources. Factors such as climate, topography, geomorphology, and lithology influence groundwater-flow interactions [1,2,3,4]. However, groundwater availability ultimately relies on the rainfall rate of the site, and the location and characteristics of the aquifer’s recharge and discharge zones [2,3,5]. A better understanding of these areas contributes to the development of groundwater-management plans and strategies that promote water-resources sustainability, which is essential in transboundary settings where water resources are shared by two or more countries. Hydrogeomorphologic studies have proven to be useful for investigating the associations between landforms and hydrological processes that affect surface-water and groundwater flow, identifying the potential impacts of changes in land-use practices, and locating possible groundwater recharge and discharge areas [1,5,6,7,8,9,10,11].
A term first introduced in 1972, hydrogeomorphology broadly described the study of landforms produced by different hydrologic processes [10]. Over the years, hydrogeomorphology evolved into “an interdisciplinary science that focuses on the interaction and linkage of hydrologic processes with landforms or earth materials and the interaction of geomorphic processes with surface and subsurface water in temporal and spatial dimensions” [11]. Frequently, geomorphologic and hydrogeomorphologic studies have focused on flood assessment and surface-water controls, landslide assessment, and in atmosphere-hydrosphere-lithosphere interactions [12,13,14,15]. On the other hand, hydrogeomorphologic studies have also been associated with the analysis of groundwater resources (e.g., [1,5,6,7]). For instance, hydrogeomorphologic mapping allowed the identification and classification of hydro-objects in Southern Italy and the modeling of catchment contribution areas [16]. Additionally, scholars in this area of study have defined the connection between landforms and hydrology and expressed the need for a holistic approach that considers the relationships among landscape engineering, agriculture, natural areas, and water-resources management [17].
The importance of using criteria including landform, elevation, lithology, and hydrology for the assessment of groundwater resources has been described by Chaminé et al. (2015) [18]. In the arid Río Asunción Basin of Sonora, México, the correlation between the basin’s groundwater-storage capacity and its lithologic units, ability to resist weathering and erosion, and presence of faults was described by Gutiérrez Anguamea (2013) [6]. The mapping methodology presented by [6] was later used for the development of a hydrogeomorphologic map for the state of Sonora, Mexico [7], an approach published in collaboration with the Mexican National Water Commission (CONAGUA) that serves as a guide for water-resources management in the region.
Located in Northwestern Mexico, the state of Sonora is bordered to the north by the state of Arizona in the United States. The United States and Mexico share history, culture, people, and water. A recent transboundary-characterization study indicates that based on geological correlations, there are 72 hydrogeologic units, or aquifers, that cross the U.S.-Mexico border [19]. One of these aquifers is the Transboundary San Pedro Aquifer (TSPA), located in the Arizona-Sonora border region. The TSPA is a Transboundary Aquifer Assessment Program (TAAP) aquifer of focus, which is a joint effort between the United States and Mexico to evaluate shared aquifers [20,21,22,23]. A number of studies and technical activities have been carried out through the TAAP in Arizona and Sonora over the last decade (i.e., [22,23,24,25]). For example, in 2016 the International Boundary and Water Commission published the Binational Study of the Transboundary San Pedro Aquifer [24]. This study, jointly developed by TAAP partners from the two countries, binationally described the physical geography, geology, hydrology, hydrogeology, and hydro-geochemistry of the TSPA. In this study, we aim to contribute to the TAAP knowledge base by using hydro-geomorphologic mapping as a tool for groundwater characterization, a novel approach that could guide land and water-management decisions in both the United States and Mexico.

2. Study Area

Located in the eastern portion of the Arizona-Sonora border, the TSPA is drained by the San Pedro River (Figure 1). The San Pedro River has its headwaters east of Cananea, Sonora, and flows northward to the United States until its confluence with the Gila River. The San Pedro River sustains hundreds of species—for example, it is an important bird habitat—and the basin contains one of the major unfragmented landscapes in the Southwest [26]. Several authors who have studied the aquifer basin have reported concerns regarding the impact of groundwater pumping on the San Pedro River, an ongoing issue that has captured the attention of scientists and stakeholders within the region [24,26,27,28].
The TSPA has an approximate area of 5000 km2 and a population of around 97,235 (Table 1). The climate in this border region is arid to semi-arid with bimodal patterns of precipitation characterized by intense summer rains associated with the North American Monsoon and winter precipitation associated with the presence of Pacific cold fronts [24,27]. The mean annual precipitation in the Mexican portion of the TSPA has been reported to be 553 mm [28]. On the other hand, 330 mm were reported in Tombstone, and 960 mm on the Huachuca mountains in Arizona [24]. Mean annual temperature was reported to range between 12 °C and 18 °C [24]. The TSPA is located in the Basin and Range Province, bordering the Sonoran and Chihuahuan Deserts [24,27].
The major economic activities within the region include tourism and military operations in the United States, and livestock, agriculture, and mining in Mexico [24]. According to data from the Mexican Public Registry of Water Rights [32], 82 wells are registered in the Mexican portion of the TSPA under the following activities: 41 for livestock activities, two for industrial uses, 14 for agricultural uses, and 25 for public, urban, residential, and miscellaneous uses. Annual groundwater extractions from these wells equals 30.67 million cubic meters per year (MCM/year) [32], and in 2015, the aquifer registered an annual groundwater deficit of −7.49 MCM [28]. In 2014, groundwater demand in the U.S. portion of the TSPA, the Sierra Vista Sub-watershed, was reported to be 38.38 MCM/year (31,119 acre-feet per year) [33]. This region also reported a groundwater deficit of 5.21 MCM (4229 acre-feet per year) during the same year [33].

3. Materials and Methods

In this study, we analyzed the satellite imagery available through the ArcGIS Online Server [34] and combined the geologic and hydrogeologic information [24], and the topographic features (Digital Terrain Model SRTM1N30W109V3, [35]) to identify composition, topographic arrangement, and the presence or the absence of structures (see Figure 2). A visual inspection of satellite imagery allows for the differentiation of rock units based on the identification of textures and tones, i.e., smoothness, roughness, and compaction [36,37,38].

3.1. Topographic Characterization (Landform Identification)

A topographic characterization was carried out using the Digital Terrain Model information available for the study site [35]. Elements within the study area were classified into four major types of terrain (described below): mountains, hills, piedmonts, and plains [39].
Mountains: Landforms with a relative height greater than 200 m, associated with endogenous folding processes, magmatism, vulcanism, and the dissection of endogenous formation structures [40]. Relative heights were considered from the base to the top of each formation analyzed in this study.
Hills: Landforms with a relative height less than 200 m. This group originates from the leveling of mountains (endogenous) or the dissection of a sloping plain (erosive exogenous). However, hills may be associated with low-elevation endogenous landforms or the product of quaternary tectonics [41].
Piedmonts: Mountainous margins or transitional zones distinguished by a change of slope and considerably lower height, ranging from 0 to 200 m depending on the behavior of the terrain. Piedmonts are composed of detrital material and present fluvial drainage [6,40].
Plains: Land surfaces with minimal slope and altitude difference. Correspond to the cumulative exogenous terrain of alluvial, wind, and coastal deposits [6,40]. The following factors were considered in the identification of a plain: land use (agricultural and urban), change in slope, and drainage pattern.

3.2. Geologic Characterization

A diverse tectonic evolution has shaped a complex geology in the TSPA with intrusive, metamorphic, volcanic-sedimentary, sedimentary, and volcanic rocks [24,30]. A Precambrian basement covered by sedimentary platform sequences—mainly carbonates—is exposed along southeastern Arizona/northeastern Sonora [24,30]. The oldest Mesozoic rocks within this region are Jurassic volcanic and sedimentary sequences covered by Cretaceous-Tertiary rocks, which are widely distributed throughout the TSPA [24,30]. The lithological units considered in this study, based on Callegary et al. (2016) [24], are shown in Table 2.

3.3. Hydrological and Hydrogeological Information

The proposed hydrogeomorphologic map includes hydrological and hydrogeological data for a better visualization of the impact groundwater extractions have on the aquifer’s distinct units. Information for this study includes a spatial layer with the hydrology, the locations of wells, and the groundwater levels for the year 2011. In addition, we identified the permeability, hydraulic conductivity of the rock units based on Gutiérrez Anguamea (2013) [6] and Freeze and Cherry (1979) [3] (Figure 3).
The primary permeability is a property directly related to the origin and formation of rock material to allow water to pass through it [42]; likewise, a secondary permeability can be interpreted based on the number and interconnection of structures that are present in a lithological unit [3,43,44]. Although it is true that the hydraulic potential of materials can be defined by direct and indirect methods—such as petrography, stratigraphy, and resistivity estimation [43]—in this study the permeability was determined based on the characteristics inherent to the formation of the rock (i.e., Figure 3 and Figure 4) and its subsequent fracturing by movements of the earth crust.
It is also presumed that the combination of permeability (primary and/or secondary) with the shape of the terrain is directly related to the potential well yield of a lithological unit (Figure 3). In other words, a portion of materials with a significant primary permeability, such as a smoothed conglomerate hill, is likely to allow water to flow through it easily, as it is composed of elements with varied granulometry and flat topography [6]. In contrast to the above, when it comes to more compact and steeper materials where the speed of surface runoff increases and the spaces between the rock crystals are smaller, the potential well yield of groundwater can be reduced [6].
The combination of the aforementioned factors allowed for the assignment of a groundwater permeability and potential well yield category to each of the elements contained in the TSPA. Categories were based on the aquifer materials, permeability, and hydraulic conductivity to describe how much and how quickly water moves within an aquifer [6].

4. Results

Based on the hydrogeomorphological analysis of the TSPA, a total of 22 units of high, medium, low, very low, and very low/null permeability/potential well yields were characterized (Table 3, Figure 5). The legend of the map is based on [46]. Extensive and highly productive (high aquifer well yield) intergranular aquifer units are shown in shades of blue. The green color range represents fissured environments. The brown areas signify the units of local extension and limited resources, as well as those that are considered to have very low well yield. Alluvial plains, cultivated plains, upper divergent plains, and unconsolidated polymictic conglomerate foothills were identified as intergranular mediums with high permeability and potential well yield.
Regarding the units of medium permeability/potential well yield, only hills of unconsolidated polymictic conglomerate were identified as such. Conglomerate hills, consolidated polymictic conglomerate hills, consolidated basalts, consolidated polymictic conglomerate foothills, consolidated sands, polymictic conglomerate mountains, and polymictic conglomerates were identified as units of limited potential well yield.
The presence of fractures and faults indicates fissured environmental units with a very low permeability and potential well yield. A second subcategory comprises medium permeability fissured units such as fissured limestone, sandstone and shale hill, polymictic conglomerate mountain and fissured sandstone, sandstone mountain, and fissured shale. A third category, the fissured volcanic mountain, was identified in a fissured medium of low permeability and potential well yield. Finally, those units whose permeability is characterized as very low/null are considered non-aquifer units. According to the flow-direction lines identified for the study area, groundwater flows from south to north and towards the San Pedro River. This information is consistent with Callegary et al. (2016) [24], who also identified cones of depressions near the cities of Sierra Vista, Tombstone, and Cananea. For this study, the hydrological discharge zones were located within the San Pedro River and its tributaries, while the recharge areas were mainly located within mountainous areas: the Huachuca, Mule, and Mustang Mountains (in the United States) and the Sierra Mariquita, Sierra Los Ajos, and Sierra San Jose (in Mexico).

5. Discussion and Conclusions

The TSPA is an aquifer shared between the United States and Mexico. It is also a TAAP aquifer of focus that has been deeply studied over the last decade. The Binational Study of the Transboundary San Pedro Aquifer [24] is one of the most relevant binational studies of the region and includes information regarding the physical geography, geology, hydrology, hydrogeology, and geochemistry of the region. According to Chaminé et al. (2015) [18], groundwater characterization must be approached based on different disciplines. These disciplines might include geology, hydrology, hydrogeology, and geochemistry, but also geomorphology and hydrogeomorphology.
Hydrogeomorphology studies describe the interactions between hydrologic processes, landforms, and lithology, and can be useful for determining potential well yields, along with recharge and recharge zones. In this study, we developed a hydrogeomorphologic map for the TSPA. This aquifer is currently experiencing groundwater deficit, with mining, military, domestic, and agricultural users competing for groundwater resources. Hydrogeomorphologic units are defined based on their formation, composition, and original texture of the different rock formations [6]. According to this study, highlands constitute potential recharge zones, and lowlands serve as groundwater-flow discharge areas.
This map makes it possible to quickly identify the functioning of the aquifer system, with the recharge and discharge zones clearly discernible. The information presented here can be used as the basis for the development of sustainable water-resources strategies that consider the hydrogeomorphic characteristics of a given aquifer region to determine how feasible it is to extract or continue extracting water in that area. Moreover, the assessment of binational aquifers needs to use consistent and harmonized methodologies to identify discharge and recharge areas in need of conservation efforts. The application of this methodology allows the locations of these areas to be identified within the framework of a pragmatic morphogenetic mapping.

Author Contributions

Conceptualization, J.I.M.S. and E.M.T.-V.; supervision, J.I.M.S. and E.M.T.-V.; methodology G.A.G.A. and J.I.M.S., validation, J.I.M.S., formal analysis, G.A.G.A., investigation, J.I.M.S. and G.A.G.A., resources, J.I.M.S., R.M., and F.J.G.N., data curation, G.A.G.A., writing—J.I.M.S., E.M.T.-V., G.A.G.A., and R.M., original draft preparation, J.I.M.S., writing—review and editing, E.M.T.-V., visualization, E.M.T.-V. and G.A.G.A., project administration, E.M.T.-V. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially funded by the U.S. Geological Survey (funding authorized by P.L. 109–448) Award Number G17AC00439 for the Transboundary Aquifer Assessment Program.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to acknowledge the Transboundary Aquifer Assessment Program (TAAP).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Teixeira, J.; Chaminé, H.I.; Marques, J.E.; Gomes, A.; Carvalho, J.M.; Albertí, A.P.; Rocha, F.T. Integrated approach of hydrogeomorphology and gis mapping to the evaluation of ground water resources: An example from the hydro-mineral system of caldas da cavaca, nw portugal. In Proceedings of the 33rd International Geological Congress, Oslo, Norway, 6–14 August 2008 . [Google Scholar]
  2. Fetter, C.W. Applied Hydrogeology; Pearson: London, UK, 2000. [Google Scholar]
  3. Freeze, R.A.; Cherry, J.A. Groundwater; Prentice-Hall: Hoboken, NJ, USA, 1979. [Google Scholar]
  4. Surrette, M.; Allen, D.M.; Journeay, M. Regional evaluation of hydraulic properties in variably fractured rock using a hydrostructural domain approach. Hydrogeol. J. 2007, 16, 11. [Google Scholar] [CrossRef]
  5. Teixeira, J.; Chaminé, H.I.; Carvalho, J.M.; Pérez-Alberti, A.; Rocha, F.T. Hydrogeomorphological mapping as a tool in groundwater exploration. J. Maps 2013, 9, 263–273. [Google Scholar] [CrossRef]
  6. Gutiérrez Anguamea, G.A. Metodología para la Elaboración del Mapa Hidrogeomorfológico de la Cuenca del Rio Asunción. Master’s Thesis, Universidad de Sonora, Hermosillo, Mexico, 2013. [Google Scholar]
  7. Gutiérrez Anguamea, G.A.; Minjárez Sosa, J.I.; Oroz Ramos, L.A. Mapa Hidrogeomorfológico del Estado de Sonora; CONAGUA: Hermosillo, Mexico, 2015. [Google Scholar]
  8. Babar, M. Hydrogeomorphology: Fundamentals, Applications and Techniques; New India Publishing Agency: New Delhi, India, 2005. [Google Scholar]
  9. Leopold, L.B. Field data: The interface between hydrology and geomorphology, scientific basis of water-resource management. In Studies in Geophysics; The National Academies Press: Cambridge, MA, USA, 1982. [Google Scholar]
  10. Scheidegger, A.E. Hydrogeomorphology: Review. J. Hydrol. 1973, 20, 193–215. [Google Scholar] [CrossRef]
  11. Sidle, R.C.; Onda, Y. Hydrogeomorphology: Overview of an emerging science. Hydrol. Processes 2004, 18, 597–602. [Google Scholar] [CrossRef]
  12. Gregory, K.J. Hydrogeomorphology: How applied should we become? Prog. Phys. Geogr. 2016, 3, 84–101. [Google Scholar] [CrossRef]
  13. Barrera, F.B.; Borja, C.; Núñez, M.F.; Lama, A. Dinámica hidrogeomorfológica e impacto antrópico en la cuenca del arroyo del partido (NW del Parque Nacional de Doñana, Huelva, España): Evaluación de procesos actuales. Cuatern. Geomorfol. Rev. Soc. Española Geomorfol. Asoc. Española Para Estud. Cuatern. 2009, 23, 45–64. [Google Scholar]
  14. Arnaud-Fassetta, G.; Cossart, E.; Fort, M. Hydro-geomorphic hazards and impact of man-made structures during the catastrophic flood of June 2000 in the Upper Guil catchment (Queyras, Southern French Alps). Geomorphology 2005, 66, 41–67. [Google Scholar] [CrossRef]
  15. Verstappen, H.T.; van Zuidam, R.A. The ITC System of Geomorphologic Survey: A Basis for the Evaluation of Natural Resources and Hazards; ITC Publication: Murarrie, Australia, 1981. [Google Scholar]
  16. Guida, D.; Cuomo, A.; Palmieri, V. Using object-based geomorphometry for hydro-geomorphological analysis in aMediterranean research catchment. Hydrol. Earth Syst. Sci. 2016, 20, 3493–3509. [Google Scholar] [CrossRef] [Green Version]
  17. Kudrna, K.; Šindelářová, M. Principles of hydrogeomorphology as a basic precondition for solution of territorial structure of unitary system of agricultural, forest and water management zákony hydrogeomorfologie jako základní předpoklad pro řešení teritoriální struktury unitární soustavy zemědělského, lesního a vodního hospodářství. J. Cent. Eur. Agric. 2006, 7, 669–676. [Google Scholar]
  18. Chaminé, H.I.; Carvalho, J.M.; Teixeira, J.; Freitas, L. Role of hydrogeological mapping in groundwater practice: Back to basics. Eur. Geol. J. 2015, 40, 34–42. [Google Scholar]
  19. Sanchez, R.; Rodriguez, L. Transboundary Aquifers between Baja California, Sonora and Chihuahua, Mexico, and California, Arizona and New Mexico, United States: Identification and Categorization. Water 2021, 13, 2878. [Google Scholar] [CrossRef]
  20. IBWC. Joint Report of the Principal Engineers Regarding the Joint Cooperative Process United States-Mexico for the Transboundary Aquifer Assessment Program; International Boundary and Water Comission: El Paso, TX, USA, 2009. [Google Scholar]
  21. Tapia-Villaseñor, E.M.; Megdal, S.B. The U.S.-Mexico Transboundary Aquifer Assessment Program as a Model for Transborder Groundwater Collaboration. Water 2021, 13, 530. [Google Scholar] [CrossRef]
  22. Callegary, J.B.; Megdal, S.B.; Tapia Villaseñor, E.M.; Petersen-Perlman, J.D.; Minjárez Sosa, I.; Monreal, R.; Gray, F.; Grijalva Noriega, F. Findings and lessons learned from the assessment of the Mexico-United States transboundary San Pedro and Santa Cruz aquifers: The utility of social science in applied hydrologic research. J. Hydrol. Reg. Stud. 2018, 20, 60–73. [Google Scholar] [CrossRef]
  23. Petersen-Perlman, J.D.; Albrecht, T.R.; Tapia-Villaseñor, E.M.; Varady, R.G.; Megdal, S.B. Science and Binational Cooperation: Bidirectionality in the Transboundary Aquifer Assessment Program in the Arizona-Sonora Border Region. Water 2021, 13, 2364. [Google Scholar] [CrossRef]
  24. Callegary, J.B.; Minjarez, I.; Tapia Villaseñor, E.M.; dos Santos, P.; Monreal Saavedra, R.; Grijalva Noriega, F.J.; Scott, C.A.; Megdal, S.B.; Oroz Ramos, L.A.; Rangel Medina, M.; et al. Binational Study of the Transboundary San Pedro Aquifer; International Boundary and Water Comission: El Paso, TX, USA, 2016. [Google Scholar]
  25. Tapia-Villaseñor, E.M.; Shamir, E.; Megdal, S.B.; Petersen-Perlman, J.D. Impacts of Variable Climate and Effluent Flows on the Transboundary Santa Cruz Aquifer. JAWRA J. Am. Water Resour. Assoc. 2020, 56, 409–430. [Google Scholar] [CrossRef]
  26. The Nature Conservancy. San Pedro River; The Nature Conservancy: Arlington, VA, USA, 2022. [Google Scholar]
  27. Minjarez Sosa, I.; Monreal Saavedra, R.; Rangel Medina, M.; Grijalva Noriega, F.; Ochoa Granillo, A.; Tapia Villaseñor, E.M.; Ramírez Días, L.L.; Cirett Galán, S.; Moncada García, G.; Montijo González, A. Actividades Hidrogeológicas en el Acuífero Río San Pedro, Estado de Sonora, Primera Etapa; Universidad de Sonora, CILA, CONAGUA: Hermosillo, Mexico, 2011. [Google Scholar]
  28. CONAGUA. Actualización de la Disponibilidad Media Anual de Agua en el Acuífero Río San Pedro (2616), Estado de Sonora; Diario Oficial de la Federación: Ciudad de México, Mexico, 2015.
  29. Instituto Nacional de Estadística y Geografía (INEGI). Censo de Población y Vivienda 2010; INEGI: Ciudad de México, Mexico, 2010; Available online: https://www.inegi.org.mx/programas/ccpv/2010/ (accessed on 30 September 2021).
  30. Instituto Nacional de Estadística y Geografía (INEGI). Censo de Población y Vivienda 2020; INEGI: Ciudad de México, Mexico, 2020; Available online: https://www.inegi.org.mx/programas/ccpv/2020/ (accessed on 30 September 2021).
  31. United States Census Bureau. Census Bureau QuickFacts: Sierra Vista City, Arizona—Tucson City, Arizona. 2022. Available online: https://www.census.gov/quickfacts/fact/table/sierravistacityarizona,tucsoncityarizona/BZA210219 (accessed on 30 September 2021).
  32. CONAGUA. Consulta a la Base de Datos del Registro Público de Derechos de Agua (REPDA). Available online: https://www.gob.mx/conagua/acciones-y-programas/consulta-la-base-de-datos-del-repda (accessed on 12 May 2021).
  33. Watershed Management Group; City of Sierra, V. A Stormwater Action Plan for Sierra Vista. 2016. Available online: https://verderiver.org/wp-content/uploads/2017/07/stormwater-action-plan-for-sierra-vista.pdf (accessed on 30 September 2021).
  34. Esri; Maxar; GeoEye; Earthstar Geographics; CNES/Airbus DS; USDA; USGS; AeroGRID; IGN; GIS User Community. World Map. 2021. Available online: https://services.arcgisonline.com/ArcGIS/rest/services/World_Imagery/MapServer (accessed on 30 September 2021).
  35. Farr, T.G.; Rosen, P.A.; Caro, E.; Crippen, R.; Duren, R.; Hensley, S.; Kobrick, M.; Paller, M.; Rodriguez, E.; Roth, L.; et al. The Shuttle Radar Topography Mission. Rev. Geophys. 2007, 45, RG2004. [Google Scholar] [CrossRef] [Green Version]
  36. Panizza, M. Environmental Geomorphology. Developments in Earth Surface Processes; Elsevier: Amsterdam, The Netherlands, 1996; Volume 4. [Google Scholar]
  37. Muñoz Jiménez, J. Paisaje y geosistema, una aproximación desde la geografía física. In Seminario “Paisaje y Medio Ambiente”; Universidad de Valladolid, Secretariado de Publicaciones e Intercambio Editorial: Valladolid, Spain, 1996. [Google Scholar]
  38. García-Romero, A. An Evaluation of Forest Deterioration in the Disturbed Mountains of Western Mexico City. Mt. Res. Dev. 2002, 22, 270–277. [Google Scholar] [CrossRef]
  39. Tapia, V.G.; Lopez, B. Mapeo geomorfológico analítico de la porción central de la Cuenca de México: Unidades morfogenéticas a escala 1:100,000. Rev. Mex. Cienc. Geológicas 2002, 19, 50–65. [Google Scholar]
  40. Lugo Hubp, J. Diccionario Geomorfológico; Instituto de Geografía, Universidad Nacional Autónoma de México: Ciudad de México, Mexico, 1989. [Google Scholar]
  41. Córdova, C. Cartografía Geomorfológica en Escala Pequeña del Occidente de la República Mexicana. Bachelor’s Thesis, Facultad de Filosofía y Letras, Universidad Autónoma de México, Ciudad de México, Mexico, 1988. [Google Scholar]
  42. Castany. Tratado Práctico de las Aguas Subterráneas. Primera Parte. Elementos de Hidrología de Superficie; Ed. Omega: Barcelona, Spain, 1971. [Google Scholar]
  43. Pedroza, V.; Gutierrez, M.; Le Roux, J.; Hernández, E.; Vasquez, J.P. Estudio Comparativo del Coeficiente de Permeabilidad K. Available online: https://biblioteca.sernageomin.cl/opac/DataFiles/14905_v2_pp_238_240.pdf (accessed on 30 September 2021).
  44. Díaz, E.J.; Contreras, N.M.; Pinto, J.E.; Velandia, F.; Morales, C.J.; Hincapie, G. Evaluación hidrogeológica preliminar de las unidades geológicas de la mesa de los santos, santander. Boletín Geol. 2009, 31, 10. [Google Scholar]
  45. LeGrand, H. Geology and Ground Water in the Statesville Area; North Carolina Department Conservation Development, Divison Mineral Resources: Raleigh, NC, USA, 1954. [Google Scholar]
  46. UNESCO. International Legend for Hidrogeological Maps; UNESCO: Paris, France, 1983. [Google Scholar]
Figure 1. The Transboundary San Pedro Aquifer (TSPA). Source: Authors’ development based on Callegary et al. (2016) [24].
Figure 1. The Transboundary San Pedro Aquifer (TSPA). Source: Authors’ development based on Callegary et al. (2016) [24].
Water 14 00906 g001
Figure 2. GIS layers used for the development of the hydrogeomorphologic map of the TSPA. (a) Digital Elevation Data [35]; (b) Satellite Image [34]; (c) Geology and Structures [24]; (d) Hydrography and Groundwater Levels [24].
Figure 2. GIS layers used for the development of the hydrogeomorphologic map of the TSPA. (a) Digital Elevation Data [35]; (b) Satellite Image [34]; (c) Geology and Structures [24]; (d) Hydrography and Groundwater Levels [24].
Water 14 00906 g002aWater 14 00906 g002b
Figure 3. Permeability and hydraulic conductivity of rock units and unconsolidated deposits. Modified from Freeze and Cherry (1979) [3].
Figure 3. Permeability and hydraulic conductivity of rock units and unconsolidated deposits. Modified from Freeze and Cherry (1979) [3].
Water 14 00906 g003
Figure 4. Relation between landform and potential well yield. Source: Freeze and Cherry (1979) [3] and LeGrand (1954) [45].
Figure 4. Relation between landform and potential well yield. Source: Freeze and Cherry (1979) [3] and LeGrand (1954) [45].
Water 14 00906 g004
Figure 5. Hydrogeomorphologic map of the San Pedro River Basin.
Figure 5. Hydrogeomorphologic map of the San Pedro River Basin.
Water 14 00906 g005
Table 1. Population centers located in the TSPA. Source: INEGI, 2010, INEGI, 2020, U.S. Census Bureau 2020 [29,30,31].
Table 1. Population centers located in the TSPA. Source: INEGI, 2010, INEGI, 2020, U.S. Census Bureau 2020 [29,30,31].
TownPopulation
Sierra Vista45,308
Tombstone1209
Naco6064
Bisbee5203
Cananea39,451
Total97,235
Table 2. Lithological units in the TSPA. Source: Authors’ development based on Callegary et al. (2016) [24].
Table 2. Lithological units in the TSPA. Source: Authors’ development based on Callegary et al. (2016) [24].
LegendLithostratigraphic UnitsDescription
Water 14 00906 i001Precambrian Igneous-Metamorphic
Complex
Igneous and metamorphic rocks
Water 14 00906 i002Early Paleozoic Sedimentary Unit Localized outcrops of detrital-carbonate rocks within the Mexican portion of the TSPA
Water 14 00906 i003Late Paleozoic Sedimentary Unit Limestone and sandstone exposed in most
topographic highs in the TSPA
Water 14 00906 i004Jurassic Felsic Volcano-Sedimentary Intercalation of volcanic rocks, sandstones,
agglomerates, basalt flows, sills, and intermediate composition
Water 14 00906 i005Jurassic Intrusive Complex Intrusive hypabyssal bodies mainly exposed on the U.S. side of the TSPA
Water 14 00906 i006Late Jurassic–Early Cretaceous
Sedimentary Unit
Conglomerate, sandstone, shale, and limestone from the Bisbee Group
Water 14 00906 i007Late Cretaceous Sedimentary Unit
(KsVs, Ks)
Sedimentary sequences
Water 14 00906 i008Cretaceous–Paleocene
Volcano-Sedimentary Unit
Rhyolitic clastic and volcanic rocks
Water 14 00906 i009Tertiary–Cretaceous Intrusive Complex Intrusive felsic rocks
Water 14 00906 i010Tertiary Felsic Volcanic Unit Rhyolitic rocks from the west-central portion of the TSPA
Water 14 00906 i011Tertiary Volcano-Sedimentary Unit Continental rocks, mainly conglomerates with intercalations of sandstone and tuff
Water 14 00906 i012Plio–Quaternary Sedimentary Unit Coarse sediments (gravels and sands) distributed within the center of the TSPA
Water 14 00906 i013AlluviumGravel, sands, silts, and clay.
Table 3. Hydrogeomorphologic Units in the Transboundary San Pedro Aquifer.
Table 3. Hydrogeomorphologic Units in the Transboundary San Pedro Aquifer.
Recharge/Non-RechargeUnitDescription
Discharge ZoneHigh permeability intergranular environment (high potential well yield)Crop plain
Upper divergent plain
Alluvial plain
Unconsolidated polymictic conglomerate piedmont
Water
Medium permeability environment (medium potential well yield)Unconsolidated polymictic conglomerate hill
Low permeability environment (low potential yield)Consolidated polymictic conglomerate hill
Consolidated polymictic conglomerate and basalt hill
Recharge ZoneMedium permeability fissured environment (medium potential well yield)Fissured limestone, sandstone, and shale hill
Fissured sandstone and shale mountain
Low permeability fissured environment (low potential well yield)Fissured volcanic mountain
Very low permeability fissured environment (very low potential well yield)Fissured polymictic conglomerate and volcanic mountain
Fissured limestone, sandstone, and shale mountain
Fissured volcanic and sandstone mountain
Fissured plutonic mountain
Fissured metamorphic mountain
Fissured volcanic hill
Fissured plutonic hill
Fissured metamorphic hill
Impervious AreasUrban zone
Volcanic hill
Plutonic hill
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sosa, J.I.M.; Gutiérrez Anguamea, G.A.; Monreal, R.; Grijalva Noriega, F.J.; Tapia-Villaseñor, E.M. Hydrogeomorphologic Mapping of the Transboundary San Pedro Aquifer: A Tool for Groundwater Characterization. Water 2022, 14, 906. https://doi.org/10.3390/w14060906

AMA Style

Sosa JIM, Gutiérrez Anguamea GA, Monreal R, Grijalva Noriega FJ, Tapia-Villaseñor EM. Hydrogeomorphologic Mapping of the Transboundary San Pedro Aquifer: A Tool for Groundwater Characterization. Water. 2022; 14(6):906. https://doi.org/10.3390/w14060906

Chicago/Turabian Style

Sosa, José Ismael Minjárez, Grisel Alejandra Gutiérrez Anguamea, Rogelio Monreal, Francisco Javier Grijalva Noriega, and Elia M. Tapia-Villaseñor. 2022. "Hydrogeomorphologic Mapping of the Transboundary San Pedro Aquifer: A Tool for Groundwater Characterization" Water 14, no. 6: 906. https://doi.org/10.3390/w14060906

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop