Climate Change Will Impact Surface Water Extents and Dynamics Across the Central United States.

Climate change is projected to impact river, lake, and wetland hydrology, with global implications for the condition and productivity of aquatic ecosystems. We integrated Sentinel-1 and Sentinel-2 based algorithms to track monthly surface water extent (2017-2021) for 32 sites across the central United States (U.S.). Median surface water extent was highly variable across sites, ranging from 3.9% to 45.1% of a site. To account for landscape-based differences (e.g., water storage capacity, land use) in the response of surface water extents to meteorological conditions, individual statistical models were developed for each site. Future changes to climate were defined as the difference between 2006-2025 and 2061-2080 using MACA-CMIP5 (MACAv2-METDATA) Global Circulation Models. Time series of climate change adjusted surface water extents were projected. Annually, 19 of the 32 sites under RCP4.5 and 22 of the 32 sites under RCP8.5 were projected to show an average decline in surface water extent, with drying most consistent across the southeast central, southwest central, and midwest central U.S. Projected declines under surface water dry conditions at these sites suggest greater impacts of drought events are likely in the future. Projected changes were seasonally variable, with the greatest decline in surface water extent expected in summer and fall seasons. In contrast, many north central sites showed a projected increase in surface water in most seasons, relative to the 2017-2021 period, likely attributable to projected increases in winter and spring precipitation exceeding increases in projected temperature.

High-frequency time series comparison of Sentinel-1 and Sentinel-2 satellites for mapping open and vegetated water across the United States (2017-2021).

Frequent observations of surface water at fine spatial scales will provide critical data to support the management of aquatic habitat, flood risk and water quality. Sentinel-1 and Sentinel-2 satellites can provide such observations, but algorithms are still needed that perform well across diverse climate and vegetation conditions. We developed surface inundation algorithms for Sentinel-1 and Sentinel-2, respectively, at 12 sites across the conterminous United States (CONUS), covering a total of >536,000 km2 and representing diverse hydrologic and vegetation landscapes. Each scene in the 5-year (2017-2021) time series was classified into open water, vegetated water, and non-water at 20 m resolution using variables from Sentinel-1 and Sentinel-2, as well as variables derived from topographic and weather datasets. The Sentinel-1 algorithm was developed distinct from the Sentinel-2 model to explore if and where the two time series could potentially be integrated into a single high-frequency time series. Within each model, open water and vegetated water (vegetated palustrine, lacustrine, and riverine wetlands) classes were mapped. The models were validated using imagery from WorldView and PlanetScope. Classification accuracy for open water was high across the 5-year period, with an omission and commission error of only 3.1% and 0.9% for the Sentinel-1 algorithm and 3.1% and 0.5% for the Sentinel-2 algorithm, respectively. Vegetated water accuracy was lower, as expected given that the class represents mixed pixels. The Sentinel-2 algorithm showed higher accuracy (10.7% omission and 7.9% commission error) relative to the Sentinel-1 algorithm (28.4% omission and 16.0% commission error). Patterns over time in the proportion of area mapped as open or vegetated water by the Sentinel-1 and Sentinel-2 algorithms were charted and correlated for a subset of all 12 sites. Our results showed that the Sentinel-1 and Sentinel-2 algorithm open water time series can be integrated at all 12 sites to improve the temporal resolution, but sensor-specific differences, such as sensitivity to vegetation structure versus pixel color, complicate the data integration for mixed-pixel, vegetated water. The methods developed here provide inundation at 5-day (Sentinel-2 algorithm) and 12-day (Sentinel-1 algorithm) time steps to improve our understanding of the short- and long-term response of surface water to climate and land use drivers in different ecoregions.

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.

Tracking rates of postfire conifer regeneration vs. deciduous vegetation recovery across the western United States.

Postfire shifts in vegetation composition will have broad ecological impacts. However, information characterizing postfire recovery patterns and their drivers are lacking over large spatial extents. In this analysis, we used Landsat imagery collected when snow cover (SCS) was present, in combination with growing season (GS) imagery, to distinguish evergreen vegetation from deciduous vegetation. We sought to (1) characterize patterns in the rate of postfire, dual-season Normalized Difference Vegetation Index (NDVI) across the region, (2) relate remotely sensed patterns to field-measured patterns of re-vegetation, and (3) identify seasonally specific drivers of postfire rates of NDVI recovery. Rates of postfire NDVI recovery were calculated for both the GS and SCS for more than 12,500 burned points across the western United States. Points were partitioned into faster and slower rates of NDVI recovery using thresholds derived from field plot data (n = 230) and their associated rates of NDVI recovery. We found plots with conifer saplings had significantly higher SCS NDVI recovery rates relative to plots without conifer saplings, while plots with ≥50% grass/forbs/shrubs cover had significantly higher GS NDVI recovery rates relative to plots with <50%. GS rates of NDVI recovery were best predicted by burn severity and anomalies in postfire maximum temperature. SCS NDVI recovery rates were best explained by aridity and growing degree days. This study is the most extensive effort, to date, to track postfire forest recovery across the western United States. Isolating patterns and drivers of evergreen recovery from deciduous recovery will enable improved characterization of forest ecological condition across large spatial scales.

Isolating Anthropogenic Wetland Loss by Concurrently Tracking Inundation and Land Cover Disturbance across the Mid-Atlantic Region, U.S.

Global trends in wetland degradation and loss have created an urgency to monitor wetland extent, as well as track the distribution and causes of wetland loss. Satellite imagery can be used to monitor wetlands over time, but few efforts have attempted to distinguish anthropogenic wetland loss from climate-driven variability in wetland extent. We present an approach to concurrently track land cover disturbance and inundation extent across the Mid-Atlantic region, United States, using the Landsat archive in Google Earth Engine. Disturbance was identified as a change in greenness, using a harmonic linear regression approach, or as a change in growing season brightness. Inundation extent was mapped using a modified version of the U.S. Geological Survey's Dynamic Surface Water Extent (DSWE) algorithm. Annual (2015-2018) disturbance averaged 0.32% (1095 km2 year-1) of the study area per year and was most common in forested areas. While inundation extent showed substantial interannual variability, the co-occurrence of disturbance and declines in inundation extent represented a minority of both change types, totaling 109 km2 over the four-year period, and 186 km2, using the National Wetland Inventory dataset in place of the Landsat-derived inundation extent. When the annual products were evaluated with permitted wetland and stream fill points, 95% of the fill points were detected, with most found by the disturbance product (89%) and fewer found by the inundation decline product (25%). The results suggest that mapping inundation alone is unlikely to be adequate to find and track anthropogenic wetland loss. Alternatively, remotely tracking both disturbance and inundation can potentially focus efforts to protect, manage, and restore wetlands.

Efficient Delineation of Nested Depression Hierarchy in Digital Elevation Models for Hydrological Analysis Using Level-Set Methods.

In terrain analysis and hydrological modeling, surface depressions (or sinks) in a digital elevation model (DEM) are commonly treated as artifacts and thus filled and removed to create a depressionless DEM. Various algorithms have been developed to identify and fill depressions in DEMs during the past decades. However, few studies have attempted to delineate and quantify the nested hierarchy of actual depressions, which can provide crucial information for characterizing surface hydrologic connectivity and simulating the fill-merge-spill hydrological process. In this paper, we present an innovative and efficient algorithm for delineating and quantifying nested depressions in DEMs using the level-set method based on graph theory. The proposed level-set method emulates water level decreasing from the spill point along the depression boundary to the lowest point at the bottom of a depression. By tracing the dynamic topological changes (i.e., depression splitting/merging) within a compound depression, the level-set method can construct topological graphs and derive geometric properties of the nested depressions. The experimental results of two fine-resolution Light Detection and Ranging-derived DEMs show that the raster-based level-set algorithm is much more efficient (~150 times faster) than the vector-based contour tree method. The proposed level-set algorithm has great potential for being applied to large-scale ecohydrological analysis and watershed modeling.

Influence of multi-decadal land use, irrigation practices and climate on riparian corridors across the Upper Missouri River headwaters basin, Montana.

The Upper Missouri River headwaters (UMH) basin (36 400 km2) depends on its river corridors to support irrigated agriculture and world-class trout fisheries. We evaluated trends (1984-2016) in riparian wetness, an indicator of the riparian condition, in peak irrigation months (June, July and August) for 158 km2 of riparian area across the basin using the Landsat normalized difference wetness index (NDWI). We found that 8 of the 19 riparian reaches across the basin showed a significant drying trend over this period, including all three basin outlet reaches along the Jefferson, Madison and Gallatin rivers. The influence of upstream climate was quantified using per reach random forest regressions. Much of the interannual variability in the NDWI was explained by climate, especially by drought indices and annual precipitation, but the significant temporal drying trends persisted in the NDWI-climate model residuals, indicating that trends were not entirely attributable to climate. Over the same period we documented a basin-wide shift from 9 % of agriculture irrigated with center-pivot irrigation to 50 % irrigated with center-pivot irrigation. Riparian reaches with a drying trend had a greater increase in the total area with center-pivot irrigation (within reach and upstream from the reach) relative to riparian reaches without such a trend (p < 0.05). The drying trend, however, did not extend to river discharge. Over the same period, stream gages (n = 7) showed a positive correlation with riparian wetness (p < 0.05) but no trend in summer river discharge, suggesting that riparian areas may be more sensitive to changes in irrigation return flows relative to river discharge. Identifying trends in riparian vegetation is a critical precursor for enhancing the resiliency of river systems and associated riparian corridors.

The potential role of very high-resolution imagery to characterise lake, wetland and stream systems across the Prairie Pothole Region, United States.

Aquatic features critical to watershed hydrology range widely in size from narrow, shallow streams to large, deep lakes. In this study we evaluated wetland, lake, and river systems across the Prairie Pothole Region to explore where pan-sharpened high-resolution (PSHR) imagery, relative to Landsat imagery, could pro-vide additional data on surface water distribution and movement, missed by Landsat. We used the monthly Global Surface Water (GSW) Landsat product as well as surface water derived from Landsat imagery using a matched filtering algorithm (MF Landsat) to help consider how including partially inundated Landsat pixels as water influenced our findings. The PSHR outputs (and MF Landsat) were able to identify ~60-90% more surface water interactions between waterbodies, relative to the GSW Landsat product. However, regardless of Landsat source, by doc-umenting many smaller (<0.2 ha), inundated wetlands, the PSHR outputs modified our interpretation of wetland size distribution across the Prairie Pothole Region.

CONNECTIVITY OF STREAMS AND WETLANDS TO DOWNSTREAM WATERS: AN INTEGRATED SYSTEMS FRAMEWORK.

Interest in connectivity has increased in the aquatic sciences, partly because of its relevance to the Clean Water Act. This paper has two objectives: (1) provide a framework to understand hydrological, chemical, and biological connectivity, focusing on how headwater streams and wetlands connect to and contribute to rivers; and (2) review methods to quantify hydrological and chemical connectivity. Streams and wetlands affect river structure and function by altering material and biological fluxes to the river; this depends on two factors: (1) functions within streams and wetlands that affect material fluxes; and (2) connectivity (or isolation) from streams and wetlands to rivers that allows (or prevents) material transport between systems. Connectivity can be described in terms of frequency, magnitude, duration, timing, and rate of change. It results from physical characteristics of a system, e.g., climate, soils, geology, topography, and the spatial distribution of aquatic components. Biological connectivity is also affected by traits and behavior of the biota. Connectivity can be altered by human impacts, often in complex ways. Because of variability in these factors, connectivity is not constant but varies over time and space. Connectivity can be quantified with field-based methods, modeling, and remote sensing. Further studies using these methods are needed to classify and quantify connectivity of aquatic ecosystems and to understand how impacts affect connectivity.

Estimating restorable wetland water storage at landscape scales.

Globally, hydrologic modifications such as ditching and subsurface drainage have significantly reduced wetland water storage capacity (i.e., volume of surface water a wetland can retain) and consequent wetland functions. While wetland area has been well documented across many landscapes and used to guide restoration efforts, few studies have directly quantified the associated wetland storage capacity. Here, we present a novel raster-based approach to quantify both contemporary and potential (i.e., restorable) storage capacities of individual depressional basins across landscapes. We demonstrate the utility of this method by applying it to the Delmarva Peninsula, a region punctuated by both depressional wetlands and drainage ditches. Across the entire peninsula, we estimated that restoration (i.e., plugging ditches) could increase storage capacity by 80%. Focusing on an individual watershed, we found that over 59% of restorable storage capacity occurs within 20 m of the drainage network, and that 93% occurs within 1 m elevation of the drainage network. Our demonstration highlights widespread ditching in this landscape, spatial patterns of both contemporary and potential storage capacities, and clear opportunities for hydrologic restoration. In Delmarva and more broadly, our novel approach can inform targeted landscape-scale conservation and restoration efforts to optimize hydrologically mediated wetland functions.

Wetlands inform how climate extremes influence surface water expansion and contraction.

Effective monitoring and prediction of flood and drought events requires an improved understanding of how and why surface-water expansion and contraction in response to climate varies across space. This paper sought to (1) quantify how interannual patterns of surface-water expansion and contraction vary spatially across the Prairie Pothole Region (PPR) and adjacent Northern Prairie (NP) in the United States, and (2) explore how landscape characteristics influence the relationship between climate inputs and surface-water dynamics. Due to differences in glacial history, the PPR and NP show distinct patterns in regards to drainage development and wetland density, together providing a diversity of conditions to examine surface-water dynamics. We mapped surface-water extent across eleven Landsat path/rows representing the PPR and NP (images spanning 1985-2015). The PPR not only experienced a 2.6-fold increase of surface-water extent under median conditions relative to the NP, but also showed a 3.4-fold greater difference in surface-water extent between drought and deluge conditions. The relationship between surface-water extent and accumulated water availability (precipitation minus potential evapotranspiration) was quantified per watershed and statistically related to variables representing hydrology-related landscape characteristics (e.g., infiltration capacity, surface storage capacity, stream density). To investigate the influence stream-connectivity has on the rate at which surface water leaves a given location, we modeled stream-connected and stream-disconnected surface water separately. Stream-connected surface water showed a greater expansion with wetter climatic conditions in landscapes with greater total wetland area. Disconnected surface water showed a greater expansion with wetter climatic conditions in landscapes with higher wetland density, lower infiltration and less anthropogenic drainage. From these findings, we can expect that shifts in precipitation and evaporative demand will have uneven effects on surface-water quantity. Accurate predictions regarding the effect of climate change on surface-water quantity will require consideration of hydrology-related landscape characteristics including wetlands.

BIOTA CONNECT AQUATIC HABITATS THROUGHOUT FRESHWATER ECOSYSTEM MOSAICS.

Freshwater ecosystems are linked at various spatial and temporal scales by movements of biota adapted to life in water. We review the literature on movements of aquatic organisms that connect different types of freshwater habitats, focusing on linkages from streams and wetlands to downstream waters. Here, streams, wetlands, rivers, lakes, ponds, and other freshwater habitats are viewed as dynamic freshwater ecosystem mosaics (FEMs) that collectively provide the resources needed to sustain aquatic life. Based on existing evidence, it is clear that biotic linkages throughout FEMs have important consequences for biological integrity and biodiversity. All aquatic organisms move within and among FEM components, but differ in the mode, frequency, distance, and timing of their movements. These movements allow biota to recolonize habitats, avoid inbreeding, escape stressors, locate mates, and acquire resources. Cumulatively, these individual movements connect populations within and among FEMs and contribute to local and regional diversity, resilience to disturbance, and persistence of aquatic species in the face of environmental change. Thus, the biological connections established by movement of biota among streams, wetlands, and downstream waters are critical to the ecological integrity of these systems. Future research will help advance our understanding of the movements that link FEMs and their cumulative effects on downstream waters.

Patterns and drivers for wetland connections in the Prairie Pothole Region, United States.

Ecosystem function in rivers, lakes and coastal waters depends on the functioning of upstream aquatic ecosystems, necessitating an improved understanding of watershed-scale interactions including variable surface-water flows between wetlands and streams. As surface water in the Prairie Pothole Region expands in wet years, surface-water connections occur between many depressional wetlands and streams. Minimal research has explored the spatial patterns and drivers for the abundance of these connections, despite their potential to inform resource management and regulatory programs including the U.S. Clean Water Act. In this study, wetlands were identified that did not intersect the stream network, but were shown with Landsat images (1990-2011) to become merged with the stream network as surface water expanded. Wetlands were found to spill into or consolidate with other wetlands within both small (2-10 wetlands) and large (>100 wetlands) wetland clusters, eventually intersecting a stream channel, most often via a riparian wetland. These surface-water connections occurred over a wide range of wetland distances from streams (averaging 90-1400 m in different ecoregions). Differences in the spatial abundance of wetlands that show a variable surface-water connection to a stream were best explained by smaller wetland-to-wetland distances, greater wetland abundance, and maximum surface-water extent. This analysis demonstrated that wetland arrangement and surface water expansion are important mechanisms for depressional wetlands to connect to streams and provides a first step to understanding the frequency and abundance of these surface-water connections across the Prairie Pothole Region.

The influence of data characteristics on detecting wetland/stream surface-water connections in the Delmarva Peninsula, Maryland and Delaware.

The dependence of downstream waters on upstream ecosystems necessitates an improved understanding of watershed-scale hydrological interactions including connections between wetlands and streams. An evaluation of such connections is challenging when, (1) accurate and complete datasets of wetland and stream locations are often not available and (2) natural variability in surface-water extent influences the frequency and duration of wetland/stream connectivity. The Upper Choptank River watershed on the Delmarva Peninsula in eastern Maryland and Delaware is dominated by a high density of small, forested wetlands. In this analysis, wetland/stream surface water connections were quantified using multiple wetland and stream datasets, including headwater streams and depressions mapped from a lidar-derived digital elevation model. Surface-water extent was mapped across the watershed for spring 2015 using Landsat-8, Radarsat-2 and Worldview-3 imagery. The frequency of wetland/stream connections increased as a more complete and accurate stream dataset was used and surface-water extent was included, in particular when the spatial resolution of the imagery was finer (i.e., <10 m). Depending on the datasets used, 12% to 60% of wetlands by count (21% to 93% of wetlands by area) experienced surface-water interactions with streams during spring 2015. This translated into a range of 50% to 94% of the watershed contributing direct surface water runoff to streamflow. This finding suggests that our interpretation of the frequency and duration of wetland/stream connections will be influenced not only by the spatial and temporal characteristics of wetlands, streams and potential flowpaths, but also by the completeness, accuracy and resolution of input datasets.

The Role of Lake Expansion in Altering the Wetland Landscape of the Prairie Pothole Region, United States.

Interannual variation in lake extent is well documented in the Prairie Pothole Region, but the role of surface-water expansion, including lake expansion, in merging with and subsuming wetlands across the landscape has been minimally considered. We examined how the expansion of surface-water extent, in particular, the expansion of lakes across parts of the Prairie Pothole Region can alter landscape-level hydrologic connectivity among substantial numbers of previously surficially disconnected wetlands. Temporally static wetland, lake, and stream datasets were fused with temporally varying Landsat-derived surface-water extent maps (1990-2011) to quantify changes in surface-water connectivity. Under deluge conditions, lakes were found to create significantly larger complexes of surficially-connected wetlands relative to non-lake surface-water connections (e.g., only wetlands or wetlands and streams). Analysis of three specific lakes showed that lakes can merge with and subsume wetlands located kilometers to tens of kilometers from the National Wetland Inventory defined lake perimeter. As climate across the Prairie Pothole Region is highly variable, understanding historic patterns of surface-water expansion and contraction under drought-to-deluge conditions will be integral to predicting future effects of climate change on wetland function, loss and influence on other aquatic systems, including downstream waters.

On the causes of rising gross ecosystem productivity in a regenerating clearcut environment: leaf area vs. species composition.

Clearcutting a forest ecosystem can result in a drastic reduction of stand productivity. Despite the severity of this disturbance type, past studies have found that the productivity of young regenerating stands can quickly rebound, approaching that of mature undisturbed stands within a few years. One of the obvious reasons is increased leaf area (LA) with each year of recovery. However, a less obvious reason may be the variability in species composition and distribution during the natural regeneration process. The purpose of this study was to investigate to what extent the increase in gross ecosystem productivity (GEP), observed during the first 4 years of recovery in a naturally regenerating clearcut stand, was due to (i) an overall expansion of leaf area and (ii) an increase in the canopy's photosynthetic capacity stemming from either species compositional shifts or drift in physiological traits within species. We found that the multi-year rise in GEP following harvest was clearly attributed to the expansion of LA rather than a change in vegetation composition. Sizeable changes in the relative abundance of species were masked by remarkably similar leaf physiological attributes for a range of vegetation types present in this early-successional environment. Comparison of upscaled leaf-chamber estimates with eddy-covariance-based estimates of light-response curves revealed a broad consistency in both maximum photosynthetic capacity and quantum yield efficiency. The approaches presented here illustrate how chamber- and ecosystem-scale measurements of gas exchange can be blended with species-level LA data to draw conclusive inferences about changes in ecosystem processes over time in a highly dynamic environment.

Post-clearcut dynamics of carbon, water and energy exchanges in a midlatitude temperate, deciduous broadleaf forest environment.

Clearcutting and other forest disturbances perturb carbon, water, and energy balances in significant ways, with corresponding influences on Earth's climate system through biogeochemical and biogeophysical effects. Observations are needed to quantify the precise changes in these balances as they vary across diverse disturbances of different types, severities, and in various climate and ecosystem type settings. This study combines eddy covariance and micrometeorological measurements of surface-atmosphere exchanges with vegetation inventories and chamber-based estimates of soil respiration to quantify how carbon, water, and energy fluxes changed during the first 3 years following forest clearing in a temperate forest environment of the northeastern US. We observed rapid recovery with sustained increases in gross ecosystem productivity (GEP) over the first three growing seasons post-clearing, coincident with large and relatively stable net emission of CO2 because of overwhelmingly large ecosystem respiration. The rise in GEP was attributed to vegetation changes not environmental conditions (e.g., weather), but attribution to the expansion of leaf area vs. changes in vegetation composition remains unclear. Soil respiration was estimated to contribute 44% of total ecosystem respiration during summer months and coarse woody debris accounted for another 18%. Evapotranspiration also recovered rapidly and continued to rise across years with a corresponding decrease in sensible heat flux. Gross short-wave and long-wave radiative fluxes were stable across years except for strong wintertime dependence on snow covered conditions and corresponding variation in albedo. Overall, these findings underscore the highly dynamic nature of carbon and water exchanges and vegetation composition during the regrowth following a severe forest disturbance, and sheds light on both the magnitude of such changes and the underlying mechanisms with a unique example from a temperate, deciduous broadleaf forest.