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.

SourceResearch connectionFrom yearThrough year

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

At least 1,297 records · Page 72Linked to original sources

USGS2021

Changes in forest connectivity from beech bark disease in Pictured Rocks National Lakeshore in the Upper Peninsula of Michigan

Within the forests of Pictured Rocks National Lakeshore, biologists are trying to understand the effects beech bark disease has on wildlife species, especially species that need forest connectivity to thrive. This project used aerial imagery collected in 2005, shortly after beech bark disease infestation, and satellite imagery from 2018. The 2018 imagery represents present day conditions and was used to locate forest canopy gaps through object-based image analysis. Forest canopy gaps were identified using the multiresolution segmentation algorithm within Trimble’s eCognition software. A time change analysis was completed to understand how the forest canopy had changed from 2005 to 2018. The analysis showed areas that had maintained forest canopy, maintained a forest canopy gap, created a new canopy gap (closed forest canopy in 2005 but open canopy gap in 2018), or created new forest canopy (open canopy gap in 2005 but closed forest canopy in 2018). There were 9,127 acres of forest canopy lost, and 72.8 percent of that lost canopy occurred in a forest type where Fagus grandifolia Ehrh. (American beech) is a common tree species. The datasets developed through this project can enhance knowledge of where canopy gaps exist and help place focus on certain areas for wildlife studies. In addition, these datasets can be used in future studies to monitor the health of the forest and conduct additional change analyses.

Michigan↗

USGS2023

Simulation of regional groundwater flow and advective transport of per- and polyfluoroalkyl substances, Joint Base McGuire-Dix-Lakehurst and vicinity, New Jersey, 2018

A three-dimensional numerical model of groundwater flow was developed and calibrated for the unconsolidated New Jersey Coastal Plain aquifers underlying Joint Base McGuire-Dix-Lakehurst (JBMDL) and vicinity, New Jersey, to evaluate groundwater flow pathways of per- and polyfluoroalkyl substances (PFAS) contamination associated with use of aqueous film forming foam (AFFF) at the base. The regional subsurface flow model spans an area of approximately 518 square miles around JBMDL and is based on a previously developed hydrogeologic framework of the area. Steady-state flow in the unconsolidated aquifers was simulated using the MODFLOW 6 groundwater flow model, which is able to account for hydrostratigraphic pinchouts and discontinuities in the Coastal Plain aquifers underlying JBMDL. To account for local patterns of fluid flow driving advective subsurface migration of PFAS, the grid was refined using quadtree meshes spanning 21 areas where historical AFFF use was identified, five off-site reconnaissance areas identified by AFCEC as areas in which the occurrence of PFAS is most likely to pose a potential danger to local drinking water supplies, and along streams that behave as drains in the base-flow-dominated Coastal Plain. Following grid refinement, four physical processes known to govern subsurface flow were introduced to the model. These included effective precipitation recharge, discharge to streams and stream-connected wetlands, regional inflows and outflows along the model bottom, and withdrawals from wells, each of which were incorporated into the model as either external or internal boundary conditions. To account for effective precipitation recharge, a specified-flow boundary was assigned along the top of the model. Similarly, regional flows predicted using the modified U.S Geological Survey’s New Jersey Coastal Plain Regional Aquifer System Analysis model were treated as specified-flow boundary conditions along the bottom of the model. Base-flow losses were treated as drains along streams delineated using a 10-foot LiDAR dataset. Drains were also assigned to cells falling within stream-connected National Hydrologic Database wetlands. Finally, well-pumpage data mined from the New Jersey Water Transfer database were added to the model to account for extraction of groundwater through pumping from industrial-supply and drinking-water-supply wells. Along model edges established at groundwater divides, where the net flux of water across the boundary is equal to zero, natural no-flow boundary conditions were imposed. The refined flow model was calibrated using the parameter-estimation (PEST) program, which adjusts model parameters by performing a gradient search over the sum-of-squared-error objective function until the parameter set that produces simulated water levels and base flows most closely matches 544 water levels and 20 estimated base flows and closely adheres to initial parameter estimates. Based on the analysis of calibration residuals, the model did not appear to be affected by significant model structural error. The MODPATH particle-tracking algorithm was used to estimate advective transport paths of PFAS in the vicinity of JBMDL. Forward tracking was used to determine paths of PFAS away from AFFF source areas to streams, wetlands, pumping wells, and geographic areas that PFAS may contaminate. Additionally, reverse tracking was used to determine particle pathlines away from off-site PFAS reconnaissance areas, or areas within which all sources of PFAS might be advectively transported into subsurface drinking-water supplies, to locations at land surface that may indicate a source of PFAS. The coupled and calibrated groundwater flow and particle-tracking transport model provide valuable tools for predicting the relative extent of PFAS contamination from onsite legacy source areas. The calibrated model also provides measures of water-level and base-flow observation influence that can help guide future data-collection efforts related to groundwater and surface water sampling for PFAS.

New Jersey↗

USGS2023

Crustal thickness and the VP/VS ratio within the Arabia Plate from P-wave receiver functions at 154 broadband seismic stations

As part of a joint Saudi Geological Survey (SGS) and United States Geological Survey project, we analyzed P-wave receiver functions from seismic stations covering most of the Kingdom of Saudi Arabia to map the thickness of the crust across the Arabia Plate. We present an update of crustal thickness estimates and fill in gaps for the western Arabian Shield and the rifted margin at the Red Sea (the coastal plain), as well as the eastern Arabian Platform. We applied a conventional H-k stacking algorithm and included careful attention to stacking weights, two forms of sedimentary corrections for stations located on the Arabian Platform, and additional processing for noisy stations. We obtained useful results at 154 stations from 898 teleseismic events over a 2-year period from 1995–1997 (for non-SGS stations) and a 6-year period from 2008–2014 (for SGS stations). Average crustal thickness (that is, depth to the Mohorovičić discontinuity [Moho] below the surface) beneath the Red Sea coastal plain (the rift margin) is 29 kilometers (km), beneath the volcanic fields (known in Arabic as harra [plural] or harrat [singular]) is 35 km, beneath the Arabian Shield (excluding harrats) is 37 km, and beneath the Arabian Platform is 38 km. Crustal thinning appears not to extend east of the rift escarpment, suggesting uniform extension that is no broader at depth than at the surface. In contrast to some previous interpretations that the Arabian Platform crust is thicker than that of the Arabian Shield, we find no statistically significant difference between their whole crustal thicknesses. However, the average sub-sedimentary crustal thickness (that is, the crystalline crust) for stations on the Arabian Platform is 34 km, 3 km thinner than the crust of the Arabian Shield. Individual station P-wave (pressure) velocity and S-wave (shear) velocity ratios (VP/VS) are highly variable for the Arabia Plate, ranging from 1.60 to 1.97 and averaging 1.75, with a standard deviation of 0.07. There are no statistically significant differences between VP/VS ratios of the different geologic regions of Saudi Arabia. Similar VP/ VS ratios, coupled with similar crustal thicknesses for harrats and the Arabian Shield, indicate that Cenozoic magmatism has contributed negligibly to crustal growth.

Open-File Report↗

USGS2024

Database and time series of nearshore waves along the Alaskan coast from the United States-Canada border to the Bering Sea

Alaska’s Arctic coast has some of the highest coastal erosion rates in the world, primarily driven by permafrost thaw and increasing wave energy. In the Arctic, a warming climate is driving sea ice cover to decrease in space and time. A lack of long-term observational wave data along Alaska’s coast challenges the ability of engineers, scientists, and planners to study and address threats and effects from wave-driven erosion and flooding. To overcome the lack of available observational wave data in the nearshore in this study by the U.S. Geological Survey, waves were downscaled with the Simulating WAves Nearshore numerical wave model (SWAN) for the hindcast period of 1979 to 2019 from the United States-Canada border to the Bering Sea utilizing nine model domains. For each domain, the model was forced at the open boundary with 2,500 representative “sea states,” which are likely combinations of significant wave heights, mean wave periods, mean wave directions, and wind speeds and directions. The sea states were obtained from the European Centre for Medium-Range Weather Forecasts “ERA5” dataset for reanalysis of winds and waves using a multivariant maximum-dissimilarity algorithm. The SWAN runs created a downscaled wave database at each grid point, which was used to reconstruct the 40-year time series in the nearshore along the 5- and 10-meter isobaths at locations approximately 400 m apart and corresponding to transects spaced approximately 50 m alongshore, as developed for USGS shoreline-change assessments. Reconstructed time series were compared to observations to validate the numerical model and the downscaled wave database method and showed overall good agreements.

Alaska↗

USGS2024

Total phosphorus and suspended-sediment concentrations and loads from two main tributaries to Upper Klamath Lake, Oregon, 2014–20

Total phosphorus (TP) and suspended-sediment concentrations (SSC) and loads were computed at two U.S. Geological Survey (USGS) streamgages in the upper Klamath River Basin on the Sprague (USGS site ID 11501000) and Williamson (USGS site ID 11502500) Rivers using high temporal resolution turbidity and streamflow data to develop surrogate regression models. Regression models were updated and validated for TP at the Williamson River site, and additional data improved a prior published TP model, increasing the coefficient of determination (R 2 ) from 0.73 to 0.88. A new TP regression model was developed for the Sprague River site using 2 years of data and showed promising results with an R 2 of 0.93. Suspended-sediment concentration (SSC) surrogate models were also updated at these sites using a longer period of record than the TP models and improved characterization of sediment transport conditions at these monitoring sites. Computations of TP loads were compared to the annual loading capacity dictated by the total maximum daily load (TMDL) for Upper Klamath Lake and showed that the combined TP load of the Williamson and Sprague Rivers approaches the annual loading capacity in water years with high annual streamflow. TP loads were also compared to loads computed by the Klamath Tribes using a long-term dataset and a regression and interpolation algorithm (RIA). The comparison showed that the two methods report similar annual loads, with the surrogate regression method generally reporting lower loads than the RIA, and the RIA annual loads falling within the range of uncertainty of the surrogate regression model results. Determining the effect of habitat and stream restoration on basin-scale TP and suspended-sediment loading is challenging using the surrogate regression method at these sites given the short period of record that TP and suspended-sediment load (SSL) data are available. However, long-term analysis by the Klamath Tribes in their larger monitoring network could provide insight into the impact of restoration at smaller spatial scales compared to the basin-wide assessment produced in this study.

Oregon↗

USGS1972

Hydrograph simulation models of the Hillsborough and Alafia Rivers, Florida: a preliminary report

Mathematical (digital) models that simulate flood hydrographs from rainfall records have been developed for the following gaging stations in the Hillsborough and Alafia River basins of west-central Florida: Hillsborough River near Tampa, Alafia River at Lithia, and north Prong Alafia River near Keysville. These models, which were developed from historical streamflow and and rainfall records, are based on rainfall-runoff and unit-hydrograph procedures involving an arbitrary separation of the flood hydrograph. These models assume the flood hydrograph to be composed of only two flow components, direct (storm) runoff, and base flow. Expressions describing these two flow components are derived from streamflow and rainfall records and are combined analytically to form algorithms (models), which are programmed for processing on a digital computing system. Most Hillsborough and Alafia River flood discharges can be simulated with expected relative errors less than or equal to 30 percent and flood peaks can be simulated with average relative errors less than 15 percent. Because of the inadequate rainfall network that is used in obtaining input data for the North Prong Alafia River model, simulated peaks are frequently in error by more than 40 percent, particularly for storms having highly variable areal rainfall distribution. Simulation errors are the result of rainfall sample errors and, to a lesser extent, model inadequacy. Data errors associated with the determination of mean basin precipitation are the result of the small number and poor areal distribution of rainfall stations available for use in the study. Model inadequacy, however, is attributed to the basic underlying theory, particularly the rainfall-runoff relation. These models broaden and enhance existing water-management capabilities within these basins by allowing the establishment and implementation of programs providing for continued development in these areas. Specifically, the models serve not only as a basis for forecasting floods, but also for simulating hydrologic information needed in flood-plain mapping and delineating and evaluating alternative flood control and abatement plans.

Florida↗

USGS1976

FORTRAN computer program for seismic risk analysis

A program for seismic risk analysis is described which combines generality of application, efficiency and accuracy of operation, and the advantage of small storage requirements. The theoretical basis for the program is first reviewed, and the computational algorithms used to apply this theory are described. The information required for running the program is listed. Published attenuation functions describing the variation with earthquake magnitude and distance of expected values for various ground motion parameters are summarized for reference by the program user. Finally, suggestions for use of the program are made, an example problem is described (along with example problem input and output) and the program is listed.

Open-File Report↗

USGS1981

A stochastic streamflow model of the Platte River at Overton, Odessa, and Grand Island, Nebraska

A stochastic model is developed to simulate flows for three seasons (September through February, March and April, May through August) at Overton, Odessa, and Grand Island, Nebraska on the Platte River. The model preserves the first and second order moment properties of the historical flow series, including significant autocorrelations within each station and cross-correlations between the stations. Higher order moments are preserved by a transformation technique that allows the residuals from the model to be approximately normal. An algorithm is given for simulation of combined-stations flows, and some simulations are carried out to determine the likelihood of specified flow shortages that could be detrimental to the wildlife habitat in the reach of the Platte River between Overton and Grand Island. (USGS)

Open-File Report↗

USGS1982

A hybrid structure for the storage and manipulation of very large spatial data sets

The map data input and output problem for geographic information systems is rapidly diminishing with the increasing availability of mass digitizing, direct spatial data capture and graphics hardware based on raster technology. Although a large number of efficient raster-based algorithms exist for performing a wide variety of common tasks on these data, there are a number of procedures which are more efficiently performed in vector mode or for which raster mode equivalents of current vector-based techniques have not yet been developed. This paper presents a hybrid spatial data structure, named the ?vaster' structure, which can utilize the advantages of both raster and vector structures while potentially eliminating, or greatly reducing, the need for raster-to-vector and vector-to-raster conversion. Other advantages of the vaster structure are also discussed.

Open-File Report↗

USGS1987

Interpretation of detailed aerial gamma-ray survey, Jabal Ashirah area, southeastern Arabian Shield, Kingdom of Saudi Arabia

A detailed aerial gamma-ray spectrometric survey of the Jabal Ashirah area in the southeastern Arabian Shield has been analyzed using computer-classification algorithms. The analysis resulted in maps that show radiometric map units and gamma-ray anomalies indicating the presence of possible concentrations of potassium and uranium. The radiometric-unit map was interpreted to 'produce a simplified radiolithic map that was correlated with the mapped geology. The gamma-ray data show uranium anomalies that coincide with a tin-bearing granite, but known gold and nickel mineralization do not have any associated gamma-ray signatures.

Open-File Report↗

USGS1989

Mass-conserving method of characteristics for streamflow modeling

A robust numerical model is presented for the computation of unsteady streamflow on steep river slopes. The one-dimensional model uses the method of characteristics on a specified space-time grid to solve the Saint-Venant equations. An additional continuity equation requirement on each space-time element provides greatly improved conservation of mass over traditional implementations of the method of characteristics on a fixed grid. The space-time geometry of the problem is described in a finite element setting. Hermite interpolation of channel parameters is used to avoid numerical difficulties that may occur with steep slopes due to discontinuities in the derivatives of data such as channel top width. Manning's equation for friction slope can be modified by a factor to make the slope more appropriate for steep rivers. The standard Manning's friction slope can also be used, if preferred. The computer model is not restricted to steep slopes, and applies as well to gently sloping streams. Two numerical examples support the mathematical approach and computational algorithm. (USGS)

Open-File Report↗

USGS1994

A comparison of change detection methods using multispectral scanner data

Change detection methods were investigated as a cooperative activity between the U.S. Geological Survey and the National Bureau of Surveying and Mapping, People's Republic of China. Subtraction of band 2, band 3, normalized difference vegetation index, and tasseled cap bands 1 and 2 data from two multispectral scanner images were tested using two sites in the United States and one in the People's Republic of China. A new statistical method also was tested. Band 2 subtraction gives the best results for detecting change from vegetative cover to urban development. The statistical method identifies areas that have changed and uses a fast classification algorithm to classify the original data of the changed areas by land cover type present for each image date.

Open-File Report↗

USGS1996

Use of particle tracking to improve numerical model calibration and to analyze ground-water flow and contaminant migration, Massachusetts Military Reservation, western Cape Cod, Massachusetts

A steady-state, three-dimensional numerical model coupled with a particle-tracking algorithm was developed to simulate the complex hydrogeologic conditions affecting ground-water flow and contaminant migration in the Cape Cod aquifer beneath the Massachusetts Military Reservation, Massachusetts. The known extents of the contaminant plumes beneath the reservation were incorporated into a particle-tracking analysis to improve model calibration. Particle tracking was used to evaluate the effects of simulated changes in hydraulic properties and in simulated hydrologic boundaries such as ponds and streams. The model simulations made during the calibration process indicated that changes in simulated hydraulic properties and hydrologic boundaries resulted in small changes in the water-table and pond altitudes and in streamflows, yet had a substantial effect on model-calculated groundwater flowpaths. Therefore, the characterization of contaminant migration using a model calibrated only on the basis of ground-water heads and flux may be inaccurate. The results of model simulations for the analysis of the effects of pumping and recharge on ground-water flow and contaminant migration indicated that ground-water flowpaths were greatly affected by subtle shifts in hydraulic gradients. These changes in the ground-water-flow system in response to hydrologic stresses such as pumping and recharge can be determined by the use of a particle-tracking analysis.

Massachusetts↗

USGS1996

Catalog of earthquake hypocenters for Augustine, Redoubt, Iliamna, and Mount Spurr volcanoes, Alaska: January 1, 1991 - December 31, 1993

The Alaska Volcano Observatory (AVO), a cooperative program of the U.S. Geological Survey, the Geophysical Institute of the University of Alaska, Fairbanks, and the Alaska Division of Geological and Geophysical Surveys, has maintained a program of seismic monitoring at potentially active volcanoes in the Cook Inlet region since 1988. The principal objectives of this program include the seismic surveillance of the Cook Inlet volcanoes and the investigation of seismic processes associated with active volcanism. This catalog reflects the status and evolution of the seismic monitoring program, and presents the basic seismic data for the time interval January 1, 1991, to December 31, 1993. For an interpretation of these data the reader should refer to several recent articles on volcano related seismicity in the Cook Inlet region (e.g. Jolly and others, 1994; Power and others, 1995; and McNutt and others, 1995). A similar catalog covers the period from October 12, 1989 to December 31, 1991 (Power and others 1993). The AVO seismic monitoring program has undergone significant changes during the catalog period. The changes included 1) new seismic stations placed at Mount Spurr and Redoubt Volcano, resulting in increased earthquake detection capability and improved earthquake locations, 2) the addition of several regional stations to the seismic data acquisition system which improved location quality near the volcano and enhanced our ability to scale eruptions and measure magnitudes of the largest volcanic earthquakes, 3) installation of a new event detection algorithm XDETECT (Rogers, 1993), which offered increased data collection capabilities , 4) modifications to the earthquake location program HYPOELLIPSE (Lahr, 1989) which now allows distinct velocity models and station corrections at each volcanic center, and 5) the addition of seismic stations at Augustine and niamna volcanoes to the data acquisition/location system. The 1992 eruptions at Mount Spurr's Crater Peak vent provided the highlight of the catalog period. The crisis included three sub-plinian eruptions, which occurred on June 27, August 18, and September 16-17, 1992. The three eruptions punctuated a complex seismic sequence which included volcano-tectonic (VT) earthquakes, tremor, and both deep and shallow long period (LP) earthquakes. The seismic sequence began on August 18, 1991, with a small swarm of volcano-tectonic events beneath Crater Peak, and spread throughout the volcanic complex by November of the same year. Elevated levels of seismicity persisted at Mount Spurr beyond the catalog time period.

Alaska↗

USGS1997

A manual for a laboratory information management system (LIMS) for light stable isotopes

The reliability and accuracy of isotopic data can be improved by utilizing database software to (i) store information about samples, (ii) store the results of mass spectrometric isotope-ratio analyses of samples, (iii) calculate analytical results using standardized algorithms stored in a database, (iv) normalize stable isotopic data to international scales using isotopic reference materials, and (v) generate multi-sheet paper templates for convenient sample loading of automated mass-spectrometer sample preparation manifolds. Such a database program is presented herein. Major benefits of this system include (i) an increase in laboratory efficiency, (ii) reduction in the use of paper, (iii) reduction in workload due to the elimination or reduction of retyping of data by laboratory personnel, and (iv) decreased errors in data reported to sample submitters. Such a database provides a complete record of when and how often laboratory reference materials have been analyzed and provides a record of what correction factors have been used through time. It provides an audit trail for stable isotope laboratories. Since the original publication of the manual for LIMS for Light Stable Isotopes, the isotopes 3 H, 3 He, and 14 C, and the chlorofluorocarbons (CFCs), CFC-11, CFC-12, and CFC-113, have been added to this program.

Open-File Report↗

USGS1998

A manual for a Laboratory Information Management System (LIMS) for light stable isotopes

Open-File Report↗

USGS1999

SutraPlot, a graphical post-processor for SUTRA, a model for ground-water flow with solute or energy transport

This report documents a graphical display post-processor (SutraPlot) for the U.S. Geological Survey Saturated-Unsaturated flow and solute or energy TRAnsport simulation model SUTRA, Version 2D3D.1. This version of SutraPlot is an upgrade to SutraPlot for the 2D-only SUTRA model (Souza, 1987). It has been modified to add 3D functionality, a graphical user interface (GUI), and enhanced graphic output options. Graphical options for 2D SUTRA (2-dimension) simulations include: drawing the 2D finite-element mesh, mesh boundary, and velocity vectors; plots of contours for pressure, saturation, concentration, and temperature within the model region; 2D finite-element based gridding and interpolation; and 2D gridded data export files. Graphical options for 3D SUTRA (3-dimension) simulations include: drawing the 3D finite-element mesh; plots of contours for pressure, saturation, concentration, and temperature in 2D sections of the 3D model region; 3D finite-element based gridding and interpolation; drawing selected regions of velocity vectors (projected on principal coordinate planes); and 3D gridded data export files. Installation instructions and a description of all graphic options are presented. A sample SUTRA problem is described and three step-by-step SutraPlot applications are provided. In addition, the methodology and numerical algorithms for the 2D and 3D finite-element based gridding and interpolation, developed for SutraPlot, are described. 1

Open-File Report↗

USGS1982

Digital classification of Landsat data for vegetation and land-cover mapping in the Blackfoot River watershed, southeastern Idaho

This paper documents the procedures, results, and final products of a digital analysis of Landsat data used to produce a vegetation and landcover map of the Blackfoot River watershed in southeastern Idaho. Resource classes were identified at two levels of detail: generalized Level I classes (for example, forest land and wetland) and detailed Levels II and III classes (for example, conifer forest, aspen, wet meadow, and riparian hardwoods). Training set statistics were derived using a modified clustering approach. Environmental stratification that separated uplands from lowlands improved discrimination between resource classes having similar spectral signatures. Digital classification was performed using a maximum likelihood algorithm. Classification accuracy was determined on a single-pixel basis from a random sample of 25-pixel blocks. These blocks were transferred to small-scale color-infrared aerial photographs, and the image area corresponding to each pixel was interpreted. Classification accuracy, expressed as percent agreement of digital classification and photo-interpretation results, was 83.0:t 2.1 percent (0.95 probability level) for generalized (Level I) classes and 52.2:t 2.8 percent (0.95 probability level) for detailed (Levels II and III) classes. After the classified images were geometrically corrected, two types of maps were produced of Level I and Levels II and III resource classes: color-coded maps at a 1:250,000 scale, and flatbed-plotter overlays at a 1:24,000 scale. The overlays are more useful because of their larger scale, familiar format to users, and compatibility with other types of topographic and thematic maps of the same scale.

Professional Paper↗