A multi‐storage groundwater concept for the SWAT model to emphasize nonlinear groundwater dynamics in lowland catchments

Hydrological models are useful tools to analyze present and future conditions of water quantity and quality. The integrated modelling of water and nutrients needs an adequate representation of the different discharge components. In common with many lowlands, groundwater contribution to the discharge in the North German lowlands is a key factor for a reasonable representation of the water balance, especially in low flow periods. Several studies revealed that the widely used Soil and Water Assessment Tool (SWAT) model performs poorly for low flow periods. This paper deals with the extension of the groundwater module of the SWAT model to enhance low flow representation. The current two‐storage concept of SWAT was further developed to a three‐storage concept. This was realized due to modification of the groundwater module by splitting the active groundwater storage into a fast and a slow contributing aquifer. The results of this study show that the groundwater module with three storages leads to a good prediction of the overall discharge especially for the recession limbs and the low flow periods. The improved performance is reflected in the signature measures for the mid‐segment (percent bias ?2.4% vs ?15.9%) and the low segment (percent bias 14.8% vs 46.8%) of the flow duration curve. The three‐storage groundwater module is more process oriented than the original version due to the introduction of a fast and a slow groundwater flow component. The three‐storage version includes a modular approach, because groundwater storages can be activated or deactivated independently for subbasin and hydrological response unit level. Copyright © 2013 John Wiley & Sons, Ltd.     

Temporal variation of suspended sediment transport in the Koga catchment,North Western Ethiopia and environmental implications

Event sediment transport and yield were studied for 45 events in the upstream part of the 260 km² agricultural Koga catchment that drains to an irrigation reservoir. Discharge and turbidity data were collected over a period of more than a year, accompanied by grab sampling. Turbidity was very well correlated with the sediment concentrations from the samples (r = 0.99), which allowed us to estimate the temporal patterns of sediment concentrations within events. The hysteresis patterns between discharge and sediment concentrations were analysed to provide insight into the different sediment sources. Anticlockwise patterns are the dominant hysteresis patterns in the area, suggesting smaller contributions of suspended sediment from the river channels than from the hillslopes and agricultural areas. Complicated types of hysteresis patterns were mostly observed for long events with multiple peaks. For a given discharge, sediment yields in August and September, when the catchment was almost completely covered with vegetation, were much smaller than during the rest of the rainy season. The hysteresis patterns and timing suggest that the sediment availability from the agricultural areas and hillslopes affects sediment yields more strongly than does peak discharge. Two distinct types of sediment rating curves were observed for the season when the agricultural land was covered with vegetation and when it was not, indicating the dominating contribution of land use/cover to sediment yields in the catchment. The rate of suspended sediment transport in the area was estimated as 25.6 t year^?1 ha^?1. Copyright © 2013 John Wiley & Sons, Ltd.     

Analysis of uncertainties in the hydrological response of a model‐based climate change impact assessment in a subcatchment of the Spree River,Germany

Climate change impact assessments form the basis for the development of suitable climate change adaptation strategies. For this purpose, ensembles consisting of stepwise coupled models are generally used [emission scenario → global circulation model → downscaling approach (DA) → bias correction → impact model (hydrological model)], in which every item is affected by considerable uncertainty. The aim of the current study is (1) to analyse the uncertainty related to the choice of the DA as well as the hydrological model and its parameterization and (2) to evaluate the vulnerability of the studied catchment, a subcatchment of the highly anthropogenically impacted Spree River catchment, to hydrological change. Four different DAs are used to drive four different model configurations of two conceptually different hydrological models (Water Balance Simulation Model developed at ETH Zürich and HBV‐light). In total, 452 simulations are carried out. The results show that all simulations compute an increase in air temperature and potential evapotranspiration. For precipitation, runoff and actual evapotranspiration, opposing trends are computed depending on the DA used to drive the hydrological models. Overall, the largest source of uncertainty can be attributed to the choice of the DA, especially regarding whether it is statistical or dynamical. The choice of the hydrological model and its parameterization is of less importance when long‐term mean annual changes are compared. The large bandwidth at the end of the modelling chain may exacerbate the formulation of suitable climate change adaption strategies on the regional scale. Copyright © 2013 John Wiley & Sons, Ltd.     

Coupled groundwater flow and heat transport simulation for estimating transient aquifer–stream exchange at the lowland River Spree (Germany)

Subsurface flow and heat transport near Freienbrink, NE Germany, was simulated in order to study groundwater–surface water exchange between a floodplains aquifer and a section of the lowland River Spree and an adjacent oxbow. Groundwater exfiltration was the dominant process, and only fast surface water level rises resulted in temporary infiltration into the aquifer. The main groundwater flow paths are identified based on a 3D groundwater flow model. To estimate mass fluxes across the aquifer–surface water interfaces, a 2D flow and heat transport modelling approach along a transect of 12 piezometers was performed. Results of steady‐state and transient water level simulations show an overall high accuracy with a Spearman coefficient ρ = 0.9996 and root mean square error (RMSE) = 0.008 m. Based on small groundwater flow velocities of about 10^?7 to 10^?6 ms^?1, mean groundwater exfiltration rates of 233 l m^?2 d^?1 are calculated. Short periods of surface water infiltration into the aquifer do not exceed 10 days, and the infiltration rates are in the same range. The heat transport was modelled with slightly less accuracy (ρ = 0.8359 and RMSE = 0.34 °C). In contrast to the predominant groundwater exfiltration, surface water temperatures determine the calculated temperatures in the upper aquifer below both surface water bodies down to 10 m during the whole simulation period. These findings emphasize prevailing of heat conduction over advection in the upper aquifer zones, which seems to be typical for lowland streams with sandy aquifer materials and low hydraulic gradients. Moreover, this study shows the potential of coupled numerical flow and heat transport modelling to understand groundwater–surface water exchange processes in detail. Copyright © 2013 John Wiley & Sons, Ltd.     

A Global Terrestrial Reference Frame from simulated VLBI and SLR data in view of GGOS

In this study, we assess the impact of two combination strategies, namely local ties (LT) and global ties (GT), on the datum realization of Global Terrestrial Reference Frames in view of the Global Geodetic Observing System requiring 1 mm-accuracy. Simulated Very Long Baseline Interferometry (VLBI) and Satellite Laser Ranging (SLR) data over a 7 year time span was used. The LT results show that the geodetic datum can be best transferred if the precision of the LT is at least 1 mm. Investigating different numbers of LT, the lack of co-located sites on the southern hemisphere is evidenced by differences of 9 mm in translation and rotation compared to the solution using all available LT. For the GT, the combination applying all Earth rotation parameters (ERP), such as pole coordinates and UT1-UTC, indicates that the rotation around the Z axis cannot be adequately transferred from VLBI to SLR within the combination. Applying exclusively the pole coordinates as GT, we show that the datum can be transferred with mm-accuracy within the combination. Furthermore, adding artificial stations in Tahiti and Nigeria to the current VLBI network results in an improvement in station positions by 13 and 12%, respectively, and in ERP by 17 and 11%, respectively. Extending to every day VLBI observations leads to 65% better ERP estimates compared to usual twice-weekly VLBI observations.     

Comparison of approaches of geographic partitioning for data anonymization

Journal of Geographical Systems - Given the large volumes of detailed data now being collected, there is a high demand for the release of this data for research purposes. In particular,...     

Stable isotopes and trace elements in modern ostracod shells: implications for reconstructing past environments on the Tibetan Plateau,China

Stable isotopes and trace elements in ostracod shells have been used widely in paleolimnological investigations of past lake hydrochemistry and climate because they provide insights into past water balance and solute evolution of lakes. Regional differences in lake characteristics and species-specific element fractionation, however, do not permit generalization of results from other regions or ostracod species to the southern Tibetan Plateau, in part because most common taxa from the southern Tibetan Plateau are endemic to the area. This study evaluated relations between present-day environmental conditions and the geochemical composition of modern ostracod shells from the southern Tibetan Plateau, to assess the suitability of using shell chemistry to infer hydrological conditions. We studied nine lakes and their catchments, located along a west–east transect in the south-central part of the Tibetan Plateau. Stable oxygen and carbon isotope values and trace element concentrations in recent shells from the four most abundant ostracod species (Leucocytherella sinensis, ?Leucocythere dorsotuberosa, Limnocythere inopinata, Tonnacypris gyirongensis) were measured, together with hydrochemical properties of host waters at the time of sampling. Results revealed significant between-species differences in stable isotope fractionation and trace element incorporation into shell calcite. Stable oxygen and carbon isotope values of ostracod shells were correlated significantly with the stable isotope composition of the respective water body \( \left( {\updelta^{18} {\text{O}}_{{{\text{H}}_{ 2} {\text{O}}}} \,{\text{and }}\updelta^{13} {\text{C}}_{{{\text{H}}_{ 2} {\text{O}}}} } \right) \), reflecting salinity and productivity, respectively. Offsets between δ¹⁸O_shell and δ¹³C_shell and inorganic calcite, the latter representing isotopic equilibrium, suggest shell formation of T. gyirongensis during spring melt. L. sinensis reproduces throughout the monsoon season until September and shows several consecutive generations, and L. inopinata molts to the adult stage after the monsoon season in August/September. The influence of pore water δ¹³C was displayed by L. inopinata, suggesting shell calcification within the sediment. Mg/Ca_shell is primarily influenced by water Mg/Ca ratios and salinity and confirms the use of this shell ratio as a proxy for precipitation-evaporation balance and lake level. In addition, Sr/Ca and Ba/Ca can be used to infer changes in salinity, at least in closed-basin lakes with calcite saturation. Observed effects of water Sr/Ca and salinity on Sr/Ca incorporation are biased by the presence of aragonite precipitation in the lakes, which removes bioavailable Sr from the host water, resulting in low Sr/Ca_shell values. Changes in carbonate mineralogy affect the bioavailability of trace elements, a process that should be considered in paleoclimate reconstructions. Oxygen isotopes and Mg/Ca_shell ratios were unaffected by water temperature. Positive correlations among Fe/Ca, Mn/Ca and U/Ca in ostracod shells, and their negative correlation with δ¹³C, which reflects organic matter decay, show the potential to infer changes in redox conditions that can be used to reconstruct past oxygen supply to bottom waters and thus past water-circulation changes within lakes. The intensity of microbial activity, associated with organic matter decomposition, can be inferred from U/Ca ratios in ostracod shells. These findings highlight the value of fossil ostracod records in lake deposits for inferring paleoenvironmental conditions on the southern Tibetan Plateau.     

Spatial variability in the chemical composition of the snowcover at high alpine sites

Summary For a comparison of snow chemistry data from different glaciers or snow fields it is important to have informations about the spatial representativeness of the data from each of the individual sites. In order to assess the representativeness of concentration data of major ions (volume weighted means of the winter accumulation) from glacier fields we investigated the variability in the average concentration of major ions from parallel samples within one snow pit and from different pits within one glacier field. The variabilities in the average concentrations of NO ₃ ^- , SO ₄ ^2- and NH ₄ ⁺ for three parallel profiles within one snow pit at Goldbergkees (Austria) were 1.2, 3.3 and 2.0% (coefficient of variation). Cl^– and Na⁺ showed larger variations (6.6 and 56.6%) possibly originating from contaminations. The variability of average concentration data from different snow pits within one glacier field was studied at La Grave (France) and at Goldbergkees (Austria). At La Grave 3 pits and at Goldbergkees 4 pits at distance of several hundred meters were sampled. Variabilities for SO ₄ ^2- and NO ₃ ^- were quite similar for the two sites (17% for both ions at La Grave, 12% and 15% at Goldbergkees). Whereas variabilities for Na⁺, NH ₄ ⁺ , Mg²⁺, Ca²⁺ and Cl^– were quite low at La Grave ( 12% and 27% for Cl^–), higher values were obtained at Goldbergkees for these ions (17–56%). Likely reasons for the higher variability observed at Goldbergkees are a) higher spatial variability of crustal aerosol species (Mg²⁺, Ca²⁺), b) problems with the detection limit of the analytical method (Ca²⁺), c) contaminations (Na⁺, Cl^–).With 4 Figures     

Petrographic and trace-element distribution studies on the dolomite-calcite in the regional metamorphic marble of the Griesscharte,Tyrol, Austria/Italy

Petrographic investigations on the light yellow marble outcrops in the Griesscharte, Tyrol, Austria/Italy revealed them to be dolomite. The white to blue-grey veins in these marbles were those of calcite. Inside these calcites and along the calcite-dolomite phase boundaries minute and microscopically unidentifiable crystallites (? 10 μ) were observed. These crystallites were recognised by X-ray diffractometry as dolomite and their form and distribution were established by scanning electron microscopy. Trace-element distributions on the calcite-dolomite neighbours were studied by neutron-activation method to elucidate the formation mechanism of these calcites in the marble.Genetic relationships of the dolomite crystallites, the calcite host and the surrounding dolomite marble have been discussed in the light of the petrographic and chemical studies. These findings substantiate a deposition of a new formed vein of calcite from the solution circulating in the fissures of the dolomite marble during prograde metamorphism rather than a simple deposition of mobilized sedimentary calcite. The dolomite crystallites in the calcite seem to have been formed post-tectonically, partly by the nucleation in the calcite cleavage planes during Mg-metasomatism and partly from the exsolutions during retrograde metamorphism.