Simulating Woodchip Bioreactor Performance Using a Dual-Porosity Model

Amended Vegetation Filters as Nature-Based Solutions for the Treatment of Pharmaceuticals: Infiltration Experiments Coupled to Reactive Transport Modelling

In small populations and scattered communities, wastewater treatment through vegetation filters (VFs), a nature-based solution, has proved to be feasible, especially for nutrient and organic matter removal. However, the presence of pharmaceuticals in wastewater and their potential to infiltrate through the vadose zone and reach groundwater is a drawback in the evaluation of VF performances. Soil amended with readily labile carbon sources, such as woodchips, enhances microbial activity and sorption processes, which could improve pharmaceutical attenuation in VFs. The present study aims to assess if woodchip amendments to a VF's soil are able to abate concentrations of selected pharmaceuticals in the infiltrating water by quantitatively describing the occurring processes through reactive transport modelling. Thus, a column experiment using soil collected from an operating VF and poplar woodchips was conducted, alongside a column containing only soil used as reference. The pharmaceuticals acetaminophen, naproxen, atenolol, caffeine, carbamazepine, ketoprofen and sulfamethoxazole were applied daily to the column inlet, mimicking a real irrigation pattern and periodically measured in the effluent. Ketoprofen was the only injected pharmaceutical that reached the column outlet of both systems within the experimental timeframe. The absence of acetaminophen, atenolol, caffeine, carbamazepine, naproxen and sulfamethoxazole in both column outlets indicates that they were attenuated even without woodchips. However, the presence of 10,11-epoxy carbamazepine and atenolol acid as transformation products (TPs) suggests that incomplete degradation also occurs and that the effect of the amendment on the infiltration of TPs is compound-specific. Modelling allowed us to generate breakthrough curves of ketoprofen in both columns and to obtain transport parameters during infiltration. Woodchip-amended columns exhibited Kd and μw values from one to two orders of magnitude higher compared to soil column. This augmentation of sorption and biodegradation processes significantly enhanced the removal of ketoprofen to over 96%.

Climate change effects on denitrification performance of woodchip bioreactors treating agricultural tile drainage

Denitrifying woodchip bioreactors (WBRs) are a nature-based technology that are increasingly used to control nonpoint source nitrate (NO3-) pollution in agricultural catchments. The treatment effectiveness of WBRs depends on temperature and hydraulic retention time (HRT), both of which are affected by climate change. Warmer temperatures will increase microbial denitrification rates, but the extent to which the resulting benefits to treatment performance may be offset by intensified precipitation and shorter HRTs is not clear. Here, we use three years of monitoring data from a WBR in Central New York State to train an integrated hydrologic-biokinetic model describing links among temperature, precipitation, bioreactor discharge, denitrification kinetics, and NO3- removal efficiencies. Effects of climate warming are assessed by first training a stochastic weather generator with eleven years of weather data from our field site, and then adjusting the distribution of precipitation intensities according to the Clausius-Clapeyron relationship between water vapor and temperature. Modeling results indicate, in our system, faster denitrification rates will outweigh the influence of intensified precipitation and discharge under warming, leading to net improvements in NO3- load reductions. Median cumulative NO3- load reductions at our study site from May - October are projected to increase from 21.7% (interquartile range 17.4%-26.1%) under baseline hydro-climate to 41.0% (interquartile range 32.6-47.1%) with a + 4 °C change in mean air temperature. This improved performance under climate warming is driven by strong nonlinear dependence of NO3- removal rates on temperature. Temperature sensitivity may increase with woodchip age and lead to stronger temperature-response in systems like this one with a highly aged woodchip matrix. While the impacts of hydro-climatic change on WBR performance will depend on site-specific properties, this hydrologic-biokinetic modeling approach provides a framework for assessing climate impacts on the effectiveness of WBRs and other denitrifying nature-based systems.

Improving accuracy of quantifying nitrate removal performance and enhancing understanding of processes in woodchip bioreactors using high-frequency data

Woodchip bioreactors have gained popularity in many countries as a conservation practice for reducing nitrate load to freshwater. However, current methods for assessing their performance may be inadequate when nitrate removal rates (RR) are determined from low-frequency (e.g., weekly) concurrent sampling at the inlet and outlet. We hypothesised that high-frequency monitoring data at multiple locations can help improve the accuracy of quantifying nitrate removal performance, enhance the understanding of processes occurring within a bioreactor, and therefore improve the design practice for bioreactors. Accordingly, the objectives of this study were to compare RRs calculated using high- and low-frequency sampling and assess the spatiotemporal variability of the nitrate removal within a bioreactor to unravel the processes occurring within a bioreactor. For two drainage seasons, we monitored nitrate concentrations at 21 locations on an hourly or two-hourly basis within a pilot-scale woodchip bioreactor in Tatuanui, New Zealand. A novel method was developed to account for the variable lag time between entry and exit of a parcel of sampled drainage water. Our results showed that this method not only enabled lag time to be accounted for but also helped quantify volumetric inefficiencies (e.g., dead zone) within the bioreactor. The average RR calculated using this method was significantly higher than the average RR calculated using conventional low-frequency methods. The average RRs of each of the quarter sections within the bioreactor were found to be different. 1-D transport modelling confirmed the effect of nitrate loading on the removal process as nitrate reduction followed Michaelis-Menten (MM) kinetics. These results demonstrate that high-frequency temporal and spatial monitoring of nitrate concentrations in the field allows improved description of bioreactor performance and better understanding of processes occurring within woodchip bioreactors. Thus, insights gained from this study can be used to optimise the design of future field bioreactors.

Field evaluation of four denitrifying woodchip bioreactors for nitrogen removal in eastern South Dakota, United States

Edge-of-field practices such as denitrifying woodchip bioreactors can be used to improve the water quality of agricultural effluents. This study evaluated the effectiveness of four field-scale woodchip bioreactors in removing nitrate‑nitrogen (nitrate-N) from subsurface drainage in eastern South Dakota. Four woodchip bioreactors were installed and monitored between 2014 and 2016 near Arlington, Baltic, Hartford, and Montrose, South Dakota. Results showed that reduction in nitrate-N concentration for the four bioreactors ranged from 7 % to 100 %, corresponding to removal rates of 5 to 27 g N/m³/day for the four bioreactors during the study period. Average Nitrate-N load reduction in the four bioreactors studied ranged from 39 % to 89 % during the study period. Reduction of nitrate-N in the four bioreactors decreased, on average, by 30 % when temperature dropped below 12 °C during the study period. Flow rate and hydraulic retention time (HRT) also influenced nitrate-N removal in the bioreactors as samples collected immediately following rainfall events showed high nitrate-N load removal compared to samples collected later after the rainfall events during the study period.

The use of recycled materials in a biofilter to polish anammox wastewater treatment plant effluent

Anammox is gaining popularity for treating wastewater containing high-strength ammonia due to lower energy demand compared to conventional nitrification-denitrification processes; however, anammox is reported to increase nitrate loads in the effluent. The objective of this study was to assess the applicability of recycled materials [recycled concrete aggregate (RCA) and rice husks (RH)] as a polishing step to improve anammox reactor effluent quality. Anammox effluents were separately passed through two single-stage columns containing RCA and RH, and one two-stage column (50% RCA, 50% RH) to quantify total N, ammonia, nitrate, nitrite, and phosphate removal efficiencies. Langmuir isotherm experiments were conducted to quantify nitrate, nitrite, and phosphate sorption capacities in the columns. The RCA column exhibited the highest phosphate sorption capacity (0.074 mg/g), while the RH column exhibited higher nitrite and nitrate adsorption (0.063 mg/g and 0.023 mg/g respectively). We created a Hydrus-1D model to estimate pseudo-first-order reaction rates in the columns. Because RCA media can form metal-phosphate precipitates, the fastest phosphate reaction rate (1.58 min-1) occurred in the RCA column. The two-stage column demonstrated the greatest overall removals for all nutrients, and removal rates were consistent throughout the experimental period. The two-stage column achieved 15% total N, 94% ammonia-N, 38% nitrate-N, 75% nitrite-N, and 27% phosphate removal. The maximum nitrite, nitrate, and phosphate adsorption capacities in the two-stage column were 0.030 mg/g, 0.017 mg/g, and 0.014 mg/g respectively. This is the first study to demonstrate that recycled materials can successfully be integrated into a biofilter as an effluent polishing step to remove nutrients from anammox wastewater.

Denitrifying bioreactor microbiome: Understanding pollution swapping and potential for improved performance

Denitrifying woodchip bioreactors are a best management practice to reduce nitrate-nitrogen (NO₃ -N) loading to surface waters from agricultural subsurface drainage. Their effectiveness has been proven in many studies, although variable results with respect to performance indicators have been observed. This paper serves the purpose of synthesizing the current state of the science in terms of the microbial community, its impact on the consistency of bioreactor performance, and its role in the production of potential harmful by-products including greenhouse gases, sulfate reduction, and methylmercury. Microbial processes other than denitrification have been observed in these bioreactor systems, including dissimilatory nitrate reduction to ammonia (DNRA) and anaerobic ammonium oxidation (anammox). Specific gene targets for denitrification, DNRA, anammox, and the production of harmful by-products are identified from bioreactor studies and other environmentally relevant systems for application in bioreactor studies. Lastly, cellulose depletion has been observed over time via increasing ligno-cellulose indices, therefore, the microbial metabolism of cellulose is an important function for bioreactor performance and management. Future work should draw from the knowledge of soil and wetland ecology to inform the study of bioreactor microbiomes.

Denitrification using permeable reactive barriers with organic substrate or zero-valent iron fillers: controlling mechanisms, challenges, and future perspectives

Nitrate as a diffusive agricultural contaminant has been causing substantial groundwater quality deterioration worldwide. In situ groundwater remediation techniques using permeable reactive barriers (PRBs) have attracted increasing interest. Particularly, PRBs based on biological denitrification, using the organic substrate as a biostimulator, and chemical nitrate reduction, using zero-valent iron (ZVI) as a reductant, are two major PRB approaches for groundwater denitrification. This review paper analyzed the published studies over the past 10 years (2010-2020) using laboratory, modeling, and field-scale approaches to explore the performance and mechanisms of these two types of PRBs. Important factors affecting the denitrification efficiencies as well as the influential mechanisms were discussed. Several research gaps have been identified and further research needs are discussed in the end.

Forecasting leachate generation from pilot woodchip stockpiles using a three-dimensional transient flow model

Leachate generation from open stockpiles of recycled woodchip materials is potentially harmful to aquatic ecosystems. There is growing interest in using numerical models to simulate leachate generation from outdoor piles, but this requires information about the hydraulic properties of the materials. The objectives of this study were to simulate leachate from woodchip piles with the numerical model HYDRUS-3D and to optimize subsets of parameters for single (SPM) and dual (DPM) pore flow models with the Bayesian Markov Chain Monte Carlo algorithm DREAM_ZS. Three experimental piles, each approximately 30 m³, were setup with mixtures of either once (coarse) or twice (fine) ground woodchips. Leachate continuously collected over a period of six months was similar across piles. As a result, subsets of optimized flow parameters for the coarse and fine woodchips were not different. Leachate predictions by the two pore flow models were similar and agreed reasonably with the field measurements, as indicated by Nash-Sutcliffe efficiency values greater than 0.6. This result suggests the simpler SPM is adequate for field predictions of leachate. However, leachate was consistently under-predicted by both pore models by 13-27% during rainfall events with more than 1 cm in 6 h. The optimized flow models can be used as a tool for studying pile management strategies.

Nitrous Oxide and Methane Dynamics in Woodchip Bioreactors: Effects of Water Level Fluctuations on Partitioning into Trapped Gas Phases

Woodchip bioreactors (WBRs) are low-cost, passive systems for nonpoint source nitrogen removal at terrestrial-aquatic interfaces. The greenhouse gases nitrous oxide (N₂O) and methane (CH₄) can be produced within WBRs, and efforts to reduce N₂O and CH₄ emissions from WBR systems require improved understanding of the biogeochemical and physical-chemical mechanisms regulating their production, transport, and release. This study evaluates the impact of trapped gas-filled void volumes as sinks of dissolved gases from water and as sources of episodic fluxes when water levels fall. Dissolved gas tracer experiments in a laboratory bioreactor were used to parameterize nonequilibrium advection-dispersion-gas transfer models and quantify trapping of gas-filled voids as a function of antecedent hydrological conditions. Experiments following a water-level rise revealed that up to 24% of the WBR pore volume was occupied by trapped gas phases, which were primarily located in pore spaces inside woodchips. This finding was confirmed with X-ray-computed microtomography. N₂O (3.3-10%) and CH₄ (4.3-14%) injected into the reactor following a water table rise partitioned into gas-filled voids and were released when water tables fell. In the case of N₂O, partitioning into trapped gas phases makes N₂O unavailable for enzymatic reduction, potentially enhancing N₂O fluxes under fluctuating water levels.

Reducing adverse side effects by seasonally lowering nitrate removal in subsurface flow constructed wetlands

Subsurface flow constructed wetlands with wood chips (SSF-CWs) have proven to effectively reduce the loss of nitrogen (N) from agricultural fields to surface water, however in some cases production of negative side effects such as methane and phosphate occur. We examined if these side effects can be avoided by decreasing the hydraulic retention time (HRT) from on average 82 h to 11 h during summer to autumn in two pilot SSF-CWs. Furthermore, we investigated the potential of the SSF-CWs to reduce phosphorus (P) loss from agricultural drainage systems. The influent and effluent concentration of total N (TN), nitrate-N, total P, phosphate-P, suspended sediment, and sulphate were monitored for five years (2013-2017). Methane concentrations were measured during two periods in 2014 and 2017. Flow was measured continuously by electromagnetic flow-meters. The nitrate-N removal was reduced from 98-100% to 27-32% and the sulphate reduction from 32-53% to 1-2% when decreasing HRT. Concurrently this resulted in a considerable decrease in the difference between the effluent and influent concentration of phosphate-P and methane concentration compared to similar periods in the preceding years. The SSF-CWs retained 67-85% of the annual loading of particulate P, but acted as both a sink and source of phosphate-P, thus further initiatives are therefore required to prevent phosphate-P release from SSF-CWs. Although during the entire monitoring period the SSF-CWs retained 29-33% of the total P loading. In summary, this study stresses how important a holistic approach is when implementing and designing new N mitigation measures.

Evaluation of Mechanistic Models for Nitrate Removal in Woodchip Bioreactors

Woodchip bioreactors (WBRs) are increasingly being applied to remove nitrate from runoff. In this study, replicate columns with aged woodchips were subjected to a range of measured flow rates and influent nitrate concentrations with an artificial stormwater matrix. Dissolved oxygen (DO), nitrate, and dissolved organic carbon (DOC) were measured along the length of the columns. A multispecies reactive transport model with Michaelis-Menten kinetics was developed to explain the concentration profiles of DO, nitrate, and DOC. Four additional models were developed based on simplifying assumptions, and all five models were tested for their ability to predict nitrate concentrations in the experimental columns. Global sensitivity analysis and constrained optimization determined the set of parameters that minimized the root-mean-squared error (RMSE) between the model and the experimental data. A k-fold validation test revealed no statistical difference in RMSE for predicting nitrate concentrations between a zero-order model and the other multispecies reactive transport models tested. Additionally, the multispecies reactive transport models demonstrated no significant differences in predicting DO and DOC concentrations. These results suggest that denitrification in an aged woodchip bioreactor at constant temperature can effectively be modeled using zero-order kinetics when nitrate concentrations are >2 mg-N L^-1. A multispecies model may be used if predicting DOC or DO concentrations is desired.

Moving Denitrifying Bioreactors beyond Proof of Concept: Introduction to the Special Section

Denitrifying bioreactors are organic carbon-filled excavations designed to enhance the natural process of denitrification for the simple, passive treatment of nitrate-nitrogen. Research on and installation of these bioreactors has accelerated within the past 10 years, particularly in watersheds concerned about high nonpoint-source nitrate loads and also for tertiary wastewater treatment. This special section, inspired by the meeting of the Managing Denitrification in Agronomic Systems Community at the 2014 Annual Meeting of the American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America, aims to firmly establish that denitrifying bioreactors for treatment of nitrate in drainage waters, groundwater, and some wastewaters have moved beyond the proof of concept. This collection of 14 papers expands the peer-reviewed literature of denitrifying bioreactors into new locations, applications, and environmental conditions. There is momentum behind the pairing of wood-based bioreactors with other media (biochar, corn cobs) and in novel designs (e.g., use within treatment trains or use of baffles) to broaden applicability into new kinds of waters and pollutants and to improve performance under challenging field conditions such as cool early season agricultural drainage. Concerns about negative bioreactor by-products (nitrous oxide and hydrogen sulfide emissions, start-up nutrient flushing) are ongoing, but this translates into a significant research opportunity to develop more advanced designs and to fine tune management strategies. Future research must think more broadly to address bioreactor impacts on holistic watershed health and greenhouse gas balances and to facilitate collaborations that allow investigation of mechanisms within the bioreactor "black box."