Tracking Dynamic Northern Surface Water Changes with High-Frequency Planet CubeSat Imagery

Quantifying the scale of erosion along major coastal aquifers of Pakistan using geospatial and machine learning approaches

Insufficient freshwater recharge and climate change resulted in seawater intrusion in most of the coastal aquifers in Pakistan. Coastal aquifers represent diverse landcover types with varying spectral properties, making it challenging to extract information about their state hence, such investigation requires a combination of geospatial tools. This study aims to monitor erosion along the major coastal aquifers of Pakistan and propose an approach that combines data fusion into the machine and deep learning image segmentation architectures for the erosion and accretion assessment in seascapes. The analysis demonstrated the image segmentation U-Net with EfficientNet backbone achieved the highest F1 score of 0.93, while ResNet101 achieved the lowest F1 score of 0.77. Resultant erosion maps indicated that Sandspit experiencing erosion at 3.14 km2 area. Indus delta is showing erosion, approximately 143 km2 of land over the past 30 years. Sonmiani has undergone substantial erosion with 52.2 km2 land. Miani Hor has experienced erosion up to 298 km2, Bhuri creek has eroded over 4.11 km2, east Phitii creek over 3.30 km2, and Waddi creek over 3.082 km2 land. Tummi creek demonstrates erosion, at 7.12 km2 of land, and East Khalri creek near Keti Bandar has undergone a measured loss of 5.2 km2 land linked with quantified reduction in the vertical sediment flow from 50 (billion cubic meters) to 10 BCM. Our analysis suggests that intense erosions are primarily a result of reduced sediment flow and climate change. Addressing this issue needs to be prioritized coastal management and climate change mitigation framework in Pakistan to safeguard communities. Leveraging emerging solutions, such as loss and damage financing and the integration of nature-based solutions (NbS), should be prioritized for the revival of the coastal aquifers.

Detecting Streamflow in Dryland Rivers Using CubeSats

Determining the flow regime of non-perennial rivers is critical in hydrology. In this study, we developed a new approach using CubeSat imagery to detect streamflow presence using differences in surface reflectance for areas within and outside of a river reach. We calibrated the approach with streamflow records in the Hassayampa River of Arizona over 3 years (2019-2021), finding good agreement in the annual fractions of flowing days at stream gages (R ² = 0.82, p < 0.0001). Subsequently, annual fractions of flowing days were derived at 90 m intervals along the Hassayampa River, finding that 12% of reaches were classified as intermittent, with the remaining as ephemeral. Using a Hovmöller diagram, streamflow presence was visualized in unprecedented spatiotemporal detail, allowing estimates of daily fraction of flowing channel and annual fractions of flowing days. This new tool opens avenues for detecting streamflow and studying hydrological and biogeochemical processes dependent on water presence in drylands.

High-Resolution Snow-Covered Area Mapping in Forested Mountain Ecosystems Using PlanetScope Imagery

Improving high-resolution (meter-scale) mapping of snow-covered areas in complex and forested terrains is critical to understanding the responses of species and water systems to climate change. Commercial high-resolution imagery from Planet Labs, Inc. (Planet, San Francisco, CA, USA) can be used in environmental science, as it has both high spatial (0.7–3.0 m) and temporal (1–2 day) resolution. Deriving snow-covered areas from Planet imagery using traditional radiometric techniques have limitations due to the lack of a shortwave infrared band that is needed to fully exploit the difference in reflectance to discriminate between snow and clouds. However, recent work demonstrated that snow cover area (SCA) can be successfully mapped using only the PlanetScope 4-band (Red, Green, Blue and NIR) reflectance products and a machine learning (ML) approach based on convolutional neural networks (CNN). To evaluate how additional features improve the existing model performance, we: (1) build on previous work to augment a CNN model with additional input data including vegetation metrics (Normalized Difference Vegetation Index) and DEM-derived metrics (elevation, slope and aspect) to improve SCA mapping in forested and open terrain, (2) evaluate the model performance at two geographically diverse sites (Gunnison, Colorado, USA and Engadin, Switzerland), and (3) evaluate the model performance over different land-cover types. The best augmented model used the Normalized Difference Vegetation Index (NDVI) along with visible (red, green, and blue) and NIR bands, with an F-score of 0.89 (Gunnison) and 0.93 (Engadin) and was found to be 4% and 2% better than when using canopy height- and terrain-derived measures at Gunnison, respectively. The NDVI-based model improves not only upon the original band-only model’s ability to detect snow in forests, but also across other various land-cover types (gaps and canopy edges). We examined the model’s performance in forested areas using three forest canopy quantification metrics and found that augmented models can better identify snow in canopy edges and open areas but still underpredict snow cover under forest canopies. While the new features improve model performance over band-only options, the models still have challenges identifying the snow under trees in dense forests, with performance varying as a function of the geographic area. The improved high-resolution snow maps in forested environments can support studies involving climate change effects on mountain ecosystems and evaluations of hydrological impacts in snow-dominated river basins.

Retrieval of Remotely Sensed Sediment Grain Size Evolution Characteristics along the Southwest Coast of Laizhou Bay Based on Support Vector Machine Learning

Grain size is the basic property of intertidal zone sediment. Grain size acts as an indicator of sedimentary processes and geomorphological evolution under human and nature interactions. The remote sensing technique provides an alternative for sediment grain-size parameter monitoring with the advantages of wide coverage and real-time surveying. This paper attempted to map the distributions of three sediment grain size contents and the mean grain size with multitemporal Landsat images along the southwestern coast of Laizhou Bay, China, from 1989 to 2015. Considering the low correlations between the measured reflectance and grain-size parameters, we used a support vector machine (SVM) to develop a nonlinear calibration model by taking several band indices as input variables. Then, the performance of the back propagation neural network (BPNN) was determined and discussed with that of the SVM. The SVM performed better than the BPNN in calibrating the four grain-size parameters based on a comparison of R2 and the root-mean-square error (RMSE). Moreover, an atmospheric correction algorithm originally proposed for case II water enabled the TM\ETM+ images to be precisely atmospherically corrected in this study. The SVM-mapped spatial-temporal grain-size variation showed a coarsening trend, which agreed with that obtained during in situ measurements in a former study. The changes in Yellow River discharge and precipitation associated with the coarsening trend were further analyzed. The yielded results showed that the coarsening trend and reduction in tidal flat area might be aggravated with overutilization. More reasonable planning would be necessary in this case.

Can Remote Sensing Fill the United States’ Monitoring Gap for Watershed Management?

Remote sensing has been heralded as the silver bullet in water quality modeling and watershed management, and yet a quantitative mapping of where its applicability is likely and most useful has not been undertaken so far. Here, we combine geospatial models of cloud cover as a proxy for the likelihood of acquiring remote scenes and the shortest time of travel to population centers as a proxy for accessibility to ground-truth remote sensing data for water quality monitoring and produce maps of the potential of remote sensing in watershed management in the United States. We generate several maps with different cost-payoff relationships to help stakeholders plan and incentivize remote sensing-based monitoring campaigns. Additionally, we combine these remote sensing potential maps with spatial indices of population, water demand, ecosystem services, pollution risk, and monitoring coverage deficits to identify where remote sensing likely has the greatest role to play. We find that the Southwestern United States and the Central plains regions are generally suitable for remote sensing for watershed management even under the most stringent costing projections, but that the potential for using remote sensing can extend further North and East as constraints are relaxed. We also find large areas in the Southern United States and sporadic watersheds in the Northeast and Northwest seaboards and the Midwest would likely benefit most from using remote sensing for watershed monitoring. Although developed herein for watershed decision support in the United States, our approach is readily generalizable to other environmental domains and across the world.

Remote Sensing of Surface Water Dynamics in the Context of Global Change—A Review

Inland surface water is often the most accessible freshwater source. As opposed to groundwater, surface water is replenished in a comparatively quick cycle, which makes this vital resource—if not overexploited—sustainable. From a global perspective, freshwater is plentiful. Still, depending on the region, surface water availability is severely limited. Additionally, climate change and human interventions act as large-scale drivers and cause dramatic changes in established surface water dynamics. Actions have to be taken to secure sustainable water availability and usage. This requires informed decision making based on reliable environmental data. Monitoring inland surface water dynamics is therefore more important than ever. Remote sensing is able to delineate surface water in a number of ways by using optical as well as active and passive microwave sensors. In this review, we look at the proceedings within this discipline by reviewing 233 scientific works. We provide an extensive overview of used sensors, the spatial and temporal resolution of studies, their thematic foci, and their spatial distribution. We observe that a wide array of available sensors and datasets, along with increasing computing capacities, have shaped the field over the last years. Multiple global analysis-ready products are available for investigating surface water area dynamics, but so far none offer high spatial and temporal resolution.

Impact of the Dates of Input Image Pairs on Spatio-Temporal Fusion for Time Series with Different Temporal Variation Patterns

Dense time series of remote sensing images with high spatio-temporal resolution are critical for monitoring land surface dynamics in heterogeneous landscapes. Spatio-temporal fusion is an effective solution to obtaining such time series images. Many spatio-temporal fusion methods have been developed for producing high spatial resolution images at frequent intervals by blending fine spatial images and coarse spatial resolution images. Previous studies have revealed that the accuracy of fused images depends not only on the fusion algorithm, but also on the input image pairs being used. However, the impact of input images dates on the fusion accuracy for time series with different temporal variation patterns remains unknown. In this paper, the impact of input image pairs on the fusion accuracy for monotonic linear change (MLC), monotonic non-linear change (MNLC), and non-monotonic change (NMC) time periods were evaluated, respectively, and the optimal selection strategies of input image dates for different situations were proposed. The 16-day composited NDVI time series (i.e., Collection 6 MODIS NDVI product) were used to present the temporal variation patterns of land surfaces in the study areas. To obtain sufficient observation dates to evaluate the impact of input image pairs on the spatio-temporal fusion accuracy, we utilized the Harmonized Landsat-8 Sentinel-2 (HLS) data. The ESTARFM was selected as the spatio-temporal fusion method for this study. The results show that the impact of input image date on the accuracy of spatio-temporal fusion varies with the temporal variation patterns of the time periods being fused. For the MLC period, the fusion accuracy at the prediction date (PD) is linearly correlated to the time interval between the change date (CD) of the input image and the PD, but the impact of the input image date on the fusion accuracy at the PD is not very significant. For the MNLC period, the fusion accuracy at the PD is non-linearly correlated to the time interval between the CD and the PD, the impact of the time interval between the CD and the PD on the fusion accuracy is more significant for the MNLC than for the MLC periods. Given the similar change of time intervals between the CD and the PD, the increments of R2 of fusion result for the MNLC is over ten times larger than those for the MLC. For the NMC period, a shorter time interval between the CD and the PD does not lead to higher fusion accuracies. On the contrary, it may lower the fusion accuracy. This study suggests that temporal variation patterns of the data must be taken into account when selecting optimal dates of input images in the fusion model.

Change Detection of Amazonian Alluvial Gold Mining Using Deep Learning and Sentinel-2 Imagery

Monitoring changes within the land surface and open water bodies is critical for natural resource management, conservation, and environmental policy. While the use of satellite imagery for these purposes is common, fine-scale change detection can be a technical challenge. Difficulties arise from variable atmospheric conditions and the problem of assigning pixels to individual objects. We examined the degree to which two machine learning approaches can better characterize change detection in the context of a current conservation challenge, artisanal small-scale gold mining (ASGM). We obtained Sentinel-2 imagery and consulted with domain experts to construct an open-source labeled land-cover change dataset. The focus of this dataset is the Madre de Dios (MDD) region in Peru, a hotspot of ASGM activity. We also generated datasets of active ASGM areas in other countries (Venezuela, Indonesia, and Myanmar) for out-of-sample testing. With these labeled data, we utilized a supervised (E-ReCNN) and semi-supervised (SVM-STV) approach to study binary and multi-class change within mining ponds in the MDD region. Additionally, we tested how the inclusion of multiple channels, histogram matching, and La*b* color metrics improved the performance of the models and reduced the influence of atmospheric effects. Empirical results show that the supervised E-ReCNN method on 6-Channel histogram-matched images generated the most accurate detection of change not only in the focal region (Kappa: 0.92 (± 0.04), Jaccard: 0.88 (± 0.07), F1: 0.88 (± 0.05)) but also in the out-of-sample prediction regions (Kappa: 0.90 (± 0.03), Jaccard: 0.84 (± 0.04), and F1: 0.77 (± 0.04)). While semi-supervised methods did not perform as accurately on 6- or 10-channel imagery, histogram matching and the inclusion of La*b* metrics generated accurate results with low memory and resource costs. These results show that E-ReCNN is capable of accurately detecting specific and object-oriented environmental changes related to ASGM. E-ReCNN is scalable to areas outside the focal area and is a method of change detection that can be extended to other forms of land-use modification.

Multi-Temporal Surface Water Classification for Four Major Rivers from the Peruvian Amazon

We describe a new minimum extent, persistent surface water classification for reaches of four major rivers in the Peruvian Amazon (i.e., Amazon, Napo, Pastaza, Ucayali). These data were generated by the Peruvian Amazon Rural Livelihoods and Poverty (PARLAP) Project which aims to better understand the nexus between livelihoods (e.g., fishing, agriculture, forest use, trade), poverty, and conservation in the Peruvian Amazon over a 35,000 km river network. Previous surface water datasets do not adequately capture the temporal changes in the course of the rivers, nor discriminate between primary main channel and non-main channel (e.g., oxbow lakes) water. We generated the surface water classifications in Google Earth Engine from Landsat TM 5, 7 ETM+, and 8 OLI satellite imagery for time periods from circa 1989, 2000, and 2015 using a hierarchical logical binary classification predominantly based on a modified Normalized Difference Water Index (mNDWI) and shortwave infrared surface reflectance. We included surface reflectance in the blue band and brightness temperature to minimize misclassification. High accuracies were achieved for all time periods (>90%).

Estimation of river flow using CubeSats remote sensing

River flow characterizes the integrated response from watersheds, so it is essential to quantify to understand the changing water cycle and underpin the sustainable management of freshwaters. However, river gauging stations are in decline with ground-based observation networks shrinking. This study proposes a novel approach of estimating river flows using the Planet CubeSats constellation with the possibility to monitor on a daily basis at the sub-catchment scale through remote sensing. The methodology relates the river discharge to the water area that is extracted from the satellite image analysis. As a testbed, a series of Surface Reflectance PlanetScope images and observed streamflow data in Araguaia River (Brazil) were selected to develop and validate the methodology. The study involved the following steps: (1) survey of measurements of water level and river discharge using in-situ data from gauge-based Conventional Station (CS) and measurements of altimetry using remote data from JASON-2 Virtual Station (JVS); (2) survey of Planet CubeSat images for dates in step 1 and without cloud cover; (3) image preparation including clipping based on different buffer areas and calculation of the Normalized Difference Vegetation Index (NDVI) per image; (4) water bodies areas calculation inside buffers in the Planet CubeSat images; and (5) correlation analysis of CubeSat water bodies areas with JVS and CS data. Significant correlations between the water bodies areas with JVS (R² = 88.83%) and CS (R² = 96.49%) were found, indicating that CubeSat images can be used as a CubeSat Virtual Station (CVS) to estimate the river flow. This newly proposed methodology using CubeSats allows for more accurate results than the JVS-based method used by the Brazilian National Water Agency (ANA) at present. Moreover, CVS requires small areas of remote sensing data to estimate with high accuracy the river flow and the height variation of the water in different timeframes. This method can be used to monitor sub-basin scale discharge and to improve water management, particularly in developing countries where the presence of conventional stations is often very limited.

Water quality challenges associated with industrial logging of a karst landscape: Guadalcanal, Solomon Islands

Human disturbance of karst landscapes in tropical volcanic islands present a unique challenge for understanding sediment transport to the coastal zone. Here we present the first evidence of urban drinking water quality impacts from industrial logging in the Solomon Islands. Despite only 6% of the Honiara's drinking water catchment being disturbed by logging, rhodamine dye tracers demonstrated complex karst sinkholes that led to high suspended sediment concentrations being transported from neighbouring Kovi catchment into the Kongulai water supply offtake point for Honiara. This has resulted in the exceedance of practical treatment thresholds of 20 NTU 9.5% of the time, leading to water supply for the majority of Honiara's residents being unavailable for 58 days in 2019. This work highlights the cost-benefit disparity between industrial logging yielding minimal short-term economic yields in comparison to on-going broader impacts of increased coastal sediment transport while restricting water supply to a developing nation's capital.

Monitoring the Transformation of Arctic Landscapes: Automated Shoreline Change Detection of Lakes Using Very High Resolution Imagery

Water bodies are a highly abundant feature of Arctic permafrost ecosystems and strongly influence their hydrology, ecology and biogeochemical cycling. While very high resolution satellite images enable detailed mapping of these water bodies, the increasing availability and abundance of this imagery calls for fast, reliable and automatized monitoring. This technical work presents a largely automated and scalable workflow that removes image noise, detects water bodies, removes potential misclassifications from infrastructural features, derives lake shoreline geometries and retrieves their movement rate and direction on the basis of ortho-ready very high resolution satellite imagery from Arctic permafrost lowlands. We applied this workflow to typical Arctic lake areas on the Alaska North Slope and achieved a successful and fast detection of water bodies. We derived representative values for shoreline movement rates ranging from 0.40–0.56 m yr−1 for lake sizes of 0.10 ha–23.04 ha. The approach also gives an insight into seasonal water level changes. Based on an extensive quantification of error sources, we discuss how the results of the automated workflow can be further enhanced by incorporating additional information on weather conditions and image metadata and by improving the input database. The workflow is suitable for the seasonal to annual monitoring of lake changes on a sub-meter scale in the study areas in northern Alaska and can readily be scaled for application across larger regions within certain accuracy limitations.

Fine-Resolution Mapping of Pan-Arctic Lake Ice-Off Phenology Based on Dense Sentinel-2 Time Series Data

The timing of lake ice-off regulates biotic and abiotic processes in Arctic ecosystems. Due to the coarse spatial and temporal resolution of available satellite data, previous studies mainly focused on lake-scale investigations of melting/freezing, hindering the detection of subtle patterns within heterogeneous landscapes. To fill this knowledge gap, we developed a new approach for fine-resolution mapping of Pan-Arctic lake ice-off phenology. Using the Scene Classification Layer data derived from dense Sentinel-2 time series images, we estimated the pixel-by-pixel ice break-up end date information by seeking the transition time point when the pixel is completely free of ice. Applying this approach on the Google Earth Engine platform, we mapped the spatial distribution of the break-up end date for 45,532 lakes across the entire Arctic (except for Greenland) for the year 2019. The evaluation results suggested that our estimations matched well with both in situ measurements and an existing lake ice phenology product. Based on the generated map, we estimated that the average break-up end time of Pan-Arctic lakes is 172 ± 13.4 (measured in day of year) for the year 2019. The mapped lake ice-off phenology exhibits a latitudinal gradient, with a linear slope of 1.02 days per degree from 55°N onward. We also demonstrated the importance of lake and landscape characteristics in affecting spring lake ice melting. The proposed approach offers new possibilities for monitoring the seasonal Arctic lake ice freeze–thaw cycle, benefiting the ongoing efforts of combating and adapting to climate change.

Satellite Flood Inundation Assessment and Forecast Using SMAP and Landsat

The capability and synergistic use of multisource satellite observations for flood monitoring and forecasts is crucial for improving disaster preparedness and mitigation. Here, surface fractional water cover (FW) retrievals derived from Soil Moisture Active Passive (SMAP) L-band (1.4 GHz) brightness temperatures were used for flood assessment over southeast Africa during the Cyclone Idai event. We then focused on five subcatchments of the Pungwe basin and developed a machine learning based approach with the support of Google Earth Engine for daily (24-h) forecasting of FW and 30-m inundation downscaling and mapping. The Classification and Regression Trees model was selected and trained using retrievals derived from SMAP and Landsat coupled with rainfall forecasts from the NOAA Global Forecast System. Independent validation showed that FW predictions over randomly selected dates are highly correlated (R = 0.87) with the Landsat observations. The forecast results captured the flood temporal dynamics from the Idai event; and the associated 30-m downscaling results showed inundation spatial patterns consistent with independent satellite synthetic aperture radar observations. The data-driven approach provides new capacity for flood monitoring and forecasts leveraging synergistic satellite observations and big data analysis, which is particularly valuable for data sparse regions.

Measuring Alpha and Beta Diversity by Field and Remote-Sensing Data: A Challenge for Coastal Dunes Biodiversity Monitoring

Combining field collected and remotely sensed (RS) data represents one of the most promising approaches for an extensive and up-to-date ecosystem assessment. We investigated the potential of the so called spectral variability hypothesis (SVH) in linking field-collected and remote-sensed data in Mediterranean coastal dunes and explored if spectral diversity provides reliable information to monitor floristic diversity, as well as the consistency of such information in altered ecosystems due to plant invasions. We analyzed alpha diversity and beta diversity, integrating floristic field and Remote-Sensing PlanetScope data in the Tyrrhenian coast (Central Italy). We explored the relationship among alpha field diversity (species richness, Shannon index, inverse Simpson index) and spectral variability (distance from the spectral centroid index) through linear regressions. For beta diversity, we implemented a distance decay model (DDM) relating field pairwise (Jaccard similarities index, Bray–Curtis similarities index) and spectral pairwise (Euclidean distance) measures. We observed a positive relationship between alpha diversity and spectral heterogeneity with richness reporting the higher R score. As for DDM, we found a significant relationship between Bray–Curtis floristic similarity and Euclidean spectral distance. We provided a first assessment of the relationship between floristic and spectral RS diversity in Mediterranean coastal dune habitats (i.e., natural or invaded). SVH provided evidence about the potential of RS for estimating diversity in complex and dynamic landscapes.

Declining greenness in Arctic-boreal lakes

Arctic and boreal regions are undergoing dramatic warming and also possess the world’s highest concentration of lakes. However, ecological changes in lakes are poorly understood. We present a continental-scale trend analysis of satellite lake color in the green wavelengths, which shows declining greenness from 1984 to 2019 in Arctic-boreal lakes across western North America. Annual 30-m Landsat composites indicate lake greenness has decreased by 15%. Our findings show a relationship between lake color, rising air temperatures, and increasing precipitation, supporting the theory that warming may be increasing connectivity between lakes and surrounding landscapes. Overall, our results bring a powerful set of observations in support of the hypothesis that lakes are sentinels for global change in rapidly warming Arctic-boreal ecosystems.

The highest concentration of the world’s lakes are found in Arctic-boreal regions [C. Verpoorter, T. Kutser, D. A. Seekell, L. J. Tranvik, Geophys. Res. Lett. 41, 6396–6402 (2014)], and consequently are undergoing the most rapid warming [J. E. Overland et al., Arctic Report Card (2018)]. However, the ecological response of Arctic-boreal lakes to warming remains highly uncertain. Historical trends in lake color from remote sensing observations can provide insights into changing lake ecology, yet have not been examined at the pan-Arctic scale. Here, we analyze time series of 30-m Landsat growing season composites to quantify trends in lake greenness for >4 × 10⁵ waterbodies in boreal and Arctic western North America. We find lake greenness declined overall by 15% from the first to the last decade of analysis within the 6.3 × 10⁶-km² study region but with significant spatial variability. Greening declines were more likely to be found in areas also undergoing increases in air temperature and precipitation. These findings support the hypothesis that warming has increased connectivity between lakes and the land surface [A. Bring et al., J. Geophys. Res. Biogeosciences 121, 621–649 (2016)], with implications for lake carbon cycling and energy budgets. Our study provides spatially explicit information linking climate to pan-Arctic lake color changes, a finding that will help target future ecological monitoring in remote yet rapidly changing regions.

Towards global flood mapping onboard low cost satellites with machine learning

Spaceborne Earth observation is a key technology for flood response, offering valuable information to decision makers on the ground. Very large constellations of small, nano satellites- 'CubeSats' are a promising solution to reduce revisit time in disaster areas from days to hours. However, data transmission to ground receivers is limited by constraints on power and bandwidth of CubeSats. Onboard processing offers a solution to decrease the amount of data to transmit by reducing large sensor images to smaller data products. The ESA's recent PhiSat-1 mission aims to facilitate the demonstration of this concept, providing the hardware capability to perform onboard processing by including a power-constrained machine learning accelerator and the software to run custom applications. This work demonstrates a flood segmentation algorithm that produces flood masks to be transmitted instead of the raw images, while running efficiently on the accelerator aboard the PhiSat-1. Our models are trained on WorldFloods: a newly compiled dataset of 119 globally verified flooding events from disaster response organizations, which we make available in a common format. We test the system on independent locations, demonstrating that it produces fast and accurate segmentation masks on the hardware accelerator, acting as a proof of concept for this approach.

Monitoring Variations in Lake Water Storage with Satellite Imagery and Citizen Science

Despite lakes being a key part of the global water cycle and a crucial water resource, there is limited understanding of whether regional or lake-specific factors control water storage variations in small lakes. Here, we study groups of small, unregulated lakes in North Carolina, Washington, Illinois, and Wisconsin, USA using lake level measurements gathered by citizen scientists and lake surface area measurements from optical satellite imagery. We show the lake level measurements to be highly accurate when compared to automated gauges (mean absolute error = 1.6 cm). We compare variations in lake water storage between pairs of lakes within these four states. On average, water storage variations in lake pairs across all study regions are moderately positively correlated (ρ = 0.49) with substantial spread in the degree of correlation. The distance between lake pairs and the extent to which their changes in volume are correlated show a weak but statistically significant negative relationship. Our results indicate that, on regional scales, distance is not a primary factor governing lake water storage patterns, which suggests that other, perhaps lakes-specific, factors must also play important roles.

Glacial Lakes Mapping Using Multi Satellite PlanetScope Imagery and Deep Learning

Glacial lakes mapping using satellite remote sensing data are important for studying the effects of climate change as well as for the mitigation and risk assessment of a Glacial Lake Outburst Flood (GLOF). The 3U cubesat constellation of Planet Labs offers the capability of imaging the whole Earth landmass everyday at 3–4 m spatial resolution. The higher spatial, as well as temporal resolution of PlanetScope imagery in comparison with Landsat-8 and Sentinel-2, makes it a valuable data source for monitoring the glacial lakes. Therefore, this paper explores the potential of the PlanetScope imagery for glacial lakes mapping with a focus on the Hindu Kush, Karakoram and Himalaya (HKKH) region. Though the revisit time of the PlanetScope imagery is short, courtesy of 130+ small satellites, this imagery contains only four bands and the imaging sensors in these small satellites exhibit varying spectral responses as well as lower dynamic range. Furthermore, the presence of cast shadows in the mountainous regions and varying spectral signature of the water pixels due to differences in composition, turbidity and depth makes it challenging to automatically and reliably extract surface water in PlanetScope imagery. Keeping in view these challenges, this work uses state of the art deep learning models for pixel-wise classification of PlanetScope imagery into the water and background pixels and compares the results with Random Forest and Support Vector Machine classifiers. The deep learning model is based on the popular U-Net architecture. We evaluate U-Net architecture similar to the original U-Net as well as a U-Net with a pre-trained EfficientNet backbone. In order to train the deep neural network, ground truth data are generated by manual digitization of the surface water in PlanetScope imagery with the aid of Very High Resolution Satellite (VHRS) imagery. The created dataset consists of more than 5000 water bodies having an area of approx. 71km2 in eight different sites in the HKKH region. The evaluation of the test data show that the U-Net with EfficientNet backbone achieved the highest F1 Score of 0.936. A visual comparison with the existing glacial lake inventories is then performed over the Baltoro glacier in the Karakoram range. The results show that the deep learning model detected significantly more lakes than the existing inventories, which have been derived from Landsat OLI imagery. The trained model is further evaluated on the time series PlanetScope imagery of two glacial lakes, which have resulted in an outburst flood. The output of the U-Net is also compared with the GLakeMap data. The results show that the higher spatial and temporal resolution of PlanetScope imagery is a significant advantage in the context of glacial lakes mapping and monitoring.

Fishpond Mapping by Spectral and Spatial-Based Filtering on Google Earth Engine: A Case Study in Singra Upazila of Bangladesh

Inland aquaculture in Bangladesh has been growing fast in the last decade. The underlying land use/land cover (LULC) change is an important indicator of socioeconomic and food structure change in Bangladesh, and fishpond mapping is essential to understand such LULC change. Previous research often used water indexes (WI), such as Normalized Difference Water Index (NDWI) and Modified Normalized Difference Water Index (MNDWI), to enhance water bodies and use shape-based metrics to assist classification of individual water features, such as coastal aquaculture ponds. However, inland fishponds in Bangladesh are generally extremely small, and little research has investigated mapping of such small water objects without high-resolution images. Thus, this research aimed to bridge the knowledge gap by developing and evaluating an automatic fishpond mapping workflow with Sentinel-2 images that is implemented on Google Earth Engine (GEE) platform. The workflow mainly includes two steps: (1) the spectral filtering phase that uses a pixel selection technique and an image segmentation method to automatically identify all-year-inundated water bodies and (2) spatial filtering phase to further classify all-year-inundated water bodies into fishponds and non-fishponds using object-based features (OBF). To evaluate the performance of the workflow, we conducted a case study in the Singra Upazila of Bangladesh, and our method can efficiently map inland fishponds with a precision score of 0.788. Our results also show that the pixel selection technique is essential in identifying inland fishponds that are generally small. As the workflow is implemented on GEE, it can be conveniently applied to other regions.

Physics-based Bathymetry and Water Quality Retrieval Using PlanetScope Imagery: Impacts of 2020 COVID-19 Lockdown and 2019 Extreme Flood in the Venice Lagoon

The recent PlanetScope constellation (130+ satellites currently in orbit) has shifted the high spatial resolution imaging into a new era by capturing the Earth’s landmass including inland waters on a daily basis. However, studies on the aquatic-oriented applications of PlanetScope imagery are very sparse, and extensive research is still required to unlock the potentials of this new source of data. As a first fully physics-based investigation, we aim to assess the feasibility of retrieving bathymetric and water quality information from the PlanetScope imagery. The analyses are performed based on Water Color Simulator (WASI) processor in the context of a multitemporal analysis. The WASI-based radiative transfer inversion is adapted to process the PlanetScope imagery dealing with the low spectral resolution and atmospheric artifacts. The bathymetry and total suspended matter (TSM) are mapped in the relatively complex environment of Venice lagoon during two benchmark events: The coronavirus disease 2019 (COVID-19) lockdown and an extreme flood occurred in November 2019. The retrievals of TSM imply a remarkable reduction of the turbidity during the lockdown, due to the COVID-19 pandemic and capture the high values of TSM during the flood condition. The results suggest that sizable atmospheric and sun-glint artifacts should be mitigated through the physics-based inversion using the surface reflectance products of PlanetScope imagery. The physics-based inversion demonstrated high potentials in retrieving both bathymetry and TSM using the PlanetScope imagery.

Timing of Landsat Overpasses Effectively Captures Flow Conditions of Large Rivers

Satellites provide a temporally discontinuous record of hydrological conditions along Earth’s rivers (e.g., river width, height, water quality). The degree to which archived satellite data effectively capture the overall population of river flow frequency is unknown. Here, we use the entire archives of Landsat 5, 7, and 8 to determine when a cloud-free image is available over the United States Geological Survey (USGS) river gauges located on Landsat-observable rivers. We compare the flow frequency distribution derived from the daily gauge record to the flow frequency distribution derived from ideally sampling gauged discharge based on the timing of cloud-free Landsat overpasses. Examining the patterns of flow frequency across multiple gauges, we find that there is not a statistically significant difference between the flow frequency distribution associated with observations contained within the Landsat archive and the flow frequency distribution derived from the daily gauge data (α = 0.05), except for hydrological extremes like maximum and minimum flow. At individual gauges, we find that Landsat observations span a wide range of hydrological conditions (97% of total flow variability observed in 90% of the study gauges) but the degree to which the Landsat sample can represent flow frequency distribution varies from location to location and depends on sample size. The results of this study indicate that the Landsat archive is, on average, representative of the temporal frequencies of hydrological conditions present along Earth’s large rivers with broad utility for hydrological, ecologic and biogeochemical evaluations of river systems.

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.

Remote Sensing of Boreal Wetlands 2: Methods for Evaluating Boreal Wetland Ecosystem State and Drivers of Change

The following review is the second part of a two part series on the use of remotely sensed data for quantifying wetland extent and inferring or measuring condition for monitoring drivers of change on wetland environments. In the first part, we introduce policy makers and non-users of remotely sensed data with an effective feasibility guide on how data can be used. In the current review, we explore the more technical aspects of remotely sensed data processing and analysis using case studies within the literature. Here we describe: (a) current technologies used for wetland assessment and monitoring; (b) the latest algorithmic developments for wetland assessment; (c) new technologies; and (d) a framework for wetland sampling in support of remotely sensed data collection. Results illustrate that high or fine spatial resolution pixels (≤10 m) are critical for identifying wetland boundaries and extent, and wetland class, form and type, but are not required for all wetland sizes. Average accuracies can be up to 11% better (on average) than medium resolution (11–30 m) data pixels when compared with field validation. Wetland size is also a critical factor such that large wetlands may be almost as accurately classified using medium-resolution data (average = 76% accuracy, stdev = 21%). Decision-tree and machine learning algorithms provide the most accurate wetland classification methods currently available, however, these also require sampling of all permutations of variability. Hydroperiod accuracy, which is dependent on instantaneous water extent for single time period datasets does not vary greatly with pixel resolution when compared with field data (average = 87%, 86%) for high and medium resolution pixels, respectively. The results of this review provide users with a guideline for optimal use of remotely sensed data and suggested field methods for boreal and global wetland studies.

Remote Sensing of Boreal Wetlands 1: Data Use for Policy and Management

Wetlands have and continue to undergo rapid environmental and anthropogenic modification and change to their extent, condition, and therefore, ecosystem services. In this first part of a two-part review, we provide decision-makers with an overview on the use of remote sensing technologies for the ‘wise use of wetlands’, following Ramsar Convention protocols. The objectives of this review are to provide: (1) a synthesis of the history of remote sensing of wetlands, (2) a feasibility study to quantify the accuracy of remotely sensed data products when compared with field data based on 286 comparisons found in the literature from 209 articles, (3) recommendations for best approaches based on case studies, and (4) a decision tree to assist users and policymakers at numerous governmental levels and industrial agencies to identify optimal remote sensing approaches based on needs, feasibility, and cost. We argue that in order for remote sensing approaches to be adopted by wetland scientists, land-use managers, and policymakers, there is a need for greater understanding of the use of remote sensing for wetland inventory, condition, and underlying processes at scales relevant for management and policy decisions. The literature review focuses on boreal wetlands primarily from a Canadian perspective, but the results are broadly applicable to policymakers and wetland scientists globally, providing knowledge on how to best incorporate remotely sensed data into their monitoring and measurement procedures. This is the first review quantifying the accuracy and feasibility of remotely sensed data and data combinations needed for monitoring and assessment. These include, baseline classification for wetland inventory, monitoring through time, and prediction of ecosystem processes from individual wetlands to a national scale.

Research Trends in the Use of Remote Sensing for Inland Water Quality Science: Moving Towards Multidisciplinary Applications

Remote sensing approaches to measuring inland water quality date back nearly 50 years to the beginning of the satellite era. Over this time span, hundreds of peer-reviewed publications have demonstrated promising remote sensing models to estimate biological, chemical, and physical properties of inland waterbodies. Until recently, most of these publications focused largely on algorithm development as opposed to implementation of those algorithms to address specific science questions. This slow evolution contrasts with terrestrial and oceanic remote sensing, where methods development in the 1970s led to publications focused on understanding spatially expansive, complex processes as early as the mid-1980s. This review explores the progression of inland water quality remote sensing from methodological development to scientific applications. We use bibliometric analysis to assess overall patterns in the field and subsequently examine 236 key papers to identify trends in research focus and scale. The results highlight an initial 30 year period where the majority of publications focused on model development and validation followed by a spike in publications, beginning in the early-2000s, applying remote sensing models to analyze spatiotemporal trends, drivers, and impacts of changing water quality on ecosystems and human populations. Recent and emerging resources, including improved data availability and enhanced processing platforms, are enabling researchers to address challenging science questions and model spatiotemporally explicit patterns in water quality. Examination of the literature shows that the past 10–15 years has brought about a focal shift within the field, where researchers are using improved computing resources, datasets, and operational remote sensing algorithms to better understand complex inland water systems. Future satellite missions promise to continue these improvements by providing observational continuity with spatial/spectral resolutions ideal for inland waters.

Monitoring Oil Exploitation Infrastructure and Dirt Roads with Object-Based Image Analysis and Random Forest in the Eastern Mongolian Steppe

Information on the spatial distribution of human disturbance is important for assessing and monitoring land degradation. In the Eastern Mongolian Steppe Ecosystem, one of the major driving factors of human-induced land degradation is the expansion of road networks mainly due to intensifications of oil exploration and exploitation. So far, neither the extents of road networks nor the extent of surrounding grasslands affected by the oil industry are monitored which is generally labor consuming. This causes that no information on the changes in the area which is affected by those disturbance drivers is available. Consequently, the study aim is to provide a cost-effective methodology to classify infrastructure and oil exploitation areas from remotely sensed images using object-based classifications with Random Forest. By combining satellite data with different spatial and spectral resolutions (PlanetScope, RapidEye, and Landsat ETM+), the product delivers data since 2005. For the classification variables, segmentation, spectral characteristics, and indices were extracted from all above mentioned imagery and used as predictors. Results show that overall accuracies of land use maps ranged 73%–93% mainly depending on satellites’ spatial resolution. Since 2005, the area of grassland disturbed by dirt roads and oil exploitation infrastructure increased by 88% with its highest expansion by 47% in the period 2005–2010. Settlements and croplands remained relatively constant throughout the 13 years. Comparison of multiscale classification suggests that, although high spatial resolutions are clearly beneficial, all datasets were useful to delineate linear features such as roads. Consequently, the results of this study provide an effective evaluation for the potential of Random Forest for extracting relatively narrow linear features such as roads from multiscale satellite images and map products that are possible to use for detailed land degradation assessments.

A Review of Earth Observation-Based Analyses for Major River Basins

Regardless of political boundaries, river basins are a functional unit of the Earth’s land surface and provide an abundance of resources for the environment and humans. They supply livelihoods supported by the typical characteristics of large river basins, such as the provision of freshwater, irrigation water, and transport opportunities. At the same time, they are impacted i.e., by human-induced environmental changes, boundary conflicts, and upstream–downstream inequalities. In the framework of water resource management, monitoring of river basins is therefore of high importance, in particular for researchers, stake-holders and decision-makers. However, land surface and surface water properties of many major river basins remain largely unmonitored at basin scale. Several inventories exist, yet consistent spatial databases describing the status of major river basins at global scale are lacking. Here, Earth observation (EO) is a potential source of spatial information providing large-scale data on the status of land surface properties. This review provides a comprehensive overview of existing research articles analyzing major river basins primarily using EO. Furthermore, this review proposes to exploit EO data together with relevant open global-scale geodata to establish a database and to enable consistent spatial analyses and evaluate past and current states of major river basins.

Winter Wheat Mapping Based on Sentinel-2 Data in Heterogeneous Planting Conditions

Monitoring and mapping the spatial distribution of winter wheat accurately is important for crop management, damage assessment and yield prediction. In this study, northern and central Anhui province were selected as study areas, and Sentinel-2 imagery was employed to map winter wheat distribution and the results were verified with Planet imagery in the 2017–2018 growing season. The Sentinel-2 imagery at the heading stage was identified as the optimum period for winter wheat area extraction after analyzing the images from different growth stages using the Jeffries–Matusita distance method. Therefore, ten spectral bands, seven vegetation indices (VI), water index and building index generated from the image at the heading stage were used to classify winter wheat areas by a random forest (RF) algorithm. The result showed that the accuracy was from 93% to 97%, with a Kappa above 0.82 and a percentage error lower than 5% in northern Anhui, and an accuracy of about 80% with Kappa ranging from 0.70 to 0.78 and a percentage error of about 20% in central Anhui. Northern Anhui has a large planting scale of winter wheat and flat terrain while central Anhui grows relatively small winter wheat areas and a high degree of surface fragmentation, which makes the extraction effect in central Anhui inferior to that in northern Anhui. Further, an optimum subset data was obtained from VIs, water index, building index and spectral bands using an RF algorithm. The result of using the optimum subset data showed a high accuracy of classification with a great advantage in data volume and processing time. This study provides a perspective for winter wheat mapping under various climatic and complicated land surface conditions and is of great significance for crop monitoring and agricultural decision-making.

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

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

Daily metre-scale mapping of water turbidity using CubeSat imagery

The potential for mapping of turbidity in inland and coastal waters using imagery from the PlanetScope (PS) and RapidEye (RE) constellations is evaluated. With >120 PS and 5 RE satellites in orbit these constellations are able to provide metre scale imagery on a daily basis and could significantly enhance high spatial resolution monitoring of turbidity worldwide. The Dark Spectrum Fitting (DSF) atmospheric correction is adapted to the PS and RE imaging systems to retrieve surface reflectances. Due to the large amount of imagery and the limited band sets on these sensors, automated pixel classification is required. This is here performed using a neural network approach, which is able to classify water pixels for clear to moderately turbid waters. Due to the limited band set and sensor performance, some issues remain with classifying extremely turbid waters and cloud shadows based on a spectral approach. Surface reflectance data compares well with in situ measurements from the AERONET-OC network. Turbidity is estimated from the Red, RedEdge (RE only) and NIR bands and is compared with measurements from autonomous stations in the San Francisco Bay area and the coastal waters around the United Kingdom. Good performance is found for Red band derived turbidity from PS data, while the NIR band performance is mediocre, likely due to calibration issues. For RE, all three turbidity products give reasonable results. A high revisit density allows for the mapping of temporal variability in water turbidity using these satellite constellations. Thanks to the RedEdge band on RE, chlorophyll a absorption can be avoided, and perhaps even estimated.

Cubesats Allow High Spatiotemporal Estimates of Satellite-Derived Bathymetry

High spatial and temporal resolution satellite remote sensing estimates are the silver bullet for monitoring of coastal marine areas globally. From 2000, when the first commercial satellite platforms appeared, offering high spatial resolution data, the mapping of coastal habitats and the extraction of bathymetric information have been possible at local scales. Since then, several platforms have offered such data, although not at high temporal resolution, making the selection of suitable images challenging, especially in areas with high cloud coverage. PlanetScope CubeSats appear to cover this gap by providing their relevant imagery. The current study is the first that examines the suitability of them for the calculation of the Satellite-derived Bathymetry. The availability of daily data allows the selection of the most qualitatively suitable images within the desired timeframe. The application of an empirical method of spaceborne bathymetry estimation provides promising results, with depth errors that fit to the requirements of the International Hydrographic Organization at the Category Zone of Confidence for the inclusion of these data in navigation maps. While this is a pilot study in a small area, more studies in areas with diverse water types are required for solid conclusions on the requirements and limitations of such approaches in coastal bathymetry estimations.

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.

Long-Term Annual Surface Water Change in the Brazilian Amazon Biome: Potential Links with Deforestation, Infrastructure Development and Climate Change

The Brazilian Amazon is one of the areas on the planet with the fastest changes in forest cover due to deforestation associated with agricultural expansion and infrastructure development. These drivers of change, directly and indirectly, affect the water ecosystem. In this study, we present a long-term spatiotemporal analysis of surface water annual change and address potential connections with deforestation, infrastructure expansion and climate change in this region. To do that, we used the Landsat Data Archive (LDA), and Earth Engine cloud computing platform, to map and analyze annual water changes between 1985 and 2017. We detected and estimated the extent of surface water using a novel sub-pixel classifier based on spectral mixture analysis, followed by a post-classification segmentation approach to isolate and classify surface water in natural and anthropic water bodies. Furthermore, we combined these results with deforestation and infrastructure development maps of roads, hydroelectric dams to quantify surface water changes linked with them. Our results showed that deforestation dramatically disrupts small streams, new hydroelectric dams inundated landmass after 2010 and that there is an overall trend of reducing surface water in the Amazon Biome and watershed scales, suggesting a potential connection to more recent extreme droughts in the 2010s.