Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Frontiers for Young Minds”

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,406 records · Page 16Linked to original sources
USGS2025

Airborne geophysical analysis to decipher salinization for coastal Louisiana

Coastal Louisiana is known for saltwater intrusion that threatens wetlands, aquifers, and rivers. However, the extent of saltwater intrusion is not well understood. This study develops an innovative framework with airborne electromagnetic (AEM) data to map chloride concentration distributions for wetlands in the Mississippi River deltaic plain and Chenier plain as well as for the Mississippi River Valley alluvial aquifer (MRVA) and Chicot aquifer. Moreover, the framework maps chloride concentrations along the Mississippi River and Atchafalaya River. Key components in the framework include the establishment of resistivity-to-chloride concentration transformation, 3D resistivity architecture building through geostatistics, and the employment of a lithologic model. The transformation functions correlate AEM resistivity data with porewater salinity measurements and groundwater and river chloride samples. The results show that AEM data reliably infers soil water chloride concentrations and correlates well with the distribution of various marsh types. AEM data reveals extensive saltwater presence at depth and near the coast, originating from salt domes and the Gulf of Mexico, respectively. The saltwater upconing pattern in the Chicot aquifer is likely due to excessive groundwater withdrawals. The AEM data also confirms a distinct tongue of saltwater intruding into the Atchafalaya Basin from the Gulf. The AEM data helps to identify faults that are obscured or eroded at the surface, which appear as leaky barriers in the subsurface where dramatic changes in chloride concentration are apparent. Finally, this study uses the AEM data to infer the presence of an extensive seawater wedge in the Mississippi River and Atchafalaya River.

Louisiana↗
USGS2024

Comparing subduction ground-motion models to observations for Cascadia

We evaluate Cascadia subduction ground-motion models (GMMs), considered for the 2023 US National Seismic Hazard Model (NSHM) update, by comparing observations to model predictions. The observations comprise regional recordings from intraslab earthquakes, including contributions from 2021 and 2022 events in southern Cascadia and global records from interface earthquakes. Since the 2018 NSHM update, new GMMs for Cascadia have been published by the Next Generation Attenuation (NGA)-Subduction Project that require independent evaluation. In the regional intraslab comparisons, we highlight a characteristic frequency dependence for Cascadia data, with short periods having lower ground motions and longer periods being comparable to other subduction zones. We evaluate differences in northern and southern Cascadia and find that the NGA-Subduction GMMs developed using southern Cascadia data perform better in this region than the model that did not consider these data. We compare ground-motion variability in Cascadia with the NGA-Subduction model predictions and find differences at short periods ( T = 0.1 s) due to the use of global versus regional data in the development of these models. Moreover, the within-event component of aleatory variability from the GMMs overpredicts the standard deviation of Cascadia recordings at very short periods ( T < 0.05 s). Using global interface earthquakes as a proxy to evaluate the Cascadia GMMs, we find long-period overprediction from a simulation-based GMM and some of the empirical GMMs. When comparing recent observations, we find a similar misfit to GMMs and the 2010 and 2022 Ferndale earthquakes. Finally, we observe different basin amplification factors arising in different subsets of the data, which indicate that differences in basin factors between empirical GMMs could arise from the data selection choices by the developers. As part of evaluating the regional basin terms, we apply basin amplification factors from the magnitude 9 Cascadia earthquake simulations to the empirical GMMs for interface earthquakes. The comparisons presented in this study indicate that the NGA-Subduction GMMs for Cascadia perform well relative to observations and older subduction GMMs.

British Columbia, California, Oregon, Washington↗
USGS2024

Constraining mean landslide occurrence rates for non-temporal landslide inventories using high-resolution elevation data

Constraining landslide occurrence rates can help to generate landslide hazard models that predict the spatial and temporal occurrence of landslides. However, most landslide inventories do not include any temporal data due to the difficulties of dating landslide deposits. Here we introduce a method for estimating the mean landslide occurrence rate of deep-seated rotational and translational slides derived solely from high-resolution (≤3 m) elevation data and globally available estimates of the diffusion coefficient for sediment flux. The method applies a linear diffusion model to the roughest landslide deposits until they reach a representative non-landslide roughness distribution. This estimates the time for a landslide deposit to be unrecognizable in high-resolution digital elevation data, which we term the mean lifetime of the landslide. Using the mean lifetime and number of landslides within an area of interest, we can estimate the mean occurrence rate of landslides over that domain. We validate this approach using a comprehensive temporal inventory of landslides in western Oregon created using age-roughness curves that are calibrated with high-resolution elevation data and radiocarbon data. We find good agreement between our diffusion method and the existing age-roughness-derived estimates, producing mean lifetimes of 4500 and 5200 years (4% difference), respectively. Hazard maps produced using the two methodologies generally agree, with the maximum differences in landslide probability reaching 0.1. Due to the relative abundance of high-resolution elevation data compared with age-dated landslides, our method could help constrain landslide occurrence rates in areas previously considered unfeasible.

Oregon↗
USGS2000

Water quality of selected springs and public-supply wells, Pine Ridge Indian Reservation, South Dakota, 1992-97

This report presents results of a water-quality study for the Pine Ridge Indian Reservation, South Dakota. The study was a cooperative effort between the U.S. Geological Survey and the Water Resources Department of the Oglala Sioux Tribe. Discharge and water-quality data were collected during 1992-97 for 14 contact springs located in the northwestern part of the Reservation. Data were collected to evaluate potential alternative sources of water supply for the village of Red Shirt, which currently obtains water of marginal quality from a well completed in the Inyan Kara aquifer. During 1995-97, water-quality data also were collected for 44 public-supply wells that serve about one-half of the Reservation's population. Quality-assurance sampling was used to evaluate the precision and accuracy of environmental samples. Ten of the springs sampled contact the White River Group, and four contact the Pierre Shale. Springs contacting the White River Group range from calcium bicarbonate to sodium bicarbonate water types. Two springs contacting the Pierre Shale have water types similar to this; however, sulfate is the dominant anion for the other two springs. In general, springs contacting the White River Group are shown to have better potential as alternative sources of water supply for the village of Red Shirt than springs contacting the Pierre Shale. Nine of the springs with better water quality were sampled repeatedly; however, only minor variability in water quality was identified. Six of these nine springs, of which five contact the White River Group, probably have the best potential for use as water supplies. Discharge from any of these six springs probably would provide adequate water supply for Red Shirt during most periods, based on a limited number of discharge measurements collected. Concentrations of lead exceeded the U.S. Environmental Protection Agency (USEPA) action level of 15 µg/L for three of these six springs. Five of these six springs also had arsenic concentrations that exceeded 10 µg/L, which could be problematic if the current maximum contaminant level (MCL) is lowered. Blending of water from one or more springs with water from the existing Inyan Kara well may be an option to address concerns regarding both quantity and quality of existing and potential sources. All nine springs that were sampled for indicator bacteria had positive detections on one or more occasions during presumptive tests. Although USEPA standards for bacteria apply only to public-water supplies, local residents using spring water for domestic purposes need to be aware of the potential health risks associated with consuming untreated water. One spring contacting the White River Group and two springs contacting the Pierre Shale exceeded 15 pCi/L for gross alpha; these values do not necessarily constitute exceedances of the MCL, which excludes radioactivity contributed by uranium and radon. Additional sampling using different analysis techniques would be needed to conclusively determine if any samples exceeded this MCL. Nine springs were sampled for selected pesticides and tritium. The pesticides atrazine, carbaryl, and 2,4-D were not detected in any of the samples. The nine springs were analyzed for tritium in order to generally assess the age of the water and to determine if concentrations exceeded the MCL established for gross beta-particle activity. Tritium results indicated two springs are composed primarily of water recharged prior to atmospheric testing of nuclear bombs and two other springs have a relatively large percentage of test-era water. The remaining five springs had tritium values that indicated some percentage of test-era water; however, additional sampling would be needed to determine whether water is predominantly pre- or post-bomb age. Of the 44 public-supply wells sampled, 42 are completed in the Arikaree aquifer, one is completed in an alluvial aquifer, and one is completed in the Inyan Kara aquifer. Water from the alluvial well is a sodium bicarbonate water type, water from Arikaree aquifer ranges from calcium bicarbonate to sodium bicarbonate types, and water from the Inyan Kara well is a calcium sulfate bicarbonate type. Of the 44 wells sampled, 28 (64 percent) tested positive for indicator bacteria in presumptive tests. Because these were single samples that generally were collected upstream from chemical treatment feeders, positive detections do not necessarily constitute exceedances of drinking-water standards. A single sample from an Arikaree well exceeded the MCL for arsenic of 50 µg/L. Arsenic exceeded 10 µg/L for six additional Arikaree wells and for the alluvial well and the Inyan Kara well, which could be problematic if the current MCL is lowered. The alluvial well also exceeded the secondary maximum contaminant level (SMCL) for dissolved solids, which is non-enforceable, and the action level for lead. The Inyan Kara well exceeded the SMCL's for iron and for manganese and the MCL of 5 pCi/L for radium-226 and 228 combined. Several Arikaree wells exceeded SMCL's for either pH, sulfate, dissolved solids, iron, or manganese. One Arikaree well exceeded the MCL of 4.0 mg/L for fluoride and another exceeded the MCL of 10 mg/L for nitrite plus nitrate. Ten Arikaree wells equalled or exceeded 15 pCi/L for gross alpha; however, these values do not necessarily constitute exceedances of the MCL, which excludes radioactivity contributed by uranium and radon. Additional sampling using different analysis techniques would be needed to conclusively determine if any samples exceeded this MCL. Eight wells, all from the Arikaree aquifer, equalled or exceeded the proposed MCL of 20 µg/L for uranium and 33 wells (75 percent) equalled or exceeded one-half of the proposed MCL. Although this standard has only been proposed, additional information regarding the extent of elevated uranium concentrations in the Arikaree aquifer, and the geochemical processes involved, may be beneficial. It was determined from analyses of uranium isotope data for five wells that the source of elevated uranium concentrations is naturally occurring, rather than anthropogenic.

South Dakota↗
USGS2025

A scientific vision and roadmap for earthquake rupture forecast developments, a USGS perspective

We articulate a scientific vision and roadmap for the development of improved Earthquake Rupture Forecast models, which are one of the two main modeling components used in modern seismic hazard and risk analysis. One primary future objective is to provide fully time-dependent models that include both elastic rebound and spatiotemporal clustering nationwide, which is particularly important for shorter-term hazard and risk considerations (e.g., earthquake insurance products). We also discuss the importance and perennial challenges associated with quantifying epistemic uncertainties, including those associated with deformation-model slip rates, un-quantified sampling errors with respect to off-fault seismicity, and any spatial covariances. The need for more physics-based approaches is also emphasized, as is the benefit of adding model valuation (quantifying usefulness) to our verification and validation protocols. Given the multidisciplinary and system-level nature of this activity, modular design is critical. Future updates will also draw from best-available science by both the United States Geological Survey and the external community. The primary goal of this paper is to highlight plans that guide research and facilitate community engagement with model development, especially with respect to lowering the entry barrier for early career scientists and engineers. The paper is written so readers can focus on the sections that interest them most (see table of contents), with the Introduction and Discussion providing a stand-alone overview and summary.

Bulletin of the Seismological Society of America↗
USGS1995

Evaluation of ground-water flow and hydrologic budget for Lake Five-O, a seepage lake in northwestern Florida

Temporal and spatial distributions of ground-water inflow to, and leakage from Lake Five-O, a softwater, seepage lake in northwestern Florida, were evaluated using hydrologic data and simulation models of the shallow ground-water system adjacent to the lake. The simulation models indicate that ground-water inflow to the lake and leakage from the lake to the ground-water system are the dominant components in the total inflow (precipitation plus ground-water inflow) and total outflow (evaporation plus leakage) budgets of Lake Five-O. Simlulated ground-water inflow and leakage were approximately 4 and 5 times larger than precipitation inputs and evaporative losses, respectively, during calendar years 1989-90. Exchanges of water between Lake Five-O and the ground-water system were consistently larger than atmospheric-lake exchanges. A consistent pattern of shallow ground-water inflow and deep leakage was also evident throughout the study period. The mean time of travel from ground-water that discharges at Lake Five-O (time from recharge at the water table to discharge at the lake) was estimated to be within a range of 3 to 6 years. Flow-path evaluations indicated that the intermediate confining unit probably has a negligible influence on the geochemistry of ground-water inflow to Lake Five-O. The hydrologic budgets and flow-path evaluations provide critical information for developing geochemical budgets for Lake Five-O and for improving the understanding of the relative importance of various processes that regulate the acid-neutralizing capacity of softwater seepage lakes in Florida.

Florida↗
USGS2025

No evidence for an active margin-spanning megasplay fault at the Cascadia Subduction Zone

It has been previously proposed that a megasplay fault within the Cascadia accretionary wedge, spanning from offshore Vancouver Island to Oregon, has the potential to slip during a future Cascadia subduction zone earthquake. This hypothetical fault has major implications for tsunami size and arrival times and is included in disaster-planning scenarios currently in use in the region. This hypothesis is evaluated in this study using CASIE21 deep-penetrating and U.S. Geological Survey high-resolution seismic reflection profiles. We map changes in wedge structural style and seismic character to identify the inner-outer wedge transition zone where a megasplay fault has been previously hypothesized to exist and evaluate evidence for active faulting within this zone. Our results indicate that there is not an active, through-going megasplay fault in Cascadia, but instead, the structure and activity of faulting at the inner-outer wedge transition zone is highly variable and segmented along strike, consistent with the segmentation of other physical and mechanical properties in Cascadia. Wedge sedimentation, plate dip, and subducting topography are proposed to play a major role in controlling megasplay fault development and evolution. Incorporating updated megasplay fault location, geometry, and activity into modeling of Cascadia earthquakes and tsunamis could help better constrain associated hazards.

British Columbia, Oregon, Washington↗
USGS2025

lasertram: A Python library for time resolved analysis of laser ablation inductively coupled plasma mass spectrometry data

Laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) data has a wide variety of uses in the geosciences for in-situ chemical analysis of complex natural materials. Improvements to instrument capabilities and operating software have drastically reduced the time required to generate large volumes of data relative to previous methodologies. Raw data from LA-ICP-MS, however, is in counts per unit time (typically counts per second), not elemental concentrations and converting these count ratesto concentrations requires additional processing. For complex materials where the ablated volume may contain a range of material compositions, a moderate amount of user input is also required if appropriate concentrations are to be accurately calculated. In geologic materials such as glasses and minerals that potentially have numerous heterogeneities (e.g., microlites or other inclusions) within them, this is typically determiningwhether the total ablation signal should be filtered to remove these heterogeneities. This necessitates that the LA-ICP-MS data processing pipeline is one that is not automated, but is also designed to enable rapid and efficient processing of large volumes of data. Here we introduce , a Python library for the time resolved analysis of LA-ICP-MS data. We outline its mathematical theory, code structure, and provide an example of how it can be used to provide the time resolved analysis necessitated by LA-ICP-MS data of complex geologic materials. Throughout the pipeline we show how metadata and data are incrementally added to the objects created such that virtually any aspect of an experiment may be interrogated and its quality assessed. We also show, that when combined with other Python libraries for building graphical user interfaces, it can be utilized outside of a pure scripting environment. can be found at https://doi.org/10.5066/P1DZUR3Z

Applied Computing and Geosciences↗
USGS2025

Groundwater age estimates for the Mississippi River Valley alluvial aquifer based on tracer data collected during 2018–20

This study characterized groundwater age across the Mississippi River Valley alluvial aquifer (MRVA). Groundwater samples from 69 MRVA wells and 19 wells in Tertiary units of the Mississippi embayment aquifer system (MEAS) were analyzed for sulfur hexafluoride (SF 6 ), tritium ( 3 H), helium (He), and (or) carbon-14 of dissolved inorganic carbon ( 14 C). The age distributions of 89 samples were estimated by fitting lumped parameter models to processed tracer concentrations with the U.S. Geological Survey software TracerLPM. Mean ages of MRVA groundwater samples ranged from 12 to 22,000 years, with a median of 140 years. Mean ages of MEAS groundwater samples ranged from 230 to 52,000 or more years, with a median of 13,500 years. The spatial distribution of MRVA groundwater ages was found to be influenced by depth, inflow of groundwater from deeper units, and soil saturated hydraulic conductivity. In parts of the MRVA, the spatial distribution of MRVA groundwater ages was found to be influenced by annual recharge and (or) annual groundwater pumpage.

Alabama, Arkansas, Illinois, Kentucky, Louisiana, ↗
USGS2024

Forecasting storm-induced coastal flooding for 21st century sea-level rise scenarios in the Hawaiian, Mariana, and American Samoan Islands

Oceanographic, coastal engineering, ecologic, and geospatial data and tools were combined to evaluate the increased risks of storm-induced coastal flooding in the populated Hawaiian, Mariana, and American Samoan Islands as a result of climate change and sea-level rise. We followed a hybrid (dynamical and statistical) downscaling approach to map flooding due to waves and storm surge at 10-square meter resolution along all 1,870 kilometers of these islands’ coastlines for annual (1-year), 20-year, and 100-year return-interval storm events and +0.00 meter (m), +0.25 m, +0.50 m, +1.00 m, +1.50 m, +2.00 m, and +3.00 m sea-level rise scenarios. We quantified the coastal flood depths and extents using the latest climate forcing from Intergovernmental Panel for Climate Change’s Sixth Assessment Report Coupled Model Intercomparison Project. The data generated using these methods provide stakeholders and decision makers with a spatially explicit, rigorous valuation of how, where, and when climate change and sea-level rise increase coastal storm-induced flooding to help identify areas where management and (or) restoration could potentially help reduce the risk to, and increase the resiliency of, the coastal communities in the populated Hawaiian, Mariana, and American Samoan Islands.

Hawaii↗
USGS2024

Improving crop-specific groundwater use estimation in the Mississippi Alluvial Plain: Implications for integrated remote sensing and machine learning approaches in data-scarce regions

Study region The Mississippi Alluvial Plain (MAP) in the United States (US). Study focus Understanding local-scale groundwater use, a critical component of the water budget, is necessary for implementing sustainable water management practices. The MAP is one of the most productive agricultural regions in the US and extracts more than 11 km 3 /year for irrigation activities. Consequently, groundwater-level declines in the MAP region pose a substantial challenge to water sustainability, and hence, we need reliable groundwater pumping monitoring solutions to manage this resource appropriately. New hydrological insights for the region We incorporate remote sensing datasets and machine learning to improve an existing lookup table-based model of groundwater use previously developed by the U.S. Geological Survey (USGS). Here, we employ Distributed Random Forests, an ensemble machine learning algorithm to predict annual and monthly groundwater use (2014–2020) throughout this region at 1-km resolution, using pumping data from existing flowmeters in the Mississippi Delta. Our model compares favorably with the existing USGS model, with higher R 2 (0.51 compared to 0.42 in the previous model), and lower root mean square error (RMSE) and mean absolute error (MAE)— 0.14 m and 0.09 m, respectively in our model, compared to 0.15 m and 0.1 m in the previous model. Therefore, this work advances our ability to predict groundwater use in regions with scarce or limited in-situ groundwater withdrawal data availability.

Journal of Hydrology Regional Studies↗
USGS2025

Applying U-Th disequilbrium for dating siliceous sinters

Continental hydrothermal systems are critical avenues for the crustal transport of heat and mass captured for geothermal energy and mineral exploration. Thus, understanding their temporal evolution and longevity is important for resource characterization. Deposits of microlaminated siliceous sinter, common surface expressions of high temperature reservoirs (> 170 °C), have the potential to trace hydrothermal histories. Geothermal reservoirs are often located within uranium-bearing silicic volcanic rock where subsurface fluid-rock interactions extract U into hydrothermal fluids. U incorporated in the surface sinter deposit has the potential to provide a U— Th disequilibrium dating option. We focus on samples from El Tatio geyser field in the Altiplano of northern Chile, the largest geothermal system in the Andes. Our resulting 230 Th/U ages, along with the water and deposit elemental compositions, suggest concentrations of U and Th vary predictably along the sinter apron. While distal facies containing the highest U concentrations (> 50 μg/g) are least affected by detrital Th corrections, they can display suspected open-system behavior. In contrast, more medial facies, where bacterial mats and other porous textures are co mmonly concentrated, have only trace amounts of U (< 0.1 μg/g), which leads to unreliable or geologically improbable dates. Proximal facies tend to date most consistently. By comparing existing 14 C ages with 230 Th/U results, 230 Th/U ages tend to be younger than the 14 C ages, supporting the presence of a 14 C-dead carbon influence. New data confirm that the onset of geothermal activity at El Tatio goes back to the late Pleistocene.

Altiplano, El Tatio geyser field↗
USGS2026

Site response models based on geometric parameters for southern California sedimentary basins

Site response in sedimentary basins is influenced by complex three-dimensional (3D) features, including trapping of seismic waves, focusing of seismic energy and basin resonance. Current ground motion models (GMMs) incorporate basin effects using one-dimensional parameters like V S30 and shear wave velocity isosurface depths, which are limited in capturing lateral and 3D effects. To address these limitations, we develop seismic site response models based on novel parameters that represent multi-dimensional properties of the Los Angeles Basin (LAB) geometry and shear wave velocity. We define a basin shape for the LAB using depth to subsurface geologic interfaces associated with the oldest sedimentary deposits (depth to a particular shear wave velocity horizon, i.e., 1.5 km/s - z 1.5 ) and the depth to the crystalline basement ( z cb ) which are determined using geologic cross sections and community seismic velocity model profiles. We explore a suite of geometric descriptors computed for the LAB and southern California, from which three parameters with the greatest predictive potential are selected and evaluated using empirical ground motion residual analyses in combination with the Boore et al. GMM. The results demonstrate that the zonal heterogeneity index ( ), standard deviation of the absolute difference between z 1.5 and z cb ( ) and standard deviation of z cb ( ) each provide a reduction in site-to-site variability ( ϕ S2S ) of empirical GMMs. The reduction in ϕ S2S is period-dependent, with average decreases of 3%, 26% and 6% for , , and , respectively. Although these reductions are modest from an engineering application perspective, they are statistically significant, underscoring the inherent difficulty in fully characterising complex basin effects. Collectively, these findings indicate that the inclusion of basin-specific geometric parameters yields measurable, albeit incremental, improvements in site response prediction and establishes a framework for the progressive refinement of seismic hazard characterisation within sedimentary basins.

California↗
USGS2026

Evaluation of submersible pressure transducers for streamflow monitoring in small streams

Compact streamgages requiring minimal infrastructure and equipped with submersible pressure transducers (PTs) are increasingly used to monitor small streams, yet disparate implementations obscure their accuracy under real-world conditions. This study isolated instrumentation-derived uncertainty in stage monitoring by co-locating various combinations of commercial vented and unvented PTs with seven U.S. Geological Survey (USGS) reference streamgages on small streams. Multi-year PT stage records from compact streamgages, collected and corrected following USGS protocols, were compared to concurrent reference observations. Vented PTs demonstrated an average measurement uncertainty of ±0.005 m and mean absolute percent error (MAPE) of ±0.2%. Unvented PTs exhibited higher uncertainty, averaging ±0.009 m and ±0.4% MAPE. Although both sensor types had stage errors up to ±0.3 m, 95% of vented PT and unvented PT observations were within 0.01 and 0.02 m of reference stage, respectively. Analysis of additional unvented PTs revealed stage errors of up to ±1.2% when using barometric sensors within 15 km of the in-water sensor. Propagation of stage error to discharge using reference rating models resulted in cumulative discharge MAPEs of ±4.5% for vented and ±5.5% for unvented PTs. These findings highlight PTs as practical alternatives to reference instrumentation when deployed with standardized procedures, potentially expanding access to reliable streamflow data.

conterminous United States↗
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↗
USGS2025

Methodology and technical input for the 2025 U.S. List of Critical Minerals—Assessing the potential effects of mineral commodity supply chain disruptions on the U.S. economy

The Secretary of the Interior, acting through the Director of the U.S. Geological Survey, is tasked by section 7002 (“Mineral Security”) of title VII (“Critical Minerals”) of the Energy Act of 2020 (Public Law 116–260, December 27, 2020, 116th Congress) with reviewing and revising the methodology used to evaluate mineral commodity supply risk and the U.S. List of Critical Minerals (LCM) no less than every 3 years. Following two previous LCM assessments, this analysis represents the latest technical input for evaluating each mineral commodity’s supply risk and determining their recommended status on the LCM. We evaluated mineral commodity supply risk using two criteria: (1) an economic effects assessment that quantified the potential effects of various trade disruption scenarios on the U.S. economy, and (2) an examination of whether the mineral commodity’s U.S. supply chain relied on a sole domestic producer that represented a single point of failure. For the first criterion, postdisruption equilibrium quantities and prices for each mineral commodity were calculated based on their price elasticities of supply and demand and the availability of excess production capacity for each yearlong foreign trade disruption scenario. Subsequently, a nonlinear optimization routine was used with detailed economic input-output tables to estimate the potential economic effects on the U.S. economy of over 1,200 scenarios for 84 mineral commodities. After accounting for the probability of each scenario’s occurrence, the overall results are presented in terms of changes in U.S. gross domestic product (GDP) by individual industry and the economy overall. The results, which ranged from a net decrease in U.S. GDP of nearly $4.5 billion to a net increase of $33 million, largely reflect U.S. import dependency and world production concentration. Using the Jenks natural breaks optimization method, a statistical classification technique, we categorized the mineral commodities into several classes based on this overall risk quantification. Mineral commodities with annualized probability-weighted net decreases in U.S. GDP greater than $2 million were recommended for inclusion on the LCM. If a mineral commodity did not meet the threshold for inclusion on the LCM under the first criterion, its domestic supply chain was examined under the second criterion, which recommended a mineral commodity for inclusion on the LCM if there was only a single domestic producer. Ultimately, the two criteria resulted in the recommendation of the addition of six mineral commodities (in descending risk order, potash, silicon, copper, silver, rhenium, and lead) to and the removal of two mineral commodities (arsenic and tellurium) from the LCM. By using an economic effects assessment, the results of this analysis provide a prioritization that can also be compared directly against other risk analyses and the cost of various risk mitigation strategies.

Open-File Report↗
USGS2012

Development of vegetation based soil quality indices for mineralized terrane in arid and semi-arid regions

Soil quality indices (SQIs) are often management driven and attempt to describe key relationships between above- and below-ground parameters. In terrestrial systems, indices that were initially developed and modified for agroecosystems have been applied to non-agricultural systems in increasing number. We develop an SQI in arid and semi-arid ecosystems of the Western US impacted by different types of geologic mineralization using the relationship between vegetation community parameters and soil abiotic and biotic properties. We analyze these relations in soils associated with three different mineralization types: podiform chromite, Cu/Mo porphyry, and acid-sulfate gold vein systems at four different sites in California and Nevada. Soil samples were collected from undisturbed soils in both mineralized and nearby unmineralized substrates as well as from waste rock and tailings. Aboveground net primary productivity (ANPP), canopy cover and shrub density were measured for the vegetative communities. Minimum data sets were developed based on correlations between the soil and vegetation parameters, refined using principal components analysis, scored using non-linear functions, and combined into an overall SQI. The indices are comprised of one or two microbial parameters and three to six abiotic parameters, the latter consisting of nutrients and metals. Given the preliminary development of this approach, the parameters and combinations to arrive at an SQI for a given site cannot at this time be correlated or compared with that of another site. This SQI approach provides a means of quantifying disturbed ecosystem recovery resulting from mining, and could be applied to other disturbances in a way that readily distills the information for potential use by land managers. However, severely disturbed areas with little to no aboveground biomass, such as unreclaimed tailings, have likely crossed an ecological threshold that precludes the use of this type of monitoring tool.

California↗
USGS2024

Groundwater flow model for the Des Moines River alluvial aquifer near Des Moines, Iowa

Des Moines Water Works (DMWW) is a regional municipal water utility that provides residential and commercial water resources to about 600,000 customers in Des Moines, Iowa, and surrounding municipalities in central Iowa. DMWW has identified a need for increased water supply and is exploring the potential for expanding groundwater production capabilities in the Des Moines River alluvial aquifer, where it operates two radial collector wells (RCWs). The U.S. Geological Survey, in cooperation with DMWW, completed a study of the Des Moines River alluvial aquifer and interactions of the RCWs with the aquifer; no previously published model has included the existing well locations, which is the focus of this model. A conceptual and numerical groundwater flow model have been developed to characterize the Des Moines River alluvial aquifer under existing conditions, to simulate water levels observed in the RCWs, and to provide publicly accessible hydrologic data and research that advance understanding of the regional hydrologic system and can potentially be used in the future to evaluate groundwater production scenarios. Model performance was assessed by comparing observed and simulated groundwater levels that included water level elevations, water level changes, water level inequality observations, surface water streamflow, and change in surface water volume from upstream to downstream. Water table elevation in the aquifer layers is on average slightly overestimated with average absolute value error less than 1.5 meters at both RCWs and less than 2.5 meters for all observation wells in the alluvial aquifer layers. The model also accurately simulated water tables greater than the RCW design minimum (a water level threshold at which RCW pumping is reduced) in all timesteps for which water level observation data existed. Water table elevation error was higher in other model layers that were not the focus of the study, and the model did not accurately match streamflow targets.

Iowa↗