Wetland Flowpaths Mediate Nitrogen and Phosphorus Concentrations across the Upper Mississippi River Basin.

Eutrophication, harmful algal blooms, and human health impacts are critical environmental challenges resulting from excess nitrogen and phosphorus in surface waters. Yet we have limited information regarding how wetland characteristics mediate water quality across watershed scales. We developed a large, novel set of spatial variables characterizing hydrological flowpaths from wetlands to streams, that is, "wetland hydrological transport variables," to explore how wetlands statistically explain the variability in total nitrogen (TN) and total phosphorus (TP) concentrations across the Upper Mississippi River Basin (UMRB) in the United States. We found that wetland flowpath variables improved landscape-to-aquatic nutrient multilinear regression models (from R2 = 0.89 to 0.91 for TN; R2 = 0.53 to 0.84 for TP) and provided insights into potential processes governing how wetlands influence watershed-scale TN and TP concentrations. Specifically, flowpath variables describing flow-attenuating environments, for example, subsurface transport compared to overland flowpaths, were related to lower TN and TP concentrations. Frequent hydrological connections from wetlands to streams were also linked to low TP concentrations, which likely suggests a nutrient source limitation in some areas of the UMRB. Consideration of wetland flowpaths could inform management and conservation activities designed to reduce nutrient export to downstream waters.

Mapping global non-floodplain wetlands.

Non-floodplain wetlands - those located outside the floodplains - have emerged as integral components to watershed resilience, contributing hydrologic and biogeochemical functions affecting watershed-scale flooding extent, drought magnitude, and water-quality maintenance. However, the absence of a global dataset of non-floodplain wetlands limits their necessary incorporation into water quality and quantity management decisions and affects wetland-focused wildlife habitat conservation outcomes. We addressed this critical need by developing a publicly available "Global NFW" (Non-Floodplain Wetland) dataset, comprised of a global river-floodplain map at 90 m resolution coupled with a global ensemble wetland map incorporating multiple wetland-focused data layers. The floodplain, wetland, and non-floodplain wetland spatial data developed here were successfully validated within 21 large and heterogenous basins across the conterminous United States. We identified nearly 33 million potential non-floodplain wetlands with an estimated global extent of over 16×106 km2. Non-floodplain wetland pixels comprised 53% of globally identified wetland pixels, meaning the majority of the globe's wetlands likely occur external to river floodplains and coastal habitats. The identified global NFWs were typically small (median 0.039 km2), with a global median size ranging from 0.018-0.138 km2. This novel geospatial Global NFW static dataset advances wetland conservation and resource-management goals while providing a foundation for global non-floodplain wetland functional assessments, facilitating non-floodplain wetland inclusion in hydrological, biogeochemical, and biological model development. The data are freely available through the United States Environmental Protection Agency's Environmental Dataset Gateway (https://gaftp.epa.gov/EPADataCommons/ORD/Global_NonFloodplain_Wetlands/, last access: 24 May 2023) and through https://doi.org/10.23719/1528331 (Lane et al., 2023a).

River Basin Simulations Reveal Wide-Ranging Wetland-Mediated Nitrate Reductions.

River basin-scale wetland restoration and creation is a primary management option for mitigating nitrogen-based water quality challenges. However, the magnitude of nitrogen reduction that will result from adding wetlands across large river basins is uncertain, partly because the areal extent, location, and physical and functional characteristics of the wetlands are unknown. We simulated over 3600 wetland restoration scenarios across the ∼450,000 km2 Upper Mississippi River Basin (UMRB) depicting varied assumptions for wetland areal extent, physical and functional characteristics, and placement strategy. These simulations indicated that restoring wetlands will reduce local nitrate yields and nitrate loads at the UMRB outlet. However, the projected magnitude of nitrate reduction varied widely across disparate scenario assumptions─e.g., restoring 4500 km2 of wetlands (i.e., 1% of UMRB area) decreased mean annual nitrate loads at the UMRB outlet between 3 and 42%. Higher magnitude nitrate reductions correlated with best-case assumptions, particularly for characteristics controlling nitrate loading rates to the wetlands. These results show that simplified claims about basin-scale wetland-mediated water quality improvements discount the breadth of possible wetland impacts across disparate wetland physical and functional conditions and highlight a need for greater clarity regarding the likelihood of these conditions at river basin scales.

Vulnerable Waters are Essential to Watershed Resilience.

Watershed resilience is the ability of a watershed to maintain its characteristic system state while concurrently resisting, adapting to, and reorganizing after hydrological (for example, drought, flooding) or biogeochemical (for example, excessive nutrient) disturbances. Vulnerable waters include non-floodplain wetlands and headwater streams, abundant watershed components representing the most distal extent of the freshwater aquatic network. Vulnerable waters are hydrologically dynamic and biogeochemically reactive aquatic systems, storing, processing, and releasing water and entrained (that is, dissolved and particulate) materials along expanding and contracting aquatic networks. The hydrological and biogeochemical functions emerging from these processes affect the magnitude, frequency, timing, duration, storage, and rate of change of material and energy fluxes among watershed components and to downstream waters, thereby maintaining watershed states and imparting watershed resilience. We present here a conceptual framework for understanding how vulnerable waters confer watershed resilience. We demonstrate how individual and cumulative vulnerable-water modifications (for example, reduced extent, altered connectivity) affect watershed-scale hydrological and biogeochemical disturbance response and recovery, which decreases watershed resilience and can trigger transitions across thresholds to alternative watershed states (for example, states conducive to increased flood frequency or nutrient concentrations). We subsequently describe how resilient watersheds require spatial heterogeneity and temporal variability in hydrological and biogeochemical interactions between terrestrial systems and down-gradient waters, which necessitates attention to the conservation and restoration of vulnerable waters and their downstream connectivity gradients. To conclude, we provide actionable principles for resilient watersheds and articulate research needs to further watershed resilience science and vulnerable-water management.

Wetland restoration yields dynamic nitrate responses across the Upper Mississippi river basin.

Wetland restoration is a primary management option for removing surplus nitrogen draining from agricultural landscapes. However, wetland capacity to mitigate nitrogen losses at large river-basin scales remains uncertain. This is largely due to a limited number of studies that address the cumulative and dynamic effects of restored wetlands across the landscape on downstream nutrient conditions. We analyzed wetland restoration impacts on modeled nitrate dynamics across 279 subbasins comprising the ∼0.5 million km2 Upper Mississippi River Basin (UMRB), USA, which covers eight states and houses ∼30 million people. Restoring ∼8,000 km2 of wetlands will reduce mean annual nitrate loads to the UMRB outlet by 12%, a substantial improvement over existing conditions but markedly less than widely cited estimates. Our lower wetland efficacy estimates are partly attributed to improved representation of processes not considered by preceding empirical studies - namely the potential for nitrate to bypass wetlands (i.e., via subsurface tile drainage) and be stored or transformed within the river network itself. Our novel findings reveal that wetlands mitigate surplus nitrogen basin-wide, yet they may not be as universally effective in tiled landscapes and because of river network processing.

Variation in stream network relationships and geospatial predictions of watershed conductivity.

Secondary salinization, the increase of anthropogenically-derived salts in freshwaters, threatens freshwater biota and ecosystems, drinking water supplies, and infrastructure. The various anthropogenic sources of salts and their locations in a watershed may result in secondary salinization of river and stream networks through multiple inputs. We developed a watershed predictive assessment to investigate the degree to which topology, land-cover, and land-use covariates affect stream specific conductivity (SC), a measure of salinity. We used spatial stream network models to predict SC throughout an Appalachian stream network in a watershed affected by surface coal mining. During high-discharge conditions, 8 to 44% of stream km in the watershed exceeded the SC benchmark of 300 µS/cm, which is meant to be protective of aquatic life in the Central Appalachian ecoregion. During low-discharge conditions, 96 to 100% of stream km exceeded the benchmark. The 2 different discharge conditions altered the spatial dependency of SC among the stream monitoring sites. During most low discharges, SC was a function of upstream-to-downstream network distances, or flow-connected distances, among the sites. Flow-connected distances are indicative of upstream dependencies affecting stream SC. During high discharge, SC was related to both flow-connected distances and flow-unconnected distances (i.e., distances between sites on different branches of the network). Flow-unconnected distances are indicative of processes on adjacent branches and their catchments affecting stream SC. With sites distributed from headwaters to the watershed outlet, the extent of impacts from secondary salinization could be better spatially predicted and assessed with spatial stream network models than with models assuming spatial independence. Importantly, the assessment also recognized the multi-scale spatial relationships that can occur between the landscape and stream network.

Spatial Convergence in Major Dissolved Ion Concentrations and Implications of Headwater Mining for Downstream Water Quality.

Spatial patterns in major dissolved solute concentrations were examined to better understand impact of surface coal mining in headwaters on downstream water chemistry. Sixty sites were sampled seasonally from 2012 to 2014 in an eastern Kentucky watershed. Watershed areas (WA) ranged from 1.6 to 400.5 km2 and were mostly forested (58%-95%), but some drained as much as 31% surface mining. Measures of total dissolved solutes and most component ions were positively correlated with mining. Analytes showed strong convergent spatial patterns with high variability in headwaters (<15 km2 WA) that stabilized downstream (WA > 75 km2), indicating hydrologic mixing primarily controls downstream values. Mean headwater solute concentrations were a good predictor of downstream values, with % differences ranging from 0.55% (Na+) to 28.78% (Mg2+). In a mined scenario where all headwaters had impacts, downstream solute concentrations roughly doubled. Alternatively, if mining impacts to headwaters were minimized, downstream solute concentrations better approximated the 300 µS/cm conductivity criterion deemed protective of aquatic life. Temporal variability also had convergent spatial patterns and mined streams were less variable due to unnaturally stable hydrology. The highly conserved nature of dissolved solutes from mining activities and lack of viable treatment options suggest forested, unmined watersheds would provide dilution that would be protective of downstream aquatic life.

Depressional wetlands affect watershed hydrological, biogeochemical, and ecological functions.

Depressional wetlands of the extensive U.S. and Canadian Prairie Pothole Region afford numerous ecosystem processes that maintain healthy watershed functioning. However, these wetlands have been lost at a prodigious rate over past decades due to drainage for development, climate effects, and other causes. Options for management entities to protect the existing wetlands, and their functions, may focus on conserving wetlands based on spatial location vis-à-vis a floodplain or on size limitations (e.g., permitting smaller wetlands to be destroyed but not larger wetlands). Yet the effects of such management practices and the concomitant loss of depressional wetlands on watershed-scale hydrological, biogeochemical, and ecological functions are largely unknown. Using a hydrological model, we analyzed how different loss scenarios by wetland size and proximal location to the stream network affected watershed storage (i.e., inundation patterns and residence times), connectivity (i.e., streamflow contributing areas), and export (i.e., streamflow) in a large watershed in the Prairie Pothole Region of North Dakota, USA. Depressional wetlands store consequential amounts of precipitation and snowmelt. The loss of smaller depressional wetlands (<3.0 ha) substantially decreased landscape-scale inundation heterogeneity, total inundated area, and hydrological residence times. Larger wetlands act as hydrologic "gatekeepers," preventing surface runoff from reaching the stream network, and their modeled loss had a greater effect on streamflow due to changes in watershed connectivity and storage characteristics of larger wetlands. The wetland management scenario based on stream proximity (i.e., protecting wetlands 30 m and ~450 m from the stream) alone resulted in considerable landscape heterogeneity loss and decreased inundated area and residence times. With more snowmelt and precipitation available for runoff with wetland losses, contributing area increased across all loss scenarios. We additionally found that depressional wetlands attenuated peak flows; the probability of increased downstream flooding from wetland loss was also consistent across all loss scenarios. It is evident from this study that optimizing wetland management for one end goal (e.g., protection of large depressional wetlands for flood attenuation) over another (e.g., protecting of small depressional wetlands for biodiversity) may come at a cost for overall watershed hydrological, biogeochemical, and ecological resilience, functioning, and integrity.

Enhancing protection for vulnerable waters.

Governments worldwide do not adequately protect their limited freshwater systems and therefore place freshwater functions and attendant ecosystem services at risk. The best available scientific evidence compels enhanced protections for freshwater systems, especially for impermanent streams and wetlands outside of floodplains that are particularly vulnerable to alteration or destruction. New approaches to freshwater sustainability - implemented through scientifically informed adaptive management - are required to protect freshwater systems through periods of changing societal needs. One such approach introduced in the US in 2015 is the Clean Water Rule, which clarified the jurisdictional scope for federally protected waters. However, within hours of its implementation litigants convinced the US Court of Appeals for the Sixth Circuit to stay the rule, and the subsequently elected administration has now placed it under review for potential revision or rescission. Regardless of its outcome at the federal level, policy and management discussions initiated by the propagation of this rare rulemaking event have potential far-reaching implications at all levels of government across the US and worldwide. At this timely juncture, we provide a scientific rationale and three policy options for all levels of government to meaningfully enhance protection of these vulnerable waters. A fourth option, a 'do-nothing' approach, is wholly inconsistent with the well-established scientific evidence of the importance of these vulnerable waters.

Do geographically isolated wetlands influence landscape functions?

Geographically isolated wetlands (GIWs), those surrounded by uplands, exchange materials, energy, and organisms with other elements in hydrological and habitat networks, contributing to landscape functions, such as flow generation, nutrient and sediment retention, and biodiversity support. GIWs constitute most of the wetlands in many North American landscapes, provide a disproportionately large fraction of wetland edges where many functions are enhanced, and form complexes with other water bodies to create spatial and temporal heterogeneity in the timing, flow paths, and magnitude of network connectivity. These attributes signal a critical role for GIWs in sustaining a portfolio of landscape functions, but legal protections remain weak despite preferential loss from many landscapes. GIWs lack persistent surface water connections, but this condition does not imply the absence of hydrological, biogeochemical, and biological exchanges with nearby and downstream waters. Although hydrological and biogeochemical connectivity is often episodic or slow (e.g., via groundwater), hydrologic continuity and limited evaporative solute enrichment suggest both flow generation and solute and sediment retention. Similarly, whereas biological connectivity usually requires overland dispersal, numerous organisms, including many rare or threatened species, use both GIWs and downstream waters at different times or life stages, suggesting that GIWs are critical elements of landscape habitat mosaics. Indeed, weaker hydrologic connectivity with downstream waters and constrained biological connectivity with other landscape elements are precisely what enhances some GIW functions and enables others. Based on analysis of wetland geography and synthesis of wetland functions, we argue that sustaining landscape functions requires conserving the entire continuum of wetland connectivity, including GIWs.

Automated riverine landscape characterization: GIS-based tools for watershed-scale research, assessment, and management.

River systems consist of hydrogeomorphic patches (HPs) that emerge at multiple spatiotemporal scales. Functional process zones (FPZs) are HPs that exist at the river valley scale and are important strata for framing whole-watershed research questions and management plans. Hierarchical classification procedures aid in HP identification by grouping sections of river based on their hydrogeomorphic character; however, collecting data required for such procedures with field-based methods is often impractical. We developed a set of GIS-based tools that facilitate rapid, low cost riverine landscape characterization and FPZ classification. Our tools, termed RESonate, consist of a custom toolbox designed for ESRI ArcGIS®. RESonate automatically extracts 13 hydrogeomorphic variables from readily available geospatial datasets and datasets derived from modeling procedures. An advanced 2D flood model, FLDPLN, designed for MATLAB® is used to determine valley morphology by systematically flooding river networks. When used in conjunction with other modeling procedures, RESonate and FLDPLN can assess the character of large river networks quickly and at very low costs. Here we describe tool and model functions in addition to their benefits, limitations, and applications.

Temporal variation in spatial sources of discharge in a large watershed.

We examined how the spatial configuration of source areas for runoff varied over time in a large watershed, in order to understand processes governing material loading to rivers. Discharge source areas within the Fox River watershed (Wisconsin, US) were mapped for two individual discharge events. The spatial distribution of source areas varied between and over the duration of individual discharge events. Relative contribution to runoff by land cover types within source areas was quantified and compared to areal abundance of land covers in the watershed. Contributions of runoff by land cover types varied over time. Moreover, the degree to which different land cover types acted as source areas differed from their abundance in the watershed. Hence, areal quantifications of land cover within a watershed may not accurately represent what land covers are source areas over given time periods. Therefore, a source-area-based approach may yield more accurate spatial analysis of material loading patterns than a watershed-based approach.