Sven Frei

Representing effects of micro-topography on runoff generation and sub-surface flow patterns by using superficial rill/depression storage height variations

An adequate representation of micro-topography in spatially explicit, physically based models can be crucial in modeling runoff generation, surface/subsurface flow interactions or subsurface flow patterns in hydrological systems with pronounced micro-topography. However, representation of micro-topography in numerical models usually requires high grid resolutions to capture relevant small scale variations in topography at the range of centimeters to meters. High grid resolutions usually result in longer simulation times, especially if fully integrated model approaches are used where the governing partial differential equations for surface and subsurface flow are solved simultaneously. This often restricts the implementation of micro-topography to plot scale models where the overall model domain is small to minimize computational cost resulting from a high grid resolution. In this study an approach is presented where a highly resolved digital elevation model (DEM) for a hummocky topography in a plot scale wetland model (10 m x 21 m x 2 m), is represented by spatially distributed rill/depression storage zones in a numerical model with a planar surface. By replacing the explicit micro-topography with spatially distributed rill/depression storage zones, important effects of micro-topography on surface flow generation and subsurface transport characteristics (e.g. residence time distributions) are being preserved, while at the same time the number of computational nodes is reduced significantly. We demonstrate that the rill/depression storage concept, which has been used for some time to represent time delays in the generation of surface runoff, can also be used to mimic subsurface flow patterns caused by micro-topography. Results further indicate that the rill/depression storage concept is an efficient tool to represent micro-topography in plot scale models because model computation times drop significantly. As important aspects of surface and subsurface flows induced by micro-topography can be mimicked adequately by applying the rill/depression storage concept on a coarser grid, it may also be a useful tool to represent micro-topography in numerical flow models beyond the plot scale.

On the value of surface saturated area dynamics mapped with thermal infrared imagery for modeling the hillslope-riparian-stream continuum

The highly dynamic processes within a hillslope-riparian-stream (HRS) continuum are known to affect streamflow generation, but are yet not fully understood. Within this study, we simulated a headwater HRS continuum in western Luxembourg with an integrated hydrologic surface subsurface model (HydroGeoSphere). The model was set up with thorough consideration of catchment-specific attributes and we performed a multi criteria model evaluation (4 years) with special focus on the temporally varying spatial patterns of surface saturation. We used a portable thermal infrared (TIR) camera to map surface saturation with a high spatial resolution and collected 20 panoramic snapshots of the riparian zone (approx. 10 m x 20 m) under different hydrologic conditions. Qualitative and quantitative comparison of the processed TIR panoramas and the corresponding model output panoramas revealed a good agreement between spatiotemporal dynamic model and field surface saturation patterns. A double logarithmic linear relationship between surface saturation extent and discharge was similar for modeled and observed data. This provided confidence in the capability of an integrated hydrologic surface subsurface model to represent temporal and spatial water flux dynamics at small (HRS continuum) scales. However, model scenarios with different parameterizations of the riparian zone showed that discharge and surface saturation were controlled by different parameters and hardly influenced each other. Surface saturation only affected very fast runoff responses with a small volumetric contribution to stream discharge, indicating that the dynamic surface saturation in the riparian zone does not necessarily imply a major control on runoff generation. This article is protected by copyright. All rights reserved.

FINIFLUX: An implicit finite element model for quantification of groundwater fluxes and hyporheic exchange in streams and rivers using radon

A quantitative understanding of groundwater-surface water interactions is vital for sustainable management of water quantity and quality. The noble gas radon-222 (Rn) is becoming increasingly used as a sensitive tracer to quantify groundwater discharge to wetlands, lakes, and rivers: a development driven by technical and methodological advances in Rn measurement. However, quantitative interpretation of these data is not trivial, and the methods used to date are based on the simplest solutions to the mass balance equation (e.g., first-order finite difference and inversion). Here we present a new implicit numerical model (FINIFLUX) based on finite elements for quantifying groundwater discharge to streams and rivers using Rn surveys at the reach scale (1–50 km). The model is coupled to a state-of-the-art parameter optimization code Parallel-PEST to iteratively solve the mass balance equation for groundwater discharge and hyporheic exchange. The major benefit of this model is that it is programed to be very simple to use, reduces nonuniqueness, and provides numerically stable estimates of groundwater fluxes and hyporheic residence times from field data. FINIFLUX was tested against an analytical solution and then implemented on two German rivers of differing magnitude, the Salzach (∼112 m3 s−1) and the Rote Main (∼4 m3 s−1). We show that using previous inversion techniques numerical instability can lead to physically impossible negative values, whereas the new model provides stable positive values for all scenarios. We hope that by making FINIFLUX freely available to the community that Rn might find wider application in quantifying groundwater discharge to streams and rivers and thus assist in a combined management of surface and groundwater systems.

Quantifying nitrate and oxygen reduction rates in the hyporheic zone using 222Rn to upscale biogeochemical turnover in rivers

Quantifying and upscaling chemical turnover in the hyporheic zone (HZ) is difficult due to limited reaction rate data, unknown carbon quality, and few methods for upscaling local measurements to river networks. Here we develop a method for quantifying reaction kinetics in situ in the HZ and upscaling biogeochemical turnover to catchment scales. Radon-222 was used to quantify water residence times in the HZ of the Roter Main River (RM), Germany. Residence times were then combined with O2, NO3−, CO2, DOC, and carbon quality (EEMs, SUVA) data to estimate Monod and first-order reaction rates. Monod parameters µmax and ksat for NO3− reduction were 11 µmol l−1 h−1 and 52 µmol l−1, respectively, while the first-order rate was 0.04 h−1. Carbon quality was highly bioavailable in the HZ and is unlikely to be limiting. Reaction kinetics was incorporated into the FINIFLUX model to upscale NO3− mass loss over a 32 km reach of the RM. The aims were to (1) to estimate hyporheic efficiency using Damkohler numbers (Da), and (2) calculate NO3− mass loss in the HZ over the reach. The Da analysis suggests that the hyporheic zone is inefficient for NO3− processing, however, this is somewhat misleading as the largest NO3− mass loss occurs at the shortest residence times where Da ≪ 1. This is due to the largest water flux occurring in the uppermost part of the sediment profile. Nitrate processing in the HZ accounted for 24 kg NO3− h−1 over the reach, which was 20% of the NO3− flux from the catchment.

Analytical modeling of hyporheic flow for in-stream bedforms: Perturbation method and implementation

Abstract Hyporheic flow and nutrient turnover in hyporheic systems are strongly influenced by in-stream bedforms. An accurate representation of topographical variations of the stream-streambed interface is therefore essential in analytical models in order to represent the couplings between hydrological and biogeochemical processes correctly. The classical Toth approach replaces the streambed surface topography by a flat surface which is identical to a truncation of the original physical flow domain into a rectangle. This simplification can lead to biased estimates of hyporheic flow and nutrient cycling within hyporheic systems. We present an alternative analytical modeling approach for solving hyporheic problems without domain truncation that explicitly accounts for topographical variations of the streambed. The presented approach is based on the application of perturbation theory. Applications of the method to hyporheic systems, ranging from the centimeter-scale of rippled bedforms to riffle structures of 10 m and larger scale, indicate a high accuracy of the approach.

Catchment Evapotranspiration and Runoff

The interplay between precipitation and evapotranspiration determines the input into the hydrological system of a catchment. Annual values of precipitation, evapotranspiration, and runoff measured at the catchment outlet for the 2002–2009 period were available. Annual precipitation clearly surmounted the sum of evapotranspiration and runoff. Part of the observed discrepancy might be due to the heterogeneity of precipitation and evapotranspiration within the catchment which has not been studied in sufficient detail. Annual evapotranspiration fluxes were remarkably constant during this period, whereas precipitation and runoff exhibited much larger interannual variability.