USGS Science⌕ Search

SEARCH · USGS Science

Results for “The Planetary Science Journal”

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.

1,665 records · Page 28Linked to original sources
USGS2025

Quaternary-active faults and the role of inherited structures in the Sacramento-San Joaquin Delta, western Central Valley, northern California

Seismic sources and their associated hazards within the Sacramento-San Joaquin Delta region of north-central California are relatively poorly characterized as compared to other, more heavily studied regions of northern California, such as the San Francisco Bay Area. Here we present a synthesis of subsurface, bedrock geology, and geodetic datasets from the Delta and from the Coast Ranges and Diablo Range to the northwest and southwest, respectively. We integrate these data and our own surface geologic and geomorphic observations to present a comprehensive review of faults in the Delta that exhibit Quaternary activity. Structural geologic data from the surrounding region highlight the significant influence that Late Cretaceous-to-Paleogene forearc structures exert on the geometry and kinematics of major Quaternary-active structures within the Delta. These inherited structures — including the Pittsburg-Kirby Hills Fault, Midland Fault, and Great Valley Fault System — exhibit a range of geometries and kinematics. Analysis of geomorphology along these structures suggests that these structures combine to accommodate Quaternary strain across the Delta region. A clearer understanding of subsurface geometries and structural relationships, built upon the regional tectonic history, provides insight into modern deformation accommodated on older structures and helps inform interpretations of seismic hazard within the Delta.

California↗
USGS2026

Origin of the Pd/Pt ratio of the J-M Reef, Stillwater Complex, Montana, USA

The J-M reef of the Stillwater Complex is characterized by a high Pd/Pt ratio (mean ~3.8 with a standard error of 0.03) with a homogeneous geospatial distribution at the deposit scale. In this contribution, we demonstrate that the Pd/Pt ratio of the reef is the product of equilibration of an immiscible sulfide liquid with a silicate melt rich in Pd relative to Pt. Despite the high tenors of the J-M reef sulfides (avg 2,700 ppm Pt and 770 ppm Pt), numerical modeling shows that the parental melts did not have extraordinary Pd and Pt concentrations. Instead, the initial composition of a plausible parental silicate melt can have Pd and Pt contents well within the expected range of a normal, mantle-derived partial melt (i.e., ~10–20 ppb for both Pd and Pt with Pd/Pt of ~1). The relative differences in the partitioning behavior of Pt and Pd between sulfide liquid and silicate melt are unlikely to produce a consistent Pd/Pt ratio across a wide range of silicate melt to sulfide liquid mass ratios (i.e., R factors). Instead, the pre-emplacement fractionation of Pt alloy from S-undersaturated silicate magma accounts for the homogeneous and high Pd/Pt ratio of the J-M reef. We show that batch equilibration of sulfide liquid with silicate melt can produce the high Pd/Pt ratios of the reef if the partition coefficients between sulfide liquid and silicate melt for Pd and Pt are extremely high (>10 6 ). In an alternative model, Pd enrichment could be achieved by sulfide upgrading in resident footwall mush even if the partition coefficients between sulfide liquid and silicate melt are relatively small (between 10 4 and 10 6 ) because the instantaneous mass ratio of silicate melt to sulfide liquid is small (R ≈ 100–700), so the partitioning behavior of Pt and Pd has little impact on the composition of sulfide liquid.

Montana↗
USGS2026

Tephra from Kīlauea’s 2008–2018 lava lake eruption—Proximal deposits and dispersal characteristics

A network of ten buckets was established early in the 2008–2018 summit eruption at Kīlauea to collect proximal tephra ejected from the new, informally named the “Overlook crater”; the buckets were emptied on most days of the eruption thereafter. This report summarizes the results of more than 2,400 different sampling intervals (most 1–3 days long) during the eruption, focusing on the physical and dispersal characteristics of the tephra deposits. The network was within about 300 meters south of the vent to capture tephra dispersed by the dominant northeast trade wind. The juvenile tephra mainly reflected spattering at the southeast (SE) sink, a downwelling area in the southeastern part of the lava lake in the Overlook crater that remained in the same area throughout the eruption, with admixtures of solid rock and secondary minerals derived from the wall of the crater. The proportion of juvenile material to lithic material ranged widely early in the eruption but was generally greater than 90 percent for the last 6 years of the eruption as lake level rose and the crater walls decreased in height and became more stable. The accumulation rate of tephra at each bucket was strongly dependent on the location of the bucket and reflects the interplay between lava lake level and wind direction and speed. The mass per unit area (m/a) of collected tephra was a maximum of about 97 kilograms per square meter for the entire eruption, equivalent to a thickness of about 75 millimeters. Thirty-two explosive events with a volcanic explosivity index of −2 to −4 deposited much of the tephra in the network. Fifteen of these rock-fall-induced events occurred within 6 days of one another. Pele’s hair, a volcanic glass predominantly associated with quieter activity at the vent in the final half of the eruption, was dispersed more than 60 kilometers downwind from the lava lake and formed a nearly continuous deposit near the Overlook crater. This eruption was probably the most frequently sampled long-lasting eruption in history, but most of the deposits are ephemeral. The collected samples, although generally of small mass, are retained by the Hawaiian Volcano Observatory and are available for detailed study.

Hawaii↗
USGS2024

Onset of aftershocks: Constraints on the Rate-and-State model

Aftershock rates typically decay with time t after the mainshock according to the Omori–Utsu law, R (t)=K(c+t) −p ⁠ , with parameters K , c , and p . The rate‐and‐state (RS) model, which is currently the most popular physics‐based seismicity model, also predicts an Omori–Utsu decay with p = 1 and a c ‐value that depends on the size of the coseismic stress change. Because the mainshock‐induced stresses strongly vary in space, the c ‐value should vary accordingly. Short‐time aftershock incompleteness (STAI) in earthquake catalogs has prevented a detailed test of this prediction so far, but the newly developed a ‐positive method for reconstructing the true earthquake rate now allows its testing. Using previously published slip models, we calculate the coseismic stress changes for the six largest mainshocks in Southern California in recent decades and estimate the maximum shear as a scalar proxy of the coseismic stress tensor. Aftershock rates reconstructed for events in different stress ranges show that the rates follow a power law with p = 1 independent of stress with no clear sign of a c ‐value. The onset of the power‐law decay is abrupt and more delayed in areas with smaller stress changes. The observations do not necessarily contradict the RS model, as STAI limits the resolution for early aftershocks, and the RS model can reproduce the observations for specific Aσ values. However, the observations lead to strong constraints, namely Aσ <10 kPa and a power‐law decay of the background rate with distance to the fault, with exponent 2.7.

Seismological Research Letters↗
USGS2024

Surficial geology and Quaternary fault map of the Las Vegas Valley, Clark County, Nevada

The Surficial Geology and Quaternary Fault Map of the Las Vegas Valley, Clark County, Nevada is a 1:50,000-scale compilation of published 1:24,000-scale geologic maps integrated with new field and desktop mapping. This geologic map compilation and GIS database are part of a broader study on the Quaternary faults in the Las Vegas Valley, the results of which are summarized in this report. The map compilation utilized a variety of digital base maps to evaluate the characteristics of surficial deposits and Quaternary fault scarps including lidar data in undeveloped parts of the valley and orthorectified historical aerial photos in the urbanized center of the valley. The map distinguishes twenty surficial map units including alluvial-fan and ephemeral-wash deposits, groundwater discharge deposits of the Las Vegas Formation, and alluvium deposits of Las Vegas Wash. Quaternary faults include the Quaternary Las Vegas Valley fault system in the center of the valley and the Frenchman Mountain fault system along the eastern margin. All Quaternary fault traces in the valley were evaluated, and the previously published fault mapping was modified as needed. The new mapping and uniform correlation of surficial deposits displaced by Quaternary faults yields a consistent valley-wide characterization of the recency of fault activity. The map and GIS database include a Quaternary geochronology compilation of 119 previously published surficial deposit ages of various types, 35 new luminescence ages, and 3 new radiocarbon ages. The new age data further bracket the ages of the Las Vegas basin Quaternary stratigraphy and provide new constraints on the timing of Quaternary fault activity.

Nevada↗
USGS2025

To heal or not to heal?: 2. The moment-recurrence time behavior of repeating earthquakes in the 2011 Prague, Oklahoma aftershock sequence is consistent with laboratory healing rates

The timing and failure conditions of an earthquake are governed by the interplay between fault reloading and restrengthening. The moment-recurrence time behavior of repeating earthquakes can give observational estimates of fault healing rates; however, it is difficult to link these observed healing rates to laboratory studies of frictional healing in part because of uncertainty in lithology. Here, we study the 2011 Prague earthquake sequence, which includes repeating earthquakes in the Arbuckle group and the granitic basement, and compare them to laboratory experiments on samples of the Arbuckle and Troy granite (representative of the basement rock) (Okamoto et al., 2025, https://doi.org/10.1029/2024JB030573 ). We find three spatially distinct groups of repeating earthquakes with different moment-recurrence behavior: (a) constant moment-recurrence time in the Arbuckle group, (b) scattered moment-recurrence time at the intersection of the foreshock-mainshock fault in the granitic basement, and (c) moment-predictable behavior outside of the foreshock-mainshock fault intersection also in the granitic basement. Our observation of stagnant healing for repeating sequences in the Arbuckle group is consistent with laboratory observations of low healing rates for moderately high pore fluid pressures in Arbuckle samples. For the moment-predictable group, the source radius that is required in order to match healing rates is consistent with source radius estimations when taking into account reasonable attenuation of the P- pulse width. Overall, we observe diverse healing behaviors in the seismic families that are consistent with laboratory healing rates, providing seismic evidence that contact-scale frictional mechanisms are relevant to large-scale earthquake dynamics.

Oklahoma↗
USGS2024

Simulation of groundwater flow and brine discharge to the Dolores River in the Paradox Valley, Montrose County, Colorado

Salinity, or total dissolved solids (TDS), of the Colorado River affects agricultural, municipal, and industrial water users and is an important concern in the Western United States. In the Paradox Valley of southwestern Colorado, natural discharge of sodium-chloride brine to the Dolores River from the underlying core of a salt-valley anticline accounts for about 6 percent of the salinity load to the Colorado River. Formation of the Paradox Valley began during the Miocene, and subsequent erosion exposed the Pennsylvania Paradox Formation in the core of the anticline where a cap rock, collapse features, breccia, and sodium-chloride saturated brine developed at the top of the exposed salt diapir. The discharge of brine to the Dolores River is affected by these dissolution features, along with seasonal hydrologic conditions and density-dependent flow between older dense brine and the younger fresh groundwater in the overlying alluvial aquifer. To reduce TDS concentrations in the Dolores River through the Paradox Valley, the Bureau of Reclamation has pumped brine from a series of shallow wells adjacent to the river since July 1996. The pumped brine is collected and piped to a deep disposal well where it is injected into the Mississippian Leadville Limestone at a depth of about 4,570-meters below land surface. The pumping and injection operation is collectively known as the Paradox Valley Unit (PVU), and by 2015, the PVU had substantially reduced TDS concentrations in the Dolores River by about 70 percent. Since 2019, injection-pressure limits and related seismic activity have constrained deep-well injection and thus brine pumping at the PVU. In cooperation with the Bureau of Reclamation, the U.S. Geological Survey developed a MODFLOW-6 three-dimensional, variable-density groundwater flow and TDS transport model of the Paradox Valley to evaluate the effects of PVU pumping operations on brine discharge to the Dolores River and to guide additional research. The finite-difference model grid consists of 76 rows and 48 columns oriented from northwest to southeast in alignment with valley topography and groundwater-flow directions in the near-surface freshwater alluvial aquifer. A 7-layer hydrogeologic framework was developed from existing datasets to represent the alluvial aquifer, cap rock, collapse breccia, and groundwater flow and TDS transport from the underlying Paradox Formation salt to the Dolores River. The model represents a 33-year transient calibration period from 1987 through 2020 that includes pre-PVU conditions from 1987 through June 1996 and post-PVU conditions from July 1996 through 2020. A 1,000-year simulation of groundwater flow and coupled TDS transport computed the initial conditions for the subsequent 33-year transient simulation. Observations of precipitation, streamflow, evaporation, agricultural land use, and PVU brine pumping rates were used to specify appropriate boundary conditions to the model representing time-varying recharge, tributary streamflow, groundwater underflow, evapotranspiration (ET), and PVU pumping. Values for average monthly streamflow and TDS concentration at the upstream streamgage, the Dolores River at Bedrock (USGS streamgage 09169500), were specified as model input where the Dolores River enters Paradox Valley. Observed pumping from the PVU, water levels and TDS concentrations in groundwater, and streamflow and estimated TDS concentrations at the downstream streamgage, the Dolores River near Bedrock (USGS streamgage 09171100), were calibration targets that constrained the manual calibration of model parameters representing aquifer hydraulic conductivity, storage, streambed conductance, recharge, and (ET). Two primary model-calibration targets were the match between observed and simulated TDS mass flux from PVU pumping wells and the match between estimated and simulated TDS mass flux to the Dolores River. The simulated TDS mass withdrawn by pumping wells is calculated by the model as the product of the assigned pumping rate and simulated groundwater TDS concentrations. Because actual pumping rates were assigned as simulated values, the total simulated PVU pumping for the 33-year calibration is within 0.5 percent of the observed values. However, simulated concentrations and thus mass flux of TDS withdrawn by the PVU pumping wells were consistently about 26 percent less than observed values for all the simulated time periods (33-year simulation, pre-PVU, and post-PVU). The representation of brine inflow was explored through additional modeling to evaluate the effect of the simulated brine source on groundwater TDS concentrations. Results indicated that a saturated-salt constant-flux brine source best replicated the magnitude and transient pattern observed for TDS mass flux from PVU pumping wells. The simulated TDS mass flux to the Dolores River is compared to estimates based on observed streamflow and specific conductance (SC) data for the downstream streamgage. The calibrated model provided a close fit of simulated to measured streamflow at the downstream streamgage, and the calibrated model fit to estimated TDS concentrations at the downstream streamgage was reasonable. The greatest differences between simulated and estimated values occurred during drought periods from June 2000 to March 2003, May 2012 to June 2013, and October 2013 to October 2014, when simulated TDS concentrations in the river were greater than estimated concentrations. In general, simulated TDS mass flux to the river for the pre-PVU period is in good agreement with estimated values (2-percent difference), but the model overestimated TDS mass flux to the river by about 41 percent during the post-PVU period. The model uncertainty with respect to TDS mass flux to the river indicates other processes or model parameters not well represented by the model are affecting the system, especially during drought. During model calibration, the most sensitive parameters were identified as vertical hydraulic conductivity of the alluvial aquifer, conductance of the Dolores River streambed, ET extinction depth and rate, and recharge rate. Five 5-year scenarios of conditions for 2021–25 were simulated to assist evaluation of alternative strategies to manage the discharge of brine into the Dolores River. The first scenario simulates no PVU pumping and serves as a base case for comparison to the other scenarios. Two scenarios simulate the effects of varying withdrawal timing at an annual rate about one-third less than during 2010 through 2018. During high-flow spring snowmelt runoff periods when brine discharge is naturally minimized, PVU pumping does not substantially affect salinity in the Dolores River, and comparison of these two scenarios indicates that scheduling brine withdrawals during times of low river stage is nearly as effective at reducing TDS mass flux to the river as pumping brine year-round. Cessation of pumping during periods of high river stage may be advantageous for system maintenance, brine injection, and seismic-risk reduction. The fourth scenario tested the effect of reducing irrigation-return flow on brine discharge and predicted a slight reduction of TDS mass flux to the Dolores River, but not as great a reduction as that of using the PVU to remove brine. The fifth scenario simulated 5 years of drought conditions without PVU pumping and indicates brine discharge during drought about 15 percent greater than during average hydrologic conditions. Results from scenario 5 are consistent with the calibrated model results and indicate that aquifer properties and ET processes and parameters may be affecting simulation results during drought. The Paradox Valley groundwater model provides a reasonable overall match to observed conditions in the Dolores River. The model is useful for evaluating relative differences between brine management scenarios to inform PVU operational decisions and to identify gaps in data and process understanding. Representation of the brine source, hydraulic-conductivity parameters, and recharge and ET processes were identified as potential areas for additional field and modeling research. Additional research in the Paradox Valley might include field-data collection that provides additional information on the hydrogeologic framework, groundwater levels, groundwater TDS concentrations, stream characteristics, and aquifer properties. Additional modeling efforts could benefit from applying advanced tools for model development, calibration, and visualization including parameter-estimation and sensitivity analysis. Statistical evaluation of known model uncertainties such as hydraulic conductivity, streambed conductance, representations of the brine source, recharge, and ET could improve the match between simulated and estimated TDS mass flux from PVU pumping wells and to the Dolores River further informing model predictions and system understanding for the Paradox Valley.

Colorado↗
USGS2025

Magnitude, depth and methodological variations of spectral stress drop within the SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest Earthquake Sequence

We present the first ensemble analysis of the 56 different sets of results submitted to the ongoing Community Stress Drop Validation Study using the 2019 Ridgecrest, California, earthquake sequence. Different assumptions and methods result in different estimation of the source contribution to recorded seismograms, and hence to the source parameters (principally corner frequency, f c ⁠ , spectral stress drop, Δσ, and seismic moment, M 0 ⁠ ) obtained from modeling calculated source spectra. For earthquakes smaller than magnitude (M) 2.5 there is negligible correlation between the f c values obtained by different studies, implying that no present method is reliable using available data. For larger magnitude events, correlation between f c measurements of different studies, within even a small M range is always higher than spectral ⁠Δσ , because the f c measurements simply reflect the underlying physical decrease in f c with increasing M. We model the observed trends of submitted f c with both magnitude and depth. Most methods report an increase in spectral Δσ with M, although a magnitude‐invariant spectral Δσ is within the confidence limits. The depth dependence is smaller and depends on whether a study allows attenuation to vary with source depth; a combination of depth‐dependent attenuation correction, and depth‐dependent shear‐wave velocity can compensate for reported depth trends. We model the submitted values to remove differing M and depth variation to investigate the relative interevent variability. We find consistent relative variation between individual events, and also lower relative spectral Δσ in the northwest of the aftershock sequence, and higher on the cross fault and in the region of main fault intersection. This large‐scale comparison implies that absolute spectral Δσ estimates are dependent on the methods used; studies of different regions or using different methods should not be directly compared and improved constraints on path and site corrections are needed to resolve these absolute spectral Δσ differences.

California↗
USGS2025

Analysis of the potential effects of Uzbekistan’s mineral endowment on the critical mineral supply of tungsten

Tungsten appears on the 2018 and 2022 U.S. Geological Survey critical mineral lists in part because of a very high global production concentration in China, which produces almost 83 percent of the world’s mined tungsten. Using known parameters and values from other tungsten mining operations, we created hypothetical scenarios in which three tungsten deposits in Uzbekistan are considered for development. Our results show that all three deposits are likely to be economically viable to develop under 2024 market conditions. If the three studied tungsten deposits were put into production, Uzbekistan could become the third-leading tungsten-producing country in the world and increase world output of tungsten by 2.7 percent. Putting these tungsten deposits in Uzbekistan into production could slightly reduce the tungsten global market concentration, therefore reducing the supply disruption potential for tungsten.

Open-File Report↗
USGS2025

Flood-Inundation Maps of the Current and Jacks Fork Rivers including the Ozark National Scenic Riverways, Southeast Missouri, 2023

Digital flood-inundation maps for a 131.8-mile reach of the Current River and a 44.6-mile reach of the Jacks Fork River, in southeast Missouri, were created by the U.S. Geological Survey (USGS) in cooperation with the Ozark Foothills Regional Planning Commission and the South Central Ozark Council of Governments. The maps also encompass the 134 miles of the Current and Jacks Fork Rivers within the Ozark National Scenic Riverways, which is the first national park area to protect a river system. The flood-inundation maps, which can be accessed through the USGS Flood Inundation Mapping Program website at https://water.usgs.gov/osw/flood_inundation/ , depict estimates of the areal extent and depth of flooding corresponding to selected water levels (hereafter referred to as “stages”) at eight reference USGS streamgages—five on the Current River (USGS station numbers 07064440, 07064533, 07066510, 07067000, and 07068000) and three on the Jacks Fork River (USGS station numbers 07065200, 07065495, and 07066000). Near-real-time stages at these streamgages may be obtained from the USGS National Water Information System at https://doi.org/10.5066/F7P55KJN or the National Weather Service National Water Prediction Service at http://water.noaa.gov/ , which also forecasts flood hydrographs at four of these sites (USGS station numbers 07067000, 07068000, 07065495, and 07066000). Flood profiles were computed for seven of the eight map reaches by means of two-dimensional hydraulic models and the remaining reach by a one-dimensional hydraulic model. The models were calibrated by using stage-streamflow relations or streamflow measurements at the USGS streamgages and from high-flow stage measurements from water-level loggers distributed throughout the reaches. The hydraulic models were used to compute water-surface profiles for flood stages at 1-foot intervals referenced to the streamgage datums. The profile stages ranged from the National Weather Service “action stage” or near bankfull, to a stage exceeding the highest recorded water level at each streamgage. The simulated water-surface profiles were then combined with a digital elevation model (derived from light detection and ranging data having a nonvegetated vertical accuracy of a maximum 10-centimeter root mean square error) to delineate the area flooded at each water level and the associated water depths. The availability of these maps, along with information regarding current stage from the USGS streamgage and forecasted high-flow stages from the National Weather Service, will provide emergency management personnel, resource managers, and residents with information that is critical for flood-response activities such as evacuations and road closures, as well as for postflood recovery efforts.

Missouri↗
USGS2025

The role of alkali bicarbonate-sulfate brines in the genesis of carbonatite REE resources at the Bear Lodge Alkaline Complex, Wyoming

Rare-earth element (REE) resources in the Bear Lodge Alkaline Complex, Wyoming, are hosted in a variably leached carbonatite dike swarm spatially related to bodies of diatreme breccia. This study examines fluid inclusions in carbonatite dikes, peripheral fluorite breccias, and smoky quartz veins to reconstruct the physiochemical conditions of REE mineralization. Results reveal a multi-stage evolution of carbonatitic fluids: 1) high-temperature (330–432°C), CO 2 -rich magmatic fluids, 2) REE-enriched alkali bicarbonate-sulfate brines, and 3) low-temperature meteoric (117–182°C) or diluted magmatic fluids. Multiphase inclusions (stage-2) contain burbankite, nahcolite, and alkali sulfate daughter crystals, linking the brines to early burbankite [(Na,Ca) 3 (Sr,Ba,REE) 3 (CO 3 ) 5 ] mineralization in carbonatite. REE mobility has been linked to the stability of aqueous complexes with ligands such as Cl-, SO 4 2- , and HCO 3 - . However, light REEs were not initially transported significant distances, rather they crystallized as burbankite from brines or brine-melts associated with the carbonatite magma. Crystallization of burbankite removed Na from the system, leaving residual fluids with high K/Na ratios and promoting potassic fenitization throughout the complex. Fractionation of light and heavy REEs and involvement of low temperature (<300°C), high K/Na bicarbonate brines may have facilitated transport and formation of peripheral occurrences that are more enriched in heavy REEs at Bear Lodge.

Wyoming↗
USGS2025

Multi-scale geophysical mapping of the brine and bedrock surfaces along the Dolores River, Paradox Valley, Colorado, December 2023

Total dissolved solids derived from salt dome–sourced brine in the underlying alluvial aquifer substantially increase with distance in the reach of the Dolores River that passes through Paradox Valley in southwestern Colorado. The area has been the site of salinity control operations since the 1990s to reduce salt loading to the downstream Colorado River. Previous airborne and ground/water-based electromagnetic (EM) geophysical data have successfully mapped the top of the brine surface, albeit with relatively coarse near-surface resolution and limited spatial coverage. This present December 2023 study used ground-based high-resolution EM and passive seismic (horizontal-to-vertical spectral ratio, HVSR) tools to map in detail the depth and thickness of the brine zone in the alluvial aquifer (top of the brine down to bedrock contact) in areas immediately surrounding the Dolores River where previous airborne EM (AEM) results indicated brine within 10 m of land surface. Results indicate the deepest bedrock is generally associated with the shallowest brine and local depressions in the collapse breccia (caprock to the Paradox Formation salt) may facilitate vertical migration of brine into the alluvial aquifer. Additionally, the ground-based EM mapping corroborated general patterns in depth to brine that were observed in previous AEM results while also revealing additional detail, including suspected focused brine discharge zones to the Dolores River. A river-based survey wherein EM data, channel depth and river water–specific conductance information were collected augmented these findings and indicated specific areas on both the western and eastern sides of the river where focused brine discharge may occur. This study comprises a large-scale, ground- and water-based geophysical mapping effort, including hundreds of HVSR soundings and 100s of kilometres of EM data, which were successfully translated into digital brine and bedrock surfaces that could be incorporated into groundwater modelling efforts, future well siting or other decision-making.

Colorado↗
USGS2024

3-D geological modeling for numerical flow simulation studies of gas hydrate reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. This study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-step gridding were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a 3-D trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

Alaska↗
USGS2024

New occurrences of the rare, REE minerals daqingshanite, törnebohmite, biraite, sahamalite, and ferriperbøeite from the Sheep Creek area, Montana, USA

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.

Montana↗
USGS1999

Comparison of two approaches for determining ground-water discharge and pumpage in the lower Arkansas River Basin, Colorado, 1997-98

In March 1994, the Colorado Division of Water Resources (CDWR) adopted “Rules Governing the Measurement of Tributary Ground Water Diversions Located in the Arkansas River Basin” (Office of the State Engineer, 1994); these initial rules were amended in February 1996 (Office of the State Engineer, 1996). The amended rules require users of wells that divert tributary ground water to annually report the water pumped monthly by each well. The rules allow a well owner to report the pumpage measured by a totalizing flowmeter (TFM) or pumpage determined from electrical power data and a power conversion coefficient (PCC) (Hurr and Litke, 1989). Opinions by representatives of the State of Kansas, presented before the Special Master hearing a court case [State of Kansas v. State of Colorado, No. 105 Original (1996)] concerning post-Compact well pumping, stated that the PCC approach does not provide the same level of accuracy and reliability as a TFM when used to determine pumpage. In 1997, the U.S. Geological Survey (USGS), in cooperation with the CDWR, began a 2-year study to compare ground-water pumpage estimates made using the TFM and the PCC approaches. The study area was along the Arkansas River between Pueblo, Colorado, and the Colorado-Kansas State line (fig. 1). The two approaches for estimating ground-water discharge and pumpage were compared for more than 100 wells completed in the alluvial aquifer of the Arkansas River Basin. The TFM approach uses an inline flowmeter to directly measure instantaneous discharge and the total volume of water pumped at a well. The PCC approach uses electrical power consumption records and a power conversion coefficient to estimate the pumpage at ground-water wells. This executive summary describes the results of the comparison of the two approaches. Specifically, (1) the differences in instantaneous discharge measured with three portable flowmeters and measured with an inline TFM are evaluated, and the statistical differences in paired instantaneous discharge between the two approaches are determined; (2) short- and long-term variations in the PCC’s are presented; (3) differences in pumpage between the two approaches are evaluated, and the statistical differences in pumpage between the two approaches are determined; (4) potential sources of discrepancy between pumpage estimates are discussed; and (5) differences in total network pumpage using the two approaches are presented. During the irrigation seasons of 1997 and 1998, instantaneous discharge and electrical power demand were measured at randomly selected wells to determine PCC’s. At more than 100 wells, the PCC’s determined during the 1998 season were applied to total electrical power consumption data that was recorded between the initial and final readings at each network well site in 1998 to estimate total ground-water pumpage. At each site, an inline TFM was installed in a full-flowing, acceptable test section of pipe on the discharge side of the pump where the measurement of discharge was made. Measurements of instantaneous ground-water discharge also were made using three different types of portable flowmeters. The average velocity multiplied by the cross-sectional area of the discharge pipe was used to compute the discharge in gallons per minute. Whenever possible, discharge measurements were made at each network site using all three types of portable flowmeters.

Colorado↗
USGS2024

Temporal, spatial, and chemical evolution of Quaternary high-silica rhyolites in the Mineral Mountains, Utah

The Mineral Mountains in southwestern Utah are a structurally controlled core complex at the confluence of the Colorado Plateau and the Basin and Range physiographic provinces. Aside from hosting Utah’s largest batholith, the Mineral Mountains host some of the State’s youngest high-silica rhyolites, which have been linked to a magma source that is presently being utilized as an enhanced geothermal system. The high-silica rhyolites take the form of effusive lavas and domes, and explosive products are rare. Previous K-Ar dating of these Pleistocene rhyolites placed eruptions between about 790 and 500 kilo-annum (ka) with contemporaneous basalts erupting in the valley to the east of the Mineral Mountains. Large uncertainties on these ages obscured the tempo of eruptions and thus hindered attempts to constrain the timescales of the petrogenetic processes that produced the rhyolites. In this study, we build on previous studies conducted in the 1970s and 1980s by using new geochronologic and geochemical data to investigate the temporal and spatial evolution of the youngest phase of volcanism in the Mineral Mountains. We identify two major eruptive periods, from approximately 850 to 750 ka and from approximately 590 to 480 ka. The older phase is characterized by the eruption of several basaltic lavas, two obsidian flows, and a series of coalescing porphyritic rhyolite domes. The younger phase included the eruption of six evolved high-silica rhyolite domes and one pyroclastic deposit, followed by the eruption of trachyandesite in the adjacent valley to the east. Whole-rock geochemical data indicate that the rhyolites can be divided into three chemical groups, with more evolved compositions erupting through time. The youngest rhyolites along the range crest have the lowest total iron and TiO 2 concentrations and the highest incompatible element concentrations, indicative of increasing differentiation with time and elevation. Improved precision on the eruption ages indicates a recurrence interval of approximately 20 thousand years. The eruptive flux for both periods of rhyolitic volcanism is about 0.01 cubic kilometers per thousand years, which is less than the magma resurgence flux rates for syn-caldera and post-caldera eruptions of the Valles Caldera and Yellowstone Caldera volcanic systems. Collectively, these geochemical, geochronological, and volumetric data may facilitate a better understanding of heat flux and the longevity of magmatic sources related to geothermal resources in similar small-volume, silicic systems.

Utah↗
USGS1997

Integrity of production wells and confining unit at the Naval Weapons Industrial Reserve Plant, Dallas, Texas, 1995

Ground water in the shallow alluvial aquifer is contaminated at the Naval Weapons Industrial Reserve Plant, Dallas, Texas. Five production wells at the site are cased through the alluvial aquifer and underlying units and are screened in either the Paluxy or Twin Mountains aquifers. Three abandoned wells, originally completed in the Twin Mountains aquifer but filled with drilling mud in 1958, also penetrate the alluvial aquifer. The Paluxy and Twin Mountains aquifers are used for drinking-water supplies in and around the Dallas-Fort Worth metroplex. Trichloroethylene and its degradation products, dichloroethylene and vinyl chloride, and the metal chromium previously have been detected in the shallow alluvial aquifer. Current (1995) analyses of water-quality samples taken from the static water column of the five production wells and one of the abandoned wells indicate no trichloroethylene, dichloroethylene, or vinyl chloride in the water column of these wells. Chromium was detected in all samples, but concentrations were less than the practical quantitation limit, which is the regulatory action level for this site. The results of borehole geophysical log analysis indicate that two of the production wells could have vertically connected intervals where cement bonding in the well annulus is poor. The other production wells have overall good bonding. Temperature logs do not indicate flow behind casing except in the screened interval of one well. Geophysical logs show the Eagle Ford Shale ranges from 147 to 185 feet thick at the site. The Eagle Ford Shale has low permeability and a high plasticity index. These physical characteristics make the Eagle Ford Shale an excellent confining unit.

Texas↗
USGS2026

Arsenic and isotope concentrations in the lower Platte River valley of eastern Nebraska, early 1970s to 2023

The City of Lincoln, Nebraska, has been monitoring concentrations of arsenic in their source water and evaluating their options for treatment and removal since at least 2002. In 2022, the City of Lincoln, Nebr., with funding assistance from the Nebraska Water Sustainability Fund, began cooperating with the U.S. Geological Survey to examine arsenic concentrations in surface water and groundwater in the lower Platte River valley and the area around City of Lincoln Water System (LWS) well field. Arsenic data collected from the Platte River since 1974 were examined using the “weighted regression on time, discharge, and season” model, which compared the streamflow (also referred to as “discharge”), time of year, and season to estimate concentrations of arsenic. Annual mean arsenic concentrations modeled for more than 49 years at the Platte River at Louisville, Nebr., U.S. Geological Survey streamgage (station 06805500), indicated a significant increasing trend. Arsenic concentrations in the Platte River were seasonal, with the highest concentrations being observed during mid- to late summer. When seasonal patterns and streamflow were combined with arsenic concentrations in the Platte River during low streamflow conditions, groundwater contributions, which can have higher arsenic concentrations, make up a larger portion of the streamflow. Arsenic samples were collected from upstream rivers in 2022 and 2023 and were paired to analyze the arsenic contributions at the U.S. Geological Survey streamgage on the Platte River near Ashland, Nebr. (station 06801000), near the City of Lincoln well field. The arsenic concentrations from the streamgage on the Platte River near Ashland, Nebr., location, were higher than the U.S. Geological Survey streamgage on the Elkhorn River at Waterloo, Nebr. (station 06800500), and significantly lower than at the U.S. Geological Survey streamgage on the Platte River near Leshara, Nebr.(station 06796500), indicating that the Platte River usually contributes a higher concentration of arsenic than does the Elkhorn River as they join near Ashland, Nebr. During 1991–2023, six groundwater monitoring wells were analyzed to identify trends in arsenic concentrations. Two of the six wells had a positive trend during the 33-year period. One monitoring well did not reveal a long-term trend during this period but showed a trend during 2019–23, correlating to a period when the island in the middle of the Platte River was connected to the east bank of the river when manganese reducing conditions were present and groundwater levels were declining in the well. Across all wells the oxidation and reduction (redox) condition during the time of sampling was assessed. Mixed anoxic and (or) oxic redox condition was the most common redox process and the highest sampled arsenic concentrations in monitoring wells were observed in anoxic conditions driven by manganese reduction. Groundwater arsenic concentrations had seasonal variation around the City of Lincoln well field, with higher arsenic concentrations tending to be further south in comparison to samples collected further north. Isotope samples were collected and analyzed in surface water and groundwater around the LWS well field. The samples indicate that the proportion of surface water present in the LWS production wells can be higher in the spring and lower in the summer. With higher arsenic concentrations observed in the stream water during the summer period, the LWS source water can be affected by these elevated arsenic concentrations even though the proportion of surface water is lower.

Nebraska↗