@article{SchreiberMunzSalzmannetal.2021, author = {Schreiber, Lisa and Munz, Matthias and Salzmann, Thomas and Oswald, Sascha E.}, title = {Coupled simulation of groundwater and drainage dynamics in a coastal fen}, series = {Grundwasser : Zeitschrift der Fachsektion Hydrogeologie in der Deutschen Gesellschaft f{\"u}r Geowissenschaften (FH-DGG)}, volume = {26}, journal = {Grundwasser : Zeitschrift der Fachsektion Hydrogeologie in der Deutschen Gesellschaft f{\"u}r Geowissenschaften (FH-DGG)}, number = {3}, publisher = {Springer}, address = {Heidelberg}, issn = {1430-483X}, doi = {10.1007/s00767-021-00486-y}, pages = {289 -- 304}, year = {2021}, abstract = {Coastal wetlands are characterized by continued human influence, e.g. with drainage ditches, coastal dikes or landscape restoration. In addition, it is important to understand the complex interactions with the sea to predict impacts of further development. In the present study the aim was to analyze surface and subsurface flow in a coastal wetland located at the Baltic Sea coastline near Warnemunde (Germany) to quantify water exchange with the Baltic Sea and analyze the effect of a storm flood event on saline intrusion. A 3-D transient groundwater model and a one-dimensional surface water model were set up and calibrated by using hydraulic head measurements. The results indicate that in addition to ditch flow, groundwater discharge to the Baltic Sea often has a significant influence on the overall water budget of the fen. From the transient modelling it became evident that water exchange between groundwater in the fen and the Baltic Sea depends on sea level and very often fluctuates between seaward and landward flow directions on daily to weekly time scales.}, language = {de} } @article{WangOswaldGraeffetal.2019, author = {Wang, Wei-shi and Oswald, Sascha Eric and Gr{\"a}ff, Thomas and Lensing, Hermann Josef and Liu, Tie and Strasser, Daniel and Munz, Matthias}, title = {Impact of river reconstruction on groundwater flow during bank filtration assessed by transient three-dimensional modelling of flow and heat transport}, series = {Hydrogeology journal : official journal of the International Association of Hydrogeologists}, volume = {28}, journal = {Hydrogeology journal : official journal of the International Association of Hydrogeologists}, number = {2}, publisher = {Springer}, address = {Berlin ; Heidelberg [u.a.]}, issn = {1431-2174}, doi = {10.1007/s10040-019-02063-3}, pages = {723 -- 743}, year = {2019}, abstract = {Bank filtration (BF) is an established indirect water-treatment technology. The quality of water gained via BF depends on the subsurface capture zone, the mixing ratio (river water versus ambient groundwater), spatial and temporal distribution of subsurface travel times, and subsurface temperature patterns. Surface-water infiltration into the adjacent aquifer is determined by the local hydraulic gradient and riverbed permeability, which could be altered by natural clogging, scouring and artificial decolmation processes. The seasonal behaviour of a BF system in Germany, and its development during and about 6 months after decolmation (canal reconstruction), was observed with a long-term monitoring programme. To quantify the spatial and temporal variation in the BF system, a transient flow and heat transport model was implemented and two model scenarios, 'with' and 'without' canal reconstruction, were generated. Overall, the simulated water heads and temperatures matched those observed. Increased hydraulic connection between the canal and aquifer caused by the canal reconstruction led to an increase of similar to 23\% in the already high share of BF water abstracted by the nearby waterworks. Subsurface travel-time distribution substantially shifted towards shorter travel times. Flow paths with travel times <200 days increased by similar to 10\% and those with <300 days by 15\%. Generally, the periodic temperature signal, and the summer and winter temperature extrema, increased and penetrated deeper into the aquifer. The joint hydrological and thermal effects caused by the canal reconstruction might increase the potential of biodegradable compounds to further penetrate into the aquifer, also by potentially affecting the redox zonation in the aquifer.}, language = {en} } @article{BarkowOswaldLensingetal.2020, author = {Barkow, Isolde S. and Oswald, Sascha Eric and Lensing, Hermann Josef and Munz, Matthias}, title = {Seasonal dynamics modifies fate of oxygen, nitrate, and organic micropollutants during bank filtration}, series = {Environmental science and pollution research : official organ of the EuCheMS Division for Chemistry and the Environment, EuCheMS DCE}, volume = {28}, journal = {Environmental science and pollution research : official organ of the EuCheMS Division for Chemistry and the Environment, EuCheMS DCE}, number = {8}, publisher = {Springer}, address = {Heidelberg}, issn = {0944-1344}, doi = {10.1007/s11356-020-11002-9}, pages = {9682 -- 9700}, year = {2020}, abstract = {Bank filtration is considered to improve water quality through microbially mediated degradation of pollutants and is suitable for waterworks to increase their production. In particular, aquifer temperatures and oxygen supply have a great impact on many microbial processes. To investigate the temporal and spatial behavior of selected organic micropollutants during bank filtration in dependence of relevant biogeochemical conditions, we have set up a 2D reactive transport model using MODFLOW and PHT3D under the user interface ORTI3D. The considered 160-m-long transect ranges from the surface water to a groundwater extraction well of the adjacent waterworks. For this purpose, water levels, temperatures, and chemical parameters were regularly measured in the surface water and groundwater observation wells over one and a half years. To simulate the effect of seasonal temperature variations on microbial mediated degradation, we applied an empirical temperature factor, which yields a strong reduction of the degradation rate at groundwater temperatures below 11 degrees C. Except for acesulfame, the considered organic micropollutants are substantially degraded along their subsurface flow paths with maximum degradation rates in the range of 10(-6) mol L-1 s(-1). Preferential biodegradation of phenazone, diclofenac, and valsartan was found under oxic conditions, whereas carbamazepine and sulfamethoxazole were degraded under anoxic conditions. This study highlights the influence of seasonal variations in oxygen supply and temperature on the fate of organic micropollutants in surface water infiltrating into an aquifer.}, language = {en} } @article{MunzOswaldSchaefferlingetal.2019, author = {Munz, Matthias and Oswald, Sascha Eric and Schaefferling, Robin and Lensing, Hermann Josef}, title = {Temperature-dependent redox zonation, nitrate removal and attenuation of organic micropollutants during bank filtration}, series = {Water research}, volume = {162}, journal = {Water research}, publisher = {Elsevier}, address = {Oxford}, issn = {0043-1354}, doi = {10.1016/j.watres.2019.06.041}, pages = {225 -- 235}, year = {2019}, abstract = {River bank filtration (RBF) is considered to efficiently remove nitrate and trace organic micropollutants (OMP) from polluted surface waters. This is essential for maintaining good groundwater quality and providing high quality drinking water. Predicting the fate of OMP during RBF is difficult as the biogeochemical factors controlling the removal efficiency are not fully understood. To determine in-situ removal efficiency and degradation rates of nitrate and OMP indicator substances we conducted a field study in a RBF system during a period of one and a half years incorporating temporally and spatially varying redox conditions and temperature changes typically occurring in temperate climates. RBF was analyzed by means of mixing ratios between infiltrated river water and groundwater as well as average residence times of surface water towards the individual groundwater observation wells. These results were used to calculate temperature dependent first order degradation rates of redox sensitive species and several OMP. Five out of ten investigated OMP were completely removed along RBF pathways. We demonstrate that degradation rates of several OMP during bank filtration were controlled by redox conditions and temperature whereby temperature itself also had a significant influence on the extent of the most reactive oxic zone. The seasonal variations in temperature alone could explain a considerable percentage of the variance in dissolved oxygen (34\%), nitrate (81\%) as well as the OMPs diclofenac (44\%) and sulfamethoxazole (76\%). Estimated in-situ degradation rates roughly varied within one order of magnitude for temperature changes between 5 degrees C and 20 degrees C. This study highlights that temporal variability in temperature and redox zonation is a significant factor for migration and degradation of nitrate and several OMPs. (C) 2019 Elsevier Ltd. All rights reserved.}, language = {en} } @article{MunzSchmidt2017, author = {Munz, Matthias and Schmidt, Christian}, title = {Estimation of vertical water fluxes from temperature time series by the inverse numerical computer program FLUX-BOT}, series = {Hydrological processes}, volume = {31}, journal = {Hydrological processes}, publisher = {Wiley}, address = {Hoboken}, issn = {0885-6087}, doi = {10.1002/hyp.11198}, pages = {2713 -- 2724}, year = {2017}, abstract = {The application of heat as a hydrological tracer has become a standard method for quantifying water fluxes between groundwater and surface water. The typical application is to estimate vertical water fluxes in the shallow subsurface beneath streams or lakes. For this purpose, time series of temperatures in the surface water and in the sediment are measured and evaluated by a vertical 1D representation of heat transport by advection and conduction. Several analytical solutions exist to calculate the vertical water flux from the measured temperatures. Although analytical solutions can be easily implemented, they are restricted to specific boundary conditions such as a sinusoidal upper temperature boundary. Numerical solutions offer higher flexibility in the selection of the boundary conditions. This, in turn, reduces the effort of data preprocessing, such as the extraction of the diurnal temperature variation from the raw data. Here, we present software to estimate water fluxes based on temperaturesFLUX-BOT. FLUX-BOT is a numerical code written in MATLAB that calculates vertical water fluxes in saturated sediments based on the inversion of measured temperature time series observed at multiple depths. FLUX-BOT applies a centred Crank-Nicolson implicit finite difference scheme to solve the one-dimensional heat advection-conduction equation. FLUX-BOT includes functions for the inverse numerical routines, functions for visualizing the results, and a function for performing uncertainty analysis. We present applications of FLUX-BOT to synthetic and to real temperature data to demonstrate its performance.}, language = {en} } @article{MunzOswaldSchmidt2017, author = {Munz, Matthias and Oswald, Sascha Eric and Schmidt, Christian}, title = {Coupled Long-Term Simulation of Reach-Scale Water and Heat Fluxes Across the River-Groundwater Interface for Retrieving Hyporheic Residence Times and Temperature Dynamics}, series = {Water resources research}, volume = {53}, journal = {Water resources research}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0043-1397}, doi = {10.1002/2017WR020667}, pages = {8900 -- 8924}, year = {2017}, abstract = {Flow patterns in conjunction with seasonal and diurnal temperature variations control ecological and biogeochemical conditions in hyporheic sediments. In particular, hyporheic temperatures have a great impact on many temperature-sensitive microbial processes. In this study, we used 3-D coupled water flow and heat transport simulations applying the HydroGeoSphere code in combination with high-resolution observations of hydraulic heads and temperatures to quantify reach-scale water and heat flux across the river-groundwater interface and hyporheic temperature dynamics of a lowland gravel bed river. The model was calibrated in order to constrain estimates of the most sensitive model parameters. The magnitude and variations of the simulated temperatures matched the observed ones, with an average mean absolute error of 0.7 degrees C and an average Nash Sutcliffe efficiency of 0.87. Our results indicate that nonsubmerged streambed structures such as gravel bars cause substantial thermal heterogeneity within the saturated sediment at the reach scale. Individual hyporheic flow path temperatures strongly depend on the flow path residence time, flow path depth, river, and groundwater temperature. Variations in individual hyporheic flow path temperatures were up to 7.9 degrees C, significantly higher than the daily average (2.8 degrees C), but still lower than the average seasonal hyporheic temperature difference (19.2 degrees C). The distribution between flow path temperatures and residence times follows a power law relationship with exponent of about 0.37. Based on this empirical relation, we further estimated the influence of hyporheic flow path residence time and temperature on oxygen consumption which was found to partly increase by up to 29\% in simulations.}, language = {en} } @article{MunzOswaldSchmidt2016, author = {Munz, Matthias and Oswald, Sascha Eric and Schmidt, Christian}, title = {Analysis of riverbed temperatures to determine the geometry of subsurface water flow around in-stream geomorphological structures}, series = {Journal of hydrology}, volume = {539}, journal = {Journal of hydrology}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0022-1694}, doi = {10.1016/j.jhydrol.2016.05.012}, pages = {74 -- 87}, year = {2016}, abstract = {The analytical evaluation of diurnal temperature variation in riverbed sediments provides detailed information on exchange fluxes between rivers and groundwater. The underlying assumption of the stationary, one-dimensional vertical flow field is frequently violated in natural systems where subsurface water flow often has a significant horizontal component. In this paper, we present a new methodology for identifying the geometry of the subsurface flow field using vertical temperature profiles. The statistical analyses are based on model optimisation and selection and are used to evaluate the shape of vertical amplitude ratio profiles. The method was applied to multiple profiles measured around in-stream geomorphological structures in a losing reach of a gravel bed river. The predominant subsurface flow field was systematically categorised in purely vertical and horizontal (hyporheic, parafluvial) components. The results highlight that river groundwater exchange flux at the head, crest and tail of geomorphological structures significantly deviated from the one-dimensional vertical flow, due to a significant horizontal component. The geometry of the subsurface water flow depended on the position around the geomorphological structures and on the river level. The methodology presented in this paper features great potential for characterising the spatial patterns and temporal dynamics of complex subsurface flow geometries by using measured temperature time series in vertical profiles. (C) 2016 Elsevier B.V. All rights reserved.}, language = {en} } @phdthesis{Munz2017, author = {Munz, Matthias}, title = {Water flow and heat transport modelling at the interface between river and aquifer}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-404319}, school = {Universit{\"a}t Potsdam}, pages = {XIII, 123}, year = {2017}, abstract = {The functioning of the surface water-groundwater interface as buffer, filter and reactive zone is important for water quality, ecological health and resilience of streams and riparian ecosystems. Solute and heat exchange across this interface is driven by the advection of water. Characterizing the flow conditions in the streambed is challenging as flow patterns are often complex and multidimensional, driven by surface hydraulic gradients and groundwater discharge. This thesis presents the results of an integrated approach of studies, ranging from the acquisition of field data, the development of analytical and numerical approaches to analyse vertical temperature profiles to the detailed, fully-integrated 3D numerical modelling of water and heat flux at the reach scale. All techniques were applied in order to characterize exchange flux between stream and groundwater, hyporheic flow paths and temperature patterns. The study was conducted at a reach-scale section of the lowland Selke River, characterized by distinctive pool riffle sequences and fluvial islands and gravel bars. Continuous time series of hydraulic heads and temperatures were measured at different depths in the river bank, the hyporheic zone and within the river. The analyses of the measured diurnal temperature variation in riverbed sediments provided detailed information about the exchange flux between river and groundwater. Beyond the one-dimensional vertical water flow in the riverbed sediment, hyporheic and parafluvial flow patterns were identified. Subsurface flow direction and magnitude around fluvial islands and gravel bars at the study site strongly depended on the position around the geomorphological structures and on the river stage. Horizontal water flux in the streambed substantially impacted temperature patterns in the streambed. At locations with substantial horizontal fluxes the penetration depths of daily temperature fluctuations was reduced in comparison to purely vertical exchange conditions. The calibrated and validated 3D fully-integrated model of reach-scale water and heat fluxes across the river-groundwater interface was able to accurately represent the real system. The magnitude and variations of the simulated temperatures matched the observed ones, with an average mean absolute error of 0.7 °C and an average Nash Sutcliffe Efficiency of 0.87. The simulation results showed that the water and heat exchange at the surface water-groundwater interface is highly variable in space and time with zones of daily temperature oscillations penetrating deep into the sediment and spots of daily constant temperature following the average groundwater temperature. The average hyporheic flow path temperature was found to strongly correlate with the flow path residence time (flow path length) and the temperature gradient between river and groundwater. Despite the complexity of these processes, the simulation results allowed the derivation of a general empirical relationship between the hyporheic residence times and temperature patterns. The presented results improve our understanding of the complex spatial and temporal dynamics of water flux and thermal processes within the shallow streambed. Understanding these links provides a general basis from which to assess hyporheic temperature conditions in river reaches.}, language = {en} } @article{MunzOswaldSchmidt2011, author = {Munz, Matthias and Oswald, Sascha Eric and Schmidt, C.}, title = {Sand box experiments to evaluate the influence of subsurface temperature probe design on temperature based water flux calculation}, series = {Hydrology and earth system sciences : HESS}, volume = {15}, journal = {Hydrology and earth system sciences : HESS}, number = {11}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1027-5606}, doi = {10.5194/hess-15-3495-2011}, pages = {3495 -- 3510}, year = {2011}, abstract = {Quantification of subsurface water fluxes based on the one dimensional solution to the heat transport equation depends on the accuracy of measured subsurface temperatures. The influence of temperature probe setup on the accuracy of vertical water flux calculation was systematically evaluated in this experimental study. Four temperature probe setups were installed into a sand box experiment to measure temporal highly resolved vertical temperature profiles under controlled water fluxes in the range of +/- 1.3 md(-1). Pass band filtering provided amplitude differences and phase shifts of the diurnal temperature signal varying with depth depending on water flux. Amplitude ratios of setups directly installed into the saturated sediment significantly varied with sand box hydraulic gradients. Amplitude ratios provided an accurate basis for the analytical calculation of water flow velocities, which matched measured flow velocities. Calculated flow velocities were sensitive to thermal properties of saturated sediment and to temperature sensor spacing, but insensitive to thermal dispersivity equal to solute dispersivity. Amplitude ratios of temperature probe setups indirectly installed into piezometer pipes were influenced by thermal exchange processes within the pipes and significantly varied with water flux direction only. Temperature time lags of small sensor distances of all setups were found to be insensitive to vertical water flux.}, language = {en} } @article{MunzKrauseTecklenburgetal.2011, author = {Munz, Matthias and Krause, Stefan and Tecklenburg, Christina and Binley, Andrew}, title = {Reducing monitoring gaps at the aquifer-river interface by modelling groundwater-surface water exchange flow patterns}, series = {Hydrological processes}, volume = {25}, journal = {Hydrological processes}, number = {23}, publisher = {Wiley-Blackwell}, address = {Malden}, issn = {0885-6087}, doi = {10.1002/hyp.8080}, pages = {3547 -- 3562}, year = {2011}, abstract = {This study investigates spatial patterns and temporal dynamics of aquifer-river exchange flow at a reach of the River Leith, UK. Observations of sub-channel vertical hydraulic gradients at the field site indicate the dominance of groundwater up-welling into the river and the absence of groundwater recharge from surface water. However, observed hydraulic heads do not provide information on potential surface water infiltration into the top 0-15 cm of the streambed as these depths are not covered by the existing experimental infrastructure. In order to evaluate whether surface water infiltration is likely to occur outside the 'window of detection', i.e. the shallow streambed, a numerical groundwater model is used to simulate hydrological exchanges between the aquifer and the river. Transient simulations of the successfully validated model (Nash and Sutcliff efficiency of 0.91) suggest that surface water infiltration is marginal and that the possibility of significant volumes of surface water infiltrating into non-monitored shallow streambed sediments can be excluded for the simulation period. Furthermore, the simulation results show that with increasing head differences between river and aquifer towards the end of the simulation period, the impact of streambed topography and hydraulic conductivity on spatial patterns of exchange flow rates decreases. A set of peak flow scenarios with altered groundwater-surface water head gradients is simulated in order to quantify the potential for surface water infiltration during characteristic winter flow conditions following the observation period. The results indicate that, particularly at the beginning of peak flow conditions, head gradients are likely to cause substantial increase in surface water infiltration into the streambed. The study highlights the potential for the improvement of process understanding of hyporheic exchange flow patterns at the stream reach scale by simulating aquifer-river exchange fluxes with a standard numerical groundwater model and a simple but robust model structure and parameterization. Copyright}, language = {en} }