Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Analytical Methods”

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 1,297 records · Page 72Linked to original sources
USGS2014

Linking aquifer spatial properties and non-Fickian transport in mobile-immobile like alluvial settings

Time-nonlocal transport models can describe non-Fickian diffusion observed in geological media, but the physical meaning of parameters can be ambiguous, and most applications are limited to curve-fitting. This study explores methods for predicting the parameters of a temporally tempered Lévy motion (TTLM) model for transient sub-diffusion in mobile–immobile like alluvial settings represented by high-resolution hydrofacies models. The TTLM model is a concise multi-rate mass transfer (MRMT) model that describes a linear mass transfer process where the transfer kinetics and late-time transport behavior are controlled by properties of the host medium, especially the immobile domain. The intrinsic connection between the MRMT and TTLM models helps to estimate the main time-nonlocal parameters in the TTLM model (which are the time scale index, the capacity coefficient, and the truncation parameter) either semi-analytically or empirically from the measurable aquifer properties. Further applications show that the TTLM model captures the observed solute snapshots, the breakthrough curves, and the spatial moments of plumes up to the fourth order. Most importantly, the a priori estimation of the time-nonlocal parameters outside of any breakthrough fitting procedure provides a reliable “blind” prediction of the late-time dynamics of subdiffusion observed in a spectrum of alluvial settings. Predictability of the time-nonlocal parameters may be due to the fact that the late-time subdiffusion is not affected by the exact location of each immobile zone, but rather is controlled by the time spent in immobile blocks surrounding the pathway of solute particles. Results also show that the effective dispersion coefficient has to be fitted due to the scale effect of transport, and the mean velocity can differ from local measurements or volume averages. The link between medium heterogeneity and time-nonlocal parameters will help to improve model predictability for non-Fickian transport in alluvial settings.

Journal of Hydrology↗
USGS2013

Analysis of changes in water-level dynamics at selected sites in the Florida Everglades

The historical modification and regulation of the hydrologic patterns in the Florida Everglades have resulted in changes in the ecosystem of South Florida and the Florida Everglades. Since the 1970s, substantial focus has been given to the restoration of the Everglades ecosystem. The U.S. Geological Survey through its Greater Everglades Priority Ecosystem Science and National Water-Quality Assessment Programs has been providing scientific information to resource managers to assist in the Everglades restoration efforts. The current investigation included development of a simple method to identify and quantify changes in historical hydrologic behavior within the Everglades that could be used by researchers to identify responses of ecological communities to those changes. Such information then could be used by resource managers to develop appropriate water-management practices within the Everglades to promote restoration. The identification of changes in historical hydrologic behavior within the Everglades was accomplished by analyzing historical time-series water-level data from selected gages in the Everglades using (1) break-point analysis of cumulative Z-scores to identify hydrologic changes and (2) cumulative water-level frequency distribution curves to evaluate the magnitude of those changes. This analytical technique was applied to six long-term water-level gages in the Florida Everglades. The break-point analysis for the concurrent period of record (1978–2011) identified 10 common periods of changes in hydrologic behavior at the selected gages. The water-level responses at each gage for the 10 periods displayed similarity in fluctuation patterns, highlighting the interconnectedness of the Florida Everglades hydrologic system. While the patterns were similar, the analysis also showed that larger fluctuations in water levels between periods occurred in Water Conservation Areas 2 and 3 in contrast to those in Water Conservation Area 1 and the Everglades National Park. Results from the analysis indicate that the cumulative Z-score curve, in conjunction with cumulative water-level frequency distribution curves, can be a useful tool in identifying and quantifying changes in historical hydrologic behavior within the Everglades. In addition to the analysis, a spreadsheet application was developed to assist in applying these techniques to time-series water-level data at gages within the Everglades and is included with this report.

Florida↗
USGS2026

dMODELS, a MATLAB software package for modeling crustal deformation near volcanic centers and active faults using Global Navigation Satellite System data—User guide

dMODELS is a MATLAB software package that implements the most common analytical models used to interpret deformation measurements near faults and active volcanic centers. This manual focuses on inversion of deformation data from the Global Navigation Satellite System (GNSS). The included case studies emphasize the GNSS inversion component of the software. Source models include pressurized spherical, spheroidal, and horizontal sill (penny-crack) magma reservoirs in a homogeneous, elastic, isotropic, flat half-space. A topography correction is available for the spherical source. Dikes and faults are described following the mathematical notation for the rectangular dislocations in a homogeneous, elastic, flat half-space. Equations have been reviewed for typographical errors present in the original literature and verified against finite-element method numerical models. GNSS data from the 2006 eruption at Augustine Volcano, Alaska; the 1998–2000 unrest at Taal Volcano, Philippines; and the 2009 earthquake in L’Aquila, Italy, are used to demonstrate the application of the software package.

Techniques and Methods↗
USGS2008

Isotope and Chemical Methods in Support of the U.S. Geological Survey Science Strategy, 2003-2008

Principal functions of the Mineral Resources Program are providing information to decision-makers related to mineral deposits on federal lands and predicting the environmental consequences of the mining or natural weathering of those deposits. Performing these functions requires that predictions be made of the likelihood of undiscovered deposits. The predictions are based on geologic and geoenvironmental models that are constructed for the various types of mineral deposits from detailed descriptions of actual deposits and detailed understanding of the processes that formed them. Over the past three decades the understanding of ore-forming processes has benefitted greatly from the integration of laboratory-based geochemical tools with field observations and other data sources. Under the aegis of the Evolution of Ore Deposits and Technology Transfer Project (EODTTP), a five-year effort that terminated in 2008, the Mineral Resources Program provided state-of-the-art analytical capabilities to support applications of several related geochemical tools.

Fact Sheet↗
USGS1996

A three-dimensional method-of-characteristics solute-transport model (MOC3D)

This report presents a model, MOC3D, that simulates three-dimensional solute transport in flowing ground water. The model computes changes in concentration of a single dissolved chemical constituent over time that are caused by advective transport, hydrodynamic dispersion (including both mechanical dispersion and diffusion), mixing (or dilution) from fluid sources, and mathematically simple chemical reactions (including linear sorption, which is represented by a retardation factor, and decay). The transport model is integrated with MODFLOW, a three-dimensional ground-water flow model that uses implicit finite-difference methods to solve the transient flow equation. MOC3D uses the method of characteristics to solve the transport equation on the basis of the hydraulic gradients computed with MODFLOW for a given time step. This implementation of the method of characteristics uses particle tracking to represent advective transport and explicit finite-difference methods to calculate the effects of other processes. However, the explicit procedure has several stability criteria that may limit the size of time increments for solving the transport equation; these are automatically determined by the program. For improved efficiency, the user can apply MOC3D to a subgrid of the primary MODFLOW grid that is used to solve the flow equation. However, the transport subgrid must have uniform grid spacing along rows and columns. The report includes a description of the theoretical basis of the model, a detailed description of input requirements and output options, and the results of model testing and evaluation. The model was evaluated for several problems for which exact analytical solutions are available and by benchmarking against other numerical codes for selected complex problems for which no exact solutions are available. These test results indicate that the model is very accurate for a wide range of conditions and yields minimal numerical dispersion for advection-dominated problems. Mass-balance errors are generally less than 10 percent, and tend to decrease and stabilize with time.

Water-Resources Investigations Report↗
USGS2009

Inference of distributional parameters from compositional samples containing nondetects

Low concentrations of elements in geochemical analyses have the peculiarity of being compositional data and, for a given level of significance, are likely to be beyond the capabilities of laboratories to distinguish between minute concentrations and complete absence, thus preventing laboratories from reporting extremely low concentrations of the analyte. Instead, what is reported is the detection limit, which is the minimum concentration that conclusively differentiates between presence and absence of the element. A spatially distributed exhaustive sample is employed in this study to generate unbiased sub-samples, which are further censored to observe the effect that different detection limits and sample sizes have on the inference of population distributions starting from geochemical analyses having specimens below detection limit (nondetects). The isometric logratio transformation is used to convert the compositional data in the simplex to samples in real space, thus allowing the practitioner to properly borrow from the large source of statistical techniques valid only in real space. The bootstrap method is used to numerically investigate the reliability of inferring several distributional parameters employing different forms of imputation for the censored data. The case study illustrates that, in general, best results are obtained when imputations are made using the distribution best fitting the readings above detection limit and exposes the problems of other more widely used practices. When the sample is spatially correlated, it is necessary to combine the bootstrap with stochastic simulation.

Conference Paper↗
USGS1971

Analytical sensitivities and energies of thermal neutron capture gamma rays II

A table of the analytical sensitivities of the principal lines in the thermal neutron capture gamma-ray spectrum from 0 to 3 MeV has been compiled for most of the elements. A tabulation of the full-energy, single-escape, and double-escape peaks has also been made according to energy. The tables are useful for spectral interpretation and calibration.

Nuclear Instruments and Methods↗
USGS2002

Determination of the common and rare alkalies in mineral analysis

Methods are described which afford a determination of each member of the alkali group and are successful in dealing with the quantities of the rare alkalies found in rocks and minerals. The procedures are relatively rapid and based chiefly on the use of chloroplatinic acid, absolute alcohol and ether, and ammonium sulfate. The percentages of all the alkalies found in a number of minerals are given.

Industrial And Engineering Chemistry Analytical Ed↗
USGS2006

U-Pb SHRIMP dating of uraniferous opals

U-Pb and U-series analyses of four U-rich opal samples using sensitive high-resolution ion microprobe (SHRIMP) demonstrate the potential of this technique for the dating of opals with ages ranging from several tens of thousand years to millions of years. The major advantages of the technique, compared to the conventional thermal ionisation mass spectrometry (TIMS), are the high spatial resolution (???20 ??m), the ability to analyse in situ all isotopes required to determine both U-Pb and U-series ages, and a relatively short analysis time which allows obtaining a growth rate of opal as a result of a single SHRIMP session. There are two major limitations to this method, determined by both current level of development of ion probes and understanding of ion sputtering processes. First, sufficient secondary ion beam intensities can only be obtained for opal samples with U concentrations in excess of ???20 ??g/g. However, this restriction still permits dating of a large variety of opals. Second, U-Pb ratios in all analyses drifted with time and were only weakly correlated with changes in other ratios (such as U/UO). This drift, which is difficult to correct for, remains the main factor currently limiting the precision and accuracy of the U-Pb SHRIMP opal ages. Nevertheless, an assumption of similar behaviour of standard and unknown opals under similar analytical conditions allowed successful determination of ages with precisions of ???10% for the samples investigated in this study. SHRIMP-based U-series and U-Pb ages are consistent with TIMS dating results of the same materials and known geological timeframes. ?? 2005 Elsevier B.V. All rights reserved.

Chemical Geology↗
USGS2022

An interactive viewer to improve operational aftershock forecasts

The U.S. Geological Survey (USGS) issues forecasts for aftershocks about 20 minutes after most earthquakes above M 5 in the United States and its territories, and updates these forecasts 75 times during the first year. Most of the forecasts are issued automatically, but some forecasts require manual intervention to maintain accuracy. It is important to identify the sequences whose forecasts will benefit from a modified approach so the USGS can provide accurate information to the public. The oaftools R package ( Paris and Michael, 2022 ) includes functions that analyze and plot earthquake sequences and their forecasts to identify which sequences require such intervention. The package includes the Operational Aftershock Forecast (OAF) Viewer, which incorporates the functions into an interactive web environment that can be used to explore aftershock sequences. The OAF Viewer starts with a global map and table of mainshocks. After a mainshock has been selected, the map and a new table show its aftershocks and the OAF Viewer generates five analytical plots: (1) magnitude–time, which is used to look for patterns in the data; (2) cumulative number, to see how the productivity of the sequence compares to a Reasenberg and Jones (1989) aftershock model over time; (3) magnitude–frequency, to compare the ratio of large to small magnitudes and extrapolate to higher magnitudes with sparse data and lower magnitudes with incomplete data; (4) forecast success, to compare the forecasts with observations for a sequence; and (5) parameter–time, which examines the temporal evolution of the forecast model parameters. The user can interact with the functions provided by the oaftools package through the OAF Viewer or by incorporating the functions into their own analysis methods. The OAF Viewer will help seismologists understand complexities in the data, communicate with the public and emergency managers, and improve the OAF system by maintaining operational awareness.

Seismological Research Letters↗
USGS2011

The effects of α-cellulose extraction and blue-stain fungus on retrospective studies of carbon and oxygen isotope variation in live and dead trees†

Tree-ring carbon and oxygen isotope ratios from live and recently dead trees may reveal important mechanisms of tree mortality. However, wood decay in dead trees may alter the δ 13 C and δ 18 O values of whole wood obscuring the isotopic signal associated with factors leading up to and including physiological death. We examined whole sapwood and α-cellulose from live and dead specimens of ponderosa pine ( Pinus ponderosa ), one-seed juniper ( Juniperous monosperma ), piñon pine ( Pinus edulis ) and white fir ( Abies concolor ), including those with fungal growth and beetle frass in the wood, to determine if α-cellulose extraction is necessary for the accurate interpretation of isotopic compositions in the dead trees. We found that the offset between the δ 13 C or δ 18 O values of α-cellulose and whole wood was the same for both live and dead trees across a large range of inter-annual and regional climate differences. The method of α-cellulose extraction, whether Leavitt-Danzer or Standard Brendel modified for small samples, imparts significant differences in the δ 13 C (up to 0.4‰) and δ 18 O (up to 1.2‰) of α-cellulose, as reported by other studies. There was no effect of beetle frass or blue-stain fungus ( Ophiostoma ) on the δ 13 C and δ 18 O of whole wood or α-cellulose. The relationships between whole wood and α-cellulose δ 13 C for ponderosa, piñon and juniper yielded slopes of ~1, while the relationship between δ 18 O of whole wood and α-cellulose was less clear. We conclude that there are few analytical or sampling obstacles to retrospective studies of isotopic patterns of tree mortality in forests of the western United States.

Rapid Communications in Mass Spectrometry↗
USGS2020

Coseismic and post-seismic gravity disturbance induced by seismic sources using a 2.5-D spectral element method

I present a prescription for computing free-air coseismic and post-seismic gravity changes induced by seismic sources in a viscoelastic earth model. I assume a spherical earth geometry and a 2.5-D calculation, that is, 3-D motions that satisfy the equations of quasi-static equilibrium on a 2-D viscoelastic structure. The prescription permits application to regional gravity computations where a 2-D structure adequately represents the structural heterogeneity. I use a hybrid approach where deformation is computed on a discretized domain and the resulting density perturbations are expanded with spherical harmonics to produce the free-air gravity field. Starting with a solution to the equations of quasi-static displacements in the Laplace transform domain for a given dislocation source, I solve Poisson’s equation using Lagrangian interpolation on spectral element nodes to compute the required deformation quantities that contribute to free-air gravity. A numerical inverse Laplace transform then yields time domain results. This methodology is tested with analytic solutions on a spherically stratified viscoelastic structure, then applied to evaluate the effect of a descending slab of relatively high viscosity on post-seismic gravity in a megathrust faulting setting.

Geophysical Journal International↗
USGS2011

Implementation and modification of a three-dimensional radiation stress formulation for surf zone and rip-current applications

Regional Ocean Modeling System (ROMS v 3.0), a three-dimensional numerical ocean model, was previously enhanced for shallow water applications by including wave-induced radiation stress forcing provided through coupling to wave propagation models (SWAN, REF/DIF). This enhancement made it suitable for surf zone applications as demonstrated using examples of obliquely incident waves on a planar beach and rip current formation in longshore bar trough morphology (Haas and Warner, 2009). In this contribution, we present an update to the coupled model which implements a wave roller model and also a modified method of the radiation stress term based on Mellor (2008, 2011a,b,in press) that includes a vertical distribution which better simulates non-conservative (i.e., wave breaking) processes and appears to be more appropriate for sigma coordinates in very shallow waters where wave breaking conditions dominate. The improvements of the modified model are shown through simulations of several cases that include: (a) obliquely incident spectral waves on a planar beach; (b) obliquely incident spectral waves on a natural barred beach (DUCK'94 experiment); (c) alongshore variable offshore wave forcing on a planar beach; (d) alongshore varying bathymetry with constant offshore wave forcing; and (e) nearshore barred morphology with rip-channels. Quantitative and qualitative comparisons to previous analytical, numerical, laboratory studies and field measurements show that the modified model replicates surf zone recirculation patterns (onshore drift at the surface and undertow at the bottom) more accurately than previous formulations based on radiation stress (Haas and Warner, 2009). The results of the model and test cases are further explored for identifying the forces operating in rip current development and the potential implication for sediment transport and rip channel development. Also, model analysis showed that rip current strength is higher when waves approach at angles of 5° to 10° in comparison to normally incident waves.

Coastal Engineering↗
USGS1999

Soil water retention and maximum capillary drive from saturation to oven dryness

This paper provides an alternative method to describe the water retention curve over a range of water contents from saturation to oven dryness. It makes two modifications to the standard Brooks and Corey [1964] (B-C) description, one at each end of the suction range. One expression proposed by Rossi and Nimmo [1994] is used in the high-suction range to a zero residual water content. (This Rossi-Nimmo modification to the Brooks-Corey model provides a more realistic description of the retention curve at low water contents.) Near zero suction the second modification eliminates the region where there is a change in suction with no change in water content. Tests on seven soil data sets, using three distinct analytical expressions for the high-, medium-, and low-suction ranges, show that the experimental water retention curves are well fitted by this composite procedure. The high-suction range of saturation contributes little to the maximum capillary drive, defined with a good approximation for a soil water and air system as H cM = ∫ 0 ∞ K rw dh c , where k rw is relative permeability (or conductivity) to water and h c is capillary suction, a positive quantity in unsaturated soils. As a result, the modification suggested to describe the high-suction range does not significantly affect the equivalence between Brooks-Corey (B-C) and van Genuchten [1980] parameters presented earlier. However, the shape of the retention curve near “natural saturation” has a significant impact on the value of the capillary drive. The estimate using the Brooks-Corey power law, extended to zero suction, will exceed that obtained with the new procedure by 25 to 30%. It is not possible to tell which procedure is appropriate. Tests on another data set, for which relative conductivity data are available, support the view of the authors that measurements of a retention curve coupled with a speculative curve of relative permeability as from a capillary model are not sufficient to accurately determine the (maximum) capillary drive. The capillary drive is a dynamic scalar, whereas the retention curve is of a static character. Only measurements of infiltration rates with time can determine the capillary drive with precision for a given soil.

Water Resources Research↗
USGS2007

Development and evaluation of a gas chromatographic method for the determination of triazine herbicides in natural water samples

A multi-residue method is described for the determination o triazine herbicides in natural water samples. The technique uses solvent extraction followed by gas chromatographic separation and detection employing nitrogen-selective devices. Seven compounds can be determined simultaneously at a nominal detection limit of 0.1 μg/L in a 1-litre sample. Three different natural water samples were used for error analysis via evaluation of recovery efficiencies and estimation of overall method precision. As an alternative to liquid-liquid partition (solvent extraction) for removal of compounds of interest from water, solid-phase extraction (SPE) techniques employing chromatographic grade silicas with chemically modified surfaces have been examined. SPE is found to provide rapid and efficient concentration with quantitative recovery of some triazine herbicides from natural water samples. Concentration factors of 500 to 1000 times are obtained readily by the SPE technique.

International Journal of Environmental Analytical ↗
USGS2012

WTAQ - A computer program for aquifer-test analysis of confined and unconfined aquifers

Computer program WTAQ was developed to implement a Laplace-transform analytical solution for axial-symmetric flow to a partially penetrating, finite-diameter well in a homogeneous and anisotropic unconfined (water-table) aquifer. The solution accounts for wellbore storage and skin effects at the pumped well, delayed response at an observation well, and delayed or instantaneous drainage from the unsaturated zone. For the particular case of zero drainage from the unsaturated zone, the solution simplifies to that of axial-symmetric flow in a confined aquifer. WTAQ calculates theoretical time-drawdown curves for the pumped well and observation wells and piezometers. The theoretical curves are used with measured time-drawdown data to estimate hydraulic parameters of confined or unconfined aquifers by graphical type-curve methods or by automatic parameter-estimation methods. Parameters that can be estimated are horizontal and vertical hydraulic conductivity, specific storage, and specific yield. A sample application illustrates use of WTAQ for estimating hydraulic parameters of a hypothetical, unconfined aquifer by type-curve methods. Copyright ASCE 2004.

Conference Paper↗
USGS1954

Spectrochemical determination of thorium in monazite by the powder-d.c. arc technique

Thorium in monazite is determined by a d.c. carbon-arc technique using zirconium as the internal standard. The analytical curve for Th II 2870.413 Å/Zr II 2844-579 Å is established by means of synthetic standards containing graduated amounts of thoria and 0.500 per cent zirconia in pegmatite base (60 parts quartz, 40 parts microchne, and 1 part ferric oxide). Monazite samples are diluted 14-fold with pegmatite base that contains 0.538 per cent ZrO 2 , so that the zirconia content of the resulting mixture is also 0.500 per cent. In addition, both the standards and the diluted monazites are mixed with one-half their weight of powdered graphite. Approximately 25 mg of the prepared samples are arced to completion at 15.5 to 17.5 amperes. With the 14-fold dilution employed, the accurate range of the method is 3 to 20 per cent thoria in the original monazite. The coefficient of variation for a single determination is 4 per cent at the 7 per cent thoria level. Tests with synthetic unknowns and chemically analyzed monazites show a maximum error of ±10 per cent of the thoria content. If niobium is substituted for zirconium as the internal standard, there is a loss of precision. Platinum as the internal standard gives results of good precision but introduces a marked sensitivity to matrix effects.

Spectrochimica Acta↗
USGS2011

Proceedings of the third USGS modeling conference, June 7-11, 2010, Broomfield, Colorado-Understanding and predicting for a changing world

The Third USGS Modeling Conference was held June 7th-11, 2010, in Broomfield, Colorado. The conference focused on the development and application of analytical and theoretical models and data availability that support managing the Nation's resources and help protect lives and property. Participants at the conference included scientists and managers from Department of the Interior (DOI) Bureaus; national and international Federal, State, and local agencies; academic institutions; and nongovernmental organizations. The conference was organized according to DOI priorities and the strategic directions of the USGS Science Strategy; the following themes were emphasized: (1) Understanding Ecosystems and Restoring America's Treasured Landscapes; (2) Climate Change and Impact; (3) New Energy Frontier and Minerals for America; (4) A National Hazards, Risk, and Resilience Assessment Program; (5) Role of Environment and Wildlife in Human Health; (6) A Water Census of the United States; and (7) New Methods of Investigation and Discovery. The conference theme-"Understanding and Predicting for a Changing World"-focused on the following goals: advance development and application of models; provide tools that address management issues; present state-of-the-art models ranging from individual phenomena to integrated systems; and foster a working community among scientists and managers.

Scientific Investigations Report↗