"Feels Like an Indie Game"-Evaluation of a Virtual Field Trip Prototype on Radioactive Waste Management Research for University Education.

This article describes the design and evaluation of a virtual field trip on the topic of radioactive waste management research for university education. We created an interactive virtual tour through the Mont Terri underground research laboratory by enhancing the virtual experiment information system, designed for domain experts, with background information, illustrations, tasks, tests, and an improved user interface. To put the tour's content into context, a conventional introductory presentation on the final disposal of radioactive waste was added. A user study with 22 participants proved a good perceived usability of the virtual tour and the virtual field trip's ability to transfer knowledge. These results suggest a benefit of employing virtual field trips in geoscientific university courses. In addition, it is conceivable to use the virtual field trip as a tool for science communication in the context of participatory processes during nuclear waste disposal site selection processes.

Modeling the impacts of plants and internal organic carbon on remediation performance in the integrated vertical flow constructed wetland.

The integrated vertical flow (IVF) constructed wetland consists of two or more chambers with heterogeneous flow patterns and strong aeration capability, possesses favorable remediation performance. The Constructed Wetland Model No.1 (CWM1) embedded in the OpenGeoSys # IPHREEQC was applied to investigate the wetland plant effects on treatment efficiency. Two fundamental functions of the plant roots (i) the radial oxygen loss (ROL) and (ii) exudation of internal organic carbon (IOC), are developed and implemented in the model to simulate the treating processes of planted laboratory-scale IVF wetlands fed by the synthetic wastewater. The good agreement between simulated results and measurements of the planted IVF wetland and the unplanted filters mimicking wetland demonstrates the combined effects of ROL and IOC and the model reliability. In summer the ammonia (NH4-N) and total nitrogen (TN) removals are high as above 90% in both IVF wetlands, and in winter they decline significantly to around 55% and 45% in unplanted wetland, contrastively to about 85% and 78% in the planted wetland. The nitrogen removal - COD/N ratio relation curves of IVF wetlands are proposed and obtained by modeling to evaluate organic carbon loading status. Based on the curves, the COD/N ratios of unplanted and planted wetlands are about 3∼7 and 3∼10 gCOD/gN for high TN removal respectively. Planted wetlands can tolerate a wider range of COD/N ratio influents than unplanted ones. The ROL in the unplanted wetland promotes COD and NH4-N removal, while may inhibit denitrification under low-temperature conditions. The single addition of IOC enhances the oxygen-consuming and restrains the nitrification under the full loaded COD condition. Summing up all organic carbon releases from substrate and roots as IOC, the quantification of IOC acts on nitrogen treatment was simulated and compared with the external organic carbon (EOC) loading from influent. IOC performs higher efficiency on TN removal than EOC at the same organic loading rates. The results provide the thoughts of the solution for low TN removal in the carbon deficient constructed wetlands.

Environmental Earth Sciences Progress Report 2020 and Outlook 2021.

The present editorial 2020 continues the series of status reports in Environmental Earth Sciences (EES) in previous years 2017 and 2019 (Kolditz et al. in Environ Earth Sci 77: 8, 2018, Kolditz et al. in Environ Earth Sci 79: 11, 2020). The year 2020 coming to an end was heavily influenced by the COVID-19 pandemic affecting all areas of life including research work and, therefore, scientific publishing as well ("Introduction"). One bright spot which shows longevity of journals that produce a quality product is that Environmental Earth Sciences (EES) is celebrating its 45th anniversary of publication. To this extent EES continues the tradition to honor the most cited papers contributing to the 2020 Impact Factor (IF) ("Highly and most cited topics") and provide information on the current status of EES as well as an outlook to 2021 ("Progress report").

The risks of long-term re-injection in supercritical geothermal systems.

Supercritical geothermal systems are appealing sources of sustainable and carbon-free energy located in volcanic areas. Recent successes in drilling and exploration have opened new possibilities and spiked interest in this technology. Experimental and numerical studies have also confirmed the feasibility of creating fluid conducting fractures in sedimentary and crystalline rocks at high temperature, paving the road towards Enhanced Supercritical Geothermal Systems. Despite their attractiveness, several important questions regarding safe exploitation remain open. We dedicate this manuscript to the first thermo-hydro-mechanical numerical study of a doublet geothermal system in supercritical conditions. Here we show that thermally-induced stress and strain effects dominate the geomechanical response of supercritical systems compared to pore pressure-related instabilities, and greatly enhance seismicity during cold water re-injection. This finding has important consequences in the design of Supercritical Geothermal Systems.

The brittle-ductile transition in active volcanoes.

Contrasting deformation mechanisms precede volcanic eruptions and control precursory signals. Density increase and high uplifts consistent with magma intrusion and pressurization are in contrast with dilatant responses and reduced surface uplifts observed before eruptions. We investigate the impact that the rheology of rocks constituting the volcanic edifice has on the deformation mechanisms preceding eruptions. We propose a model for the pressure and temperature dependent brittle-ductile transition through which we build a strength profile of the shallow crust in two idealized volcanic settings (igneous and sedimentary basement). We have performed finite element analyses in coupled thermo-hydro-mechanical conditions to investigate the influence of static diking on the local brittle-ductile transition. Our results show that in active volcanoes: (i) dilatancy is an appropriate indicator for the brittle-ductile transition; (ii) the predicted depth of the brittle-ductile transition agrees with the observed attenuated seismicity; (iii) seismicity associated with diking is likely to be affected by ductile deformation mode caused by the local temperature increase; (iv) if failure occurs within the edifice, it is likely to be brittle-dilatant with strength and stiffness reduction that blocks stress transfers within the volcanic edifice, ultimately damping surface uplifts.

Tensor Field Visualization using Fiber Surfaces of Invariant Space.

Scientific visualization developed successful methods for scalar and vector fields. For tensor fields, however, effective, interactive visualizations are still missing despite progress over the last decades. We present a general approach for the generation of separating surfaces in symmetric, second-order, three-dimensional tensor fields. These surfaces are defined as fiber surfaces of the invariant space, i.e. as pre-images of surfaces in the range of a complete set of invariants. This approach leads to a generalization of the fiber surface algorithm by Klacansky et al. [16] to three dimensions in the range. This is due to the fact that the invariant space is three-dimensional for symmetric second-order tensors over a spatial domain. We present an algorithm for surface construction for simplicial grids in the domain and simplicial surfaces in the invariant space. We demonstrate our approach by applying it to stress fields from component design in mechanical engineering.

Effect of solution pH on aging dynamics and surface structural evolution of mZVI particles: H2 production and spectroscopic/microscopic evidence.

A microscale zero-valent iron (mZVI)-based in situ reactive zone is a promising technology for contaminated groundwater remediation. Estimation of mZVI aging behavior after its injection into the subsurface is essential for efficiency and longevity assessments. In this study, batch tests were conducted to investigate the effect of initial pH on mZVI aging dynamics, as well as the formation and evolution of aging products over 112 days. Results indicated that mZVI aging accelerated with decreasing initial pH. Corrosion rates of mZVI particles under pH 6.0 and 7.5 were approximately two orders of magnitude higher than those observed at pH 9.0. The morphological, structural, and compositional evolution of mZVI particles in three systems (pH = 6.0, 7.5, and 9.0) were investigated using scanning electron microscopy, X-ray diffraction, and X-ray photoelectron spectroscopy. In acidic and neutral solutions, a thick passivation layer with loosely and unevenly distributed aging precipitates was observed, and Fe3O4 was the final aging precipitate. Nevertheless, in alkaline solutions, minute aging precipitates were detected on the mZVI surface at 112 day. Characterization results suggested that mZVI was oxidized via the Fe0-Fe(OH)2-Fe3O4 route. These findings shed new light on mZVI aging mechanisms, particularly its physicochemical characteristics and the structural evolution of mZVI in field-scale groundwater remediation applications.

Development and application of a novel method for regional assessment of groundwater contamination risk in the Songhua River Basin.

The main objective of this study is to quantify the groundwater contamination risk of Songhua River Basin by applying a novel approach of integrating public datasets, web services and numerical modelling techniques. To our knowledge, this study is the first to establish groundwater risk maps for the entire Songhua River Basin, one of the largest and most contamination-endangered river basins in China. Index-based groundwater risk maps were created with GIS tools at a spatial resolution of 30arc sec by combining the results of groundwater vulnerability and hazard assessment. Groundwater vulnerability was evaluated using the DRASTIC index method based on public datasets at the highest available resolution in combination with numerical groundwater modelling. As a novel approach to overcome data scarcity at large scales, a web mapping service based data query was applied to obtain an inventory for potential hazardous sites within the basin. The groundwater risk assessment demonstrated that <1% of Songhua River Basin is at high or very high contamination risk. These areas were mainly located in the vast plain areas with hotspots particularly in the Changchun metropolitan area. Moreover, groundwater levels and pollution point sources were found to play a significantly larger impact in assessing these areas than originally assumed by the index scheme. Moderate contamination risk was assigned to 27% of the aquifers, predominantly associated with less densely populated agricultural areas. However, the majority of aquifer area in the sparsely populated mountain ranges displayed low groundwater contamination risk. Sensitivity analysis demonstrated that this novel method is valid for regional assessments of groundwater contamination risk. Despite limitations in resolution and input data consistency, the obtained groundwater contamination risk maps will be beneficial for regional and local decision-making processes with regard to groundwater protection measures, particularly if other data availability is limited.

Distinct kinetics and mechanisms of mZVI particles aging in saline and fresh groundwater: H2 evolution and surface passivation.

Application of microscale zero-valent iron (mZVI) is a promising technology for in-situ contaminated groundwater remediation; however, its longevity is negatively impacted by surface passivation, especially in saline groundwater. In this study, the aging behavior of mZVI particles was investigated in three media (milli-Q water, fresh groundwater and saline groundwater) using batch experiments to evaluate their potential corrosion and passivation performance under different field conditions. The results indicated that mZVI was reactive for 0-7 days of exposure to water and then gradually lost H2-generating capacity over the next hundred days in all of the tested media. In comparison, mZVI in saline groundwater exhibited the fastest corrosion rate during the early phase (0-7 d), followed by the sharpest kinetic constant decline in the latter phases. The SEM-EDS and XPS analyses demonstrated that in the saline groundwater, a thin and compact oxide film was immediately formed on the surface and significantly shielded the iron reactive site. Nevertheless, in fresh groundwater and milli-Q water, a passive layer composed of loosely and unevenly distributed precipitates slowly formed, with abundant reactive sites available to support continuous iron corrosion. These findings provide insight into the molecular-scale mechanism that governs mZVI passivation and provide implications for long-term mZVI application in saline contaminated groundwater.

Remediation of trichloroethylene-contaminated groundwater by three modifier-coated microscale zero-valent iron.

Building a microscale zero-valent iron (mZVI) reaction zone is a promising in situ remediation technology for restoring groundwater contaminated by trichloroethylene (TCE). In order to determine a suitable modifier that could not only overcome gravity sedimentation of mZVI but also improve its remediation efficiency for TCE, the three biopolymers xanthan gum (XG), guargum (GG), and carboxymethyl cellulose (CMC) were employed to coat mZVI for surface modification. The suspension stability of the modified mZVI and its TCE removal efficiency were systematically investigated. The result indicated that XG as a shear-thinning fluid showed the most remarkable efficiency of preventing mZVI from gravity sedimentation and enhancing the TCE removal efficiency by mZVI. In a 480-h experiment, the presence of XG (3 g L(-1)) increased the TCE removal efficiency by 31.85 %, whereas GG (3 g L(-1)) and CMC (3 g L(-1)) merely increased by 15.61 and 9.69 % respectively. The pH value, Eh value, and concentration of ferrous ion as functions of the reaction time were recorded in all the reaction systems, which indicated that XG worked best in buffering the pH value of the solution and inhibiting surface passivation of mZVI.

Transport and retention of xanthan gum-stabilized microscale zero-valent iron particles in saturated porous media.

Microscale zero valent iron (mZVI) is a promising material for in-situ contaminated groundwater remediation. However, its usefulness has been usually inhibited by mZVI particles' low mobility in saturated porous media for sedimentation and deposition. In our study, laboratory experiments, including sedimentation studies, rheological measurements and transport tests, were conducted to investigate the feasibility of xanthan gum (XG) being used as a coating agent for mZVI particle stabilization. In addition, the effects of XG concentration, flow rate, grain diameter and water chemistry on XG-coated mZVI (XG-mZVI) particle mobility were explored by analyzing its breakthrough curves and retention profiles. It was demonstrated that XG worked efficiently to enhance the suspension stability and mobility of mZVI particles through the porous media as a shear thinning fluid, especially at a higher concentration level (3 g/L). The results of the column study showed that the mobility of XG-mZVI particles increased with an increasing flow rate and larger grain diameter. At the highest flow rate (2.30 × 10(-3) m/s) within the coarsest porous media (0.8-1.2 mm), 86.52% of the XG-mZVI flowed through the column. At the lowest flow rate (0.97 × 10(-4) m/s) within the finest porous media (0.3-0.6 mm), the retention was dramatically strengthened, with only 48.22% of the particles flowing through the column. The XG-mZVI particles appeared to be easily trapped at the beginning of the column especially at a low flow rate. In terms of two representative water chemistry parameters (ion strength and pH value), no significant influence on XG-mZVI particle mobility was observed. The experimental results suggested that straining was the primary mechanism of XG-mZVI retention under saturated condition. Given the above results, the specific site-related conditions should be taken into consideration for the design of a successful delivery system to achieve a compromise between maximizing the radius of influence of the injection and minimizing the injection pressure.

MEVA--An Interactive Visualization Application for Validation of Multifaceted Meteorological Data with Multiple 3D Devices.

BACKGROUND: To achieve more realistic simulations, meteorologists develop and use models with increasing spatial and temporal resolution. The analyzing, comparing, and visualizing of resulting simulations becomes more and more challenging due to the growing amounts and multifaceted character of the data. Various data sources, numerous variables and multiple simulations lead to a complex database. Although a variety of software exists suited for the visualization of meteorological data, none of them fulfills all of the typical domain-specific requirements: support for quasi-standard data formats and different grid types, standard visualization techniques for scalar and vector data, visualization of the context (e.g., topography) and other static data, support for multiple presentation devices used in modern sciences (e.g., virtual reality), a user-friendly interface, and suitability for cooperative work. METHODS AND RESULTS: Instead of attempting to develop yet another new visualization system to fulfill all possible needs in this application domain, our approach is to provide a flexible workflow that combines different existing state-of-the-art visualization software components in order to hide the complexity of 3D data visualization tools from the end user. To complete the workflow and to enable the domain scientists to interactively visualize their data without advanced skills in 3D visualization systems, we developed a lightweight custom visualization application (MEVA - multifaceted environmental data visualization application) that supports the most relevant visualization and interaction techniques and can be easily deployed. Specifically, our workflow combines a variety of different data abstraction methods provided by a state-of-the-art 3D visualization application with the interaction and presentation features of a computer-games engine. Our customized application includes solutions for the analysis of multirun data, specifically with respect to data uncertainty and differences between simulation runs. In an iterative development process, our easy-to-use application was developed in close cooperation with meteorologists and visualization experts. The usability of the application has been validated with user tests. We report on how this application supports the users to prove and disprove existing hypotheses and discover new insights. In addition, the application has been used at public events to communicate research results.

Mechanism insights into enhanced trichloroethylene removal using xanthan gum-modified microscale zero-valent iron particles.

This report focuses on the enhancement in trichloroethylene (TCE) removal from contaminated groundwater using xanthan gum (XG)-modified, microscale, zero-valent iron (mZVI). Compared with bare mZVI, XG-coated mZVI increased the TCE removal efficiency by 30.37% over a 480-h experimental period. Because the TCE removal is attributed to both sorption and reduction processes, the contributions from sorption and reduction were separately investigated to determine the mechanism of XG on TCE removal using mZVI. The results showed that the TCE sorption capacity of mZVI was lower in the presence of XG, whereas the TCE reduction capacity was significantly increased. The FTIR spectra confirmed that XG, which is rich in hydrophilic functional groups, was adsorbed onto the iron surface through intermolecular hydrogen bonds, which competitively repelled the sorption and mass transfer of TCE toward reactive sites. The variations in the pH, Eh, and Fe(2+) concentration as functions of the reaction time were recorded and indicated that XG buffered the solution pH, inhibited surface passivation, and promoted TCE reduction by mZVI. Overall, the XG-modified mZVI was considered to be potentially effective for the in-situ remediation of TCE contaminated groundwater due to its high stability and dechlorination reactivity.

Description and verification of a novel flow and transport model for silicate-gel emplacement.

We present a novel approach for the numerical simulation of the gelation of silicate solutions under density-dependent flow conditions. The method utilizes an auxiliary, not density-dependent solute that is subject to a linear decay function to provide temporal information that is used to describe the viscosity change of the fluid. By comparing the modeling results to experimental data, we are able to simulate the behavior and the gelation process of the injected solute for three different compositions, including long-term stability of the gelated area, and non-gelation of low concentrations due to hydro-dynamic dispersion. This approach can also be used for other types of solutes with this gelling property and is useful in a variety of applications in geological, civil and environmental engineering.

Enhanced biodegradation by hydraulic heterogeneities in petroleum hydrocarbon plumes.

In case of dissolved electron donors and acceptors, natural attenuation of organic contaminant plumes in aquifers is governed by hydrodynamic mixing and microbial activity. Main objectives of this work were (i) to determine whether aerobic and anaerobic biodegradation in porous sediments is controlled by transverse dispersion, (ii) to elucidate the effect of sediment heterogeneity on mixing and biodegradation, and (iii) to search for degradation-limiting factors. Comparative experiments were conducted in two-dimensional sediment microcosms. Aerobic toluene and later ethylbenzene degradation by Pseudomonas putida strain F1 was initially followed in a plume developing from oxic to anoxic conditions and later under steady-state mixing-controlled conditions. Competitive anaerobic degradation was then initiated by introduction of the denitrifying strain Aromatoleum aromaticum EbN1. In homogeneous sand, aerobic toluene degradation was clearly controlled by dispersive mixing. Similarly, under denitrifying conditions, microbial activity was located at the plume's fringes. Sediment heterogeneity caused flow focusing and improved the mixing of reactants. Independent from the electron accepting process, net biodegradation was always higher in the heterogeneous setting with a calculated efficiency plus of 23-100% as compared to the homogeneous setup. Flow and reactive transport model simulations were performed in order to interpret and evaluate the experimental results.

Determination of first-order degradation rate constants from monitoring networks.

In this article, different strategies for estimating first-order degradation rate constants from measured field data are compared by application to multiple, synthetic, contaminant plumes. The plumes were generated by numerical simulation of contaminant transport and degradation in virtual heterogeneous aquifers. These sites were then individually and independently investigated on the computer by installation of extensive networks of observation wells. From the data measured at the wells, that is, contaminant concentrations, hydraulic conductivities, and heads, first-order degradation rates were estimated by three 1D centerline methods, which use only measurements located on the plume axis, and a two-dimensional method, which uses all concentration measurements available downgradient from the contaminant source. Results for both strategies show that the true rate constant used for the numerical simulation of the plumes in general tends to be overestimated. Overestimation is stronger for narrow plumes from small source zones, with an average overestimation factor of about 5 and single values ranging from 0.5 to 20, decreasing for wider plumes, with an average overestimation factor of about 2 and similar spread. Reasons for this overestimation are identified in the velocity calculation, the dispersivity parameterization, and off-centerline measurements. For narrow plumes, the one- and the two-dimensional strategies show approximately the same amount of overestimation. For wider plumes, however, incorporation of all measurements in the two-dimensional approach reduces the estimation error. No significant relation between the number of observation wells in the monitoring network and the quality of the estimated rate constant is found for the two-dimensional approach.

Influence of small-scale heterogeneities on contaminant transport in fractured crystalline rock.

We present a sequence of purely advective transport models that demonstrate the influence of small-scale geometric inhomogeneities on contaminant transport in fractured crystalline rock. Special weight is placed on the role of statistically generated variable fracture apertures. The fracture network geometry and the aperture distribution are based on information from an in situ radionuclide retardation experiment performed at Grimsel test site (Swiss Alps). The obtained breakthrough curves are fitted with the advection dispersion equation and continuous-time random walks (CTRW). CTRW is found to provide superior fits to the late-arrival tailing and is also found to show a good correlation with the velocity distributions obtained from the hydraulic models. The impact of small-scale heterogeneities, both in fracture geometry and aperture, on transport is shown to be considerable.

Uncertainty assessment of contaminant plume length estimates in heterogeneous aquifers.

The Virtual Aquifer approach is used in this study to assess the uncertainty involved in the estimation of contaminant plume lengths in heterogeneous aquifers. Contaminant plumes in heterogeneous two-dimensional conductivity fields and subject to first order and Michaelis-Menten (MM) degradation kinetics are investigated by the center line method. First order degradation rates and plume lengths are estimated from point information obtained along the plume center line. Results from a Monte-Carlo investigation show that the estimated rate constant is highly uncertain and biased towards overly high values. Uncertainty and bias amplify with increasing heterogeneity up to maximum values of one order of magnitude. Calculated plume lengths reflect this uncertainty and bias. On average, plume lengths are estimated to about 50% of the true plume length. When plumes subject to MM degradation kinetics are investigated by using a first order rate law, an additional error is introduced and uncertainty as well as bias increase, causing plume length estimates to be less than 40% of the true length. For plumes with MM degradation kinetics, therefore, a regression approach is used which allows the determination of the MM parameters from center line data. Rate parameters are overestimated by a factor of two on average, while plume length estimates are about 80% of the true length. Plume lengths calculated using the MM parameters are thus closer to the correct length, as compared to the first order approximation. This approach is therefore recommended if field data collected along the center line of a plume give evidence of MM kinetics.

Numerical simulation of reactive processes in an experiment with partially saturated bentonite.

Bentonites are preferred materials for use as engineered barriers for high-level nuclear waste repositories. Simulation of geochemical processes in bentonite is therefore important for long-term safety assessment of those repositories. In this work, the porewater chemistry of a bentonite sample subject to simultaneous heating and hydration, as studied by Cuevas et al. [Cuevas, J., Villar, M., Fernández, A., Gómez, P., Martín, P., 1997. Porewaters extracted from compacted bentonite subjected to simultaneous heating and hydration. Applied Geochemistry 12, 473-481.], was assessed with a non-isothermal reactive transport model by coupling the geochemical software PHREEQC2 with the object-oriented FEM simulator GeoSys/RockFlow. Reactive transport modelling includes heat transport, two-phase flow, multicomponent transport and geochemical reactions in the liquid phase, i.e. ion exchange, mineral dissolution/precipitation and equilibrium reactions. Simulations show that the easily soluble minerals in bentonite determine the porewater chemistry. Temperature affects both two-phase flow and geochemical reactions. Porosity change due to dissolution/precipitation is low during the experiment. However, changes of the effective porosity caused by bentonite swelling can be very large. The simulated results agree well with the experimental data.