SEARCH · Geology Reports

Results for “Remote Sensing of the Environment”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11Linked to original sources

USGS2021

Investigation of land surface phenology detections in shrublands using multiple scale satellite data

Shrublands occupy about 13% of the global land surface, contain about one-third of the biodiversity, store about half of the global terrestrial carbon, and provide many ecosystem services to a large amount of world's human population and livestock. Because phenology is a sensitive indicator of the response of shrubland ecosystems to climate change, the alteration of ecosystems following species invasions, and the dynamics of shrubland ecosystem function, biodiversity, and carbon budgets, it is critical to monitor and assess phenological dynamics in shrubland ecosystems. However, most current land surface phenology (LSP) products derived from satellite data do not provide phenology detections in some semiarid shrublands where the amplitude of seasonal vegetation greenness is small. Therefore, we investigated the LSP detection using multiple spatial resolution satellite data and examined the impacts of spatial scales and shrubland ecosystem components (shrub and herb cover) on LSP detections over the western United States. Specifically, greenup onset date (GUD) in shrublands was detected from 30 m Harmonized Landsat and Sentinel-2 (HLS) data and 500 m Visible Infrared Imaging Radiometer Suite (VIIRS) data to quantify scale effects. The GUD spatial patterns were explored with 30 m pixel variations in shrubland ecosystem components. The results show that GUD values varied with percent vegetation cover and shifted to earlier dates with increasing vegetation cover, demonstrating that satellite observations were not able to capture greenup onset until a threshold of green vegetation cover is reached. GUD was mostly undetectable from both HLS and VIIRS pixels with vegetation cover less than 10% and became fully detectable with vegetation covers larger than 50%. Similarly, the differences of GUD between HLS and VIIRS detections also decreased with increased vegetation cover. As a result of high shrubland heterogeneity, GUD from 30 m HLS pixels could be partially detected within a 500 m pixel despite GUD being undetectable from VIIRS time series. Moreover, vegetation cover heterogeneity also made it difficult for GUD at 30 m to be aggregated to coarse scales (such as to 500 m VIIRS pixels). These findings have significant implications to the detection and characterization of shrubland LSP responses to environmental and climate changes.

Nevada, Utah↗

USGS2021

Quantifying plant-soil-nutrient dynamics in rangelands: Fusion of UAV hyperspectral-LiDAR, UAV multispectral-photogrammetry, and ground-based LiDAR-digital photography in a shrub-encroached desert grassland

Rangelands cover 70% of the world's land surface, and provide critical ecosystem services of primary production, soil carbon storage, and nutrient cycling. These ecosystem services are governed by very fine-scale spatial patterning of soil carbon, nutrients, and plant species at the centimeter-to-meter scales, a phenomenon known as “islands of fertility”. Such fine-scale dynamics are challenging to detect with most satellite and manned airborne platforms. Remote sensing from unmanned aerial vehicles (UAVs) provides an alternative option for detecting fine-scale soil nutrient and plant species changes in rangelands tn0020 smaller extents. We demonstrate that a model incorporating the fusion of UAV multispectral and structure-from-motion photogrammetry classifies plant functional types and bare soil cover with an overall accuracy of 95% in rangelands degraded by shrub encroachment and disturbed by fire. We further demonstrate that employing UAV hyperspectral and LiDAR fusion greatly improves upon these results by classifying 9 different plant species and soil fertility microsite types (SFMT) with an overall accuracy of 87%. Among them, creosote bush and black grama, the most important native species in the rangeland, have the highest producer's accuracies at 98% and 94%, respectively. The integration of UAV LiDAR-derived plant height differences was critical in these improvements. Finally, we use synthesis of the UAV datasets with ground-based LiDAR surveys and lab characterization of soils to estimate that the burned rangeland potentially lost 1474 kg/ha of C and 113 kg/ha of N owing to soil erosion processes during the first year after a prescribed fire. However, during the second-year post-fire, grass and plant-interspace SFMT functioned as net sinks for sediment and nutrients and gained approximately 175 kg/ha C and 14 kg/ha N, combined. These results provide important site-specific insight that is relevant to the 423 Mha of grasslands and shrublands that are burned globally each year. While fire, and specifically post-fire erosion, can degrade some rangelands, post-fire plant-soil-nutrient dynamics might provide a competitive advantage to grasses in rangelands degraded by shrub encroachment. These novel UAV and ground-based LiDAR remote sensing approaches thus provide important details towards more accurate accounting of the carbon and nutrients in the soil surface of rangelands.

New Mexico↗

USGS2021

Detecting subtle change from dense landsat time series: Case studies of mountain pine beetle and spruce beetle disturbance

In contrast to abrupt changes caused by land cover conversion, subtle changes driven by a shift in the condition, structure, or other biological attributes of land often lead to minimal and slower alterations of the terrestrial surface. Accurate mapping and monitoring of subtle change are crucial for an early warning of long-term gradual change that may eventually result in land cover conversion. Freely accessible moderate-resolution datasets such as the Landsat archive have great potential to characterize subtle change by capturing low-magnitude spectral changes in long-term observations. However, past studies have reported limited success in accurately extracting subtle changes from satellite-based time series analysis. In this study, we introduce a supervised framework named ‘PIDS’ to detect subtle forest disturbance from a comprehensive Landsat data archive by leveraging disturbance-based calibration sites. PIDS consists of four components: (1) Parameter optimization; (2) Index selection; (3) Dynamic stratified monitoring; and (4) Spatial consideration. PIDS was applied to map the early stage of bark beetle infestations (i.e., a lower per-pixel fraction of trees cover that show visual signs of infestation), which are a typical example of subtle change in conifer forests. Landsat Analysis Ready Data were used as the time series inputs for mapping mountain pine beetle and spruce beetle disturbance between 2001 and 2019 in Colorado, USA. PIDS-detection map assessment showed that the overall performance of PIDS (namely ‘F1 score’) was 0.86 for mountain pine beetle and 0.73 for spruce beetle, making a substantial improvement (> 0.3) compared to other approaches/products including COntinuous monitoring of Land Disturbance, LandTrendr, and the National Land Cover Database forest disturbance product. A sub-pixel analysis of tree canopy mortality percentage was performed by linking classified high-resolution (0.3- and 1-m) aerial imagery and 30-m PIDS-detection maps. Results show that PIDS typically detects mountain pine beetle infestation when ≥56% of a Landsat pixel is occupied by red-stage canopy mortality (one year after initial infestation), and spruce beetle infestation when ≥55% is occupied by gray-stage mortality (two years after initial infestation). This study addresses an important methodological goal pertinent to the utility of event-based reference samples for detecting subtle forest change, which could be potentially applied to other types of subtle land change.

Colorado↗

USGS2021

Determination of burn severity models ranging from regional to continental scales for the conterminous United States

Identifying meaningful measures of ecological change over large areas is dependent on the quantification of robust relationships between ecological metrics and remote sensing products. Over the past several decades, ground observations of wildfire and prescribed fire severity have been acquired across hundreds of wildland fires in the United States, primarily utilizing the Composite Burn Index (CBI) plot protocol. These observations have been coupled to spaceborne passive spectral reflectance indices (e.g. Landsat-derived variations of the Normalized Burn Ratio [NBR]) to produce regression models describing their relationship. Here we develop regression models by vegetation type for multiple vegetation classification systems representing a range of spatial scales, and a decision tree framework for evaluating these regression models. Our overall goals were to determine which scale of ecological classifications provided the best estimate of burn severity from Landsat data and how to choose the best regression model. We aggregated a total of 6280 CBI plots for 234 wildland fires that burned between 1994 and 2017 and produced Landsat-derived NBR and differenced NBR (dNBR) values for each plot. We then calculated best fit linear or higher order regression equations between CBI and NBR/dNBR for each landcover classification system from smallest to largest scale: LANDFIRE Biophysical Settings (BPS), National Vegetation Classification macrogroup (NVC) landcover classifications, Omernick III, II, and I ecoregions, LANDFIRE Fire Regime Groups (FRG), and the entire conterminous United States (CONUS) dataset. The CONUS regression model goodness of fit was moderate (R 2 = 0.55, P < 0.001) for dNBR and poor (R 2 = 0.30, P < 0.001) for NBR. Within landcover classifications, CBI was better fit by dNBR than NBR. Finer scale regional regression models including BPS (dNBR R 2 ¯ "> R2¯ = 0.56 and 0.00–0.83 R 2 range; NBR R 2 ¯ "> R(bar) 2 = 0.43 and 0.00–0.82 R 2 range) and NVC (dNBR R 2 ¯ "> R(bar) 2 = 0.55 and 0.15–0.78 R 2 range; NBR R 2 ¯ "> R(bar) 2 = 0.41 and 0.00–0.79 R 2 range) were on average the same or better than the CONUS models for dNBR and NBR, with the strongest fit models exhibiting R 2 ≥ 0.70, whereas larger scale regional models R 2 ¯ "> R(bar) 2 ranged from 0.28 to 0.5. However, variation in accuracy among landcover types indicate that dNBR and NBR regression models could be used to effectively estimate CBI for future fires in certain regions, while for other regions models may require additional field observations or alternative spectral transformations. Our decision tree schema can be used to help users determine which scale is likely to produce the most accurate results using our models. The CBI regression models developed here, paired with the decision tree, provide users with a simple method to estimate burn severity in units of CBI for any fire within CONUS with moderate to high levels of confidence and provide a template for further development of models with new data going forward.

Conterminous United States↗

USGS2024

Mapping and classification of volcanic deposits using multi-sensor unoccupied aerial systems

The deposits from volcanic eruptions represent the record of activity at a volcano. Identification, classification, and interpretation of these deposits are crucial to the understanding of volcanic processes and assessing hazards. However, deposits often cover large areas and can be difficult or dangerous to access, making field mapping dangerous and time-consuming. Remote sensing techniques are often used to map and identify the deposits of volcanic eruptions, though these techniques present their own trade-offs in terms of image resolution, wavelength, and observation frequency. Here, we present a new approach for mapping and classifying volcanic deposits using a multi-sensor unoccupied aerial system (UAS), and demonstrate its application on lava and tephra deposits associated with the 2018 eruption of Sierra Negra volcano (Galápagos Archipelago, Ecuador). We surveyed the study area and collected visible and thermal infrared (TIR) images. We used structure-from-motion photogrammetry to create a digital elevation model (DEM) from the visual images and calculated the solar heating rate of the surface from temperature maps based on the TIR images. We find that the solar heating rate is highest for tephra deposits and lowest for ʻaʻā lava, with pāhoehoe lava having intermediate values. This is consistent with the solar heating rate correlating to the density and particle size of the surface. The solar heating rate for the lava flow also decreases with increasing distance from the vent, consistent with an increase in density as the lava degasses. We combined the surface roughness (calculated from the DEM) and the solar heating rate of the surface to remotely classify tephra deposits and different lava morphologies. We applied both supervised and unsupervised machine learning algorithms and demonstrate that supervised methods can replicate the manual classification while the unsupervised method can identify major surface units with no ground truth information. These methods allow for remote mapping and classification at high spatial resolution (< 1 meter) of a variety of volcanic deposits, with potential for application to deposits from other processes (e.g., fluvial, glacial) and deposits on other planetary bodies.

Galápagos Archipelago, Sierra Negra volcano↗

USGS2021

Validation of the U.S. Geological Survey’s Land Change Monitoring, Assessment and Projection (LCMAP) collection 1.0 annual land cover products 1985–2017

The U.S. Geological Survey Land Change Monitoring, Assessment and Projection (USGS LCMAP) has released a suite of annual land cover and land cover change products for the conterminous United States (CONUS). The accuracy of these products was assessed using an independently collected land cover reference sample dataset produced by analysts interpreting Landsat data, high-resolution aerial photographs, and other ancillary data. The reference sample of nearly 25,000 pixels and the accompanying 33-year time series of annual land cover reference labels allowed for a comprehensive assessment of accuracy of the LCMAP land cover and land cover change products. Overall accuracy (± standard error) for the per-pixel assessment across all years for the eight land cover classes was 82.5% (±0.2%). Overall accuracy was consistent year-to-year within a range of 1.5% but varied regionally with lower accuracy in the eastern United States. User’s accuracy (UA) and producer’s accuracy (PA) for CONUS ranged from the higher accuracies of Water (UA=96%, PA=93%) and Tree Cover (UA=90%, PA=83%) to the lower accuracies of Wetland (UA=69%, PA=74%) and Barren (UA=43%, PA=57%). For a binary change / no change classification, UA of change was 13% (±0.5%) and PA was 16% (±0.6%) for CONUS when agreement was defined as a match by the exact year of change. UA and PA improved to 28% and 34% when agreement was defined as the change being detected by the map and reference data within a ±2-year window. Change accuracy was higher in the eastern United States compared to the western US. UA was 49% (±0.3) and PA was 54% (±0.3) for the footprint of change (defined as the area experiencing at least one land cover change from 1985–2017). For class-specific loss and gain when agreement was defined as an exact year match, UA and PA were generally below 30%, with Tree Cover loss being the most accurately mapped change (UA=25%, PA=31%). These accuracy results provide users with information to assess the suitability of LCMAP data and information to guide future research for improving LCMAP products, particularly focusing on the challenges of accurately mapping annual land cover change.

Remote Sensing of Environment↗

USGS2021

Hydrological control shift from river level to rainfall in the reactivated Guobu slope besides the Laxiwa hydropower station in China

Landslides are common geohazards associated with natural drivers such as precipitation, land degradation, toe erosion by rivers and wave attack, and ground shaking. On the other hand, human alterations such as inundation by water impoundment or rapid drawdown may also destabilize the surrounding slopes. The Guobu slope is an ancient rockslide on the banks of the Laxiwa hydropower station reservoir (China), which reactivated during the reservoir impoundment in 2009. We extracted three-dimensional surface displacements with azimuth and range radar interferometry using European Space Agency's Copernicus Sentinel-1 and German Aerospace Center's TerraSAR-X data during 20152019. The upper part of the Guobu rockslide is characterized by toppling and is mostly subsiding with maximum rates over 0.4 m/yr and 0.7 m/yr in the vertical and horizontal directions, respectively. During filling of the reservoir prior to 2014, there was a long-wavelength in-phase response between rising reservoir level and GPS-observed increased slope movements. After the reservoir water level stabilized from 2015 to 2019, the slide movement became seasonal and we see a correlation between rainfall and landslide movement. These observations suggest that the slide motion is now primarily controlled by rainfall. The spatiotemporal landslide displacements allow us to estimate the hydraulic diffusivity of the rock mass, to be on the order (~1.05 × 10 ‐7 m 2 /s) and the thickness of the moving rock mass (~200 m). Our results demonstrate that InSAR is a useful tool for monitoring the rockslide movement as a function of seasonal precipitation.

Remote Sensing of Environment↗

USGS2021

A novel automatic phenology learning (APL) method of training sample selection using multiple datasets for time-series land cover mapping

The long record of Landsat imagery, which is the cornerstone of Earth observation, provides an opportunity to monitor land use and land cover (LULC) change and understand the interactions between the climate and earth system through time. A few change detection algorithms such as Continuous Change Detection and Classification (CCDC) have been developed to utilize all available Landsat images for change detection and characterization at local or global scales. However, the reliable, rapid, and reproducible collection of training samples have become a challenge for time series land cover classification at a large scale. To meet the challenge, we proposed an automatic phenology learning (APL) method with the assumption that the temporal profiles of samples within the same land cover type are the same or similar at a local scale to generate evenly distributed training samples automatically. We designed the method to build land cover patterns for each category based on consensus samples derived from multiple existing scientific datasets including LANDFIRE's (LF) Existing Vegetation Type (EVT), USGS National Land Cover Database (NLCD), National Agricultural Statistics Service (NASS) Cropland Data Layer (CDL), and National Wetlands Inventory (NWI). Then we calculated the Time-Weighted Dynamic Time Warping (twDTW) distance between any undefined samples and land cover patterns in the same geographical region as prior knowledge. Finally, we selected the optimal land cover category for each undefined sample from the land cover products based on the designed criteria iteratively using the twDTW distance as an indicator. The method was applied in the footprint of 10 selected Landsat Analysis Ready Data (ARD) tiles in the eastern and western conterminous United States (CONUS) to produce annual land cover maps from 1985 to 2017. The accuracy assessment and visual comparison revealed that the APL method can generate reliable training samples without any manual interpretation, producing better land cover results especially for the grass/shrub and wetland land cover classes. Applying the APL method, the overall accuracy of the annual land cover maps was improved by 2% over the accuracy of Land Change Monitoring, Assessment, and Projection (LCMAP) Collection 1.0 Science Products in the research regions. Our results also indicate that the APL method provides an approach for best use of different land cover products and meets the requirement of intensive sampling for training data collection.

Remote Sensing of Environment↗

USGS2022

Monitoring and characterizing multi-decadal variations of urban thermal condition using time-series thermal remote sensing and dynamic land cover data

Urban development and associated land cover and land use change alter the thermal, hydrological, and physical properties of the land surface. Assessments of surface urban heat island (UHI) usually focused on using remote sensing and land cover data to quantify UHI intensity and spatial distribution within a certain period. However, the mechanisms and complex interactions in landscape dynamics and land surface thermal features are still being assessed. In this study, we developed and implemented a novel approach to characterize landscape thermal conditions by focusing on UHI intensity and its spatiotemporal variation using the recently available time series of Landsat land surface temperature and land cover change products. We analyzed land surface temperature changes in urban and surrounding non-urban lands to quantify the UHI intensity and landscape thermal conditions in the Atlanta and Minneapolis metropolitan areas of the United States. Our results revealed that UHI intensities had averages of 3.4 °C and 3.3 °C in the Atlanta and Minneapolis metropolitan areas, respectively. The dominant land cover type in rural areas and urban imperviousness cover determines the UHI intensity. Increasing trends of 0.04 °C/year and 0.01 °C/year in UHI intensity between 1985 and 2018 were found in Atlanta and Minneapolis, respectively. The UHI intensity variations in 1985 and 2018 suggest that the magnitudes and temporal variations of UHI intensity averaged from all urban land cover classes are close to the UHI intensity estimated from the low intensity urban area only while the UHI intensities are more than 2 °C larger in medium to high and high intensity urban areas. The UHI intensities estimated from the maximum temperature that have statistically significant increasing trends suggest that the maximum temperature is a good element for measuring UHI effect. Urban land cover dynamics play an important role in controlling temporal variation of UHI and the UHI hotspots. Our findings support the scientific value of implementing the prototype approach as an objective framework to quantify and monitor UHI intensity at a large geographic extent.

Remote Sensing of Environment↗

USGS2022

Incorporating interpreter variability into estimation of the total variance of land cover area estimates under simple random sampling

Area estimates of land cover and land cover change are often based on reference class labels determined by analysts interpreting satellite imagery and aerial photography. Different interpreters may assign different reference class labels to the same sample unit. This interpreter variability is typically not accounted for in variance estimators applied to area estimates of land cover. A simple measurement model provides the basis for an estimator of the total variance ( V Total ) that takes into account both sampling variance and interpreter variance. This method requires two or more reference class interpretations (i.e., repeated measurements) obtained by analysts, working independently of each other, for the full sample or a random subsample of the full sample. Estimators of the total variance ( V ̂ Total "> V̂Total ) and the variance component attributable to interpreters ( V ̂ 1 "> V̂1 ) were obtained for the case of two reference class interpretations per repeated sample unit. To evaluate the effect of interpreter variability on variance estimation, we used land cover reference data interpreted by seven analysts who each interpreted the same 300 sample pixels from a region of the Pacific Northwest of the United States. From these data, we estimated the contribution of interpreter variance to the total variance (i.e., V ̂ 1 / V ̂ Total "> V̂1/V̂Total ) and the relative bias of the standard simple random sampling variance estimator ( V ̂ stand "> V̂stand ) as an estimator of V Total , defined as 100%*( V ̂ stand − V ̂ Total "> V̂stand−V̂Total )/ V ̂ Total "> V̂Total . For each of five land cover classes, we computed V ̂ 1 "> V̂1 , V ̂ Total "> V̂Total , and V ̂ stand "> V̂stand using the sample data from each of the 21 possible pairwise combinations of the seven interpreters, and then calculated the mean of V ̂ 1 / V ̂ Total "> V̂1/V̂Total and the mean of the estimated relative bias of V ̂ stand "> V̂stand over these 21 pairs. Based on the mean of V ̂ 1 / V ̂ Total "> V̂1/V̂Total per class, interpreter variance contributed from 25% (cropland) to 76% (grass/shrub) of the total variance, indicating that interpreter variance was a non-negligible component of the total variance. Typically, the standard variance estimator, V ̂ stand "> V̂stand , underestimated the total variance with the mean estimated relative bias ranging from −3% (cropland) to −33% (grass/shrub). Classes with greater inconsistency between pairs of interpreters had larger contributions of interpreter variance to the total variance ( V ̂ 1 / V ̂ Total "> V̂1/V̂Total ) and larger negative estimated relative bias of V ̂ stand "> V̂stand . Given that interpreter variance can contribute substantially to the total variance, the repeated measurements approach offers a practical way to incorporate this variability into an estimator of the total variance.

Washington↗

USGS2022

Hyperspectral remote sensing of white mica: A review of imaging and point-based spectrometer studies for mineral resources, with spectrometer design considerations

Over the past ~30 years, hyperspectral remote sensing of chemical variations in white mica have proven to be useful for ore deposit studies in a range of deposit types. To better understand mineral deposits and to guide spectrometer design, this contribution reviews relevant papers from the fields of remote sensing, spectroscopy, and geology that have utilized spectral changes caused by chemical variation in white micas. This contribution reviews spectral studies conducted at the following types of mineral deposits: base metal sulfide, epithermal, porphyry, sedimentary rock hosted gold deposits, orogenic gold, iron oxide copper gold, and unconformity-related uranium. The structure, chemical composition, and spectral features of white micas, in this contribution defined as muscovite, paragonite, celadonite, phengite, illite, and sericite, are given. Reviewed laboratory spectral studies determined that shifts in the position of the white mica 2200 nm combination feature of 1 nm correspond to a change in Al oct content of approximately ±1.05%. Many of the reviewed spectral studies indicated that a shift in the position of the white mica 2200 nm combination feature of 1 nm was geologically significant. A sensitivity analysis of spectrometer characteristics; bandpass, sampling interval, and channel position, is conducted using spectra of 19 white micas with deep absorption features to determine minimum characteristics required to accurately measure a shift in the position of the white mica 2200 nm combination feature. It was determined that a sampling interval < 16.3 nm and bandpass <17.5 nm are needed to achieve a root mean square error (RMSE) of 2 nm, whereas a sampling interval < 8.8 nm and bandpass <9.8 nm are needed to achieve a RMSE of 1 nm. For comparison, commonly used imaging spectrometers HyMap, AVIRIS-Classic, SpecTIR®'s AisaFENIX 1K, and HySpex tm SWIR 384 have 2.1, 1.2, 0.96, and 0.95 nm RMSE in determining the position of the 2200 nm white mica combination feature, respectively. An additional sensitivity analysis is conducted to determine the effect of signal to noise ratio (SNR) on the RMSE of the position of the white mica 2200 nm combination feature, using spectra of 18 white micas with deep absorption features. For a spectrometer with sampling interval and bandpass of 1 nm, we estimate that RMSEs of 1 and 1.5 nm are achievable with spectra having a minimum SNR of approximately 246 and 64, respectively. For a spectrometer with sampling interval and bandpass of 5 nm, we estimate that RMSEs of 1 and 1.5 nm are attainable with spectra having a minimum SNR of approximately 431 and 84, respectively. When using a spectrometer with a sampling interval 8.8 nm and a bandpass of 9.8 nm, a RMSE of 1 is only achievable with convolved, noiseless reference spectra. For the 8.8_9.8 nm spectrometer, spectra with SNR of 250 and 100 result in RMSE of 1.1 and 1.3, respectively. Therefore, fine spectral resolution characteristics achieve RMSEs better than 1 nm for high SNR spectra while spectrometers with coarse spectral resolution have larger RMSE, perform well with noisy data, and are useful for white mica studies if RMSE of 1.1 to 1.5 nm is acceptable.

Remote Sensing of Environment↗

USGS2022

Mapping actual evapotranspiration using Landsat for the conterminous United States: Google Earth Engine implementation and assessment of the SSEBop model

The estimation and mapping of actual evapotranspiration (ETa) is an active area of applied research in the fields of agriculture and water resources. Thermal remote sensing-based methods, using coarse resolution satellites, have been successful at estimating ETa over the conterminous United States (CONUS) and other regions of the world. In this study, we present CONUS-wide ETa from Landsat thermal imagery-using the Operational Simplified Surface Energy Balance (SSEBop) model in the Google Earth Engine (GEE) cloud computing platform. Over 150,000 Landsat satellite images were used to produce 10 years of annual ETa (2010–2019) at unprecedented scale. The accuracy assessment of the SSEBop results included point-based evaluation using monthly Eddy Covariance (EC) data from 25 AmeriFlux stations as well as basin-scale comparison with annual Water Balance ETa (WBET) for more than 1000 sub-basins. Evaluations using EC data showed generally mixed performance with weaker (R 2 < 0.6) correlation on sparsely vegetated surfaces such as grasslands or woody savanna and stronger correlation (R 2 > 0.7) over well-vegetated surfaces such as croplands and forests, but location-specific conditions rather than cover type were attributed to the variability in accuracy. Croplands performed best with R 2 of 0.82, root mean square error of 29 mm/month, and average bias of 12%. The WBET evaluation indicated that the SSEBop model is strong in explaining the spatial variability (up to R 2 > 0.90) of ETa across large basins, but it also identified broad hydro-climatic regions where the SSEBop ETa showed directional biases, requiring region-specific model parameter improvement and/or bias correction with an overall 7% bias nationwide. Annual ETa anomalies over the 10-year period captured widely reported drought-affected regions, for the most part, in different parts of the CONUS, indicating their potential applications for mapping regional- and field-scale drought and fire effects. Due to the coverage of the Landsat Path/Row system, the availability of cloud-free image pixels ranged from less than 12 (mountainous cloud-prone regions and U.S. Northeast) to more than 60 (U.S. Southwest) per year. However, this study reinforces a promising application of Landsat satellite data with cloud-computing for quick and efficient mapping of ETa for agricultural and water resources assessments at the field scale.

Remote Sensing of Environment↗

USGS2022

Detection and characterization of coastal tidal wetland change in the northeastern US using Landsat time series

Coastal tidal wetlands are highly altered ecosystems exposed to substantial risk due to widespread and frequent land-use change coupled with sea-level rise, leading to disrupted hydrologic and ecologic functions and ultimately, significant reduction in climate resiliency. Knowing where and when the changes have occurred, and the nature of those changes, is important for coastal communities and natural resource management . Large-scale mapping of coastal tidal wetland changes is extremely difficult due to their inherent dynamic nature. To bridge this gap, we developed an automated algorithm for DEtection and Characterization of cOastal tiDal wEtlands change (DECODE) using dense Landsat time series. DECODE consists of three elements, including spectral break detection, land cover classification and change characterization. DECODE assembles all available Landsat observations and introduces a water level regressor for each pixel to flag the spectral breaks and estimate harmonic time-series models for the divided temporal segments. Each temporal segment is classified (e.g., vegetated wetlands, open water, and others – including unvegetated areas and uplands) based on the phenological characteristics and the synthetic surface reflectance values calculated from the harmonic model coefficients, as well as a generic rule-based classification system. This harmonic model-based approach has the advantage of not needing the acquisition of satellite images at optimal conditions (i.e., low tide status) to avoid underestimating coastal vegetation caused by the tidal fluctuation. At the same time, DECODE can also characterize different kinds of changes including land cover change and condition change (i.e., land cover modification without conversion). We used DECODE to track status of coastal tidal wetlands in the northeastern United States from 1986 to 2020. The overall accuracy of land cover classification and change detection is approximately 95.8% and 99.8%, respectively. The vegetated wetlands and open water were mapped with user's accuracy of 94.6% and 99.0%, and producer's accuracy of 98.1% and 93.5%, respectively. The cover change and condition change were mapped with user's accuracy of 68.0% and 80.0%, and producer's accuracy of 80.5% and 97.1%, respectively. Approximately 3283 km 2 of the coastal landscape within our study area in the northeastern United States changed at least once (12% of the study area), and condition changes were the dominant change type (84.3%). Vegetated coastal tidal wetland decreased consistently (~2.6 km 2 per year) in the past 35 years, largely due to conversion to open water in the context of sea-level rise.

Remote Sensing of Environment↗

USGS2022

Spectral mixture analysis for surveillance of harmful algal blooms (SMASH): A field-, laboratory-, and satellite-based approach to identifying cyanobacteria genera from remotely sensed data

Algal blooms around the world are increasing in frequency and severity, often with the possibility of adverse effects on human and ecosystem health. The health and economic impacts associated with harmful algal blooms, or HABs, provide compelling rationale for developing new methods for monitoring these events via remote sensing. Although concentrations of chlorophyll- a and key pigments like phycocyanin are routinely estimated from satellite images and used to infer algal or cyanobacterial cell counts, current methods are unable to provide information on the taxonomic composition of a bloom. This study introduced a new approach capable of differentiating among genera based on their reflectance characteristics: Spectral Mixture Analysis for Surveillance of HABs, or SMASH. The foundation of SMASH is a multiple endmember spectral mixture analysis (MESMA) algorithm that takes a library of cyanobacteria endmembers and a hyperspectral image as input and estimates the fractional abundance of each genus, plus water, on a per-pixel basis. Importantly, we assume that the water column consists of only pure water and cyanobacteria, implying that our linear spectral unmixing models do not account for other optically active constituents such as suspended sediment and colored dissolved organic matter (CDOM). We used reflectance spectra for 12 genera measured under a microscope to populate an algal spectral library and applied the SMASH workflow to satellite images from four waterbodies across the United States. Normalized spectral separability scores indicated that the 12 genera were distinct from one another and the MESMA algorithm reproduced known input fractions for simulated mixtures that included all pairwise combinations of genera and water. We used Upper Klamath Lake as an example to illustrate data products generated via SMASH: maps of the normalized difference chlorophyll index and cyanobacterial index, a MESMA-based classification of algal genera, fraction images for each endmember, and a root mean square error (RMSE) image that summarizes uncertainty. For Upper Klamath Lake, these outputs highlighted a complex algal bloom featuring several genera, primarily Aphanizomenon , and intricate spatial patterns associated with gyres. The maximum RMSE constraint imposed on the MESMA algorithm provided a means of avoiding false positive detection of genera not present in a waterbody but must not be set so low as to leave much of an image unclassified in cases where genera included in the library are present. Comparison of endmember fractions with relative biovolumes calculated from field samples indicated that taxonomic information from SMASH was consistent with field observations. For example, the algorithm successfully identified Microcystis in Owasco Lake but avoided misclassifying Asterionella , a genus not yet included in our library, in Detroit Lake. This proof-of-concept investigation demonstrates the potential of SMASH to enhance our understanding of algal blooms, particularly with respect to their spatial and temporal dynamics.

New York, Oregon, Texas↗

USGS2022

Fifty years of Landsat science and impacts

Since 1972, the Landsat program has been continually monitoring the Earth, to now provide 50 years of digital, multispectral, medium spatial resolution observations. Over this time, Landsat data were crucial for many scientific and technical advances. Prior to the Landsat program, detailed, synoptic depictions of the Earth's surface were rare, and the ability to acquire and work with large datasets was limited. The early years of the Landsat program delivered a series of technological breakthroughs, pioneering new methods, and demonstrating the ability and capacity of digital satellite imagery, creating a template for other global Earth observation missions and programs. Innovations driven by the Landsat program have paved the way for subsequent science, application, and policy support activities. The economic and scientific value of the knowledge gained through the Landsat program has been long recognized, and despite periods of funding uncertainty, has resulted in the program's 50 years of continuity, as well as substantive and ongoing improvements to payload and mission performance. Free and open access to Landsat data, enacted in 2008, was unprecedented for medium spatial resolution Earth observation data and substantially increased usage and led to a proliferation of science and application opportunities. Here, we highlight key developments over the past 50 years of the Landsat program that have influenced and changed our scientific understanding of the Earth system. Major scientific and programmatic impacts have been realized in the areas of agricultural crop mapping and water use, climate change drivers and impacts, ecosystems and land cover monitoring, and mapping the changing human footprint. The introduction of Landsat collection processing, coupled with the free and open data policy, facilitated a transition in Landsat data usage away from single images and towards time series analyses over large areas and has fostered the widespread use of science-grade data. The launch of Landsat-9 on September 27, 2021, and the advanced planning of its successor mission, Landsat-Next, underscore the sustained institutional support for the program. Such support and commitment to continuity is recognition of both the historic impact the program, and the future potential to build upon Landsat's remarkable 50-year legacy.

Remote Sensing of Environment↗

USGS2022

Spaceborne InSAR mapping of landslides and subsidence in rapidly deglaciating terrain, Glacier Bay National Park and Preserve and vicinity, Alaska and British Columbia

The Glacier Bay area in southeastern Alaska and British Columbia, encompassing Glacier Bay National Park and Preserve, has experienced rapid glacier retreat since the end of the Little Ice Age in the mid-1800s. The impact that rapid deglaciation has had on the slope stability of valley walls and on the sedimentation of fans and deltas adjacent to fjords and inlets is an ongoing research topic. Using 3-year (2018–2020) Sentinel-1 datasets, and an automated time-series persistent scatterer interferometric synthetic aperture radar (PSInSAR) processing method, we detected landslides or delta subsidence at 27 sites within a vast 180 × 180 km remote region encompassing Glacier Bay proper. Most of the sites that we identified had not been previously identified. We categorized the hazard source areas that we identified into three general types:1) slow-moving landslides on steep rocky slopes not near (> 2 km away from) present-day glacier termini (e.g., Tidal Inlet landslide), 2) slow-moving landslides directly adjacent to (< 2 km away from), and associated with glacier thinning and retreat, and 3) subsidence of glacial outwash fan deltas. In categories 1 and 2, we observed 22 landslides moving at velocities ranging from 0.5 to 4 cm/yr. In category 3, we detected five fan deltas subsiding at velocities ranging from 0.5 to 6 cm/yr. Within our measurement error, these velocities were consistent during the monitoring period. Because acceleration was not observed, the issuance of warnings of imminent rapid failure is not warranted, however, continued remote monitoring is warranted. Our interferometric synthetic aperture radar (InSAR) results could be combined with other data sets including field observations, subaerial and submarine landslide inventories, bedrock fabric mapping from newly available light detection and ranging (lidar) data, and geologic maps to produce an inherent susceptibility map for landslides in bedrock and fan deltas. This map could be used to forecast susceptibility for both earthquake and climatically induced landslides.

Alaska, British Columbia↗

USGS2022

Uncertainty of ICESat-2 ATL06- and ATL08-derived snow depths for glacierized and vegetated mountain regions

Seasonal snow melt dominates the hydrologic budget across a large portion of the globe. Snow accumulation and melt vary over a broad range of spatial scales, preventing accurate extrapolation of sparse in situ observations to watershed scales. The lidar onboard the Ice, Cloud, and land Elevation, Satellite (ICESat-2) was designed for precise mapping of ice sheets and sea ice, and here we assess the feasibility of snow depth-mapping using ICESat-2 data in more complex and rugged mountain landscapes. We explore the utility of ATL08 Land and Vegetation Height and ATL06 Land Ice Height differencing from reference elevation datasets in two end member study sites. We analyze ∼3 years of data for Reynolds Creek Experimental Watershed in Idaho's Owyhee Mountains and Wolverine Glacier in southcentral Alaska's Kenai Mountains. Our analysis reveals decimeter-scale uncertainties in derived snow depth and glacier mass balance at the watershed scale. Both accuracy and precision decrease as slope increases: the magnitudes of the median and median of the absolute deviation of elevation errors (MAD) vary from ∼0.2 m for slopes <5° to >1 m for slopes >20°. For glacierized regions, failure to account for intra- and inter-annual evolution of glacier surface elevations can strongly bias ATL06 elevations, resulting in under-estimation of the mass balance gradient with elevation. Based on these results, we conclude that ATL08 and ATL06 observations are best suited for characterization of watershed-scale snow depth and mass balance gradients over relatively shallow slopes with thick snowpacks . In these regions, ICESat-2 elevation residual-derived snow depth and mass balance transects can provide valuable watershed scale constraints on terrain parameter- and model-derived estimates of snow accumulation and melt.

Alaska, Idaho↗

USGS2023

Toward consistent change detection across irregular remote sensing time series observations

The use of remote sensing in time series analysis enables wall-to-wall monitoring of the land surface and is critical for assessing and understanding land cover and land use change and for understanding the Earth system as a whole. However, variability in remote sensing observation frequency through time and across space presents challenges for producing consistent change detection results throughout the available satellite record using approaches such as the Continuous Change Detection and Classification (CCDC) change detection methodology. Here we investigate new modifications to this methodology with the goal of improving accuracy and consistency in results and increasing flexibility for operational usage and future development. The modified method (Band-First Probability, or CCD-BFP) change detection procedure works by calculating a test for each band through time before summarizing between bands. We evaluate the CCD-BFP method compared to an existing implementation of CCDC using a variety of approaches, including a validation dataset of human-interpreted locations, comparison with data from fire events, use of simulated remote sensing data, and qualitative inspection of areas of interest. We find CCD-BFP improves consistency across time and space compared to the existing implementation of CCDC, with more similarity in rates of change across Landsat swath boundaries and before and after the launch of Landsat 7. Also, we find that CCD-BFP detects more of the change events in the validation dataset while reducing the overall number of change detections, indicating that it is able to more accurately capture the most notable land surface change events.

Remote Sensing of Environment↗