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.     

Assessment of transmission losses and groundwater recharge from runoff events in a wadi under shortage of data on lateral inflow,Negev, Israel

A hydrological–lithostratigraphical model was developed for assessment of transmission losses and groundwater recharge from runoff events in arid water courses where hydrological and meteorological records are incomplete. Water balance equations were established for reaches between hydrometric stations. Because rainfall and tributary flow data are scarce, lateral inflow, which is an essential component of the water balance equation, could not be estimated directly. The solution was obtained by developing a method which includes a hydrological–lithostratigraphical analogy. This is based on the following assumptions: (a) runoff resulting from a given rainfall event is related to the watershed surface lithology; (b) for a given event, the spatial distribution of runoff reflects the distribution of rainfall: and (c) transmission losses are uniquely related to the total inflow to the reach. The latter relationship, called the loss function, and the water balance equation comprise a model which simultaneously assesses lateral inflow and transmission losses for runoff events recorded at the terminal stations. The model was applied to three reaches of the arid Nahal Tsin in Israel. In this case study, the transmission losses were of the same order of magnitude as the flow at the major hydrometric stations. The losses were subdivided into channel moistening, which subsequently evaporates, and deep percolation, which recharges groundwater. For large runoff events, evaporation was substantially smaller than the losses. The mean annual recharge of groundwater from runoff events in the Tsin watershed was 4·1×10⁶ m³, while the mean annual flow volume at the major stations ranged from 0·6 to 1·5×10⁶ m³. Once in 100 years, the annual recharge may be seven times higher than the mean annual value, but the recharge during most years is very small. Copyright © 1999 John Wiley & Sons, Ltd.     

A statistically based runoff‐yield model coupling infiltration excess and saturation excess mechanisms

A statistically based runoff‐yield model is proposed in this paper. The model considers spatial heterogeneities of rainfall, soil infiltration capacity and soil water storage capacity that are main factors controlling runoff‐yield process. It assumes that the spatial variation of rainfall intensity at each time step can be characterized by a probability density function, which is estimated by matching the hyetograph through goodness‐of‐fit measure, whereas the spatial heterogeneities of soil infiltration capacity and soil water storage capacity are described by parabola‐type functions. Surface runoff is calculated according to infiltration excess mechanism; the statistical distribution of surface runoff rate can be deduced with the joint distribution of rainfall intensity and soil infiltration rate, thus obtaining a quasi‐analytical solution for surface runoff. Based on saturation excess mechanism, the groundwater flow (flows below the ground are collectively referred to as groundwater flow) is calculated by infiltration and the probability distribution of soil water storage capacity. Consequently, the total runoff is composed of infiltration excess and saturation excess runoff components. As an example, this model is applied to flood event simulation in Dongwan catchment, a semi‐humid region and a tributary of Yellow River in China. It indicates that the proposed runoff‐yield model could achieve acceptable accuracy. Copyright © 2011 John Wiley & Sons, Ltd.     

Role of epikarst in near‐surface hydrological processes in a soil mantled subtropical dolomite karst slope: implications of field rainfall simulation experiments

Epikarst exerts a strong control on run‐off generation in karst regions, but it is still unclear in karst regions. Our study aimed to demonstrate the effect of epikarst on near‐surface hydrological processes in a subtropical cockpit karst region of southwest China, using plot‐scale rainfall simulation experiments with different rainfall intensities (low and high) and antecedent moisture conditions (dry and wet). A trench excavated to the epikarst lower boundary allowed identification of flow pathways in the entire soil–epikarst architecture system, thus facilitating the water balance calculations using a conceptual model with the assumption of a two‐stage hydrological evolution. More than 70% of the total rainfall water moved vertically through the shallow soil layer and then was redistributed by the epikarst as subsurface flow occurring on the soil–epikarst interface, depression filling on epikarst surface, water held by epikarst and deep percolation. Epikarst water regulation capacity, defined as the sum of depression filling on epikarst surface, water held by epikarst, epikarst seepage flow and deep percolation, was 58 mm (wet antecedent condition) and 223 mm (dry antecedent condition). Total run‐off from the soil–epikarst system was dominated by saturated subsurface flow showing a threshold process controlled by epikarst storage capacity (storing as much as 181 mm of rainfall water under dry antecedent condition). Our study proved that despite the epikarst being relatively poorly developed and covered by a soil mantle, it still exerted a strong influence on near‐surface hydrological processes and thus should be adequately considered in future modelling of water recharge and depletion dynamics in this integrated soil–epikarst system. Copyright © 2015 John Wiley & Sons, Ltd.     

Tracer hydrology of the data-scarce and heterogeneous Central American Isthmus

Numerous socio-economic activities depend on the seasonal rainfall and groundwater recharge cycle across the Central American Isthmus. Population growth and unregulated land use changes resulted in extensive surface water pollution and a large dependency on groundwater resources. This work combines stable isotope variations in rainfall, surface water, and groundwater of Costa Rica, Nicaragua, El Salvador, and Honduras to develop a regionalized rainfall isoscape, isotopic lapse rates, spatial–temporal isotopic variations, and air mass back trajectories determining potential mean recharge elevations, moisture circulation patterns, and surface water–groundwater interactions. Intra-seasonal rainfall modes resulted in two isotopically depleted incursions (W-shaped isotopic pattern) during the wet season and two enriched pulses during the mid-summer drought and the months of the strongest trade winds. Notable isotopic sub-cloud fractionation and near-surface secondary evaporation were identified as common denominators within the Central American Dry Corridor. Groundwater and surface water isotope ratios depicted the strong orographic separation into the Caribbean and Pacific domains, mainly induced by the governing moisture transport from the Caribbean Sea, complex rainfall producing systems across the N-S mountain range, and the subsequent mixing with local evapotranspiration, and, to a lesser degree, the eastern Pacific Ocean fluxes. Groundwater recharge was characterized by (a) depleted recharge in highland areas (72.3%), (b) rapid recharge via preferential flow paths (13.1%), and enriched recharge due to near-surface secondary fractionation (14.6%). Median recharge elevation ranged from 1,104 to 1,979 m a.s.l. These results are intended to enhance forest conservation practices, inform water protection regulations, and facilitate water security and sustainability planning in the Central American Isthmus.     

Assessment of groundwater recharge processes through karst vadose zone by cave percolation monitoring

Recharge processes of karst aquifers are difficult to assess given their strong heterogeneity and the poorly known effect of vadose zone on infiltration. However, recharge assessment is crucial for the evaluation of groundwater resources. Moreover, the vulnerability of karst aquifers depends on vadose zone behaviour because it is the place where most contamination takes place. In this work, an in situ experimental approach was performed to identify and quantify flow and storage processes occurring in karst vadose zone. Cave percolation monitoring and dye tracing were used to investigate unsaturated zone hydrological processes. Two flow components (diffuse and quick) were identified and, respectively, account for 66% and 34% of the recharge. Quickflow was found to be the result of bypass phenomenon in vadose zone related to water saturation. We identify the role of epikarst as a shunting area, most of the storage in the vadose zone occurring via the diffuse flow component in low permeability zones. Relationship between rainfall intensity and transit velocity was demonstrated, with 5 times higher velocities for the quick recharge mode than the diffuse mode. Modelling approach with KarstMod software allowed to simulate the hybrid recharge through vadose zone and shows promising chances to properly assess the recharge processes in karst aquifer based on simple physical models.     

Groundwater Recharge at Five Representative Sites in the Hebei Plain,China

Accurate estimates of groundwater recharge are essential for effective management of groundwater, especially when supplies are limited such as in many arid and semiarid areas. In the Hebei Plain, China, water shortage is increasingly restricting socioeconomic development, especially for agriculture, which heavily relies on groundwater. Human activities have greatly changed groundwater recharge there during the past several decades. To obtain better estimates of recharge in the plain, five representative sites were selected to investigate the effects of irrigation and water table depth on groundwater recharge. At each site, a one‐dimensional unsaturated flow model (Hydrus‐1D) was calibrated using field data of climate, soil moisture, and groundwater levels. A sensitivity analysis of evapotranspirative fluxes and various soil hydraulic parameters confirmed that fine‐textured surface soils generally generate less recharge. Model calculations showed that recharge on average is about 175 mm/year in the piedmont plain to the west, and 133 mm/year in both the central alluvial and lacustrine plains and the coastal plain to the east. Temporal and spatial variations in the recharge processes were significant in response to rainfall and irrigation. Peak time‐lags between infiltration (rainfall plus irrigation) and recharge were 18 to 35 days in the piedmont plain and 3 to 5 days in the central alluvial and lacustrine plains, but only 1 or 2 days in the coastal plain. This implies that different time‐lags corresponding to different water table depths must be considered when estimating or modeling groundwater recharge.     

Hedgerow impacts on soil-water transfer due to rainfall interception and root-water uptake

Hedgerow is one of the most important rural landscapes in the world, especially in Europe. Knowledge about the hydrological role of hedgerows is useful in many fields of study, such as hydrological modelling and rural landscape management. The aim of this study was to investigate the impact of a hedgerow on rainfall distribution, soil-water potential gradient, lateral water transfer and water balance. A hillslope with a hedgerow perpendicular to the slope was monitored. To evaluate hedgerow rainfall interception, rainfall was measured (hourly, daily, and by rainfall event) both next to and up to 16 m upslope and 12 m downslope perpendicularly away from the hedgerow. The strongest correlation between rainfall next to the hedgerow and rainfall at more distant points was obtained using data measured per rainfall event. The average percentage of rainfall intercepted equalled 28% for the leafed period and 12% for the leafless period. The impact of the hedgerow on spatial rainfall distribution was related to distance from the hedgerow and rainfall amount. Annual distribution of soil-water potential showed that the hedgerow influenced it up to 9 m upslope and 6 m downslope, the area in which most of the hedgerow's roots were observed. The soil was driest at the end of summer, which delayed soil rewetting during autumn in areas surrounding the hedgerow. Annual groundwater dynamics exhibited three distinct periods due to temporal rainfall distribution and, especially at the end of summer, root-water uptake. In addition, the total potential gradient showed that unsaturated flow was directed towards the hedgerow in summer and autumn. These results indicate that at the local scale hedgerows influences (1) spatial rainfall distribution, (2) soil rewetting, and (3) groundwater recharge, often at distances well beyond the hedgerow's drip line. Consequently, the processes involved in soil-water dynamics around hedgerows should be integrated into relevant hydrological models, especially for catchments with a dense hedgerow network. Copyright © 2008 John Wiley & Sons, Ltd.     

Effects of Linking a Soil‐Water‐Balance Model with a Groundwater‐Flow Model

A previously published regional groundwater‐flow model in north‐central Nebraska was sequentially linked with the recently developed soil‐water‐balance (SWB) model to analyze effects to groundwater‐flow model parameters and calibration results. The linked models provided a more detailed spatial and temporal distribution of simulated recharge based on hydrologic processes, improvement of simulated groundwater‐level changes and base flows at specific sites in agricultural areas, and a physically based assessment of the relative magnitude of recharge for grassland, nonirrigated cropland, and irrigated cropland areas. Root‐mean‐squared (RMS) differences between the simulated and estimated or measured target values for the previously published model and linked models were relatively similar and did not improve for all types of calibration targets. However, without any adjustment to the SWB‐generated recharge, the RMS difference between simulated and estimated base‐flow target values for the groundwater‐flow model was slightly smaller than for the previously published model, possibly indicating that the volume of recharge simulated by the SWB code was closer to actual hydrogeologic conditions than the previously published model provided. Groundwater‐level and base‐flow hydrographs showed that temporal patterns of simulated groundwater levels and base flows were more accurate for the linked models than for the previously published model at several sites, particularly in agricultural areas.     

Is the connectivity function a good indicator of soil infiltrability distribution and runoff flow dimension?

Among the studies on runoff connectivity of soils with heterogeneous properties, the need to understand the relationships between soil heterogeneity and the associated runoff organization and amount is frequently mentioned. In this study, we simulate the stationary runoff–runon process on bi‐dimensional (2D) flat slopes for five infiltrability distributions, one of them correlated, as a function of rainfall intensity and flow dimension. We define flow dimension by 1 + ε, where ε is the outflow fraction transferred from one pixel to each of the two lateral downslope pixels. Our aim is to assess the effect of ε and soil heterogeneity on the connectivity function compared to the mean runoff flow rate, the wet area and the number of runoff patterns. The analysis of connectivity is carried within the percolation framework. The results show that the integral connectivity scale is more sensitive to the flow dimension and soil heterogeneity compared to the other variables. The wet area fraction does not depend on ε. Unlike previous studies, we find that increased runoff production is not necessarily related to increased connectivity. The use of the connectivity function within the percolation framework appears to be a valuable method for assessing the impact of soil heterogeneity and flow dimension on the runoff organization during a rainfall event. Copyright © 2014 John Wiley & Sons, Ltd.     

Effects of meteorological and hydrogeological factors on gross recharge percentage at unconfined sandy aquifers with an equatorial climate

A thorough understanding of rainfall recharge processes and their controlling factors is essential for management of groundwater systems. This study investigates the effects of various meteorological and hydrogeological factors on the gross recharge percentages, the rainfall–recharge relationships and the recharge threshold values for unconfined sandy aquifers under an equatorial climate. Among the meteorological factors investigated, rainfall intensity was found to have the most significant impact on the gross recharge rate. The effects of potential evaporation rate, relative humidity and air temperature on the gross recharge percentage were significant when the vadose zone thickness is larger than 2·5 m. The recharge threshold values were found to depend strongly on the vadose zone thickness. The rainfall–recharge relationships could generally be well defined by a normal–log relationship. The rainfall–recharge relationships derived here are applicable to yield estimates of gross recharge percentages for unconfined sandy aquifers under an equatorial climate, using rainfall intensity and vadose zone thickness as input variables. In this study, a theory was developed and validated to provide physical explanations for the observations, based on the residence time of the percolated rainwater within the vadose zone. Among the soil hydraulic parameters tested, porosity and saturated hydraulic conductivity were found to have the most pronounced effects on the gross recharge percentage. Utilizing the sensitivity results and the theory derived, an approach was developed for extending the application of the derived rainfall–recharge relationships to other sand textures. The approach was found to be capable of producing rough and fast estimations of gross recharge percentage for other sand textures. Copyright © 2007 John Wiley & Sons, Ltd.     

Grey degree and grey prediction of groundwater head

Boundary types and geologic conditions, which possess random and obscure characteristics, produce variations in groundwater heads. In this study, a groundwater flow system was regarded as a grey system and a grey degree defined and quantified those characteristics. Data were obtained from field records in Gao-Shu and De-Xie Stations, which are located in the upper and middle sites of Pingtung Plain, Taiwan, respectively. An estimated interval was used to represent the upper and lower groundwater head limits. The analyses showed that the grey degree in the wet season was smaller than in the dry season. This was due to the greater uncertainty in groundwater pumping and rainfall recharge in the dry season. With high confidence levels, the grey degree of the groundwater heads decreased. Under spatial distribution, the uncertainty of the groundwater pumping and rainfall recharge in the middle zone was greater than in the upper zone because the middle zone included residential districts and agricultural regions. Thus, the grey degree in the upper zone was less than the middle zone. The grey model can be used to predict a groundwater head based on the observed groundwater head data and the random degree in the groundwater head each month can be judged by the grey interval.     

Water percolation process studies in a Mediterranean karst area

To investigate processes of water percolation, the drip response of stalactites in a karstic cave below a 143 m² sprinkling plot was measured. The experiment was conducted in Mount Carmel, Israel, at the end of the dry season and intended to simulate a series of two high‐intensity storms on dry and wet soils. In addition to hydrometric measurements (soil moisture, surface runoff, stalactite dripping rates), two types of tracers (electrical conductivity and bromide) were used to study recharge processes, water origin and mixing inside a 28‐m vadose zone. Results suggested that slow, continuous percolation through the rock matrix is of minor importance and that percolating water follows a complicated pattern including vertical and horizontal flow directions. While bromide tracing allowed identification of quick direct flow paths at all drips with maximum flow velocities of 4·3 m/h, mixing analysis suggested that major water fractions were mobilized by piston flow, pushing out water stored in the unsaturated zone above the cave. Under dry preconditions, 80 mm of artificial rainfall applied in less than 7 h was not enough to initiate significant downward water percolation. Most water was required to fill uppermost soil and rock storages. Under wet preconditions during the second day sprinkling, higher water contents in soils and karst cavities facilitated piston flow effects and a more intense response of the cave drips. Results indicate that in Mediterranean karst regions, filling of the unsaturated zone, including soil and rock storages, is an important precondition for the onset of significant water percolation and recharge. This results in a higher seasonal threshold for water percolation than for the generation of surface runoff. Copyright © 2010 John Wiley & Sons, Ltd.     

The mechanism of unsaturated flow through a volcanic ash layer under humid climatic conditions

It is very interesting and meaningful to investigate the rainfall-groundwater recharge process under the humid climatic condition of Japan, where mean annual precipitation is about 1600 mm. The present study has investigated soil water movement in the unsaturated zones of a volcanic ash layer, called the ‘Kanto Loam formation’, using environmental tritium as a tracer. The site selected is a flat ground surface on a terraced upland which has a deep unsaturated zone (about 20 m) with a relatively high water content (about 70 per cent) consisting of nearly uniform Kanto Loam formation. The tritium concentrations in groundwater, soil waters having different matric potentials, precipitation, and the seepage water moving through the formation into a man-made cave were measured to characterize the rainfall-groundwater recharge process and the effect of large pore spaces in the formation mentioned by previous studies. Because of the humid climate of Japan, there appears to be a unique soil water flow characteristic which may involve percolation through large pore spaces during heavy rainfall. However, in a fine grained and high water content soil like the Kanto Loam formation, the existence of this flow through large pore spaces does not have a significant effect upon the whole recharge process. The recharge model of displacement flow with dispersion is useful in estimating the tritium concentration profile of soil water. The calculated result shows a recharge rate of 2.5 mm/day. The value obtained reflects the hydrological characteristics of the uplands covered with volcanic ash.     

The impact of climate change on groundwater recharge and runoff in a humid,equatorial catchment: sensitivity of projections to rainfall intensity

Abstract

Projected warming in equatorial Africa, accompanied by greater evaporation and more frequent heavy precipitation events, may have substantial but uncertain impacts on terrestrial hydrology. Quantitative analyses of climate change impacts on catchment hydrology require high-resolution (<50 km) climate data provided by regional climate models (RCMs). We apply validated precipitation and temperature data from the RCM PRECIS (Providing Regional Climates for Impact Studies) to a semi-distributed soil moisture balance model (SMBM) in order to quantify the impacts of climate change on groundwater recharge and runoff in a medium-sized catchment (2098 km²) in the humid tropics of southwestern Uganda. The SMBM explicitly accounts for changes in soil moisture, and partitions effective precipitation into groundwater recharge and runoff. Under the A2 emissions scenario (2070–2100), climate projections from PRECIS feature not only rises in catchment precipitation and modelled potential evapotranspiration by 14% and 53%, respectively, but also increases in rainfall intensity. We show that the common application of the historical rainfall distribution using delta factors to the SMBM grossly underestimates groundwater recharge (i.e. 55% decrease relative to the baseline period of 1961–1990). By transforming the rainfall distribution to account for changes in rainfall intensity, we project increases in recharge and runoff of 53% and 137%, respectively, relative to the baseline period.     

Analytic solutions to transient groundwater flow under time-dependent sources in a heterogeneous aquifer bounded by fluctuating river stage

Analytical solutions for the water table and lateral discharge in a heterogeneous unconfined aquifer with time-dependent source and fluctuating river stage were derived and compared with those in an equivalent homogeneous aquifer. The heterogeneous aquifer considered consists of a number of sections of different hydraulic conductivity values. The source term and river stage were assumed to be time-dependent but spatially uniform. The solutions derived is useful in studying various groundwater flow problems in a horizontally heterogeneous aquifer since the spatially piecewise-constant hydraulic conductivity and temporally piecewise-constant recharge and lateral discharge can be used to quantify variations in these processes commonly observed in reality. Applying the solutions derived to an aquifer of three sections of different hydraulic conductivity values shown that (1) the aquifer heterogeneity significantly increases the spatial variation of the water table and thus its gradient but it has little effect on lateral discharge in the case of temporally and spatially uniform recharge, (2) the time-dependent but spatially uniform recharge increases the temporal variation of groundwater table over the entire aquifer but its effect on lateral discharge is limited in the zone near the river, and (3) the effect of river stage fluctuation on the water table and lateral discharge is limited in the zone near the river and the effect of the heterogeneity is to increase lateral discharge to or recharge from the river.     

The groundwater–land-surface–atmosphere connection: Soil moisture effects on the atmospheric boundary layer in fully-coupled simulations

This study combines a variably-saturated groundwater flow model and a mesoscale atmospheric model to examine the effects of soil moisture heterogeneity on atmospheric boundary layer processes. This parallel, integrated model can simulate spatial variations in land-surface forcing driven by three-dimensional (3D) atmospheric and subsurface components. The development of atmospheric flow is studied in a series of idealized test cases with different initial soil moisture distributions generated by an offline spin-up procedure or interpolated from a coarse-resolution dataset. These test cases are performed with both the fully-coupled model (which includes 3D groundwater flow and surface water routing) and the uncoupled atmospheric model. The effects of the different soil moisture initializations and lateral subsurface and surface water flow are seen in the differences in atmospheric evolution over a 36-h period. The fully-coupled model maintains a realistic topographically-driven soil moisture distribution, while the uncoupled atmospheric model does not. Furthermore, the coupled model shows spatial and temporal correlations between surface and lower atmospheric variables and water table depth. These correlations are particularly strong during times when the land-surface temperatures trigger shifts in wind behavior, such as during early morning surface heating.     

What is groundwater and what does this mean to fauna? - An opinion

Many subsurface waters are considered groundwater but are influenced in shallow depths by hyporheic, parafluvial and/or soil interception water to such a degree that groundwater fauna (stygofauna) communities may be significantly altered. Recharge, even if spatially and temporally distinct, delivers input of dissolved oxygen, organic matter (OM), and nutrients that caters sustainably for ubiquists such as stygophiles and hyporheic fauna, but renders the life of uncompetitive stygobites difficult or impossible. The impact of recharge at shallow groundwater thus needs to be taken into account when determining groundwater fauna reference communities and when evaluating monitoring studies.One of the main characteristics of groundwater is low OM concentration. In contrast, high OM concentrations are typical of hyporheic or parafluvial waters, which are enriched by OM from the river, the riparian soils and from interflow, and which contribute significantly to river OM balance. Consequently, for ecological studies on subsurface waters, both the origin of the water and OM, and the intensity of surface water interactions should be considered. Here, we discuss how groundwater spatial and temporal heterogeneity translates into faunal distribution patterns. In terms of the origin of water and OM, and from an ecological point of view, we need to distinguish between (i) shallow groundwater characterized by infiltrating precipitation and soil recharge, (ii) shallow groundwater interacting with surface water bodies such as continuously flowing and ephemeral streams and rivers, and (iii) “old” groundwater which has no recent connections to the surface and is thus largely secluded from input of nutrients and carbon. Water in the first two groups is characterized by high amounts of OM of varying quality, while water in the third group is characterized by low amounts of low quality OM. Consequently, stygophiles dominate in groups 1 and 2, with hyporheic fauna taking up a considerable proportion in group 2, while stygobites only dominate in group 3. Thus, for studies aiming to assess impacts on groundwater, only sampling sites of the third group should be used for reference sites as these are the most likely sites to have little surface impact and a stygofauna representative of the deeper aquifer.     

The groundwater recharge regime of some slightly metamorphosed neoproterozoic sedimentary rocks: an application of natural environmental tracers

Isotope data of precipitation and groundwater in parts of the Voltaian Basin in Northern Ghana were used to explain the groundwater recharge regime in the area. Groundwater recharge is an important parameter in the development of a decision support system for the management and efficient utilization of groundwater resources in the area. It is therefore important to establish the processes and sources of groundwater recharge. δ¹⁸O and δ²H data for local precipitation suggest enrichment relative to the Global Meteoric Water Line (GMWL) and indicate that precipitation takes place at a relative humidity less than 100%. The groundwater data plot on an evaporation line with a slope of 5, suggesting a high degree of evaporative enrichment of the precipitation in the process of vertical infiltration and percolation through the unsaturated zone into the saturated zone. This finding is consistent with the observation of high evapotranspiration rates in the area and ties in with the fact that significant clay fraction in the unsaturated zone limits vertical percolation and thus exposes the percolating rainwater to the effects of high temperatures and low humidities resulting in high evapotranspiration rates. Groundwater recharge estimates from the chloride mass balance, CMB, method suggest recharge in the range of 1.8–32% of the annual average precipitation in the form of rainfall. The highest rates are associated with areas where open wells encourage significant amount of groundwater recharge from precipitation in the area. In the northern parts of the study area, groundwater recharge is lower than 12%. The recharge so computed through the application of the CMB methodology takes on a spatial distribution akin to the converse of the spatial pattern of both δ¹⁸O and δ²H in the area. As such, the locations of the highest recharge are associated with the most depleted values of the two isotopes. This observation is consistent with the assertion that low vertical hydraulic conductivities slow down vertical percolation of precipitation down to the groundwater water. The percolating precipitation water thus gets enriched in the heavier isotopes through high evapotranspiration rates. At the same time, the amount of water that finally reaches the water table is considerably reduced. Copyright © 2013 John Wiley & Sons, Ltd.     

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).