A Comparison of Modeling Approaches in Simulating Chlorinated Ethene Removal in a Constructed Wetland by a Microbial Consortia

The purpose of this study is to compare different approaches to modeling the reductive dechlorination of chlorinated ethenes in the anaerobic region of an upward flow constructed wetland by microbial consortia. A controlled simulation experiment that compares three different approaches to modeling the degradation of chlorinated ethenes in wetland environments is conducted and investigates how each of the modeling approaches affect simulation results. Concepts like microbial growth in the form of a biofilm and spatially varying contaminant concentrations bring the validity of the CSTR assumption into question. These concepts are incorporated into the different modeling approaches to evaluate the CSTR assumption. Model simulations show that spatially varying contaminant concentrations have a significant effect on contaminant effluent concentrations. Additionally, the significance of the incorporation of a biofilm concept depends on the time characteristics of both diffusive mass transport and reaction kinetics

Landfill Leachate Production and Gas Generation Numerical Model

Numerous processes occur in landfills which lend themselves to modeling. Many of the processes are mutually interdependent. An unsteady numerical model is developed combining the major processes. The three-dimensional moisture transport equations and boundary conditions are solved using an implicit finite difference scheme. The boundaries are determined through a two-dimensional runoff model for the landfill surface and a one-dimensional leachate liner flow model at the bottom of the landfill. The runoff model accounts for evapotranspiration, runoff, infiltration, and leachate recirculation. Richard\u27s equation is solved for saturated and unsaturated vertical flows and Darcy\u27s Equation is solved for lateral flow between adjacent saturated landfill cells. Results of the moisture flow are used to solve contaminant production and transport equations. Contaminant production uses moisture flow and previous leaching history to generate source terms. The source terms and recirculated contaminants are used to implicitly solve contaminant transport equations which account for advection, diffusion, and dispersion of the contaminant. Landfill temperatures are predicted by solving an energy equation implicitly. Temperatures are combined with moisture content and gas production history to determine gas generation. The model is applied to three Wisconsin lysimeters and a Kentucky landfill to demonstrate the simulation of leachate and contaminant production and transport. Comparison to the HELP water balance model is also done for a Wisconsin lysimeter. The model is also applied to an existing landfill to demonstrate the gas generation portions of the model

Arsenic Speciation and Phytoremediation Modeling for Environmental Management

Arsenic has been used throughout recorded history but during the industrial revolution widespread use led to global environmental impact. The two forms that should be considered in environment management are arsenate and arsenite. The calculations of environmental risk for arsenic exposure relies the toxicity of arsenite however, in well aeriated surface soils arsenate may be the predominate form. Ecological risk assessments based on arsenite studies will lead to restrictive remediation requirements that do not adequately reflect the level of risk. Arsenate resembles phosphate and as such has a greater affinity for phytoremediation. Phytoremediation is one of the most viable and cost effective cleanup techniques developed. Different mathematical approaches have been implemented to characterize phytoremediation systems to address concerns with performance. A system dynamic model is presented to describe solute transport in groundwater coupled to sorption by plant roots, translocation into plant stems, and evapotranspiration. The model was tested and assessed using published and peer-reviewed experimental data, to assess its capability to mimic phytoremediation processes. The model is consistent with previous research establishing the extraction process as a constringent factor for this cleanup technique. The model included modules that can estimate rainfall, seasonal temperature and growth. The modules allow for the independent verification of data before input into the model. The implementation of phytoremediation model can provide information about: pollutant-media-plant interaction, pollutant concentration and flow rate through the plant

Is the Hyporheic Zone Relevant beyond the Scientific Community?

Rivers are important ecosystems under continuous anthropogenic stresses. The hyporheic zone is a ubiquitous, reactive interface between the main channel and its surrounding sediments along the river network. We elaborate on the main physical, biological, and biogeochemical drivers and processes within the hyporheic zone that have been studied by multiple scientific disciplines for almost half a century. These previous efforts have shown that the hyporheic zone is a modulator for most metabolic stream processes and serves as a refuge and habitat for a diverse range of aquatic organisms. It also exerts a major control on river water quality by increasing the contact time with reactive environments, which in turn results in retention and transformation of nutrients, trace organic compounds, fine suspended particles, and microplastics, among others. The paper showcases the critical importance of hyporheic zones, both from a scientific and an applied perspective, and their role in ecosystem services to answer the question of the manuscript title. It identifies major research gaps in our understanding of hyporheic processes. In conclusion, we highlight the potential of hyporheic restoration to efficiently manage and reactivate ecosystem functions and services in river corridors. View Full-Tex

Modeling contaminant transport and fate and subsequent impacts on ecosystems

Assessing risks associated with the release of metals into the environment and managing remedial activities requires simulation tools that depict speciation and risk with accurate mechanistic models and well-defined transport parameters. Such tools need to address the following processes: (1) aqueous speciation, (2) distribution mechanisms, (3) transport, and (4) ecological risk. The primary objective of this research is to develop a simulation tool that accounts for these processes. Speciation in the aqueous phase can be assessed with geochemical equilibrium models, such as MINEQL+. Furthermore, metal distribution can be addressed mechanistically. Studies with Pb sorption to amorphous aluminum (HAG), iron (HFO), and manganese (HMO) oxides, as well as oxide coatings, demonstrated that intraparticle diffusion is the rate-limiting mechanism in the sorption process, where best-fit surface diffusivities ranged from 10-18 to 10-15 cm2 s-1 Intraparticle surface diffusion was incorporated into the Groundwater Modeling System (GMS) to accurately simulate metal contaminant mobility where oxides are present. In the model development, the parabolic concentration layer approximation and the operator split technique were used to solve the microscopic diffusion equation coupled with macroscopic advection and dispersion. The resulting model was employed for simulating Sr90 mobility at the U.S. Department of Energy (DOE) Hanford Site. The Sr90 plume is observed to be migrating out of the 100-N area extending into other areas of the Hanford Site and beyond. Once bioavailability is understood, static or dynamic ecological risk assessments can be conducted. Employing the ERA model, a static ecological risk assessment for exposure to depleted uranium (DU) at Aberdeen and Yuma Proving Grounds (APG and YPG) revealed that a reduction in plant root weight is considered likely to occur. For most terrestrial animals at YPG, the predicted DU dose is less than that which would result in a decrease in offspring. However, for the lesser long-nosed bat, reproductive effects are expected to occur through the reduction in size and weight of offspring. At APG, based on very limited data, it is predicted that uranium uptake will not likely affect survival of terrestrial animals and aquatic species. In model validation, sampling of pocket mice, kangaroo rat, white-throated woodrat, deer, and milfoil showed that body burden concentrations fall into the distributions simulated at both sites. This static risk assessment provides a solid background for applying the dynamic approach. Overall, this research contributes to a holistic approach in developing accurate mechanistic models for simulating metal contaminant mobility and bioavailability in subsurface environments

Dynamic p-enrichment schemes for multicomponent reactive flows

We present a family of p-enrichment schemes. These schemes may be separated into two basic classes: the first, called \emph{fixed tolerance schemes}, rely on setting global scalar tolerances on the local regularity of the solution, and the second, called \emph{dioristic schemes}, rely on time-evolving bounds on the local variation in the solution. Each class of $p$-enrichment scheme is further divided into two basic types. The first type (the Type I schemes) enrich along lines of maximal variation, striving to enhance stable solutions in "areas of highest interest." The second type (the Type II schemes) enrich along lines of maximal regularity in order to maximize the stability of the enrichment process. Each of these schemes are tested over a pair of model problems arising in coastal hydrology. The first is a contaminant transport model, which addresses a declinature problem for a contaminant plume with respect to a bay inlet setting. The second is a multicomponent chemically reactive flow model of estuary eutrophication arising in the Gulf of Mexico.Comment: 29 pages, 7 figures, 3 table

Modeling Vertical Flow Treatment Wetland Hydraulics to Optimize Treatment Efficiency

An upward Vertical Flow Treatment Wetland (uVFTW) at Wright Patterson AFB designed to bioremediate contaminated groundwater exhibits hydraulic short-circuiting. Prior studies estimated that groundwater flowed through less than 50% of the wetland’s volume, and that the mean residence time was significantly less than the nominal residence time, which was calculated assuming flow through the entire wetland volume. The objective of this research was to investigate how uVFTW hydraulics affects treatment efficiency, and to propose design strategies to maximize treatment efficiency. A groundwater flow and contaminant transport model of a uVFTW that couples hydraulics and degradation kinetics was built and applied to estimate the effectiveness of engineering solutions aimed at improving treatment efficiency. Model simulations indicate that the engineering solutions improve hydraulic residence times, volumetric utilization, and treatment efficiency over the existing wetland, but also that increasing hydraulic residence time only has a significant impact on treatment efficiency when the time scale for the biodegradation process is similar to the wetland residence time. Degradation kinetics must be quantitatively understood to determine an optimum range for hydraulic residence time, and to ensure that resources are not wasted in an attempt to improve hydraulic performance where no improvement in degradation performance is possible

Modeling water resources management at the basin level: review and future directions

Water quality / Water resources development / Agricultural production / River basin development / Mathematical models / Simulation models / Water allocation / Policy / Economic aspects / Hydrology / Reservoir operation / Groundwater management / Drainage / Conjunctive use / Surface water / GIS / Decision support systems / Optimization methods / Water supply

Investigating the migration of immiscible contaminant fluid flow in homogeneous and heterogeneous aquifers with high-precision numerical simulations

Numerical modeling of the migration of three-phase immiscible fluid flow in variably saturated zones is challenging due to the different behavior of the system between unsaturated and saturated zones. This behavior results in the use of different numerical methods for the numerical simulation of the fluid flow depending on whether it is in the unsaturated or saturated zones. This paper shows that using a high-resolution shock-capturing conservative method to resolve the nonlinear governing coupled partial differential equations of a three-phase immiscible fluid flow allows the numerical simulation of the system through both zones providing a unitary vision (and resolution) of the migration of an immiscible contaminant problem within a porous medium. In particular, using different initial scenarios (including impermeable “lenses” in heterogeneous aquifers), three-dimensional numerical simulation results are presented on the temporal evolution of the contaminant migration following the saturation profiles of the three-phases fluids flow in variably saturated zones. It is considered either light nonaqueous phase liquid with a density less than the water, or dense nonaqueous phase liquid, which has densities greater than the water initially released in unsaturated dry soil. Our study shows that the fate of the migration of immiscible contaminants in variably saturated zones can be accurately described, using a unique mathematical conservative model, with different evolution depending on the value of the system’s physical parameters, including the contaminant density, and accurately tracking the evolution of the sharp (shock) contaminant front