USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geosciences”

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 361 records · Page 20Linked to original sources
USGS1988

Normalization of oxygen and hydrogen isotope data

To resolve confusion due to expression of isotopic data from different laboratories on non-corresponding scales, oxygen isotope analyses of all substances can be expressed relative to VSMOW or VPDB (Vienna Peedee belemnite) on scales normalized such that the δ 18 O of SLAP is −55.5% relative to VSMOW. H 3 + contribution in hydrogen isotope ratio analysis can be easily determined using two gaseous reference samples that differ greatly in deuterium content.

Chemical Geology: Isotope Geoscience Section↗
USGS1989

Paleozoic age of the Capo Spartivento Orthogneiss, Sardinia, Italy

Zircon U Pb "> UPb isotope dating of the Capo Spartivento Orthogneiss, proposed as a possible Precambrian basement of southern Sardinia, shows that this rock is Caledonian in age. Conventional multi-grain analyses yield an imprecise age of roughly 480 Ma, and ion-microprobe analyses of cores of single grains yield a consistent age of 449 Ma. Though some inherited grains of Proterozoic age are present in the zircon population, they are neither abundant nor consistent with Caledonian growth of new zircons within an older protolith.

Chemical Geology: Isotope Geoscience Section↗
USGS1989

Contrasting zircon morphology and UPb systematics in peralkaline and metaluminous post-orogenic granite complexes of the Arabian Shield, Kingdom of Saudi Arabia

Uzircon ages are reported for seven metaluminous-to-peralkaline post-orogenic granites from the Late Proterozoic Arabian Shield of Saudi Arabia. Zircons from the metaluminous rocks are prismatic, with length-to-width ratios of ∼ 2–4: 1 and small pyramidal terminations. In contrast, zircons from three of the four peralkaline complexes either lack well-developed prismatic faces (are pseudo-octahedral) or are anhedral. Some zircons from the peralkaline granites contain inherited radiogenic Pb and have very high common Pb contents ( 206 Pb/ 204 Pb < 150), making the UPb method poorly suited for determining the age of these rocks. Zircons in the metaluminous granites do not contain inheritance and yield well-defined concordia intercepts. The span of ages of the seven complexes (670-470 Ma) indicates that post-orogenic granitic magmatism was not a singular event in the Arabian Shield but rather occurred as multiple intrusive episodes from the Late Proterozoic to the Middle Ordovician.

Chemical Geology: Isotope Geoscience Section↗
USGS1990

UPb ages of zircon rims: A new analytical method using the air-abrasion technique

We present a new technique for directly dating, by conventional techniques, the rims of zircons. Several circumstances, such as a xenocrystic or inherited component in igneous zircon and metamorphic overgrowths on igneous cores, can result in grains with physically distinct age components. Pneumatic abrasion has been previously shown by Krogh to remove overgrowths and damaged areas of zircon, leaving more resistant and isotopically less disturbed parts available for analysis. A new abrader design, which is capable of very gently grinding only tips and interfacial edges of even needle-like grains, permits easy collection of abraded material for dating. Five examples demonstrate the utility of the “dust-collecting” technique, including two studies that compare conventional, ion microprobe and abrader data. Common Pb may be strongly concentrated in the outermost zones of many zircons and this Pb is not easily removed by leaching (even in weak HF). Thus, the benefit of removing only the outermost zones (and avoiding mixing of age components) is somewhat compromised by the much higher common Pb contents which result in less precise age determinations. A very brief abrasion to remove the high common Pb zones prior to collection of material for dating is selected.

Chemical Geology: Isotope Geoscience Section↗
USGS1991

Resetting of RbSr ages of volcanic rocks by low-grade burial metamorphism

We report a nine-point RbSr whole-rock isochron age of 70±3 Ma (MSWD 3.97) for Mid-Jurassic volcanic rocks. The same rocks have also been dated by the UThPb method on zircon, giving a crystallization age of 166 ± 11 Ma, over twice as old as the RbSr age. The data demonstrate that whole-rock RbSr ages of volcanic rocks, even lava flows with SiO 2 content as low as 57 wt.%, are susceptible to complete resetting. The rocks range in composition from rhyodacite tuffs to andesite lavas. The complete breakdown of all major minerals that contain Rb and Sr resulted in an alteration mineral assemblage consisting of phengite, albite, secondary quartz, and minor amounts of chlorite and epidote. Phengite is the K-bearing product of the breakdown of biotite and K-feldspar. Pressure during low-grade metamorphism of the volcanic rocks, estimated from phengite composition to have been in the range of 4 to 6 kbar, points to thrust-related burial as the main cause of resetting. Consequently, such reset isochrons may date large-scale events such as regional thrusting and metamorphism. The coherent resetting of the RbSr isochron suggests large-scale pervasive fluid movement during thrust-related burial metamorphism.

Chemical Geology: Isotope Geoscience Section↗
USGS1991

Caution on the use of Viton® or FETFE® O-rings in carbon dioxide sample containers for δ 18 0 analysis

After 10 days, 3-&mu;mol CO 2 samples in containers having glass stopcocks with Viton &reg; or FETFE O-rings were enriched in 18 O by 1.5% as a result of absorption by the elastomer; this amount of enrichment is &sim;20 times greater than the precision of &delta; 18 O measurements of most laboratories. No change in 13 C content was observed. Increasing the sample size to 100 &mu;mol resulted in an 18 O enrichment of 0.2% and did not affect the 13 C content. Caution needs to be exercised in selecting sample containers for CO 2 isotope-ratio samples of < 200 &mu;mol. If stopcocks are used in construction of containers for such samples, the use of all-glass stopcocks with Apiezon N &reg; hydrocarbon-based grease will eliminate the fractionation of oxygen isotopes.

Chemical Geology: Isotope Geoscience Section↗
USGS2023

Rapid estimation of minimum depth-to-bedrock from lidar leveraging deep-learning-derived surficial material maps

Previously glaciated landscapes often share similar surficial characteristics, including large areas of exposed bedrock, blankets of till deposits, and alluvium-floored valleys. These materials play significant roles in geologic and hydrologic resources, geohazards, and landscape evolution; however, the vast extents of many previously glaciated landscapes have rendered comprehensive, detailed field mapping difficult. While recent advances in remote sensing have facilitated mapping of surficial materials and landforms, manual map creation has remained a time-intensive task. The development of convolutional neural networks (CNNs) for image classification has provided a new opportunity for rapid characterization of digital elevation models, thus enabling efficient mapping of surficial materials and landforms. We have developed a methodology that leverages existing geologic maps and high-resolution (1–3 m) lidar data to train a U-Net CNN to classify alluvium and exposed bedrock in previously glaciated regions. Coupled with U.S. Geological Survey-developed geomorphometry tools capable of approximating stream incision depths, these classifications can be used to estimate the minimum thicknesses of stream-proximal hillslope sediments in areas where streams have undergone minimal incision into bedrock. We validate this approach in the context of the Neversink River watershed, a subbasin of the Delaware River Basin and significant water source for New York City. Evaluation of deep learning model performance demonstrates substantial agreement with manually drawn maps of alluvium and exposed bedrock. Validation of the minimum sediment thickness map using borehole data and passive seismic measurements shows the greatest performance for shallow materials and decreased performance in deep sediments, as well as in areas where bedrock exposures were too small to be resolved by lidar. To resolve these issues and create more accurate surficial maps, we are training new CNNs with additional geologic data and exploring advanced approaches for estimating depths of stream incision.

Delaware, New Jersey, New York, Pennsylvania↗
USGS2023

Near surface sediments introduce low frequency noise into gravity models

3D geologic modeling and mapping often relies on gravity modeling to identify key geologic structures, such as basin depth, fault offset, or fault dip. Such gravity models generally assume either homogeneous or spatially uncorrelated densities within modeled rock bodies and overlying sediments, with average densities typically derived from surface and drill-hole sampling. The noise contributed to the gravity anomaly by these density assumptions is zero in the homogeneous case and typically <200 μGal in the uncorrelated case. Rock bodies and sediments, however, show both a range of densities and spatial correlation of these densities, in both surface and drill-hole samples, and this correlation causes an increase in power in the low frequency content of the resulting gravity anomaly. Spatial correlation of densities can be modeled as a Gaussian random field (GRF), with the random field parameters derived from drill-hole and geologic map data. Data from alluvial fan sediments in southern Nevada indicate correlation lengths of up to 300 m in the vertical dimension and kilometers in the horizontal dimension. The resulting GRF density models show that the noise contributed to the measured gravity anomaly is of low frequency and can be several mGal in amplitude, contradicting the common attribution of lower frequencies to deeper sources. This low-frequency noise increases in power with an increase in sediment thickness. Its presence increases the ambiguity of interpretations of subsurface geologic body shape, such as basin analyses that attempt to quantify concealed basement fault depths, offsets, and dip angles. In the southwestern United States, where basin analyses are important for natural resource applications, such ambiguity increases the uncertainty of subsequent process modeling.

Applied Computing and Geosciences↗
USGS2025

Predictive regressive models of recent marsh sediment thickness improve the quantification of coastal marsh sediment budgets

Coastal marsh wetlands experience variations in vertical gains and losses through time, which have allowed them to infill relict topography and record variations in drivers. The stratigraphic unit associated with the development of the marsh also reflects the long-term importance of key ecosystem services supplied by the marsh environment, including carbon storage and storm mitigation. Mapping these coastal wetland sediments and the marsh unit thickness is challenging as traditional coastal geophysical tools are not easily deployable (acoustic methods) or are unreliable in saline-soil environments (e.g., ground-penetrating radar), leaving core-based methods the most viable mapping method. In the present study, we utilized prior information on the geologic architecture of the region to select spatial and physical metrics that likely persisted throughout evolution of the marsh during the late Holocene. We then assessed the individual and collective power of these metrics to predict marsh thickness observed from cores. Employing regressive predictive models powered by these data, we improve the quantification of marsh thickness for a coastal fringing marsh within the Grand Bay estuary in Mississippi and Alabama (USA). The information gained from this approach yields improved estimates of the carbon stocks in this environment. Additionally, the stored sediment masses reflect the past, and potential future, persistence of the Grand Bay marsh under historical and present marsh-estuarine sediment exchange fluxes. Such improvements to both the sediment budget of recent marsh stratigraphic units and the spatial extent provide new resources for comparison with large-scale landscape models, the latter of which may be used, when validated, to predict future change and ecosystem transformations.

Alabama. Mississippi↗
USGS2025

Statistical approaches for modeling correlated grade and tonnage distributions and applications for mineral resource assessments

Correlations between grade and tonnage exist in mineral resource data compiled from published reports, but they are not always addressed during quantitative assessment of undiscovered mineral resources. Failure to account for correlated grade and tonnage distributions can result in geologically unrealistic assessment results. Current software tools simulate univariate ore tonnage and multivariate resource grades of undiscovered deposits independently. As a result, analysts are forced to rely on ad-hoc solutions to minimize the correlation issues by: 1) creating subsets of data with restricted criteria; 2) truncating grade and tonnage distributions; and 3) testing model robustness using exploratory data analysis. While these methods represent pragmatic solutions, the statistical solutions presented here provide additional options to address real correlations in grade and tonnage data used for mineral resource assessments. We present a modified version of the MapMark4 package in R that introduces two alternatives for modeling grade and tonnage distributions, consisting of a multivariate solution that accounts for correlations between ore tonnage and metal grades and an empirical solution that utilizes simple random sampling with replacement to reproduce coupled grades and tonnages from the input data. We present simulations for contained ore and metal for three case studies representing tungsten skarn, komatiite-hosted nickel, and sediment-hosted carbonate amagmatic zinc-lead (Mississippi Valley-type) deposits. Employing the methods presented here yields quantitative mineral resource assessment results that more closely reflect the empirical distributions of grades and tonnages observed in nature and expands the applicability of these tools for ongoing critical mineral resource assessments.

Applied Computing and Geosciences↗
USGS2025

Extracting data from maps: Lessons learned from the artificial intelligence for critical mineral assessment competition

The U.S. Geological Survey (USGS), Defense Advanced Projects Research Agency (DARPA), NASA Jet Propulsion Laboratory (JPL), and MITRE ran a 12-week machine learning competition aimed at accelerating development of AI tools for critical mineral assessments. The Artificial Intelligence for Critical Mineral Assessment Competition solicited innovative solutions for two challenges: 1) automated georeferencing of historical maps, and 2) automated feature extraction from historical maps. Competitors used a new dataset of historical map images to train, validate, and evaluate their models. Automated georeferencing pipelines attained a median root-mean square error of 1.1 km. Prompt-based extraction (i.e., with user input) of polygons, polylines, and points from geologic maps yielded median F1-scores of 0.77, 0.56, 0.35, respectively. Geologic maps pose numerous challenges for AI workflows because they vary significantly. However, despite its short duration, the competition yielded promising results that have since spurred further innovation in this area and led to the development of new AI tools to semi-automate key, time-consuming parts of the assessment workflow.

Applied Computing and Geosciences↗
USGS2004

CORRELATOR 5.2 - A program for interactive lithostratigraphic correlation of wireline logs

The limited radius of investigation of petrophysical measurements made in boreholes and the relatively large distances between wells result in an incomplete sensing of the subsurface through well logging. CORRELATOR is a program for estimating geological properties between logged boreholes. An initial and fundamental step is the lithostratigraphic correlation of logs in different wells. The method employed by the program closely emulates the process of visual inspection used by experienced subsurface geologists in manual correlation. Mathematically, the determination of lithostratigraphical equivalence is based on the simultaneous assessment of similarity in shale content, similarity in the patterns of vertical variation in a petrophysical property that is measured with high vertical resolution, and spatial consistency of stratigraphic relationships as determined by an expert system. Multiple additional options for processing log readings allow maximization in the extraction of information from pairs of logs per well and great flexibility in the final display of results in the form of cross sections and dip diagrams. ?? 2004 Elsevier Ltd. All rights reserved.

Computers & Geosciences↗
USGS2004

A visual basic program to generate sediment grain-size statistics and to extrapolate particle distributions

Measures that describe and summarize sediment grain-size distributions are important to geologists because of the large amount of information contained in textural data sets. Statistical methods are usually employed to simplify the necessary comparisons among samples and quantify the observed differences. The two statistical methods most commonly used by sedimentologists to describe particle distributions are mathematical moments (Krumbein and Pettijohn, 1938) and inclusive graphics (Folk, 1974). The choice of which of these statistical measures to use is typically governed by the amount of data available (Royse, 1970). If the entire distribution is known, the method of moments may be used; if the next to last accumulated percent is greater than 95, inclusive graphics statistics can be generated. Unfortunately, earlier programs designed to describe sediment grain-size distributions statistically do not run in a Windows environment, do not allow extrapolation of the distribution's tails, or do not generate both moment and graphic statistics (Kane and Hubert, 1963; Collias et al., 1963; Schlee and Webster, 1967; Poppe et al., 2000) 1 . Owing to analytical limitations, electro-resistance multichannel particle-size analyzers, such as Coulter Counters, commonly truncate the tails of the fine-fraction part of grain-size distributions. These devices do not detect fine clay in the 0.6–0.1 μm range (part of the 11-phi and all of the 12-phi and 13-phi fractions). Although size analyses performed down to 0.6 μm microns are adequate for most freshwater and near shore marine sediments, samples from many deeper water marine environments (e.g. rise and abyssal plain) may contain significant material in the fine clay fraction, and these analyses benefit from extrapolation. The program (GSSTAT) described herein generates statistics to characterize sediment grain-size distributions and can extrapolate the fine-grained end of the particle distribution. It is written in Microsoft Visual Basic 6.0 and provides a window to facilitate program execution. The input for the sediment fractions is weight percentages in whole-phi notation (Krumbein, 1934; Inman, 1952), and the program permits the user to select output in either method of moments or inclusive graphics statistics (Fig. 1). Users select options primarily with mouse-click events, or through interactive dialogue boxes.

Computers & Geosciences↗
USGS2005

Use of XML and Java for collaborative petroleum reservoir modeling on the Internet

The GEMINI (Geo-Engineering Modeling through INternet Informatics) is a public-domain, web-based freeware that is made up of an integrated suite of 14 Java-based software tools to accomplish on-line, real-time geologic and engineering reservoir modeling. GEMINI facilitates distant collaborations for small company and academic clients, negotiating analyses of both single and multiple wells. The system operates on a single server and an enterprise database. External data sets must be uploaded into this database. Feedback from GEMINI users provided the impetus to develop Stand Alone Web Start Applications of GEMINI modules that reside in and operate from the user's PC. In this version, the GEMINI modules run as applets, which may reside in local user PCs, on the server, or Java Web Start. In this enhanced version, XML-based data handling procedures are used to access data from remote and local databases and save results for later access and analyses. The XML data handling process also integrates different stand-alone GEMINI modules enabling the user(s) to access multiple databases. It provides flexibility to the user to customize analytical approach, database location, and level of collaboration. An example integrated field-study using GEMINI modules and Stand Alone Web Start Applications is provided to demonstrate the versatile applicability of this freeware for cost-effective reservoir modeling. ?? 2005 Elsevier Ltd. All rights reserved.

Computers & Geosciences↗
USGS2005

Dynamic modeling of Tampa Bay urban development using parallel computing

Urban land use and land cover has changed significantly in the environs of Tampa Bay, Florida, over the past 50 years. Extensive urbanization has created substantial change to the region's landscape and ecosystems. This paper uses a dynamic urban-growth model, SLEUTH, which applies six geospatial data themes (slope, land use, exclusion, urban extent, transportation, hillside), to study the process of urbanization and associated land use and land cover change in the Tampa Bay area. To reduce processing time and complete the modeling process within an acceptable period, the model is recoded and ported to a Beowulf cluster. The parallel-processing computer system accomplishes the massive amount of computation the modeling simulation requires. SLEUTH calibration process for the Tampa Bay urban growth simulation spends only 10 h CPU time. The model predicts future land use/cover change trends for Tampa Bay from 1992 to 2025. Urban extent is predicted to double in the Tampa Bay watershed between 1992 and 2025. Results show an upward trend of urbanization at the expense of a decline of 58% and 80% in agriculture and forested lands, respectively.

Computers & Geosciences↗
USGS2005

Statistical analysis of water-quality data containing multiple detection limits: S-language software for regression on order statistics

Trace contaminants in water, including metals and organics, often are measured at sufficiently low concentrations to be reported only as values below the instrument detection limit. Interpretation of these "less thans" is complicated when multiple detection limits occur. Statistical methods for multiply censored, or multiple-detection limit, datasets have been developed for medical and industrial statistics, and can be employed to estimate summary statistics or model the distributions of trace-level environmental data. We describe S-language-based software tools that perform robust linear regression on order statistics (ROS). The ROS method has been evaluated as one of the most reliable procedures for developing summary statistics of multiply censored data. It is applicable to any dataset that has 0 to 80% of its values censored. These tools are a part of a software library, or add-on package, for the R environment for statistical computing. This library can be used to generate ROS models and associated summary statistics, plot modeled distributions, and predict exceedance probabilities of water-quality standards. ?? 2005 Elsevier Ltd. All rights reserved.

Computers & Geosciences↗
USGS2006

Nonlinear inversion of potential-field data using a hybrid-encoding genetic algorithm

Using a genetic algorithm to solve an inverse problem of complex nonlinear geophysical equations is advantageous because it does not require computer gradients of models or "good" initial models. The multi-point search of a genetic algorithm makes it easier to find the globally optimal solution while avoiding falling into a local extremum. As is the case in other optimization approaches, the search efficiency for a genetic algorithm is vital in finding desired solutions successfully in a multi-dimensional model space. A binary-encoding genetic algorithm is hardly ever used to resolve an optimization problem such as a simple geophysical inversion with only three unknowns. The encoding mechanism, genetic operators, and population size of the genetic algorithm greatly affect search processes in the evolution. It is clear that improved operators and proper population size promote the convergence. Nevertheless, not all genetic operations perform perfectly while searching under either a uniform binary or a decimal encoding system. With the binary encoding mechanism, the crossover scheme may produce more new individuals than with the decimal encoding. On the other hand, the mutation scheme in a decimal encoding system will create new genes larger in scope than those in the binary encoding. This paper discusses approaches of exploiting the search potential of genetic operations in the two encoding systems and presents an approach with a hybrid-encoding mechanism, multi-point crossover, and dynamic population size for geophysical inversion. We present a method that is based on the routine in which the mutation operation is conducted in the decimal code and multi-point crossover operation in the binary code. The mix-encoding algorithm is called the hybrid-encoding genetic algorithm (HEGA). HEGA provides better genes with a higher probability by a mutation operator and improves genetic algorithms in resolving complicated geophysical inverse problems. Another significant result is that final solution is determined by the average model derived from multiple trials instead of one computation due to the randomness in a genetic algorithm procedure. These advantages were demonstrated by synthetic and real-world examples of inversion of potential-field data. ?? 2005 Elsevier Ltd. All rights reserved.

Computers & Geosciences↗
USGS2006

An efficient Matlab script to calculate heterogeneous anisotropically elastic wave propagation in three dimensions

We have created a second-order finite-difference solution to the anisotropic elastic wave equation in three dimensions and implemented the solution as an efficient Matlab script. This program allows the user to generate synthetic seismograms for three-dimensional anisotropic earth structure. The code was written for teleseismic wave propagation in the 1-0.1 Hz frequency range but is of general utility and can be used at all scales of space and time. This program was created to help distinguish among various types of lithospheric structure given the uneven distribution of sources and receivers commonly utilized in passive source seismology. Several successful implementations have resulted in a better appreciation for subduction zone structure, the fate of a transform fault with depth, lithospheric delamination, and the effects of wavefield focusing and defocusing on attenuation. Companion scripts are provided which help the user prepare input to the finite-difference solution. Boundary conditions including specification of the initial wavefield, absorption and two types of reflection are available. ?? 2005 Elsevier Ltd. All rights reserved.

Computers & Geosciences↗