Denitrification potential in the subsurface environment in the Manawatu River catchment, New Zealand: Indications from oxidation-reduction conditions, hydrogeological factors, and implications for nutrient management

Coping with groundwater pollution in high-nitrate leaching areas: The efficacy of denitrification

The Ningxia Yellow River irrigation area, characterized by an arid climate and high leaching of NO3--N, exhibits complex and unique groundwater nitrate (NO3--N) pollution, with denitrification serving as the principal mechanism for NO3--N removal. The characteristics of N leaching from paddy fields and NO3--N removal by groundwater denitrification were investigated through a two-year field observation. The leaching losses of total nitrogen (TN) and NO3--N accounted for 10.81-27.34% and 7.59-12.74%, respectively, of the N input. The linear relationship between NO3--N leaching and N input indicated that the fertilizer-induced emission factor (EF) of NO3--N leaching in direct dry seeding and seedling-raising and transplanting paddy fields was 8.2% (2021, R2 = 0.992) and 6.7% (2022, R2 = 0.994), respectively. The study highlighted that the quadratic relationship between the NO3--N leaching loss and N input (R2 = 0.999) significantly outperformed the linear relationship. Groundwater denitrification capacity was characterized by monitoring the concentrations of dinitrogen (N2) and nitrous oxide (N2O). The results revealed substantial seasonal fluctuations in excess N2 and N2O concentrations in groundwater, particularly following fertilization and irrigation events. The removal efficiency of NO3--N via groundwater denitrification ranged from 42.70% to 74.38%, varying with depth. Groundwater denitrification capacity appeared to be linked to dissolved organic carbon (DOC) concentration, redox conditions, fertilization, irrigation, and soil texture. The anthropogenic-alluvial soil with limited water retention accelerated the leaching of NO3--N into groundwater during irrigation. This process enhances the groundwater recharge capacity and alters the redox conditions of groundwater, consequently impacting groundwater denitrification activity. The DOC concentration emerged as the primary constraint on the groundwater denitrification capacity in this region. Hence, increasing carbon source concentration and enhancing soil water retention capacity are vital for improving the groundwater denitrification capacity and NO3--N removal efficiency. This study provides practical insights for managing groundwater NO3--N pollution in agricultural areas, optimizing fertilization strategies and improving groundwater quality.

Determining the likelihood and cost of detecting reductions of nitrate‑nitrogen concentrations in groundwater across New Zealand

Nitrate‑nitrogen (NO3-N) is a contaminant of concern in groundwater worldwide. Stakeholders need information on the ability to detect changes in NO3-N concentrations to prove that land management practices are meeting water quality aims. We created a database of quarterly to monthly NO3-N measurements in 948 sites across New Zealand; 186 of those sites had mean residence time (MRT) data. New Zealand has set a target of sufficient land use mitigations in the next 30 years to ensure steady state surface water concentrations do not exceed 2.4 mg L-1. Here we assess whether the current monitoring network could identify the impacts of these mitigations, assuming that the mitigations are successfully implemented at the source. Only 41 % of the network could detect statistically significant reductions with the current standard quarterly sampling after 30 years of monitoring. The percentage of sites increased to 60 % with increased monitoring frequency (often weekly) but this required a 100-300 % increase in monitoring costs. However, policy makers and stakeholders typically require information on policy and mitigation effectiveness within 5-10 years. Detection within 5-10 years was very unlikely (0-20 % of sites) regardless of the sampling frequency. Importantly, these analyses include the impacts of groundwater lag and temporal dispersion on the likelihood of detecting change, ignoring these impacts, incorrectly, yields a much higher likelihood of detecting reductions. We conclude that the current monitoring network is unlikely to be fit for the purpose of detecting NO3-N reductions within practical timeframes or budgets. Furthermore, we conclude that lag and temporal dispersion effects must be included in detection power calculations; we therefore recommend that MRT data is regularly collected. We also provide a python package to enable easy detection power calculations with lag and temporal dispersion impacts, thereby supporting the development of robust change-detection monitoring networks.

Assessing the environmental factors affecting the sustainability of Aini Falaj system

This study investigates the spatial distribution patterns and environmental factors influencing the Aini Falaj system in a specific study area. The research findings are presented through the lens of the following four categories: collinearity diagnostics, spatial autocorrelation analysis, kernel density (KD) findings, and multivariate geographically weighted regression (MGWR) analysis. The collinearity diagnostics were applied to examine the interrelationships among 18 independent environmental variables. The results indicate the absence of significant multicollinearity concerns, with most variables showing values below the critical threshold of five for variance inflation factors (VIFs). The selected variables indicate minimal intercorrelation, suggesting that researchers should be confident utilizing them in subsequent modelling or regression analyses. A spatial autocorrelation analysis using Moran's Index revealed positive spatial autocorrelation and significant clustering patterns in the distribution of live and non-functional Aini Falajs. High concentrations of live or dead Falajs tended to be surrounded by neighbouring areas with similar characteristics. These findings provide insights into the ecological preferences and habitat associations of Aini Falajs, thereby aiding conservation strategies and targeted studies. The kernel density (KD) analysis depicted distribution patterns of live and dry Aini Falajs through hotspots and cold spots. Specific regions exhibited high-density areas of live Falajs, indicating favourable environmental conditions or historical factors contributing to their concentrated distribution. Identifying these high-density zones can enhance our understanding of the spatial patterns and potential factors influencing the prevalence and sustainability of Aini Falajs. The multivariate geographically weighted regression (MGWR) models revealed strong associations between the live or dead status of Aini Falajs and environmental factors. The precipitation, topographic wetness index (TWI), aspect and slope exerted positive impacts on the live status, while evaporation, solar radiation, distance to drains and drain density exerted negative influences. Similar associations were observed for the dead status, emphasising the importance of controlling evaporation, shading mechanisms, proper drainage planning and sustainable land-use practices. This study provides valuable insights into the spatial distributions and factors influencing the live and dead status of Aini Falajs, thereby contributing to our understanding of their ecological dynamics and guiding conservation efforts and management strategies.

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.

Quantifying MCPA load pathways at catchment scale using high temporal resolution data

Detection of the agricultural acid herbicide MCPA (2-methyl-4-chlorophenoxyacetic acid) in drinking water source catchments is of growing concern, with economic and environmental implications for water utilities and wider ecosystem services. MCPA is poorly adsorbed to soil and highly mobile in water, but hydrological pathway processes are relatively unknown at the catchment scale and limited by coarse resolution data. This understanding is required to target mitigation measures and to provide a framework to monitor their effectiveness. To address this knowledge gap, this study reports findings from river discharge and synchronous MCPA concentration datasets (continuous 7 hour and with additional hourly sampling during storm events) collected over a 7 month herbicide spraying season. The study was undertaken in a surface (source) water catchment (384 km²-of which 154 km² is agricultural land use) in the cross-border area of Ireland. Combined into loads, and using two pathway separation techniques, the MCPA data were apportioned into event and baseload components and the former was further separated to quantify a quickflow (QF) and other event pathways. Based on the 7 hourly dataset, 85.2 kg (0.22 kg km^-2 by catchment area, or 0.55 kg km^-2 by agricultural area) of MCPA was exported from the catchment in 7 months. Of this load, 87.7 % was transported via event flow pathways with 72.0 % transported via surface dominated (QF) pathways. Approximately 12 % of the MCPA load was transported via deep baseflows, indicating a persistence in this delayed pathway, and this was the primary pathway condition monitored in a weekly regulatory sampling programme. However, overall, the data indicated a dominant acute, storm dependent process of incidental MCPA loss during the spraying season. Reducing use and/or implementing extensive surface pathway disconnection measures are the mitigation options with greatest potential, the success of which can only be assessed using high temporal resolution monitoring techniques.

Upscaling of Surface Water and Groundwater Interactions in Hyporheic Zone from Local to Regional Scale

The groundwater (GW) and surface water (SW) interaction (SW-GW) through the hyporheic zone is a significant component in sustainable water resource management. The complexities in SW-GW interactions increase from a local to a regional scale and are affected by variation in hydraulic, hydrologic, and hydrogeologic (3H) processes. Controlling factors and their upscaling of these processes to assess SW-GW interaction have not been addressed sufficiently in previous studies. Additionally, it is unclear what the effective factors are at different scales during the upscaling. Therefore, the present review focused on controlling factors of 3H processes in SW-GW interaction and their upscaling techniques. Relevancy of controlling factors was identified at different scales. Applications of different approaches and their uncertainties were also discussed for the characterization of SW-GW interactions. The study revealed that the improved data from different approaches is crucial for machine learning training and its application in the SW and GW assessment at local, sub-catchment, and catchment scales. Based on the outcomes, a framework has been proposed to execute modalities of controlling factors using remote sensing, geophysics, and artificial intelligence. The proposed framework could help in handling big data and accurate upscaling for water resource management.

Contrasting subsurface denitrification characteristics under temperate pasture lands and its implications for nutrient management in agricultural catchments

Subsurface denitrification plays a key role in the reduction or 'attenuation' of nitrate contamination of groundwater and surface waters. We investigated subsurface denitrification characteristics in the vadose zone and shallow groundwater at four sites under pastoral farming in the Manawatū River catchment, located in the lower part of North Island, New Zealand. The denitrification potential of the vadose zone was determined by the laboratory incubation assays measuring the denitrifying enzyme activity (DEA) in soil samples collected from different soil horizons (up to 2.1 m below ground surface), whereas denitrification rates in shallow groundwaters were measured in situ by single-well, push-pull tests conducted in piezometers installed at multiple depths at the study sites. Soils and underlying geology, defining hydrogeologic settings, appear to influence the spatial variability of subsurface denitrification characteristics at the study sites. Where the vadose zone is thin and composed of coarse-textured soils, percolation of nitrate was evident in observed high nitrate-nitrogen concentrations (>5 mg L^-1) in oxic and young shallow groundwaters, but low nitrate-nitrogen concentrations (<0.05 mg L^-1) were observed in the reduced shallow groundwater found underneath the fine textured soils and/or a thick vadose zone. This was confirmed by the push-pull tests measuring denitrification rates from 0.08 to 1.07 mg N L^-1 h^-1 in the reduced shallow groundwaters (dissolved oxygen or DO < 0.5 mg L^-1), while negligible in the oxic groundwaters (DO > 5 mg L^-1) found at the study sites. These contrasting subsurface denitrification characteristics determine the ultimate delivery of nitrate losses from agricultural soils to receiving waters, where the fine textured thick vadose zone and reducing groundwater conditions offer nitrate reduction in the subsurface environment.

Soil Salinity and Moisture Control the Processes of Soil Nitrification and Denitrification in a Riparian Wetlands in an Extremely Arid Regions in Northwestern China

Soil nitrification and denitrification are key nitrogen (N) removal processes in riparian wetlands in extremely arid regions, but the driving factors of the two N processes in these wetlands are still unclear. We measured soil nitrification and denitrification rates and related environmental properties in a typical riparian wetland in the middle reaches of the Heihe River, northwestern China. Our results showed that rates of soil nitrification and denitrification exhibited moderate variability, ranging from 52.77 to 221.18 μg kg−1 h−1 and 91.25 to 428.26 μg kg−1 h−1, respectively. Soil salinity was high, with mean electrical conductivity (EC) of 6.8 mS cm−1. Soil salinity and moisture were the key factors influencing nitrification and denitrification in this riparian wetland in an extremely arid region. Soil salinity exerted significant inhibitory impact on soil nitrification when EC was > 4.05 mS cm−1. Soil nitrification increased with an increase in soil moisture when soil water content < 27.03% and decreased with an increase in soil moisture when soil water content > 27.03%. Denitrification had a significantly negative relationship with soil salinity, and significantly positive relationship with soil moisture. The interaction of soil salinity and moisture played a central role in regulating soil denitrification. Based on these results, we propose that water consumption of riparian wetlands, and the planting of halophytes, should be increased to reduce soil salinity and increase soil moisture, which is essential for sustaining soil N removal function in riparian wetlands in extremely arid regions.

Spatial Distribution of Integrated Nitrate Reduction across the Unsaturated Zone and the Groundwater Body in Germany

Nitrate pollution in groundwater and its mitigation strategies is currently a topic of controversial debate in Germany, and the demand for harmonised approaches for the implementation of regulations is increasing. Important factors that need to be considered when planning mitigation measures are the nitrogen inputs into water bodies and the natural nitrate reduction capacity. The present study introduces a nationwide, harmonised and simplified approach for estimating nitrate reduction as an integral quantity across the unsaturated zone and the groundwater body. The nitrate reduction rates vary from 0% to 100%, and are on average 57%, with high values in the north of Germany and low values in the south. Hydrogeological characteristics are associated with the estimated nitrate reduction rates, whereby the influence of aquifer type and redox conditions are particularly relevant. The nitrate reduction rates are substantially higher in porous aquifers and under anaerobic conditions than in fractured, consolidated aquifers and under aerobic conditions. This contribution presents a harmonised conceptual approach to derive the nitrate reduction rate at a 1 km × 1 km resolution. This information can be used when planning and designing mitigation measures to meet the groundwater nitrate limits.

Controls on species distribution and biogeochemical cycling in nitrate-contaminated groundwater and surface water, southeastern Australia

A detailed study of groundwater and surface water nitrate over four seasons across an area of varied landuse provided insights into the mechanisms that underlie accumulation and transport of nitrate. High nitrate concentrations found in a significant percentage of surface water and shallow groundwater samples are due to anthropogenic contamination. Statistics (PCA, ANOVA, parsimonious model and general linear regression) were used to explore the relationship between NO₃^- and land use, and confirmed that areas of high NO₃^- concentration are associated with dairy pasture and horticulture. Seasonally, NO₃^- levels are greater during winter, the wettest part of the year. Values of δ¹⁵N showed that most nitrate is sourced from livestock waste, with a smaller contribution from synthetic fertilizer. Direct wash-off of animal waste from dairy farms results in higher NO₃^- concentrations in surface water than in groundwater. Denitrification is an important NO₃^- attenuation mechanism which reduces NO₃^- to NH₄, as demonstrated by the PCA analysis, which showed positive correlation of NO₃^- concentrations with dissolved oxygen and negative correlations with NH₄⁺, Fe²⁺and Mn²⁺; the latter two species may act as the electron donors necessary for reduction of NO₃^-. The often high NO₃^- concentrations in shallow groundwater are decreased by denitrification, which can occur at relatively shallow depths (<3 m). The relatively small NO₃^- concentrations in deeper groundwater are due partly to denitrification, but more to originally lower NO₃^- concentrations, as the age of deeper groundwater shows that it was recharged before agriculture was established in the study area. Overall, the study demonstrates the usefulness of hydrogeochemical characterisation and multivariate statistics in the evaluation of impacts of agricultural land-use on regional N cycling. In particular, the results show that efforts to mitigate NO₃^- pollution from farms should concentrate more on wash-off of animal waste than the contribution of nitrogenous synthetic fertilizer.

Achieving unbiased predictions of national-scale groundwater redox conditions via data oversampling and statistical learning

An important policy consideration for integrated land and water management is to understand the spatial distribution of nitrate attenuation in the groundwater system, for which redox condition is the key indicator. This paper proposes a methodology to accommodate the computational demands of large datasets, and presents national-scale predictions of groundwater redox class for New Zealand. Our approach applies statistical learning methods to relate the redox class determined on groundwater samples to spatially varying attributes. The trained model uses these spatial variables to predict redox status in areas without sample data. We assembled the groundwater sample data from regional authority databases, and assigned each sample a redox class. A key achievement was to overcome the influence of sample selection bias on model training via oversampling. We removed additional bias imposed by imbalances in the predictor variables by applying a conditional inference random forest classifier. The unbiased trained model uses eight predictors, and achieves a high validation performance (accuracy 0.81, kappa 0.71), providing good confidence in model predictions. National maps are provided for redox class and probability at specified depths. Feature importance rankings indicate that reducing conditions are associated with poorly-drained soils, and to a lesser extent, high hydrological variability, low elevation, and low-permeability lithology. These conditions are common in New Zealand's coastal and lowland plains, where artificial drainage is required to make land suitable for production. The spatial extent of reduced groundwater increases with depth, suggesting a shallow influence of soil infiltration or mobile organic carbon, and a deeper influence of lithological electron donors. Our model provides unbiased predictions at a scale relevant for environmental policy development and legislation. Identifying where the ecosystem service provided by denitrification can be utilised will enable spatially targeted interventions that can achieve the desired environmental outcome in a more cost-effective manner than non-targeted interventions.

Water-quality issues facing dairy farming: potential natural and built attenuation of nitrate losses in sensitive agricultural catchments

Context Dairy farming will be increasingly scrutinised for its environmental impacts, in particular for its impacts on freshwater quality in New Zealand and elsewhere. Management and mitigation of high nitrate losses is one of the greatest water-quality challenges facing dairy farming in New Zealand and other countries. Management of critical flow pathways and nitrate-attenuation capacity could offer potential solutions to this problem and help maintain dairy-farming productivity, while reducing its water-quality impacts. Aims The present paper reviewed the key water-quality issues faced by dairy farming and assessed potential of emerging edge-of-paddock technologies, and catchment-scale nutrient-attenuation practices, to reduce nitrate losses from dairy farming to receiving water bodies. Methods We developed a conceptual catchment-scale modelling analysis assessing potential natural and built attenuation of nitrate losses from dairy farming in the Tararua and Rangitikei catchments (located in the lower part of the North Island, New Zealand). Key results This exploratory analysis suggests that a reduction of greater than 25% in the river nitrate loads from dairy-farming areas could potentially be achieved by spatially aligning dairy land with areas of high subsurface nitrate-attenuation capacity, and by managing critical flow pathways using innovative edge-of-field technologies such as controlled drainage, drainage-water harvesting for supplemental irrigation, woodchip bioreactors, and constructed wetlands in the study catchments. Conclusions The research findings highlighted the potential to better understand, map and effectively utilise existing natural and new built-in nitrate-attenuation capacity to significantly reduce water-quality impacts from dairy farming across environmentally sensitive agricultural catchments. This knowledge and tools could help farmers close the gap between what can be achieved with current, in-field mitigation practises and the nitrogen-loss allocation imposed by regulatory authorities. Implications However, the research findings presented here are based on a coarse-scale, conceptual modelling analysis, and therefore further research is recommended to develop tools and practices to better understand, map and effectively utilise existing natural and new built-in nitrogen attenuation capacity at farm-scale to achieve productive and environmentally friendly pastoral dairy farming across agricultural landscapes.

3D characterization of the subsurface redox architecture in complex geological settings

Nitrogen (N) leaching caused by agricultural activities is one of the major threats to the aquatic ecosystems and public health. Moving from the agricultural soils through the subsurface and reemerging to the surface water, N undergoes various biogeochemical reactions along pathways in the subsurface, which occur heterogeneously in space and time. Thus to improve our understanding on the fate and distribution of N in the aquatic environment, detailed knowledge about the subsurface hydrogeological and biogeochemical conditions, especially the redox conditions, are essential. In this study, 3D information of the redox conditions termed the redox architecture was investigated in two Danish catchments with intensive agriculture underlain by glacial deposits. Towed transient electromagnetic (tTEM) resistivity was interpreted which reveals the subsurface geological structures at a few hectare scale. These geophysical data were integrated with sediment and water chemistry for the redox architecture interpretations. The top soils of both catchments are characterized as clay-till, but the tTEM showed that the subsurface hydrogeological structures are distinctively different. We identified three types of redox architectures in the studied catchments: 1) a planar redox architecture with a single redox interface; 2) a geological-window redox architecture with local complexity; and 3) a glaciotectonic-thrusted redox architecture with high complexity. The baseflow N load at the catchment outlets reflect the contributions of N via oxic pathways through the complex redox architectures of the subsurface. We conclude that in some landscapes, the redox architecture cannot be simplified as a single interface that roughly follows the terrain; hence, thorough investigations of the structural heterogeneity of the local redox architectures will be necessary to improve simulations of N evolution along pathways and quantifications of N attenuation under various mitigation scenarios.

The influence of unsaturated zone drainage status on denitrification and the redox succession in shallow groundwater

Since nitrate is a major agricultural freshwater contaminant, denitrification is the environmentally most important step in the ecological succession of redox processes that can occur in groundwater. Understanding where and to what extent denitrification occurs would enable spatially differentiated land management and regulation. We investigated in a dairy farming catchment in the North Island of New Zealand the influence of the unsaturated zone's drainage status on the redox succession in shallow groundwater along a well transect spanning drainage conditions from well drained to very poorly drained. Groundwater samples were analysed for a variety of parameters including nitrate, tritium, dinitrogen, argon, methane and nitrous oxide. The redox classification based on measured redox-sensitive parameters broadly matched the a priori assessed drainage status of the overlying unsaturated zone. Only the groundwater underlying the well-drained soil was oxic and reflected the N losses from the intensive pastoral land use, with nitrate nitrogen concentrations up to 9.6 mg L^-1. All other sites had mildly to strongly reduced groundwater and concomitantly decreasing or low nitrate concentrations, even at the water table. The tritium-derived mean residence time (MRT) estimates for the oxic groundwater (12 and 14 y) were within the range found in mildly reduced groundwater from the imperfectly drained sites (6-24 y), with the exception of one sample from below an aquitard (105 y). In contrast, the strongly reduced groundwater observed at the poorly and very poorly drained sites was relatively immobile (55 to >110 y). Denitrification was confirmed by the nitrate dual isotope signatures, and by the occurrence of excess dinitrogen, and likely occurred in both the unsaturated and saturated zones. A coherent sampling scheme throughout the unsaturated zone - saturated zone continuum should be used in future studies to allow ascertaining the exact location of denitrification activity in vertical profiles.

Management of slurry in Gran Canaria Island with full-scale natural treatment systems for wastewater (NTSW). One year experience in livestock farms

The aim of this work is to describe the performance of three full-scale natural treatment systems for wastewater, which operated in an integrated manner in livestock pig farms (1000-1500pigsintotal) over one year. Slurry management was performed with these natural treatment systems operating under the normal waste loading conditions of the livestock farms in which were integrated. The systems were comprised of elements such as first generation digesters, subsurface flow constructed wetlands and facultative ponds. The facilities, located on the island of Gran Canaria (Spain), enabled the study of viable alternatives for effluent management characterized by low-cost treatments. The systems were evaluated in terms of chemical oxygen demand removal efficiency, operating with variable organic loading. Values of between 80% and 90% were obtained. A comparison was also made of first-generation cascade flow digester operation (<70% removal efficiency), with complete-mix digesters (<20% removal efficiency), and finally with facultative ponds combined with subsurface flow constructed wetlands (<91% removal efficiency). It was also verified that when natural treatment systems for wastewater combine different elements they have better removal efficiency and better response to load and/or flow changes.

Do Wet-Dry Ratio and Fe-Mn System Affect Oxidation-Reduction Potential Nonlinearly in the Subsurface Wastewater Infiltration Systems?

To understand characteristics of on-line oxidation-reduction potential (ORP) in a subsurface wastewater infiltration system (SWIS) under different intermittent influent conditions, ORP among five matrix depths at wet-dry ratios (R_wds) of 2:1, 1:1 and 1:2 with a hydraulic load of 0.10 m³·(m²·d)^-1 were monitored. Results showed that the optimal R_wd for the SWIS was 1:1. In that case, ORP at 40 and 65 cm depths changed significantly, by 529 mV and 261 mV, respectively, from the inflow period to the dry period, which was conducive to the recovery of the oxidation environment. It was concluded that ORP varied nonlinearly in strongly aerobic and hypoxic environment. Wastewater was fed into the SWIS at 80 cm and dissolved oxygen diffused at the initial period of one cycle. As a consequence, ORP at 65 cm increased with water content increasing. However, ORP at 40 and 95 cm displayed inverse trends. Moreover, results showed that ORP decreased with Fe²⁺ and Mn²⁺ increasing under aerobic conditions (p < 0.05) because Fe²⁺ and Mn²⁺ moved with wastewater flow. Effluent met reuse requirements and no clogging was found in the SWIS during the operation.

Vertical stratification of redox conditions, denitrification and recharge in shallow groundwater on a volcanic hillslope containing relict organic matter

Natural denitrification in groundwater systems has been recognised as an ecosystem service that reduces the impact of agriculturally-derived nitrate inputs to surface waters. Identification of this ecosystem service within the landscape would permit spatially differentiated land management and legislation. However, spatial variation in groundwater redox conditions poses a significant challenge to such a concept. To gain understanding of the small-scale mosaic of biogeochemical and hydrological controls on denitrification, we established a well field consisting of 11 multilevel well (MLW) clusters on a hillslope containing relict organic matter buried by volcanic deposits 1.8 ka before present. Based on site-specific redox classification thresholds, vertical redox gradients and denitrification potentials were detected at 7 of the 11 sites. Palaeosols or woody debris, which had previously been identified in laboratory experiments as resident electron donors fuelling denitrification, were visually recognisable at 4 of the 7 MLW sites with vertical redox gradients. Moderately enhanced groundwater dissolved organic carbon (DOC) concentrations occurred where resident electron donors were evident. DOC concentrations were lower where anoxic and nitrate-depleted groundwater was found but with an absence of resident electron donors. In these instances, it was assumed that nitrate reduction had occurred somewhere upgradient of the sampled well screen along the lateral groundwater flow path, with the proximate electron donor (DOC) largely consumed in the process, since no evidence was found for denitrification being fuelled by inorganic electron donors. Due to high variability in the isotopic signature of nitrate in oxidised groundwater, the nitrate dual isotope method did not yield firm evidence for denitrification. However, realistic vertical patterns were obtained using the excess N₂ method. Tritium-based age dating revealed that oxic conditions were restricted to young groundwater (mean residence time ≤ 3 y), while anoxic conditions were observed across a wider age range (3-25 y).

Transport and potential attenuation of nitrogen in shallow groundwaters in the lower Rangitikei catchment, New Zealand

Intensive agricultural activities are generally associated with nitrogen leaching from agricultural soils, and this nitrogen has the potential to percolate and contaminate groundwater and surface waters. We assessed surface water and groundwater interactions, and nitrogen leaching and its potential attenuation in shallow groundwater in the lower Rangitikei River catchment (832km²), New Zealand. We combined regional- and local-scale field surveys and experiments, nutrient budget modelling, and hydraulic and geochemical methods, to gain an insight into leaching, transformation and transport of nitrogen via groundwaters to the river in the study area. Concurrent river flow gaugings (in January 2015) and a piezometric map, developed from measured depths to groundwater in 110 bores (in October 2014), suggest groundwater discharges to the Rangitikei River in the upper parts of the study area, while there is groundwater recharge near the coast. The groundwater redox characterisation, based on sampling and analysis of 15 mostly shallow bores (<30m below ground level (bgl)), suggests groundwater across the lower Rangitikei catchment in general is under anoxic/reduced conditions. The groundwater typically has low dissolved oxygen (DO) concentrations (<1mg/L), suggesting the subsurface environment is conducive to potential attenuation by 'denitrification' of NO₃-N in groundwater. We further measured NO₃-N attenuation in shallow groundwater piezometers (3-6mbgl) using single-well push-pull tests. We found generally low levels (<0.5mg/L) of NO₃-N in shallow groundwater piezometers (>5mbgl), despite being installed under intensive land uses, such as dairying and cropping. Our in-field push-pull tests showed NO₃-N reduction at four shallow groundwater piezometers, with the rates of reduction varying from 0.04mgNL^-1h^-1 to 1.57mgNL^-1h^-1. This highlights the importance of a sound understanding of not only the sources, but also transport and transformation, or fate, of nutrients leached from farms, to mitigate the likely impacts of land use on water quality and ecosystem health in agricultural catchments.