USGS Science⌕ Search

SEARCH · USGS Science

Results for “Algorithms”

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 541 records · Page 30Linked to original sources
USGS2021

Surface flow velocities from space: Particle image velocimetry of satellite video of a large, sediment-laden river

Conventional, field-based streamflow monitoring in remote, inaccessible locations such as Alaska poses logistical challenges. Safety concerns, financial considerations, and a desire to expand water-observing networks make remote sensing an appealing alternative means of collecting hydrologic data. In an ongoing effort to develop non-contact methods for measuring river discharge, we evaluated the potential to estimate surface flow velocities from satellite video of a large, sediment-laden river in Alaska via particle image velocimetry (PIV). In this setting, naturally occurring sediment boil vortices produced distinct water surface features that could be tracked from frame to frame as they were advected by the flow, obviating the need to introduce artificial tracer particles. In this study, we refined an end-to-end workflow that involved stabilization and geo-referencing, image preprocessing, PIV analysis with an ensemble correlation algorithm, and post-processing of PIV output to filter outliers and scale and geo-reference velocity vectors. Applying these procedures to image sequences extracted from satellite video allowed us to produce high resolution surface velocity fields; field measurements of depth-averaged flow velocity were used to assess accuracy. Our results confirmed the importance of preprocessing images to enhance contrast and indicated that lower frame rates (e.g., 0.25 Hz) lead to more reliable velocity estimates because longer capture intervals allow more time for water surface features to translate several pixels between frames, given the relatively coarse spatial resolution of the satellite data. Although agreement between PIV-derived velocity estimates and field measurements was weak ( R 2 = 0.39) on a point-by-point basis, correspondence improved when the PIV output was aggregated to the cross-sectional scale. For example, the correspondence between cross-sectional maximum velocities inferred via remote sensing and measured in the field was much stronger ( R 2 = 0.76), suggesting that satellite video could play a role in measuring river discharge. Examining correlation matrices produced as an intermediate output of the PIV algorithm yielded insight on the interactions between image frame rate and sensor spatial resolution, which must be considered in tandem. Although further research and technological development are needed, measuring surface flow velocities from satellite video could become a viable tool for streamflow monitoring in certain fluvial environments.

Alaska↗
USGS2018

Probabilistic substrate classification with multispectral acoustic backscatter: A comparison of discriminative and generative models

We propose a probabilistic graphical model for discriminative substrate characterization, to support geological and biological habitat mapping in aquatic environments. The model, called a fully-connected conditional random field (CRF), is demonstrated using multispectral and monospectral acoustic backscatter from heterogeneous seafloors in Patricia Bay, British Columbia, and Bedford Basin, Nova Scotia. Unlike previously proposed discriminative algorithms, the CRF model considers both the relative backscatter magnitudes of different substrates and their relative proximities. The model therefore combines the statistical flexibility of a machine learning algorithm with an inherently spatial treatment of the substrate. The CRF model predicts substrates such that nearby locations with similar backscattering characteristics are likely to be in the same substrate class. The degree of allowable proximity and backscatter similarity are controlled by parameters that are learned from the data. CRF model results were evaluated against a popular generative model known as a Gaussian Mixture model (GMM) that doesn’t include spatial dependencies, only covariance between substrate backscattering response over different frequencies. Both models are used in conjunction with sparse bed observations/samples in a supervised classification. A detailed accuracy assessment, including a leave-one-out cross-validation analysis, was performed using both models. Using multispectral backscatter, the GMM model trained on 50% of the bed observations resulted in a 75% and 89% average accuracies in Patricia Bay and Bedford Basin, respectively. The same metrics for the CRF model were 78% and 95%. Further, the CRF model resulted in a 91% mean cross-validation accuracy across four substrate classes at Patricia Bay, and a 99.5% mean accuracy across three substrate classes at Bedford Basin, which suggest that the CRF model generalizes extremely well to new data. This analysis also showed that the CRF model was much less sensitive to the specific number and locations of bed observations than the generative model, owing to its ability to incorporate spatial autocorrelation in substrates. The CRF therefore may prove to be a powerful ‘spatially aware’ alternative to other discriminative classifiers.

Geosciences↗
USGS2011

Comparison of Bayesian clustering and edge detection methods for inferring boundaries in landscape genetics

Recently, techniques available for identifying clusters of individuals or boundaries between clusters using genetic data from natural populations have expanded rapidly. Consequently, there is a need to evaluate these different techniques. We used spatially-explicit simulation models to compare three spatial Bayesian clustering programs and two edge detection methods. Spatially-structured populations were simulated where a continuous population was subdivided by barriers. We evaluated the ability of each method to correctly identify boundary locations while varying: (i) time after divergence, (ii) strength of isolation by distance, (iii) level of genetic diversity, and (iv) amount of gene flow across barriers. To further evaluate the methods’ effectiveness to detect genetic clusters in natural populations, we used previously published data on North American pumas and a European shrub. Our results show that with simulated and empirical data, the Bayesian spatial clustering algorithms outperformed direct edge detection methods. All methods incorrectly detected boundaries in the presence of strong patterns of isolation by distance. Based on this finding, we support the application of Bayesian spatial clustering algorithms for boundary detection in empirical datasets, with necessary tests for the influence of isolation by distance

International Journal of Molecular Sciences↗
USGS2019

Monitoring landscape dynamics in central U.S. grasslands with harmonized Landsat-8 and Sentinel-2 time series data

Remotely monitoring changes in central U.S. grasslands is challenging because these landscapes tend to respond quickly to disturbances and changes in weather. Such dynamic responses influence nutrient cycling, greenhouse gas contributions, habitat availability for wildlife, and other ecosystem processes and services. Traditionally, coarse-resolution satellite data acquired at daily intervals have been used for monitoring. Recently, the harmonized Landsat-8 and Sentinel-2 (HLS) data increased the temporal frequency of the data. Here we investigated if the increased data frequency provided adequate observations to characterize highly dynamic grassland processes. We evaluated HLS data available for 2016 to (1) determine if data from Sentinel-2 contributed to an improvement in characterizing landscape processes over Landsat-8 data alone, and (2) quantify how observation frequency impacted results. Specifically, we investigated into estimating annual vegetation phenology, detecting burn scars from fire, and modeling within-season wetland hydroperiod and growth of aquatic vegetation. We observed increased sensitivity to the start of the growing season (SOST) with the HLS data. Our estimates of the grassland SOST compared well with ground estimates collected at a phenological camera site. We used the Continuous Change Detection and Classification (CCDC) algorithm to assess if the HLS data improved our detection of burn scars following grassland fires and found that detection was considerably influenced by the seasonal timing of the fires. The grassland burned in early spring recovered too quickly to be detected as change events by CCDC; instead, the spectral characteristics following these fires were incorporated as part of the ongoing time-series models. In contrast, the spectral effects from late-season fires were detected both by Landsat-8 data and HLS data. For wetland-rich areas, we used a modified version of the CCDC algorithm to track within-season dynamics of water and aquatic vegetation. The addition of Sentinel-2 data provided the potential to build full time series models to better distinguish different wetland types, suggesting that the temporal density of data was sufficient for within-season characterization of wetland dynamics. Although the different data frequency, in both the spatial and temporal dimensions, could cause inconsistent model estimation or sensitivity sometimes; overall, the temporal frequency of the HLS data improved our ability to track within-season grassland dynamics and improved results for areas prone to cloud contamination. The results suggest a greater frequency of observations, such as from harmonizing data across all comparable Landsat and Sentinel sensors, is still needed. For our study areas, at least a 3-day revisit interval during the early growing season (weeks 14–17) is required to provide a >50% probability of obtaining weekly clear observations.

Minnesota↗
USGS2019

Overall methodology design for the United States National Land Cover Database 2016 products

The National Land Cover Database (NLCD) 2016 provides a suite of data products, including land cover and land cover change of the conterminous United States from 2001 to 2016, at two- to three-year intervals. The development of this product is part of an effort to meet the growing demand for longer temporal duration and more frequent, accurate, and consistent land cover and change information. To accomplish this, we designed a new land cover strategy and developed comprehensive methods, models, and procedures for NLCD 2016 implementation. Major steps in the new procedures consist of data preparation, land cover change detection and classification, theme-based postprocessing, and final integration. Data preparation includes Landsat imagery selection, cloud detection, and cloud filling, as well as compilation and creation of more than 30 national-scale ancillary datasets. Land cover change detection includes single-date water and snow/ice detection algorithms and models, two-date multi-index integrated change detection models, and long-term multi-date change algorithms and models. The land cover classification includes seven-date training data creation and 14-run classifications. Pools of training data for change and no-change areas were created before classification based on integrated information from ancillary data, change-detection results, Landsat spectral and temporal information, and knowledge-based trajectory analysis. In postprocessing, comprehensive models for each land cover theme were developed in a hierarchical order to ensure the spatial and temporal coherence of land cover and land cover changes over 15 years. An initial accuracy assessment on four selected Landsat path/rows classified with this method indicates an overall accuracy of 82.0% at an Anderson Level II classification and 86.6% at the Anderson Level I classification after combining the primary and alternate reference labels. This methodology was used for the operational production of NLCD 2016 for the Conterminous United States, with final produced products available for free download.

Remote Sensing↗
USGS2020

Gap fill of Land surface temperature and reflectance products in Analysis Ready Data

The recently released Landsat Analysis Ready Data (ARD) over the United States provides the opportunity to investigate landscape dynamics using dense time series observations at 30-m resolution. However, the dataset often contains data gaps (or missing data) because of cloud contamination or data acquisition strategy. We present a new algorithm that focuses on data gap filling using clear observations from orbit overlap regions. Multiple linear regression models were established for each pixel time series to estimate stable predictions and uncertainties. The model's training data came from stratified random samples based on the time series similarity between the pixel and data from the overlap regions. The algorithm was evaluated using four tiles (5,000 × 5,000 30-m pixels for each tile) from 2018 land surface temperature data (LST) in Atlanta, Georgia. The accuracy was assessed using 1,000 randomly masked pixels and daily air temperature from eight ground stations. Both assessments showed the r2 value above 0.75, except two stations with mixed Landsat pixels. We also compared our results with the eMODIS LST product in terms of annual mean temperature. The two maps showed a similar spatial pattern at the region level, but our results showed more spatial detail in the urban area that matched the pattern of impervious surface. We also applied the method on ARD surface reflectance bands at Fairbanks, Alaska, to illustrate its improvements in surface reflectance products and in land change modeling. This approach can also be applied to other datasets, vegetation indexes, or spectral reflectance bands of other sensors.

Remote Sensing↗
USGS2020

Inferring surface flow velocities in sediment-laden Alaskan rivers from optical image sequences acquired from a helicopter

The remote, inaccessible location of many rivers in Alaska creates a compelling need for remote sensing approaches to streamflow monitoring. Motivated by this objective, we evaluated the potential to infer flow velocities from optical image sequences acquired from a helicopter deployed above two large, sediment-laden rivers. Rather than artificial seeding, we used an ensemble correlation particle image velocimetry (PIV) algorithm to track the movement of boil vortices that upwell suspended sediment and produce a visible contrast at the water surface. This study introduced a general, modular workflow for image preparation (stabilization and geo-referencing), preprocessing (filtering and contrast enhancement), analysis (PIV), and postprocessing (scaling PIV output and assessing accuracy via comparison to field measurements). Applying this method to images acquired with a digital mapping camera and an inexpensive video camera highlighted the importance of image enhancement and the need to resample the data to an appropriate, coarser pixel size and a lower frame rate. We also developed a Parameter Optimization for PIV (POP) framework to guide selection of the interrogation area (IA) and frame rate for a particular application. POP results indicated that the performance of the PIV algorithm was highly robust and that relatively large IAs (64–320 pixels) and modest frame rates (0.5–2 Hz) yielded strong agreement ( R 2 > 0 . 9 "> R 2 > 0 . 9 ) between remotely sensed velocities and field measurements. Similarly, analysis of the sensitivity of PIV accuracy to image sequence duration showed that dwell times as short as 16 s would be sufficient at a frame rate of 1 Hz and could be cut in half if the frame rate were doubled. The results of this investigation indicate that helicopter-based remote sensing of velocities in sediment-laden rivers could contribute to noncontact streamgaging programs and enable reach-scale mapping of flow fields.

Alaska↗
USGS2021

Exploring VIIRS continuity with MODIS in an expedited capability for monitoring drought-related vegetation conditions

Vegetation has been effectively monitored using remote sensing time-series vegetation index (VI) data for several decades. Drought monitoring has been a common application with algorithms tuned to capturing anomalous temporal and spatial vegetation patterns. Drought stress models, such as the Vegetation Drought Response Index (VegDRI), often use VIs like the Normalized Difference Vegetation Index (NDVI). The EROS expedited Moderate Resolution Imaging Spectroradiometer (eMODIS)-based, 7-day NDVI composites are integral to the VegDRI. As MODIS satellite platforms (Terra and Aqua) approach mission end, the Visible Infrared Imaging Radiometer Suite (VIIRS) presents an alternate NDVI source, with daily collection, similar band passes, and moderate spatial resolution. This study provides a statistical comparison between EROS expedited VIIRS (eVIIRS) 375-m and eMODIS 250-m and tests the suitability of replacing MODIS NDVI with VIIRS NDVI for drought monitoring and vegetation anomaly detection. For continuity with MODIS NDVI, we calculated a geometric mean regression adjustment algorithm using 375-m resolution for an eMODIS-like NDVI (eVIIRS’) eVIIRS’ = 0.9887 × eVIIRS − 0.0398. The resulting statistical comparisons (eVIIRS’ vs. eMODIS NDVI) showed correlations consistently greater than 0.84 throughout the three years studied. The eVIIRS’ VegDRI results characterized similar drought patterns and hotspots to the eMODIS-based VegDRI, with near zero bias.

Kansas↗
USGS2021

Estimating actual evapotranspiration over croplands using vegetation index methods and dynamic harvested area

Advances in estimating actual evapotranspiration (ETa) with remote sensing (RS) have contributed to improving hydrological, agricultural, and climatological studies. In this study, we evaluated the applicability of Vegetation-Index (VI) -based ETa (ET-VI) for mapping and monitoring drought in arid agricultural systems in a region where a lack of ground data hampers ETa work. To map ETa (2000–2019), ET-VIs were translated and localized using Landsat-derived 3- and 2-band Enhanced Vegetation Indices (EVI and EVI2) over croplands in the Zayandehrud River Basin (ZRB) in Iran. Since EVI and EVI2 were optimized for the MODerate Imaging Spectroradiometer (MODIS), using these VIs with Landsat sensors required a cross-sensor transformation to allow for their use in the ET-VI algorithm. The before- and after- impact of applying these empirical translation methods on the ETa estimations was examined. We also compared the effect of cropping patterns’ interannual change on the annual ETa rate using the maximum Normalized Difference Vegetation Index (NDVI) time series. The performance of the different ET-VIs products was then evaluated. Our results show that ETa estimates agreed well with each other and are all suitable to monitor ETa in the ZRB. Compared to ETc values, ETa estimations from MODIS-based continuity corrected Landsat-EVI (EVI2) (EVI MccL and EVI2 MccL ) performed slightly better across croplands than those of Landsat-EVI (EVI2) without transformation. The analysis of harvested areas and ET-VIs anomalies revealed a decline in the extent of cultivated areas and a loss of corresponding water resources downstream. The findings show the importance of continuity correction across sensors when using empirical algorithms designed and optimized for specific sensors. Our comprehensive ETa estimation of agricultural water use at 30 m spatial resolution provides an inexpensive monitoring tool for cropping areas and their water consumption.

Remote Sensing↗
USGS2022

Deep and machine learning image classification of coastal wetlands using unpiloted aircraft system multispectral images and lidar datasets

The recent developments of new deep learning architectures create opportunities to accurately classify high-resolution unoccupied aerial system (UAS) images of natural coastal systems and mandate continuous evaluation of algorithm performance. We evaluated the performance of the U-Net and DeepLabv3 deep convolutional network architectures and two traditional machine learning techniques (support vector machine (SVM) and random forest (RF)) applied to seventeen coastal land cover types in west Florida using UAS multispectral aerial imagery and canopy height models (CHM). Twelve combinations of spectral bands and CHMs were used. Our results using the spectral bands showed that the U-Net (83.80–85.27% overall accuracy) and the DeepLabV3 (75.20–83.50% overall accuracy) deep learning techniques outperformed the SVM (60.50–71.10% overall accuracy) and the RF (57.40–71.0%) machine learning algorithms. The addition of the CHM to the spectral bands slightly increased the overall accuracy as a whole in the deep learning models, while the addition of a CHM notably improved the SVM and RF results. Similarly, using bands outside the three spectral bands, namely, near-infrared and red edge, increased the performance of the machine learning classifiers but had minimal impact on the deep learning classification results. The difference in the overall accuracies produced by using UAS-based lidar and SfM point clouds, as supplementary geometrical information, in the classification process was minimal across all classification techniques. Our results highlight the advantage of using deep learning networks to classify high-resolution UAS images in highly diverse coastal landscapes. We also found that low-cost, three-visible-band imagery produces results comparable to multispectral imagery that do not risk a significant reduction in classification accuracy when adopting deep learning models.

Florida↗
USGS2024

Physics-based satellite-derived bathymetry (SDB) using Landsat OLI images

The estimation of depth in optically shallow waters using satellite imagery can be efficient and cost-effective. Active sensors measure the distance traveled by an emitted laser pulse propagating through the water with high precision and accuracy if the bottom peak intensity of the waveform is greater than the noise level. However, passive optical imaging of optically shallow water involves measuring the radiance after the sunlight undergoes downward attenuation on the way to the sea floor, and the reflected light is then attenuated while moving back upward to the water surface. The difficulty of satellite-derived bathymetry (SDB) arises from the fact that the measured radiance is a result of a complex association of physical elements, mainly the optical properties of the water, bottom reflectance, and depth. In this research, we attempt to apply physics-based algorithms to solve this complex problem as accurately as possible to overcome the limitation of having only a few known values from a multispectral sensor. Major analysis components are atmospheric correction, the estimation of water optical properties from optically deep water, and the optimization of bottom reflectance as well as the water depth. Specular reflection of the sky radiance from the water surface is modeled in addition to the typical atmospheric correction. The physical modeling of optically dominant components such as dissolved organic matter, phytoplankton, and suspended particulates allows the inversion of water attenuation coefficients from optically deep pixels. The atmospheric correction and water attenuation results are used in the ocean optical reflectance equation to solve for the bottom reflectance and water depth. At each stage of the solution, physics-based models and a physically valid, constrained Levenberg–Marquardt numerical optimization technique are used. The physics-based algorithm is applied to Landsat Operational Land Imager (OLI) imagery over the shallow coastal zone of Guam, Key West, and Puerto Rico. The SDB depths are compared to airborne lidar depths, and the root mean squared error (RMSE) is mostly less than 2 m over water as deep as 30 m. As the initial choice of bottom reflectance is critical, along with the bottom reflectance library, we describe a pure bottom unmixing method based on eigenvector analysis to estimate unknown site-specific bottom reflectance.

Remote Sensing↗
USGS2024

An ensemble mean method for remote sensing of actual evapotranspiration to estimate water budget response across a restoration landscape

Estimates of actual evapotranspiration (ETa) are valuable for effective monitoring and management of water resources. In areas that lack ground-based monitoring networks, remote sensing allows for accurate and consistent estimates of ETa across a broad scale—though each algorithm has limitations (i.e., ground-based validation, temporal consistency, spatial resolution). We developed an ensemble mean ETa (EMET) product to incorporate advancements and reduce uncertainty among algorithms (e.g., energy-balance, optical-only), which we use to estimate vegetative water use in response to restoration practices being implemented on the ground using management interventions (i.e., fencing pastures, erosion control structures) on a private ranch in Baja California Sur, Mexico. This paper describes the development of a monthly EMET product, the assessment of changes using EMET over time and across multiple land use/land cover types, and the evaluation of differences in vegetation and water distribution between watersheds treated by restoration and their controls. We found that in the absence of a ground-based monitoring network, the EMET product is more robust than using a single ETa data product and can augment the efficacy of ETa-based studies. We then found increased ETa within the restored watershed when compared to the control sites, which we attribute to increased plant water availability.

Baja California Sur↗
USGS2025

Advancements in satellite observations of inland and coastal waters: Building towards a global validation network

The use of satellite-based remote sensing imagery for water quality monitoring of inland and coastal waters has become widespread over the last few decades, with the expansion of, and investment in, operational Earth-observing missions. Satellite-based sensors are uniquely suited to provide synoptic, system-wide water quality parameter estimates that supplement traditional field-based sampling methods. The remote sensing of water quality parameter estimates is particularly valuable in systems with high temporal and spatial variability, as well as in areas that are difficult to access, or where agencies lack funding for routine monitoring. However, optically complex inland and coastal waters pose additional challenges for developing robust remote sensing retrieval models for optical properties and water quality parameters. One of the biggest challenges is collecting high quality field measurements that are used to calibrate and validate the retrieval algorithms. Here, we present the current status of satellite missions, field methods that include instruments used and commonly measured parameters, and repositories of historical field data that are relevant to inland and coastal water studies. We then present data requirements for model validation and highlight gaps in validation coverage. Finally, we provide suggestions for future field campaigns to improve coordination with remote sensing data collection and to ensure that field data is well suited for use in model or algorithm development.

Remote Sensing↗
USGS2015

Categorisation of northern California rainfall for periods with and without a radar brightband using stable isotopes and a novel automated precipitation collector

During landfall of extratropical cyclones between 2005 and 2011, nearly 1400 precipitation samples were collected at intervals of 30-min time resolution with novel automated collectors at four NOAA sites in northern California [Alta (ATA), Bodega Bay (BBY), Cazadero (CZD) and Shasta Dam (STD)] during 43 events. Substantial decreases were commonly followed hours later by substantial increases in hydrogen isotopic composition ( δ 2 H VSMOW where VSMOW is Vienna Standard Mean Ocean Water) and oxygen isotopic composition ( δ 18 O VSMOW ) of precipitation. These variations likely occur as pre-cold frontal precipitation generation transitions from marine vapour masses having low rainout to cold cloud layers having much higher rainout (with concomitant brightband signatures measured by an S-band profiling radar and lower δ 2 H VSMOW values of precipitation), and finally to shallower, warmer precipitating clouds having lower rainout (with non-brightband signatures and higher δ 2 H VSMOW values of precipitation), in accord with ‘seeder–feeder’ precipitation. Of 82 intervals identified, a remarkable 100.5 ‰ decrease in δ 2 H VSMOW value was observed for a 21 January 2010 event at BBY. Of the 61 intervals identified with increases in δ 2 H VSMOW values as precipitation transitioned to shallower, warmer clouds having substantially less rainout (the feeder part of the seeder–feeder mechanism), a remarkable increase in δ 2 H VSMOW value of precipitation of 82.3 ‰ was observed for a 10 February 2007 event at CZD. All CZD and ATA events having δ 2 H VSMOW values of precipitation below −105 ‰ were atmospheric rivers (ARs), and of the 13 events having δ 2 H VSMOW values of precipitation below −80 ‰, 77 % were ARs. Cloud echo-top heights (a proxy for atmospheric temperature) were available for 23 events. The mean echo-top height is greater for higher rainout periods than that for lower rainout periods in 22 of the 23 events. The lowest δ 2 H VSMOW of precipitation of 28 CZD events was −137.9 ‰ on 16 February 2009 during an AR with cold precipitating clouds and very high rainout with tops >6.5 km altitude. An altitude effect of −2.5 ‰ per 100 m was measured from BBY and CZD δ 2 H VSMOW data and of −1.8 ‰ per 100 m for CZD and ATA δ 2 H VSMOW data. We present a new approach to categorise rainfall intervals using δ 2 H VSMOW values of precipitation and rainfall rates. We term this approach the algorithmic-isotopic categorisation of rainfall, and we were able to identify higher rainout and/or lower rainout periods during all events in this study. We conclude that algorithmic-isotopic categorisation of rainfall can enable users to distinguish between tropospheric vapour masses having relatively high rainout (typically with brightband rain and that commonly are ARs) and vapour masses having lower rainout (commonly with non-brightband rain).

California↗
USGS2015

A comparative study of surface waves inversion techniques at strong motion recording sites in Greece

Surface wave method was used for the estimation of Vs vs depth profile at 10 strong motion stations in Greece. The dispersion data were obtained by SASW method, utilizing a pair of electromechanical harmonic-wave source (shakers) or a random source (drop weight). In this study, three inversion techniques were used a) a recently proposed Simplified Inversion Method (SIM), b) an inversion technique based on a neighborhood algorithm (NA) which allows the incorporation of a priori information regarding the subsurface structure parameters, and c) Occam's inversion algorithm. For each site constant value of Poisson's ratio was assumed (ν=0.4) since the objective of the current study is the comparison of the three inversion schemes regardless the uncertainties resulting due to the lack of geotechnical data. A penalty function was introduced to quantify the deviations of the derived Vs profiles. The Vs models are compared as of Vs(z), Vs30 and EC8 soil category, in order to show the insignificance of the existing variations. The comparison results showed that the average variation of SIM profiles is 9% and 4.9% comparing with NA and Occam's profiles respectively whilst the average difference of Vs30 values obtained from SIM is 7.4% and 5.0% compared with NA and Occam's.

Conference Paper↗
USGS2009

Detection and characterization of benthic filamentous algal stands ( Cladophora sp.) on rocky substrata using a high-frequency echosounder

A high-frequency echosounder was used to detect and characterize percent cover and stand height of the benthic filamentous green alga Cladophora sp. on rocky substratum of the Laurentian Great Lakes. Comparisons between in situ observations and estimates of the algal stand characteristics (percent cover, stand height) derived from the acoustic data show good agreement for algal stands that exceeded the height threshold for detection by acoustics (~7.5 cm). Backscatter intensity and volume scattering strength were unable to provide any predictive power for estimating algal biomass. A comparative analysis between the only current commercial software (EcoSAV™) and an alternate method using a graphical user interface (GUI) written in MATLAB® confirmed previous findings that EcoSAV functions poorly in conditions where the substrate is uneven and bottom depth changes rapidly. The GUI method uses a signal processing algorithm similar to that of EcoSAV but bases bottom depth classification and algal stand height classification on adjustable thresholds that can be visualized by a trained analyst. This study documents the successful characterization of nuisance quantities of filamentous algae on hard substrate using an acoustic system and demonstrates the potential to significantly increase the efficiency of collecting information on the distribution of nuisance macroalgae. This study also highlights the need for further development of more flexible classification algorithms that can be used in a variety of aquatic ecosystems.

Great Lakes↗
USGS2026

Origins, evolutions, and future directions of Landsat science products for advancing global inland water and coastal ocean observations

In April 2020, the U.S. Geological Survey (USGS) Earth Resources Observation and Science (EROS) Center introduced a Level 2 provisional Aquatic Reflectance (AR) product for the Landsat 8 Operational Land Imager (OLI), marking the initial phase in developing a standardized global product for Landsat-derived surface water measurements. The goal of USGS EROS aquatic product research and development is to prepare for an operational processing architecture for Landsat Collection 3 in the late 2020s that will enable use of quality-controlled data for emerging Landsat aquatic science applications. To achieve this, we released a subset of the Landsat 8/9 provisional AR products (Crawford et al., 2025, https://doi.org/10.5066/P14MBBRM) and examined its general performance through the Science Algorithms to Operations (SATO) framework alongside quantitative assessment using community made inland water data records (GLObal Reflectance community dataset for Imaging and optical sensing of Aquatic environments, GLORIA) and radiometric coastal validation platforms (NASA’s Ocean Color component of the Aerosol Robotic Network, AERONET-OC). Variability within the validation datasets indicate that the performance of the Landsat 8/9 provisional AR retrieval is highly context-dependent; errors are minimal in optically simple waters (e.g., clear to moderately turbid coastal waters) but increase considerably in optically complex waters where factors such as elevated levels of turbidity, chlorophyll (Chl a ) concentrations, or colored dissolved organic matter (CDOM) dominate the water column. Additionally, this paper examines key algorithmic considerations for atmospheric correction, highlighting factors that influence accuracy, scalability, and computational efficiency necessary for collection processing in the operational Landsat Product Generation System (LPGS). This paper is intended to communicate with aquatic scientists, satellite oceanographers, and the broader Earth observation community on the origins, requirements, challenges, successes, and future objectives for operationalizing global AR data products for Landsat satellite missions.

Earth System Science Data↗
USGS2025

Review and synthesis of the applications of machine learning to coalbed methane recovery

Over the last 30 years, a substantial literature has evolved on the use of machine learning (ML) to assess, predict, and improve the efficiency of coalbed methane (CBM) recovery. In the United States, the production of CBM declined as shale gas production matured, but CBM continues to be an important energy resource in other parts of the world. ML applications that have the potential to improve CBM reservoir management and production forecasts, and to increase exploration and operational efficiency, are still of significant interest. The integration of geostatistical techniques into the CBM ML applications has been largely absent but represents an opportunity for improvement. The literature demonstrates the widespread interest in, and applicability of, ML algorithms applied to CBM problems, and that they continue to result in improvements in predictive performance. However, (1) much of the research is more academic than operational, (2) many results are based on simulations, or small or proprietary datasets, (3) ML performance information can be inconsistent and sometimes entirely omitted, (4) most methodologies are unique to the specific CBM situation and likely not generalizable, (5) no standard data repositories are available to directly compare the performance of competing algorithms, and (6) the spatial component is often omitted. Finally, relatively new ML protocols involving causality analysis and reinforced learning, as well as hybrid workflows combining both supervised and unsupervised learning, are anticipated to dominate the future investigations. Integration of geostatistical and geospatial analysis with ML should enhance performance.

Book chapter↗