scholarly journals Application of pore water stable isotope method to characterise a wetland system

2018 ◽  
Author(s):  
Katarina David ◽  
Wendy Timms ◽  
Cath E. Hughes ◽  
Jagoda Crawford ◽  
Dayna McGeeney

Abstract. Three naturally intact wetland systems (swamps) were characterized based on sediment cores, analysis of surface water, groundwater and porewater stable isotopes. These swamps are classified as temperate highland peat swamps on sandstone (THPSS) and in Australia they are listed as threatened ecological communities. This study is the first application of the stable isotope direct vapour equilibration method in a wetland, enabling quantification of the contributions of evaporation, rainfall and groundwater to swamp water balance. This technique enables understanding of the depth of evaporative losses and the relative importance of groundwater flow within the swamp environment without the need for intrusive piezometer installation at multiple locations and depths. Additional advantages of the stable isotope direct vapour equilibration technique include detailed spatial and vertical depth profiles of δ18O and δ2H, with good accuracy comparable to the porewater compression technique. Depletion of δ18O and δ2H in porewater with increasing depth (to around 40–60 cm depth) was observed in two swamps, but remained uniform with depth in the third swamp. Within the upper surficial zone, the measurements respond to seasonal trends and are subject to evaporation in the capillary zone. Below this depth the pore water δ18O and δ2H signature approaches that of groundwater indicating lateral groundwater contribution. Significant differences were found in stable pore water isotopes for samples collected after dry weather period compared to wet periods where recharge of depleted rainfall was apparent. The organic rich soil in the upper 40–60 cm retains significant saturation following precipitation events and maintains moisture necessary for ecosystem functioning. An important finding for wetland and ecosystem response to changing groundwater conditions (and potential ground movement) are the observations that basal sands underlay the swamps, allowing relatively rapid drainage at the base of the swamp and interaction with lateral groundwater contribution. Based on the novel stable isotope direct vapour equilibration analysis of swamp sediment, our study identified the following important processes: rapid infiltration of rainfall to the water table with longer retention of moisture in the upper 40–60 cm and lateral groundwater flow contribution at the base. This study also found, that evaporation estimated using stable isotope direct vapour equilibration method is more realistic compared to reference evapotranspiration (ET). Importantly, if swamp discharge data were available in combination with pore water isotope profiles, an appropriate transpiration could be determined for these swamps. Based on the results, the groundwater contribution to the swamp is a significant component of the water balance during dry period. Our methods could complement other monitoring studies and numerical water balance models to improve prediction of the hydrological response of the swamp to changes in water conditions due to natural or anthropogenic influences.

2018 ◽  
Vol 22 (11) ◽  
pp. 6023-6041 ◽  
Author(s):  
Katarina David ◽  
Wendy Timms ◽  
Catherine E. Hughes ◽  
Jagoda Crawford ◽  
Dayna McGeeney

Abstract. Three naturally intact wetland systems (swamps) were characterised based on sediment cores, analysis of surface water, swamp groundwater, regional groundwater and pore water stable isotopes. These swamps are classified as temperate highland peat swamps on sandstone (THPSS) and in Australia they are listed as threatened endangered ecological communities under state and federal legislation. This study applies the stable isotope direct vapour equilibration method in a wetland, aiming at quantification of the contributions of evaporation, rainfall and groundwater to swamp water balance. This technique potentially enables understanding of the depth of evaporative losses and the relative importance of groundwater flow within the swamp environment without the need for intrusive piezometer installation at multiple locations and depths. Additional advantages of the stable isotope direct vapour equilibration technique include detailed spatial and vertical depth profiles of δ18O and δ2H, with good accuracy comparable to other physical and chemical extraction methods. Depletion of δ18O and δ2H in pore water with increasing depth (to around 40–60 cm depth) was observed in two swamps but remained uniform with depth in the third swamp. Within the upper surficial zone, the measurements respond to seasonal trends and are subject to evaporation in the capillary zone. Below this depth the pore water δ18O and δ2H signature approaches that of regional groundwater, indicating lateral groundwater contribution. Significant differences were found in stable pore water isotope samples collected after the dry weather period compared to wet periods where recharge of depleted rainfall (with low δ18O and δ2H values) was apparent. The organic-rich soil in the upper 40 to 60 cm retains significant saturation following precipitation events and maintains moisture necessary for ecosystem functioning. An important finding for wetland and ecosystem response to changing swamp groundwater conditions (and potential ground movement) is that basal sands are observed to underlay these swamps, allowing relatively rapid drainage at the base of the swamp and lateral groundwater contribution. Based on the novel stable isotope direct vapour equilibration analysis of swamp sediment, our study identified the following important processes: rapid infiltration of rainfall to the water table with longer retention of moisture in the upper 40–60 cm and lateral groundwater flow contribution at the base. This study also found that evaporation estimated using the stable isotope direct vapour equilibration method is more realistic compared to reference evapotranspiration (ET). Importantly, if swamp discharge data were available in combination with pore water isotope profiles, an appropriate transpiration rate could be determined for these swamps. Based on the results, the groundwater contribution to the swamp is a significant and perhaps dominant component of the water balance. Our methods could complement other monitoring studies and numerical water balance models to improve prediction of the hydrological response of the swamp to changes in water conditions due to natural or anthropogenic influences.


Water ◽  
2021 ◽  
Vol 13 (4) ◽  
pp. 563
Author(s):  
Wiebe Förster ◽  
Jan C. Scholten ◽  
Michael Schubert ◽  
Kay Knoeller ◽  
Nikolaus Classen ◽  
...  

The eutrophic Lake Eichbaumsee, a ~1 km long and 280 m wide (maximum water depth 16 m) dredging lake southeast of Hamburg (Germany), has been treated for water quality improvements using various techniques (i.e., aeration plants, removal of dissolved phosphorous by aluminum phosphorous precipitation, and by Bentophos® (Phoslock Environmental Technologies, Sydney, Australia), adsorption) during the past ~15 years. Despite these treatments, no long-term improvement of the water quality has been observed and the lake water phosphorous content has continued to increase by e.g., ~670 kg phosphorous between autumn 2014 and autumn 2019. As no creeks or rivers drain into the lake and hydrological groundwater models do not suggest any major groundwater discharge into the lake, sources of phosphorous (and other nutrients) are unknown. We investigated the phosphorous fluxes from sediment pore water and from groundwater in the water body of the lake. Sediment pore water was extracted from sediment cores recovered by divers in August 2018 and February 2019. Diffusive phosphorous fluxes from pore water were calculated based on phosphorus gradients. Stable water isotopes (δ2H, δ18O) were measured in the lake water, in interstitial waters in the banks surrounding the lake, in the Elbe River, and in three groundwater wells close to the lake. Stable isotope (δ2H, δ18O) water mass balance models were used to compute water inflow/outflow to/from the lake. Our results revealed pore-water borne phosphorous fluxes between 0.2 mg/m2/d and 1.9 mg/m2/d. Assuming that the measured phosphorous fluxes are temporarily and spatially representative for the whole lake, about 11 kg/a to 110 kg/a of phosphorous is released from sediments. This amount is lower than the observed lake water phosphorous increase of ~344 kg between April 2018 and November 2018. Water stable isotope (δ2H, δ18O) compositions indicate a water exchange between an aquifer and the lake water. Based on stable isotope mass balances we estimated an inflow of phosphorous from the aquifer to the lake of between ~150 kg/a and ~390 kg/a. This result suggests that groundwater-borne phosphorous is a significant phosphorous source for the Eichbaumsee and highlights the importance of groundwater for lake water phosphorous balances.


2015 ◽  
Vol 19 (6) ◽  
pp. 2617-2635 ◽  
Author(s):  
M. Sprenger ◽  
T. H. M. Volkmann ◽  
T. Blume ◽  
M. Weiler

Abstract. Determining the soil hydraulic properties is a prerequisite to physically model transient water flow and solute transport in the vadose zone. Estimating these properties by inverse modelling techniques has become more common within the last 2 decades. While these inverse approaches usually fit simulations to hydrometric data, we expanded the methodology by using independent information about the stable isotope composition of the soil pore water depth profile as a single or additional optimization target. To demonstrate the potential and limits of this approach, we compared the results of three inverse modelling strategies where the fitting targets were (a) pore water isotope concentrations, (b) a combination of pore water isotope concentrations and soil moisture time series, and (c) a two-step approach using first soil moisture data to determine water flow parameters and then the pore water stable isotope concentrations to estimate the solute transport parameters. The analyses were conducted at three study sites with different soil properties and vegetation. The transient unsaturated water flow was simulated by solving the Richards equation numerically with the finite-element code of HYDRUS-1D. The transport of deuterium was simulated with the advection-dispersion equation, and a modified version of HYDRUS was used, allowing deuterium loss during evaporation. The Mualem–van Genuchten and the longitudinal dispersivity parameters were determined for two major soil horizons at each site. The results show that approach (a), using only the pore water isotope content, cannot substitute hydrometric information to derive parameter sets that reflect the observed soil moisture dynamics but gives comparable results when the parameter space is constrained by pedotransfer functions. Approaches (b) and (c), using both the isotope profiles and the soil moisture time series, resulted in good simulation results with regard to the Kling–Gupta efficiency and good parameter identifiability. However, approach (b) has the advantage that it considers the isotope data not only for the solute transport parameters but also for water flow and root water uptake, and thus increases parameter realism. Approaches (b) and (c) both outcompeted simulations run with parameters derived from pedotransfer functions, which did not result in an acceptable representation of the soil moisture dynamics and pore water stable isotope composition. Overall, parameters based on this new approach that includes isotope data lead to similar model performances regarding the water balance and soil moisture dynamics and better parameter identifiability than the conventional inverse model approaches limited to hydrometric fitting targets. If only data from isotope profiles in combination with textural information is available, the results are still satisfactory. This method has the additional advantage that it will not only allow us to estimate water balance and response times but also site-specific time variant transit times or solute breakthrough within the soil profile.


2010 ◽  
Vol 7 (5) ◽  
pp. 7809-7838 ◽  
Author(s):  
M. Larocque ◽  
V. Fortin ◽  
M. C. Pharand ◽  
C. Rivard

Abstract. Groundwater contribution to river flows, generally called base flows, often accounts for a significant proportion of total flow rate, especially during the dry season. The objective of this work is to test simple approaches requiring limited data to understand groundwater contribution to river flows. The Noire river basin in southern Quebec is used as a case study. A lumped conceptual hydrological model (the MOHYSE model), a groundwater flow model (MODFLOW) and hydrograph separation are used to provide estimates of base flow for the study area. Results show that the methods are complementary. Hydrograph separation and the MOHYSE surface flow model provide similar annual estimates for the groundwater contribution to river flow, but monthly base flows can vary significantly between the two methods. Both methods have the advantage of being easily implemented. However, the distinction between aquifer contribution and shallow subsurface contribution to base flow can only be made with a groundwater flow model. The aquifer renewal rate estimated with the MODFLOW model for the Noire River is 30% of the recharge estimated from base flow values. This is a significantly difference which can be crucial for regional-scale water management.


2021 ◽  
Author(s):  
Ivan Vorobevskii ◽  
Rico Kronenberg

<p>‘Just drop a catchment and receive reasonable model output’ – still stays as motto and main idea of the ‘Global BROOK90’ project. The open-source R-package is build-up on global land cover, soil, topographical, meteorological datasets and the lumped hydrological model as a core to simulate water balance components on HRU scale all over the world in an automatic mode. First introduced in EGU2020 and followed by GitHub code release including an publication of methodology with few examples we want to continue with the insights on the current state and highlight the future steps of the project.</p><p>A global validation of discharge and evapotranspiration components of the model showed promising results. We used 190 small (median size of 64 km<sup>2</sup>) catchments and FLUXNET data which represent a wide range of relief, vegetation and soil types within various climate zones. The model performance was evaluated with NSE, KGE, KGESS and MAE. In more than 75 % of the cases the framework performed better than the mean of the observed discharge. On a temporal scale the performance is significantly better on a monthly vs daily scale. Cluster analysis revealed that some of the site characteristics have a significant influence on the performance. Additionally, it was found that Global BROOK90 outperforms GloFAS ERA5 discharge reanalysis (for the category with smallest catchments).</p><p>A cross-combination of three different BROOK90 setups and three forcing datasets was set up to reveal uncertainties of the Global BROOK90 package using a small catchment in Germany as a case study. Going from local to regional and finally global scale we compared mixtures of model parameterization schemes (original calibrated BROOK90, EXTRUSO and Global BROOK90) and meteorological datasets (local gauges, RaKlida and ERA5). Besides high model performances for a local dataset plus a calibrated model and weaker results for ERA5 and the Global BROOK90, it was found that the ERA5 dataset is still able to provide good results when combined with a regional and local parameterization. On the other side, the combination of a global parameterization with local and regional forcings gives still adequate, but much worse results. Furthermore, a hydrograph separation revealed that the Global BROOK90 parameterization as well as ERA5 discharge data perform weaker especially within low flow periods.</p><p>Currently, some new features are added to the original package. First, with the recent release of the ERA5 extension, historical simulations with the package now are expanded to 1950-2021 period. Additionally, an alternative climate reanalysis dataset is included in the framework (Merra-2, 0.5x0.625-degree spatial resolution, starting from 1980). A preliminary validation shows insignificant differences between both meteorological datasets with respect to the discharge based model performance.</p><p>Further upgrades of the framework will include the following core milestones: recognition of forecast and climate projections and parameter optimization features. In the nearest future we plan to utilize full power of the Climate Data Store for easy access to seasonal forecasts (i.e. ECMWF, DWD, NCEP) as well as climate projections (CMIP5) to extend the package’s scope to predict near and far future water balance components.</p>


Water ◽  
2020 ◽  
Vol 12 (6) ◽  
pp. 1775
Author(s):  
Jozsef Szilagyi

Watershed-scale annual evapotranspiration (ET) is routinely estimated by a simplified water balance as the difference in catchment precipitation (P) and stream discharge (Q). With recent developments in ET estimation by the calibration-free generalized complementary relationship, the water balance equation is employed to estimate watershed/basin P at an annual scale as ET + Q on the United States (US) Geological Survey’s Hydrologic Unit Code (HUC) 2- and 6-level watersheds over the 1979–2015 period. On the HUC2 level, mean annual PRISM P was estimated with a correlation coefficient (R) of 0.99, relative bias (RB) of zero, root-mean-squared-error (RMSE) of 54 mm yr−1, ratio of standard deviations (RS) of 1.08, and Nash–Sutcliffe efficiency (NSE) of 0.98. On the HUC6 level, R, RS, and NSE hardly changed, RB remained zero, while RMSE increased to 90 mm yr−1. Even the long-term linear trend values were found to be fairly consistent between observed and estimated values with R = 0.97 (0.81), RMSE = 0.63 (1.63) mm yr−1, RS = 0.99 (1.05), NSE = 0.92 (0.59) on the HUC2 and HUC6 (in parentheses) levels. This calibration-free water-balance method demonstrates that annual watershed precipitation can be estimated with an acceptable accuracy from standard atmospheric/radiation and stream discharge data.


Sign in / Sign up

Export Citation Format

Share Document