Temporal variation of stable isotopes in a precipitation–groundwater system: implications for determining the mechanism of groundwater recharge in high mountain–hills of the Loess Plateau,China

The groundwater in shallow loess aquifers in high mountain–hills in the western Loess Plateau in China is almost the sole water resource for local residents. However, the question about how the loess groundwater naturally circulates in these high mountain–hills, characterized by low precipitation and high potential evaporation, remains unclear. The objectives of this study are to evaluate the application of hydrogen and oxygen isotopes to (1) examine temporal variations of the isotopic composition of precipitation and shallow groundwater and (2) uncover the mechanism of groundwater recharge in high mountain–hills. Results from 2 years of monitoring data show a difference in the stable isotopes for groundwater and local precipitation between the winter and summer periods. Similar to precipitation, stable isotopes in groundwater are observed to be depleted in winter and enriched in summer, particularly in oxygen isotope. A prominent characteristic is that H and O isotopes of groundwater show a very clear response to strong precipitation in the rainy season in 2013. The results highlight that local precipitation is the likely recharge source for groundwater in shallow loess aquifers. Annual recharge from local precipitation maintains the groundwater resource in the shallower loess aquifer. The mechanisms governing shallow loess groundwater recharge in high mountain–hills were evaluated. In addition to possible vertical slow percolation of soil water through the unsaturated zone, rapid groundwater recharge mechanisms have been identified as temporal preferential infiltration through sinkholes, slip surface or landslide surface and through the interface of loess layer and palaeo‐soils. Most groundwater can be recharged after a heavy rainy season. Copyright © 2015 John Wiley & Sons, Ltd.     

Groundwater-supported evapotranspiration within glaciated watersheds under conditions of climate change

This paper analyzes the effects of geology and geomorphology on surface-water/-groundwater interactions, evapotranspiration, and recharge under conditions of long-term climatic change. Our analysis uses hydrologic data from the glaciated Crow Wing watershed in central Minnesota, USA, combined with a hydrologic model of transient coupled unsaturated/saturated flow (HYDRAT2D). Analysis of historical water-table (1970–1993) and lake-level (1924–2002) records indicates that larger amplitude and longer period fluctuations occur within the upland portions of watersheds due to the response of the aquifer system to relatively short-term climatic fluctuations. Under drought conditions, lake and water-table levels fell by as much as 2–4 m in the uplands but by 1 m in the lowlands. The same pattern can be seen on millennial time scales. Analysis of Holocene lake-core records indicates that Moody Lake, located near the outlet of the Crow Wing watershed, fell by as much as 4 m between about 4400 and 7000 yr BP. During the same time, water levels in Lake Mina, located near the upland watershed divide, fell by about 15 m. Reconstructed Holocene climate as represented by HYDRAT2D gives somewhat larger drops (6 and 24 m for Moody Lake and Lake Mina, respectively). The discrepancy is probably due to the effect of three-dimensional flow. A sensitivity analysis was also carried out to study how aquifer hydraulic conductivity and land-surface topography can influence water-table fluctuations, wetlands formation, and evapotranspiration. The models were run by recycling a wet year (1985, 87 cm annual precipitation) over a 10-year period followed by 20 years of drier and warmer climate (1976, 38 cm precipitation). Model results indicated that groundwater-supported evapotranspiration accounted for as much as 12% (10 cm) of evapotranspiration. The aquifers of highest hydraulic conductivity had the least amount of groundwater-supported evapotranspiration owing to a deep water table. Recharge was even more sensitive to aquifer hydraulic conductivity, especially in the lowland regions. These findings have important implications for paleoclimatic studies, because the hydrologic response of a surface-water body will vary across the watershed to a given climate signal.     

Distinguishing Mountain Front and Mountain Block Recharge in an Intermontane Basin of the Himalayan Region

Intermontane basin aquifers worldwide, particularly in the Himalayan region, are recharged largely by the adjoining mountains. Recharge in these basins can occur either by water infiltrating from streams near mountain fronts (MFs) as mountain front recharge (MFR) or by sub-surface mountain block infiltration as mountain block recharge (MBR). MFR and MBR recharge are challenging to distinguish and are least quantified, considering the lack of extensive understanding of the hydrological processes in the mountains. This study used oxygen and hydrogen isotopes (δ¹⁸O and δ²H), electrical conductivity (EC) data, hydraulic head, and groundwater level data to differentiate MFR and MBR. Groundwater level data provide information about the groundwater-surface water interactions and groundwater flow directions, whereas isotopes and EC data are used to distinguish and quantify different recharge sources. The present methodology is tested in an intermontane basin of the Himalayan region. The results suggest that karst springs (KS) and deep groundwater (DGW) recharge are dominated by snowmelt (47% ± 10% and 46% ± 9%) as MBR from adjacent mountains, insignificantly affected by evaporation. The hydraulic head data and isotopes indicate Quaternary shallow groundwater (SGW) aquifer system recharge as MFR of local meteoric water with significant evaporation. The results indicate several flow paths in the aquifer system, a local flow for KS, intermediate flow for SGW, and regional flow for DGW. The findings will significantly impact water resource management in the area and provide vital baseline knowledge for sustainable groundwater management in other Himalayan intermontane basins.     

Groundwater recharge from an oxbow lake-wetland system in the Mississippi Alluvial Plain

The Mississippi River Valley Alluvial Aquifer ranks among the most overdrafted aquifers in the United States due to intensive irrigation. Concern over declining water levels has increased focus on understanding the sources of recharge. Numerous oxbow lakes overlie the aquifer that are often considered hydraulically disconnected from the groundwater system due to fine-grained bottom sediments. In the current study, groundwater levels in and around a 445-ha oxbow lake-wetland in Mississippi were monitored for a 2-year period that included an unusually long low-water condition in the lake (>17 months), followed by a high-water event lasting over 4 months before returning to earlier low-water levels. The high-water pulse (>4 m rise) provided a unique opportunity to track the impact in the underlying alluvial aquifer. During low-water conditions, groundwater flowed westward beneath the lake. Following the lake rise, groundwater beneath and near the perimeter responded as quickly as the same day, with more delayed responses moving away from the lake. Within 2 months, a groundwater mound formed near the centre of the oxbow (>3 m increase), with a reversal in the local hydraulic gradient towards the east. Flow returned to a westward gradient when the lake level dropped back below 0.3 m. Analysis of precipitation and nearby river stage could not account for the observed behavior. Recharge to the aquifer is attributed to rising water levels spreading over point bar deposits and into the surrounding forested wetlands where preferential flow pathways are likely to exist due to buried and decomposing tree remains. An earlier study in the wetland demonstrated an increasing redox potential in isolated zones, consistent with the existence of preferential flow pathways through the bottom sediments (Lahiri & Davidson, 2020). Retaining high-water levels in oxbow lakes could be a relatively low-cost water management practice for enhancing aquifer recharge.     

Modelling regional effects of artificial groundwater recharge in a multilayer aquifer characterized by perched water tables

This paper presents an approach to estimate the effects of a managed recharge experiment in a multilayer aquifer characterized by the presence of perched water tables in the Medina del Campo groundwater body, Douro basin, central Spain. A numerical model was developed to evaluate the effect of artificial recharge on the shallow sector of a regional-scale aquifer and on formerly active wetlands. The model was developed in the Visual MODFLOW Pro v.2011.1 environment in order to represent and analyse the regional impact of this artificial recharge event. Results suggest that the assumption of a single perched system may prove useful in regional contexts where data is limited. From a study site perspective, managed recharge is observed to increase shallow storage along the riverbanks, which is considered valuable for environmental purposes. However, downstream wetlands are unlikely to experience a significant recovery. Furthermore, only a small percentage of artificial recharge is expected to reach the deep regional aquifer. This method can be exported to settings characterized by the presence of perched aquifers and associated groundwater dependent ecosystems.     

Origin of sulphate in the unsaturated zone and groundwater of a loess aquifer

The use of the sulphate mass balance (SMB) between precipitation and soil water as a supplementary method to estimate the diffuse recharge rate assumes that the sulphate in soil water originated entirely from atmospheric deposition; however, the origin of sulphate in soil and groundwater is often unclear, especially in loess aquifers. This study analysed the sulphur (δ³⁴S-SO₄) and oxygen (δ¹⁸O-SO₄) isotopes of sulphate in precipitation, water-extractable soil water, and shallow groundwater samples and used these data along with hydrochemical data to determine the sources of sulphate in the thick unsaturated zone and groundwater of a loess aquifer. The results suggest that sulphate in groundwater mainly originated from old precipitation. When precipitation percolates through the unsaturated zone to recharge groundwater, sulphates were rarely dissolved due to the formation of CaCO₃ film on the surface of sulphate minerals. The water-extractable sulphate in the deep unsaturated zone (>10 m) was mainly derived from the dissolution of evaporite minerals and there was no oxidation of sulphide minerals during the extraction of soil water by elutriating soil samples with deionized water. The water-extractable concentration of SO₄ was not representative of the actual SO₄ concentration in mobile soil water. Therefore, the recharge rate cannot be estimated by the SMB method using the water-extractable concentration of SO₄ in the loess areas. This study is important for identifying sulphate sources and clarifying the proper method for estimating the recharge rate in loess aquifers.     

Characterization of Flow Dynamics and Vulnerability in a Coastal Aquifer System

Traditional aquifer vulnerability techniques primarily rely on spatial property data for a region and are limited by their ability to directly or indirectly assess flow and transport processes occurring from the surface to depth within an aquifer system. The main objective of this study was to investigate groundwater vulnerability in terms of aquifer interconnectivity and flow dynamics. A combination of stable isotopes, groundwater age‐dating (radiocarbon), and geomorphic/geogenic spatial analyses was applied to a regional, highly developed coastal aquifer to explain the presence of nitrate at depth. The average δ¹³C value (?17.3 ± 2‰ VPDB, n = 27) is characteristic of groundwater originating from locally infiltrated precipitation through extensively cultivated soils. The average δ¹⁸O and δD values (?4.0 ± 0.1‰ VSMOW, n = 27; δD: ?19.3 ± 1‰ VSMOW, n = 27, respectively) are similar to precipitation water derived from maritime sources feeding the region's surface water and groundwater. Stable and radioactive isotopes reveal significant mixing between shallow and deep aquifers due to high velocities, hydraulic connection, and input of local recharge water to depths. Groundwater overdevelopment enhances deeper and faster modern water downward flux, amplifying aquifer vulnerability. Therefore, aquifer vulnerability is a variable, dependent on the type and degree of stress conditions experienced by a groundwater system as well as the geospatial properties at the near surface.     

Assessing coupling between lakes and layered aquifers in a complex Pleistocene landscape based on water level dynamics

The biosphere reserve Schorfheide-Chorin is a scenic region with many lakes. Hydraulic coupling between lakes and groundwater is difficult to assess due to the very heterogeneous Pleistocene deposits with a complex layering of different aquifers, part of them being confined. Thus, a principal component analysis of time series of groundwater and lake water levels was performed. The first two principal components provided a quantitative measure of damping of the input signal, i.e., the extent to which time series of groundwater pressure heads or lake water levels are smoothed and delayed with respect to the input signal, i.e., groundwater recharge or precipitation minus evapotranspiration, respectively. The lakes differed substantially with respect to damping behaviour, indicating different impacts of deep groundwater contribution. For most of the groundwater wells, damping increased linearly with mean depth to water table. In contrast, some wells exhibited nearly identical behaviour independent of depth. High-pass filtered data of water table level from these wells were strongly and inversely correlated with those of barometric pressure fluctuations, pointing to a confined aquifer which was evidently not connected to the adjacent lake.     

Geochemistry and Provenance of Springs in a Baja California Sur Mountain Catchment

Fractured rock aquifers cover much of Earth's surface and are important mountain sites for groundwater recharge but are poorly understood. To investigate groundwater systematics of a fractured-dominated aquifer in Baja California Sur, Mexico, we examined the spatial patterns of aquifer recharge and connectivity using the geochemistry of springs. We evaluate a range of geochemical data within the context of two endmember hypotheses describing spatial recharge patterns and fracture connectivity. Hypothesis 1 is that the aquifer system is segmented, and springs are fed by local recharge. Hypothesis 2 is that the aquifer system is well connected, with dominant recharge occurring in the higher elevations. The study site is a small <15 km² catchment. Thirty-four distinct springs and two wells were identified in the study area, and 24 of these sites were sampled for geochemical analyses along an elevation gradient and canyon transect. These analyses included major ion composition, trace element and strontium isotopes, δ¹⁸O and δ²H isotopes, radiocarbon, and tritium. δ¹⁸O and δ²H isotopes suggest that the precipitation feeding the groundwater system has at least two distinct sources. Carbon isotopes showed a change along the canyon transect, suggesting that shorter flowpaths feed springs in the top of the transect, and longer flowpaths discharge near the bottom. Geochemical interpretations support a combination of the two proposed hypotheses. Understanding of the connectivity and provenance of these springs is significant as they are the primary source of water for the communities that inhabit this region and may be impacted by changes in recharge and use.     

Hydrodynamics of floodplain wetlands in a chalk catchment: The River Lambourn, UK

Variations in floodplain channel water levels and valley floor groundwater levels (measured in piezometers and boreholes) are examined at selected points along the course of the River Lambourn, a chalk river in southern England. A local alluvial gravel aquifer in the valley bottom is associated with numerous small wetlands that extend over much of the river's perennial profile. Variations in hydraulic gradient between local borehole levels and/or floodplain channel water levels are described for three sites in the seasonal section of the channel at Bockhampton, East Garston and West Shefford. The results indicate that observed groundwater levels are closely associated with flows from discrete springs at the margins of the channel and floodplain. However, as the floodplain widens and the alluvial gravel aquifer increases in size, the gravel aquifer accounts for a substantial down-valley component of groundwater flow with a diffuse vertical water flux. In the lower catchment, the exchange of flows between the gravel aquifer and the river enables some attenuation of floodplain water-table variability, providing a stable hydrological regime for valley-bottom wetlands. Catchment controls upon the local, valley-bottom, wetland regime are demonstrated with the application of a simple groundwater model developed using MODFLOW. The model is used to simulate groundwater discharge to the river in the upper and lower catchment, in addition to the water level regime at selected points in the valley bottom in the lower catchment. The results demonstrate the importance of taking catchment-scale water flow into account when managing isolated wetlands in a permeable catchment.     

Determination of groundwater flow-through regimes of shallow lakes and wetlands from numerical analysis of stable isotope and chloride tracer distribution patterns

A numerical model previously developed to systematically examine groundwater flow in vertical section near shallow surface water bodies such as lakes, wetlands and ponds is further developed to include simulation of the distribution patterns of hydrogeochemical and stable isotopic tracers in relation to the surface water body and the geometry of distribution patterns of the tracers in the groundwater release zone of the lake. Many different possible flow regimes are identified, however, in this paper attention is focused on flow-through water bodies, as these are the flow regimes observed in field validation. Two shallow lakes on the Swan Coastal Plain of south-west Western Australia are the subject of field studies where hydrogeochemical and stable isotopes are used to validate the flow-through groundwater flow regime predicted by the modelling confirming the validity of the approach. The flow regime transition diagrams introduced in earlier papers are extended to include consideration of the hydrogeochemical and stable isotopic indices Cl_L/Cl₊ and (1000+δ_L)/(1000+δ₊). These ratios are introduced as an additional two of nine non-dimensional ratios that are necessary to analyse the problem. The ratios represent the chloride and isotopic composition (Cl_L and δ_L), respectively, in the groundwater release zones of the lakes, relative to these parameters in the groundwater capture zone (Cl₊ and δ₊) for the lake. Field data from the case studies plotted on appropriately configured transition diagrams demonstrates the overall validity of the modelling approach and its underlying assumptions. It is concluded that isotopic and hydrogeochemical data are invaluable in interpreting the interaction between lakes or wetlands and regional aquifers as it is very difficult to make physical or hydraulic measurements in the field that allow an understanding of lake–aquifer interaction. The tools and concepts developed that are summarized in the presented transition diagrams are invaluable starting points for the consideration and analysis of other case-specific examples of groundwater–surface water interaction and will improve the scientific basis of decision-making concerning lake and wetland management and groundwater interaction by water resource and environmental managers.     

Semi‐Arid Aquifer Responses to Forest Restoration Treatments and Climate Change

The purpose of this study was to develop an interpretive groundwater‐flow model to assess the impacts that planned forest restoration treatments and anticipated climate change will have on large regional, deep (>400 m), semi‐arid aquifers. Simulations were conducted to examine how tree basal area reductions impact groundwater recharge from historic conditions to 2099. Novel spatial analyses were conducted to determine areas and rates of potential increases in groundwater recharge. Changes in recharge were applied to the model by identifying zones of basal area reduction from planned forest restoration treatments and applying recharge‐change factors to these zones. Over a 10‐year period of forest restoration treatment, a 2.8% increase in recharge to one adjacent groundwater basin (the Verde Valley sub‐basin) was estimated, compared to conditions that existed from 2000 to 2005. However, this increase in recharge was assumed to quickly decline after treatment due to regrowth of vegetation and forest underbrush and their associated increased evapotranspiration. Furthermore, simulated increases in groundwater recharge were masked by decreases in water levels, stream baseflow, and groundwater storage resulting from surface water diversions and groundwater pumping. These results indicate that there is an imbalance between water supply and demand in this regional, semi‐arid aquifer. Current water management practices may not be sustainable into the far future and comprehensive action should be taken to minimize this water budget imbalance.     

Water-table rise due to infiltration from canals

Analytical solutions are presented for transient water-table profiles in finite unconfined aquifers due to infiltration occurring from a single canal or more than one canal. The aquifer boundaries are assumed to be parallel to the alignment of the canal(s) and the lower impervious bottom is horizontal. Further, the aquifer is assumed to be homogeneous and isotropic. Losses such as evapotranspiration are ignored; likewise replenishment from other sources, for example, rainfall, is also neglected.     

Recharge characteristics of an unconfined aquifer from the rainfall-water table relationship

The determination of recharge levels of unconfined aquifers, recharged entirely by rainfall, is done by developing a model for the aquifer that estimates the water-table levels from the history of rainfall observations and past water-table levels. In the present analysis, the model parameters that influence the recharge were not only assumed to be time dependent but also to have varying dependence rates for various parameters. Such a model is solved by the use of a recursive least-squares method. The variable-rate parameter variation is incorporated using a random walk model. From the field tests conducted at Tomago Sandbeds, Newcastle, Australia, it was observed that the assumption of variable rates of time dependency of recharge parameters produced better estimates of water-table levels compared to that with constant-recharge parameters. It was observed that considerable recharge due to rainfall occurred on the very same day of rainfall. The increase in water-table level was insignificant for subsequent days of rainfall. The level of recharge very much depends upon the intensity and history of rainfall. Isolated rainfalls, even of the order of 25 mm day⁻¹, had no significant effect on the water-table levels.     

Water use strategies of two co‐occurring tree species in a semi‐arid karst environment

The flow of precipitation from the surface through to groundwater in karst systems is a complex process involving storage in the unsaturated zone and diffuse and preferential recharge pathways. The processes associated with this behaviour are not well understood, despite the prevalence of karst aquifers being used as freshwater supplies. As a result, uncertainty regarding the ecohydrological processes in this geological setting remains large. In response to the need to better understand the impact of woody vegetation on groundwater recharge, annual evapotranspiration (ET) rates and tree water sources were measured for two years above a shallow, fresh karst aquifer. Water use strategies of the co‐occurring Eucalyptus diversifolia subsp. diversifolia Bonpl. and Allocasuarina verticillata (Lam.) L. Johnson were investigated using a monthly water balance approach, in conjunction with measurement of the stable isotopes of water, leaf water potentials and soil matric potentials. The results suggest that it is unlikely groundwater resources are required to sustain tree transpiration, despite its shallow proximity to the soil surface, and that similarities exist between ET losses and the estimated long‐term average rainfall for this area. Irrespective of stand and morphological differences, E. diversifolia and A. verticillata ET rates showed remarkable convergence, demonstrating the ability of these co‐occurring species to maximise their use of the available precipitation, which avoids the requirement to differentiate between these species when estimating ET at a landscape scale. We conclude that the water holding capacity of porous geological substrates, such as those associated with karst systems, will play an important role in equilibrating annual rainfall variability and should be considered when assessing ecohydrological links associated with karst systems. Copyright © 2013 John Wiley & Sons, Ltd.     

Estimation of the percentage of annual groundwater recharge with bomb tritium using a cumulative mass balance method

The bomb tritium (³H) distribution patterns in the aquifer beneath an abandoned landfill at the Canadian Forces Base (CFB) Borden, Ontario, and in a sandy aquifer at Whiteshell Nuclear Research Establishments (WNRE) Pinawa, Manitoba, all in Canada, were delineated in great detail. A sampling and monitoring network of multilevel samplers and bundle piezometers were used. The directions of groundwater flow were established, and the boundary between the tritiated and non-tritiated zones of the two aquifers were closely demarcated. Using a cumulative mass balance method, the³H input mass into the aquifers was compared with the³H mass in groundwater storage to estimate the percentages of annual groundwater recharge from 1953 to 1978. Two recharge calculations for theeffective recharge zone and thetotal recharge area of the aquifers as established from the flownet analysis, and the distributions of dissolved geochemical constitutents showed that theeffective recharge zone calculations gave higher values of 30.6 cm/yr for CFB Borden and 20.1 cm/yr for WNRE while thetotal recharge areas gave lower values of 19.1 and 10.1 cm/yr for the Borden and WNRE aquifers respectively. The two recharge values provide possible minimum and maximum recharge estimates for the two study areas.     

Effects of variations of river stage and hydraulic conductivity on temporal scaling of groundwater levels: numerical simulations

It was found in previous studies that groundwater levels may fluctuate as a temporal fractal. In this study numerical simulations of groundwater level fluctuations in an unconfined aquifer near a river were conducted to investigate the effects of aquifer heterogeneity and river stage variations on the fractal behavior of the water levels, h(t). Groundwater recharge was taken to be a white-noise process. The aquifer heterogeneity was simulated with a second-order stationary field of hydraulic conductivity (K) with an exponential variogram model. The results showed that groundwater levels fluctuate as a temporal fractal in both homogeneous and heterogeneous aquifers as long as K is less than 10 m/d. Most aquifers may indeed act as a fractal filter which takes a random non-fractal recharge inputs and produces a fractal responses of groundwater level fluctuations. A crossover in temporal scaling of h(t) may appear in more permeable aquifers. Fluctuations of the groundwater level in a homogeneous aquifer are dominated by the recharge process when the river stage is constant or by the river stage variations when the river stage varies in highly permeable aquifers. Heterogeneity plays an important role in the temporal scaling of h(t) in more permeable aquifers: the stronger the heterogeneity, the stronger the temporal scaling of h(t).     

Influence of Seasonal Variations in Sea Level on the Salinity Regime of a Coastal Groundwater–Fed Wetland

Seasonal variations in sea level are often neglected in studies of coastal aquifers; however, they may have important controls on processes such as submarine groundwater discharge, sea water intrusion, and groundwater discharge to coastal springs and wetlands. We investigated seasonal variations in salinity in a groundwater‐fed coastal wetland (the RAMSAR listed Piccaninnie Ponds in South Australia) and found that salinity peaked during winter, coincident with seasonal sea level peaks. Closer examination of salinity variations revealed a relationship between changes in sea level and changes in salinity, indicating that sea level–driven movement of the fresh water‐sea water interface influences the salinity of discharging groundwater in the wetland. Moreover, the seasonal control of sea level on wetland salinity seems to override the influence of seasonal recharge. A two‐dimensional variable density model helped validate this conceptual model of coastal groundwater discharge by showing that fluctuations in groundwater salinity in a coastal aquifer can be driven by a seasonal coastal boundary condition in spite of seasonal recharge/discharge dynamics. Because seasonal variations in sea level and coastal wetlands are ubiquitous throughout the world, these findings have important implications for monitoring and management of coastal groundwater–dependent ecosystems.     

The Effect of Modeled Recharge Distribution on Simulated Groundwater Availability and Capture

Simulating groundwater flow in basin‐fill aquifers of the semiarid southwestern United States commonly requires decisions about how to distribute aquifer recharge. Precipitation can recharge basin‐fill aquifers by direct infiltration and transport through faults and fractures in the high‐elevation areas, by flowing overland through high‐elevation areas to infiltrate at basin‐fill margins along mountain fronts, by flowing overland to infiltrate along ephemeral channels that often traverse basins in the area, or by some combination of these processes. The importance of accurately simulating recharge distributions is a current topic of discussion among hydrologists and water managers in the region, but no comparative study has been performed to analyze the effects of different recharge distributions on groundwater simulations. This study investigates the importance of the distribution of aquifer recharge in simulating regional groundwater flow in basin‐fill aquifers by calibrating a groundwater‐flow model to four different recharge distributions, all with the same total amount of recharge. Similarities are seen in results from steady‐state models for optimized hydraulic conductivity values, fit of simulated to observed hydraulic heads, and composite scaled sensitivities of conductivity parameter zones. Transient simulations with hypothetical storage properties and pumping rates produce similar capture rates and storage change results, but differences are noted in the rate of drawdown at some well locations owing to the differences in optimized hydraulic conductivity. Depending on whether the purpose of the groundwater model is to simulate changes in groundwater levels or changes in storage and capture, the distribution of aquifer recharge may or may not be of primary importance.     

Application of an integrated surface water‐groundwater model to multi‐aquifers modeling in Choushui River alluvial fan,Taiwan

This research proposes a combination of SWAT and MODFLOW, MD‐SWAT‐MODFLOW, to address the multi‐aquifers condition in Choushui River alluvial fan, Taiwan. The natural recharge and unidentified pumping/recharge are separately estimated. The model identifies the monthly pumping/recharge rates in multi‐aquifers so that the daily streamflow can be simulated correctly. A multi‐aquifers condition means a subsurface formation composed of at least the unconfined aquifer, the confined aquifer, and an in‐between aquitard. In such a case, the variation of groundwater level is related to pumping/recharge activities in vertically adjacent aquifer and the river‐aquifer interaction. Both factors in turn affect the streamflow performance. Results show that MD‐SWAT‐MODFLOW performs better than SWAT alone in terms of simulated streamflow, especially during low flow period, when pumping/recharge rates are properly estimated. A sensitivity analysis of individual parameter suggests that the vertical leakance may be the most sensitive among all investigated MODFLOW parameters in terms of the estimated pumping/recharge among aquifers, and the Latin‐Hypercube‐One‐factor‐At‐a‐Time sensitivity analysis indicates that the hydraulic conductivity of channel is the most sensitive to the model performance. It also points out the necessity to simultaneously estimate pumping/recharge rates in multi‐aquifers. The estimated net pumping rate can be treated as a lower bound of the actual local pumping rate. As a whole, the model provides the spatio‐temporal groundwater use, which gives the authorities insights to manage groundwater resources. Copyright © 2012 John Wiley & Sons, Ltd.