Nutrient Removal Potential of Headwater Wetlands in Coastal Plains of Alabama, USA.

Headwater streams drain over 70% of the land in the United States with headwater wetlands covering 6.59 million hectares. These ecosystems are important landscape features in the southeast United States, with underlying effects on ecosystem health, water yield, nutrient cycling, biodiversity, and water quality. However, little is known about the relationship between headwater wetlands' nutrient function (i.e., nutrient load removal (RL) and removal efficiency (ER)) and their physical characteristics. Here, we investigate this relationship for 44 headwater wetlands located within the Upper Fish River watershed (UFRW) in coastal Alabama. To accomplish this objective, we apply the process-based watershed model SWAT (Soil and Water Assessment Tool) to generate flow and nutrient loadings to each study wetland and subsequently quantify the wetland-level nutrient removal efficiencies using the process-based wetland model WetQual. Results show that the calculated removal efficiencies of the headwater wetlands in the UFRW are 75-84% and 27-35% for nitrate (NO3-) and phosphate (PO4+), respectively. The calculated nutrient load removals are highly correlated with the input loads, and the estimated PO4+ ER shows a significant decreasing trend with increased input loadings. The relationship between NO3-ER and wetland physical characteristics such as area, volume, and residence time is statistically insignificant (p > 0.05), while for PO4+, the correlation is positive and statistically significant (p < 0.05). On the other hand, flashiness (flow pulsing) and baseflow index (fraction of inflow that is coming from baseflow) have a strong effect on NO3- removal but not on PO4+ removal. Modeling results and statistical analysis point toward denitrification and plant uptake as major NO3- removal mechanisms, whereas plant uptake, diffusion, and settling of sediment-bound P were the main mechanisms for PO4+ removal. Additionally, the computed nutrient ER is higher during the driest year of the simulated period compared to during the wettest year. Our findings are in line with global-level studies and offer new insights into wetland physical characteristics affecting nutrient removal efficiency and the importance of headwater wetlands in mitigating water quality deterioration in coastal areas. The regression relationships for NO3- and PO4+ load removals in the selected 44 wetlands are then used to extrapolate nutrient load removals to 348 unmodeled non-riverine and non-riparian wetlands in the UFRW (41% of UFRW drains to them). Results show that these wetlands remove 51-61% of the NO3- and 5-10% of the PO4+ loading they receive from their respective drainage areas. Due to geographical proximity and physiographic similarity, these results can be scaled up to the coastal plains of Alabama and Northwest Florida.

Sub-monthly time scale forecasting of harmful algal blooms intensity in Lake Erie using remote sensing and machine learning.

Harmful algal blooms of cyanobacteria (CyanoHAB) have emerged as a serious environmental concern in large and small water bodies including many inland lakes. The growth dynamics of CyanoHAB can be chaotic at very short timescales but predictable at coarser timescales. In Lake Erie, cyanobacteria blooms occur in the spring-summer months, which, at annual timescale, are controlled by the total spring phosphorus (TP) load into the lake. This study aimed to forecast CyanoHAB cell count at sub-monthly (e.g., 10-day) timescales. Satellite-derived cyanobacterial index (CI) was used as a surrogate measure of CyanoHAB cell count. CI was related to the in-situ measured chlorophyll-a and phycocyanin concentrations and Microcystis biovolume in the lake. Using available data on environmental and lake hydrodynamics as predictor variables, four statistical models including LASSO (Least Absolute Shrinkage and Selection Operator), artificial neural network (ANN), random forest (RF), and an ensemble average of the three models (EA) were developed to forecast CI at 10-, 20- and 30-day lead times. The best predictions were obtained by using the RF and EA algorithms. It was found that CyanoHAB growth dynamics, even at sub-monthly timescales, are determined by coarser timescale variables. Meteorological, hydrological, and water quality variations at sub-monthly timescales exert lesser control over CyanoHAB growth dynamics. Nutrients discharged into the lake from rivers other than the Maumee River were also important in explaining the variations in CI. Surprisingly, to forecast CyanoHAB cell count, average solar radiation at 30 to 60 days lags were found to be more important than the average solar radiation at 0 to 30 days lag. Other important variables were TP discharged into the lake during the previous 10 years, TP and TKN discharged into the lake during the previous 120 days, the average water level at 10-day lag and 60-day lag.

A Machine Learning Approach to Predict Watershed Health Indices for Sediments and Nutrients at Ungauged Basins.

Effective water quality management and reliable environmental modeling depend on the availability, size, and quality of water quality (WQ) data. Observed stream water quality data are usuallEEy sparse in both time and space. Reconstruction of water quality time series using surrogate variables such as streamflow have been used to evaluate risk metrics such as reliability, resilience, vulnerability, and watershed health (WH) but only at gauged locations. Estimating these indices for ungauged watersheds has not been attempted because of the high-dimensional nature of the potential predictor space. In this study, machine learning (ML) models, namely random forest regression, AdaBoost, gradient boosting machines, and Bayesian ridge regression (along with an ensemble model), were evaluated to predict watershed health and other risk metrics at ungauged hydrologic unit code 10 (HUC-10) basins using watershed attributes, long-term climate data, soil data, land use and land cover data, fertilizer sales data, and geographic information as predictor variables. These ML models were tested over the Upper Mississippi River Basin, the Ohio River Basin, and the Maumee River Basin for water quality constituents such as suspended sediment concentration, nitrogen, and phosphorus. Random forest, AdaBoost, and gradient boosting regressors typically showed a coefficient of determination R2>0.8 for suspended sediment concentration and nitrogen during the testing stage, while the ensemble model exhibited R2>0.95. Watershed health values with respect to suspended sediments and nitrogen predicted by all ML models including the ensemble model were lower for areas with larger agricultural land use, moderate for areas with predominant urban land use, and higher for forested areas; the trained ML models adequately predicted WH in ungauged basins. However, low WH values (with respect to phosphorus) were predicted at some basins in the Upper Mississippi River Basin that had dominant forest land use. Results suggest that the proposed ML models provide robust estimates at ungauged locations when sufficient training data are available for a WQ constituent. ML models may be used as quick screening tools by decision makers and water quality monitoring agencies for identifying critical source areas or hotspots with respect to different water quality constituents, even for ungauged watersheds.

Two-Layer numerical model of soil moisture dynamics: Model assessment and Bayesian uncertainty estimation.

A two-layer model based on the integrated form of Richards' equation (RE) was recently developed to simulate the soil water movement in the roots layer and the vadose zone with a relatively shallow and dynamic water table. The model simulates thickness-averaged volumetric water content and matric suction as opposed to point values and was numerically verified for three soil textures using HYDRUS as a benchmark. However, the strengths and limitations of the two-layer model and its performance in stratified soils and under actual field conditions have not been tested. This study further examined the two-layer model using two numerical verification experiments and, most importantly, tested its performance at site-level under actual, highly variable hydroclimate conditions. Moreover, model parameters were estimated and uncertainty and sources of errors were quantified using a Bayesian framework. First, the two-layer model was evaluated for 231 soil textures under varying soil layer thicknesses with a uniform soil profile. Second, the two-layer model was assessed for stratified conditions where the top and bottom soil layers have contrasting hydraulic conductivities. The model was evaluated by comparing soil moisture and flux estimates to those from the HYDRUS model. Last, a case study of model application using data from a Soil Climate Analysis Network (SCAN) site was presented. Bayesian Monte Carlo (BMC) method was implemented for model calibration and quantifying sources of uncertainty under real hydroclimate and soil conditions. For a homogeneous soil profile, the two-layer model generally had excellent performance in estimating volumetric water content and fluxes, while the model performance slightly declined with increasing layer thickness and coarser textured soils. The model configurations regarding layer thicknesses and soil textures that generate accurate soil moisture and flux estimations were further suggested. With the two layers of contrasting permeability, model-simulated soil moisture contents and fluxes agreed well with those computed by HYDRUS, indicating that the two-layer model accurately handles the water flow dynamics around the layer interface. In the field application, given the highly variable hydroclimate conditions, the two-layer model combined with the BMC method showed good agreement with the observed average soil moisture of the root zone and the vadose zone below (RMSE <0.021 during calibration and <0.023 during validation periods). The contribution of parametric uncertainty to the total model uncertainty was too small compared to other sources. The numerical tests and the site level application showed that the two-layer model can reliably simulate thickness-averaged soil moisture and estimate fluxes in the vadose zone under various soil and hydroclimate conditions. Results also indicated that the BMC method could be a robust framework for vadose zone hydraulic parameters identification and model uncertainty estimation.

A secondary assessment of sediment trapping effectiveness by vegetated buffers.

Vegetated buffers and filter strips are a widely used Best Management Practice (BMP) for enhancing streamside ecosystem quality and water quality improvement through nonpoint source pollutant removal. Most existing studies are either site-specific, rely on limited data points, or evaluate buffer width and slope as the only design variables for predicting sediment reduction, not considering other parameters such as soil texture, vegetation types, and runoff loads that can significantly influence the buffer efficiency. In this paper, we carry out a meta-analysis of published studies and fit regression models to explore the sediment removal capacity of riparian buffers. We compiled 905 data points from over 90 studies (including data from an online BMP database) documenting sediment trapping by vegetated buffers and recorded data regarding buffer characteristics such as buffer width, slope, area, vegetation type, sediment loading, water flow rates, and sediment removal efficiency. We found that an exponential regression model describing the relationship between sediment removal efficiency by the buffer and water inflow/outflow volume ratio explained 44% of the variance. Adding the square root of roughness increased the R 2 to 0.50. The model performance was compared with other sediment reduction regression models reported in the literature. The results point towards the importance of considering flow parameters in vegetative buffer design. The improved empirical relationships derived here can be used at local scales to understand sediment trapping potential by vegetated buffers for water quality mitigation purposes and can be built into extant hydrologic models for improved watershed-scale assessments.

A two-layer numerical model of soil moisture dynamics: Model development.

Simulating water moisture flow in variably saturated soils with a relatively shallow water table is challenging due to the high nonlinear behavior of Richards' equation (RE). A two-layer approximation of RE was derived in this paper, which describes vertically-averaged soil moisture content and flow dynamics in the root zone and the unsaturated soil below. To this end, the partial differential equation (PDE) describing RE was converted into two-coupled ordinary differential equations (ODEs) describing dynamic vertically-averaged soil moisture variations in the two soil zones subject to a deep or shallow water table in addition to variable soil moisture flux and pressure conditions at the surface. The coupled ODEs were solved numerically using the iterative Huen's method for a variety of flux and pressure-controlled top and bottom boundary conditions (BCs). The numerical model was evaluated for three typical soil textures with free-drainage and mixed flux-pressure head at the bottom boundary under various atmospheric conditions. The results of soil water contents and fluxes were validated using HYDRUS-1D as a benchmark. Simulated values showed that the new model is numerically stable and generally accurate in simulating vertically-averaged soil moisture in the two layers under various flux and prescribed pressure BCs. A hypothetical simulation scenario involving desaturation of initially saturated soil profile caused by exponentially declining water table demonstrated the robustness of the numerical model in tracking vertically-averaged moisture contents in the roots layer and the lower vadose soil as the water table continued to fall. The two-layer model can be used by researchers to simulate variably saturated soils in wetlands and by water resources planners for efficient coupling of land-surface systems to groundwater and management of conjunctive use of surface and groundwater.

Uncertainty quantification in reconstruction of sparse water quality time series: Implications for watershed health and risk-based TMDL assessment.

Despite the plethora of methods available for uncertainty quantification, their use has been limited in the practice of water quality (WQ) modeling. In this paper, a decision support tool (DST) that yields a continuous time series of WQ loads from sparse data using streamflows as predictor variables is presented. The DST estimates uncertainty by analyzing residual errors using a relevance vector machine. To highlight the importance of uncertainty quantification, two applications enabled within the DST are discussed. The DST computes (i) probability distributions of four measures of WQ risk analysis- reliability, resilience, vulnerability, and watershed health- as opposed to single deterministic values and (ii) concentration/load reduction required in a WQ constituent to meet total maximum daily load (TMDL) targets along with the associated risk of failure. Accounting for uncertainty reveals that a deterministic analysis may mislead about the WQ risk and the level of compliance attained with established TMDLs.

Composite measures of watershed health from a water quality perspective.

Water quality data at gaging stations are typically compared with established federal, state, or local water quality standards to determine if violations (concentrations of specific constituents falling outside acceptable limits) have occurred. Based on the frequency and severity of water quality violations, risk metrics such as reliability, resilience, and vulnerability (R-R-V) are computed for assessing water quality-based watershed health. In this study, a modified methodology for computing R-R-V measures is presented, and a new composite watershed health index is proposed. Risk-based assessments for different water quality parameters are carried out using identified national sampling stations within the Upper Mississippi River Basin, the Maumee River Basin, and the Ohio River Basin. The distributional properties of risk measures with respect to water quality parameters are reported. Scaling behaviors of risk measures using stream order, specifically for the watershed health (WH) index, suggest that WH values increased with stream order for suspended sediment concentration, nitrogen, and orthophosphate in the Upper Mississippi River Basin. Spatial distribution of risk measures enable identification of locations exhibiting poor watershed health with respect to the chosen numerical standard, and the role of land use characteristics within the watershed.

Bayesian Monte Carlo and maximum likelihood approach for uncertainty estimation and risk management: Application to lake oxygen recovery model.

Model uncertainty estimation and risk assessment is essential to environmental management and informed decision making on pollution mitigation strategies. In this study, we apply a probabilistic methodology, which combines Bayesian Monte Carlo simulation and Maximum Likelihood estimation (BMCML) to calibrate a lake oxygen recovery model. We first derive an analytical solution of the differential equation governing lake-averaged oxygen dynamics as a function of time-variable wind speed. Statistical inferences on model parameters and predictive uncertainty are then drawn by Bayesian conditioning of the analytical solution on observed daily wind speed and oxygen concentration data obtained from an earlier study during two recovery periods on a eutrophic lake in upper state New York. The model is calibrated using oxygen recovery data for one year and statistical inferences were validated using recovery data for another year. Compared with essentially two-step, regression and optimization approach, the BMCML results are more comprehensive and performed relatively better in predicting the observed temporal dissolved oxygen levels (DO) in the lake. BMCML also produced comparable calibration and validation results with those obtained using popular Markov Chain Monte Carlo technique (MCMC) and is computationally simpler and easier to implement than the MCMC. Next, using the calibrated model, we derive an optimal relationship between liquid film-transfer coefficient for oxygen and wind speed and associated 95% confidence band, which are shown to be consistent with reported measured values at five different lakes. Finally, we illustrate the robustness of the BMCML to solve risk-based water quality management problems, showing that neglecting cross-correlations between parameters could lead to improper required BOD load reduction to achieve the compliance criteria of 5 mg/L.

Aggregate Measures of Watershed Health from Reconstructed Water Quality Data with Uncertainty.

Risk-based measures such as reliability, resilience, and vulnerability (R-R-V) have the potential to serve as watershed health assessment tools. Recent research has demonstrated the applicability of such indices for water quality (WQ) constituents such as total suspended solids and nutrients on an individual basis. However, the calculations can become tedious when time-series data for several WQ constituents have to be evaluated individually. Also, comparisons between locations with different sets of constituent data can prove difficult. In this study, data reconstruction using a relevance vector machine algorithm was combined with dimensionality reduction via variational Bayesian noisy principal component analysis to reconstruct and condense sparse multidimensional WQ data sets into a single time series. The methodology allows incorporation of uncertainty in both the reconstruction and dimensionality-reduction steps. The R-R-V values were calculated using the aggregate time series at multiple locations within two Indiana watersheds. Results showed that uncertainty present in the reconstructed WQ data set propagates to the aggregate time series and subsequently to the aggregate R-R-V values as well. This data-driven approach to calculating aggregate R-R-V values was found to be useful for providing a composite picture of watershed health. Aggregate R-R-V values also enabled comparison between locations with different types of WQ data.

Applying a statewide geospatial leaching tool for assessing soil vulnerability ratings for agrochemicals across the contiguous United States.

A large-scale leaching assessment tool not only illustrates soil (or groundwater) vulnerability in unmonitored areas, but also can identify areas of potential concern for agrochemical contamination. This study describes the methodology of how the statewide leaching tool in Hawaii modified recently for use with pesticides and volatile organic compounds can be extended to the national assessment of soil vulnerability ratings. For this study, the tool was updated by extending the soil and recharge maps to cover the lower 48 states in the United States (US). In addition, digital maps of annual pesticide use (at a national scale) as well as detailed soil properties and monthly recharge rates (at high spatial and temporal resolutions) were used to examine variations in the leaching (loads) of pesticides for the upper soil horizons. Results showed that the extended tool successfully delineated areas of high to low vulnerability to selected pesticides. The leaching potential was high for picloram, medium for simazine, and low to negligible for 2,4-D and glyphosate. The mass loadings of picloram moving below 0.5 m depth increased greatly in northwestern and central US that recorded its extensive use in agricultural crops. However, in addition to the amount of pesticide used, annual leaching load of atrazine was also affected by other factors that determined the intrinsic aquifer vulnerability such as soil and recharge properties. Spatial and temporal resolutions of digital maps had a great effect on the leaching potential of pesticides, requiring a trade-off between data availability and accuracy. Potential applications of this tool include the rapid, large-scale vulnerability assessments for emerging contaminants which are hard to quantify directly through vadose zone models due to lack of full environmental data.

Watershed reliability, resilience and vulnerability analysis under uncertainty using water quality data.

A method for assessment of watershed health is developed by employing measures of reliability, resilience and vulnerability (R-R-V) using stream water quality data. Observed water quality data are usually sparse, so that a water quality time-series is often reconstructed using surrogate variables (streamflow). A Bayesian algorithm based on relevance vector machine (RVM) was employed to quantify the error in the reconstructed series, and a probabilistic assessment of watershed status was conducted based on established thresholds for various constituents. As an application example, observed water quality data for several constituents at different monitoring points within the Cedar Creek watershed in north-east Indiana (USA) were utilized. Considering uncertainty in the data for the period 2002-2007, the R-R-V analysis revealed that the Cedar Creek watershed tends to be in compliance with respect to selected pesticides, ammonia and total phosphorus. However, the watershed was found to be prone to violations of sediment standards. Ignoring uncertainty in the water quality time-series led to misleading results especially in the case of sediments. Results indicate that the methods presented in this study may be used for assessing the effects of different stressors over a watershed. The method shows promise as a management tool for assessing watershed health.

An auxiliary method to reduce potential adverse impacts of projected land developments: subwatershed prioritization.

An index based method is developed that ranks the subwatersheds of a watershed based on their relative impacts on watershed response to anticipated land developments, and then applied to an urbanizing watershed in Eastern Pennsylvania. Simulations with a semi-distributed hydrologic model show that computed low- and high-flow frequencies at the main outlet increase significantly with the projected landscape changes in the watershed. The developed index is utilized to prioritize areas in the urbanizing watershed based on their contributions to alterations in the magnitude of selected flow characteristics at two spatial resolutions. The low-flow measure, 7Q10, rankings are shown to mimic the spatial trend of groundwater recharge rates, whereas average annual maximum daily flow, QAMAX, and average monthly median of daily flows, QMMED, rankings are influenced by both recharge and proximity to watershed outlet. Results indicate that, especially with the higher resolution, areas having quicker responses are not necessarily the more critical areas for high-flow scenarios. Subwatershed rankings are shown to vary slightly with the location of water quality/quantity criteria enforcement. It is also found that rankings of subwatersheds upstream from the site of interest, which could be the main outlet or any interior point in the watershed, may be influenced by the time scale of the hydrologic processes.