Geology Reports⌕ Search

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 667 records · Page 37Linked to original sources
USGS1991

The effects of overburden stress on the specific storage and hydraulic conductivity of artesian aquifers

A mathematical algorithm is developed to determine the depth-dependent profiles of specific storage and hydraulic conductivity resulting from overburden stress in horizontally isotropic artesian aquifers. Vertical variations in the void ratio of the aquifer matrix brought about by overburden stresses determine the pore-volume compressibility and matrix permeability at specific depths within the aquifers which, in turn, determine the depth-dependent profiles of specific storage and hydraulic conductivity. Time-drawdown curves are obtained for two sands subjected to different overburden stresses. For shallow artesian aquifers with low overburden stress and high aquifer matrix compressibility, noticeable vertical gradients in specific storage occur. These vertical gradients cause deviations from the classical time-drawdown curves defined by the Theis solution. These deviations are negligible for deep artesian aquifers.

Journal of Hydrology↗
USGS1982

An oilspill trajectory analysis model with a variable wind deflection angle

The oilspill trajectory movement algorithm consists of a vector sum of the surface drift component due to wind and the surface current component. In the U.S. Geological Survey oilspill trajectory analysis model, the surface drift component is assumed to be 3.5% of the wind speed and is rotated 20 degrees clockwise to account for Coriolis effects in the Northern Hemisphere. Field and laboratory data suggest, however, that the deflection angle of the surface drift current can be highly variable. An empirical formula, based on field observations and theoretical arguments relating wind speed to deflection angle, was used to calculate a new deflection angle at each time step in the model. Comparisons of oilspill contact probabilities to coastal areas calculated for constant and variable deflection angles showed that the model is insensitive to this changing angle at low wind speeds. At high wind speeds, some statistically significant differences in contact probabilities did appear.

Ocean Engineering↗
USGS1987

Moment tensor solutions estimated using optimal filter theory for 51 selected earthquakes, 1980-1984

The 51 global events that occurred from January 1980 to March 1984, which were chosen by the convenors of the Symposium on Seismological Theory and Practice, have been analyzed using a moment tensor inversion algorithm (Sipkin). Many of the events were routinely analyzed as part of the National Earthquake Information Center's (NEIC) efforts to publish moment tensor and first-motion fault-plane solutions for all moderate- to large-sized (mb>5.7) earthquakes. In routine use only long-period P-waves are used and the source-time function is constrained to be a step-function at the source (??-function in the far-field). Four of the events were of special interest, and long-period P, SH-wave solutions were obtained. For three of these events, an unconstrained inversion was performed. The resulting time-dependent solutions indicated that, for many cases, departures of the solutions from pure double-couples are caused by source complexity that has not been adequately modeled. These solutions also indicate that source complexity of moderate-sized events can be determined from long-period data. Finally, for one of the events of special interest, an inversion of the broadband P-waveforms was also performed, demonstrating the potential for using broadband waveform data in inversion procedures. ?? 1987.

Physics of the Earth and Planetary Interiors↗
USGS1988

Use of airborne imaging spectrometer data to map minerals associated with hydrothermally altered rocks in the northern Grapevine Mountains, Nevada, and California

Three flightlines of Airborne Imaging Spectrometer (AIS) data, acquired over the northern Grapevine Mountains, Nevada, and California, were used to map minerals associated with hydrothermally altered rocks. The data were processed to remove vertical striping, normalized using an equal area normalization, and reduced to reflectance relative to an average spectrum derived from the data. An algorithm was developed to automatically calculate the absorption band parameters band position, band depth, and band width for the strongest absorption feature in each pixel. These parameters were mapped into an intensity, hue, saturation (IHS) color system to produce a single color image that summarized the absorption band information, This image was used to map areas of potential alteration based upon the predicted relationships between the color image and mineral absorption band. Individual AIS spectra for these areas were then examined to identify specific minerals. Two types of alteration were mapped with the AIS data. Areas of quartz-sericite-pyrite alteration were identified based upon a strong absorption feature near 2.21 μm, a weak shoulder near 2.25 μm, and a weak absorption band near 2.35 μm caused by sericite (fine-grained muscovite). Areas of argillic alteration were defined based on the presence of montmorillonite, identified by a weak to moderate absorption feature near 2.21 μm and the absence of the 2.35 μm band. Montmorillonite could not be identified in mineral mixtures. Calcite and dolomite were identified based on sharp absorption features near 2.34 and 2.32 μm, respectively. Areas of alteration identified using the AIS data corresponded well with areas mapped using field mapping, field reflectance spectra, and laboratory spectral measurements.

California, Nevada↗
USGS1993

Mapping playa evaporite minerals with AVIRIS data: A first report from Death Valley, California

Efflorescent salt crusts in Death Valley, California, were mapped by using Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) data and a recently developed least-squares spectral band-fitting algorithm. Eight different saline minerals were remotely identified, including three borates, hydroboracite, pinnoite, and rivadavite, that have not been previously reported from the Death Valley efflorescent crusts. The three borates are locally important phases in the crusts, and at least one of the minerals, rivadavite, appears to be forming directly from brine. Borates and other evaporite minerals provide a basis for making remote chemical measurements of desert hydrologic systems. For example, in the Eagle Borax Spring area, the AVIRIS mineral maps pointed to elevated magnesium and boron levels in the ground waters, and to the action of chemical divides causing subsurface fractionation of calcium. Many other chemical aspects of playa brines should have an expression in the associated evaporite assemblages. Certain anhydrous evaporites, including anhydrite, glauberite, and thenardite, lack absorption bands in the visible and near-infrared wavelength range, and crusts composed of these minerals could not be characterized by using AVIRIS. In these situations, thermal-infrared remote sensing data may complement visible and near-infrared data for mapping evaporites. Another problem occurred in wet areas of Death Valley, where water absorption caused low signal levels in the 2.0-2.5 ??m wavelength region that obscured any spectral features of evaporite minerals. Despite these difficulties, the results of this study demonstrate the potential for using AVIRIS and other imaging spectrometer data to study playa chemistry. Such data can be useful for understanding chemical linkages between evaporites and ground waters, and will facilitate studies of how desert ground-water regimes change through time in response to climatic and other variables.

California↗
USGS1995

Use of semipermeable membrane devices (SPMDS) to determine bioavailable organochlorine pesticide residues in streams receiving irrigation drainwater

The semipermeable membrane device (SPMD), consisting of a neutral lipid (triolein) enclosed in polyethylene layflat tubing, is very effective in sequestering bioavailable organochlorine (OC) pesticides in the environment. We used SPMDs to sequester OC pesticide residues in streams receiving irrigation drainwater and found toxaphene, the DDT complex, dieldrin, and endrin. Ambient water concentrations of the OC pesticides were calculated using an algorithm developed previously. Toxaphene residues were estimated to range from 300 to 7,000 ng/L.

Chemosphere↗
USGS1977

A method and fortran program for the computerized reconstruction of three-dimensional objects from serial sections

Three-dimensional objects can be reconstructed automatically from sets of serial sections by using a suitable algorithm. For two-phase systems it is sufficient to match phase boundary loops between adjacent sections. When all possible neighboring pairs of sections have been compared in this manner, the bounding surfaces of the object are described by sets of triangular and nonplanar quadrilateral patches. ?? 1977.

Computers & Geosciences↗
USGS1992

A spatial model to aggregate point-source and nonpoint-source water-quality data for large areas

More objective and consistent methods are needed to assess water quality for large areas. A spatial model, one that capitalizes on the topologic relationships among spatial entities, to aggregate pollution sources from upstream drainage areas is described that can be implemented on land surfaces having heterogeneous water-pollution effects. An infrastructure of stream networks and drainage basins, derived from 1:250,000-scale digital-elevation models, define the hydrologic system in this spatial model. The spatial relationships between point- and nonpoint pollution sources and measurement locations are referenced to the hydrologic infrastructure with the aid of a geographic information system. A maximum-branching algorithm has been developed to simulate the effects of distance from a pollutant source to an arbitrary downstream location, a function traditionally employed in deterministic water quality models. ?? 1992.

Computers & Geosciences↗
USGS1981

Significant results from using earth observation satellites for mineral and energy resource exploration

A large number of Earth-observation satellites orbit our world several times each day, providing new information about the land and sea surfaces and the overlying thin layer of atmosphere that makes our planet unique. Meteorological satellites have had the longest history of experimental use and most are now considered operational. The geologic information collected by the Landsat, Polar Orbiting Geophysical Observatory (POGO), Magsat, Heat Capacity Mapping Mission (HCMM) and Seasat land and ocean observation systems is being thoroughly tested, and some of these systems are now approaching operational use. Landsat multispectral images provide views of large areas of the Earth under uniform lighting conditions and can be obtained at a variety of scales and formats. Not only do the Landsat data provide highly useful images showing surficial materials and structures such as folds and faults, but also measurements and computer-derived ratios of the brightness of different rock types, alteration zones, and mineral associations. These data have led to the finding of a variety of new ore deposits. In addition, the combination of Landsat digital data and aeromagnetic data has extended the use of Landsat as an exploration tool which can be used to readily relate surface features to subsurface anomalies. Magsat data, now being collected, are helping refine information on major crustal anomalies that were first recognized during the analysis of POGO data. The more nearly circular orbit, lower altitude, and increased sophistication of its vector magnetometer enable Magsat to provide more precise information than POGO. Information of this type is required to develop crustal models. Although Magsat is designed to operate for only 4–8 months, the number of orbits that it should be able to make will be sufficient to accomplish its mission and to record a major magnetic storm expected in 1980. HCMM is a two-band visible to near-IR (0.55–1.1 μm) and thermal infrared (10.2–12.5 μm) system designed to measure reflected solar energy, determine the heat capacity of rocks and to monitor soil moisture, thermal effluents, plant canopy temperatures and snow cover. Launched in April 1978, it is in sun-synchronous, circular orbit at an altitude of 620 km. It is a relatively low-resolution system with an instantaneous field of view (IFOV) of 500–600 m and a swath width of 716 km. However, the system is designed to detect objects in the range of 260°–340° K with a sensitivity (NEδT) of 0.4°K at 280°. Recording the thermal radiation of urban heat islands and high thermal inertia of quartzite strata in the Appalachian region are two examples of its land applications. Launched in June 1978, Seasat operated for only 100 days, but successfully acquired much information over both sea and land. The collection of synthetic aperture radar (SAR) imagery and radar altimetry was particularly important to geologists. Although there are difficulties in processing and distributing these data in a timely manner, initial evaluations indicate that the radar imagery supplements Landsat data by increasing the spectral range and offering a different look angle. The radar altimeter provides accurate profiles over narrow strips of land (1 km wide) and has demonstrated usefulness in measuring icecap surfaces (Greenland, Iceland, and Antarctica). The Salar of Uyuni in southern Bolivia served as a calibration site for the altimeter and has enabled investigators to develop a land-based smoothing algorithm that is believed to increase the accuracy of the system to 10 cm. Data from the altimeter are currently being used to measure subsidence resulting from ground water withdrawal in the Phoenix-Tucson area.

Advances in Space Research↗
USGS1979

A two-dimensional hydrodynamic model of a tidal estuary

A finite element model is described which is used in the computation of tidal currents in an estuary. This numerical model is patterned after an existing algorithm and has been carefully tested in rectangular and curve-sided channels with constant and variable depth. One of the common uncertainties in this class of two-dimensional hydrodynamic models is the treatment of the lateral boundary conditions. Special attention is paid specifically to addressing this problem. To maintain continuity within the domain of interest, ‘smooth’ curve-sided elements must be used at all shoreline boundaries. The present model uses triangular, isoparametric elements with quadratic basis functions for the two velocity components and a linear basis function for water surface elevation. An implicit time integration is used and the model is unconditionally stable. The resultant governing equations are nonlinear owing to the advective and the bottom friction terms and are solved iteratively at each time step by the Newton-Raphson method. Model test runs have been made in the southern portion of San Francisco Bay, California (South Bay) as well as in the Bay west of Carquinez Strait. Owing to the complex bathymetry, the hydrodynamic characteristics of the Bay system are dictated by the generally shallow basins which contain deep, relict river channels. Great care must be exercised to ensure that the conservation equations remain locally as well as globally accurate. Simulations have been made over several representative tidal cycles using this finite element model, and the results compare favourably with existing data. In particular, the standing wave in South Bay and the progressive wave in the northern reach are well represented.

California↗
USGS1989

A comparison of two finite element models of tidal hydrodynamics using a North Sea data set

Using the region of the English Channel and the southern bight of the North Sea, we systematically compare the results of two independent finite element models of tidal hydrodynamics. The model intercomparison provides a means for increasing our understanding of the relevant physical processes in the region in question as well as a means for the evaluation of certain algorithmic procedures of the two models.

Advances in Water Resources↗
USGS1992

Allozyme evidence for crane systematics and polymorphisms within populations of Sandhill, Sarus, Siberian and whooping cranes

Electrophoretic analysis of proteins yielded evidence on the relationships of species of cranes and on genetic diversity within populations of some species. Diversity within the Greater Sandhill crane and a Florida population of the Florida Sandhill crane was similar to that of most other vertebrates, but diversity was low in the Mississippi Sandhill crane, in the Okefenokee population of the Florida Sandhill crane, and within the Siberian and Sarus cranes. Diversity was surprisingly high among whooping cranes, whose number dropped to less than 25 early in this century. Phylogenetic analysis, using both character state and distance algorithms, yielded highly concordant trees for the 15 species. The African crowned cranes ( Balearica ) were widely divergent from all other cranes. Species of Anthropoides, Bugeranus , and Grus clustered closely but sorted into two lineages: a Whooper Group consisted of the whooping, common, hooded, black-necked, white-naped, and red-crowned cranes of genus Grus ; and a Sandhill Group included the Sandhill, Siberian, Sarus, and Brolga cranes of genus Grus , the wattled crane of genus Bugeranus , and the Demoiselle and blue cranes of genus Anthropoides .

Florida↗
USGS2016

Regional flow duration curves: Geostatistical techniques versus multivariate regression

A period-of-record flow duration curve (FDC) represents the relationship between the magnitude and frequency of daily streamflows. Prediction of FDCs is of great importance for locations characterized by sparse or missing streamflow observations. We present a detailed comparison of two methods which are capable of predicting an FDC at ungauged basins: (1) an adaptation of the geostatistical method, Top-kriging, employing a linear weighted average of dimensionless empirical FDCs, standardised with a reference streamflow value; and (2) regional multiple linear regression of streamflow quantiles, perhaps the most common method for the prediction of FDCs at ungauged sites. In particular, Top-kriging relies on a metric for expressing the similarity between catchments computed as the negative deviation of the FDC from a reference streamflow value, which we termed total negative deviation (TND). Comparisons of these two methods are made in 182 largely unregulated river catchments in the southeastern U.S. using a three-fold cross-validation algorithm. Our results reveal that the two methods perform similarly throughout flow-regimes, with average Nash-Sutcliffe Efficiencies 0.566 and 0.662, (0.883 and 0.829 on log-transformed quantiles) for the geostatistical and the linear regression models, respectively. The differences between the reproduction of FDC's occurred mostly for low flows with exceedance probability (i.e. duration) above 0.98.

Advances in Water Resources↗
USGS2023

Operationalizing crop model data assimilation for improved on-farm situational awareness

The ability of ‘digital agriculture’ to support on-farm decision making is predicated on the real-time combination of observations and prior knowledge into an integrated digital environment. The mathematical discipline that seeks to provide this integration is known as model data assimilation (DA), with demonstrated benefits including improved predictive reliability, and the capacity to identify unexpected changes in field conditions and potential measurement errors. Despite routine adoption in other fields, the delayed adoption of DA in agriculture is due to the need to express end-of-season outcomes such as yield, update forecasts of these outcomes throughout the growing season as data become available, and enhance forecast reliability. To overcome these challenges, three guiding principles are introduced, providing a means to operationalize crop model DA for robust on-farm decision support. We apply the guiding principles using a South Australian viticulture case study. Our case study involves application of an iterative form of a widely used DA algorithm (ensemble Kalman filter) to dynamically update both static parameters and states associated with a grapevine simulation model. Daily weather data as well as fortnightly ground-based leaf area index (LAI) data are used for assimilation. It is shown how crop model DA can lead to not only significant improvements in forecasts of LAI but also to forecasts of end-of-season yield. The guiding principles also enable observations of greatest value to be identified throughout the season. This study highlights the role that formal crop model DA can play in agricultural decision support through enhancing situational awareness in real time.

Agricultural and Forest Meteorology↗
USGS2014

Characterizing recent and projecting future potential patterns of mountain pine beetle outbreaks in the Southern Rocky Mountains

The recent widespread mountain pine beetle (MPB) outbreak in the Southern Rocky Mountains presents an opportunity to investigate the relative influence of anthropogenic, biologic, and physical drivers that have shaped the spatiotemporal patterns of the outbreak. The aim of this study was to quantify the landscape-level drivers that explained the dynamic patterns of MPB mortality, and simulate areas with future potential MPB mortality under projected climate-change scenarios in Grand County, Colorado, USA. The outbreak patterns of MPB were characterized by analysis of a decade-long Landsat time-series stack, aided by automatic attribution of change detected by the Landsat-based Detection of Trends in Disturbance and Recovery algorithm (LandTrendr). The annual area of new MPB mortality was then related to a suite of anthropogenic, biologic, and physical predictor variables under a general linear model (GLM) framework. Data from years 2001–2005 were used to train the model and data from years 2006–2011 were retained for validation. After stepwise removal of non-significant predictors, the remaining predictors in the GLM indicated that neighborhood mortality, winter mean temperature anomaly, and residential housing density were positively associated with MPB mortality, whereas summer precipitation was negatively related. The final model had an average area under the curve (AUC) of a receiver operating characteristic plot value of 0.72 in predicting the annual area of new mortality for the independent validation years, and the mean deviation from the base maps in the MPB mortality areal estimates was around 5%. The extent of MPB mortality will likely expand under two climate-change scenarios (RCP 4.5 and 8.5) in Grand County, which implies that the impacts of MPB outbreaks on vegetation composition and structure, and ecosystem functioning are likely to increase in the future.

Colorado↗
USGS2008

Satellite remotely-sensed land surface parameters and their climatic effects for three metropolitan regions

By using both high-resolution orthoimagery and medium-resolution Landsat satellite imagery with other geospatial information, several land surface parameters including impervious surfaces and land surface temperatures for three geographically distinct urban areas in the United States – Seattle, Washington, Tampa Bay, Florida, and Las Vegas, Nevada, are obtained. Percent impervious surface is used to quantitatively define the spatial extent and development density of urban land use. Land surface temperatures were retrieved by using a single band algorithm that processes both thermal infrared satellite data and total atmospheric water vapor content. Land surface temperatures were analyzed for different land use and land cover categories in the three regions. The heterogeneity of urban land surface and associated spatial extents were shown to influence surface thermal conditions because of the removal of vegetative cover, the introduction of non-transpiring surfaces, and the reduction in evaporation over urban impervious surfaces. Fifty years of in situ climate data were integrated to assess regional climatic conditions. The spatial structure of surface heating influenced by landscape characteristics has a profound influence on regional climate conditions, especially through urban heat island effects.

Advances in Space Research↗
USGS2010

Optimal control of native predators

We apply decision theory in a structured decision-making framework to evaluate how control of raccoons ( Procyon lotor ), a native predator, can promote the conservation of a declining population of American Oystercatchers ( Haematopus palliatus ) on the Outer Banks of North Carolina. Our management objective was to maintain Oystercatcher productivity above a level deemed necessary for population recovery while minimizing raccoon removal. We evaluated several scenarios including no raccoon removal, and applied an adaptive optimization algorithm to account for parameter uncertainty. We show how adaptive optimization can be used to account for uncertainties about how raccoon control may affect Oystercatcher productivity. Adaptive management can reduce this type of uncertainty and is particularly well suited for addressing controversial management issues such as native predator control. The case study also offers several insights that may be relevant to the optimal control of other native predators. First, we found that stage-specific removal policies (e.g., yearling versus adult raccoon removals) were most efficient if the reproductive values among stage classes were very different. Second, we found that the optimal control of raccoons would result in higher Oystercatcher productivity than the minimum levels recommended for this species. Third, we found that removing more raccoons initially minimized the total number of removals necessary to meet long term management objectives. Finally, if for logistical reasons managers cannot sustain a removal program by removing a minimum number of raccoons annually, managers may run the risk of creating an ecological trap for Oystercatchers.

North Carolina↗
USGS2021

Prioritizing restoration areas to conserve multiple sagebrush-associated wildlife species

Strategic restoration of altered habitat is one method for addressing worldwide biodiversity declines. Within the sagebrush steppe of western North America, habitat degradation has been linked to declines in many species, making restoration a priority for managers; however, limited funding, spatiotemporal variation in restoration success, and the need to manage for diverse wildlife species make decision-making regarding restoration actions challenging. To address the challenge of spatial conservation prioritization, we developed the Prioritizing Restoration of Sagebrush Ecosystems Tool (PReSET). This decision support tool utilizes the prioritizr package in program R and an integer linear programming algorithm to select parcels representing both high biodiversity value and high probability of restoration success. We tested PReSET on a sagebrush steppe system within southwestern Wyoming using distributional data for six species with diverse life histories and a spatial layer of predicted sagebrush recovery times to identify restoration targets at both broad and local scales. While the broad-scale portion of our tool outputs can inform policy, the local-scale results can be applied directly to on-the-ground restoration. We identified restoration priority areas with greater precision than existing spatial prioritizations and incorporated range differences among species. We noted tradeoffs, including that restoring for habitat connectivity may require restoration actions in areas with lower probability of success. Future applications of PReSET will draw from emerging datasets, including spatially-varying economic costs of restoration, animal movement data, and additional species, to further improve our ability to target effective sagebrush restoration.

Wyoming↗