Water percolation in a thick unsaturated loess layer considering the ground‐atmosphere interaction

The key objective of this paper is to advance our present understanding of how surface water infiltrates in thick unsaturated loess, which is found in arid and semiarid regions of the world, considering the ground‐atmosphere interaction. In situ data for a period of 1 year in thick loess layer at a site in the Loess Plateau of China that has groundwater table at 97.5 m depth were collected for achieving this objective. Climate factors, mainly rainfall and actual evaporation, were measured. In addition, variations of soil temperature and water content at different depths in the unsaturated zone were also measured. The data were used to interpret the water percolation characteristics by dividing the thick unsaturated zone into three zones; namely, (i) surface zone, which constitutes the top 1.0 m, (ii) unsteady zone, which is from 1.0 to 7.0 m, and (iii) steady zone, which is below 7.0 m. In the surface zone, soil temperature and water content are sensitive to climate factors. There is a variation of water content associated with the cumulative influence of infiltration and evaporation in the precipitation and nonprecipitation periods, respectively. In the unsteady zone, the water content is relatively constant; however, temperature varies in different seasons. Water percolation in this zone is both in liquid and vapour phases. In the steady zone, both soil temperature and water content are constant during the entire investigation period. The percolation velocity in this zone is approximately 1.23 × 10^?8 m/s or 0.39 m/year, which suggests that it will take approximately 230.8 years for surface water to pass through the thick unsaturated zone and recharge the groundwater.     

Water balance of a tropical woodland ecosystem, Northern Australia: a combination of micro-meteorological, soil physical and groundwater chemical approaches

A combination of micro-meteorological, soil physical and groundwater chemical methods enabled the water balance of a tropical eucalypt savanna ecosystem in Northern Australia to be estimated. Heat pulse and eddy correlation were used to determine overstory and total evapotranspiration, respectively. Measurements of soil water content, matric suction and water table variations were used to determine changes in soil moisture storage throughout the year. Groundwater dating with chlorofluorocarbons was used to estimate net groundwater recharge rates, and stream gauging was used to determine surface runoff. The wet season rainfall of 1585 mm is distributed as: evapotranspiration 810 mm, surface runoff (and shallow subsurface flow) into the river 410 mm, groundwater recharge 200 mm and increase in soil store 165 mm. Of the groundwater recharge, 160 mm enters the stream as baseflow in the wet season, 20 mm enters as baseflow in the dry season, and the balance (20 mm) is distributed to and used by minor vegetation types within the catchment or discharges to the sea. In the dry season, an evapotranspiration of 300 mm comprises 135 mm rainfall and 165 mm from the soil store. Because of the inherent errors of the different techniques, the water balance surplus (estimated at 20 mm) cannot be clearly distinguished from zero. It may also be as much as 140 mm. To our knowledge, this is the first time that such diverse methods have been combined to estimate all components of a catchment's water balance.     

Seasonal soil water dynamics in the jarrah forest,Western Australia. I: Results from a hillslope transect with coarse-textured soil profiles

Seasonal soil water dynamics were measured on a hillslope transect in the jarrah forest of southwest Western Australia over the period 1984-86 using mercury manometer tensiometers, gypsum blocks, and a neutron moisture meter. The soil water potential gradients indicated downward vertical drainage flux through winter and spring. There was generally a change to an upwards flux in early summer which was sustained through to autumn. A shallow ephemeral saturation zone was identified in and above a duricrust layer, lasting up to three days after heavy, late winter rainfall. The annual maximum to minimum unsaturated soil water storage on the hillslope was approximately 400 mm to 6 m depth and 480 mm to 15 m depth. This did not change significantly in years of substantially different winter rainfall. The magnitude of seasonal soil water storage was similar to other forested areas with deep soil profiles. The depth of observable infiltration was dependent on annual rainfall. This was consistent with the observation that groundwater levels responded to rainfall over the whole hillslope in wet years but only responded on the lower slopes in dry years. The average summer drying rate of the soil profile to 6 m depth of 3.5 mm day^?1 was within the range of values reported for forests elsewhere. In late summer, following an extended drought period, the drying rate decreased downslope but increased midslope.     

Modelling the effects of soil type and root distribution on shallow groundwater resources

Vegetated, shallow groundwater environments typically have high environmental and economic value. A sound understanding of the complex interactions and feedbacks between surface vegetation and groundwater resources is crucial to managing and maintaining healthy ecosystems while responding to human needs. A vegetated shallow groundwater environment was modelled using the software HYDRUS 2D to investigate the effects of several combinations of soil type and root distributions on shallow groundwater resources. Three rainfall regimes coupled to both natural and anthropogenically affected groundwater conditions were used to investigate the effect that combinations of four soil types and five root distributions can have on (a) groundwater level drops, (b) groundwater depletion, (c) groundwater recharge and (d) water stress conditions. Vegetation with roots distributed across the whole unsaturated zone and vegetation with dimorphic root systems (i.e. roots having larger concentrations both near the surface and the capillary fringe) behaved differently from vegetation growing roots mainly near the saturated zone. Specifically, vegetation with roots in the unsaturated zone caused water‐table drops and groundwater depletions that were half the amount due to deep‐rooted vegetation. Vegetation with a large portion of roots near the soil surface benefited from rainfall and was less vulnerable to water‐table lowering; as such, the fraction of the total area of roots affected by water stress conditions could be 40% smaller than in the case with deep‐rooted vegetation. However, roots uniformly distributed in the unsaturated zone could halve groundwater recharge rates observed in bare soils. Our analysis provided insights that can enable the formulation of site‐ and purpose‐specific management plans to respond to both human and ecosystem water requirements. Copyright © 2015 John Wiley & Sons, Ltd.     

Climate,vegetation, and soil controls on hydraulic redistribution in shallow tree roots

Hydraulic redistribution defined as the translocation of soil moisture by plant root systems in response to water potential gradients is a phenomenon widely documented in different climate, vegetation, and soil conditions. Past research has largely focused on hydraulic redistribution in deep tree roots with access to groundwater and/or winter rainfall, while the case of relatively shallow (i.e., ≈1–2 m deep) tree roots has remained poorly investigated. In fact, it is not clear how hydraulic redistribution in shallow root zones is affected by climate, vegetation, and soil properties. In this study, we developed a model to investigate the climate, vegetation, and soil controls on the net direction and magnitude of hydraulic redistribution in shallow tree root systems at the growing season to yearly timescale. We used the model to evaluate the effect of hydraulic redistribution on the water stress of trees and grasses. We found that hydraulic lift increases with decreasing rainfall frequency, depth of the rooting zone, root density in the deep soil and tree leaf area index; at the same time for a given rainfall frequency, hydraulic lift increases with increasing average rainstorm depth and soil hydraulic conductivity. We propose that water drainage into deeper soil layers can lead to the emergence of vertical water potential gradients sufficient to explain the occurrence of hydraulic lift in shallow tree roots without invoking the presence of a shallow water table or winter precipitation. We also found that hydraulic descent reduces the water stress of trees and hydraulic lift reduces the water stress of grass with important implications on tree–grass interactions.     

A conceptual hydrologic model for a forested Carolina bay depressional wetland on the Coastal Plain of South Carolina,USA

This paper describes how climate influences the hydrology of an ephemeral depressional wetland. Surface water and groundwater elevation data were collected for 7 years in a Coastal Plain watershed in South Carolina USA containing depressional wetlands, known as Carolina bays. Rainfall and temperature data were compared with water‐table well and piezometer data in and around one wetland. Using these data a conceptual model was created that describes the hydrology of the system under wet, dry, and drought conditions. The data suggest this wetland operates as a focal point for groundwater recharge under most climate conditions. During years of below‐normal to normal rainfall the hydraulic gradient indicated the potential for groundwater recharge from the depression, whereas during years of above‐normal rainfall, the hydraulic gradient between the adjacent upland, the wetland margin, and the wetland centre showed the potential for groundwater discharge into the wetland. Using high‐resolution water‐level measurements, this groundwater discharge condition was found to hold true even during individual rainfall events, especially under wet antecedent soil conditions. The dynamic nature of the hydrology in this Carolina bay clearly indicates it is not an isolated system as previously believed, and our groundwater data expand upon previous hydrologic investigations at similar sites which do not account for the role of groundwater in estimating the water budget of such systems. Copyright © 2007 John Wiley & Sons, Ltd.     

IHMS—Integrated Hydrological Modelling System. Part 1. Hydrological processes and general structure

A newly Integrated Hydrological Modelling System (IHMS) has been developed to study the impact of changes in climate, land use and water management on groundwater and seawater intrusion (SWI) into coastal areas. The system represents the combination of three models, which can, if required, be run separately. It has been designed to assess the combined impact of climate, land use and groundwater abstraction changes on river, drainage and groundwater flows, groundwater levels and, where appropriate, SWI. The approach is interdisciplinary and reflects an integrated water management approach. The system comprises three packages: the Distributed Catchment Scale Model (DiCaSM), MODFLOW (96 and 2000) and SWI models. In addition to estimating all water balance components, DiCaSM, produces the recharge data that are used as input to the groundwater flow model of the US Geological Survey, MODFLOW. The latter subsequently generates the head distribution and groundwater flows that are used as input to the SWI model, SWI. Thus, any changes in land use, rainfall, water management, abstraction, etc. at the surface are first handled by DiCaSM, then by MODFLOW and finally by the SWI. The three models operate at different spatial and temporal scales and a facility (interface utilities between models) to aggregate/disaggregate input/output data to meet a desired spatial and temporal scale was developed allowing smooth and easy communication between the three models. As MODFLOW and SWI are published and in the public domain, this article focuses on DiCaSM, the newly developed unsaturated zone DiCaSM and equally important the interfacing utilities between the three models. DiCaSM simulates a number of hydrological processes: rainfall interception, evapotranspiration, surface runoff, infiltration, soil water movement in the root zone, plant water uptake, crop growth, stream flow and groundwater recharge. Input requirements include distributed data sets of rainfall, land use, soil types and digital terrain; climate data input can be either distributed or non‐distributed. The model produces distributed and time series output of all water balance components including potential evapotranspiration, actual evapotranspiration, rainfall interception, infiltration, plant water uptake, transpiration, soil water content, soil moisture (SM) deficit, groundwater recharge rate, stream flow and surface runoff. This article focuses on details of the hydrological processes and the various equations used in DiCaSM, as well as the nature of the interface to the MODFLOW and SWI models. Furthermore, the results of preliminary tests of DiCaSM are reported; these include tests related to the ability of the model to predict the SM content of surface and subsurface soil layers, as well as groundwater levels. The latter demonstrates how the groundwater recharge calculated from DiCaSM can be used as input into the groundwater model MODFLOW using aggregation and disaggregation algorithms (built into the interface utility). SWI has also been run successfully with hypothetical examples and was able to reproduce the results of some of the original examples of Bakker and Schaars ( 2005 ). In the subsequent articles, the results of applications to different catchments will be reported. Copyright © 2010 John Wiley & Sons, Ltd.     

Modelling non-stationary water ages in a tropical rainforest: A preliminary spatially distributed assessment

Pristine tropical forests play a critical role in regional and global climate systems. For a better understanding of the eco-hydrology of tropical “evergreen” vegetation, it is essential to know the partitioning of water into transpiration and evaporation, runoff and associated water ages. For this purpose, we evaluated how topography and vegetation influence water flux and age dynamics at high temporal (hourly) and spatial (10 m) resolution using the Spatially Distributed Tracer-Aided Rainfall-Runoff model for the tropics (STARRtropics). The model was applied in a tropical rainforest catchment (3.2 km²) where data were collected biweekly to monthly and during intensive monitoring campaigns from January 2013 to July 2018. The STARRtropics model was further developed, incorporating an isotope mass balance for evapotranspiration partitioning into transpiration and evaporation. Results exhibited a rapid streamflow response to rainfall inputs (water and isotopes) with limited mixing and a largely time-invariant baseflow isotope composition. Simulated soil water storage showed a transient response to rainfall inputs with a seasonal component directly resembling the streamflow dynamics which was independently evaluated using soil water content measurements. High transpiration fluxes (max 7 mm/day) were linked to lower slope gradients, deeper soils and greater leaf area index. Overall water partitioning resulted in 65% of the actual evapotranspiration being driven by vegetation with high transpiration rates over the drier months compared to the wet season. Time scales of water age were highly variable, ranging from hours to a few years. Stream water ages were conceptualized as a mixture of younger soil water and slightly older, deeper soil water and shallow groundwater with a maximum age of roughly 2 years during drought conditions (722 days). The simulated soil water ages ranged from hours to 162 days and for shallow groundwater up to 1,200 days. Despite the model assumptions, experimental challenges and data limitation, this preliminary spatially distributed model study enhances knowledge about the water ages and overall young water dominance in a tropical rainforest with little influence of deeper and older groundwater.     

How important is groundwater availability and stream perenniality to riparian and floodplain tree growth?

Riparian vegetation is important for stream functioning and as a major landscape feature. For many riparian plants, shallow groundwater is an important source of water, particularly in areas where rainfall is low, either annually or seasonally, and when extended dry conditions prevail for all or part of the year. The nature of tree water relationships is highly complex. Therefore, we used multiple lines of evidence to determine the water sources used by the dominant tree species Eucalyptus camaldulensis (river red gum), growing in riparian and floodplain areas with varying depth to groundwater and stream perenniality. Dendrometer bands were used to measure diel, seasonal, and annual patterns of tree water use and growth. Water stable isotopes (δ²H and δ¹⁸O) in plant xylem, soil water, and groundwater were measured to determine spatial and temporal patterns in plant water source use. Our results indicated riparian trees located on relatively shallow groundwater had greater growth rates, larger diel responses in stem diameter, and were less reactive to extended dry periods, than trees in areas of deep groundwater. These results were supported by isotope analysis that suggested all trees used groundwater when soil water stores were depleted at the end of the dry season, and this was most pronounced for trees with shallow groundwater. Trees may experience more frequent periods of water deficit stress and undergo reduced productivity in scenarios where water table accessibility is reduced, such as drawdown from groundwater pumping activities or periods of reduced rainfall recharge. The ability of trees to adapt to changing groundwater conditions may depend on the speed of change, the local hydrologic and soil conditions as well as the species involved. Our results suggest that E. camaldulesis growing at our study site is capable of utilizing groundwater even to depths >10 m, and stream perenniality is likely to be a useful indicator of riparian tree use of groundwater.     

Methods of groundwater recharge estimation in eastern Mediterranean—a water balance model application in Greece,Cyprus and Jordan

Groundwater recharge studies in semi‐arid areas are fundamental because groundwater is often the only water resource of importance. This paper describes the water balance method of groundwater recharge estimation in three different hydro‐climatic environments in eastern Mediterranean, in northwest Greece (Aliakmonas basin/Koromilia basin), in Cyprus (Kouris basin and Larnaka area) and in Jordan (northern part of Jordan). For the Aliakmonas basin, groundwater recharge was calculated for different sub‐catchments. For the Upper Aliakmonas basin (Koromilia basin), a watershed‐distributed model was developed and recharge maps were generated on a daily basis. The mean annual recharge varied between 50 and 75 mm/year (mean annual rainfall 800 mm/year). In Cyprus, the mean groundwater recharge estimates yielded 70 mm/year in the Kouris basin. In the Larnaka area, groundwater recharge ranged from 30 mm/year (lowland) to 200 mm/year (mountains). In Jordan, the results indicated recharge rates ranging from 80 mm/year for very permeable karstified surfaces in the upper part of the Salt basin, where rainfall reaches 500 mm/year to less than 10 mm/year and to only about 1 mm/year in the southernmost part of the basin. For the north part of Jordan, a watershed‐distributed model was developed and recharge maps were generated. This water balance model was used for groundwater recharge estimations in many regions with different climatic conditions and has provided reliable results. It has turned out to be an important tool for the management of the limited natural water resources, which require a detailed understanding of regional hydro(geo)logical processes. Copyright © 2007 John Wiley & Sons, Ltd.     

Soil water dynamics and water balance on a tropical coral island

Understanding soil water dynamics and the water balance of tropical coral islands is important for the utilization and management of their limited freshwater resources, which is only from rainfall. However, there is a significant knowledge gap in the influence of soil water on the water cycle of coral islands. Soil water dynamics and the water balance of Zhaoshu Island, Xisha Archipelago were thus investigated using soil moisture measurements and the Hydrus-1D model from October 2018 to September 2019. Over the study period, vegetation transpiration, soil evaporation, groundwater recharge and storage in the vadose zone were approximately 196, 330, 365 and 20 mm, occupying 22%, 36%, 40% and 2% of annual rainfall total (911 mm), respectively. For the wet season (from May to October) these values became 75, 202, 455 and 40 mm, occupying 10%, 26% and 59% and 5% of the seasonal rainfall total (772 mm), respectively. During the dry season (from November to April), a dry soil layer between 40 and 120 cm depth of the soil profile was identified that prevented water exchange between the upper soil layers and the groundwater resulting in the development of deep roots so that vegetation could extract groundwater to supplement their water requirements. Vegetation not only consumes all dry season rainfall (140 mm) but extracts water deeply from groundwater (90 mm) as well as from the vadose layer (20 mm). As such, the vegetation appears to be groundwater-dependent ecosystems. The research results aid us to better understand the process of water dynamics on coral islands and to protect coral island ecosystems.     

Groundwater recharge estimation in arid hardrock‐alluvium aquifers using combined water‐table fluctuation and groundwater balance approaches

This paper proposes an approach to estimate groundwater recharge using an optimization‐based water‐table fluctuation method combined with a groundwater balance model in an arid hardrock‐alluvium region, located at the Oman–United Arab Emirates border. We introduce an “effective hardrock thickness” term to identify the percentage of the considered hardrock thickness in which effective groundwater flow takes place. The proposed method is based upon a Thiessen polygon zoning approach. The method includes subpolygons to represent specific geologic units and to enhance the confidence of the estimated groundwater recharge. Two linear and 1 nonlinear submodels were developed to evaluate the model components for the calibration (October 1996 to September 2008) and validation (October 2008 to September 2013) periods. Long‐term annual groundwater recharge from rainfall and return flow over the model domain are estimated as 24.62 and 5.71 Mm³, respectively, while the effective groundwater flow circulation is found to occur in the upper 7% of the known hardrock thickness (42 m), confirming conclusions of previous field studies. Considering a total difference in groundwater levels between eastern and western points of the study area of the order of 220 m and a 12‐year monthly calibration period, a weighted root mean squared error in predicted groundwater elevation of 2.75 m is considered quite reasonable for the study area characterized by remarkable geological and hydrogeological diversity. The proposed approach provides an efficient and robust method to estimate groundwater recharge in regions with a complex geological setting in which interaction between fractured and porous media cannot be easily assessed.     

Quantifying groundwater discharge to a small perennial stream in southern Ontario, Canada

A study of the interaction between groundwater and surface water was undertaken within a small agricultural watershed in southern Ontario, Canada. Groundwater contributions to streamflow were measured along a section of stream during baseflow conditions and during rainfall events. Four techniques were used to estimate the contribution of groundwater to the stream along a 450 m reach (three during baseflow and one during stormflow conditions). Under baseflow conditions, streamflow measurements using the velocity–area technique indicated that the net groundwater flux to the stream during the summer months was 10 ml s⁻¹ m⁻¹. Hydrometric measurements (i.e. hydraulic gradient and hydraulic conductivity) taken using mini-piezometers installed in the sediments beneath the stream resulted in net groundwater flux estimates that were four to five times lower. Seepage meters failed to provide any measurements of water flux into or out of the stream. Therefore, based on these results, the velocity–area technique gives the best estimate of groundwater discharge. Hydrograph separations were conducted using isotopic ratios and electrical conductivity on two large rainfall events with different antecedent moisture conditions in the catchment. Both events showed that pre-event water (generally considered groundwater) dominated streamflow and tile drain flow with 64%–80% of the total discharge contributed by pre-event water. High water table conditions within the catchment resulted in greater stream discharge and a greater contribution of event water in the streamflow than that observed under low water table conditions for similar intensity storm events. The results also showed that differences in riparian zone width, vegetation and surface saturation conditions between the upper and lower catchment can influence the relative magnitude of streamflow response from the two catchment areas.     

Temporal dynamics of soil water balance components in a karst range in southeastern Spain: estimation of potential recharge

Abstract

This paper analyses the temporal dynamics of soil water balance components in a representative recharge area of the Sierra de Gádor (Almeria, southeastern Spain) in two hydrological years. Two approaches are used to estimate daily potential recharge (PR): Approach 1 based on deriving PR from the water balance as the difference between measurements of rainfall (P) and actual evapotranspiration (E) obtained by eddy covariance; and Approach 2 with PR obtained from the dynamic pattern of the soil moisture (θ) recorded at two depths in the site's thin soil (average 0.35 m thickess). For the hydrological year 2003/04, which was slightly drier than the 30-year average, E accounted for 64% of rainfall and occurred mainly in late spring and early summer. The PR estimated by Approach 1 was 181 ± 18 mm year^-1 (36% of rainfall), suggesting an effective groundwater recharge in the study area. In the unusually dry hydrological year 2004/05, E was about 215 mm year^-1, close to the annual rainfall input, and allowing very little (8 ± 12 mm year^-1) PR according to Approach 1. Estimation of PR based on Approach 2 resulted in PR rates lower than those found by Approach 1, because Approach 2 does not take into account the recharge that occurs through preferential flow pathways (cracks, joints and fissures) which were not monitored with the θ probes. Moreover, using Approach 2, the PR estimates differed widely depending on the time scale considered: with daily mean θ data, PR estimation was lower, especially in late spring, while θ data at 30 min resolution yielded a more reliable prediction of the fraction of total PR resulting from the downward movement of soil water by gravity.

Citation Cantón, Y., Villagarcía, L., Moro, M. J., Serrano-Ortíz, P., Were, A., Alcalá, F. J., Kowalski, A. S., Solé-Benet, A., Lázaro, R. & Domingo, F. (2010) Temporal dynamics of soil water balance components in a karst range in southeastern Spain: estimation of potential recharge. Hydrol. Sci. J. 55(5), 737–753.     

An integrated water balance model for assessing water scarcity in a data-sparse interfluve in eastern India

Abstract

The objective of this study is to measure the balance of water demand versus water resource availability in an interfluve of West Bengal, India to support water resource planning, particularly of inter-basin transfers. Surface water availability was modelled using the US Soil Conservation Service curve number (SCS-CN) approach, whilst groundwater availability was modelled based on water-level fluctuations and the rainfall infiltration method. Water use was modelled separately for the agricultural, industrial, and domestic sectors using a predominantly normative approach and water use to availability ratios calculated for different administrative areas within the interfluve. Overall, the approach suggested that the interfluve receives 327 × 10⁶ m³ year^-1 of excess water after satisfying these sectoral demands, but that the eastern part of the study area is in deficit. However, a sensitivity analysis carried on the approach to several assumptions in the model suggested changed circumstances would produce surplus/deficit ranging from ?215 × 10⁶ to 435 × 10⁶ m³ year^-1 . The approach could have potential for localised water balance modelling in other Indian catchments.

Editor D. Koutsoyiannis; Associate editor D. Hughes     

Probabilistic risk assessment of contaminant transport in groundwater and vapour intrusion following remediation of a contaminant source

Community-scale simulations were performed to investigate the risk to groundwater and indoor air receptors downgradient of a contaminated site following the remediation of a long-term source. Six suites of Monte Carlo simulations were performed using a numerical model that accounted for groundwater flow, reactive solute transport, soil gas flow, and vapour intrusion in buildings. The model was applied to a three-dimensional, community-scale (250 m × 1000 m × 14 m) domain containing heterogeneous, spatially correlated distributions of the hydraulic conductivity, fraction of organic carbon, and biodegradation rate constant, which were varied between realizations. Analysis considered results from both individual realizations as well as the suite of Monte Carlo simulations expressed through several novel, integrated parameters, such as the probability of exceeding a regulatory standard in either groundwater or indoor air. Results showed that exceedance probabilities varied considerably with the consideration of biodegradation in the saturated zone, and were less sensitive to changes in the variance of hydraulic conductivity or the incorporation of heterogeneous distributions of organic carbon at this spatial scale. A sharp gradient in exceedance probability existed at the lateral edges of the plumes due to variability in lateral dispersion, which defined a narrow region of exceedance uncertainty. Differences in exceedance probability between realizations (i.e., due to heterogeneity uncertainty) were similar to differences attributed to changes in the variance of hydraulic conductivity or fraction of organic carbon. Simulated clean-up times, defined by reaching an acceptable exceedance probability, were found to be on the order of decades to centuries in these community-scale domains. Results also showed that the choice of the acceptable exceedance probability level (e.g., 1 vs. 5 %) would likely affect clean up times on the order of decades. Moreover, in the scenarios examined here, the risk of exceeding indoor air standards was greater than that of exceeding groundwater standards at all times and places. Overall, simulations of coupled transport processes combined with novel spatial and temporal quantification metrics for Monte Carlo analyses, provide practical tools for assessing risk in wider communities when considering site remediation.     

Environmental isotopic and hydrochemical characteristics of groundwater systems in Daying and Qicun geothermal fields,Xinzhou Basin,Shanxi, China

The conceptual hydrogeological model of the low to medium temperature Daying and Qicun geothermal fields has been proposed, based on hydrochemical characteristics and isotopic compositions. The two geothermal fields are located in the Xinzhou basin of Shanxi, China and exhibit similarities in their broad‐scale flow patterns. Geothermal water is derived from the regional groundwater flow system of the basin and is characterized by Cl·SO₄‐Na type. Thermal water is hydrochemically distinct from cold groundwater having higher total dissolved solids (TDS) (>0·8 g/l) and Sr contents, but relatively low Ca, Mg and HCO₃ contents. Most shallow groundwater belongs to local flow systems which are subject to evaporation and mixing with irrigation returns. The groundwater residence times estimated by tritium and ¹⁴C activities indicate that deep non‐thermal groundwater (130–160 m) in the Daying region range from modern (post‐1950s) in the piedmont area to more than 9·4 ka BP (Before Present) in the downriver area and imply that this water belong to an intermediate flow system. Thermal water in the two geothermal fields contains no detectable active ¹⁴C, indicating long residence times (>50 ka), consistent with this water being part of a large regional flow system. The mean recharge elevation estimated by using the obtained relationship Altitude (m) = ? 23·8 × δ²H (‰ ) ? 121·3, is 1980 and 1880 m for the Daying and Qicun geothermal fields, respectively. The annual infiltration rates in the Daying and Qicun geothermal fields can be estimated to be 9029 × 10³ and 4107 × 10³ m³/a, respectively. The variable ⁸⁶Sr/⁸⁷Sr values in the thermal and non‐thermal groundwater in the two fields reflect different lithologies encountered along the flow path(s) and possibly different extents of water‐rock interaction. Based on the analysis of groundwater flow systems in the two geothermal fields, hydrogeochemical inverse modelling was performed to indicate the possible water‐rock interaction processes that occur under different scenarios. Copyright © 2010 John Wiley & Sons, Ltd.     

Evaluating the interactions between surface water and groundwater in the arid mid-eastern Yanqi Basin,northwestern China

The environment of Bosten Lake in the Mid-Eastern Yanqi Basin (MEYB), an arid inland area in northwest China, has deteriorated greatly due to increasing groundwater exploitation and changes in the interactions between groundwater and surface water. This study intended to simulate the spatio-temporal variability of groundwater and surface water across the entire MEYB over the period 2000–2013. The applicable groundwater flow model and mass balance calculation method for river water were constructed to evaluate the change in groundwater recharged by and discharged to different segments of the Kaidu River. Simulation results show that the entire river seepage in the MEYB increased from 1.05 to 6.17 × 10⁸ m³/year between 2000 and 2013. The increasing river seepage, induced by increasing groundwater exploitation, plays the most important role in the water level decline in the downstream reaches of the Kaidu River and in Bosten Lake. This implies that the current utilization of groundwater resources in the MEYB is unsustainable.     

Coupled surface water–groundwater model and its application in the arid Shiyang River basin,China

The planning and management of water resources in the Shiyang River basin, China require a tool for assessing the impact of groundwater and stream use on water supply reliabilities and improving many environment‐related problems such as soil desertification induced by recent water‐related human activities. A coupled model, integrating rule‐based lumped surface water model and distributed three‐dimensional groundwater flow model, has been established to investigate surface water and groundwater management scenarios that may be designed to restore the deteriorated ecological environment of the downstream portion of the Shiyang River basin. More than 66% of the water level among 24 observation wells have simulation error less than 1·0 m. The overall trend of the temporal changes of simulated and observed surface runoff at the Caiqi gauging station remains almost the same. The calibration was considered satisfactory. Initial frameworks for water allocation, including agricultural water‐saving projects, water diversion within the basin and inter‐basin water transfer, reducing agricultural irrigation area and surface water use instead of groundwater exploitation at the downstream were figured out that would provide a rational use of water resources throughout the whole basin. Sixteen scenarios were modelled to find out the most appropriate management strategies. The results showed that in the two selected management options, the groundwater budget at the Minqin basin was about 1·4 × 10⁸ m³/a and the ecological environment would be improved significantly, but the deficit existed at the Wuwei basin and the number was about 0·8 × 10⁸ m³/a. Water demand for domestic, industry and urban green area would be met in the next 30 years, but the water shortage for meeting the demand of agricultural water use in the Shiyang River basin was about 2·2 × 10⁸ m³/a. It is suggested that more inter‐basin water transfer should be required to obtain sustainable water resource use in the Shiyang River basin. Copyright © 2009 John Wiley & Sons, Ltd.