Significance of stemflow in groundwater recharge. 1: Evaluation of the stemflow contribution to recharge using a mass balance approach

Stemflow was evaluated in a water balance and its contribution to groundwater recharge determined. Gross precipitation, throughfall and stemflow were measured for one year in a pine forest (Tsukuba, Japan) to determine each component of the water balance in the forest. Groundwater recharge rates by stemflow and throughfall were calculated from a mass balance method using chloride in subsurface waters. The stemflow in the water balance was relatively small when estimated as a value per canopy projected area of the tree in the forest. However, the results for the mass balance of chloride in subsurface waters indicated that it was impossible to disregard the stemflow in determining groundwater recharge. Although the ratio of stemflow to the net precipitation was small in the water balance, the effect of stemflow on groundwater recharge was relatively large.     

Hydrological modelling for assessing spatio-temporal groundwater recharge variations in the water-stressed Amathole Water Supply System,Eastern Cape,South Africa

To increase the resilience of regional water supply systems in South Africa in the face of anticipated climatic changes and a constant increase in water demand, water supply sources require diversification. Many water-stressed metropolitan regions in South Africa depend largely on surface water to cover their water demand. While climatic and river discharge data is widely available in these regions, information on groundwater resources – which could support supply source diversification – is scarce. Groundwater recharge is a key parameter that is used to estimate groundwater amounts that can be sustainably exploited at a sub-watershed level. Therefore, the objective of this study was to develop a reliable hydrological modelling routine that enables the assessment of regional spatio-temporal variations of groundwater recharge to discern the most promising areas for groundwater development. Accordingly, we present a semi-distributed hydrological modelling approach that incorporates water balance routines coupled with baseflow modelling techniques to yield spatio-temporal variations of groundwater recharge on a regional level. The approach is demonstrated for the actively managed catchment areas of the Amathole Water Supply System situated in a semi-arid part of the Eastern Cape of South Africa. In the investigated study area, annual groundwater recharge exhibits a high spatio-temporal heterogeneity and is estimated to vary between ~0.5% and 8% of annual rainfall. Despite some uncertainties induced by limited data availability, calibration and validation of the model were found to be satisfactory and yielded model results similar to (point) data of annual groundwater recharge reported in earlier studies. Our approach is therefore found to derive crucial information for efficiently targeting more detailed groundwater exploration studies and could work as a blueprint for orientating groundwater potential exploration in similar environments.     

Impacts of land-use and climate changes on surface runoff in a tropical forest watershed (Brazil)

ABSTRACT

Surface runoff generation capacity can be modified by land-use and climate changes. Annual runoff volumes have been evaluated in a small watershed of tropical forest (Brazil), using the Soil and Water Assessment Tool (SWAT) model. Firstly, the accuracy of SWAT in runoff predictions has been assessed by default input parameters and improved by automatic calibration, using 20-year observations. Then, the hydrological response under land uses (cropland, pasture and deforested soil) alternative to tropical forest and climate change scenarios has been simulated. SWAT application has showed that, if forest was replaced by crops or pasture, the watershed’s hydrological response would not significantly be affected. Conversely, a complete deforestation would slightly increase its runoff generation capacity. Under forecasted climate scenarios, the runoff generation capacity of the watershed will tend to decrease and will not be noticeably different among the representative concentration pathways. Pasture and bare soil will give the lowest and highest runoff coefficients, respectively.     

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.     

Evaluating the use of a gridded climate surface for modelling groundwater recharge in a semi‐arid region (Okanagan Basin,Canada)

Spatially distributed groundwater recharge was simulated for a segment of a semi‐arid valley using three different treatments of meteorological input data and potential evapotranspiration (PET). For the same area, timeframe, land cover characteristics and soil properties, groundwater recharge was estimate using (i) single‐station climate data with monthly PET calculated by the Thornthwaite method; (ii) single‐station climate data with daily PET calculated by the Penman–Monteith method; and (iii) daily gridded climate data with spatially distributed PET calculated using the Penman–Monteith method. For each treatment, the magnitude and distribution of actual evapotranspiration (AET) for summer months compared well with those estimated for a 5‐year crop study, suggesting that the near‐surface hydrological processes were replicated and that subsequent groundwater recharge rates are realistic. However, for winter months, calculated AET was near zero when using the Thornthwaite PET method. Mean annual groundwater recharge varied from ～3·2 to 10·0 mm when PET was calculated by the Thornthwaite method, and from ～1·8 to 7·5 mm when PET was calculated by the Penman–Monteith method. Comparisons of bivariate plots of seasonal recharge rates estimated from single‐station versus gridded surface climate reveal that there is greater variability between the different methods for spring months, which is the season of greatest recharge. Furthermore, these seasonal differences are shown to provide different results when compared to the depth to water table, which could lead to different results of evaporative extinction depth. These findings illustrate potential consequences of using different approaches for representing spatial meteorological input data, which could provide conflicting predictions when modelling the influence of climate change on groundwater recharge. Copyright © 2010 John Wiley & Sons, Ltd.     

Water uptake by trees in a riparian hardwood forest (Rhine floodplain,France)

Water flow in the soil–root–stem system was studied in a flooded riparian hardwood forest in the upper Rhine floodplain. The study was undertaken to identify the vertical distribution of water uptake by trees in a system where the groundwater is at a depth of less than 1 m. The three dominant ligneous species (Quercus robur, Fraxinus excelsior and Populus alba) were investigated for root structure (vertical extension of root systems), leaf and soil water potential (Ψ_m), isotopic signal (¹⁸O) of soil water and xylem sap. The root density of oak and poplar was maximal at a depth of 20 to 60 cm, whereas the roots of the ash explored the surface horizon between 0 and 30 cm, which suggests a complementary tree root distribution in the hardwood forest. The flow density of oak and poplar was much lower than that of the ash. However, in the three cases the depth of soil explored by the roots reached 1·2 m, i.e. just above a bed of gravel. The oak roots had a large lateral distribution up to a distance of 15 m from the trunk. The water potential of the soil measured at 1 m from the trunk showed a zone of strong water potential between 20 and 60 cm deep. The vertical profile of soil water content varied from 0·40 to 0·50 cm³ cm^?3 close to the water table, and 0·20 to 0·30 cm³ cm^?3 in the rooting zone. The isotopic signal of stem water was constant over the whole 24‐h cycle, which suggested that the uptake of water by trees occurred at a relatively constant depth. By comparing the isotopic composition of water between soil and plant, it was concluded that the water uptake occurred at a depth of 20 to 60 cm, which was in good agreement with the root and soil water potential distributions. The riparian forest therefore did not take water directly from the water table but from the unsaturated zone through the effect of capillarity. Copyright © 2007 John Wiley & Sons, Ltd.     

A comparison of groundwater recharge estimation methods in a semi‐arid,coastal avocado and citrus orchard (Ventura County,California)

We assess the relative merits of application of the most commonly used field methods (soil‐water balance (SWB), chloride mass balance (CMB) and soil moisture monitoring (NP)) to determine recharge rates in micro‐irrigated and non‐irrigated areas of a semi‐arid coastal orchard located in a relatively complex geological environment. Application of the CMB method to estimate recharge rates was difficult owing to the unusually high, variable soil‐water chloride concentrations. In addition, contrary to that expected, the chloride concentration distribution at depths below the root zone in the non‐irrigated soil profiles was greater than that in the irrigated profiles. The CMB method severely underestimated recharge rates in the non‐irrigated areas when compared with the other methods, although the CMB method estimated recharge rates for the irrigated areas, that were similar to those from the other methods, ranging from 42 to 141 mm/year. The SWB method, constructed for a 15‐year period, provided insight into the recharge process being driven by winter rains rather than summer irrigation and indicated an average rate of 75 mm/year and 164 mm/year for the 1984 – 98 and 1996 – 98 periods, respectively. Assuming similar soil‐water holding capacity, these recharge rates applied to both irrigated and non‐irrigated areas. Use of the long period of record was important because it encompassed both drought and heavy rainfall years. Successful application of the SWB method, however, required considerable additional field measurements of orchard ET_c, soil‐water holding capacity and estimation of rainfall interception – runoff losses. Continuous soil moisture monitoring (NP) was necessary to identify both daily and seasonal seepage processes to corroborate the other recharge estimates. Measured recharge rates during the 1996 – 1998 period in both the orchards and non‐irrigated site averaged 180 mm/year. The pattern of soil profile drying during the summer irrigation season, followed by progressive wetting during the winter rainy season was observed in both irrigated and non‐irrigated soil profiles, confirming that groundwater recharge was rainfall driven and that micro‐irrigation did not ‘predispose’ the soil profile to excess rainfall recharge. The ability to make this recharge assessment, however, depended on making multiple field measurements associated with all three methods, suggesting that any one should not be used alone. Copyright © 2000 John Wiley & Sons, Ltd.     

Water quality in the Dickinson Bayou watershed (Texas, Gulf of Mexico) and health issues

The Dickinson Bayou watershed (near Houston, Texas, Gulf of Mexico) provides habitat for numerous coastally influenced communities of wildlife, including scores of birds and fish. Encroaching development and impervious surfaces are altering the habitat and degrading water quality. Herein we have defined the current health of the bayou using water quality data collected between 2000 and 2006. Elevated bacteria (fecal coliform, Escherichia coli and Enterococcus) and depressed dissolved oxygen concentrations (often <3 mg l⁻¹) are the two major impairments to this ecosystem. While nutrient ratios indicate primary productivity may be nitrogen limited, concerns of eutrophication persist because the bayou has a low intrinsic flushing rate. Consistent with this is the magnitude of algal blooms (ca. 100 μg chl l⁻¹) which often occur in spring/summer. The findings of this study will assist with the understanding of the influence of urban development on small watersheds.     

Hydrogeological characterization of groundwater storage and drainage in an alpine karst aquifer (the Kanin massif,Julian Alps)

The Kanin massif is an important trans‐boundary aquifer, which stretches between Slovenia and Italy. The groundwater is only partially exploited, mainly for water supply, but the aquifer exhibits great potential for future exploitation. Since no consistent regional overview of the hydrogeological functioning of the Kanin massif was available, the decision was made to perform a study of this area, using a pragmatic approach based on 3D geological and hydrogeological modelling. The so‐called KARSYS approach was applied, with the aim of characterizing the groundwater reserves within this karst massif and of locating the main drainage axes that carry groundwater from the recharge areas to the respective springs. Delineation of the catchment areas of the corresponding springs was carried out, and some new explanations were obtained, especially with regard to the Mo?nica spring, which is located in Slovenia and forms a potential source of drinking water. It was found that this spring's catchment area extends as far as the Italian ski resort of Sella Nevea. The conceptual model also provides a possible explanation about the underground drainage towards the Boka spring and waterfall, which has been a challenge for decades. This new explanation is based on the existence of a perched groundwater body that feeds the Boka spring via a system of conduits. Despite some limitations, the results, which consist of a visualization of the underground drainage and groundwater storage within the Kanin massif, can be used as a basis for planning the sustainable management of karst waters in the studied area. Copyright © 2014 John Wiley & Sons, Ltd.     

Distribution and sources of sedimentary organic matter in a tropical estuary,south west coast of India (Cochin estuary): A baseline study

Surface sediments samples were collected from 9 stations of the Cochin estuary during the monsoon, post-monsoon and pre-monsoon seasons and were analyzed for grain size, total organic carbon (OC), total nitrogen (TN) and stable isotopic ratios of carbon (δ¹³C) and nitrogen (δ¹⁵N) to identify major sources of organic matter in surface sediments. Sediment grain size is found to be the key factor influencing the organic matter accumulation in surface sediments. The δ¹³C values ranges from ?27.5‰ to ?21.7‰ in surface sediments with a gradual increase from inner part of the estuary to the seaward side that suggest an increasing contribution of marine autogenous organic matter towards the seaward side. The δ¹⁵N value varies between 3.1‰ and 6.7‰ and it exhibits complex spatial and seasonal distributions in the study area. It is found that the dynamic cycling of nitrogen through various biogeochemical and organic matter degradation processes modifies the OC/TN ratios and δ¹⁵N to a considerable degree. The fraction of terrestrial organic matter in the total organic matter pool ranges from 13% to 74% in the surface sediments as estimated by δ¹³C based two end member mixing model.     

Negligible soil erosion in a burned mountain watershed,Canadian Rockies: field and modelling investigations considering the role of duff

Increased soil erosion in immediate post‐wildfire years has been well documented in the literature, but many unanswered questions remain about the factors controlling erosional responses in different regional settings. The field site for the present study was located in a closed canopy, subalpine forest in Kootenay National Park, British Columbia that was subjected to a high‐intensity crown fire in the summer of 2003. Low soil erosion values were documented at the study site in the years immediately following the 2003 wildfire, with estimates ranging from approximately 10^‐1 up to 10⁰ t ha^‐1. Following the wildfire, notable duff coverage (the duff layer is the combined fermentation and humus soil organic layers) remained above the mineral soil. This finding supports earlier studies documenting only partial duff consumption by high‐intensity wildfires in the boreal forest of Canada. It is postulated that remnant duff coverage after many high‐intensity wildfires impacts the hydrological and soil erosional response to rainstorm events in post‐wildfire years. In particular, duff provides detention storage for infiltrating rainfall and, therefore, may inhibit the generation of overland flow. Furthermore, duff also provides a physical barrier to soil erosion. The Green–Ampt model of rainfall infiltration is employed to better assess how interactions between rainfall duration/intensity and soil/duff properties affect hydrological response and the generation of overland flow. Model results show that duff provides an effective zone for detention storage and that duff accommodates all rainfall intensities to which it was subjected without the occurrence of surface ponding. In addition, the penetration of the wetting front is relatively slow in duff due to its high porosity and water storage potential. Copyright © 2011 John Wiley & Sons, Ltd.     

Hydraulic subsurface measurements and hydrodynamic modelling as indicators for groundwater flow systems in the Rotondo granite,Central Alps (Switzerland)

Regional groundwater flow in high mountainous terrain is governed by a multitude of factors such as geology, topography, recharge conditions, structural elements such as fracturation and regional fault zones as well as man‐made underground structures. By means of a numerical groundwater flow model, we consider the impact of deep underground tunnels and of an idealized major fault zone on the groundwater flow systems within the fractured Rotondo granite. The position of the free groundwater table as response to the above subsurface structures and, in particular, with regard to the influence of spatial distributed groundwater recharge rates is addressed. The model results show significant unsaturated zones below the mountain ridges in the study area with a thickness of up to several hundred metres. The subsurface galleries are shown to have a strong effect on the head distribution in the model domain, causing locally a reversal of natural head gradients. With respect to the position of the catchment areas to the tunnel and the corresponding type of recharge source for the tunnel inflows (i.e. glaciers or recent precipitation), as well as water table elevation, the influence of spatial distributed recharge rates is compared to uniform recharge rates. Water table elevations below the well exposed high‐relief mountain ridges are observed to be more sensitive to changes in groundwater recharge rates and permeability than below ridges with less topographic relief. In the conceptual framework of the numerical simulations, the model fault zone has less influence on the groundwater table position, but more importantly acts as fast flow path for recharge from glaciated areas towards the subsurface galleries. This is in agreement with a previous study, where the imprint of glacial recharge was observed in the environmental isotope composition of groundwater sampled in the subsurface galleries. Copyright © 2012 John Wiley & Sons, Ltd.     

Water balance simulation for resource evaluation at the catchment scale: application to the Nema (Sudano-Sahelian zone,Senegal)

Abstract

The Nema is a river in western Senegal where only a minority of inhabitants has access to drinking water. Rainfall has been decreasing in this region since the 1960s. It is crucial to understand how this change affects groundwater recharge. The objective of this research is to determine the current proportions of groundwater recharge, runoff, subsurface runoff and evapotranspiration using a simulation approach. The Nash criterion and water balance error were used to evaluate the quality of the simulations. The following results were obtained: the Nash criterion was 0.73 for calibration (0.73 for validation), and the water balance error was ?0.35% and 0.005%, respectively, for the hydrological years 1995/96 and 1997/98. Evapotranspiration and groundwater recharge are the main processes involved.

Editor Z.W. Kundzewicz; Associate editor D. Hughes     

Parameter transferability under changing climate: case study with a land surface model in the Durance watershed,France

Abstract

In physically-based land surface models, the parameters can all be prescribed a priori but calibration can be used to enhance the realism of the simulations in well instrumented domains. In such a case, the transferability of calibrated parameters under non-stationary conditions needs to be addressed, especially in the context of climate change. To this end, we used the Catchment Land Surface Model (CLSM) in the Upper Durance watershed located in the French Alps, which experienced a significant increase in temperature over the last century. The CLSM is forced by a 50-year meteorological dataset of good quality. Four parameters of the CLSM (one related to snow processes and three to soil properties) are calibrated against discharge observations with a multi-objective algorithm. First, the robustness of the CLSM parameterizations is tested by the Differential Split Sample Test (DSST). The simulations show good performances over a wide range of retrospective climatic conditions, except when the parameters are calibrated over a period with a large contribution of snowmelt to annual mean discharge. Then, the use of a climate change scenario reveals that the parameterizations of soil moisture processes in the CLSM are responsible for an increasing dispersion among simulations when facing dry and warm conditions. However, the differences between the simulated changes of river discharge remain very small. This work shows that calibration conveys some uncertainties, but they are moderate in the studied case, and pertain to the most conceptual parameterizations of this physically-based model.     

Hydrological footprints of urban developments in the Lake Simcoe watershed,Canada: a combined paired‐catchment and change detection modelling approach

Urban sprawl and regional climate variability are major stresses on surface water resources in many places. The Lake Simcoe watershed (LSW) Ontario, Canada, is no exception. The LSW is predominantly agricultural but is experiencing rapid population growth because of its proximity to the Greater Toronto area. This has led to extensive land use changes that have impacted its water resources and altered run‐off patterns in some rivers draining to the lake. Here, we use a paired‐catchment approach, hydrological change detection modelling and remote sensing analysis of satellite images to evaluate the impacts of land use change on the hydrology of the LSW (1994 to 2008). Results show that urbanization increased up to 16% in Lovers Creek, the most urban‐impacted catchment. Annual run‐off from Lovers Creek increased from 239 to 442 mm/year in contrast to the reference catchment (Black River at Washago) where run‐off was relatively stable with an annual mean of 474 mm/year. Increased annual run‐off from Lovers Creek was not accompanied by an increase in annual precipitation. Discriminant function analysis suggests that early (1992–1997; pre‐major development) and late (2004–2009; fully urbanized) periods for Lovers Creek separated mainly based on model parameter sets related to run‐off flashiness and evapotranspiration. As a result, parameterization in either period cannot be used interchangeably to produce credible run‐off simulations in Lovers Creek because of greater scatter between the parameters in canonical space. Separation of early and late‐period parameter sets for the reference catchment was based on climate and snowmelt‐related processes. This suggests that regional climatic variability could be influencing hydrologic change in the reference catchment, whereas urbanization amplified the regional natural hydrologic changes in urbanizing catchments of the LSW. Copyright © 2014 John Wiley & Sons, Ltd.