Estimating ground water recharge using flow models of perched karstic aquifers

The fraction of rain that is annually recharged to ground water is a function of the transient quantities of precipitation (wet vs. dry years) as well as other meteorological and geologic factors, and thus it is very difficult to estimate. In this study, we have used long records (20 to 30 years) of precipitation and spring discharge to reconstruct the transient character of yearly recharge. These data sets were used to calibrate numerical ground water flow models on the less than 3 km(2) scale for four separate perched karstic aquifers in the Judean and Samarian Mountains of Israel. The stratification and karstic character of the local carbonate rock aquifers cause ground water to flow through discrete dissolution channels and to discharge at isolated springs. An innovative, dual-porosity approach was used where a finite-difference solution simulates flow in the rock matrix, while the karstic channels are simulated using computationally simple drains. Perched conditions are also simulated innovatively using MODFLOW by treating the bottom unsaturated layer as if it is saturated, but by assuming zero pressure head throughout the "unsaturated" layer. Best fitting between measured and computed spring hydrograph data has allowed us to develop a set of empirical functions relating measured precipitation to recharge to the aquifer. The generic methodology presented gives insight into the suspected changes in aquifer recharge rates between particularly wet or dry years.     

Spatial evolution of nitrogen and phosphorus loads along a small Mediterranean river: implication of bed sediments

This study addresses the spatial variations in water quality along the River Vène (France). The Vène drains a 67 km² rural basin, with a large karstic area, located in a Mediterranean context. A 1 day sampling campaign was conducted along the river, in winter low‐flow conditions (February 2003). Physico‐chemical parameters and water flow discharge were measured in situ during the sampling campaign. Water quality was evaluated by determining the concentrations of nitrogen and phosphorus in water and bed‐sediment samples. Nitrogen and phosphorus loads were evaluated taking into account the measured concentrations and discharge. The campaign included 18 sampling points and concerned the whole river from the spring to the outlet, plus the main inputs, i.e. sewage treatment works, main tributaries and karstic springs. The spatial evolution of nitrogen and phosphorus loads along the river allowed the significant role of point‐source inputs to be demonstrated. The decrease in nutrient loads along the river occurred mainly in specific reaches where fine sediments had accumulated. In these zones, phosphorus is trapped in the bed sediments in calcium‐bound phosphates due to precipitation processes. Copyright © 2005 John Wiley & Sons, Ltd.     

Conceptualization of Karstic Aquifer with Multiple Outlets Using a Dual Porosity Model

In this study, two conceptual models, the classic reservoir (CR) model and exchange reservoirs model embedded by dual porosity approach (DPR) are developed for simulation of karst aquifer functioning drained by multiple outlets. The performances of two developed models are demonstrated at a less developed karstic aquifer with three spring outlets located in Zagros Mountain in the south‐west of Iran using 22‐years of daily data. During the surface recharge, a production function based on water mass balance is implemented for computing the time series of surface recharge to the karst formations. The efficiency of both models has been assessed for simulation of daily spring discharge during the recession and also surface recharge periods. Results indicate that both CR and DPR models are capable of simulating the ordinates of spring hydrographs which drainage less developed karstic aquifer. However, the goodness of fit criteria indicates outperformance of DPR model for simulation of total hydrograph ordinates. In addition, the DPR model is capable of quantifying hydraulic properties of two hydrologically connected overlapping continua conduits network and fissure matrix which lays important foundations for the mining operation and water resource management whereas homogeneous model representations of the karstic subsurface (e.g., the CR) do not work accurately in the karstic environment.     

Importance of unsaturated zone flow for simulating recharge in a humid climate

Transient recharge to the water table is often not well understood or quantified. Two approaches for simulating transient recharge in a ground water flow model were investigated using the Trout Lake watershed in north-central Wisconsin: (1) a traditional approach of adding recharge directly to the water table and (2) routing the same volume of water through an unsaturated zone column to the water table. Areas with thin (less than 1 m) unsaturated zones showed little difference in timing of recharge between the two approaches; when water was routed through the unsaturated zone, however, less recharge was delivered to the water table and more discharge occurred to the surface because recharge direction and magnitude changed when the water table rose to the land surface. Areas with a thick (15 to 26 m) unsaturated zone were characterized by multimonth lags between infiltration and recharge, and, in some cases, wetting fronts from precipitation events during the fall overtook and mixed with infiltration from the previous spring snowmelt. Thus, in thicker unsaturated zones, the volume of water infiltrated was properly simulated using the traditional approach, but the timing was different from simulations that included unsaturated zone flow. Routing of rejected recharge and ground water discharge at land surface to surface water features also provided a better simulation of the observed flow regime in a stream at the basin outlet. These results demonstrate that consideration of flow through the unsaturated zone may be important when simulating transient ground water flow in humid climates with shallow water tables.     

Locating the zone of saline intrusion in a coastal karst aquifer using springflow data

Coastal fresh water aquifers are an increasingly desirable resource. In a karstic aquifer, sea water intrusion occurs as a salt water wedge, like in porous media. However, preferential flow conduits may alter the spatial and temporal distribution of the salt water. This is typically the case when the outlet of the aquifer is a brackish spring. This paper shows that salinity and flow rate variations at a spring, where salinity is inversely proportional to discharge, can help to understand the hydrodynamic functioning of the aquifer and to locate the fresh water-sea water mixing zone deep inside the aquifer. The volume of water-filled conduit between the sea water intrusion zone and the spring outlet is calculated by the integral over time of the flow rate during the time lag between the flow rate increase and the salinity decrease as measured at the spring. In the example of the spring at Almyros of Heraklio (Crete, Greece), this time lag is variable, depending on the discharge, but the volume of water-filled conduit appears to be constant, which shows that the processes of salt water intrusion and mixing in the conduit are constant throughout the year. The distance between the spring and the zone where sea water enters the conduit is estimated and provides an indication of the position where only fresh water is present in the conduit.     

Brackish karstic springs model: application to Almiros spring in Crete

A mathematical model is proposed to simulate brackish karstic springs. Rainfall data constitutes model input information while output information is the discharge and the chloride concentration of the water versus time. The model was constructed by considering the mass and mechanical energy balance on the hydrodynamic analog, which includes three reservoirs outflowing in a tube that lies adjacent to the spring. Two reservoirs emulate the karstic system, and the third one emulates the sea. The discharge of the spring is given by the sum of the discharge of the reservoirs, and the chloride concentration by the solution of the mixing problem between the fresh and the salty water, which exists in the tube leading to the spring. The model is applied to the spring of Almiros at Heraklion, Crete, Greece. The agreement between model values and field measurements is very good for depletion periods and satisfactory for recharge periods.     

Impact of water transfer on interaction between surface water and groundwater in the lowland area of North China Plain

The contradiction between the freshwater shortage and the large demand of freshwater by irrigation was the key point in cultivated lowland area of North China Plain. Water transfer project brings fresh water from water resource‐rich area to water shortage area, which can in turn change the hydrological cycle in this region. Major ions and stable isotopes were used to study the temporal variations of interaction between surface water and groundwater in a hydrological year after a water transfer event in November 2014. Irrigation canal received transferred Yellow River, with 2.9% loss by evaporation during water transfer process. The effect of transferred water on shallow groundwater decreased with increasing distance from the irrigation canal. Pit pond without water transfer receives groundwater discharge. During dry season after water transfer event, shallow groundwater near the irrigation canal was recharged by lateral seepage and deep percolation of irrigation, whereas shallow groundwater far from irrigation canal was recharged by deep percolation of deep groundwater irrigation. Canal water lost by evaporation was 2.7–17.4%. Influence of water transfer gradually disappeared until March as the water usage of agricultural irrigation increased. In the dry season, groundwater discharged to irrigation canal and pond; 2.2–31.6% canal water and 11.3–20.0% pond water were lost by evaporation. In the rainy season (June to September), surface water was fed mainly by precipitation and surface run‐off, whereas groundwater was recharged by infiltration of precipitation. The two‐end member mix model showed that the mixing ratio of precipitation in pond and irrigation canal were 73–83.4% (except one pond with 28.1%) and 77.3–99.9%, respectively. Transferred water and precipitation were the important recharge sources for shallow groundwater, which decreased groundwater salinity in cultivated lowland area of North China Plain. With the temporary and spatial limitation of water transfer effects, increased water transfer amounts and frequency may be an effective way of mitigating regional water shortage. In addition, reducing the evaporation of surface water is also an important way to increase the utilization of transfer water.     

Variations of phytoplankton community structure related to water quality trends in a tropical karstic coastal zone

Phytoplankton community structure in coastal areas is a result of various environmental factors such as nutrients, light, grazing, temperature, and salinity. The Yucatan Peninsula is a karstic tropical region that is strongly influenced by submerged groundwater discharge (SGD) into the coastal zone. Phytoplankton community structure and its relationship with regional and local water quality variables were studied in four ports of the northwestern Yucatan Peninsula. Water quality was strongly related to SGD, and variations in phytoplankton community structure were related to local nutrient loading and hydrographic conditions, turbulence, and human impacts. Our study provides an ecological baseline for the Yucatan Peninsula and serves as a basis for establishing monitoring programs to predict changes at sites with high hydrological variation and in developing an early alert system for harmful toxic algal blooms.     

Climate change impact on the hydrological budget of a large Mediterranean island

ABSTRACT

Crete is a Mediterranean, karst-dominated island, characterized by long drought periods. The Karst-SWAT model, combined with 11 climate change scenarios, was run to assess climate change impacts on the island under two set-ups, both using the auto-irrigation function of the model: (1) with water drawn from the shallow or deep aquifer, and (2) with irrigated water derived from an unlimited outside source. The first set-up provided insight into the fluctuation of future irrigation needs, and when compared to the second set-up, enabled quantification of the future water deficit. The Water Exploitation Index was used to describe the spatial variability of future water stress on Crete. A decrease in both surface and karstic spring flows is foreseen, especially after 2060 (24.2 and 16.5%, respectively). Simulated irrigation water demand and water deficit show continuous increase throughout the projection period (2020–2098).     

Multifractal behaviour of long‐term karstic discharge fluctuations

Karstic watersheds are highly complex hydrogeological systems that are characterized by a multiscale behaviour corresponding to the different pathways of water in these systems. The main issue of karstic spring discharge fluctuations consists in the presence and the identification of characteristic time scales in the discharge time series. To identify and characterize these dynamics, we acquired, for many years at the outlet of two karstic watersheds in South of France, discharge data at 3‐mn, 30‐mn and daily sampling rate. These hydrological records constitute to our knowledge the longest uninterrupted discharge time series available at these sampling rates. The analysis of the hydrological records at different levels of detail leads to a natural scale analysis of these time series in a multifractal framework. From a universal class of multifractal models based on cascade multiplicative processes, the time series first highlights two cut‐off scales around 1 and 16 h that correspond to distinct responses of the aquifer drainage system. Then we provide estimates of the multifractal parameters α and C₁ and the moment of divergence q_D corresponding to the behaviour of karstic systems. These results constitute the first estimates of the multifractal characteristics of karstic spingflows based on 10 years of high‐resolution discharge time series and should lead to several improvements in rainfall‐karstic springflow simulation models. Copyright © 2012 John Wiley & Sons, Ltd.     

Stable isotope‐based mean catchment altitudes of springs in the Lebanon Mountains

Stable isotopes of water are known to provide information on mean altitudes of spring recharge areas, which is an important parameter for groundwater resources management especially in karstic environments. Very often, a lack of precipitation input data limits the possibility for an appropriate estimation of mean catchment altitudes. In the Jeita spring catchment, Lebanon, a characterization of precipitation input was possible with samples collected at six stations at varying altitudes (88 amount‐weighted monthly samples). A local meteoric water line for the Jeita spring catchment was characterized as δ²H = 6.04 * δ¹⁸O + 8.45 (R ² = .92) for a 2‐year observation period between October 2012 and September 2014. Integral samples from the snow layer were collected at 22 sites at altitudes ranging from 1,000 to 2,300 m above sea level at the end of February 2012 and February 2013, when snow height reached a maximum of more than 6 m at the highest peak in the catchment. Water samples were continuously collected from six springs (Jeita, Kashkoush, Labbane, Assal, Afqa, and Rouaiss). Jeita spring water samples were collected additionally in daily time steps during the snowmelt season in 2012. Mean isotope values of the sampled springs range from ?6.8‰ to ?8.2‰, and from ?33‰ to ?44‰, for δ¹⁸O and δ²H, respectively. The stable isotope data show that input variability (space and time, snow cover, and rainfall) has direct impacts on mean altitude estimates of spring catchments. A more profound interpretation of spring response to rainfall for six local springs in the Lebanon Mountains was possible in comparison to four earlier described springs collected in the Anti‐Lebanon Mountains in Syria.     

A MODFLOW Infiltration Device Package for Simulating Storm Water Infiltration

This article describes a MODFLOW Infiltration Device (INFD) Package that can simulate infiltration devices and their two‐way interaction with groundwater. The INFD Package relies on a water balance including inflow of storm water, leakage‐like seepage through the device faces, overflow, and change in storage. The water balance for the device can be simulated in multiple INFD time steps within a single MODFLOW time step, and infiltration from the device can be routed through the unsaturated zone to the groundwater table. A benchmark test shows that the INFD Package's analytical solution for stage computes exact results for transient behavior. To achieve similar accuracy by the numerical solution of the MODFLOW Surface‐Water Routing (SWR1) Process requires many small time steps. Furthermore, the INFD Package includes an improved representation of flow through the INFD sides that results in lower infiltration rates than simulated by SWR1. The INFD Package is also demonstrated in a transient simulation of a hypothetical catchment where two devices interact differently with groundwater. This simulation demonstrates that device and groundwater interaction depends on the thickness of the unsaturated zone because a shallow groundwater table (a likely result from storm water infiltration itself) may occupy retention volume, whereas a thick unsaturated zone may cause a phase shift and a change of amplitude in groundwater table response to a change of infiltration. We thus find that the INFD Package accommodates the simulation of infiltration devices and groundwater in an integrated manner on small as well as large spatial and temporal scales.     

Transport of suspended solids from a karstic to an alluvial aquifer: the role of the karst/alluvium interface

This study focuses on the coupled transport of dissolved constituents and particulates, from their infiltration on a karst plateau to their discharge from a karst spring and their arrival at a well in an alluvial plain. Particulate markers were identified and the transport of solids was characterised in situ in porous and karstic media, based on particle size analyses, SEM, and traces. Transport from the sinkhole to the spring appeared to be dominated by flow through karst: particulate transport was apparently conservative between the two sites, and there was little difference in the overall character of the particle size distribution of the particulates infiltrating the sinkhole and of those discharging from the spring. Qualitatively, the mineralogy of the infiltrating and discharging material was similar, although at the spring an autochthonous contribution from the aquifer was noted (chalk particles eroded from the parent rock by weathering). In contrast, transport between the spring and the well appears to be affected by the overlying alluvium: particles in the water from the well, showed evidence of considerable size-sorting. Additionally, SEM images of the well samples showed the presence of particles originating from the overlying alluvial system; these particles were not found in samples from the sinkhole or the spring. The differences between the particulates discharging from the spring and the well indicate that the water pumped from the alluvial plain is coming from the karst aquifer via the very transmissive, complex geologic interface between the underlying chalk formation and the gravel at the base of the overlying alluvial system.     

Improvement of the solute transfer in a conceptual unsaturated zone scheme: a case study of the Seine River basin

For predicting the evolution of solute concentrations in groundwater and testing the impact of remediation policies, a coupling between the agronomical model STICS and the hydrogeological model MODCOU was implemented. When applied to the Seine River basin, this model accurately represents the temporal evolution of average nitrate concentrations in the aquifer, but with large local errors. We propose an improvement to the simple unsaturated zone (UZ) scheme NonsatSW used in STICS–MODCOU. The modifications are based on a comparison with the mechanistic model Metis considered as a reference as it solves Richards' equation. A more realistic saturation profile and a varying percolation rate are integrated in NonsatSW. This new model, named NonsatVG, is assessed by comparing it with NonsatSW and Metis. In an ideal case, NonsatVG generates a solute transfer and a dispersion closer to that of Metis than of NonsatSW. In real cases, without additional calibration, NonsatVG and Metis simulate better the average transfer velocities of the observed nitrate profiles. Furthermore, modifications in NonsatVG give a direct relationship between the depth of the water table and the saturation profile. We obtain, therefore, as in Metis, an evolution of the solute transfer velocity depending on the piezometric level. These dynamics are not simulated in NonsatSW. Despite a modified water transfer through the UZ, NonsatVG is also as valid as NonsatSW in the modelling of water transfer to the saturated zone. Finally, an application to the Seine basin shows that solute transfer velocities are lower with NonsatVG than with NonsatSW, but are in better agreement with literature. Copyright © 2010 John Wiley & Sons, Ltd.     

Neural network simulation of spring flow in karst environments

Daily discharges of two springs lying in a karstic environment were simulated for a period of 2.5 years with the use of a multi-layer perceptron back-propagation neural network. Two models were developed for the springs, one relying on the original data and another where the missing discharge values were supplemented by assuming linear relationships during base flow conditions. For both springs the mean square error of the two models did not differ significantly, with an improvement exhibited at the extremes, during the network’s training phase, by the model that utilized the extended data set, the results of which are reported here. The time lag between precipitation and spring discharge differed significantly for the two springs indicating that in karstic environments hydraulic behavior is dominated, even within a few hundred meters, by local conditions. Optimum training results were attained with a Levenberg–Marquardt algorithm resulting in a network architecture consisting of two input layer neurons, four hidden layer neurons, and one output layer neuron, the spring’s discharge. The neural network’s predictions captured the behavior for both springs and followed very closely the discontinuities in the discharge time series. Under-/over-estimation of observed discharges for the two springs remained below 3 %, with the exception of a few local maxima where the predicted discharges diverged more strongly from observed values. Inclusion of temperature data did not add to the improvement of predictions. Finally, optimum predictions were attained when past discharge data were added to the input record and discharge differentials rather than direct discharges were calculated resulting in elimination of any local maximum discrepancy between observed and predicted discharge values.     

鄱阳湖典型洲滩湿地水分补排关系

湿地水分在地下水含水层-土壤-植物-大气界面的运移和转换是维持能量和营养物平衡的重要环节,水分运移是湿地生态水文过程研究的关键.数值模型模拟已成为水分运移研究的重要手段,然而限于复杂的湿地自然条件及有限的监测手段,部分界面水分通量连续动态变化数据的获取及定量化工作较为困难,目前应用数值模拟法于湿地水分运移研究的案例仍不多见.本文以鄱阳湖典型湿地为研究区,构建垂向一维数值模型,阐释了湖泊水位显著季节性变化条件下,湿地水分在不同界面的传输过程,量化了湿地水分的补排关系.结果表明:(1)界面水分通量季节性差异大,降雨入渗地面和根系层水分渗漏均对降雨变化响应敏感,主要集中在4—6月,分别占年总量(1450和1053 mm)的65%和73%.土面蒸发和植物蒸腾年总量为176和926 mm,土面蒸发主要受气候条件影响,植物蒸腾还与植物生长特征有关,均集中在7—8月,分别占年总量的30%和47%.深层土壤向浅层根系层的水分补给集中发生在地下水浅埋时段6—8月,占年总量(609 mm)的76%;(2)湿地植物根系层水分补排受鄱阳湖水位季节性波动影响显著.除丰水期(7—9月)主要补给为深层土壤水外,退、枯、涨水期的主要补给均为降水入渗.涨水期(4—6月)和枯水期(12—3月)的主要排泄为根系层水分渗漏,丰水期以植物蒸腾排泄为主,退水期(10—11月),土面蒸发与植物蒸腾为主要排泄,且比重相当.本文定量了鄱阳湖典型湿地不同界面水分连续交换关系,区分了土面蒸发和植物蒸腾,辨析了各界面水分的主要影响因子,研究结果有助于深入理解水分在湿地生态系统地下水含水层-土壤-植物-大气界面的相互作用机制,认识湖泊洲滩湿地水量平衡,为揭示湖泊水情变化对湿地生态的可能影响提供依据,为湿地生态水文过程研究提供重要方法和理论参考.     

The discovery of surface runoff in the megadunes of Badain Jaran Desert,China, and its significance

The Badain Jaran Desert exhibits the greatest difference in altitude of all of the world’s deserts. On the slopes of megadunes in the desert, there are physical and chemical deposits produced by surface runoff. In addition, we have observed rarely-seen infiltration-excess surface runoff in the megadune depressions as well as spring streams at the base of megadunes. We used electron microscopy, energy spectrum analysis, infiltration experiments, moisture content determinations and grain-size analysis to study the mineral and chemical composition of the runoff precipitates, and grain-size of the deposits associated with the runoff, together with the hydrological balance in the megadune area, and the atmospheric precipitation mechanism responsible for groundwater recharge and for supplying water to lakes. The observations of shallow runoff and infiltration-excess surface runoff indicate the occurrence of strong and effective precipitation in summer, which would provide an important source for groundwater recharge. Several lines of evidence, such as the physical and chemical deposits resulting from shallow subsurface runoff, spring streams, infiltration-excess runoff, and gravity capillary water with a moisture content of 3–6%, demonstrate that precipitation reaches the base of the megadunes through infiltration and subsequently becomes groundwater. The chemical deposits, such as newly-formed calcite and gypsum, and gray-black physical deposits, as well as different stages in the development of fan-shaped landforms resulting from shallow subsurface runoff, indicate that groundwater recharge in the area is the result of long-term precipitation, rather than intermittent individual major rainfall events. Fine sand layers with a low infiltration capacity lead to subsurface runoff emerging at the ground surface. Five factors play an important role in maintaining a positive water balance and in replenishing groundwater via rainfall: effective rainfall as a water source, the high infiltration capacity of the sands enabling rainfall to rapidly become capillary water in the dunes, low evapotranspiration rates due to the sparse vegetation, the fact that the depth of the sand layer influenced by evaporation is shallow enough to maximize the deep infiltration of rainfall, and rapidly-moving gravity capillary water in the sandy dunes. These five factors together constitute a mechanism for groundwater recharge from rainfall, and explain the origin of the groundwater and lakes in the area. Our findings represent a significant advance in research on the hydrological cycle, including groundwater recharge conditions and recharge mechanisms, in this desert region.     

Water regime of Vrana Lake in Dalmatia (Croatia): changes,risks and problems

Abstract

Vrana Lake in Dalmatia is a karstic kryptodepression connected to the nearby sea through the karstic subsoil and a canal. Due to interactions with the sea, lake water salinity increases greatly during severe dry periods, seriously endangering the ecosystem. Trend analysis (1961–2010) reveals a decrease in precipitation and surface inflow, but an increase in air temperature, and in sea and lake water levels. Lake inflow and water losses are only partially monitored. Average annual inflow from the monitored part of the catchment is 1722 m³ s^-1, but total inflow is significantly greater; the average difference between total inflow and cumulative water losses is 3072 m³ s^-1. The paper uses modelling to evaluate total inflow into the lake system, taking into consideration projected climate changes/variations till 2100 from the RegCM3 and ALADIN climate models. The analysis indicates marked decrease in discharge values by the end of this century, by as much as 60%.

Editor Z.W. Kundzewicz     

Prediction of flow rate of karstic springs using support vector machines

Complex void space structure and flow patterns in karstic aquifers render behaviour prediction of karstic springs difficult. Four support vector regression-based models are proposed to predict flow rates from two adjacent karstic springs in Greece (Mai Vryssi and Pera Vryssi). Having no accurate estimates of the groundwater flow pattern, we used four kernels: linear, polynomial, Gaussian radial basis function and exponential radial basis function (ERBF). The data used for training and testing included daily and mean monthly precipitation, and spring flow rates. The support vector machine (SVM) performance depends on hyper-parameters, which were optimized using a grid search approach. Model performance was evaluated using root mean square error and correlation coefficient. Polynomial kernel performed better for Mai Vryssi and the ERBF for Pera Vryssi. All models except one performed better for Pera Vryssi. Our models performed better than generalized regression neural network, radial basis function neural network and ARIMA models.     

Delineation of regional arid karstic aquifers: an integrative data approach

This research integrates data procedures for the delineation of regional ground water flow systems in arid karstic basins with sparse hydrogeologic data using surface topography data, geologic mapping, permeability data, chloride concentrations of ground water and precipitation, and measured discharge data. This integrative data analysis framework can be applied to evaluate arid karstic aquifer systems globally. The accurate delineation of ground water recharge areas in developing aquifer systems with sparse hydrogeologic data is essential for their effective long-term development and management. We illustrate the use of this approach in the Cuatrociénegas Basin (CCB) of Mexico. Aquifers are characterized using geographic information systems for ground water catchment delineation, an analytical model for interbasin flow evaluation, a chloride balance approach for recharge estimation, and a water budget for mapping contributing catchments over a large region. The test study area includes the CCB of Coahuila, Mexico, a UNESCO World Biosphere Reserve containing more than 500 springs that support ground water-dependent ecosystems with more than 70 endemic organisms and irrigated agriculture. We define recharge areas that contribute local and regional ground water discharge to springs and the regional flow system. Results show that the regional aquifer system follows a topographic gradient that during past pluvial periods may have linked the Río Nazas and the Río Aguanaval of the Sierra Madre Occidental to the Río Grande via the CCB and other large, currently dry, upgradient lakes.