Search

Extreme flow conditions in river discharge have far-reaching environmental and economic consequences. The retention of surface water in lakes, wetlands, and floodplains can potentially moderate these extreme flows by modifying the timing, duration, and magnitude of flow generation. However, efforts to characterize the impact of surface water storage on river discharge have been limited in geographic extent. In this analysis, a suite of hydrologic signatures, quantifying components of watershed flow regimes, was calculated from daily discharge at 72 gaged watersheds across the conterminous United States. Random forest models were developed to explain variability in six hydrologic signatures related to flashiness and high and low flow conditions. In addition to traditionally considered variables such as climate, land cover, topography, and geology, a novel remote sensing (Sentinel-1 & 2) approach was used to study the contribution of surface water storage dynamics to each signature's variability. While climate variables explained much of the variability in the hydrologic signatures, models for five of the six signatures showed an improvement in explanatory power when landscape characteristics were added. Automated variable selection is part of the modeling process and can be indicative of the relative importance of certain variables over others. When all variables were considered, four of the six signature models selected remotely sensed inundation variables. The amount of semi-permanent and permanent floodplain inundation, for example, was both negatively correlated with, and showed the greatest variable importance for wet season flashiness. Further, increases in seasonal floodplain inundation were positively correlated with increases in peak flows. This suggests that the storage of surface water on floodplains is relevant to both flashiness and high flow signatures. In addition, spatial variability in the amount of semi-permanent and permanent non-floodplain water helped explain variability in the baseflow index. These findings suggest that watershed surface water storage dynamics explain a portion of streamflow signature variability. The results underscore the need for protection and restoration of surface water storage systems, such as wetlands, across watersheds. [ABSTRACT FROM AUTHOR]

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. Plain Language Summary: Climate change is expected to impact rivers, lakes, and wetlands. In this effort we used multiple satellites to track monthly surface water extent (2017–2021) for 32 sites across the central United States. The average amount of surface water was highly variable across sites. Individual statistical models, relating meteorological variables to surface water extent, were developed for each site. The models were then updated with climate change adjusted variables. Most sites were projected to show a decline in surface water extent, with drying most consistent across the southeast central, southwest central, and midwest central U.S. Projected declines under dry conditions at these sites suggest greater impacts of drought events are likely. 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, likely attributable to projected increases in winter and spring precipitation. Key Points: Surface water across the central U.S. responds to episodic, seasonal, and interannual variability in water availabilityMost sites showed projected declines in surface water extents under RCP4.5 and 8.5, peaking in summer‐fall and increasing drought impactsSelect sites, concentrated in the north central U.S., projected increases in surface water, associated with greater precipitation [ABSTRACT FROM AUTHOR]

Wildfires and housing development have increased since the 1990s, presenting unique challenges for wildfire management. However, it is unclear how the relative influences of housing growth and changing wildfire occurrence have altered risk to homes, or the potential for wildfire to threaten homes. We used a random forests model to predict burn probability in relation to weather variables at 1‐km resolution and monthly intervals from 1990 through 2019 in the Southern Rocky Mountains ecoregion. We quantified risk by combining the predicted burn probabilities with decadal housing density. We then compared the predicted burn probabilities and risk across the study area with observed values and quantified trends. Finally, we evaluated how housing growth and changes in burn probability influenced risk individually and combined. Fires burned 9055 km2 and exposed more than 8500 homes from 1990 to 2019. Observed burned area increased 632% from the 1990s to the 2000s, which combined with housing growth, resulted in a 1342% increase in homes exposed. Increases continued in the 2010s but at lower rates; burned area by 65% and exposure by 32%. The random forests model had excellent fit and high correlation with observations (AUC = 0.88 and r = 0.9). Observed values were within the 95% uncertainty interval for all years except 2016 (burned area) and 2000 (exposure). However, our model overpredicted in years with low observed burned area and underpredicted in years with high observed burned area. Overpredictions in risk resulted in lower rates of change in predicted risk compared with change in observed exposure. Increases in risk between the 1990s and 2000s were primarily due to warmer and drier weather conditions and secondarily because of housing growth. However, increases between the 2000s and 2010s were primarily due to housing growth. Our modeling approach identifies spatial and temporal patterns of wildfire potential and risk, which is critical information to guide decision‐making. Because the drivers behind risk shift over time, strategies to mitigate risk may need to account for multiple drivers simultaneously. [ABSTRACT FROM AUTHOR]

Background: Remotely sensed burned area products are critical to support fire modelling, policy, and management but often require further processing before use. Aim: We calculated fire history metrics from the Landsat Burned Area Product (1984–2020) across the conterminous U.S. (CONUS) including (1) fire frequency, (2) time since last burn (TSLB), (3) year of last burn, (4) longest fire-free interval, (5) average fire interval length, and (6) contemporary fire return interval (cFRI). Methods: Metrics were summarised by ecoregion and land ownership, and related to historical and cheatgrass datasets to demonstrate further applications of the products. Key results: The proportion burned ranged from 0.7% in the Northeast Mixed Woods to 74.1% in the Kansas Flint Hills. The Flint Hills and Temperate Prairies showed the highest burn frequency, while the Flint Hills and the Sierra Nevada and Klamath Mountains showed the shortest TSLB. Compared to private, public land had greater burned area (19 of 31 ecoregions) and shorter cFRI (25 of 31 ecoregions). Conclusions: Contemporary fire history metrics can help characterise recent fire regimes across CONUS. Implications: In regions with frequent fire, comparison of contemporary with target fire regimes or invasive species datasets enables the efficient incorporation of burned area data into decision-making. We present contemporary fire history metrics for the conterminous United States (CONUS) derived from 37 years of the Landsat Burned Area Product (1984–2020) and provide examples of how these metrics can inform decision-making. Fire regimes are diverse across CONUS, but most ecoregions showed more burning on public than private land. [ABSTRACT FROM AUTHOR]

Prescribed fires and wildfires are common in wetland ecosystems across the Southeastern United States. However, the wetland burned area has been chronically underestimated across the region due to (1) spectral confusion between open water and burned area, (2) rapid post-fire vegetation regrowth, and (3) high annual precipitation limiting clear-sky satellite observations. We developed a machine learning algorithm specifically for burned area in wetlands, and applied the algorithm to the Sentinel-2 archive (2016-2019) across the Southeastern US (>290,000 km²). Combining Landsat-8 imagery with Sentinel-2 increased the annual clear-sky observation count from 17 to 46 in 2016 and from 16 to 78 in 2019. When validated withWorldView imagery, the Sentinel-2 burned area had a 29% and 30% omission and commission rates of error for burned area, respectively, compared to the US Geological Survey Landsat-8 Burned Area Product (L8 BA), which had a 47% and 8% omission and commission rate of error, respectively. The Sentinel-2 algorithm and the L8 BA mapped burned area within 78% and 60% of wetland fire perimeters (n = 555) compiled from state and federal agencies, respectively. This analysis demonstrated the potential of Sentinel-2 to support efforts to track the burned area, especially across challenging ecosystem types, such as wetlands. [ABSTRACT FROM AUTHOR]

Over the past 2 decades, several global burned area products have been produced and released to the public. However, the accuracy assessment of such products largely depends on the availability of reliable reference data that currently do not exist on a global scale or whose production require a high level of dedication of project resources. The important lack of reference data for the validation of burned area products is addressed in this paper. We provide the Burned Area Reference Database (BARD), the first publicly available database created by compiling existing reference BA (burned area) datasets from different international projects. BARD contains a total of 2661 reference files derived from Landsat and Sentinel-2 imagery. All those files have been checked for internal quality and are freely provided by the authors. To ensure database consistency, all files were transformed to a common format and were properly documented by following metadata standards. The goal of generating this database was to give BA algorithm developers and product testers reference information that would help them to develop or validate new BA products. BARD is freely available at 10.21950/BBQQU7 (Franquesa et al., 2020). [ABSTRACT FROM AUTHOR]

This study explores the viability of using simulated monthly runoff as a proxy for landscape‐scale surface‐depression storage processes simulated by the United States Geological Survey's National Hydrologic Model (NHM) infrastructure across the conterminous United States (CONUS). Two different temporal resolution model codes (daily and monthly) were run in the NHM with the same spatial discretization. Simulated values of daily surface‐depression storage (treated as a decimal fraction of maximum volume) as computed by the daily Precipitation‐Runoff Modeling System (NHM‐PRMS) and normalized runoff (0 to 1) as computed by the Monthly Water Balance Model (NHM‐MWBM) were aggregated to monthly and annual values for each hydrologic response unit (HRU) in the CONUS geospatial fabric (HRU; n = 109,951) and analyzed using Spearman's rank correlation test. Correlations between simulated runoff and surface‐depression storage aggregated to monthly and annual values were compared to identify where which time scale had relatively higher correlation values across the CONUS. Results show Spearman's rank values >0.75 (highly correlated) for the monthly time scale in 28,279 HRUs (53.35%) compared to the annual time scale in 41,655 HRUs (78.58%). The geographic distribution of HRUs with highly correlated monthly values show areas where surface‐depression storage features are known to be common (e.g., Prairie Pothole Region, Florida). Research Impact Statement: Correlation between simulated runoff and surface depression storage using continental‐extent models to improve understanding and representation of surface‐depression storage processes. [ABSTRACT FROM AUTHOR]

The Upper Missouri River headwaters (UMH) basin (36 400 km 2) 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 km 2 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. [ABSTRACT FROM AUTHOR]

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 provide 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 documenting many smaller (<0.2 ha), inundated wetlands, the PSHR outputs modified our interpretation of wetland size distribution across the Prairie Pothole Region. [ABSTRACT FROM AUTHOR]

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. Research Impact Statement: An efficient level‐set method for delineating and characterizing nested depressions in digital elevation models (DEMs) for terrain analysis and hydrological modeling. [ABSTRACT FROM AUTHOR]

Abstract. Landsat Normalised Difference Vegetation Index (NDVI) is commonly used to monitor post-fire green-up; however, most studies do not distinguish new growth of conifer from deciduous or herbaceous species, despite potential consequences for local climate, carbon and wildlife. We found that dual season (growing and snow cover) NDVI improved our ability to distinguish conifer tree presence and density. We then examined the post-fire pattern (1984–2017) in Landsat NDVI for fires that occurred a minimum of 20 years ago (1986–1997). Points were classified into four categories depending on whether NDVI, 20 years post-fire, had returned to pre-fire values in only the growing season, only under snow cover, in both seasons or neither. We found that each category of points showed distinct patterns of NDVI change that could be used to characterise the average pre-fire and post-fire vegetation condition Of the points analysed, 43% showed a between-season disagreement if NDVI had returned to pre-fire values, suggesting that using dual-season NDVI can modify our interpretations of post-fire conditions. We also found an improved correlation between 5- and 20-year NDVI change under snow cover, potentially attributable to snow masking fast-growing herbaceous vegetation. This study suggests that snow-cover Landsat imagery can enhance characterisations of forest recovery following fire. We used a time series approach (1985–2017) to evaluate the post-fire (1986–1997) change in greenness for conifer forests across the United States Rocky Mountains. We found that examining trends in post-fire greenness in both the growing season and snow-cover season allowed us to isolate conifer species, improving our understanding of both the pre-fire and post-fire vegetation condition. [ABSTRACT FROM AUTHOR]

Landscape carbon (C) flux estimates help assess the ability of terrestrial ecosystems to buffer further increases in anthropogenic carbon dioxide (CO2) emissions. Advances in remote sensing have led to coarse‐scale estimates of gross primary productivity (GPP; e.g., MODIS 17), yet efforts to develop spatial respiration products are lacking. Here we demonstrate a method to predict growing season soil respiration at a regional scale in a mixed subalpine ecosystem. We related field measurements (n = 396) of growing season soil respiration mostly from subalpine forests in the Southern Rocky Mountains ecoregion to a suite of biophysical predictors using a Random Forest model (30‐m pixel size). We found that Landsat Enhanced Vegetation Index, growing season aridity index, temperature, precipitation, elevation, and slope aspect explained spatiotemporal variability in soil respiration. Our model had a psuedo‐r2 of 0.45 and root‐mean‐square error of roughly one quarter of the mean value of respiration. Predicted growing season soil respiration across the region was remarkably consistent across 2004, 2005, and 2006 (150‐day sums of 542.8, 544.3, and 536.5 g C/m2, respectively). Yet we observed substantial variability in spatial patterns of soil respiration predictions that varied among years, suggesting that our method is sensitive to changes in respiration drivers. Mean predicted growing season soil respiration was 73% of MODIS GPP, while predicted soil respiration was generally within 20% of nocturnal net ecosystem exchange from nearby eddy covariance towers. Thus, geospatial and remotely sensed data sets can be used to estimate soil respiration at landscape scales. Plain Language Summary: Soil respiration returns carbon dioxide back to the atmosphere and is an important part of the carbon cycle, but estimates of soil respiration across large landscapes are difficult to come by. Soil respiration is sensitive to changes in climate and vegetation, which are available as mapped data products, thanks to remote sensing and geospatial technology. We developed a statistical model that mapped soil respiration across three forests and an entire region based on climate and vegetation spatial data. While this work was limited to subalpine forests in the Southern Rocky Mountains, our method can be used in other ecosystems to better understand how ecosystems interact with atmospheric carbon dioxide. Key Points: Subalpine soil respiration was estimated across the Southern Rocky Mountains using 396 point measurements, Landsat Enhanced Vegetation Index, climate, and terrainPredicted soil respiration compared reasonably well to eddy covariance nocturnal respiration and MODIS GPPThis method shows promise for large‐scale estimates of soil respiration, a large component of the terrestrial carbon cycle [ABSTRACT FROM AUTHOR]

Interpretations of post-fire condition and rates of vegetation recovery can influence management priorities, actions and perception of latent risks from landslides and floods. In this study, we used the Waldo Canyon fire (2012, Colorado Springs, Colorado, USA) as a case study to explore how a time series (2011–2016) of high-resolution images can be used to delineate burn extent and severity, as well as quantify post-fire vegetation recovery. We applied an object-based approach to map burn severity and vegetation recovery using Worldview-2, Worldview-3 and QuickBird-2 imagery. The burned area was classified as 51% high, 20% moderate and 29% low burn-severity. Across the burn extent, the shrub cover class showed a rapid recovery, resprouting vigorously within 1 year, whereas 4 years post-fire, areas previously dominated by conifers were divided approximately equally between being classified as dominated by quaking aspen saplings with herbaceous species in the understorey or minimally recovered. Relative to using a pixel-based Normalised Difference Vegetation Index (NDVI), our object-based approach showed higher rates of revegetation. High-resolution imagery can provide an effective means to monitor post-fire site conditions and complement more prevalent efforts with moderate- and coarse-resolution sensors. We classified burn extent, severity, and vegetation recovery for a time series of high-resolution images over the Waldo Canyon fire, Colorado, USA. High-resolution (2 m) imagery is increasingly available and can be used to enhance our ability, relative to coarser resolution sources of imagery, to map fire characteristics and recovery. [ABSTRACT FROM AUTHOR]

Abstract: 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. [ABSTRACT FROM AUTHOR]

Abstract: 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) briefly 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. [ABSTRACT FROM AUTHOR]

Abstract: Connectivity is a fundamental but highly dynamic property of watersheds. Variability in the types and degrees of aquatic ecosystem connectivity presents challenges for researchers and managers seeking to accurately quantify its effects on critical hydrologic, biogeochemical, and biological processes. However, protecting natural gradients of connectivity is key to protecting the range of ecosystem services that aquatic ecosystems provide. In this featured collection, we review the available evidence on connections and functions by which streams and wetlands affect the integrity of downstream waters such as large rivers, lakes, reservoirs, and estuaries. The reviews in this collection focus on the types of waters whose protections under the U.S. Clean Water Act have been called into question by U.S. Supreme Court cases. We synthesize 40+ years of research on longitudinal, lateral, and vertical fluxes of energy, material, and biota between aquatic ecosystems included within the Act's frame of reference. Many questions about the roles of streams and wetlands in sustaining downstream water integrity can be answered from currently available literature, and emerging research is rapidly closing data gaps with exciting new insights into aquatic connectivity and function at local, watershed, and regional scales. Synthesis of foundational and emerging research is needed to support science‐based efforts to provide safe, reliable sources of fresh water for present and future generations. [ABSTRACT FROM AUTHOR]

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 used Landsat imagery to characterize variability in surface water extent across 11 Landsat path/rows representing the PPR and NP (images spanned 1985-2015). The PPR not only experienced a 2.6-fold greater surface water extent under median conditions relative to the NP, but also showed a 3.4-fold greater change 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, but lower total wetland density. 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 wetland storage and arrangement. [ABSTRACT FROM AUTHOR]

Abstract: In this study, the Precipitation‐Runoff Modelling System (PRMS) was used to simulate changes in surface‐water depression storage in the 1,126‐km2 Upper Pipestem Creek basin located within the Prairie Pothole Region of North Dakota, USA. The Prairie Pothole Region is characterized by millions of small water bodies (or surface‐water depressions) that provide numerous ecosystem services and are considered an important contribution to the hydrologic cycle. The Upper Pipestem PRMS model was extracted from the U.S. Geological Survey's (USGS) National Hydrologic Model (NHM), developed to support consistent hydrologic modelling across the conterminous United States. The Geospatial Fabric database, created for the USGS NHM, contains hydrologic model parameter values derived from datasets that characterize the physical features of the entire conterminous United States for 109,951 hydrologic response units. Each hydrologic response unit in the Geospatial Fabric was parameterized using aggregated surface‐water depression area derived from the National Hydrography Dataset Plus, an integrated suite of application‐ready geospatial datasets. This paper presents a calibration strategy for the Upper Pipestem PRMS model that uses normalized lake elevation measurements to calibrate the parameters influencing simulated fractional surface‐water depression storage. Results indicate that inclusion of measurements that give an indication of the change in surface‐water depression storage in the calibration procedure resulted in accurate changes in surface‐water depression storage in the water balance. Regionalized parameterization of the USGS NHM will require a proxy for change in surface‐storage to accurately parameterize surface‐water depression storage within the USGS NHM. [ABSTRACT FROM AUTHOR]

Abstract: 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. [ABSTRACT FROM AUTHOR]

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 used Landsat imagery to characterize variability in surface-water extent across eleven Landsat path/rows representing the PPR and NP (images spanned 1985-2015). The PPR not only experienced a 2.6-fold greater surface-water extent under median conditions relative to the NP, but also showed a 3.4-fold greater change 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. [ABSTRACT FROM AUTHOR]

The U.S. Geological Survey has produced the Landsat Burned Area Essential Climate Variable (BAECV) product for the conterminous United States (CONUS), which provides wall-to-wall annual maps of burned area at 30 m resolution (1984-2015). Validation is a critical component in the generation of such remotely sensed products. Previous efforts to validate the BAECV relied on a reference dataset derived from Landsat, which was effective in evaluating the product across its timespan but did not allow for consideration of inaccuracies imposed by the Landsat sensor itself. In this effort, the BAECV was validated using 286 high-resolution images, collected from GeoEye-1, QuickBird-2,Worldview-2 and RapidEye satellites. A disproportionate sampling strategy was utilized to ensure enough burned area pixels were collected. Errors of omission and commission for burned area averaged 22±4% and 48±3%, respectively, across CONUS. Errors were lowest across the western U.S. The elevated error of commission relative to omission was largely driven by patterns in the Great Plains which saw low errors of omission (13±13%) but high errors of commission (70±5%) and potentially a region-growing function included in the BAECV algorithm. While the BAECV reliably detected agricultural fires in the Great Plains, it frequently mapped tilled areas or areas with low vegetation as burned. Landscape metrics were calculated for individual fire events to assess the influence of image resolution (2 m, 30 m and 500 m) on mapping fire heterogeneity. As the spatial detail of imagery increased, fire events were mapped in a patchier manner with greater patch and edge densities, and shape complexity, which can influence estimates of total greenhouse gas emissions and rates of vegetation recovery. The increasing number of satellites collecting high-resolution imagery and rapid improvements in the frequency with which imagery is being collected means greater opportunities to utilize these sources of imagery for Landsat product validation. [ABSTRACT FROM AUTHOR]

Natural variability in surface-water extent and associated characteristics presents a challenge to gathering timely, accurate information, particularly in environments that are dominated by small and/or forested wetlands. This study mapped inundation extent across the Upper Choptank River Watershed on the Delmarva Peninsula, occurring within both Maryland and Delaware. We integrated six quad-polarized Radarsat-2 images, Worldview-3 imagery, and an enhanced topographic wetness index in a random forest model. Output maps were filtered using light detection and ranging (lidar)-derived depressions to maximize the accuracy of forested inundation extent. Overall accuracy within the integrated and filtered model was 94.3%, with 5.5% and 6.0% errors of omission and commission for inundation, respectively. Accuracy of inundation maps obtained using Radarsat-2 alone were likely detrimentally affected by less than ideal angles of incidence and recent precipitation, but were likely improved by targeting the period between snowmelt and leaf-out for imagery collection. Across the six Radarsat-2 dates, filtering inundation outputs by lidar-derived depressions slightly elevated errors of omission for water (+1.0%), but decreased errors of commission (-7.8%), resulting in an average increase of 5.4% in overall accuracy. Depressions were derived from lidar datasets collected under both dry and average wetness conditions. Although antecedent wetness conditions influenced the abundance and total area mapped as depression, the two versions of the depression datasets showed a similar ability to reduce error in the inundation maps. Accurate mapping of surface water is critical to predicting and monitoring the effect of human-induced change and interannual variability on water quantity and quality. [ABSTRACT FROM AUTHOR]

Warmer conditions over the past two decades have contributed to rapid expansion of bark beetle outbreaks killing millions of trees over a large fraction of western United States ( US) forests. These outbreaks reduce plant productivity by killing trees and transfer carbon from live to dead pools where carbon is slowly emitted to the atmosphere via heterotrophic respiration which subsequently feeds back to climate change. Recent studies have begun to examine the local impacts of bark beetle outbreaks in individual stands, but the full regional carbon consequences remain undocumented for the western US. In this study, we quantify the regional carbon impacts of the bark beetle outbreaks taking place in western US forests. The work relies on a combination of postdisturbance forest regrowth trajectories derived from forest inventory data and a process-based carbon cycle model tracking decomposition, as well as aerial detection survey ( ADS) data documenting the regional extent and severity of recent outbreaks. We find that biomass killed by bark beetle attacks across beetle-affected areas in western US forests from 2000 to 2009 ranges from 5 to 15 Tg C yr−1 and caused a reduction of net ecosystem productivity ( NEP) of about 6.1-9.3 Tg C y−1 by 2009. Uncertainties result largely from a lack of detailed surveys of the extent and severity of outbreaks, calling out a need for improved characterization across western US forests. The carbon flux legacy of 2000-2009 outbreaks will continue decades into the future (e.g., 2040-2060) as committed emissions from heterotrophic respiration of beetle-killed biomass are balanced by forest regrowth and accumulation. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

In this paper, we address various challenges in multi-pedestrian and vehicle tracking in high-resolution aerial imagery by intensive evaluation of a number of traditional and Deep Learning based Single- and Multi-Object Tracking methods. We also describe our proposed Deep Learning based Multi-Object Tracking method AerialMPTNet that fuses appearance, temporal, and graphical information using a Siamese Neural Network, a Long Short-Term Memory, and a Graph Convolutional Neural Network module for more accurate and stable tracking. Moreover, we investigate the influence of the Squeeze-and-Excitation layers and Online Hard Example Mining on the performance of AerialMPTNet. To the best of our knowledge, we are the first to use these two for regression-based Multi-Object Tracking. Additionally, we studied and compared the L 1 and Huber loss functions. In our experiments, we extensively evaluate AerialMPTNet on three aerial Multi-Object Tracking datasets, namely AerialMPT and KIT AIS pedestrian and vehicle datasets. Qualitative and quantitative results show that AerialMPTNet outperforms all previous methods for the pedestrian datasets and achieves competitive results for the vehicle dataset. In addition, Long Short-Term Memory and Graph Convolutional Neural Network modules enhance the tracking performance. Moreover, using Squeeze-and-Excitation and Online Hard Example Mining significantly helps for some cases while degrades the results for other cases. In addition, according to the results, L 1 yields better results with respect to Huber loss for most of the scenarios. The presented results provide a deep insight into challenges and opportunities of the aerial Multi-Object Tracking domain, paving the way for future research. [ABSTRACT FROM AUTHOR]

The Visible Infrared Imaging Radiometer Suite (VIIRS) onboard the Suomi National Polar-orbiting Partnership (SNPP) satellite has been a reliable source of ocean color data products, including five moderate (M) bands and one imagery (I) band normalized water-leaving radiance spectra nLw(λ). The spatial resolutions of the M-band and I-band nLw(λ) are 750 m and 375 m, respectively. With the technique of convolutional neural network (CNN), the M-band nLw(λ) imagery can be super-resolved from 750 m to 375 m spatial resolution by leveraging the high spatial resolution features of I1-band nLw(λ) data. However, it is also important to enhance the spatial resolution of VIIRS-derived chlorophyll-a (Chl-a) concentration and the water diffuse attenuation coefficient at the wavelength of 490 nm (Kd(490)), as well as other biological and biogeochemical products. In this study, we describe our effort to derive high-resolution Kd(490) and Chl-a data based on super-resolved nLw(λ) images at the VIIRS five M-bands. To improve the network performance over extremely turbid coastal oceans and inland waters, the networks are retrained with a training dataset including ocean color data from the Bohai Sea, Baltic Sea, and La Plata River Estuary, covering water types from clear open oceans to moderately turbid and highly turbid waters. The evaluation results show that the super-resolved Kd(490) image is much sharper than the original one, and has more detailed fine spatial structures. A similar enhancement of finer structures is also found in the super-resolved Chl-a images. Chl-a filaments are much sharper and thinner in the super-resolved image, and some of the very fine spatial features that are not shown in the original images appear in the super-resolved Chl-a imageries. The networks are also applied to four other coastal and inland water regions. The results show that super-resolution occurs mainly on pixels of Chl-a and Kd(490) features, especially on the feature edges and locations with a large spatial gradient. The biases between the original M-band images and super-resolved high-resolution images are small for both Chl-a and Kd(490) in moderately to extremely turbid coastal oceans and inland waters, indicating that the super-resolution process does not change the mean values of the original images. [ABSTRACT FROM AUTHOR]

Accurate spatial information of agricultural fields is important for providing actionable information to farmers, managers, and policymakers. On the other hand, the automated detection of field boundaries is a challenging task due to their small size, irregular shape and the use of mixed-cropping systems making field boundaries vaguely defined. In this paper, we propose a strategy for field boundary detection based on the fully convolutional network architecture called ResU-Net. The benefits of this model are two-fold: first, residual units ease training of deep networks. Second, rich skip connections within the network could facilitate information propagation, allowing us to design networks with fewer parameters but better performance in comparison with the traditional U-Net model. An extensive experimental analysis is performed over the whole of Denmark using Sentinel-2 images and comparing several U-Net and ResU-Net field boundary detection algorithms. The presented results show that the ResU-Net model has a better performance with an average F1 score of 0.90 and average Jaccard coefficient of 0.80 in comparison to the U-Net model with an average F1 score of 0.88 and an average Jaccard coefficient of 0.77. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The Delmarva Peninsula in the eastern United States is partially characterized by thousands of small, forested, depressional wetlands that are highly sensitive to weather variability and climate change, but provide critical ecosystem services. Due to the relatively small size of these depressional wetlands and their occurrence under forest canopy cover, it is very challenging to map their inundation status based on existing remote sensing data and traditional classification approaches. In this study, we applied a state-of-the-art U-Net semantic segmentation network to map forested wetland inundation in the Delmarva area by integrating leaf-off WorldView-3 (WV3) multispectral data with fine spatial resolution light detection and ranging (lidar) intensity and topographic data, including a digital elevation model (DEM) and topographic wetness index (TWI). Wetland inundation labels generated from lidar intensity were used for model training and validation. The wetland inundation map results were also validated using field data, and compared to the U.S. Fish and Wildlife Service National Wetlands Inventory (NWI) geospatial dataset and a random forest output from a previous study. Our results demonstrate that our deep learning model can accurately determine inundation status with an overall accuracy of 95% (Kappa = 0.90) compared to field data and high overlap (IoU = 70%) with lidar intensity-derived inundation labels. The integration of topographic metrics in deep learning models can improve the classification accuracy for depressional wetlands. This study highlights the great potential of deep learning models to improve the accuracy of wetland inundation maps through use of high-resolution optical and lidar remote sensing datasets. [ABSTRACT FROM AUTHOR]

Degradation of streams and associated riparian habitat across the Missouri River Headwaters Basin has motivated several stream restoration projects across the watershed. Many of these projects install a series of beaver dam analogues (BDAs) to aggrade incised streams, elevate local water tables, and create natural surface water storage by reconnecting streams with their floodplains. Satellite imagery can provide a spatially continuous mechanism to monitor the effects of these in-stream structures on stream surface area. However, remote sensing-based approaches to map narrow (e.g., <5 m wide) linear features such as streams have been under-developed relative to efforts to map other types of aquatic systems, such as wetlands or lakes. We mapped pre- and post-restoration (one to three years post-restoration) stream surface area and riparian greenness at four stream restoration sites using Worldview-2 and 3 images as well as a QuickBird-2 image. We found that panchromatic brightness and eCognition-based outputs (0.5 m resolution) provided high-accuracy maps of stream surface area (overall accuracy ranged from 91% to 99%) for streams as narrow as 1.5 m wide. Using image pairs, we were able to document increases in stream surface area immediately upstream of BDAs as well as increases in stream surface area along the restoration reach at Robb Creek, Alkali Creek and Long Creek (South). Although Long Creek (North) did not show a net increase in stream surface area along the restoration reach, we did observe an increase in riparian greenness, suggesting increased water retention adjacent to the stream. As high-resolution imagery becomes more widely collected and available, improvements in our ability to provide spatially continuous monitoring of stream systems can effectively complement more traditional field-based and gage-based datasets to inform watershed management. [ABSTRACT FROM AUTHOR]

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.

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.