Dissolution in a variably confined carbonate platform: effects of allogenic runoff,hydraulic damming of groundwater inputs,and surface–groundwater exchange at the basin scale

In variably confined carbonate platforms, impermeable confining units collect rainfall over large areas and deliver runoff to rivers or conduits in unconfined portions of platforms. Runoff can increase river stage or conduit heads in unconfined portions of platforms faster than local infiltration of rainfall can increase groundwater heads, causing hydraulic gradients between rivers, conduits and the aquifer to reverse. Gradient reversals cause flood waters to flow from rivers and conduits into the aquifer where they can dissolve limestone. Previous work on impacts of gradient reversals on dissolution has primarily emphasized individual caves and little research has been conducted at basin scales. To address this gap in knowledge, we used legacy data to assess how a gradient of aquifer confinement across the Suwannee River Basin, north‐central Florida affected locations, magnitudes and processes of dissolution during 2005–2007, a period with extreme ranges of discharge. During intense rain events, runoff from the confining unit increased river stage above groundwater heads in unconfined portions of the platform, hydraulically damming inputs of groundwater along a 200 km reach of river. Hydraulic damming allowed allogenic runoff with SI_CAL < ?4 to fill the entire river channel and flow into the aquifer via reversing springs. Storage of runoff in the aquifer decreased peak river discharges downstream and contributed to dissolution within the aquifer. Temporary storage of allogenic runoff in karst aquifers represents hyporheic exchange at a scale that is larger than found in streams flowing over non‐karst aquifers because conduits in karst aquifers extend the area available for exchange beyond river beds deep into aquifers. Post‐depositional porosity in variably confined carbonate platforms should thus be enhanced along rivers that originate on confining units. This distribution should be considered in models of porosity distribution used to manage water and hydrocarbon resources in carbonate rocks. Copyright © 2013 John Wiley & Sons, Ltd.     

Airborne Thermal Remote Sensing for Estimation of Groundwater Discharge to a River

Traditional methods for studying surface water and groundwater interactions have usually been limited to point measurements, such as geochemical sampling and seepage measurement. A new methodology is presented for quantifying groundwater discharge to a river, by using river surface temperature data obtained from airborne thermal infrared remote sensing technology. The Hot Spot Analysis toolkit in ArcGIS was used to calculate the percentage of groundwater discharge to a river relative to the total flow of the river. This methodology was evaluated in the midstream of the Heihe River in the arid and semiarid northwest China. The results show that the percentage of groundwater discharge relative to the total streamflow was as high as 28%, which is in good agreement with the results from previous geochemical studies. The data analysis methodology used in this study is based on the assumption that the river water is fully mixed except in the areas of extremely low flow velocity, which could lead to underestimation of the amount of groundwater discharge. Despite this limitation, this remote sensing‐based approach provides an efficient means of quantifying the surface water and groundwater interactions on a regional scale.     

Quantifying River‐Groundwater Interactions of New Zealand's Gravel‐Bed Rivers: The Wairau Plain

New Zealand's gravel‐bed rivers have deposited coarse, highly conductive gravel aquifers that are predominantly fed by river water. Managing their groundwater resources is challenging because the recharge mechanisms in these rivers are poorly understood and recharge rates are difficult to predict, particularly under a more variable future climate. To understand the river‐groundwater exchange processes in gravel‐bed rivers, we investigate the Wairau Plain Aquifer using a three‐dimensional groundwater flow model which was calibrated using targeted field observations, “soft” information from experts of the local water authority, parameter regularization techniques, and the model‐independent parameter estimation software PEST. The uncertainty of simulated river‐aquifer exchange flows, groundwater heads, spring flows, and mean transit times were evaluated using Null‐space Monte‐Carlo methods. Our analysis suggests that the river is hydraulically perched (losing) above the regional water table in its upper reaches and is gaining downstream where marine sediments overlay unconfined gravels. River recharge rates are on average 7.3 m³/s, but are highly dynamic in time and variable in space. Although the river discharge regularly hits 1000 m³/s, the net exchange flow rarely exceeds 12 m³/s and seems to be limited by the physical constraints of unit‐gradient flux under disconnected rivers. An important finding for the management of the aquifer is that changes in aquifer storage are mainly affected by the frequency and duration of low‐flow periods in the river. We hypothesize that the new insights into the river‐groundwater exchange mechanisms of the presented case study are transferable to other rivers with similar characteristics.     

Formation waters discharge to rivers near oil sands projects

Groundwater discharges in the western Canadian oil sands region impact river water quality. Mapping groundwater discharges to rivers in the oil sands region is important to target water quality monitoring efforts and to ensure injected wastewater and steam remain sequestered rather than eventually resurfacing. Saline springs composed of Pleistocene‐aged glacial meltwater exist in the region, but their spatial distribution has not been mapped comprehensively. Here we show that formation waters discharge into 3 major rivers as they flow through the Athabasca Oil Sands Region adjacent to many active oil sands projects. These discharges increase river chloride concentrations from river headwaters to downstream reaches by factors of ~23 in the Christina River, ~4 in the Clearwater River, and ~5 in the Athabasca River. Our survey provides further evidence for the substantial impact of formation water discharges on river water quality, even though they comprise less than ~2% of total streamflow. Geochemical evidence supporting formation water discharges as the leading control on river salinity include increases in river chloride concentrations, Na/(Na + Ca) ratios, Cl/(Cl + SO₄) ratios and decreases in ⁸⁷Sr/⁸⁶Sr ratios; each mixing trend is consistent with saline groundwater discharges sourced from Cretaceous or Devonian aquifers. These regional subsurface‐to‐surface connections signify that injected wastewater or steam may potentially resurface in the future, emphasizing the critical importance of mapping groundwater flow paths to understand present‐day streamflow quality and to predict the potential for injected fluids to resurface.     

Formation of phreatic caves in an eogenetic karst aquifer by CO2 enrichment at lower water tables and subsequent flooding by sea level rise

Formation of extensive phreatic caves in eogenetic karst aquifers is widely believed to require mixing of fresh and saltwater. Extensive phreatic caves also occur, however, in eogenetic karst aquifers where fresh and saltwater do not mix, for example in the upper Floridan aquifer. These caves are thought to have formed in their modern settings by dissolution from sinking streams or by convergence of groundwater flow paths on springs. Alternatively, these caves have been hypothesized to have formed at lower water tables during sea level low‐stands. These hypotheses have not previously been tested against one another. Analyzing morphological data and water chemistry from caves in the Suwannee River Basin in north‐central Florida and water chemistry from wells in the central Florida carbonate platform indicates that phreatic caves within the Suwannee River Basin most likely formed at lower water tables during lower sea levels. Consideration of the hydrological and geochemical constraints posed by the upper Floridan aquifer leads to the conclusion that cave formation was most likely driven by dissolution of vadose CO₂ gas into the groundwater. Sea level rise and a wetter climate during the mid‐Holocene lifted the water table above the elevation of the caves and placed the caves tens of meters below the modern water table. When rising water tables reached the land surface, surface streams formed. Incision of surface streams breached the pre‐existing caves to form modern springs, which provide access to the phreatic caves. Phreatic caves in the Suwannee River Basin are thus relict and have no causal relationship with modern surficial drainage systems. Neither mixing dissolution nor sinking streams are necessary to form laterally extensive phreatic caves in eogenetic karst aquifers. Dissolution at water tables, potentially driven by vadose CO₂ gas, offers an underappreciated mechanism to form cavernous porosity in eogenetic carbonate rocks. Copyright © 2012 John Wiley & Sons, Ltd.     

Evaluation of riverbed magnetic susceptibility for mapping biogeochemical hot spots in groundwater-impacted rivers

Redox hot spots occurring as metal-rich anoxic groundwater discharges through oxic wetland and river sediments commonly result in the formation of iron (Fe) oxide precipitates. These redox-sensitive precipitates influence the release of nutrients and metals to surface water and can act as ‘contaminant sponges’ by absorbing toxic compounds. We explore the feasibility of a non-invasive, high-resolution magnetic susceptibility (MS) technique to efficiently map the spatial variations of magnetic Fe oxide precipitates in the shallow bed of three rivers impacted by anoxic groundwater discharge. Laboratory analyses on Mashpee River (MA, USA) sediments demonstrate the sensitivity of MS to sediment Fe concentrations. Field surveys in the Mashpee and Quashnet rivers (MA, USA) reveal several discrete high MS zones, which are associated with likely anoxic groundwater discharge as evaluated by riverbed temperature, vertical head gradient, and groundwater chemistry measurements. In the East River (CO, USA), widespread cobbles/rocks exhibit high background MS from geological ferrimagnetic minerals, thereby obscuring the relatively small enhancement of MS from groundwater induced Fe oxide precipitates. Our study suggests that, in settings with low geological sources of magnetic minerals such as lowland rivers and wetlands, MS may serve as a complementary tool to temperature methods for efficiently mapping Fe oxide accumulation zones due to anoxic groundwater discharges that may function as biogeochemical hot spots and water quality control points in gaining systems.     

Incubation,hatching and survival of eggs of Atlantic salmon (Salmo salar) in spawning redds influenced by groundwater

Many west coastal and northern Norwegian rivers run through deep, confined valleys with permeable layers of glacial and alluvial deposits. Groundwater flows through these permeable layers and enter lakes and rivers as underwater seepage and springs. Groundwater inflow to inland Norwegian rivers may constitute 40–100% of total water discharge during low flow periods in late summer and winter. Juvenile salmonids may take advantage of groundwater upwellings and actively seek out such patches. In regulated rivers groundwater influx may create refuges during low flow or hydropeaking episodes. The importance of groundwater for salmon redd site selection and egg survival is also clear, although less known and documented in regulated rivers.Eggs of Atlantic salmon (Salmo salar) are deposited in redds in river bed gravels lacking fine sediments and with high oxygen levels. Egg development is therefore dependent on the interaction of a number of environmental factors such as groundwater influx, oxygen and temperature. Atlantic salmon in the regulated River Suldalslågen, Western Norway, spawn relatively late compared to other Norwegian rivers, with a peak in early January. Newly emerged fry are found from the end of May to the beginning of June, i.e. “swim up” one month earlier than expected using models for egg and alevin development and river water temperatures. The most plausible explanation is that groundwater has a higher and more stable temperature than surface river water. In field experiments, fertilized salmon eggs were placed in boxes close to natural spawning redds in the river bed at sites influenced and those not influenced by groundwater. A difference of up to 40 days in 50% hatching was found, and “swim up” occurred at the end of May in boxes influenced by groundwater.Preliminary studies have revealed that groundwater also plays an important role in survival of salmon eggs in the River Suldalslågen when dewatered in winter. Eggs placed in boxes in groundwater seepage areas during winter in the dewatered river bed survived even when covered by ice and snow. The survival from fertilization until 30 April, one month before hatching, was 91%, the same survival as found for eggs placed in boxes in the wetted river bed. However, mortality from fertilization to hatching was higher compared to the eggs placed in wetted river bed, 57 and 91% respectively.Groundwater creates a horizontal and vertical mosaic of temperatures in spawning redd areas leading to potentially greater variation in spawning sites, time of hatching and “swim up”. This is likely to increase egg survival during low flow periods in regulated rivers. In conclusion, the interaction between groundwater and surface river water should therefore be considered when managing fish populations in regulated rivers.     

Exploring landscape and geologic controls on spatial patterning of streambank groundwater discharge in a mixed land use watershed

Preferential groundwater discharge features along stream corridors are ecologically important at local and stream network scales, yet we lack quantification of the multiscale controls on the spatial patterning of groundwater discharge. Here we identify physical attributes that best explain variation in the presence and lateral extent of preferential groundwater discharges along two 5th order streams, the Housatonic and Farmington Rivers, and 32 1st to 4th order reaches across the Farmington River network. We mapped locations of preferential groundwater discharge exposed along streambanks using handheld thermal infrared cameras paired with high-resolution topographic and land use land cover datasets, surficial soil characteristic maps, and depth-to-bedrock geophysical measurements. The unconfined Housatonic River, MA, USA (12 km) had fewer discharge locations and less lateral extent (41 discharge locations with 38 m of active discharge/km of river) compared to the partially confined Farmington River, CT, USA (26 km; 169 discharge locations with 129 m of active discharge/km of river). Using a moving window analysis, we found along both rivers that discharge was more likely to occur where bank slopes were steeper, floodplain extent was narrower, and degree of confinement was higher. Along the Farmington River, groundwater discharge was more likely to occur where saturated hydraulic conductivity was higher and depth-to-bedrock was shallower. Among the 32 stream reaches surveyed (33.2 km of total stream length) within the Farmington River watershed, preferential discharge was observed in all but two stream reaches, varied from 0 to 25% of lateral extent along stream banks (mean = 6%), and was more likely to occur where stream reach slopes were steep, saturated hydraulic conductivity was high, and watershed urbanization was low. Our results show that, though both surface (e.g., topographic, land use land cover) and subsurface (e.g., soil characteristics, bedrock depth) factors control the prevalence of streambank preferential groundwater discharge, the dominant controls vary across valley settings and stream sizes.     

Combining isotopic tracers (222Rn and δ13C) for improved modelling of groundwater discharge to small rivers

In regions where aquifers sustain rivers, the location and quantification of groundwater discharge to surface water are important to prevent pollution hazards, to quantify and predict low flows and to manage water supplies. ²²²Rn is commonly used to determine groundwater discharge to rivers. However, using this isotopic tracer is challenging because of the high diffusion capacity of ²²²Rn in open water. This study illustrates how a combination of isotopic tracers can contribute to an enhanced understanding of groundwater discharge patterns in small rivers. The aim of this paper is to combine ²²²Rn and δ¹³C_DIC to better constrain the physical parameters related to the degassing process of these tracers in rivers. The Hallue River (northern France) was targeted for this study because it is sustained almost exclusively by a fractured chalk aquifer. The isotopes ²²²Rn, δ¹³C_DIC, δ²H and δ¹⁸O were analysed along with other natural geochemical tracers. A mass balance model was used to simulate ²²²Rn and δ¹³C_DIC. The results of δ²H and δ¹⁸O analyses prove that evaporation did not occur in the river. The calibration of a numerical model to reproduce ²²²Rn and δ¹³C_DIC provides a best‐fit diffusive layer thickness of 3.21 × 10^?5 m. This approach is particularly useful for small rivers flowing over carbonate aquifers with high groundwater DIC where the evolution of river DIC reflects the competing processes of groundwater inflow and CO₂ degassing. This approach provides a means to evaluate groundwater discharge in small ungauged rivers. Copyright © 2014 John Wiley & Sons, Ltd.     

Studies of a regulated dryland river: surface–groundwater interactions

Unlike rivers in humid regions, dryland rivers typically exhibit reduced flow in the downstream direction as a result of transmission losses, which include seepage of streamflow into the aquifer, evaporation, and transpiration. However, much remains to be learned about the nature of the exchange between surface water and groundwater in these landscapes, especially in terms of spatial and temporal variability. Our study focused on streambank seepage and groundwater flow in the alluvial aquifer, specifically on answering questions such as: Is there seasonal variability in seepage losses? Is seepage permanently lost? Can losses be reduced by killing riparian vegetation? To better understand the magnitude, variability, and fate of streambank seepage, we assessed river stages, groundwater hydraulic gradients, and groundwater flow paths at two sites along a reach of the Pecos River, a dryland perennial river in West Texas. We found that along this reach the river was losing water to the aquifer even under low‐flow conditions; but seepage was controlled by a number of different mechanisms. Seepage increased not only during high‐flow events but also when the groundwater level was declining owing to long periods of no irrigation release. Tamarix (saltcedar) control did not affect hydraulic gradients nor reduce streambank seepage and given that this reach of the Pecos River is a losing one, streamflow will not be enhanced by controlling saltcedar. These findings can be used to improve basic conceptual models of dryland river systems and to predict hydrologic responses to changes in the timing and magnitude of streamflows and to riparian vegetation management. Copyright © 2012 John Wiley & Sons, Ltd.     

Three-dimensional semi-analytical solution to groundwater flow in confined and unconfined wedge-shaped aquifers

The Laplace domain solutions have been obtained for three-dimensional groundwater flow to a well in confined and unconfined wedge-shaped aquifers. The solutions take into account partial penetration effects, instantaneous drainage or delayed yield, vertical anisotropy and the water table boundary condition. As a basis, the Laplace domain solutions for drawdown created by a point source in uniform, anisotropic confined and unconfined wedge-shaped aquifers are first derived. Then, by the principle of superposition the point source solutions are extended to the cases of partially and fully penetrating wells. Unlike the previous solution for the confined aquifer that contains improper integrals arising from the Hankel transform [Yeh HD, Chang YC. New analytical solutions for groundwater flow in wedge-shaped aquifers with various topographic boundary conditions. Adv Water Resour 2006;26:471–80], numerical evaluation of our solution is relatively easy using well known numerical Laplace inversion methods. The effects of wedge angle, pumping well location and observation point location on drawdown and the effects of partial penetration, screen location and delay index on the wedge boundary hydraulic gradient in unconfined aquifers have also been investigated. The results are presented in the form of dimensionless drawdown-time and boundary gradient-time type curves. The curves are useful for parameter identification, calculation of stream depletion rates and the assessment of water budgets in river basins.     

Hydrogeochemical Evolution Along Groundwater Flow Paths in the Manas River Basin,Northwest China

The impacts of long-term pumping on groundwater chemistry remain unclear in the Manas River Basin, Northwest China. In this study, major ions within five surface water and 105 groundwater samples were analyzed to identify hydrogeochemical processes affecting groundwater composition and evolution along the regional-scale groundwater flow paths using the multivariate techniques of hierarchical cluster analysis (HCA) and principal components analysis (PCA) and traditional graphical methods for analyzing groundwater geochemistry. HCA classified the groundwater samples into four clusters (C1 to C4). PCA reduced the dimensionality of geochemical data into three PCs, which explained 86% of the total variance. The results of HCA and PCA were used to identify three zones: “recharge,” “transition,” and “discharge.” In the recharge zone the groundwater type is Ca-HCO₃-SO₄ and is primarily impacted by the dissolution of calcite and silicate weathering. In the transition zone the groundwater type is Ca-HCO₃-SO₄-Cl and is impacted by rock dissolution and reverse ion exchange. In the discharge zone the groundwater type is Na-Cl and is impacted by evaporation and reverse ion exchange. In addition, anthropogenic activities impact the groundwater chemistry in the study area. The groundwater type generally changes from Ca-HCO₃-SO₄ in the recharge area to Na-Cl in the discharge area along the regional-scale groundwater flow paths. This study provides a process-based knowledge for understanding the interaction of groundwater flow patterns and geochemical evolution within the Manas River Basin.     

Analytical investigation of the exact groundwater divide between rivers beyond the Dupuit–Forchheimer approximation

The groundwater divide is a key feature of river basins and significantly influenced by subsurface hydrological processes. For an unconfined aquifer between two parallel rivers or ditches, it has long been defined as the top of the water table based on the Dupuit–Forchheimer approximation. However, the exact groundwater divide is subject to the interface between two local flow systems transporting groundwater to rivers from the infiltration recharge. This study contributes a new analytical model for two-dimensional groundwater flow between rivers of different water levels. The flownet is delineated in the model to identify groundwater flow systems and the exact groundwater divide. Formulas with two dimensionless parameters are derived to determine the distributed hydraulic head, the top of the water table and the groundwater divide. The locations of the groundwater divide and the top of the water table are not the same. The distance between them in horizontal can reach up to 8.9% of the distance between rivers. Numerical verifications indicate that simplifications in the analytical model do not significantly cause misestimates in the location of the groundwater divide. In contrast, the Dupuit–Forchheimer approximation yields an incorrect water table shape. The new analytical model is applied to investigate groundwater divides in the Loess Plateau, China, with a Monte Carlo simulation process taking into account the uncertainties in the parameters.     

Riverine groundwater discharge estimation in a dynamic river corridor using 222Rn

Groundwater discharge flux into rivers (riverine groundwater discharge or RGD) is essential information for the conservation and management of aquatic ecosystems and resources. One way to estimate area-integrated groundwater discharge into surface water bodies is to measure the concentration of a groundwater tracer within the water body. We assessed groundwater discharge using ²²²Rn, a tracer common in many surface water studies, through field measurements, surface water ²²²Rn mass balance model, and groundwater flow simulation, for the seldom studied but ubiquitous setting of a flooding river corridor. The investigation was conducted at the dam-regulated Lower Colorado River (LCR) in Austin, Texas, USA. We found that ²²²Rn in both the river water and groundwater in the river bank changed synchronously over a 12-hour flood cycle. A ²²²Rn mass balance model allowed for estimation of groundwater discharge into a 500-m long reach of the LCR over the flood. The groundwater discharge ranged between negative values (indicating recharge) to 1570 m³/h; groundwater discharge from groundwater flow simulations corroborated these estimates. However, for the dynamic groundwater discharge estimated by the ²²²Rn box model, assuming whether the groundwater ²²²Rn endmember was constant or dynamic led to notably different results. The resultant groundwater discharge estimates are also highly sensitive to river ²²²Rn values. We thus recommend that when using this approach to accurately characterize dynamic groundwater discharge, the ²²²Rn in near-stream groundwater should be monitored at the same frequency as river ²²²Rn. If this is not possible, the ²²²Rn method can still provide reasonable but approximate groundwater discharge given background information on surface water-groundwater exchange time scales.     

Tracking Groundwater Discharge to a Large River using Tracers and Geophysics

Few studies have investigated large reaches of rivers in which multiple sources of groundwater are responsible for maintaining baseflow. This paper builds upon previous work undertaken along the Fitzroy River, one of the largest perennial river systems in north‐western Australia. Synoptic regional‐scale sampling of both river water and groundwater for a suite of environmental tracers (⁴He, ⁸⁷Sr/⁸⁶Sr, ²²²Rn and major ions), and subsequent modeling of tracer behavior in the river, has enabled definition and quantification of groundwater input from at least three different sources. We show unambiguous evidence of both shallow “local” groundwater, possibly recharged to alluvial aquifers beneath the adjacent floodplain during recent high‐flow events, and old “regional” groundwater introduced via artesian flow from deep confined aquifers. We also invoke hyporheic exchange and either bank return flow or parafluvial flow to account for background ²²²Rn activities and anomalous chloride trends along river reaches where there is no evidence of the local or regional groundwater inputs. Vertical conductivity sections acquired through an airborne electromagnetic (AEM) survey provide insights to the architecture of the aquifers associated with these sources and general groundwater quality characteristics. These data indicate fresh groundwater from about 300 m below ground preferentially discharging to the river, at locations consistent with those inferred from tracer data. The results demonstrate how sampling rivers for multiple environmental tracers of different types—including stable and radioactive isotopes, dissolved gases and major ions—can significantly improve conceptualization of groundwater—surface water interaction processes, particularly when coupled with geophysical techniques in complex hydrogeological settings.     

Dams in the last large free-flowing rivers of Patagonia,the Santa Cruz River,environmental features,and macroinvertebrate community

Three large rivers have their headwaters in the Patagonian Ice Fields (PIFs) in the Andes Mountains, the largest mid-latitude ice masses on Earth: Santa Cruz, Baker and Pascua. They are the last large free flowing rivers in Patagonia, but plans are advanced for building dams for hydroelectric power generation. The three PIF rivers, with a discharge dominated by ice melt, share a common, unique hydrograph compared to that of the other eight large rivers in the region: a distinct seasonal cycle, and an extremely stable discharge, with much lower variability than other rivers. In this study we present the first extensive survey of habitats and benthic macroinvertebrates in the least studied system, the Santa Cruz River. We assess how much of the natural capital provided and sustained by benthic invertebrates are expected to be lost by flooding and discuss how dams would affect riverine habitat and biota. In the Santa Cruz River, we conducted an intensive field survey during September 2010; a total of 52 sites located at regular 6 km intervals were sampled along the 310 river-km for macroinvertebrates and seventeen habitat variables. Although some habitat structure is apparent at the local scale, the Santa Cruz River could be described as very homogeneous. Macroinvertebrate density and the richness (38 genera) found in the Santa Cruz River resulted to be one of the lowest in comparison with 42 other Patagonian rivers. Albeit weak, the structure of the macroinvertebrates assemblages was successfully described by a reduced set of variables. The reduced flow variation and the lack of bed scouring flows have a direct and negative effect on the heterogeneity of riverbeds and banks. The high turbidity of the Santa Cruz River may also contribute to shorter food webs, by affecting autotrophic production, general trophic structure, and overall macroinvertebrate productivity and diversity. Dams will obliterate 51% of the lotic environment, including the most productive sections of the river according to our macroinvertebrate data. Since Santa Cruz River has a naturally homogeneous flow cycle, dams may provide more variable flows and more diverse habitat. Our data provide critically valuable baseline information to understand the effects of dams on the unique set of glacial driven large rivers of Patagonia.     

Analysis of steady ground water flow toward wells in a confined-unconfined aquifer

A confined aquifer may become unconfined near the pumping wells when the water level falls below the confining unit in the case where the pumping rate is great and the excess hydraulic head over the top of the aquifer is small. Girinskii's potential function is applied to analyze the steady ground water flow induced by pumping wells with a constant-head boundary in a mixed confined-unconfined aquifer. The solution of the single-well problem is derived, and the critical radial distance at which the flow changes from confined to unconfined condition is obtained. Using image wells and the superposition method, an analytic solution is presented to study steady ground water flow induced by a group of pumping wells in an aquifer bounded by a river with constant head. A dimensionless function is introduced to determine whether a water table condition exists or not near the pumping wells. An example with three pumping wells is used to demonstrate the patterns of potentiometric surface and development of water table around the wells.     

Hydrogeological and hydrogeochemical control of groundwater salinity in an arid inland basin: Dunhuang Basin,northwestern China

High groundwater salinity has become a major concern in the arid alluvial plain of the Dunhuang Basin in northwestern China because it poses a significant challenge to water resource management. Isotopic and geochemical analyses were conducted on 55 water samples from springs, boreholes and surface water to identify potential sources of groundwater salinity and analyse the processes that control increasing salinity. The total dissolved solid (TDS) content in the groundwater ranged from 400 to 41 000 mg/l, and high TDS values were commonly associated with shallow water tables and flow‐through and discharge zones in unconfined aquifers. Various groundwater contributions from rainwater, agricultural irrigation, river water infiltration and lateral inflows from mountains were identified by major ions and δD and δ¹⁸O. In general, HCO₃^? and SO₄^2? were the dominant anions in groundwater with a salinity of <2500 mg/l, whereas Cl^? and SO₄^2? were the dominant anions in groundwater with a salinity of >2500 mg/l. The major ion concentrations indicated that mineral weathering, including carbonate and evaporite dissolution, primarily affected groundwater salinity in recharge areas. Evapotranspiration controlled the major ion concentration evolution and salinity distribution in the unconfined groundwaters in the flow‐through and discharge areas, although it had a limited effect on groundwater in the recharge areas and confined aquifers. Agricultural irrigation increased the water table and enhanced evapotranspiration in the oasis areas of the basin. TDS and Cl became more concentrated, but H and O isotopes were not enriched in the irrigation district, indicating that transpiration dominated the increasing salinity. For other places in the basin, as indicated by TDS, Cl, δD and δ¹⁸O characteristics, evaporation, transpiration and water–rock interactions dominated at different hydrogeological zones, depending on the plant coverage and hydrogeological conditions. Groundwater ages of ³H, and δD and δ¹⁸O compositions and distributions suggest that most of the groundwaters in Dunhuang Basin have a paleometeoric origin and experienced a long residence time. These results can contribute to groundwater management and future water allocation programmes in the Dunhuang Basin. Copyright © 2015 John Wiley & Sons, Ltd.     

Using principal component analysis in the investigation of groundwater hydrochemistry of Upper Jezireh Basin,Syria

Principal component analysis (PCA) was applied to hydrochemical and isotopic data of 34 groundwater samples. This allowed the reduction of 20 variables to four significant PCs that explain 81.9% of the total variance; F1 (47.1%) explains the groundwater mineralization, whereas F2 (17%) shows isotopic enrichment and nitrate pollution. Based on an iso-factor scores map of F1, three water zones were delineated: Zone A (F1 < ?1), with fresh groundwater from the unconfined aquifer; Zone B (1 > F1 > ?1), with moderate mineralization from the confined–unconfined aquifer boundary; and Zone C (F1 > 1), with the most mineralized hot water from the confined aquifer. The iso-factor scores map of F2 delineates positive values representing samples from the unconfined aquifer, with freshwater and nitrate contamination associated with stable isotope enrichment, whereas negative values represent samples from the confined aquifer. The results clearly demonstrate the usefulness of PCA in groundwater hydrochemistry investigations.     

Stable isotopic and geochemical identification of groundwater evolution and recharge sources in the arid Shule River Basin of Northwestern China

Stable isotopic (δD_VSMOW and δ¹⁸O_VSMOW) and geochemical signatures were employed to constrain the geochemical evolution and sources of groundwater recharge in the arid Shule River Basin, Northwestern China, where extensive groundwater extraction occurs for agricultural and domestic supply. Springs in the mountain front of the Qilian Mountains, the Yumen‐Tashi groundwater (YTG), and the Guazhou groundwater (GZG) were Ca‐HCO₃, Ca‐Mg‐HCO₃‐SO₄ and Na‐Mg‐SO₄‐Cl type waters, respectively. Total dissolved solids (TDS) and major ion (Mg²⁺, Na⁺, Ca²⁺, K⁺, SO₄^2?, Cl^? and NO₃^?) concentrations of groundwater gradually increase from the mountain front to the lower reaches of the Guazhou Basin. Geochemical evolution in groundwater was possibly due to a combination of mineral dissolution, mixing processes and evapotranspiration along groundwater flow paths. The isotopic and geochemical variations in melt water, springs, river water, YTG and GZG, together with the end‐member mixing analysis (EMMA) indicate that the springs in the mountain front mainly originate from precipitation, the infiltration of melt water and river in the upper reaches; the lateral groundwater from the mountain front and river water in the middle reaches are probably effective recharge sources for the YTG, while contribution of precipitation to YTG is extremely limited; the GZG is mainly recharged by lateral groundwater flow from the Yumen‐Tashi Basin and irrigation return flow. The general characteristics of groundwater in the Shule River Basin have been initially identified, and the results should facilitate integrated management of groundwater and surface water resources in the study area. Copyright © 2015 John Wiley & Sons, Ltd.