USGS Science⌕ Search

SEARCH · USGS Science

Results for “Remote Sensing of Environment”

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

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

At least 235 records · Page 13Linked to original sources
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↗
USGS2023

Evaluation of Landsat image compositing algorithms

We proposed a new image compositing algorithm (MAX-RNB) based on the maximum ratio of Near Infrared (NIR) to Blue band (RNB), and evaluated it together with nine other compositing algorithms: MAX-NDVI (maximum Normalized Difference Vegetation Index), MED-NIR (median NIR band), WELD (conterminous United States Web-Enabled Landsat Data), BAP (Best Available Pixel), PAC (Phenology Adaptive Composite), WPS (Weighted Parametric Scoring), MEDOID (medoid measurement), COSSIM (cosine similarity), and NLCD (National Land Cover Database). Each algorithm was applied to time series of Landsat observations collected within two separate years at six locations around the world, to produce monthly (July 1 ± 15 days), seasonal (July 1 ± 45 days), and annual (July 1 ± 180 days) composite images free of cloud, cloud shadow, and snow/ice. By comparing the composite images to reference Landsat images acquired in the growing season (closest to July 1 within ±15 days) for each year, we evaluated the performance of the algorithms in preserving the spectral and spatial fidelity (hereafter referred to as spectral and spatial evaluation, respectively), as well as land cover classification and land change detection (hereafter referred to as application evaluation). The results demonstrated that no single algorithm outperformed all other algorithms in all the evaluations, but that performance depended on compositing intervals and cloud cover. For monthly composites, the MAX-RNB algorithm generally produced the best results in the spectral and application evaluations. For seasonal composites, the NLCD algorithm produced the best results in the spectral and application evaluations. For annual composites, the PAC algorithm produced the best results in the spectral evaluation and change detection, whereas BAP produced the best results in land cover classification. The BAP algorithm also produced the best results in the spatial evaluation for all the compositing periods. This study provides a comprehensive guidance for selecting the most appropriate image compositing algorithm for different Landsat-based applications.

Remote Sensing of Environment↗
USGS2023

Historical Structure from Motion (HSfM): Automated processing of historical aerial photographs for long-term topographic change analysis

Precisely measuring the Earth’s changing surface on decadal to centennial time scales is critical for many science and engineering applications, yet long-term records of quantitative landscape change are often temporally and geographically sparse. Archives of scanned historical aerial photographs provide an opportunity to augment these records with accurate elevation measurements that capture the historical state of the Earth surface. Structure from Motion (SfM) photogrammetry workflows produce high-quality digital elevation models (DEMs) and orthoimage mosaics from these historical images, but time-intensive tasks like manual image preprocessing (e.g., fiducial marker identification) and ground control point (GCP) selection impede processing at scale. We developed an automated method to process historical images and generate self-consistent time series of high-resolution (0.5–2 m) DEMs and orthomosaics, without manual GCP selection. The method relies on SfM to correct camera interior and exterior orientation and a robust multi-stage co-registration approach using modern reference terrain datasets for geolocation refinement. We demonstrate the method using scanned images from the North American Glacier Aerial Photography (NAGAP) archive collected between 1967 and 1997. We present results for two sites with variable photo acquisition geometry and overlap — Mount Baker and South Cascade Glacier in Washington State, USA. The automated method corrects initial camera position errors of several kilometers and produces accurately georeferenced, high-resolution DEMs and orthoimages, regardless of camera configuration, acquisition geometry, terrain characteristics, and reference DEM properties. The average RMS reprojection error following bundle adjustment optimization was 0.67 px (0.15 m) for the 261 images contributing to 10 final DEM mosaics between 1970 and 1992 at Mount Baker, and 0.65 px (0.13 m) for the 243 images contributing to 18 individual DEMs between 1967 and 1997 at South Cascade Glacier. The relative accuracy of elevation values in the historical time series stacks was 0.68 m at Mount Baker and 0.37 m at South Cascade Glacier. Our products have reduced systematic error and improved accuracy compared to DEM products generated using SfM with manual GCP selection. Final elevation change measurement precision was &#x223C; "> ∼ 0.7–1.0 m over a 30-year period, enabling the study of processes with rates as low as &#x223C; "> ∼ 1-3 cm/yr. Our results demonstrate the potential of this scalable method to rapidly process archives of historical imagery and deliver new quantitative insights on long-term geodetic change and Earth surface processes.

Washington↗
USGS2023

Elevation-based probabilistic mapping of irregularly flooded wetlands along the northern Gulf of Mexico coast

Irregularly flooded wetlands are found above the mean high water tidal datum and are exposed to tides and saltwater less frequently than daily. These wetlands provide important ecosystem services, such as providing habitat for fish and wildlife, enhancing water quality, ameliorating flooding impacts, supporting coastal food webs, and protecting upslope areas from erosion. Mapping irregularly flooded wetlands is challenging given their expansive coverage and dynamic nature. Furthermore, coastal wetlands are expected to change over the coming century due to sea-level rise and changes in the frequency and intensity of extreme storms. Consequently, coastal managers need baseline information on the spatial distribution of wetlands along with efficient and repeatable methods for observing changes. In this study, we used coastal wetlands from existing land use land cover data, best available lidar-derived digital elevation models, and Monte Carlo simulations to incorporate elevation uncertainty to create a probabilistic map of irregularly flooded wetlands along the northern Gulf of Mexico coast (USA). Our approach integrated findings from a review of coastal wetland elevation error in lidar datasets and an analysis of spatial autocorrelations of wetland elevation. We found a positive correlation ( r = 0.563, p < 0.0001) when comparing the probability estimated from a digital elevation model and in situ elevation observations. The differences in probability had a mean bias error of −0.04 (i.e., digital elevation model-based probability tends to be slightly lower), a mean absolute error of 0.20, and a root mean square error of 0.26. Beyond this overall validation, we explored error metrics for land cover classes and lidar collection details. To quantify areal coverage of the probabilistic output, we classified the probability values into equal bins using an interval of 0.33. The areal coverage of the lowest probability bin (“unlikely”; probability ≤0.33) was separated into the upper and lower portions of the irregularly flooded wetland zone. Of the coastal wetlands along the northern Gulf of Mexico coast about 38% were classified as unlikely and low with the greatest coverage in south Louisiana and the Everglades and around 33% were classified as unlikely and high with the greatest coverage in the Everglades and Texas. The relative coverage within the highest probability bin (“likely”; probability >0.66) covered around 13%, with the greatest coverage in south Florida, south Louisiana, and Texas. The framework developed in this study can be transferred to other coastal wetland areas and updated to observe changes with sea-level rise.

Gulf of Mexico coastline↗
USGS2023

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

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

continental United States↗
USGS2023

Evaluating and mitigating the impact of systematic geolocation error on canopy height measurement performance of GEDI

NASA's Global Ecosystem Dynamics Investigation (GEDI) is designed to provide high-resolution measurements of forest structure and topography between 52° N and S. However, current geolocation accuracy may limit further science applications of footprint-level products as early adopters have found it difficult to align with in-situ forestry inventory data and high-resolution imagery for calibration and validation purpose. Here we developed a new means to rapidly evaluate and mitigate the impact of systematic geolocation error on the performance of GEDI's forest height estimates in the US. By integrating nationwide high-resolution airborne lidar data collected through the 3D Elevation Program of the USGS, we provided optimal geolocation adjustments of GEDI at per beam level and tracked their performances over the first 18-mo. Our results suggest that the first release of GEDI product (R01) can have large systematic geolocation errors at beam level (i.e., 50.5% of beams with an error > 20 m). Its impact on canopy height measurement could drastically vary in space and time, which in turn also offers a separate indirect method to evaluate and track geolocation performance. The second release of GEDI data (R02) has achieved a much-improved systematic geolocation accuracy which is shown to meet the mission requirement (0.2% beams >20 m and 80.8% beams <10 m) and should be able to meet requirements from many practical science applications tolerant to moderate geolocation errors. In sum, our approach has provided a short-term solution for an enhanced Cal/Val strategy for GEDI. While further improvements will certainly be made in future releases, it can potentially create an alternative pathway to generate and validate biomass products by linking GEDI footprint samples directly with in-situ data collections.

Remote Sensing of Environment↗
USGS2023

The EnMAP imaging spectroscopy mission towards operations

EnMAP (Environmental Mapping and Analysis Program) is a high-resolution imaging spectroscopy remote sensing mission that was successfully launched on April 1st, 2022. Equipped with a prism-based dual-spectrometer, EnMAP performs observations in the spectral range between 418.2 nm and 2445.5 nm with 224 bands and a high radiometric and spectral accuracy and stability. EnMAP products, with a ground instantaneous field-of-view of 30 m x 30 m at a swath width of 30 km, allow for the qualitative and quantitative analysis of surface variables from frequently and consistently acquired observations on a global scale. This article presents the EnMAP mission and details the activities and results of the Launch and Early Orbit and Commissioning Phases until November 1st, 2022. The mission capabilities and expected performances for the operational Routine Phase are provided for existing and future EnMAP users.

Remote Sensing of Environment↗
USGS2023

Quantification of wetland vegetation communities features with airborne AVIRIS-NG, UAVSAR, and UAV LiDAR data in Peace-Athabasca Delta

Arctic-boreal wetlands, important ecosystems for biodiversity and ecological services, are experiencing hydrological changes including permafrost thaw, earlier snowmelt, and increased wildfire susceptibility. These changes are affecting wetland productivity, species diversity, and biogeochemical cycles. However, given the diverse forms and structures of wetland vegetation communities, traditional wetland maps generated from lower spatial and spectral resolution satellite imagery lack community-level vegetation classification and miss spatially complex patterns. In this study, we built a cloud-based workflow to map wetland vegetation community of the Peace-Athabasca Delta (PAD), Canada, by leveraging high-resolution (5-m) airborne multi-sensor datasets, namely NASA's Airborne Visible/Infrared Imaging Spectrometer-Next Generation (AVIRIS-NG) and Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR), and a historical LiDAR archive. Validation of our classifications using ground references indicates that classifications derived from AVIRIS-NG have higher accuracies (≥87.9%) than either UAVSAR (65.6%) or LiDAR (75.9%) for mapping wetland vegetation communities. We also show improved classification accuracy when combining information from multiple sensors. In particular, incorporating AVIRIS-NG and UAVSAR datasets substantially reduced omission errors of wet graminoid and wet shrub classes from 29.6% to 20.5% and from 10.8% to 7.5%, respectively. Combining AVIRIS-NG and LiDAR datasets further improves overall accuracy (+2.2%) for most classifications, especially emergent vegetation, wet graminoid, and wet shrub. The best performing model, using features derived from all three sensors, achieved an overall accuracy of 93.5%. The framework established here can be used to leverage extensive airborne AVIRIS-NG and UAVSAR datasets collected across Alaska and northwest Canada to understand the spatial distribution of Arctic-Boreal wetland vegetation communities.

Remote Sensing of Environment↗