Synopsis of basin-scale 3D framework modeling at the USGS Mapping Program
No abstract available.
SEARCH · USGS Science
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.
No abstract available.
Redox conditions in the Numerous Sand Channels Zone beneath a petrochemical reclamation site in Harris County, Texas, range from sulfate reducing to methanogenic as indicated by the presence of methane in ground water and the range of molecular hydrogen concentrations. Assessment of the potential for reductive dechlorination using BIOCHLOR as a screening tool indicated conditions favoring anaerobic degradation of chlorinated organic compounds in the Numerous Sand Channels Zone. Evidence supporting reductive dechlorination includes apparently biogenic cis -1,2-dichloroethene; an increased ratio of 1,2-dichloroethane to 1,1,2-trichloroethane downgradient from the assumed contaminant source area; ethene and methane concentrations greater than background concentrations within the area of the contaminant plume; and a positive correlation of the ratio of ethene to vinyl chloride as a function of methane concentrations. The body of evidence presented in this report argues for hydrogenolysis of trichloroethene to cis -1,2-dichloroethene; of 1,1,2-trichloroethane to 1,2-dichloroethane; and of vinyl chloride to ethene within the Numerous Sand Channels Zone. Simulations using BIOCHLOR yielded apparent first-order decay constants for reductive dechlorination in the sequence tetrachloroethene trichloroethene cis -1,2-dichloroethene vinyl chloride ethene within the range of literature values reported for each compound and apparent first-order decay constants for reductive dechlorination in the sequence 1,1,2-trichloroethane 1,2-dichloroethane slightly greater than literature values reported for each compound along the upgradient segment of a simulated ground-water flowpath. Except for vinyl chloride, apparent rates of reductive dechlorination for all simulated species show a marked decrease along the downgradient segment of the simulated ground-water flowpath. Evidence for reductive dechlorination of chlorinated ethenes within the Numerous Sand Channels Zone indicates potential for natural attenuation of chlorinated ethenes. Reductive dechlorination of chlorinated ethanes apparently occurs to a lesser extent, indicating relatively less potential for natural attenuation of chlorinated ethanes. Additional data are needed on the concentrations and distribution of chlorinated ethenes and ethanes in individual fine sand intervals of the Numerous Sand Channels Zone. This information, combined with lower minimum reporting levels for future chloroethane analyses, might enable a more complete and quantitative assessment of the potential for natural attenuation at the site.
There is increasing interest in extracting critical minerals (CM), including rare earth elements (REE), from coals in the United States to address the overreliance on imported REE. The Gulf Coast and the Fort Union regions are the two major lignite-bearing regions within the United States. Recent REE and CM studies of the lignite in these regions have indicated that the coals may be a viable source material for REE and CM extraction. To evaluate in-place coal as a potential source of REE and CM, the coal resources need to be quantified. This study presents the results of a regional analysis of Gulf Coast Coal Region lignite and bituminous coal resources that might be available as potential sources of REE and CM. The resource analysis used kriging methods to develop isopleth maps of cumulative coal thickness throughout the region using data from drill holes and measured mine and outcrop sections for a total of 31,181 data points. The estimated total remaining coal resource in the Gulf Coast Coal Region is about 83 billion metric tons in the upper 90 meters (approximately 300 feet) of the subsurface. The current estimate of remaining coal resources includes approximately 1.88 billion metric tons of coal that have been extracted in the region since 1983. Texas accounted for 39 percent (32 billion metric tons) of the remaining total resource, followed by Mississippi (24 percent, 20 billion metric tons ), Louisiana (14 percent, 12 billion metric tons), Tennessee (10 percent, 8.5 billion metric tons), and Arkansas (6 percent, 5.1 billion metric tons). Alabama, Georgia, Kentucky, and Missouri each accounted for less than 5 percent of the total resource; Georgia had the smallest resource estimated at 7 million metric tons. Here we report the first known statewide lignite resource estimates for Georgia, Kentucky (820 million metric tons), and Missouri (1,800 million metric tons). A comparison of the results of this study with those of previous Gulf Coast Coal Region and Fort Union Coal Region resource studies is difficult because each study used different data sources, assessment methodologies, overburden depths, and qualifying coal thicknesses.
Groundwater is heavily used for public supply and industrial uses in the Baton Rouge, Louisiana, area. Lowered water levels resulting from groundwater withdrawals have induced the movement of saltwater towards wells in East Baton Rouge and West Baton Rouge Parishes. Saltwater intrusion has the potential to affect water supply infrastructure, reduce water availability for some uses, and increase treatment costs. To document current conditions, samples were collected from 161 wells screened in 10 aquifers of the Southern Hills regional aquifer system during November 2021 through February 2022. The results were compared with historical data to identify where chloride concentrations are increasing, which could indicate that saltwater intrusion is occurring. Saltwater intrusion, to varying degrees and areal extents, was observed in most of the 10 aquifers. The limited availability of monitoring wells near or within some of the known saltwater plume areas restricts tracking of the movement or delineation of the plumes’ current extents.
Abstract Groundwater flow between 2014 through 2040 was simulated in the New Jersey Coastal Plain based on three withdrawal scenarios. Two of the scenarios were based on projected population trends and the assumption of water conservation; the nominal water-loss scenario projected a status quo in the efficiency of water loss in the delivery systems whereas the optimal water-loss scenario projected a better water-loss efficiency resulting in less withdrawals. The third scenario assumes that all wells will withdraw water at their full allocation level which is generally much more than reported withdrawals in 2013 or projected under the other two scenarios. Maps and summaries of heads and drawdowns are presented for nine confined aquifers. All the aquifers have areas with heads below sea level by 2040. Of the three scenarios, the drawdowns are most extreme in the full allocation scenarios; there are large areas of head decline greater than 20 feet in 5 of the 9 confined aquifers. The exceptions are the Vincentown aquifer, despite some areas of large drawdown in the vicinity of wells, and the three Potomac-Raritan-Magothy (PRM) aquifers where withdrawals are regulated by Critical Area restrictions. The nominal and optimal water-loss scenarios have some areas of head declines; most are less than 15 feet. The simulation of these scenarios shows some extensive areas of head recovery as well—especially in the aquifers that are regulated by the Critical Area restrictions. Budgets of inflow and outflow components were calculated for 44 hydrologic budget areas (HBAs). The budget analysis shows that the water movement is complex and varies based on the aquifer geometry and location of pumping wells. Flow components between the unconfined and confined parts of the system were summarized by HUC11 (hydrologic unit code 11) basins.
Estimating the evolving state of stress along active fault systems can provide insight into the conditions that generated past ground‐rupturing earthquakes and influenced their ability to propagate through areas of geometric complexity, such as fault branches and stepovers. We use quasi‐static forward numerical models that incorporate the 3D complex configuration of active faults in southern California to estimate shear tractions on the geometrically complex southern San Andreas and San Jacinto faults from 1000 to 1900 C.E. These tractions include interseismic accumulation of traction due to tectonic loading, viscoelastic relaxation of shear stress within the upper crust between earthquakes, and effects of other earthquakes on the fault network. We simulate ground‐rupturing earthquakes based on the along‐strike earthquake extents modeled by Scharer and Yule (2020) , assuming that stress drop is complete in each earthquake. We use Monte Carlo simulations to estimate uncertainty in evolving shear tractions due to uncertainties in earthquake timing and in upper‐crustal viscosity. Pre‐earthquake shear tractions typically do not exceed ∼2 MPa. Although ruptures with length <200 km have pre‐earthquake shear tractions that range from near zero to ∼1.75 MPa, these tractions are not less than ∼0.4 MPa for earthquakes with rupture length >200 km. Earthquakes with long (>200 km) ruptures occur only in the single‐stranded part of the system, whereas those with short (<125 km) rupture length and high pre‐earthquake shear traction occur near fault stepovers and branches. This suggests that high accumulated shear traction encourages longer rupture propagation, but may not be sufficient to overcome geometric complexities. This modeling approach informs our understanding of rupture propagation and provides estimates of fault shear tractions that are unavailable from direct measurements.
We obtained 12 core samples for physical and chemical characterization from three serpentinite mud volcanoes (Yinazao, Asùt Tesoru, and Fantangisña) located on the forearc of the Mariana subduction system, that were drilled during International Ocean Discovery Program Expedition 366. Two samples from the Fantangisña mud volcano are interpreted to be clay-rich fault gouges derived from the subduction channel. Their bulk compositions are intermediate between the serpentinites and oceanic basalts. The oceanic crustal materials in the gouges have been thoroughly metasomatized and the serpentinites extensively altered to the trioctahedral, Mg-rich smectite clays saponite and corrensite. The only relict phases in clasts of crustal rock are accessory Ti- and P-bearing minerals. The two fault gouge samples have lower frictional strengths ( μ < 0.2) than the serpentinites ( μ = 0.2–0.4), and their measured permeabilities are also somewhat lower. Their physical and compositional properties correspond to saponitic gouges from other faults that juxtapose serpentinite against crustal rocks, in particular gouges from the two creeping traces of the San Andreas Fault recovered in the core from the San Andreas Fault Observatory at Depth. The décollement beneath Fantangisña mud volcano is thus expected to be very weak and likely characterized by stable slip.
The Basin and Range Province is host to many important low-sulfidation epithermal deposits. Within this broad zone of extension, epithermal deposits are hosted by specific areas of Miocene and younger bimodal volcanism. In northern Nevada, rifting and related volcanic activity occurred in response to thermal bulging during the development of the Yellowstone hotspot. The Colorado River Extensional Corridor in southern Nevada, adjoining eastern California and northwestern Arizona, is a major zone of Miocene crustal extension formed during the transformation of the western margin of North America from a convergent to a transform plate boundary. The style of mineralization in the low-sulfidation epithermal deposits in these areas is strongly controlled by the nature of the volcanic successions. Highgrade, low-tonnage deposits are commonly found in flow-dominated volcanic successions where ore deposition occurred through short-lived periods of fluid flashing. In contrast, low-grade, large-tonnage deposits are more commonly located in clastic-dominated successions where fluid infiltration of the permeable hosts and cooling predominated.
The Ma'millah gold mine is about 90 km southeast of At Ta'if at lat 21°03'N., long 41°18'E., in southwestern Saudi Arabia. The deposit, which was worked extensively by ancient miners, consists of several veins and massive lenses of quartz along a regional fault. The extensive ancient mine workings, large size of some quartz outcrops, and significant gold content of several samples indicated that subsurface exploration was warranted. Accordingly, two holes were drilled into the deposit during the interval December 1972 to February 1973. The drill findings were disappointing; the grade of vein material was found to be too low to be of economic interest and the quartz bodies were found to pinch out both along strike and downdip. The drilling program therefore was terminated. The deposit is not worth further exploration
Over 30 small, discontinuous, tabular carbonatite bodies are located in the Sheep Creek area, Ravalli County, southwest Montana. The age and origin of these REE-Nb-rich deposits are currently being investigated. The purpose of this paper is to document the occurrence of several rare minerals, including daqingshanite, törnebohmite, biraite, sahamalite, and ferriperbøeite, in two of the carbonatite bodies. These minerals are found in association with monazite, hydroxylbastnäsite, ferriallanite, calcite, dolomite, baryte, quartz, actinolite, apatite, celsian, and Sr-rich aragonite. Automated SEM-EDS was used to target the areas of interest in polished specimens for more detailed spot SEM-EDS and electron probe microanalysis. Raman spectra were also acquired for each of the rare minerals. The complex mineralogy of the Sheep Creek carbonatites is most likely due to several overlapping thermal events, including primary magmatic, overprinting hydrothermal, and supergene weathering stages. The rare minerals described in this study are believed to be hydrothermal and/or carbothermal in origin, although no estimates of temperature are available at this time.
Logistic regression was used to relate anthropogenic (man-made) and natural factors to the occurrence of elevated concentrations of nitrite plus nitrate as nitrogen in ground water in the Columbia Basin Ground Water Management Area, eastern Washington. Variables that were analyzed included well depth, depth of well casing, ground-water recharge rates, presence of canals, fertilizer application amounts, soils, surficial geology, and land-use types. The variables that best explain the occurrence of nitrate concentrations above 3 milligrams per liter in wells were the amount of fertilizer applied annually within a 2-kilometer radius of a well and the depth of the well casing; the variables that best explain the occurrence of nitrate above 10 milligrams per liter included the amount of fertilizer applied annually within a 3-kilometer radius of a well, the depth of the well casing, and the mean soil hydrologic group, which is a measure of soil infiltration rate. Based on the relations between these variables and elevated nitrate concentrations, models were developed using logistic regression that predict the probability that ground water will exceed a nitrate concentration of either 3 milligrams per liter or 10 milligrams per liter. Maps were produced that illustrate the predicted probability that ground-water nitrate concentrations will exceed 3 milligrams per liter or 10 milligrams per liter for wells cased to 78 feet below land surface (median casing depth) and the predicted depth to which wells would need to be cased in order to have an 80-percent probability of drawing water with a nitrate concentration below either 3 milligrams per liter or 10 milligrams per liter. Maps showing the predicted probability for the occurrence of elevated nitrate concentrations indicate that the irrigated agricultural regions are most at risk. The predicted depths to which wells need to be cased in order to have an 80-percent chance of obtaining low nitrate ground water exceed 600 feet in the irrigated agricultural regions, whereas wells in dryland agricultural areas generally need a casing in excess of 400 feet. The predicted depth to which wells need to be cased to have at least an 80-percent chance to draw water with a nitrate concentration less than 10 milligrams per liter generally did not exceed 800 feet, with a 200-foot casing depth typical of the majority of the area.
The project was to provide information on the quality of ground water from rural-domestic-water wells within the Osage Reservation and compare the water-quality to proximity to oil wells. About 38,500 oil wells have been drilled in the Reservation since drilling began in 1896. About 1,480 square miles or 64 percent of the Reservation is within a quarter mile of an oil well. The unconfined Quaternary sand aquifer covers about 315 square miles or about 14 percent of the Reservation and the confined Ada-Vamoosa sandstone aquifer covers about 800 square miles or about 35 percent of the Reservation. Fifty-eight percent of the Quaternary aquifer and 69 percent of the outcrop area of the Ada-Vamoosa aquifer are within a quarter mile of an oil well. One hundred twenty domestic ground-water wells were sampled from the Quaternary and Ada-Vamoosa aquifers. Forty-nine percent of the Reservation is underlain by the aquifers. Ground-water quality is good on most of the Reservation, but the use of domestic-water-supply wells tend to minimize water-quality problems. Existing water-supply wells commonly are located in areas that produce usable volumes of potable water. Several constituents in samples from the Ada-Vamoosa aquifer within a quarter mile of an oil well were significantly greater than from the aquifer not near oil wells. The constituents include specific conductance, dissolved solids, sodium, sulfate, chloride, bromide, and silica. These ions are probably derived from brine water. In the Ada-Vamoosa aquifer subgroups, 57 percent of the samples near oil wells and 24 percent of the samples not near oil wells had dissolved-solids concentrations greater than 500 milligrams per liter. The water quality in the Quaternary and Ada-Vamoosa aquifers is similar in areas where no oil wells have been drilled but is significantly different for several constituents. Median concentrations of major constituents from the Ada-Vamoosa aquifer not near oil wells were less than or equal to values from the Quaternary aquifer. Sixty-four percent of the water-quality samples from the Quaternary and 51 percent from the Ada-Vamoosa aquifers have dissolved-solids concentrations less than the secondary drinking water regulations of 500 milligrams per liter. Fifty-nine percent of the aquifer samples in the Quaternary aquifer subgroups not near oil wells and 70 percent of the samples near oil wells had dissolved solids less than 500 milligrams per liter. Areas in the Ada-Vamoosa aquifer near Hominy, Pershing, and Hula Lake have dissolved-solids concentrations greater than the secondary drinking water regulations. Water-quality samples from the Quaternary aquifer in these areas also have dissolved-solids concentrations greater than 500 milligrams per liter.
Executive Summary ModelMuse is a graphical user interface for several groundwater modeling programs. ModelMuse was updated to generate the input files for the parameter estimation software suite PEST. The software is used with MODFLOW or SUTRA models to run PEST-based parameter estimation and display the updated model inputs after parameter estimation. The PEST input files can also be used with the PEST++ version 5 software suite. Parameter estimation typically requires defining the parameters being adjusted during calibration and observations for assessing calibration quality. After a parameter is defined in ModelMuse, it can be applied to all or part of a model dataset. Pilot points—a parameterization device that facilitates higher levels of parameterization—can be used to assign spatially variable distributions of model inputs. Parameters can be applied to temporally varying features, such as boundary conditions, by either applying them to all the values in a series in one step or by applying separate parameters to individual members of a series. ModelMuse allows the definition of many observation types from various model output files. For MODFLOW 6 and SUTRA models, new options were added to ModelMuse to allow it to display the changed input after parameter estimation is complete. For MODFLOW–2005 and MODFLOW–NWT models, ModelMuse can import an entire model for visualization. An example illustrates the use of PEST with a MODFLOW 6 model in ModelMuse.
Earthquake emergency response and loss modeling require accurate models of ground‐motion intensity metrics (IMs), such as those produced by ShakeMap. In the critical minutes and hours after an earthquake, however, knowledge of the earthquake source geometry limits the accuracy of these models, especially in areas where seismic instrumentation is insufficient to constrain the shaking intensity pattern. Previous attempts to overcome this limitation have centered on computing the mean distance metrics from a suite of assumed source geometries. That approach is inadequate with modern ground‐motion models (GMMs) because they typically include terms for the hanging wall effect, depth to top of rupture, and other nonlinear terms. Such approaches also propagate the distance uncertainty through the GMMs and may not properly allocate the additional uncertainty into the between‐ and within‐event components. Here we present Finite‐Fault SIMulation with Mixed Effects Regression (FFSIMMER), a method to address these limitations by computing IM distributions from a suite of rupture realizations that fully accounts for all geometry‐related GMM terms and allows for the separation of the between‐ and within‐event components of the total uncertainty. The accuracy of this method is limited by the available constraints on the source parameters. To address this, FFSIMMER can incorporate constraints based on knowledge of the source, such as the approximate strike and dip based on knowledge of the plate tectonics and focal mechanism of the earthquake. We demonstrate the improvement in accuracy that can be achieved with the 2025 M 7.7 Mandalay, Burma (Myanmar), earthquake by incorporating constraints on the range of plausible source parameters.
In this chapter I discuss the use of so-called macroseismic data, i.e., reports of damage and other effects of shaking on humans and the built environment, to improve the characterization of earthquakes and the ground motions they produce. Macroseismic data are critical not only to investigate earthquakes that occurred before the start of the instrumental era in seismology, but are also, by virtue of the spatially rich nature of the data, of tremendous potential importance to investigate early instrumental and even recent earthquakes. Fueled by both the advent of on-line systems that now produce unprecedented volumes of macroseismic intensity data and by recognition of the importance of understanding key historical earthquakes, there has been a growing appreciation for the value of macroseismic data. It is, however, critical to understand both the unique challenges and the unique opportunities associated with these data. I review the evolution of thought and practices regarding analysis of macroseismic data, and the use of such data to improve the characterization of historic and early instrumental earthquakes.
In recent years, there have been several advancements related to the modelling of near-source effects of earthquake rupture on strong ground shaking, leading to an improved characterization of ground motions and resulting seismic hazard. Some of these modifications have stemmed from physics-based numerical modelling of the earthquake rupture process, using physics-based dynamic rupture simulations. These contributions have led to a better understanding of how fault rupture characteristics, geometry, and the style of faulting can interact with the hypocenter-dependence on the path from source to site that may ultimately guide the development of seismic directivity models. Moving forward, the application of modern techniques can be used to incorporate these source characteristics and near-fault ground motion behavior that contribute to the azimuthally varying effects that result in rupture directivity. One example is the application of machine learning methods to support more automated integration of new predictor variables in model development and open more evaluation opportunities to access residuals. Here, we utilize several techniques to take advantage of the plethora of synthetic data and its ability to supplement preexisting trends observed in data. We showcase two examples of how models can be either developed, expanded upon, or constrained using artificial neural network model (ANNs). We evaluate the performance of the ANN with existing methods, comparing misfit, potential limitations, and ability to continue to improve upon these methods in the future. One approach uses a set of simulations with corresponding synthetic ground motions from the Southern California Earthquake Center (SCEC) CyberShake study to develop a ground motion model adapted to incorporate seismic directivity information using an ANN. This large database (TBs) enables us to train the model to capture magnitude, period, and distance variations and how these parameters relate to amplification from hypocenters located along finite-faults. In some cases, there is reduced misfit from better representing source features that aren’t included in base ground motion models that neglect hypocenter location (e.g. azimuthal variation, source-to-site terms). Another ANN method uses a shallow-layered neural network model to better fit a hypocenter-independent model. This method adjusts the median and aleatory variability to account for the averaged impact of various hypocenter distributions to fit the underlying directivity adjustment model. This method serves as a template to apply to other directivity models, improving computational efficiency and more readily enabling integration in hazard codes.
We investigated sediment core records from the Cascadia Margin (Ocean Drilling Program Sites 1249 and 1252 at Hydrate Ridge; Integrated Ocean Drilling Program Site U1325 offshore Vancouver Island) using a Zr/Rb heavy mineral proxy from X-ray fluorescence (XRF) core scanning to identify intervals of primary detrital magnetic susceptibility (κ) and predict intervals where diagenesis caused magnetite dissolution by hydrogen sulfide. We also measured total sulfur (TS) content, grain size distributions, total organic carbon (TOC) content, and the magnetic mineral assemblage to further constrain the role of diagenesis on κ. Understanding how κ can be used to better characterize the varied effects of detrital and diagenetic signals in marine settings is important for understanding biogeochemical cycling and records of paleoenvironmental change. The upper 100 m of slope basin Site 1252 contains multiple intervals (> 90 m total) of decreased κ correlated with elevated TS content, consistent with dissolution of magnetite and precipitation of pyrite, iron monosulfides, and/or elemental sulfur. Similarly at the other slope basin site, Site U1325, κ is lower and TS is elevated in the interval between 24 and 51 mbsf, due to sulfide formation. At both slope basin sites, these low κ intervals correspond with high TOC, suggesting the possibility that organoclastic sulfate reduction (OSR) is likely a major driver of diagenetic alteration of κ at these sites. High TS:TOC ratios at Site U1325 suggests anaerobic oxidation of methane (AOM) during sulfate-methane transition zone (SMTZ) migration may have contributed to alteration of κ. In contrast, within the upper 90 m of Site 1249, a methane seep site at the summit of Hydrate Ridge, κ is almost entirely altered by diagenetic processes, with much of the low κ explained by a high degree of iron sulfide formation, while some intervals are affected by precipitation of magnetic iron sulfides that maintain or even increase κ. The presence of abundant methane seepage and gas hydrate as well as chemosynthetic seafloor fauna at this site, suggests that sulfide is released to the water column and AOM, rather than OSR, drives diagenetic alteration of κ at this site. Overall, the slope basin sites show episodic variation of κ that is influenced by TOC content, likely driven by changes in marine primary productivity and sedimentation rate, while the seep site shows consistently altered κ with lower TS content and no correlation with TOC. Methane seep environments likely experience loss of hydrogen sulfide to the water column and oxidation of hydrogen sulfide by seafloor seep fauna, which limits the amount of solid phase sulfur (pyrite, iron monosulfides, elemental sulfur) that can be precipitated within the sediments. In contrast, the migration of a buried SMTZ at slope basin sites results in enhanced sulfur precipitation within the sediments. This integrated magnetic and geochemical approach reveals the diagenetic production pathway and residence time of sulfide with the sediment column ultimately controls the style and degree of diagenetic loss of κ in marine sediments. This approach works best in environments with unaltered reference intervals/sites, consistent magnetic mineralogy, and clay-to-silt grain sizes.
No abstract available.