Automated Extraction of Surface Water Extent from Sentinel-1 Data

Automated in-season rice crop mapping using Sentinel time-series data and Google Earth Engine: A case study in climate-risk prone Bangladesh

High-resolution mapping of rice fields is crucial for understanding and managing rice cultivation in countries like Bangladesh, particularly in the face of climate change. Rice is a vital crop, cultivated in small scale farms that contributes significantly to the economy and food security in Bangladesh. Accurate mapping can facilitate improved rice production, the development of sustainable agricultural management policies, and formulation of strategies for adapting to climatic risks. To address the need for timely and accurate rice mapping, we developed a framework specifically designed for the diverse environmental conditions in Bangladesh. We utilized Sentinel-1 and Sentinel-2 time-series data to identify transplantation and peak seasons and employed the multi-Otsu automatic thresholding approach to map rice during the peak season (April-May). We also compared the performance of a random forest (RF) classifier with the multi-Otsu approach using two different data combinations: D1, which utilizes data from the transplantation and peak seasons (D1 RF) and D2, which utilizes data from the transplantation to the harvest seasons (D2 RF). Our results demonstrated that the multi-Otsu approach achieved an overall classification accuracy (OCA) ranging from 61.18% to 94.43% across all crop zones. The D2 RF showed the highest mean OCA (92.15%) among the fourteen crop zones, followed by D1 RF (89.47%) and multi-Otsu (85.27%). Although the multi-Otsu approach had relatively lower OCA, it proved effective in accurately mapping rice areas prior to harvest, eliminating the need for training samples that can be challenging to obtain during the growing season. In-season rice area maps generated through this framework are crucial for timely decision-making regarding adaptive management in response to climatic stresses and forecasting area-wide productivity. The scalability of our framework across space and time makes it particularly suitable for addressing field data scarcity challenges in countries like Bangladesh and offers the potential for future operationalization.

Anthropogenic impacts on urban blue space and its reciprocal effect on human and socio-ecological health

Quantifying anthropogenic impacts on blue space (BS) and its effect on human and socio-ecological health was least explored. The present study aimed to do this in reference to the urban BS transformation scenario of Eastern India. To measure BS transformation, Landsat image-based water indices were run from 1990 to 2021. Anthropogenic impact score (AIS) and 7 components scores of 78 selected BS on 70 parameters related data driven from the field. Total 345 respondents were taken for human and socio-ecological health assessment. For this, depression (DEP), anxiety (ANX), stress (STR), physical activities (PA), social capital (SC), therapeutic landscape (TL) and environment building (EB) parameters were taken. The result exhibited that BS was reduced. About 50% of urban core BS was reported highly impacted. Human and socio-ecological health was identified as good in proximity to BS, but it was observed better in the cases of larger peripheral BS. AIS on BS was found to be positively associated with mental health (0.47-0.63) and negatively associated with PA, SC, TL and EB (-0.50 to -0.90). Standard residual in ordinary least square was reported low (-1.5 to 1.5) in 95% BS. Therefore, BS health restoration and management is crucial for sustaining the living environment.

Use of remote sensing-based pressure-state-response framework for the spatial ecosystem health assessment in Langfang, China

Land use/land cover changes are occurring at an unprecedented rate and scale because of the economic development that has jeopardized the ecosystem's health. Ecosystem health should be studied and monitored at spatiotemporal scale to promote sustainable development and ecological civilization. The goal of this study was to assess the spatial ecosystem health of Langfang at the city and administrative levels using city's regional characteristics. Remote sensing-based pressure-state-response (PSR) framework, analytical hierarchy process (AHP), and principal component analysis (PCA) were utilized for spatial ecosystem health index (SEHI) formulation, indicator weighting, and indicator selection in several epochs (1990, 2003, 2013, and 2021), respectively. SEHI was formulated by combining subindices of pressure, state and response. The spatial ecosystem pressure index (SEIP) identified that the pressure was increasing on the ecosystem. In contrast, the spatial ecosystem state index (SEIS) pointed out an improvement in the state of the ecosystem since 1990. The worst state of the ecosystem was observed for the year 2013. The spatial ecosystem response index (SEIR) indicated that the response of the ecosystem towards the exerted pressures and states remained variable; however, it was reasonably good in 1990. All the administrative units of Langfang were associated with a healthy score for the spatial ecosystem health index (SEHI) for 1990 (pre-industrialization epoch), while the SEHI significantly reduced in 2013 (industrialization epoch) however improved for the later epochs (circular economy and ecological civilization epoch). The SEHI was moderately healthy for Dachang, Dacheng, Guan, Guangyang, and Yongqing while relatively healthy for the remaining administrative units in 2021. SEHI identified that spatial health has been improving since 2003 though not reaching the 1990's level for Langfang. Therefore, efforts should be focused on minimizing pressure and stabilizing the state to improve the spatial ecosystem health of Langfang. The developed SEHI can assist policymakers in analyzing regional health, identifying development strategies, driving environmental restoration, and quantifying needed changes.

Effects of optical and radar satellite observations within Google Earth Engine on soil organic carbon prediction models in Spain

The modeling and mapping of soil organic carbon (SOC) has advanced through the rapid growth of Earth observation data (e.g., Sentinel) collection and the advent of appropriate tools such as the Google Earth Engine (GEE). However, the effects of differing optical and radar sensors on SOC prediction models remain uncertain. This research aims to investigate the effects of different optical and radar sensors (Sentinel-1/2/3 and ALOS-2) on SOC prediction models based on long-term satellite observations on the GEE platform. We also evaluate the relative impact of four synthetic aperture radar (SAR) acquisition configurations (polarization mode, band frequency, orbital direction and time window) on SOC mapping with multiband SAR data from Spain. Twelve experiments involving different satellite data configurations, combined with 4027 soil samples, were used for building SOC random forest regression models. The results show that the synthesis mode and choice of satellite images, as well as the SAR acquisition configurations, influenced the model accuracy to varying degrees. Models based on SAR data involving cross-polarization, multiple time periods and "ASCENDING" orbits outperformed those involving copolarization, a single time period and "DESCENDING" orbits. Moreover, combining information from different orbital directions and polarization modes improved the soil prediction models. Among the SOC models based on long-term satellite observations, the Sentinel-3-based models (R² = 0.40) performed the best, while the ALOS-2-based model performed the worst. In addition, the predictive performance of MSI/Sentinel-2 (R² = 0.35) was comparable with that of SAR/Sentinel-1 (R² = 0.35); however, the combination (R² = 0.39) of the two improved the model performance. All the predicted maps involving Sentinel satellites had similar spatial patterns that were higher in northwest Spain and lower in the south. Overall, this study provides insights into the effects of different optical and radar sensors and radar system parameters on soil prediction models and improves our understanding of the potential of Sentinels in developing soil carbon mapping.

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

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.

The urgent need to develop a new grassland map in China: based on the consistency and accuracy of ten land cover products

Grasslands are the most dominant terrestrial ecosystem in China, but few national grassland maps have been generated. The grassland resource map produced in the 1980s is widely used as background data, but it has not been updated for almost 40 years. Therefore, a reliable map depicting the current spatial distribution of grasslands across the country is urgently needed. In this study, we evaluated the grassland consistency and accuracy of ten land cover datasets (GLC2000, GlobCover, CCI-LC, MCD12Q1, CLUD, GlobeLand30, GLC-FCS30, CGLS-LC100, CLCD, and FROM-GLC) for 2000, 2010, and 2020 based on extensive fieldwork. We concluded that the area of these ten grassland products ranges from 107.80×10⁴ to 332.46×10⁴ km², with CLCD and MCD12Q1 having the highest area consistency. The spatial and sample consistency is highest in the regions of east-central Inner Mongolia, the Qinghai-Tibet Plateau and northern Xinjiang, while the distribution of southern grasslands is scattered and differs considerably among the ten products. MCD12Q1 is significantly more accurate than the other nine products, with an overall accuracy (OA) reaching 77.51% and a kappa coefficient of 0.51; CLCD is slightly less accurate than MCD12Q1 (OA=73.02%, kappa coefficient=0.45) and is more conducive to the fine monitoring and management of grassland because of its 30-meter resolution. The highest accuracy of grassland was found in the Inner Mongolia-Ningxia region and Qinghai-Tibet Plateau, while the accuracy was worst in the southeastern region. In the future grassland mapping, cartographers should improve the accuracy of the grassland distribution in South China and regions where grassland is confused with forest, cropland and bare land. We specify the availability of valuable data in existing land cover datasets for China's grasslands and call for researchers and the government to actively produce a new generation of grassland maps.

Developing an integrated approach based on geographic object-based image analysis and convolutional neural network for volcanic and glacial landforms mapping

Rapid detection and mapping of landforms are crucially important to improve our understanding of past and presently active processes across the earth, especially, in complex and dynamic volcanoes. Traditional landform modeling approaches are labor-intensive and time-consuming. In recent years, landform mapping has increasingly been digitized. This study conducted an in-depth analysis of convolutional neural networks (CNN) in combination with geographic object-based image analysis (GEOBIA), for mapping volcanic and glacial landforms. Sentinel-2 image, as well as predisposing variables (DEM and its derivatives, e.g., slope, aspect, curvature and flow accumulation), were segmented using a multi-resolution segmentation algorithm, and relevant features were selected to define segmentation scales for each landform category. A set of object-based features was developed based on spectral (e.g., brightness), geometrical (e.g., shape index), and textural (grey level co-occurrence matrix) information. The landform modelling networks were then trained and tested based on labelled objects generated using GEOBIA and ground control points. Our results show that an integrated approach of GEOBIA and CNN achieved an ACC of 0.9685, 0.9780, 0.9614, 0.9767, 0.9675, 0.9718, 0.9600, and 0.9778 for dacite lava, caldera, andesite lava, volcanic cone, volcanic tuff, glacial circus, glacial valley, and suspended valley, respectively. The quantitative evaluation shows the highest performance (Accuracy > 0.9600 and cross-validation accuracy > 0.9400) for volcanic and glacial landforms and; therefore, is recommended for regional and large-scale landform mapping. Our results and the provided automatic workflow emphasize the potential of integrated GEOBIA and CNN for fast and efficient landform mapping as a first step in the earth's surface management.

Satellite Imagery-Based Identification of High-Risk Areas of Schistosome Intermediate Snail Hosts Spread after Flood

Snail intermediate host monitoring and control are essential for interrupting the parasitic disease schistosomiasis. Identifying large-scale high-risk areas of snail spread after floods has been greatly facilitated by remote sensing imagery. However, previous studies have usually assumed that all inundation areas carry snails and may have overestimated snail spread areas. Furthermore, these studies only used a single environmental factor to estimate the snail survival risk probability, failing to analyze multiple variables, to accurately distinguish the snail survival risk in the snail spread areas. This paper proposes a systematic framework for early monitoring of snail diffusion to accurately map snail spread areas from remote sensing imagery and enhance snail survival risk probability estimation based on the snail spread map. In particular, the flooded areas are extracted using the Sentinel-1 Dual-Polarized Water Index based on synthetic aperture radar images to map all-weather flooding areas. These flood maps are used to extract snail spread areas, based on the assumption that only inundation areas that spatially interacted with (i.e., are close to) the previous snail distribution regions before flooding are identified as snail spread areas, in order to reduce the misclassification in snail spread area identification. A multiple logistic regression model is built to analyze how various types of snail-related environmental factors, including the normalized difference vegetation index (NDVI), wetness, river and channel density, and landscape fractal dimension impact snail survival, and estimate its risk probabilities in snail spread area. An experiment was conducted in Jianghan Plain, China, where snails are predominantly linearly distributed along the tributaries and water channels of the middle and lower reaches of the Yangtze River. The proposed method could accurately map floods under clouds, and a total area of 231.5 km2 was identified as the snail spread area. The snail survival risk probabilities were thus estimated. The proposed method showed a more refined snail spread area and a more reliable degree of snail survival risk compared with those of previous studies. Thus, it is an efficient way to accurately map all-weather snail spread and survival risk probabilities, which is helpful for schistosomiasis interruption.

Rapid Extreme Tropical Precipitation and Flood Inundation Mapping Framework (RETRACE): Initial Testing for the 2021–2022 Malaysia Flood

The 2021–2022 flood is one of the most serious flood events in Malaysian history, with approximately 70,000 victims evacuated daily, 54 killed and total losses up to MYR 6.1 billion. From this devastating event, we realized the lack of extreme precipitation and flood inundation information, which is a common problem in tropical regions. Therefore, we developed a Rapid Extreme TRopicAl preCipitation and flood inundation mapping framEwork (RETRACE) by utilizing: (1) a cloud computing platform, the Google Earth Engine (GEE); (2) open-source satellite images from missions such as Global Precipitation Measurement (GPM), Sentinel-1 SAR and Sentinel-2 optical satellites; and (3) flood victim information. The framework was demonstrated with the 2021–2022 Malaysia flood. The preliminary results were satisfactory with an optimal threshold of five for flood inundation mapping using the Sentinel-1 SAR data, as the accuracy of inundated floods was up to 70%. Extreme daily precipitation of up to 230 mm/day was observed and resulted in an inundated area of 77.43 km2 in Peninsular Malaysia. This framework can act as a useful tool for local authorities and scientists to retrace the extreme precipitation and flood information in a relatively short period for flood management and mitigation strategy development.

Water Body Mapping Using Long Time Series Sentinel-1 SAR Data in Poyang Lake

Mapping water bodies with a high accuracy is necessary for water resource assessment, and mapping them rapidly is necessary for flood monitoring. Poyang Lake is the largest freshwater lake in China, and its wetland is one of the most important in the world. Poyang Lake is affected by floods from the Yangtze River basin every year, and the fluctuation of the water area and water level directly or indirectly affects the ecological environment of Poyang Lake. Synthetic Aperture Radar (SAR) is particularly suitable for large-scale water body mapping, as SAR allows data acquisition regardless of illumination and weather conditions. The two-satellite Sentinel-1 constellation, providing C-Band SAR data, passes over the Poyang Lake about five times a month. With its high temporal-spatial resolution, the Sentinel-1 SAR data can be used to accurately monitor the water body. After acquiring all the Sentinel-1 (1A and 1B) SAR data, to ensure the consistency of data processing, we propose the use of a Python and SeNtinel Application Platform (SNAP)-based engine (SARProcMod) to process the data and construct a Poyang Lake Sentinel-1 SAR dataset with a 10 m resolution. To extract water body information from Sentinel-1 SAR data, we propose an automatic classification engine based on a modified U-Net convolutional neural network (WaterUNet), which classifies all data using artificial sample datasets with a high validation accuracy. The results show that the maximum and minimum water areas in our study area were 2714.08 km2 on 20 July 2020, and 634.44 km2 on 4 January 2020. Compared to the water level data from the Poyang gauging station, the water area was highly correlated with the water level, with the correlation coefficient being up to 0.92 and the R2 from quadratic polynomial fitting up to 0.88; thus, the resulting relationship results can be used to estimate the water area or water level of Poyang Lake. According to the results, we can conclude that Sentinel-1 SAR and WaterUNet are very suitable for water body monitoring as well as emergency flood mapping.

Analysis of Environmental and Atmospheric Influences in the Use of SAR and Optical Imagery from Sentinel-1, Landsat-8, and Sentinel-2 in the Operational Monitoring of Reservoir Water Level

In this work, we aim to evaluate the feasibility and operational limitations of using Sentinel-1 synthetic aperture radar (SAR) data to monitor water levels in the Poço da Cruz reservoir from September 2016–September 2020, in the semi-arid region of northeast Brazil. To segment water/non-water features, SAR backscattering thresholding was carried out via the graphical interpretation of backscatter coefficient histograms. In addition, surrounding environmental effects on SAR polarization thresholds were investigated by applying wavelet analysis, and the Landsat-8 and Sentinel-2 normalized difference water index (NDWI) and modified normalized difference water index (MNDWI) were used to compare and discuss the SAR results. The assessment of the observed and estimated water levels showed that (i) SAR accuracy was equivalent to that of NDWI/Landsat-8; (ii) optical image accuracy outperformed SAR image accuracy in inlet branches, where the complexity of water features is higher; and (iii) VV polarization outperformed VH polarization. The results confirm that SAR images can be suitable for operational reservoir monitoring, offering a similar accuracy to that of multispectral indices. SAR threshold variations were strongly correlated to the normalized difference vegetation index (NDVI), the soil moisture variations in the reservoir depletion zone, and the prior precipitation quantities, which can be used as a proxy to predict cross-polarization (VH) and co-polarization (VV) thresholds. Our findings may improve the accuracy of the algorithms designed to automate the extraction of water levels using SAR data, either in isolation or combined with multispectral images.

Monitoring Post-Flood Recovery of Croplands Using the Integrated Sentinel-1/2 Imagery in the Yangtze-Huai River Basin

The increasingly frequent flooding imposes tremendous and long-lasting damages to lives and properties in impoverished rural areas. Rapid, accurate, and large-scale flood mapping is urgently needed for flood management, and to date has been successfully implemented benefiting from the advancement in remote sensing and cloud computing technology. Yet, the effects of agricultural emergency response to floods have been limitedly evaluated by satellite-based remote sensing, resulting in biased post-flood loss assessments. Addressing this challenge, this study presents a method for monitoring post-flood agricultural recovery using Sentinel-1/2 imagery, tested in three flood-affected main grain production areas, in the middle and lower Yangtze and Huai River, China. Our results indicated that 33~72% of the affected croplands were replanted and avoided total crop failures in summer 2020. Elevation, flood duration, crop rotation scheme, and flooding emergency management affect the post-flood recovery performance. The findings also demonstrate rapid intervention measures adjusted to local conditions could reduce the agricultural failure cost from flood disasters to a great extent. This study provides a new alternative for comprehensive disaster loss assessment in flood-prone agricultural regions, which will be insightful for worldwide flood control and management.

Integrating SAR and Optical Remote Sensing for Conservation-Targeted Wetlands Mapping

The Prairie Pothole Region (PPR) contains numerous depressional wetlands known as potholes that provide habitats for waterfowl and other wetland-dependent species. Mapping these wetlands is essential for identifying viable waterfowl habitat and conservation planning scenarios, yet it is a challenging task due to the small size of the potholes, and the presence of emergent vegetation. This study develops an open-source process within the Google Earth Engine platform for mapping the spatial distribution of wetlands through the integration of Sentinel-1 C-band SAR (synthetic aperture radar) data with high-resolution (10-m) Sentinel-2 bands. We used two machine-learning algorithms (random forest (RF) and support vector machine (SVM)) to identify wetlands across the study area through supervised classification of the multisensor composite. We trained the algorithms with ground truth data provided through field studies and aerial photography. The accuracy was assessed by comparing the predicted and actual wetland and non-wetland classes using statistical coefficients (overall accuracy, Kappa, sensitivity, and specificity). For this purpose, we used four different out-of-sample test subsets, including the same year, next year, small vegetated, and small non-vegetated test sets to evaluate the methods on different spatial and temporal scales. The results were also compared to Landsat-derived JRC surface water products, and the Sentinel-2-derived normalized difference water index (NDWI). The wetlands derived from the RF model (overall accuracy 0.76 to 0.95) yielded favorable results, and outperformed the SVM, NDWI, and JRC products in all four testing subsets. To provide a further characterization of the potholes, the water bodies were stratified based on the presence of emergent vegetation using Sentinel-2-derived NDVI, and, after excluding permanent water bodies, using the JRC surface water product. The algorithm presented in the study is scalable and can be adopted for identifying wetlands in other regions of the world.

Effects of Meteorological Parameters on Surface Water Loss in Burdur Lake, Turkey over 34 Years Landsat Google Earth Engine Time-Series

The current work aims to examine the effect of meteorological parameters as well as the temporal variation on the Burdur Lake surface water body areas in Turkey. The data for the surface area of Burdur Lake was collected over 34 years between 1984 and 2019 by Google Earth Engine. The seasonal variation in the water bodies area was collected using our proposed extraction method and 570 Landsat images. The reduction in the total area of the lake was observed between 206.6 km2 in 1984 to 125.5 km2 in 2019. The vegetation cover and meteorological parameters collected that affect the observed variation of surface water bodies for the same area include precipitation, evapotranspiration, albedo, radiation, and temperature. The selected meteorological variables influence the reduction in lake area directly during various seasons. Correlations showed a strong positive or negative significant relationship between the reduction and the selected meteorological variables. A factor analysis provided three components that explain 82.15% of the total variation in the data. The data provide valuable references for decision makers to develop sustainable agriculture and industrial water use policies to preserve water resources as well as cope with potential climate changes.

Semi-Automated Semantic Segmentation of Arctic Shorelines Using Very High-Resolution Airborne Imagery, Spectral Indices and Weakly Supervised Machine Learning Approaches

Precise coastal shoreline mapping is essential for monitoring changes in erosion rates, surface hydrology, and ecosystem structure and function. Monitoring water bodies in the Arctic National Wildlife Refuge (ANWR) is of high importance, especially considering the potential for oil and natural gas exploration in the region. In this work, we propose a modified variant of the Deep Neural Network based U-Net Architecture for the automated mapping of 4 Band Orthorectified NOAA Airborne Imagery using sparsely labeled training data and compare it to the performance of traditional Machine Learning (ML) based approaches—namely, random forest, xgboost—and spectral water indices—Normalized Difference Water Index (NDWI), and Normalized Difference Surface Water Index (NDSWI)—to support shoreline mapping of Arctic coastlines. We conclude that it is possible to modify the U-Net model to accept sparse labels as input and the results are comparable to other ML methods (an Intersection-over-Union (IoU) of 94.86% using U-Net vs. an IoU of 95.05% using the best performing method).

Estimation of Infiltration Volumes and Rates in Seasonally Water-Filled Topographic Depressions Based on Remote-Sensing Time Series

In semi-arid ecoregions of temperate zones, focused snowmelt water infiltration in topographic depressions is a key, but imperfectly understood, groundwater recharge mechanism. Routine monitoring is precluded by the abundance of depressions. We have used remote-sensing data to construct mass balances and estimate volumes of temporary ponds in the Tambov area of Russia. First, small water bodies were automatically recognized in each of a time series of high-resolution Planet Labs images taken in April and May 2021 by object-oriented supervised classification. A training set of water pixels defined in one of the latest images using a small unmanned aerial vehicle enabled high-confidence predictions of water pixels in the earlier images (Cohen’s Κ = 0.99). A digital elevation model was used to estimate the ponds’ water volumes, which decreased with time following a negative exponential equation. The power of the exponent did not systematically depend on the pond size. With adjustment for estimates of daily Penman evaporation, function-based interpolation of the water bodies’ areas and volumes allowed calculation of daily infiltration into the depression beds. The infiltration was maximal (5–40 mm/day) at onset of spring and decreased with time during the study period. Use of the spatially variable infiltration rates improved steady-state shallow groundwater simulations.

Remote Sensing of Wetlands in the Prairie Pothole Region of North America

The Prairie Pothole Region (PPR) of North America is an extremely important habitat for a diverse range of wetland ecosystems that provide a wealth of socio-economic value. This paper describes the ecological characteristics and importance of PPR wetlands and the use of remote sensing for mapping and monitoring applications. While there are comprehensive reviews for wetland remote sensing in recent publications, there is no comprehensive review about the use of remote sensing in the PPR. First, the PPR is described, including the wetland classification systems that have been used, the water regimes that control the surface water and water levels, and the soil and vegetation characteristics of the region. The tools and techniques that have been used in the PPR for analyses of geospatial data for wetland applications are described. Field observations for ground truth data are critical for good validation and accuracy assessment of the many products that are produced. Wetland classification approaches are reviewed, including Decision Trees, Machine Learning, and object versus pixel-based approaches. A comprehensive description of the remote sensing systems and data that have been employed by various studies in the PPR is provided. A wide range of data can be used for various applications, including passive optical data like aerial photographs or satellite-based, Earth-observation data. Both airborne and spaceborne lidar studies are described. A detailed description of Synthetic Aperture RADAR (SAR) data and research are provided. The state of the art is the use of multi-source data to achieve higher accuracies and hybrid approaches. Digital Surface Models are also being incorporated in geospatial analyses to separate forest and shrub and emergent systems based on vegetation height. Remote sensing provides a cost-effective mechanism for mapping and monitoring PPR wetlands, especially with the logistical difficulties and cost of field-based methods. The wetland characteristics of the PPR dictate the need for high resolution in both time and space, which is increasingly possible with the numerous and increasing remote sensing systems available and the trend to open-source data and tools. The fusion of multi-source remote sensing data via state-of-the-art machine learning is recommended for wetland applications in the PPR. The use of such data promotes flexibility for sensor addition, subtraction, or substitution as a function of application needs and potential cost restrictions. This is important in the PPR because of the challenges related to the highly dynamic nature of this unique region.

Monitoring Drought through the Lens of Landsat: Drying of Rivers during the California Droughts

Water scarcity during severe droughts has profound hydrological and ecological impacts on rivers. However, the drying dynamics of river surface extent during droughts remains largely understudied. Satellite remote sensing enables surveys and analyses of rivers at fine spatial resolution by providing an alternative to in-situ observations. This study investigates the seasonal drying dynamics of river extent in California where severe droughts have been occurring more frequently in recent decades. Our methods combine the use of Landsat-based Global Surface Water (GSW) and global river bankful width databases. As an indirect comparison, we examine the monthly fractional river extent (FrcSA) in 2071 river reaches and its correlation with streamflow at co-located USGS gauges. We place the extreme 2012–2015 drought into a broader context of multi-decadal river extent history and illustrate the extraordinary change between during- and post-drought periods. In addition to river extent dynamics, we perform statistical analyses to relate FrcSA with the hydroclimatic variables obtained from the National Land Data Assimilation System (NLDAS) model simulation. Results show that Landsat provides consistent observation over 90% of area in rivers from March to October and is suitable for monitoring seasonal river drying in California. FrcSA reaches fair (>0.5) correlation with streamflow except for dry and mountainous areas. During the 2012–2015 drought, 332 river reaches experienced their lowest annual mean FrcSA in the 34 years of Landsat history. At a monthly scale, FrcSA is better correlated with soil water in more humid areas. At a yearly scale, summer mean FrcSA is increasingly sensitive to winter precipitation in a drier climate; and the elasticity is also reduced with deeper ground water table. Overall, our study demonstrates the detectability of Landsat on the river surface extent in an arid region with complex terrain. River extent in catchments of deficient water storage is likely subject to higher percent drop in a future climate with longer, more frequent droughts.

Sentinel-1&2 Multitemporal Water Surface Detection Accuracies, Evaluated at Regional and Reservoirs Level

Water stock monitoring is a major issue for society on a local and global scale. Sentinel-1&2 satellites provide frequent acquisitions to track water surface dynamics, proxy variables to enable water surface volume monitoring. How do we combine such observations along time for each sensor? What advantages and disadvantages of single-date, monthly or time-windowed estimations? In this context, we analysed the impact of merging information through different types and lengths of time-windows. Satellite observations were processed separately on optical (Sentinel-2) and radar (Sentinel-1) water detectors at 10 m resolution. The analysis has been applied at two scales. First, validating with 26 large scenes (110 × 110 km) in different climatic zones in France, time-windows yielded an improvement on radar detection (F1-score improved from 0.72 to 0.8 for 30 days on average logic) while optical performances remained stable (F1-score 0.89). Second, validating reservoir area estimations with 29 instrumented reservoirs (20–1250 ha), time-windows presented in all cases an improvement on both optical and radar error for any window length (5–30 days). The mean relative absolute error in optical area detection improved from 16.9% on single measurements to 12.9% using 15 days time-windows, and from 22.15% to 15.1% in radar detection). Regarding reservoir filling rates, we identified an increased negative bias for both sensors when the reservoir is nearly full. This work helped to compare accuracies of separate optical and radar capabilities, where optical statistically outperforms radar at both local and large scale to the detriment of less frequent measurements. Furthermore, we propose a geomorphological indicator of reservoirs to predict the quality of radar area monitoring (R2 = 0.58). In conclusion, we suggest the use of time-windows on operational water mapping or reservoir monitoring systems, using 10–20 days time-windows with average logic, providing more frequent and faster information to water managers in periods of crisis (e.g., water shortage) compared to monthly estimations.

Impervious Surfaces Mapping at City Scale by Fusion of Radar and Optical Data through a Random Forest Classifier

Urbanization increases the amount of impervious surfaces, making accurate information on spatial and temporal expansion trends essential; the challenge is to develop a cost- and labor-effective technique that is compatible with the assessment of multiple geographical locations in developing countries. Several studies have identified the potential of remote sensing and multiple source information in impervious surface quantification. Therefore, this study aims to fuse datasets from the Sentinel 1 and 2 Satellites to map the impervious surfaces of nine Pakistani cities and estimate their growth rates from 2016 to 2020 utilizing the random forest algorithm. All bands in the optical and radar images were resampled to 10 m resolution, projected to same coordinate system and geometrically aligned to stack into a single product. The models were then trained, and classifications were validated with land cover samples from Google Earth’s high-resolution images. Overall accuracies of classified maps ranged from 85% to 98% with the resultant quantities showing a strong linear relationship (R-squared value of 0.998) with the Copernicus Global Land Services data. There was up to 9% increase in accuracy and up to 12 % increase in kappa coefficient from the fused data with respect to optical alone. A McNemar test confirmed the superiority of fused data. Finally, the cities had growth rates ranging from 0.5% to 2.5%, with an average of 1.8%. The information obtained can alert urban planners and environmentalists to assess impervious surface impacts in the cities.

Estimation of the Monthly Dynamics of Surface Water in Wetlands from Satellite and Secondary Hydro-Climatological Data

Satellites produce valuable information for studying the surface water in wetlands, but in many cases the period covered, the spatial resolution and/or the revisit frequency is not enough to produce long historical series. In this paper we propose a novel method which uses regression models that include climatic and hydrological variables to complete the satellite information. We used this method in the Lagunas de Ruidera wetland (Spain). We approached the monthly dynamic of the surface water for a long period (1984–2015). Information from LANDSAT (30-m resolution) and MODIS (250-m resolution) satellites were tested but, due to the size of some lagoons, only the LANDSAT approach produced satisfactory results. An ensemble of regression models based on hydro-climatological explanatory variables was defined to complete the gaps in the monthly surface water. It showed a root mean squared error of around 476 pixels (0.4 Km2) in the cross-validation analysis. Our analysis showed that the explanatory variables with a more significant participation in the regression ensemble are the aquifer discharge, the effective precipitation and the surface water from the previous month. From January to June, the mean surface water in Lagunas de Ruidera is around 4.3 Km2. In summer a reduction of around 13% of the surface water can be observed, which is recovered during the autumn.

Sentinel-1 and 2 Time-Series for Vegetation Mapping Using Random Forest Classification: A Case Study of Northern Croatia

Land-cover (LC) mapping in a morphologically heterogeneous landscape area is a challenging task since various LC classes (e.g., crop types in agricultural areas) are spectrally similar. Most research is still mostly relying on optical satellite imagery for these tasks, whereas synthetic aperture radar (SAR) imagery is often neglected. Therefore, this research assessed the classification accuracy using the recent Sentinel-1 (S1) SAR and Sentinel-2 (S2) time-series data for LC mapping, especially vegetation classes. Additionally, ancillary data, such as texture features, spectral indices from S1 and S2, respectively, as well as digital elevation model (DEM), were used in different classification scenarios. Random Forest (RF) was used for classification tasks using a proposed hybrid reference dataset derived from European Land Use and Coverage Area Frame Survey (LUCAS), CORINE, and Land Parcel Identification Systems (LPIS) LC database. Based on the RF variable selection using Mean Decrease Accuracy (MDA), the combination of S1 and S2 data yielded the highest overall accuracy (OA) of 91.78%, with a total disagreement of 8.22%. The most pertinent features for vegetation mapping were GLCM Mean and Variance for S1, NDVI, along with Red and SWIR band for S2, whereas the digital elevation model produced major classification enhancement as an input feature. The results of this study demonstrated that the aforementioned approach (i.e., RF using a hybrid reference dataset) is well-suited for vegetation mapping using Sentinel imagery, which can be applied for large-scale LC classifications.

Mapping Outburst Floods Using a Collaborative Learning Method Based on Temporally Dense Optical and SAR Data: A Case Study with the Baige Landslide Dam on the Jinsha River, Tibet

Outburst floods resulting from giant landslide dams can cause devastating damage to hundreds or thousands of kilometres of a river. Accurate and timely delineation of flood inundated areas is essential for disaster assessment and mitigation. There have been significant advances in flood mapping using remote sensing images in recent years, but little attention has been devoted to outburst flood mapping. The short-duration nature of these events and observation constraints from cloud cover have significantly challenged outburst flood mapping. This study used the outburst flood of the Baige landslide dam on the Jinsha River on 3 November 2018 as an example to propose a new flood mapping method that combines optical images from Sentinel-2, synthetic aperture radar (SAR) images from Sentinel-1 and a Digital Elevation Model (DEM). First, in the cloud-free region, a comparison of four spectral indexes calculated from time series of Sentinel-2 images indicated that the normalized difference vegetation index (NDVI) with the threshold of 0.15 provided the best separation flooded area. Subsequently, in the cloud-covered region, an analysis of dual-polarization RGB false color composites images and backscattering coefficient differences of Sentinel-1 SAR data were found an apparent response to ground roughness’s changes caused by the flood. We carried out the flood range prediction model based on the random forest algorithm. Training samples consisted of 13 feature vectors obtained from the Hue-Saturation-Value color space, backscattering coefficient differences/ratio, DEM data, and a label set from the flood range prepared from Sentinel-2 images. Finally, a field investigation and confusion matrix tested the prediction accuracy of the end-of-flood map. The overall accuracy and Kappa coefficient were 92.3%, 0.89 respectively. The full extent of the outburst floods was successfully obtained within five days of its occurrence. The multi-source data merging framework and the massive sample preparation method with SAR images proposed in this paper, provide a practical demonstration for similar machine learning applications using remote sensing.

An Optical and SAR Based Fusion Approach for Mapping Surface Water Dynamics over Mainland China

Earth Observation (EO) data is a critical information source for mapping and monitoring water resources over large inaccessible regions where hydrological in-situ networks are sparse. In this paper, we present a simple yet robust method for fusing optical and Synthetic Aperture Radar (SAR) data for mapping surface water dynamics over mainland China. This method uses a multivariate logistic regression model to estimate monthly surface water extent over a four-year period (2017 to 2020) from the combined usages of Sentinel-1, Sentinel-2 and Landsat-8 imagery. Multi-seasonal high-resolution images from the Chinese Gaofen satellites are used as a reference for an independent validation showing a high degree of agreement (overall accuracy 94%) across a diversity of climatic and physiographic regions demonstrating potential scalability beyond China. Through inter-comparison with similar global scale products, this paper further shows how this new mapping technique provides improved spatio-temporal characterization of inland water bodies, and for better capturing smaller water bodies (< 0.81 ha in size). The relevance of the results is discussed, and we find this new enhanced monitoring approach has the potential to advance the use of Earth observation for water resource management, planning and reporting.

Characterizing Wetland Inundation and Vegetation Dynamics in the Arctic Coastal Plain Using Recent Satellite Data and Field Photos

Arctic wetlands play a critical role in the global carbon cycle and are experiencing disproportionate impacts from climate change. Even though Alaska hosts 65% of U.S. wetlands, less than half of the wetlands in Alaska have been mapped by the U.S. Fish and Wildlife Service National Wetlands Inventory (NWI) or other high-resolution wetlands protocols. The availability of time series satellite data and the development of machine learning algorithms have enabled the characterization of Arctic wetland inundation dynamics and vegetation types with limited ground data input. In this study, we built a semi-automatic process to generate sub-pixel water fraction (SWF) maps across the Coastal Plain of the Arctic National Wildlife Refuge (ANWR) in Alaska using random forest regression and 139 Sentinel-2 images taken in ice-free seasons from 2016 to 2019. With this, we characterized the seasonal dynamics of wetland inundation and explored their potential usage in determining NWI water regimes. The highest levels of surface water expression were detected in June, resulting from seasonal active layer thaw and snowmelt. Inundation was most variable in riverbeds, lake and pond margins, and depressional wetlands, where water levels fluctuate substantially between dry and wet seasons. NWI water regimes that indicate frequent inundation, such as permanently flooded wetlands, had high SWF values (SWF ≥ 90%), while those with infrequent inundation, such as temporarily flooded wetlands, had low SWF values (SWF < 10%). Vegetation types were also classified through the synergistic use of a vegetation index, water regimes, synthetic-aperture radar (SAR) data, topographic data, and a random forest classifier. The random forest classification algorithms demonstrated good performance in classifying Arctic wetland vegetation types, with an overall accuracy of 0.87. Compared with NWI data produced in the 1980s, scrub-shrub wetlands appear to have increased from 91 to 258 km2 over the last three decades, which is the largest percentage change (182%) among all vegetation types. However, additional field data are needed to confirm this shift in vegetation type. This study demonstrates the potential of using time series satellite data and machine learning algorithms in characterizing inundation dynamics and vegetation types of Arctic wetlands. This approach could aid in the creation and maintenance of wetland inventories, including the NWI, in Arctic regions and enable an improved understanding of long-term wetland dynamics.

Using remote sensing to assess peatland resilience by estimating soil surface moisture and drought recovery

Peatland areas provide a range of ecosystem services, including biodiversity, carbon storage, clean water, and flood mitigation, but many areas of peatland in the UK have been degraded through human land use including drainage. Here, we explore whether remote sensing can be used to monitor peatland resilience to drought. We take resilience to mean the rate at which a system recovers from perturbation; here measured literally as a recovery timescale of a soil surface moisture proxy from drought lowering. Our objectives were (1) to assess the reliability of Sentinel-1 Synthetic Aperture Radar (SAR) backscatter as a proxy for water table depth (WTD); (2) to develop a method using SAR to estimate below-ground (hydrological) resilience of peatlands; and (3) to apply the developed method to different sites and consider the links between resilience and land management. Our inferences of WTD from Sentinel-1 SAR data gave results with an average Pearson's correlation of 0.77 when compared to measured WTD values. The 2018 summer drought was used to assess resilience across three different UK peatland areas (Dartmoor, the Peak District, and the Flow Country) by considering the timescale of the soil moisture proxy recovery. Results show clear areas of lower resilience within all three study sites, which often correspond to areas of high drainage and may be particularly vulnerable to increasing drought severity/events under climate change. This method is applicable to monitoring peatland resilience elsewhere over larger scales, and could be used to target restoration work towards the most vulnerable areas.

Flood Mapping in Vegetated Areas Using an Unsupervised Clustering Approach on Sentinel-1 and -2 Imagery

The European Space Agency’s Sentinel-1 constellation provides timely and freely available dual-polarized C-band Synthetic Aperture Radar (SAR) imagery. The launch of these and other SAR sensors has boosted the field of SAR-based flood mapping. However, flood mapping in vegetated areas remains a topic under investigation, as backscatter is the result of a complex mixture of backscattering mechanisms and strongly depends on the wave and vegetation characteristics. In this paper, we present an unsupervised object-based clustering framework capable of mapping flooding in the presence and absence of flooded vegetation based on freely and globally available data only. Based on a SAR image pair, the region of interest is segmented into objects, which are converted to a SAR-optical feature space and clustered using K-means. These clusters are then classified based on automatically determined thresholds, and the resulting classification is refined by means of several region growing post-processing steps. The final outcome discriminates between dry land, permanent water, open flooding, and flooded vegetation. Forested areas, which might hide flooding, are indicated as well. The framework is presented based on four case studies, of which two contain flooded vegetation. For the optimal parameter combination, three-class F1 scores between 0.76 and 0.91 are obtained depending on the case, and the pixel- and object-based thresholding benchmarks are outperformed. Furthermore, this framework allows an easy integration of additional data sources when these become available.

A Multi-Scale Deep Neural Network for Water Detection from SAR Images in the Mountainous Areas

Water detection from Synthetic Aperture Radar (SAR) images has been widely utilized in various applications. However, it remains an open challenge due to the high similarity between water and shadow in SAR images. To address this challenge, a new end-to-end framework based on deep learning has been proposed to automatically classify water and shadow areas in SAR images. This end-to-end framework is mainly composed of three parts, namely, Multi-scale Spatial Feature (MSF) extraction, Multi-Level Selective Attention Network (MLSAN) and the Improvement Strategy (IS). Firstly, the dataset is input to MSF for multi-scale low-level feature extraction via three different methods. Then, these low-level features are fed into the MLSAN network, which contains the Encoder and Decoder. The Encoder aims to generate different levels of features using residual network of 101 layers. The Decoder extracts geospatial contextual information and fuses the multi-level features to generate high-level features that are further optimized by the IS. Finally, the classification is implemented with the Softmax function. We name the proposed framework as MSF-MLSAN, which is trained and tested using millimeter wave SAR datasets. The classification accuracy reaches 0.8382 and 0.9278 for water and shadow, respectively; while the overall Intersection over Union (IoU) is 0.9076. MSF-MLSAN demonstrates the success of integrating SAR domain knowledge and state-of-the-art deep learning techniques.

High Spatiotemporal Resolution Mapping of Surface Water in the Southwest Poyang Lake and Its Responses to Climate Oscillations

Accurately quantifying spatiotemporal changes in surface water is essential for water resources management, nevertheless, the dynamics of Poyang Lake surface water areas with high spatiotemporal resolution, as well as its responses to climate change, still face considerable uncertainties. Using the time series of Sentinel-1 images with 6- or 12-day intervals, the Sentinel-1 water index (SWI), and SWI-based water extraction model (SWIM) from 2015 to 2020 were used to document and study the short-term characteristics of southwest Poyang Lake surface water. The results showed that the overall accuracy of surface water area was satisfactory with an average of 91.92%, and the surface water area ranged from 129.06 km² on 2 March 2017 to 1042.57 km² on 17 July 2016, with significant intra- and inter-month variability. Within the 6-day interval, the maximum change of lake area was 233.42 km² (i.e., increasing from 474.70 km² up to 708.12 km²). We found that the correlation coefficient between the water area and the 45-day accumulated precipitation reached to 0.75 (p < 0.001). Moreover, a prediction model was built to predict the water area based on climate records. These results highlight the significance of high spatiotemporal resolution mapping for surface water in the erratic southwest Poyang Lake under a changing climate. The automated water extraction algorithm proposed in this study has potential applications in delineating surface water dynamics at broad geographic scales.

Comparing Sentinel-1 Surface Water Mapping Algorithms and Radiometric Terrain Correction Processing in Southeast Asia Utilizing Google Earth Engine

Satellite remote sensing plays an important role in the monitoring of surface water for historical analysis and near real-time applications. Due to its cloud penetrating capability, many studies have focused on providing efficient and high quality methods for surface water mapping using Synthetic Aperture Radar (SAR). However, few studies have explored the effects of SAR pre-processing steps used and the subsequent results as inputs into surface water mapping algorithms. This study leverages the Google Earth Engine to compare two unsupervised histogram-based thresholding surface water mapping algorithms utilizing two distinct pre-processed Sentinel-1 SAR datasets, specifically one with and one without terrain correction. The resulting surface water maps from the four different collections were validated with user-interpreted samples from high-resolution Planet Scope data. It was found that the overall accuracy from the four collections ranged from 92% to 95% with Cohen’s Kappa coefficients ranging from 0.7999 to 0.8427. The thresholding algorithm that samples a histogram based on water edge information performed best with a maximum accuracy of 95%. While the accuracies varied between methods it was found that there is no statistical significant difference between the errors of the different collections. Furthermore, the surface water maps generated from the terrain corrected data resulted in a intersection over union metrics of 95.8%–96.4%, showing greater spatial agreement, as compared to 92.3%–93.1% intersection over union using the non-terrain corrected data. Overall, it was found that algorithms using terrain correction yield higher overall accuracy and yielded a greater spatial agreement between methods. However, differences between the approaches presented in this paper were not found to be significant suggesting both methods are valid for generating accurate surface water maps. High accuracy surface water maps are critical to disaster planning and response efforts, thus results from this study can help inform SAR data users on the pre-processing steps needed and its effects as inputs on algorithms for surface water mapping applications.

Wetland Surface Water Detection from Multipath SAR Images Using Gaussian Process-Based Temporal Interpolation

Wetlands provide society with a myriad of ecosystem services, such as water storage, food sources, and flood control. The ecosystem services provided by a wetland are largely dependent on its hydrological dynamics. Constant monitoring of the spatial extent of water surfaces and the duration of flooding of a wetland is necessary to understand the impact of drought on the ecosystem services a wetland provides. Synthetic aperture radar (SAR) has the potential to reveal wetland dynamics. Multitemporal SAR image analysis for wetland monitoring has been extensively studied based on the advances of modern SAR missions. Unfortunately, most previous studies utilized monopath SAR images, which result in limited success. Tracking changes in individual wetlands remains a challenging task because several environmental factors, such as wind-roughened water, degrade image quality. In general, the data acquisition frequency is an important factor in time series analysis. We propose a Gaussian process-based temporal interpolation (GPTI) method that enables the synergistic use of SAR images taken from multiple paths. The proposed model is applied to a series of Sentinel-1 images capturing wetlands in Okanogan County, Washington State. Our experimental analysis demonstrates that the multiple path analysis based on the proposed method can extract seasonal changes more accurately than a single path analysis.

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.

A Pathway to the Automated Global Assessment of Water Level in Reservoirs with Synthetic Aperture Radar (SAR)

Global measurements of reservoir water levels are crucial for understanding Earth’s hydrological dynamics, especially in the context of global industrialization and climate change. Although radar altimetry has been used to measure the water level of some reservoirs with high accuracy, it is not yet feasible unless the water body is sufficiently large or directly located at the satellite’s nadir. This study proposes a gauging method applicable to a wide range of reservoirs using Sentinel–1 Synthetic Aperture Radar data and a digital elevation model (DEM). The method is straightforward to implement and involves estimating the mean slope–corrected elevation of points along the reservoir shoreline. We test the model on six case studies and show that the estimated water levels are accurate to around 10% error on average of independently verified values. This study represents a substantial step toward the global gauging of lakes and reservoirs of all sizes and in any location where a DEM is available.

Tropical Wetland (TropWet) Mapping Tool: The Automatic Detection of Open and Vegetated Waterbodies in Google Earth Engine for Tropical Wetlands

Knowledge of the location and extent of surface water and inundated vegetation is vital for a range of applications including flood risk management, biodiversity monitoring, quantifying greenhouse gas emissions, and mapping water-borne disease risk. Here, we present a new tool, TropWet, which enables users of all abilities to map wetlands in herbaceous dominated regions based on simple unmixing of optical Landsat satellite imagery in the Google Earth Engine. The results demonstrate transferability throughout the African continent with a high degree of accuracy (mean 91% accuracy, st. dev 2.6%, n = 10,800). TropWet demonstrated considerable improvements over existing globally available surface water datasets for mapping the extent of important wetlands like the Okavango, Botswana. TropWet was able to provide frequency inundation maps as an indicator of malarial mosquito aquatic habitat extent and persistence in Barotseland, Zambia. TropWet was able to map flood extent comparable to operational flood risk mapping products in the Zambezi Region, Namibia. Finally, TropWet was able to quantify the effects of the El Niño/Southern Oscillation (ENSO) events on the extent of photosynthetic vegetation and wetland extent across Southern Africa. These examples demonstrate the potential for TropWet to provide policy makers with crucial information to help make national, regional, or continental scale decisions regarding wetland conservation, flood/disease hazard mapping, or mitigation against the impacts of ENSO.

Automatic Surface Water Mapping Using Polarimetric SAR Data for Long-Term Change Detection

Mapping the distribution and persistence of surface water in a timely fashion has broad value for tracking dynamic events like flooding, and for monitoring the effects of climate and human activities on natural resource values and biodiversity. Traditionally, surface water is mapped from optical imagery using semi-automatic approaches. However, this process is time-consuming and the accuracy of results can vary among image interpreters. In recent years, Synthetic Aperture Radar (SAR) images have been increasingly used. Microwave signals sensitive to water content make SAR systems useful for mapping surface water, saturated soils, and flooded vegetation. In this study, a fully automatic method based on robust stepwise thresholding was developed to map and track the change in the extent of surface water using Polarimetric SAR data. The application of this method in both Radarsat-2 and Sentinel-1 data in central Ontario, Canada demonstrates that the developed robust stepwise thresholding approach could facilitate rapid mapping of open water areas with a promising accuracy of over 95%. In addition, the time-series extent of surface water extracted from May 2008 to August 2016 reveals the dynamic nature of surface inundation, and the trend was consistent with the local precipitation data.

Mapping Forested Wetland Inundation in the Delmarva Peninsula, USA Using Deep Convolutional Neural Networks

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.

Seasonal cycles of lakes on the Tibetan Plateau detected by Sentinel-1 SAR data

The long-term evolution of lakes on the Tibetan Plateau (TP) could be observed from Landsat series of satellite data since the 1970s. However, the seasonal cycles of lakes on the TP have received little attention due to high cloud contamination of the commonly-used optical images. In this study, for the first time, the seasonal cycle of lakes on the TP was detected using Sentinel-1 Synthetic Aperture Radar (SAR) data with a high repeat cycle. A total of approximately 6000 Level-1 scenes were obtained that covered all large lakes (>50 km²) in the study area. The images were extracted from stripmap (SM) and interferometric wide swath (IW) modes that had a pixel spacing of 40 m in the range and azimuth directions. The lake boundaries extracted from Sentinel-1 data using the algorithm developed in this study were in good agreement with in-situ measurements of lake shoreline, lake outlines delineated from the corresponding Landsat images in 2015 and lake levels for Qinghai Lake. Upon analysis, it was found that the seasonal cycles of lakes exhibited drastically different patterns across the TP. For example, large size lakes (>100 km²) reached their peaks in August-September while lakes with areas of 50-100 km² reached their peaks in early June-July. The peaks of seasonal cycles for endorheic lakes were more pronounced than those for exorheic lakes with flat peaks, and glacier-fed lakes with additional supplies of water exhibited delayed peaks in their seasonal cycles relative to those of non-glacier-fed lakes. Large-scale atmospheric circulation systems, such as the westerlies, Indian summer monsoon, transition in between, and East Asian summer monsoon, were also found to affect the seasonal cycles of lakes. The results of this study suggest that Sentinel-1 SAR data are a powerful tool that can be used to fill gaps in intra-annual lake observations.

Automatic Extraction of Water Inundation Areas Using Sentinel-1 Data for Large Plain Areas

Accurately quantifying water inundation dynamics in terms of both spatial distributions and temporal variability is essential for water resources management. Currently, the water map is usually derived from synthetic aperture radar (SAR) data with the support of auxiliary datasets, using thresholding methods and followed by morphological operations to further refine the results. However, auxiliary datasets may lose efficacy on large plain areas, whilst the parameters of morphological operations are hard to be decided in different situations. Here, a heuristic and automatic water extraction (HAWE) method is proposed to extract the water map from Sentinel-1 SAR data. In the HAWE, we integrate tile-based thresholding and the active contour model, in which the former provides a convincing initial water map used as a heuristic input, and the latter refines the initial map by using image gradient information. The proposed approach was tested on the Dongting Lake plain (China) by comparing the extracted water map with the reference data derived from the Sentinel-2 dataset. For the two selected test sites, the overall accuracy of water classification is between 94.90% and 97.21% whilst the Kappa coefficient is within the range of 0.89 and 0.94. For the entire study area, the overall accuracy is between 94.32% and 96.7% and the Kappa coefficient ranges from 0.80 to 0.90. The results show that the proposed method is capable of extracting water inundations with satisfying accuracy.

A Sentinel-1 Based Processing Chain for Detection of Cyclonic Flood Impacts

In the future, climate change will induce even more severe hurricanes. Not only should these be better understood, but there is also a necessity to improve the assessment of their impacts. Flooding is one of the most common powerful impacts of these storms. Analyzing the impacts of floods is essential in order to delineate damaged areas and study the economic cost of hurricane-related floods. This paper presents an automated processing chain for Sentinel-1 synthetic aperture radar (SAR) data. This processing chain is based on the S1-Tiling algorithm and the normalized difference ratio (NDR). It is able to download and clip S1 images on Sentinel-2 tiles footprints, perform multi-temporal filtering, and threshold NDR images to produce a mask of flooded areas. Applied to two different study zones, subject to hurricanes and cyclones, this chain is reliable and simple to implement. With the rapid mapping product of EMS Copernicus (Emergency Management Service) as reference, the method confers up to 95% accuracy and a Kappa value of 0.75.

Multi-Method Tracking of Monsoon Floods Using Sentinel-1 Imagery

Floods cause great losses in terms of human life and damages to settlements. Since the exposure is a proxy of the risk, it is essential to track flood evolution. The increasing availability of Synthetic Aperture Radar (SAR) imagery extends flood tracking capabilities because of its all-water and day/night acquisition. In this paper, in order to contribute to a better evaluation of the potential of Sentinel-1 SAR imagery to track floods, we analyzed a multi-pulse flood caused by a typhoon in the Camarines Sur Province of Philippines between the end of 2018 and the beginning of 2019. Multiple simple classification methods were used to track the spatial and temporal evolution of the flooded area. Our analysis indicates that Valley Emphasis based manual threshold identification, Otsu methodology, and K-Means Clustering have the potential to be used for tracking large and long-lasting floods, providing similar results. Because of its simplicity, the K-Means Clustering algorithm has the potential to be used in fully automated operational flood monitoring, also because of its good performance in terms of computation time.

Automatic Inundation Mapping Using Sentinel-2 Data Applicable to Both Camargue and Doñana Biosphere Reserves

Flooding periodicity is crucial for biomass production and ecosystem functions in wetland areas. Local monitoring networks may be enriched by spaceborne derived products with a temporal resolution of a few days. Unsupervised computer vision techniques are preferred, since human interference and the use of training data may be kept to a minimum. Recently, a novel automatic local thresholding unsupervised methodology for separating inundated areas from non-inundated ones led to successful results for the Doñana Biosphere Reserve. This study examines the applicability of this approach to Camarque Biosphere Reserve, and proposes alternatives to the original approach to enhance accuracy and applicability for both Camargue and Doñana wetlands in a scientific quest for methods that may serve accurately biomes at both protected areas. In particular, it examines alternative inputs for automatically estimating thresholds while applying various algorithms for estimating the splitting thresholds. Reference maps for Camargue are provided by local authorities, and generated using Sentinel-2 Band 8A (NIR) and Band 12 (SWIR-2). The alternative approaches examined led to high inundation mapping accuracy. In particular, for the Camargue study area and 39 different dates, the alternative approach with the highest overall Kappa coefficient is 0.84, while, for the Doñana Biosphere Reserve and Doñana marshland (a subset of Doñana Reserve) and 7 different dates, is 0.85 and 0.94, respectively. Moreover, there are alternative approaches with high overall Kappa for all areas, i.e., at 0.79 for Camargue, over 0.91 for Doñana marshland, and over 0.82 for Doñana Reserve. Additionally, this study identifies the alternative approaches that perform better when the study area is extensively covered by temporary flooded and emergent vegetation areas (i.e., Camargue Reserve and Doñana marshland) or when it contains a large percentage of dry areas (i.e., Doñana Reserve). The development of credible automatic thresholding techniques that can be applied to different wetlands could lead to a higher degree of automation for map production, while enhancing service utilization by non-trained personnel.

Impacts of Radiometric Uncertainty and Weather-Related Surface Conditions on Soil Moisture Retrievals with Sentinel-1

The radiometric uncertainty of Synthetic Aperture Radar (SAR) observations and weather-related surface conditions caused by frozen conditions, snow and intercepted rain affect the backscatter ( σ 0 ) observations and limit the accuracy of soil moisture retrievals. This study estimates Sentinel-1’s radiometric uncertainty, identifies the effects of weather-related surface conditions on σ 0 and investigates their impact on soil moisture retrievals for various conditions regarding soil moisture, surface roughness and incidence angle. Masking rules for the surface conditions that disturb σ 0 were developed based on meteorological measurements and timeseries of Sentinel-1 observations collected over five forests, five meadows and five cultivated fields in the eastern part of the Netherlands. The Sentinel-1 σ 0 observations appear to be affected by frozen conditions below an air temperature of 1 ∘ C , snow during Sentinel-1’s morning overpasses on meadows and cultivated fields and interception after more than 1.8 m m of rain in the 12 h preceding a Sentinel-1 overpass, whereas dew was not found to be of influence. After the application of these masking rules, the radiometric uncertainty was estimated by the standard deviation of the seasonal anomalies timeseries of the Sentinel-1 forest σ 0 observations. By spatially averaging the σ 0 observations, the Sentinel-1 radiometric uncertainty improves from 0.85 dB for a surface area of 0.25 h a to 0.30 dB for 10 h a for the VV polarization and from 0.89 dB to 0.36 dB for the VH polarization, following approximately an inverse square root dependency on the surface area over which the σ 0 observations are averaged. Deviations in σ 0 were combined with the σ 0 sensitivity to soil moisture as simulated with the Integral Equation Method (IEM) surface scattering model, which demonstrated that both the disturbing effects by the weather-related surface conditions (if not masked) and radiometric uncertainty have a significant impact on the soil moisture retrievals from Sentinel-1. The soil moisture retrieval uncertainty due to radiometric uncertainty ranges from 0.01 m 3 m − 3 up to 0.17 m 3 m − 3 for wet soils and small surface areas. The impacts on soil moisture retrievals are found to be weakly dependent on the surface roughness and the incidence angle, and strongly dependent on the surface area (or the σ 0 disturbance caused by a weather-related surface condition for a specific land cover type) and the soil moisture itself.

Integrating LiDAR data and multi-temporal aerial imagery to map wetland inundation dynamics using Google Earth Engine

The Prairie Pothole Region of North America is characterized by millions of depressional wetlands, which provide critical habitats for globally significant populations of migratory waterfowl and other wildlife species. Due to their relatively small size and shallow depth, these wetlands are highly sensitive to climate variability and anthropogenic changes, exhibiting inter- and intra-annual inundation dynamics. Moderate-resolution satellite imagery (e.g., Landsat, Sentinel) alone cannot be used to effectively delineate these small depressional wetlands. By integrating fine spatial resolution Light Detection and Ranging (LiDAR) data and multi-temporal (2009-2017) aerial images, we developed a fully automated approach to delineate wetland inundation extent at watershed scales using Google Earth Engine. Machine learning algorithms were used to classify aerial imagery with additional spectral indices to extract potential wetland inundation areas, which were further refined using LiDAR-derived landform depressions. The wetland delineation results were then compared to the U.S. Fish and Wildlife Service National Wetlands Inventory (NWI) geospatial dataset and existing global-scale surface water products to evaluate the performance of the proposed method. We tested the workflow on 26 watersheds with a total area of 16,576 km2 in the Prairie Pothole Region. The results showed that the proposed method can not only delineate current wetland inundation status but also demonstrate wetland hydrological dynamics, such as wetland coalescence through fill-spill hydrological processes. Our automated algorithm provides a practical, reproducible, and scalable framework, which can be easily adapted to delineate wetland inundation dynamics at broad geographic scales.

Hydrologic Lag Effects on Wetland Greenhouse Gas Fluxes

Hydrologic margins of wetlands are narrow, transient zones between inundated and dry areas. As water levels fluctuate, the dynamic hydrology at margins may impact wetland greenhouse gas (GHG) fluxes that are sensitive to soil saturation. The Prairie Pothole Region of North America consists of millions of seasonally-ponded wetlands that are ideal for studying hydrologic transition states. Using a long-term GHG database with biweekly flux measurements from 88 seasonal wetlands, we categorized each sample event into wet to wet (W→W), dry to wet (D→W), dry to dry (D→D), or wet to dry (W→D) hydrologic states based on the presence or absence of ponded water from the previous and current event. Fluxes of methane were 5-times lower in the D→W compared to W→W states, indicating a lag ‘ramp-up’ period following ponding. Nitrous oxide fluxes were highest in the W→D state and accounted for 20% of total emissions despite accounting for only 5.2% of wetland surface area during the growing season. Fluxes of carbon dioxide were unaffected by transitions, indicating a rapid acclimation to current conditions by respiring organisms. Results of this study highlight how seasonal drying and re-wetting impact GHGs and demonstrate the importance of hydrologic transitions on total wetland GHG balance.

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.

Preface: Remote Sensing for Flood Mapping and Monitoring of Flood Dynamics

This Special Issue is a collection of papers that focus on the use of remote sensing data and describe methods for flood monitoring and mapping. These articles span a wide range of topics; present novel processing techniques and review methods; and discuss limitations and challenges. This preface provides a brief overview of the content.

Automatic Detection of Open and Vegetated Water Bodies Using Sentinel 1 to Map African Malaria Vector Mosquito Breeding Habitats

Providing timely and accurate maps of surface water is valuable for mapping malaria risk and targeting disease control interventions. Radar satellite remote sensing has the potential to provide this information but current approaches are not suitable for mapping African malarial mosquito aquatic habitats that tend to be highly dynamic, often with emergent vegetation. We present a novel approach for mapping both open and vegetated water bodies using serial Sentinel-1 imagery for Western Zambia. This region is dominated by the seasonally inundated Upper Zambezi floodplain that suffers from a number of public health challenges. The approach uses open source segmentation and machine learning (extra trees classifier), applied to training data that are automatically derived using freely available ancillary data. Refinement is implemented through a consensus approach and Otsu thresholding to eliminate false positives due to dry flat sandy areas. The results indicate a high degree of accuracy (mean overall accuracy 92% st dev 3.6) providing a tractable solution for operationally mapping water bodies in similar large river floodplain unforested environments. For the period studied, 70% of the total water extent mapped was attributed to vegetated water, highlighting the importance of mapping both open and vegetated water bodies for surface water mapping.

Improved Automated Detection of Subpixel-Scale Inundation—Revised Dynamic Surface Water Extent (DSWE) Partial Surface Water Tests

In order to produce useful hydrologic and aquatic habitat data from the Landsat system, the U.S. Geological Survey has developed the “Dynamic Surface Water Extent” (DSWE) Landsat Science Product. DSWE will provide long-term, high-temporal resolution data on variations in inundation extent. The model used to generate DSWE is composed of five decision-rule based tests that do not require scene-based training. To allow its general application, required inputs are limited to the Landsat at-surface reflectance product and a digital elevation model. Unlike other Landsat-based water products, DSWE includes pixels that are only partially covered by water to increase inundation dynamics information content. Previously published DSWE model development included one wetland-focused test developed through visual inspection of field-collected Everglades spectra. A comparison of that test’s output against Everglades Depth Estimation Network (EDEN) in situ data confirmed the expectation that omission errors were a prime source of inaccuracy in vegetated environments. Further evaluation exposed a tendency toward commission error in coniferous forests. Improvements to the subpixel level “partial surface water” (PSW) component of DSWE was the focus of this research. Spectral mixture models were created from a variety of laboratory and image-derived endmembers. Based on the mixture modeling, a more “aggressive” PSW rule improved accuracy in herbaceous wetlands and reduced errors of commission elsewhere, while a second “conservative” test provides an alternative when commission errors must be minimized. Replication of the EDEN-based experiments using the revised PSW tests yielded a statistically significant increase in mean overall agreement (4%, p = 0.01, n = 50) and a statistically significant decrease (11%, p = 0.009, n = 50) in mean errors of omission. Because the developed spectral mixture models included image-derived vegetation endmembers and laboratory spectra for soil groups found across the US, simulations suggest where the revised DSWE PSW tests perform as they do in the Everglades and where they may prove problematic. Visual comparison of DSWE outputs with an unusual variety of coincidently collected images for locations spread throughout the US support conclusions drawn from Everglades quantitative analyses and highlight DSWE PSW component strengths and weaknesses.