Flow and Transport Modeling in Heterogeneous Sediments Using an Integral Approach

An integral approach which can simultaneously model turbulent flow and transport at the sediment-water interface has been recently developed and validated for homogeneous sediment which was achieved by comparing numerical results to flume experiments on flow and transport over a rippled streambed and through the sediment for neutral, gaining, and losing conditions. In the present study, we validated the approach for heterogeneous conditions by comparing numerical simulations of flow and transport in heterogeneous sediment to analytical solutions as well as flume experiments on flow and transport through rippled streambed consisting of heterogeneous sediment. For this complex setup, simulation and experimental results agree well showing that flow and tracer transport prefer paths through areas with bigger grain diameters and higher porosities. The effect of flow redirections under losing and gaining conditions on hyporheic flow and residence times is discussed.

High-Resolution Integrated Transport Model for Studying Surface Water-Groundwater Interaction

Transport processes that lead to exchange of mass between surface water and groundwater play a significant role for the ecological functioning of aquatic systems, for hydrological processes and for biogeochemical transformations. In this study, we present a novel integral modeling approach for flow and transport at the sediment-water interface. The model allows us to simultaneously simulate turbulent surface and subsurface flow and transport with the same conceptual approach. For this purpose, a conservative transport equation was implemented to an existing approach that uses an extended version of the Navier-Stokes equations. Based on previous flume studies which investigated the spreading of a dye tracer under neutral, losing and gaining flow conditions the new solver is validated. Tracer distributions of the experiments are in close agreement with the simulations. The simulated flow paths are significantly affected by in- and outflowing groundwater flow. The highest velocities within the sediment are found for losing condition, which leads to shorter residence times compared to neutral and gaining conditions. The largest extent of the hyporheic exchange flow is observed under neutral condition. The new solver can be used for further examinations of cases that are not suitable for the conventional coupled models, for example, if Reynolds numbers are larger than 10. Moreover, results gained with the integral solver provide high-resolution information on pressure and velocity distributions at the rippled streambed, which can be used to improve flow predictions. This includes the extent of hyporheic exchange under varying ambient groundwater flow conditions.

Using computational fluid dynamics to describe H2S mass transfer across the water-air interface in sewers

For the past 70 years, researchers have dealt with the investigation of odour in sewer systems caused by hydrogen sulphide formations and the development of approaches to describe it. The state-of-the-art models are one-dimensional. At the same time, flow and transport phenomena in sewers can be three-dimensional, for example the air flow velocities in circular pipes or flow velocities of water and air in the reach of drop structures. Within the past years, increasing computational capabilities enabled the development of more complex models. This paper uses a three-dimensional two-phase computational fluid dynamics model to describe mass transfer phenomena between the two phases: water and air. The solver has been extended to be capable of accounting account for temperature dependency, the influence of pH value and a conversion to describe simulated air phase concentrations as partial pressure. Its capabilities are being explored in different application examples and its advantages compared to existing models are demonstrated in a highly complex three-dimensional test case. The resulting interH2SFoam solver is a significant step in the direction of describing and analysing H₂S emissions in sewers.

The fate of polar trace organic compounds in the hyporheic zone

The hyporheic zone (HZ) is often considered to efficiently remove polar trace organic compounds (TrOCs) from lotic systems, mitigating potential adverse effects of TrOCs on ecosystem functioning and drinking water production. Predicting the fate of TrOCs in the hyporheic zone (HZ) is difficult as the in-situ removal rate constants are not known and the biogeochemical factors as well as hydrological conditions controlling the removal efficiency are not fully understood. To determine the in-situ removal efficiency of the HZ for a variety of TrOCs as a function of the biogeochemical milieu, we conducted a field study in an urban river near Berlin, Germany. Subsurface flow was studied by time series of temperature depth profiles and the biogeochemical milieu of the HZ by concentration depth profiles. These results, in conjunction with a 1D advection-dispersion transport model, were used to calculate first-order removal rate constants of several polar TrOCs in the HZ. For the majority of TrOCs investigated, removal rate constants were strongly dependent on redox conditions, with significantly higher removal rates observed under predominantly suboxic (i.e. denitrifying) compared to anoxic (i.e. Fe and Mn reducing) conditions. Compared to previous studies on the fate of TrOCs in saturated sediments, half-lives within oxic/suboxic sections of the HZ were relatively low, attributable to the site-specific characteristics of the HZ in a stream dominated by wastewater treatment plant effluent. For nine out of thirteen investigated TrOCs, concentrations decreased significantly in the HZ with relative removal percentages ranging from 32% for primidone to 77% for gabapentin. For many TrOCs, removal efficiency decreased drastically as redox conditions became anoxic. For the majority of compounds investigated here, the HZ indeed acts as an efficient bioreactor that is capable of removing TrOCs along relatively short flow paths. Depending on the TrOC, removal capacity may be enhanced by either increasing the magnitude of groundwater-surface exchange fluxes, by increasing the total residence time in the HZ or the exposure time to suboxic zones, respectively.

Investigating groundwater-lake interactions by hydraulic heads and a water balance

Discharge of groundwater into lakes (lacustrine groundwater discharge, LGD) can play a major role in water balances of lakes. Unfortunately, studies often neglect this input path because of methodological difficulties in its determination. Direct measurements of LGD are labor-consuming and prone to error. The present study uses both spatially variable hydraulic-head data and meteorological data to estimate groundwater input by LGD and lake water output through infiltration. The study sites are two shallow, groundwater-fed lakes without any surface inflows or outflows. Horizontally interpolated groundwater heads were combined with lake water levels to obtain vertical hydraulic gradients between the aquifer and the lake, which are separated by a thick layer of lake bed sediment which has a much lower hydraulic conductivity than the underlying aquifer. By fitting the hydraulic gradient to the results of a simple mass balance and considering the process of clogging, we were able to estimate the hydraulic conductivity of the lake bed sediments. We calculated groundwater inputs by LGD and lake water outputs by infiltration on an annual basis. Although our method requires several assumptions, the results are reasonable and provide useful information about the exchange between the aquifer and the lake, which can, for example, be used for the calculation of nutrient mass balances.

Seeking a compromise between pharmaceutical pollution and phosphorus load: Management strategies for Lake Tegel, Berlin

Lake Tegel (Berlin, Germany) is controlled by two main inflows: inflow #1 (River Havel) is heavily phosphorus-laden, whereas inflow #2 is an artificial confluence that includes discharge from a municipal wastewater treatment plant distinguished by high levels of phosphorus and pharmaceuticals. To reduce the phosphorus load on the lake, a phosphorus elimination plant (PEP) is situated at inflow #2. Moreover, the two inflows are short-circuited by a pipeline that transfers part of the inflow #1 water to the PEP and finally releases it into inflow #2. The pipeline and the PEP have contributed to a continuous reduction in the total phosphorus concentration of Lake Tegel in the past 25 years. We investigate the question of whether the existing lake pipeline can also be used to reduce the amount of pharmaceuticals in Lake Tegel originating from inflow #2 by dilution with water from River Havel, by diverting part of inflow #2 around the lake, or by a combination of both strategies. The circulation pattern of Lake Tegel is complicated by complex bathymetry and numerous islands and is therefore highly sensitive to winds. We tested seven different management scenarios by hydrodynamic modeling for a period of 16 years with the two-dimensional version of the Princeton Ocean Model (POM). None of the scenarios provided a strategy optimal for both pharmaceuticals and phosphorus. Nonetheless, compound regimes, such as alternating the pipe flow direction or adding another pipeline, allowed the most abundant pharmaceutical (carbamazepine) to be reduced while maintaining the current phosphorus level. This study demonstrates the ability of immediate lake regulation measures to maintain water quality. In the case of Lake Tegel, the pipeline can be fully effective with regard to pharmaceuticals only in combination with additional efforts such as advanced pharmaceutical treatment of wastewater and/or phosphorus reduction in the River Havel catchment.

Exposure to human pharmaceuticals Carbamazepine, Ibuprofen and Bezafibrate causes molecular effects in Dreissena polymorpha

Carbamazepine (CBZ), Ibuprofen (IBU) and Bezafibrate (BEZ) were tested for their potential to bioaccumulate and provoke molecular changes in the non-target organism Dreissena polymorpha. mRNA changes of enzymes and other proteins involved in the prevention from protein damage (heat shock protein 70, hsp70) and oxidative stress (superoxide dismutase, SOD; catalase, CAT; metallothionein, MT), biotransformation (pi-class glutathione S-transferase, piGST; aryl hydrocarbon receptor, AH-R), elimination (P-glycoprotein, P-gp) and reversible protein posttranslational modification (protein phosphatase 2A, PP2A) served as molecular biomarkers. Mussels were exposed in a flow-through system to increasing concentrations of the three substances (1, 10, 100 and 1000 nM). The two lower concentrations correspond to environmentally relevant concentrations detected in surface and effluent waters, respectively. Measuring tissue concentration after one, four and seven days the uptake of CBZ and IBU by the mussels could be evidenced, whereas no accumulation data could be achieved for BEZ. The bioconcentration factor was highest for mussels exposed to the lowest CBZ and IBU concentrations, with 90 and 460-fold higher tissue concentration, respectively, after seven days. CBZ was the only substance tested which caused a significant increase in gill mRNA level of hsp70 after only one day exposure, evidencing the potential of CBZ to immediately provoke a stress condition and assumingly protein damage in gills. After longer exposure, mussels displayed down-regulated mRNA levels of hsp70 and SOD in gills, as well as of MT and P-gp in the digestive gland, hinting on an inhibitory character of CBZ. In IBU exposed mussels increased oxidant stress conditions were evidenced by induced mRNA levels in the digestive gland of CAT and MT, as well as SOD after one and four days, respectively. A concentration as found at sewage treatment plant effluents provoked an increase in transcript levels of piGST, suggesting enhanced need for biotransformation of IBU or by-products derived from oxidative stress. Also exposure to an environmentally relevant BEZ concentration provoked an immediate increase in piGST transcript level in the digestive gland followed by up-regulated hsp70 after four and seven days evidencing a chronic stress condition for the mussels.

The synthetic gestagen levonorgestrel impairs metamorphosis in Xenopus laevis by disruption of the thyroid system

Synthetic gestagens, including levonorgestrel (LNG), are active compounds in contraceptives, and several studies report their occurrence in surface waters. However, information about endocrine-disrupting effects in nontarget organisms is scarce. The present study investigated effects of LNG exposure on thyroid hormone-dependent metamorphosis of Xenopus laevis. Premetamorphic X. laevis tadpoles at Nieuwkoop and Faber (NF) stage 48 were exposed in a flow-through culture system to four LNG concentrations (10(-11), 10(-10), 10(-9), and 10(-8)M) over the period of metamorphosis. At NF 58 and 66, tadpoles were examined sex specifically. Developmental time and organismal responses were recorded and correlated with molecular and histopathological endpoints. Exposure to 10(-8)M LNG caused an inhibition of metamorphosis resulting in developmental arrest at early climax stages as giant tadpoles or tailed frogs. In brain-pituitary tissue of NF 58 tadpoles, gene expression of thyroid-stimulating hormone (β-subunit; TSHβ), TH receptor β (TRβ), and deiodinase type 3 (D3) was not changed. Instead, prolactin (PRL) messenger RNA (mRNA) was significantly increased by 10(-9)M LNG in females and by 10(-8)M LNG in both sexes. In NF 66 tadpoles, mRNA levels of TSHβ mRNA were significantly increased in the 10(-9) and 10(-8)M LNG treatment groups indicating a hypothyroid state. No changes of TRβ, D3, and PRL gene expression were detected. Histopathological evaluation of thyroid gland sections revealed no typical sign of hypothyroidism but rather an inactivated appearance of the thyroid. In conclusion, our data demonstrate for the first time a completely new aspect of thyroid system disruption caused by synthetic gestagens in developing amphibians.

Molecular effects and bioaccumulation of levonorgestrel in the non-target organism Dreissena polymorpha

Bioaccumulation and effects of the contraceptive hormone levonorgestrel were examined in the non-target organism Dreissena polymorpha. Molecular biomarkers of biotransformation, elimination, antioxidant defence and protein damage were analyzed after exposure to increasing concentrations of levonorgestrel in a flow-through system. The lowest concentration (0.312 μg L(-1)) was 100-fold bioconcentrated within four days. A decrease of the bioconcentration factor was observed within one week for the highest test concentrations (3.12 and 6.24 μg L(-1)) suggesting enhanced excretory processes. The immediate mRNA up-regulation of pi class glutathione S-transferase proved that phase II biotransformation processes were induced. Disturbance of fundamental cell functions was assumed since the aryl hydrocarbon receptor has been permanently down-regulated. mRNA up-regulation of P-glycoprotein, superoxide dismutase and metallothioneine suggested enhanced elimination processes and ongoing oxidative stress. mRNA up-regulation of heat shock protein 70 in mussels exposed to the two highest concentrations clearly indicated impacts on protein damage.

The beta-receptor blocker metoprolol alters detoxification processes in the non-target organism Dreissena polymorpha

Due to increasing amounts of pharmaceutically active compounds (PhACs) in the aquatic environment, their largely unknown effects to non-target organisms need to be assessed. This study examined physiological changes in the freshwater mussel Dreissena polymorpha exposed to increasing concentrations (0.534, 5.34, 53.4 and 534 microg L(-1)) of the beta-blocker metoprolol in a flow-through system for seven days. The two lower concentrations represent the environmentally relevant range. Surprisingly, metallothionein mRNA was immediately up-regulated in all treatments. For the two higher concentrations mRNA up-regulation in gills was found for P-glycoprotein after one day, and after four days for pi class glutathione S-transferase, demonstrating elimination and biotransformation processes, respectively. Additionally, catalase and superoxide dismutase were up-regulated in the digestive gland indicating oxidative stress. In all treated mussels a significant up-regulation of heat shock protein mRNA was observed in gills after four days, which suggests protein damage and the requirement for repair processes. Metoprolol was 20-fold bioaccumulated for environmentally relevant concentrations.

Model based verification and prognosis of acidification and sulphate releasing processes downstream of a former sewage field in Berlin (Germany)

An ammonium contamination plume originating from sewage field management practices over several decades is affecting the water quality at the well fields of the Friedrichshagen waterworks in Berlin, Germany. Because hydraulic measures were unsuccessful due to the fixation of ammonium on the aquifer matrix by cation exchange, an in situ nitrification measure by injection of oxygen gas was chosen to protect the extraction wells. In order to assess the hydro chemical processes accompanying this in situ measure, reactive transport modelling was performed. The relevant processes are the dissolution of oxygen gas and the nitrification of ammonium which initiate secondary geochemical processes like sulphate release, acidification and hardening. The reactive transport modelling began with the deduction of a reaction network, followed by the mathematical formulation and incorporation of reactive terms into a reactive transport solver. Two model versions were set up: (1) a simplified large scale model to evaluate the long-term reaction zoning to be expected due to permanent oxygen gas injection, and (2) a verification of the monitored hydrochemistry during a first field test performed near the contamination source. The results of reactive transport modelling demonstrate that in situ injection of oxygen gas will be effective in reducing the ammonium load from the well fields, and that acidification processes near the production wells can be minimized. Finally, a line of gas injection wells extending over the whole width of the ammonium contamination plume will be constructed to protect the well fields from further ammonium load.

Transient leakance and infiltration characteristics during lake bank filtration

Infiltration capacity of bank filtration systems depends on water extraction and hydraulic resistance of the bed sediments. Lakebed hydraulics may be especially affected by clogging, which is dependent on settlement of fine particles, redox potential, and other factors. In the field, most of these processes are difficult to quantify, and thus, when calculating response to pumping the water flux across the sediment surface is assumed to be linearly dependent on the hydraulic gradient. However, this assumption was not adequate to describe conditions at a bank filtration site located at Lake Tegel, Berlin, Germany. Hence, we first assumed the leakage coefficient (or leakance) is spatially distributed and also temporally variant. Furthermore, observations show that the leakance is considerably higher in shallow than in deeper areas; hence, leakance was assumed to be dependent on the existence and thickness of an unsaturated zone below the lake. The proposed explanation of spatial and temporal variability in leakance involves a hypothesis for redox dependent and reversible biogeochemical clogging, supported by geochemical observations in surface water and ground water. Four leakance approaches are implemented in the ground water flow code MODFLOW2000 and calibrated by inverse modeling using the parameter estimation software PEST. These concepts are evaluated by examining the fit to the hydraulic heads, to infiltration measurements, transport modeling results, and considering the degrees of freedom due to the number of calibration parameters. The leakage concept based on the assumption of the influence of an unsaturated zone on clogging processes best explains the field data.

Ecological scenarios analyzed and evaluated by a shallow lake model

We applied the complex ecosystem model EMMO, which was adopted to the shallow lake Müggelsee (Germany), in order to evaluate a large set of ecological scenarios. By means of EMMO, 33 scenarios and 17 indicators were defined to characterize their effects on the lake ecosystem. The indicators were based on model outputs of EMMO and can be separated into biological indicators, such as chlorophyll-a and cyanobacteria, and hydro-chemical indicators, such as phosphorus. The question to be solved was, what is the ranking of the scenarios based on their characterization by these 17 indicators? And how can we handle high quantities of complex data within evaluation procedures? The scenario evaluation was performed by partial order theory which, however, did not provide a clear result. By subsequently applying the hierarchical cluster analysis (complete linkage) it was possible to reduce the data matrix to indicator and scenario representatives. Even though this step implies losses of information, it simplifies the application of partial order theory and the post processing by METEOR. METEOR is derived from partial order theory and allows the stepwise aggregation of indicators, which subsequently leads to a distinct and clear decision. In the final evaluation result the best scenario was the one which defines a minimum nutrient input and no phosphorus release from the sediment while the worst scenario is characterized by a maximum nutrient input and extensive phosphorus release from the sediment. The reasonable and comprehensive results show that the combination of partial order, cluster analysis and METEOR can handle big amounts of data in a very clear and transparent way, and therefore is ideal in the context of complex ecosystem models, like that we applied.

Modeling seasonal redox dynamics and the corresponding fate of the pharmaceutical residue phenazone during artificial recharge of groundwater

Reactive multicomponent transport modeling was used to investigate and quantify the factors that affect redox zonation and the fate of the pharmaceutical residue phenazone during artificial recharge of groundwater at an infiltration site in Berlin, Germany. The calibrated model and the corresponding sensitivity analysis demonstrated thattemporal and spatial redox zonation at the study site was driven by seasonally changing, temperature-dependent organic matter degradation rates. Breakthrough of phenazone at monitoring wells occurred primarily during the warmer summer months, when anaerobic conditions developed. Assuming a redox-sensitive phenazone degradation behavior the model results provided an excellent agreement between simulated and measured phenazone concentrations. Therefore, the fate of phenazone was shown to be indirectly controlled by the infiltration water temperature through its effect on the aquifer's redox conditions. Other factors such as variable residence times appeared to be of less importance.