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.

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 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.

Predicting water quality in unmonitored watersheds using artificial neural networks.

Land use and land cover (LULC) play a central role in fate and transport of water quality (WQ) parameters in watersheds. Developing relationships between LULC and WQ parameters is essential for evaluating the quality of water resources. In this paper, we present an artificial neural network (ANN)-based methodology to predict WQ parameters in watersheds with no prior WQ data. The model relies on LULC percentages, temperature, and stream discharge as inputs. The approach is applied to 18 watersheds in west Georgia, United States, having a LULC gradient and varying in size from 2.96 to 26.59 km2. Out of 18 watersheds, 12 were used for training, 3 for validation, and 3 for testing the ANN model. The WQ parameters tested are total dissolved solids (TDS), total suspended solids (TSS), chlorine (Cl), nitrate (NO3), sulfate (SO4), sodium (Na), potassium (K), total phosphorus (TP), and dissolved organic carbon (DOC). Model performances are evaluated on the basis of a performance rating system whereby performances are categorized as unsatisfactory, satisfactory, good, or very good. Overall, the ANN models developed using the training data performed quite well in the independent test watersheds. Based on the rating system TDS, Cl, NO3, SO4, Na, K, and DOC had a performance of at least "good" in all three test watersheds. The average performance for TSS and TP in the three test watersheds were "good." Overall the model performed better in the pastoral and forested watersheds with an average rating of "very good." The average model performance at the urban watershed was "good." This study showed that if WQ and LULC data are available from multiple watersheds in an area with relatively similar physiographic properties, then one can successfully predict the impact of LULC changes on WQ in any nearby watershed.