Abstract: Elevated groundwater nitrate levels are common in drylands, often in excess of WHO guidelines, with concern for human and animal health. In light of recent attempts to identify nitrate sources in the Kalahari this paper presents the first unsaturated zone (USZ) nitrate profiles and recharge rate estimates for the important transboundary Stampriet Basin, alongside the first rainfall chemistry records. Elevated subsurface nitrate reaches 100–250 and 250–525 mg/L NO3–N, with NO3–N/Cl of 4–12, indicating input above evapotranspiration. Chloride mass balance recharge rates range from 4 to 27 mm/y, indicating a vertical movement of these nitrate pulses toward the water table over multi-decadal timescales. These profiles are sampled from dune crests, away from high concentrations of animals and without termite mounds. Given low-density animal grazing is unlikely to contribute consistent spot-scale nitrate over decades, these profiles give an initial estimate of naturally-produced concentrations. This insight is important for the management of the Stampriet Basin and wider Kalahari groundwater. This study expands our knowledge about elevated nitrate in dryland USZs, demonstrating that it can occur as pulses, probably in response to transient vegetation cover and that it is not limited to long-residence time USZs with very limited downward moisture flux (recharge).

Abstract: In coastal salinized groundwater systems, contamination from various nitrate (NO3) inputs combined with complex hydrogeochemical processes make it difficult to distinguish NO3 sources and identify potential NO3 transformtation processes. Effective field-based NO3 studies in coastal areas are needed to improve the understanding of NO3 contamination dynamics in groundwater of such complex coastal systems. This study focuses on a typical Mediterranean coastal agricultural area, located in Tunisia, experiencing substantial NO3 contamination from multiple anthropogenic sources. Here, multiple isotopic tracers (δ18OH2O, δ2HH2O, δ15NNO3, δ18ONO3, and δ11B) combined with a Bayesian isotope MixSIAR model are used (i) to identify the major NO3 sources and their contributions, and (ii) to describe the potential NO3 transformation processes. The measured NO3 concentrations in groundwater are above the natural baseline threshold, suggesting anthropogenic influence. The measured isotopic composition of NO3 indicates that manure, soil organic matter, and sewage are the potential sources of NO3, while δ11B values constrain the NO3 contamination to manure; a finding that is supported by the results of MixSIAR model revealing that manure-derived NO3 dominates over other likely sources. Nitrate derived from manure in the study area is attributed to organic fertilizers used to promote crop growth, and livestock that deposit manure directly on the ground surface. Evidence for ongoing denitrification in groundwaters of the study area is supported by an enrichment in both 15N and 18O in the remaining NO3, although isotopic mass balances between the measured and the theoretical δ18ONO3 values also suggest the occurrence of nitrification. The simultaneous occurrence of these biogeochemical processes with heterogeneous distribution across the study area reflect the complexity of interactions within the investigated coastal aquifer. The multiple isotopic tracer approach used here can identify the effect of multiple NO3 anthropogenic activities in coastal environments, which is fundamental for sustainable groundwater resources management.

Abstract: Drawing on advances in nonstationary frequency analysis and the science of causation and attribution, this study employs a newly developed nonstationary stochastic paired watershed approach to determine the effect of forest harvesting on snowmelt-generated floods. Moreover, this study furthers the application of stochastic physics to evaluate the environmental controls and drivers of flood response. Physically-based climate and time-varying harvesting data are used as covariates to drive the nonstationary flood frequency distribution parameters to detect, attribute, and quantify the effect of harvesting on floods in the snow-dominated Deadman River (878 km2) and nested Joe Ross Creek (99 km2) watersheds. Harvesting only 21% of the watershed caused a 38% and 84% increase in the mean but no increase in variability around the mean of the frequency distribution in the Deadman River and Joe Ross Creek, respectively. Consequently, the 7-year, 20-year, 50-year, and 100-year flood events became approximately two, four, six, and ten times more frequent in both watersheds. An increase in the mean is posited to occur from an increase in moisture availability following harvest from suppressed snow interception and increased net radiation reaching the snowpack. Variability was not increased because snowmelt synchronization was inhibited by the buffering capacity of abundant lakes, evenly distributed aspects, and widespread spatial distribution of cutblocks in the watersheds, preventing any potential for harvesting to increase the efficiency of runoff delivery to the outlet. Consistent with similar recent studies, the effect of logging on floods is controlled not only by the harvest rate but most importantly the physiographic characteristics of the watershed and the spatial distribution of the cutblocks. Imposed by the probabilistic framework to understanding and predicting the relation between extremes and their environmental controls, commonly used in the general sciences but not forest hydrology, it is the inherent nature of snowmelt-driven flood regimes which cause even modest increases in magnitude, especially in the upper tail of the distribution, to translate into surprisingly large changes in frequency. Contrary to conventional wisdom, harvesting influenced small, medium, and very large flood events, and the sensitivity to harvest increased with increasing flood event size and watershed area.

Abstract: Qanat is a passive (unpumped) horizontal well (or a slant well with a very mild inclined angle) that is capable of extracting water from aquifers by gravity. Many qanats are constructed along the radius of the alluvial fan wedge-shaped aquifers. Analytical modeling of such a qanat-aquifer system provides great benefit for quickly screening different designs of qanats and accessing the performance of qanat discharge in the field. The previous analytical modeling of discharge of qanats, however, did not consider the wedge-shaped aquifers. Thus, the goal of this research is to obtain semi-analytical solutions of discharge variations of qanats in alluvial fan aquifers with nearby pumping wells, subjected to areal recharges due to rainfall. The uniform head boundary is considered inside the qanat (because of its enormous permeability in respect to the background aquifer). The influences of the aquifer lateral boundaries on discharge of qanat and its sensitivity to hydraulic and geometric parameters are explored. The influences of the lateral boundaries on the discharge of qanat due to areal recharge and nearby pumping wells discharge are also explored. The results of this study can be utilized for multiple purposes: 1) to predict the discharge of qanat in an alluvial fan aquifer and explore the influences of the areal recharge and nearby pumping well discharge; 2) to estimate the hydraulic parameters of the alluvial fan aquifer depleted by a qanat; 3) to determine the location of the nearby pumping well to minimize its influences on the discharge of a qanat; 4) to calculate the water budgets of aquifers depleted by qanats and pumping wells and replenished by areal recharge among other applications. This paper is an extension to the work presented by Sedghi and Zhan (2020) (which concerns an infinite unconfined aquifer) for an unconfined alluvial fan aquifer setting.