Identification and reliability of microbial aerobic respiration and denitrification kinetics using a single-well push-pull field test

Prediction of post-Darcy flow based on the spatial non-local distribution of hydraulic gradient: Preliminary assessment of wastewater management

Understanding high-velocity pollutant transport dependent on the large hydraulic gradient and/or heterogeneity of the aquifer and criteria for the onset of post-Darcy flow have attracted considerable attention in water resources and environmental engineering applications. In this study, a parameterized model is established based on the equivalent hydraulic gradient (EHG) which affected by spatial nonlocality of nonlinear head distribution due to the inhomogeneity at a wide range of scales. Two parameters relevant to the spatially non-local effect were selected to predict the development of post-Darcy flow. Over 510 sets of laboratory one-dimensional (1-D) steady hydraulic experimental data were used to validate the performance of this parameterized EHG model. The results show that (1) the spatial nonlocal effect of the whole upstream is related to the mean grain size of the medium, and the anomalous variation due to the small grain size implies the existence of the particle size threshold. (2) The parameterized EHG model can effectively capture the nonlinear trend that fails to be described by the traditional local form of nonlinear models, even if the specific discharge stabilizes at the later stages. (3) The Sub-Darcy flow distinguished by the parameterized EHG model can be equated to the post-Darcy flow, and then the criteria for the post-Darcy flow will be strictly distinguished under the premise of determining the hydraulic conductivity. The results of this study facilitate the identification and prediction of high-velocity non-Darcian flow in wastewater management and provide insight into mass transport by advection at the fine-scale.

Degradation of seven pesticides and two metabolites before and during aquifer storage transfer and recovery operation

Degradation of 7 common pesticides (bentazon, boscalid, chloridazon, fluopyram, flutolanil, imidacloprid, and methoxyfenozide) and 2 metabolites of chloridazon (desphenyl-chloridazon, and methyl-desphenyl-chloridazon) was studied in an anoxic and brackish sandy aquifer before and during Aquifer Storage Transfer and Recovery (ASTR) operation. Fresh tile drainage water was injected and stored for later re-use as irrigation water. We hypothesized that electron acceptors (O₂, NO₃), dissolved organic carbon (∼24.7 mg/L), nutrients (NO₃: ∼14.1 mg/L, NH₄: ∼0.13 mg/L, PO₄: ∼5.2 mg/L), and biodegrading bacteria in tile drainage water could stimulate degradation of the pesticides and metabolites (ranging between 0.013 and 10.8 μg/L) introduced in the aquifer. Pesticide degradation was studied at 6 depths in the aquifer using push-pull tests lasting ±18 days before the onset of ASTR operation. Degradation was too limited to quantify and/or could not be assessed because of the potential occurrence of pesticide retardation. Utilizing push-pull tests to obtain degradation constants should only be considered in future studies for non-retarding pesticides with relative low half-lives (here <20 days). During ASTR operation, pesticide degradation was studied at the same depths during 3 storage periods equally spread over 1.5 years of ASTR operation. Overall, trends of degradation were observed, although with relatively high half-lives of at least 53 days. Microbial adaptation of the aquifer and/or bioaugmentation by the injected biodegrading bacteria did not result in enhanced degradation during consecutive storage periods. Operational monitoring data over longer periods and distances yielded half-lives of at least 141 days. The slow degradation mostly agrees with previous studies. The injected tile drainage water composition did therefore not notably stimulate pesticide degradation. The relatively persistent behavior of the studied pesticides/metabolites implies that ASTR abstracted water will have generally high pesticide concentrations, and non-abstracted water may form a contamination risk for the surrounding native brackish groundwater.

A new solution to transient single-well push-pull test with low-permeability non-Darcian leakage effects

An accurate solute transport model is critical to the interpretation of single-well push-pull (SWPP) test. Previous studies of SWPP test generally consider solitary aquifer that is confined by impermeable layers. Also, existing solutions for solute transport in aquifer-aquitard systems only consider the injection phase and over-simplify the flow field by assuming uniformly distributed velocity in the aquitard. In this study, we developed a numerical model with Dirichlet boundary condition for SWPP test affected by leakage described by a low-permeability non-Darcian expression involving a threshold pressure gradient (I₀). Our SWPP test model considered transient flow in multi-phases, which include injection, chase, rest and extraction phases. Finite-difference scheme was adopted to solve the models of flow and solute transport. The results indicate that an increasing hydraulic diffusivity leads to a greater peak value of breakthrough curve (BTC) while a medium with larger grain size results in less estimation error when using steady-state flow model to interpret the transient SWPP test. Additionally, A greater I₀ makes the solute stored in aquitard more difficult to be extracted out due to dispersion dominance, which results in higher BTC values. For the purpose of application, a lumped dimensionless index called the non-Darcian index (NDI) was proposed to quantify the overestimation degree by neglecting leakage, and the underestimation degree by accounting for Darcian leakage, when interpreting the SWPP test with low-permeability non-Darcian leakage. The long-term slope of breakthrough curve coupled with the NDI can be employed to determine the cases in which the low-permeability non-Darcian leakage should be considered.

Reaction chain modeling of denitrification reactions during a push-pull test

Field quantitative estimation of reaction kinetics is required to enhance our understanding of biogeochemical reactions in aquifers. We extended the analytical solution developed by Haggerty et al. (1998) to model an entire 1st order reaction chain and estimate the kinetic parameters for each reaction step of the denitrification process. We then assessed the ability of this reaction chain to model biogeochemical reactions by comparing it with experimental results from a push-pull test in a fractured crystalline aquifer (Ploemeur, French Brittany). Nitrates were used as the reactive tracer, since denitrification involves the sequential reduction of nitrates to nitrogen gas through a chain reaction (NO3(-)→NO2(-)→NO→N2O→N2) under anaerobic conditions. The kinetics of nitrate consumption and by-product formation (NO2(-), N2O) during autotrophic denitrification were quantified by using a reactive tracer (NO3(-)) and a non-reactive tracer (Br(-)). The formation of reaction by-products (NO2(-), N2O, N2) has not been previously considered using a reaction chain approach. Comparison of Br(-) and NO3(-) breakthrough curves showed that 10% of the injected NO3(-) molar mass was transformed during the 12 h experiment (2% into NO2(-), 1% into N2O and the rest into N2 and NO). Similar results, but with slower kinetics, were obtained from laboratory experiments in reactors. The good agreement between the model and the field data shows that the complete denitrification process can be efficiently modeled as a sequence of first order reactions. The 1st order kinetics coefficients obtained through modeling were as follows: k1=0.023 h(-1), k2=0.59 h(-1), k3=16 h(-1), and k4=5.5 h(-1). A next step will be to assess the variability of field reactivity using the methodology developed for modeling push-pull tracer tests.

Simulating MODFLOW-based reactive transport under radially symmetric flow conditions

Radially symmetric flow and solute transport around point sources and sinks is an important specialized topic of groundwater hydraulics. Analysis of radial flow fields is routinely used to determine heads and flows in the vicinity of point sources or sinks. Increasingly, studies also consider solute transport, biogeochemical processes, and thermal changes that occur in the vicinity of point sources/sinks. Commonly, the analysis of hydraulic processes involves numerical or (semi-) analytical modeling methods. For the description of solute transport, analytical solutions are only available for the most basic transport phenomena. Solving advanced transport problems numerically is often associated with a significant computational burden. However, where axis-symmetry applies, computational cost can be decreased substantially in comparison with full three-dimensional (3D) solutions. In this study, we explore several techniques of simulating conservative and reactive transport within radial flow fields using MODFLOW as the flow simulator, based on its widespread use and ability to be coupled with multiple solute and reactive transport codes of different complexity. The selected transport simulators are MT3DMS and PHT3D. Computational efficiency and accuracy of the approaches are evaluated through comparisons with full 2D/3D model simulations, analytical solutions, and benchmark problems. We demonstrate that radial transport models are capable of accurately reproducing a wide variety of conservative and reactive transport problems provided that an adequate spatial discretization and advection scheme is selected. For the investigated test problems, the computational load was substantially reduced, with the improvement varying, depending on the complexity of the considered reaction network.

Submersible microbial fuel cell sensor for monitoring microbial activity and BOD in groundwater: focusing on impact of anodic biofilm on sensor applicability

A sensor, based on a submersible microbial fuel cell (SUMFC), was developed for in situ monitoring of microbial activity and biochemical oxygen demand (BOD) in groundwater. Presence or absence of a biofilm on the anode was a decisive factor for the applicability of the sensor. Fresh anode was required for application of the sensor for microbial activity measurement, while biofilm-colonized anode was needed for utilizing the sensor for BOD content measurement. The current density of SUMFC sensor equipped with a biofilm-colonized anode showed linear relationship with BOD content, to up to 250 mg/L (∼233 ± 1 mA/m(2)), with a response time of <0.67 h. This sensor could, however, not measure microbial activity, as indicated by the indifferent current produced at varying active microorganisms concentration, which was expressed as microbial adenosine-triphosphate (ATP) concentration. On the contrary, the current density (0.6 ± 0.1 to 12.4 ± 0.1 mA/m(2)) of the SUMFC sensor equipped with a fresh anode showed linear relationship, with active microorganism concentrations from 0 to 6.52 nmol-ATP/L, while no correlation between the current and BOD was observed. It was found that temperature, pH, conductivity, and inorganic solid content were significantly affecting the sensitivity of the sensor. Lastly, the sensor was tested with real contaminated groundwater, where the microbial activity and BOD content could be detected in <3.1 h. The microbial activity and BOD concentration measured by SUMFC sensor fitted well with the one measured by the standard methods, with deviations ranging from 15% to 22% and 6% to 16%, respectively. The SUMFC sensor provides a new way for in situ and quantitative monitoring contaminants content and biological activity during bioremediation process in variety of anoxic aquifers.