Geology ReportsSearch

SEARCH · Geology Reports

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 127 records · Page 7Linked to original sources

The automated reference toolset: A soil-geomorphic ecological potential matching algorithm

Ecological inventory and monitoring data need referential context for interpretation. Identification of appropriate reference areas of similar ecological potential for site comparison is demonstrated using a newly developed automated reference toolset (ART). Foundational to identification of reference areas was a soil map of particle size in the control section (PSCS), a theme in US Soil Taxonomy. A 30-m resolution PSCS map of the Colorado Plateau (366,000 km 2 ) was created by interpolating ∼5000 field soil observations using a random forest model and a suite of raster environmental spatial layers representing topography, climate, general ecological community, and satellite imagery ratios. The PSCS map had overall out of bag accuracy of 61.8% (Kappa of 0.54, p < 0.0001), and an independent validation accuracy of 93.2% at a set of 356 field plots along the southern edge of Canyonlands National Park, Utah. The ART process was also tested at these plots, and matched plots with the same ecological sites (ESs) 67% of the time where sites fell within 2-km buffers of each other. These results show that the PSCS and ART have strong application for ecological monitoring and sampling design, as well as assessing impacts of disturbance and land management action using an ecological potential framework. Results also demonstrate that PSCS could be a key mapping layer for the USDA-NRCS provisional ES development initiative.

Soil Science Society of America Journal

Accuracy of EAARL lidar ground elevations using a bare-earth algorithm in marsh and beach grasses on the Chandeleur Islands, Louisiana

The NASA Experimental Advanced Airborne Lidar (EAARL) is an airborne lidar (light detection and ranging) instrument designed to map coastal topography and bathymetry. The EAARL system has the capability to capture each laser-pulse return over a large signal range and can digitize the full waveform of the backscattered energy. Because of this ability to capture the full waveform, the EAARL system can map features such as coral reefs, beaches, coastal vegetation, and trees, where extreme variations in the laser backscatter are caused by different physical and optical characteristics. Post-processing of the EAARL data is accomplished using the Airborne Lidar Processing System (ALPS) (Nayegandhi and others, 2009). In ALPS, the waveform of the lidar is analyzed and split into first and last returns. The 'first returns' are indicative of vegetation-canopy height, or bare ground in the absence of vegetation, whereas 'last returns' typically represent 'bare-earth' elevations under vegetation. To test the accuracy of the first-return and bare-earth EAARL data, topographic and vegetation height surveys were conducted in the Chandeleur Islands, concurrent with an EAARL lidar survey and an aerial oblique-photographic survey from September 20 to 27, 2006. The Chandeleur Islands are a north-south-oriented chain of low-lying islands located approximately 100 kilometers east of the city of New Orleans, Louisiana. The islands are narrow north-south strips of land with marsh on the landward (west sides) and sandy beaches on their gulfward (east sides). Prior to Hurricane Katrina, which made landfall at Buras, Louisiana, as a Category 3 storm on August 29, 2005, prominent, 3- to 4-meter-high sand dunes were present in the northern Chandeleurs. The storm removed them, leaving post-storm island elevations of generally less than 2 meters above 0.0 NAVD88. This report is part of a study of the impact of Hurricane Katrina on the Chandeleur Islands using pre-storm and post-storm lidar surveys to detect morphological changes. The islands lost over 80 percent of their land area during Hurricane Katrina, and in the first 2 years following Katrina, many of the islands experienced continued shoreline retreat (Sallenger and others, 2007). In addition to land-area losses, the loss of dunes made the islands increasingly vulnerable to future storm impacts. The U.S. Geological Survey, along with partners in the Louisiana Department of Natural Resources and the U.S. Army Corps of Engineers, continues to monitor changes in shoreline position, land area, and elevation in the Chandeleur Islands.

Open-File Report

Improving temporal frequency of Landsat surface temperature products using the gap-filling algorithm

Remotely sensed surface temperature (ST) has been widely used to monitor and assess landscape thermal conditions, hydrologic modeling, and surface energy balance. Landsat thermal sensors have continuously measured the Earth surface thermal radiance since August 1982. The thermal radiance measurements are atmospherically compensated and converted to Landsat STs and delivered as part of the U.S. Geological Survey Landsat Collection 1 U.S. Analysis Ready Data; however, the low satellite revisit cycles combined with the presence of clouds and cloud shadows reduce the number of valid retrievals. This reduction can limit the ability to monitor annual or seasonal variations in the surface thermal budget. These factors reduce the ability to use the temperature data to fit time series for historical trend analysis to match background climate variations. In this study, we implemented an approach that uses linear harmonic least absolute shrinkage and selection operator regression models to fill gaps because of clouds, shadows, and coarse temporal resolution. The gap-filled data provide increased temporal density of Landsat ST records. The gap-filled Landsat ST, therefore, can allow for an improved monitoring of annual, seasonal, or even monthly landscape thermal conditions.

Georgia

Application of Gauss algorithm and Monte Carlo simulation to the identification of aquifer parameters

The Gauss optimization technique can be used to identify the parameters of a model of a groundwater system for which the parameter identification problem is formulated as a least squares comparison between the response of the prototype and the response of the model. Unavoidable uncertainty in the true stress on the prototype and in the true response of the prototype to that stress will introduce errors into the parameter identification problem. A method for evaluating errors in the predictions of future water levels due to errors in recharge estimates was demonstrated. The method involves a Monte Carlo simulation of the parameter identification problem and of the prediction problem. The steps in the method are: (1) to prescribe the distribution of the recharge estimates; (2) to use this distribution to generate random sets of recharge estimates; (3) to use the Gauss optimization technique to identify the corresponding set of parameter estimates for each set of recharge estimates; (4) to make the corresponding set of hydraulic head predictions for each set of parameter estimates; and (5) to examine the distribution of hydraulic head predictions and to draw appropriate conclusions. Similarly, the method can be used independently or simultaneously to estimate the effect on hydraulic head predictions of errors in the measured water levels that are used in the parameter identification problem. The fit of the model to the data that are used to identify parameters is not a good indicator of these errors. A Monte Carlo simulation of the parameter identification problem can be used, however, to evaluate the effects on water level predictions of errors in the recharge (and pumpage) data used in the parameter identification problem. (Lantz-PTT)

Open-File Report

Nominal 30-m cropland extent map of continental Africa by integrating pixel-based and object-based algorithms using Sentinel-2 and Landsat-8 Data on Google Earth Engine

A satellite-derived cropland extent map at high spatial resolution (30-m or better) is a must for food and water security analysis. Precise and accurate global cropland extent maps, indicating cropland and non-cropland areas, are starting points to develop higher-level products such as crop watering methods (irrigated or rainfed), cropping intensities (e.g., single, double, or continuous cropping), crop types, cropland fallows, as well as for assessment of cropland productivity (productivity per unit of land), and crop water productivity (productivity per unit of water). Uncertainties associated with the cropland extent map have cascading effects on all higher-level cropland products. However, precise and accurate cropland extent maps at high spatial resolution over large areas (e.g., continents or the globe) are challenging to produce due to the small-holder dominant agricultural systems like those found in most of Africa and Asia. Cloud-based geospatial computing platforms and multi-date, multi-sensor satellite image inventories on Google Earth Engine offer opportunities for mapping croplands with precision and accuracy over large areas that satisfy the requirements of broad range of applications. Such maps are expected to provide highly significant improvements compared to existing products, which tend to be coarser in resolution, and often fail to capture fragmented small-holder farms especially in regions with high dynamic change within and across years. To overcome these limitations, in this research we present an approach for cropland extent mapping at high spatial resolution (30-m or better) using the 10-day, 10 to 20-m, Sentinel-2 data in combination with 16-day, 30-m, Landsat-8 data on Google Earth Engine (GEE). First, nominal 30-m resolution satellite imagery composites were created from 36,924 scenes of Sentinel-2 and Landsat-8 images for the entire African continent in 2015–2016.

Remote Sensing

Investigating hydrologic alteration under historical and future scenarios in the Mobile River and Perdido River basins using the cubist algorithm

This study investigates the impacts of human activities and climate variability on hydrologic alterations in the Mobile River and Perdido River Basins of Alabama. The research uses a machine learning approach, specifically cubist models, to quantify and predict changes in flow duration curves (FDCs) under both historical (1980–2009) and future climate scenarios. Future climate projections include the Representative Concentration Pathways (RCP 4.5 and RCP 8.5) and the Shared Socioeconomic Pathways (SSP2 4.5 and SSP5 8.5), evaluated for two future periods: 1980–2069 and 1980–2099. The models incorporate a wide range of covariates, including basin geomorphology, aquifer characteristics, land cover, water storage, environmental factors, solar radiation, census data, and water use data. Under the baseline period (1980–2009), most level 12 hydrologic unit codes (HUC12s) in both basins showed alterations, with substantial differences observed between pre- and post-alteration FDCs. The model performance varied, with a Nash–Sutcliffe Efficiency between 0.91 and 0.95 for testing and between 0.98 and 0.99 for training during the baseline period. Future projections under the RCP 4.5 and RCP 8.5 scenarios generally differed significantly from baseline conditions across all flow regimes ( p < 0.05). In contrast, SSP2 4.5 showed comparatively limited statistical significance, while SSP5 8.5 exhibited significant departures from baseline conditions across all flow regimes, reflecting the greater influence of high-emissions climate forcing on projected hydrologic alterations. Overall, the RCP scenarios projected more widespread statistically significant changes than the corresponding SSP scenarios at the same forcing level, particularly when comparing RCP4.5 with SSP2-4.5, while both RCP8.5 and SSP5-8.5 consistently indicated greater hydrologic alterations than their moderate-emissions counterparts. These findings highlight the importance of considering different flow regimes when assessing the impacts of climate variability on streamflow. This study contributes to the understanding of hydrologic alterations in the Mobile River and Perdido River Basins, providing insights for water resource management and ecological conservation efforts in the region.

Alabama

Seasonal comparisons of sea ice concentration estimates derived from SSM/I, OKEAN, and RADARSAT data

The Special Sensor Microwave Imager (SSM/I) microwave satellite radiometer and its predecessor SMMR are primary sources of information for global sea ice and climate studies. However, comparisons of SSM/I, Landsat, AVHRR, and ERS-1 synthetic aperture radar (SAR) have shown substantial seasonal and regional differences in their estimates of sea ice concentration. To evaluate these differences, we compared SSM/I estimates of sea ice coverage derived with the NASA Team and Bootstrap algorithms to estimates made using RADARSAT, and OKEAN-01 satellite sensor data. The study area included the Barents Sea, Kara Sea, Laptev Sea, and adjacent parts of the Arctic Ocean, during October 1995 through October 1999. Ice concentration estimates from spatially and temporally near-coincident imagery were calculated using independent algorithms for each sensor type. The OKEAN algorithm implemented the satellite's two-channel active (radar) and passive microwave data in a linear mixture model based on the measured values of brightness temperature and radar backscatter. The RADARSAT algorithm utilized a segmentation approach of the measured radar backscatter, and the SSM/I ice concentrations were derived at National Snow and Ice Data Center (NSIDC) using the NASA Team and Bootstrap algorithms. Seasonal and monthly differences between SSM/I, OKEAN, and RADARSAT ice concentrations were calculated and compared. Overall, total sea ice concentration estimates derived independently from near-coincident RADARSAT, OKEAN-01, and SSM/I satellite imagery demonstrated mean differences of less than 5.5% (S.D.<9.5%) during the winter period. Differences between the SSM/I NASA Team and the SSM/I Bootstrap concentrations were no more than 3.1% (S.D.<5.4%) during this period. RADARSAT and OKEAN-01 data both yielded higher total ice concentrations than the NASA Team and the Bootstrap algorithms. The Bootstrap algorithm yielded higher total ice concentrations than the NASA Team algorithm. Total ice concentrations derived from OKEAN-01 and SSM/I satellite imagery were highly correlated during winter, spring, and fall, with mean differences of less than 8.1% (S.D.<15%) for the NASA Team algorithm, and less than 2.8% (S.D.<13.8%) for the Bootstrap algorithm. Respective differences between SSM/I NASA Team and SSM/I Bootstrap total concentrations were less than 5.3% (S.D.<6.9%). Monthly mean differences between SSM/I and OKEAN differed annually by less than 6%, with smaller differences primarily in winter. The NASA Team and Bootstrap algorithms underestimated the total sea ice concentrations relative to the RADARSAT ScanSAR no more than 3.0% (S.D.<9%) and 1.2% (S.D.<7.5%) during cold months, and no more than 12% and 7% during summer, respectively. ScanSAR tended to estimate higher ice concentrations for ice concentrations greater than 50%, when compared to SSM/I during all months. ScanSAR underestimated total sea ice concentration by 2% compared to the OKEAN-01 algorithm during cold months, and gave an overestimation by 2% during spring and summer months. Total NASA Team and Bootstrap sea ice concentration estimates derived from coincident SSM/I and OKEAN-01 data demonstrated mean differences of no more than 5.3% (S.D.<7%), 3.1% (S.D.<5.5%), 2.0% (S.D.<5.5%), and 7.3% (S.D.<10%) for fall, winter, spring, and summer periods, respectively. Large disagreements were observed between the OKEAN and NASA Team results in spring and summer for estimates of the first-year (FY) and multiyear (MY) age classes. The OKEAN-01 algorithm and data tended to estimate, on average, lower concentrations of young or FY ice and higher concentrations of total and MY ice for all months and seasons. Our results contribute to the growing body of documentation about the levels of disparity obtained when seasonal sea ice concentrations are estimated using various types of satellite data and algorithms.

Remote Sensing of Environment

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

Exploring the factors controlling the error characteristics of the Surface Water and Ocean Topography mission discharge estimates

The Surface Water and Ocean Topography (SWOT) satellite mission will measure river width, water surface elevation, and slope for rivers wider than 50-100 m. SWOT observations will enable estimation of river discharge by using simple flow laws such as the Manning-Strickler equation, complementing in-situ streamgages. Several discharge inversion algorithms designed to compute unobserved flow law parameters ( e.g. friction coefficient, bathymetry) have been proposed, but to date, a systematic assessment of factors controlling algorithm performance has not been conducted. Here, we assess the performance of the five algorithms that are expected to be used in the construction of the SWOT product. To perform this assessment, we used synthetic SWOT observations created with hydraulic model output corrupted with SWOT-like error. Prior information provided to the algorithms was purposefully limited to an estimate of mean annual flow (MAF), designed to produce a “worst case” benchmark. Prior MAF error was an important control on algorithm performance, but discharge estimates produced by the algorithms are less biased than the MAF; thus, the discharge algorithms improve on the prior. We show for the first time that accuracy and frequency of remote sensing observations are less important than prior bias, hydraulic variability among reaches, and flow law accuracy in governing discharge algorithm performance. The discharge errors and error sensitivities reported herein are a bounding benchmark, representing worst possible expected errors and error sensitivities. This study lays the groundwork to develop predictive power of algorithm performance, and thus map the global distribution of worst-case SWOT discharge accuracy.

Water Resources Research

A two-dimensional, reach-scale implementation of space-time image velocimetry (STIV) and comparison to particle image velocimetry (PIV)

Image-based algorithms have become a powerful tool for estimating flow velocities in rivers. In this study, we generalize the space-time image velocimetry (STIV) framework for reach-scale application rather than along a cross section. The new algorithm provides information on both the magnitude and orientation of velocity vectors, and we refer to the algorithm as two-dimensional STIV, or 2D-STIV. The workflow involves setting up a grid, using centreline tangent vectors as initial estimates of flow direction, and then extracting space-time images (STIs) along search lines radiating from each grid node. The autocorrelation function is used to infer the inclination of streak lines present in STIs, which represents the advection of water surface features. Information on flow direction is obtained by evaluating various candidate search lines and identifying that which yields the highest velocity. This search can be performed exhaustively or via optimization. We applied the new 2D-STIV algorithm to three test cases, one simulated data set and two natural channels, and compared image-derived velocities to modelled or measured values. We also applied two established particle image velocimetry (PIV) algorithms to the same data sets. 2D-STIV performed as well as the two PIV algorithms for simulated images. For a natural river with distinct water surface features, 2D-STIV was effective for much of the channel but also led to a more patchy, irregular velocity field than the two PIV algorithms. For a site lacking obvious surface features, exhaustive 2D-STIV led to velocity estimates uncorrelated with field data while the optimization-based version produced erratic flow directions. 2D-STIV also required greater image sequence durations, higher frame rates, and generally longer computational run times. Overall, ensemble PIV was the most reliable algorithm.

Earth Surface Processes and Landforms

Fish body midline segmentation using binary search

Body and caudal fin locomotion is ubiquitous in aquatic vertebrates, and kinematic models describing it are used in robotics, biomechanics and fisheries research. This paper presents a new algorithm to translate continuous body midlines of fish into a series of interconnected segments by identifying favorable joint positions along the body. The algorithm employs binary search to generate parsimonious kinematic models, aiming at minimizing the number of segments yet keeping approximation error below a user-defined threshold. To achieve this, the algorithm maximizes the length of each segment by determining the most distal joint position through repetitive shrinking of the search space. Theoretical and empirical analysis using two different datasets show that the binary search algorithm is substantially faster when compared to segment growing algorithm, which employs linear search to generate its models. There is four-fold improvement in computation time when generating models with less than 10 segments, which are typically sufficient to describe fish and fish-inspired robot movements. Furthermore, the multi-segment models generated by the binary search algorithm matched the ground truth models obtained through dynamic programming in over 97% of cases, and on average, contained one fewer segment than those produced by the Ramer–Douglas–Peucker algorithm, which is widely used in curvature simplification tasks. Our findings suggest that the binary search algorithm provides a computationally efficient approach for generating compact kinematic models and may facilitate the analysis of large datasets with high temporal and spatial resolution.

Computers and Electronics in Agriculture

Simplification of polylines by segment collapse: Minimizing areal displacement while preserving area

This paper reports on a new Area Preserving Segment Collapse (APSC) algorithm for simplifying polygonal boundaries while preserving polygonal area at simplified target scales and minimizing areal displacement. A general segment collapse algorithm is defined by iteratively collapsing segments to Steiner points in priority order, guided by placement and displacement functions. The algorithm is specified by defining functions that minimize areal displacement under the constraint that the areas of adjoining polygons are preserved exactly. Self-intersections can be avoided by testing for intersections with two new line segments associated with each segment collapse operation. The paper demonstrates simplification results for a sample of 10 lakes formed from alpine, Karst, glacial and arid desert processes as well as artificial dams. APSC results are compared with three other simplification routines and evaluated for area preservation, linear and areal displacement, complexity and introduction of boundary self-intersections. Results confirm that the APSC algorithm preserves area exactly and indicate that it outperforms the other tested algorithms for minimizing areal displacement while producing reasonably low measures of linear displacement. Self-intersections can occur more commonly with APSC than other algorithms but are avoided with the proposed topology check. The APSC algorithm additionally preserves polygon complexity better than other tested algorithms.

International Journal of Cartography