Model-based analysis of nutrient retention and management for a lowland river

International audienceIn the context of the European Water Framework Directive options for improving the water quality of the lowland river Havel (Germany) were assessed. The lower section of this river is actually a polytrophic river-lake system suffering from high external nutrient loading and exhibiting significant in-river turnover. In order to gain a better understanding of present conditions and to allow integrated scenarios of nutrient management to be evaluated the catchment models SWIM and ArcEGMO-Urban were coupled with a simple, newly developed nutrient TRAnsport Model (TraM). Using the TraM model, the retention of nitrogen and phosphorus in a 55 km reach of the Lower Havel River was quantified and its temporal variation was analyzed. It was examined that about 30% of the external nitrogen input to the Lower Havel is retained within the surveyed river section. A comparison of simulation results generated with and without consideration of phosphorus retention/release revealed that summer TP concentrations are currently increased by 100?200% due to internal loading. Net phosphorus release rates of about 20 mg P m?2 d-1 in late summer were estimated for the Havel lakes. Scenario simulations with lowered external nutrient inputs revealed that persistent phosphorus limitation of primary production cannot be established within the next decade. It was shown that a further reduction in nitrogen concentrations requires emissions to be reduced in all inflows. Though the TraM model needs further extension it proved to be appropriate for conducting integrated catchment and river modeling

Use of soil moisture dynamics and patterns at different spatio-temporal scales for the investigation of subsurface flow processes

Spatial patterns as well as temporal dynamics of soil moisture have a major influence on runoff generation. The investigation of these dynamics and patterns can thus yield valuable information on hydrological processes, especially in data scarce or previously ungauged catchments. The combination of spatially scarce but temporally high resolution soil moisture profiles with episodic and thus temporally scarce moisture profiles at additional locations provides information on spatial as well as temporal patterns of soil moisture at the hillslope transect scale. This approach is better suited to difficult terrain (dense forest, steep slopes) than geophysical techniques and at the same time less cost-intensive than a high resolution grid of continuously measuring sensors. Rainfall simulation experiments with dye tracers while continuously monitoring soil moisture response allows for visualization of flow processes in the unsaturated zone at these locations. Data was analyzed at different spacio-temporal scales using various graphical methods, such as space-time colour maps (for the event and plot scale) and binary indicator maps (for the long-term and hillslope scale). Annual dynamics of soil moisture and decimeter-scale variability were also investigated. The proposed approach proved to be successful in the investigation of flow processes in the unsaturated zone and showed the importance of preferential flow in the Malalcahuello Catchment, a data-scarce catchment in the Andes of Southern Chile. Fast response times of stream flow indicate that preferential flow observed at the plot scale might also be of importance at the hillslope or catchment scale. Flow patterns were highly variable in space but persistent in time. The most likely explanation for preferential flow in this catchment is a combination of hydrophobicity, small scale heterogeneity in rainfall due to redistribution in the canopy and strong gradients in unsaturated conductivities leading to self-reinforcing flow paths

Use of soil moisture dynamics and patterns for the investigation of runoff generation processes with emphasis on preferential flow

International audienceSpatial patterns as well as temporal dynamics of soil moisture have a major influence on runoff generation. The investigation of these dynamics and patterns can thus yield valuable information on hydrological processes, especially in data scarce or previously ungauged catchments. The combination of spatially scarce but temporally high resolution soil moisture profiles with episodic and thus temporally scarce moisture profiles at additional locations provides information on spatial as well as temporal patterns of soil moisture at the hillslope transect scale. This approach is better suited to difficult terrain (dense forest, steep slopes) than geophysical techniques and at the same time less cost-intensive than a high resolution grid of continuously measuring sensors. Rainfall simulation experiments with dye tracers while continuously monitoring soil moisture response allows for visualization of flow processes in the unsaturated zone at these locations. Data was analyzed at different spacio-temporal scales using various graphical methods, such as space-time colour maps (for the event and plot scale) and indicator maps (for the long-term and hillslope scale). Annual dynamics of soil moisture and decimeter-scale variability were also investigated. The proposed approach proved to be successful in the investigation of flow processes in the unsaturated zone and showed the importance of preferential flow in the Malalcahuello Catchment, a data-scarce catchment in the Andes of Southern Chile. Fast response times of stream flow indicate that preferential flow observed at the plot scale might also be of importance at the hillslope or catchment scale. Flow patterns were highly variable in space but persistent in time. The most likely explanation for preferential flow in this catchment is a combination of hydrophobicity, small scale heterogeneity in rainfall due to redistribution in the canopy and strong gradients in unsaturated conductivities leading to self-reinforcing flow paths

Three perceptions of the evapotranspiration landscape: comparing spatial patterns from a distributed hydrological model, remotely sensed surface temperatures, and sub-basin water balances

A problem encountered by many distributed hydrological modelling studies is high simulation errors at interior gauges when the model is only globally calibrated at the outlet. We simulated river runoff in the Elbe River basin in central Europe (148 268 km2) with the semi-distributed eco-hydrological model SWIM (Soil and Water Integrated Model). While global parameter optimisation led to Nash–Sutcliffe efficiencies of 0.9 at the main outlet gauge, comparisons with measured runoff series at interior points revealed large deviations. Therefore, we compared three different strategies for deriving sub-basin evapotranspiration: (1) modelled by SWIM without any spatial calibration, (2) derived from remotely sensed surface temperatures, and (3) calculated from long-term precipitation and discharge data. The results show certain consistencies between the modelled and the remote sensing based evapotranspiration rates, but there seems to be no correlation between remote sensing and water balance based estimations. Subsequent analyses for single sub-basins identify amongst others input weather data and systematic error amplification in inter-gauge discharge calculations as sources of uncertainty. The results encourage careful utilisation of different data sources for enhancements in distributed hydrological modelling

Analysis and simulation of nutrient retention and management for a lowland river-lake system

International audienceIn the context of the European Water Framework Directive, we studied the possible impact of reduced emissions on phosphorus and nitrogen concentrations in a lowland river-lake system (Havel River, Germany). As a prerequisite, we quantified the retention of nutrients in the river from mass balances and deduced its seasonal variation. We detected that about 30% of the total nitrogen input is retained within the surveyed river section. In contrast, phosphorus release from sediments was shown to cause a considerable increase in present P concentrations. Average net phosphorus release rates of about 20 mg P m?2 d?1 in late summer were estimated for the Havel Lakes. Based on the observed patterns of N retention and P release we parametrized a newly developed water quality simulation program (TRAM), which allows alternative model approaches of different complexity to be implemented and tested. To account for the future trend of internal P loading, the phosphorus excess in lake sediments was estimated from core samples and included in the model as a state variable. For analyzing scenarios of reduced nutrient emissisions, the water quality simulation program was linked to mesoscale hydrological catchment models for the first time. From scenario simulations we conclude that internal P loading is likely to counteract efforts of emission control for decades. Even by significant reductions in external P loads, a persistent phosphorus limitation of primary production can hardly be established in the analyzed time frame of 13 years. Though in the short run a continued reduction in nitrogen loads appears to be the more promising approach of eutrophication management, we recommend enhanced efforts to diminish both N and P emissions

Plot and field scale soil moisture dynamics and subsurface wetness control on runoff generation in a headwater in the Ore Mountains

This study presents an application of an innovative sampling strategy to assess soil moisture dynamics in a headwater of the Weißeritz in the German eastern Ore Mountains. A grassland site and a forested site were instrumented with two Spatial TDR clusters (STDR) that consist of 39 and 32 coated TDR probes of 60 cm length. Distributed time series of vertically averaged soil moisture data from both sites/ensembles were analyzed by statistical and geostatistical methods. Spatial variability and the spatial mean at the forested site were larger than at the grassland site. Furthermore, clustering of TDR probes in combination with long-term monitoring allowed identification of average spatial covariance structures at the small field scale for different wetness states. The correlation length of soil water content as well as the sill to nugget ratio at the grassland site increased with increasing average wetness and but, in contrast, were constant at the forested site. As soil properties at both the forested and grassland sites are extremely variable, this suggests that the correlation structure at the forested site is dominated by the pattern of throughfall and interception. We also found a very strong correlation between antecedent soil moisture at the forested site and runoff coefficients of rainfall-runoff events observed at gauge Rehefeld. Antecedent soil moisture at the forest site explains 92% of the variability in the runoff coefficients. By combining these results with a recession analysis we derived a first conceptual model of the dominant runoff mechanisms operating in this catchment. Finally, we employed a physically based hydrological model to shed light on the controls of soil- and plant morphological parameters on soil average soil moisture at the forested site and the grassland site, respectively. A homogeneous soil setup allowed, after fine tuning of plant morphological parameters, most of the time unbiased predictions of the observed average soil conditions observed at both field sites. We conclude that the proposed sampling strategy of clustering TDR probes is suitable to assess unbiased average soil moisture dynamics in critical functional units, in this case the forested site, which is a much better predictor for event scale runoff formation than pre-event discharge. Long term monitoring of such critical landscape elements could maybe yield valuable information for flood warning in headwaters. We thus think that STDR provides a good intersect of the advantages of permanent sampling and spatially highly resolved soil moisture sampling using mobile rods

The role of spatial variability of soil moisture for modelling surface runoff generation at the small catchment scale

International audienceA Bayesian approach is described for dealing with the problem of infilling and generating stochastic flow sequences using rainfall data to guide the flow generation process, and including bounded (censored) observed flow and rainfall data to provide additional information. Solutions are obtained using a Gibbs sampling procedure. Particular problems discussed include developing new procedures for fitting transformations when bounded values are available, coping with additional information in the form of values, or bounds, for totals of flows across several sites, and developing relationships between annual flow and rainfall data. Examples are shown of both infilled values of unknown past river flows, with assessment of uncertainty, and realisations of flows representative of what might occur in the future. Several procedures for validating the model output are described and the central estimates of flows, taken as a surrogate for historical observed flows, are compared with long term regional flow and rainfall data

Application of two phosphorus models with different complexities in a mesoscale river catchment

The water balance and phosphorus inputs of surface waters of the Weiße Elster catchment, Germany, have been quantified using the models GROWA/MEPhos and SWAT. A comparison of the model results shows small differences in the mean long-term total runoff for the entire study area. All relevant pathways of phosphorus transport were considered in MEPhos with phosphorus inputs resulting to about 65% from point sources. SWAT focuses on agricultural areas and estimates a phosphorus input of about 60% through erosion. The mean annual phosphorus input from erosion calculated with SWAT is six times higher than the estimation with MEPhos due to the differing model concepts. This shows the uncertainty contributed by the modelling description of phosphorus pathways