Contaminant degradation in physically and chemically heterogeneous aquifers

The value of simplified models of radionuclide transport for the safety assessment of nuclear waste repositories: A benchmark study

In order to assess sites for a deep geological repository for storing high-level nuclear waste safely in Germany, various numerical models and tools will be in use. For their interaction within one workflow, their reproducibility, and reliability version-controlled open-source solutions and careful documentation of model setups, results and verifications are of special value. However, spatially fully resolved models including all relevant physical and chemical processes are neither computationally feasible for large domains nor is the data typically available to parameterize such models. Thus, simplified models are crucial for the pre-assessment of possible sites to narrow down the list of suitable candidates for which detailed site investigations and fully resolved models will be done at a later stage. Still, the accuracy of these simplified models is of importance as the pre-assessment of suitable sites will be based on them. In this study, we compare the modelling capabilities of TransPyREnd, a one-dimensional transport code based on finite differences, specifically developed for the fast estimation of radionuclide transport by the German federal company for radioactive waste disposal (BGE), with OpenGeoSys, which is a modelling platform based on finite elements in up to three spatial dimensions. Both codes are used in the site selection procedure for the German nuclear waste repository. The comparison of the model results of TransPyREnd and OpenGeoSys is augmented by comparisons with an analytical solution for a homogeneous material. For the purpose of numerical benchmarking, we consider a geological profile located in southern Germany as an example where the hypothetical repository is located in a clay-stone formation. TransPyREnd and OpenGeoSys yield overall similar results. However, both codes use different discretizations which impact is the highest for strongly sorbing compounds, while the difference gets negligible for less sorbing and more diffusive compounds as higher diffusion tends to blur the initial conditions. Overall, the OpenGeoSys model is more exact whereas the TransPyREnd model has considerable faster run times. We found in our example, that significant substance amounts are only leaving the host rock formation, if apparent diffusion is high, for which case both codes give similar results, while relative differences are considerable for strongly sorbing compounds. However, in the latter case no significant substance amount of radionuclides leaves the host-rock formation, thus deeming the differences in the model results minor for the overall safety assessments of sites.

Virtual experiments to assess opportunities and pitfalls of CSIA in physical-chemical heterogeneous aquifers

Degradation of chlorinated ethenes (CEs) in low conductivity layers of aquifers reduces pollution plume tailing and accelerates natural attenuation timeframes. The degradation pathways involved are often different from those in the higher conductive layers and might go undetected when only highly conductive layers are targeted in site assessments. Reactive transport model simulations (PHT3D in FloPy) were executed to assess the performance of dual carbon and chlorine compound specific stable isotope analysis (CSIA) in degradation pathway identification and quantification in a coupled physical-chemical heterogeneous virtual aquifer. Degradation rate constants were assumed correlated to the hydraulic conductivity: positively for oxidative transformation (higher oxygen availability in coarser sands) and negatively for chemical reduction (higher content of reducing solids in finer sediments). Predicted carbon isotope ratios were highly heterogeneous. They generally increased downgradient of the pollution source but the large variation across depth illustrates that monotonously increasing isotope ratios downgradient, as were associated with the oxidative component, are not necessarily a common situation when degradation is favorable in low conductivity layers. Dual carbon-chlorine CSIA performed well in assessing the occurrence of the spatially separated degradation pathways and the overall degradation, provided appropriate enrichment factors were known and sufficiently different. However, pumping to obtain groundwater samples especially from longer well screens causes a bias towards overestimation of the contribution of oxidative transformation associated with the higher conductive zones. As degradation was less intense in these highly conductive zones under the simulated conditions, overall degradation was underestimated. In contrast, in the usual case of limited CSIA data, dual CSIA plots may rather indicate dominance of chemical reduction, while oxidative transformation could go unnoticed, despite being an equally important degradation pathway.

Effects of local transverse dispersion on macro-scale coefficients of oxygen-limited biodegradation in a stratified formation

The correct characterization of macro-scale contaminant transport and transformation rates is an important issue for modeling reactive transport in heterogeneous aquifers. While previous studies have investigated field-scale heterogeneity of transport and biochemical properties, the effects of local transverse dispersion on macro-scale transport and transformation rates have not been well understood. In this paper, the process of oxygen-limited biodegradation in a stratified aquifer is analysed by spectral perturbation approach, and longitudinal macrodispersivity, effective biodegradation rate, effective retardation factor and effective velocity are derived for the coupled transport equations of a system consisting of a contaminant and an oxidizing agent (oxygen). The effects of local transverse dispersion on these macro-scale coefficients are studied. It is shown that local transverse dispersion can smooth the heterogeneity in biodegradation and sorption processes and enlarge effective biodegradation rate and retardation factor. The local transverse dispersion can also limit the effects of heterogeneity in biodegradation process on longitudinal macrodispersivities and effective velocities for the contaminant and dissolved oxygen. But the effects of heterogeneity in sorption process on the contaminant macrodispersivity is likely to be magnified by local transverse dispersion.

Exposure-time based modeling of nonlinear reactive transport in porous media subject to physical and geochemical heterogeneity

Transport of reactive solutes in groundwater is affected by physical and chemical heterogeneity of the porous medium, leading to complex spatio-temporal patterns of concentrations and reaction rates. For certain cases of bioreactive transport, it could be shown that the concentrations of reactive constituents in multi-dimensional domains are approximately aligned with isochrones, that is, lines of identical travel time, provided that the chemical properties of the matrix are uniform. We extend this concept to combined physical and chemical heterogeneity by additionally considering the time that a water parcel has been exposed to reactive materials, the so-called exposure time. We simulate bioreactive transport in a one-dimensional domain as function of time and exposure time, rather than space. Subsequently, we map the concentrations to multi-dimensional heterogeneous domains by means of the mean exposure time at each location in the multi-dimensional domain. Differences in travel and exposure time at a given location are accounted for as time difference. This approximation simplifies reactive-transport simulations significantly under conditions of steady-state flow when reactions are restricted to specific locations. It is not expected to be exact in realistic applications because the underlying assumption, such as neglecting transverse mixing altogether, may not hold. We quantify the error introduced by the approximation for the hypothetical case of a two-dimensional, binary aquifer made of highly-permeable, non-reactive and low-permeable, reactive materials releasing dissolved organic matter acting as electron donor for aerobic respiration and denitrification. The kinetically controlled reactions are catalyzed by two non-competitive bacteria populations, enabling microbial growth. Even though the initial biomass concentrations were uniform, the interplay between transport, non-uniform electron-donor supply, and bio-reactions led to distinct spatial patterns of the two types of biomass at late times. Results obtained by mapping the exposure-time based results to the two-dimensional domain are compared with simulations based on the two-dimensional, spatially explicit advection-dispersion-reaction equation. Once quasi-steady state has been reached, we find a good agreement in terms of the chemical-compound concentrations between the two approaches inside the reactive zones, whereas the exposure-time based model is not able to capture reactions occurring in the zones with zero electron-donor release. We conclude that exposure-time models provide good approximations of nonlinear bio-reactive transport when transverse mixing is not the overall controlling process and all reactions are essentially restricted to distinct reactive zones.

The effects of physical and geochemical heterogeneities on hydro-geochemical transport and effective reaction rates

The role of coupled physical and geochemical heterogeneities in hydro-geochemical transport is investigated by simulating three-dimensional transport in a heterogeneous system with kinetic mineral reactions. Ensembles of 100 physically heterogeneous realizations were simulated for three geochemical conditions: 1) spatially homogeneous reactive mineral surface area, 2) reactive surface area positively correlated to hydraulic heterogeneity, and 3) reactive surface area negatively correlated to hydraulic heterogeneity. Groundwater chemistry and the corresponding effective reaction rates were calculated at three transverse planes to quantify differences in plume evolution due to heterogeneity in mineral reaction rates and solute residence time (τ). The model is based on a hypothetical CO2 intrusion into groundwater from a carbon capture utilization and storage (CCUS) operation where CO2 dissolution and formation of carbonic acid created geochemical dis-equilibrium between fluids and the mineral galena that resulted in increased aqueous lead (Pb(2+)) concentrations. Calcite dissolution buffered the pH change and created conditions of galena oversaturation, which then reduced lead concentrations along the flow path. Near the leak kinetic geochemical reactions control the release of solutes into the fluid, but further along the flow path mineral solubility controls solute concentrations. Simulation results demonstrate the impact of heterogeneous distribution of geochemical reactive surface area in coordination with physical heterogeneity on the effective reaction rate (Krxn,eff) and Pb(2+) concentrations within the plume. Dissimilarities between ensemble Pb(2+) concentration and Krxn,eff are attributed to how geochemical heterogeneity affects the time (τeq) and therefore advection distance (Leq) required for the system to re-establish geochemical equilibrium. Only after geochemical equilibrium is re-established, Krxn,eff and Pb(2+) concentrations are the same for all three geochemical conditions. Correlation between reactive surface area and hydraulic conductivity, either positive or negative, results in variation in τeq and Leq.

General quantification of catchment-scale nutrient and pollutant transport through the subsurface to surface and coastal waters

This study develops a general quantification framework for consistent intermodel and intercatchment comparison of the nutrient and pollutant mass loading from multiple sources in a catchment area to downstream surface and coastal waters. The framework accounts for the wide spectrum of different transport pathways and travel times through the subsurface (soil, groundwater, sediment) and the linked surface (streams, lakes, wetlands) water systems of a catchment. The account is based on key flow partitioning and mass delivery fractions, which can be quantified differently by different flow and transport and reaction models. The framework application is exemplified for two Swedish catchment cases with regard to the transport of phosphorus and of a generic attenuating solute. The results show essential differences in model quantifications of transport pathways and temporal spreading, with important implications for our understanding of cause and effect in the catchment-scale nutrient and pollutant loading to downstream waters.