SEARCH · USGS Science

Results for “Journal of Analytical and Applied Pyrolysis”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,688 records · Page 19Linked to original sources

USGS2024

Conceptualization and simulation of groundwater flow and groundwater availability in the Boone and Roubidoux aquifers in northeastern Oklahoma, 1980–2017

Oklahoma Groundwater Law (Oklahoma Statute § 82-1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations to determine the maximum annual yield for the State’s groundwater basins. The Boone and Roubidoux aquifers (also known as the Springfield Plateau aquifer and Ozark aquifer, respectively) are bedrock aquifers that extend from northeastern Oklahoma into Kansas, Arkansas, and Missouri. At present (2024), the Oklahoma Water Resources Board has yet to legally issue orders for the final determination of maximum annual yields for the Boone and Roubidoux aquifers. To support determination of a maximum annual yield, the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, developed a hydrogeologic framework, a conceptual groundwater-flow model, and a calibrated numerical groundwater-flow model for the Boone and Roubidoux aquifers. Three types of groundwater-availability scenarios were simulated by using the calibrated numerical model. These scenarios were used to (1) estimate equal-proportionate-share groundwater withdrawal rates (groundwater withdrawal applied equally over the aquifer), (2) quantify the potential effects of projected groundwater withdrawals on groundwater storage over a 50-year period, and (3) simulate the potential effects of a hypothetical 10-year drought. For the Boone aquifer, equal-proportionate-share groundwater withdrawal rates were 1.10, 0.98, and 0.96 acre-feet per acre per year for the 20-, 40-, and 50-year scenarios, respectively. For the Roubidoux aquifer, equal-proportionate-share groundwater withdrawal rates were 1.76, 1.34, and 1.25 acre-feet per acre per year for the 20-, 40-, and 50-year simulations, respectively. For the 50-year scenarios, stream seepage was minimally affected. Over the 10-year drought scenario, groundwater storage in the Boone and Roubidoux aquifers decreased by 660,451 acre-feet (6.7 percent) and 508,472 acre-feet (1.0 percent), respectively.

Oklahoma↗

USGS2026

The 2025 Puerto Rico and Virgin Islands U.S. National Seismic Hazard Model Update: Ground motion model selection and comparison

We evaluate, select, and describe the ground-motion models (GMMs) used in the 2025 update of the U.S. National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). We identify the most appropriate models that align with GMM selection criteria for use in the PRVI region to improve the accuracy of seismic hazard assessments. The update incorporates globally applicable GMMs suited for the active crustal and subduction earthquakes in the region. We include region-specific adjustments to these GMMs derived from local site response analyses derived from ground motion records. The unadjusted and regionally-corrected GMMs are combined to create a robust model for predicting median ground motion. The model integrates epistemic uncertainty through a median ground motion logic tree that accounts for variations in magnitude and distance. This study compares the GMMs selected for the 2025 PRVI NSHM, including both as-provided and regionally adjusted NGA-West2 and NGA-Subduction models, with those used in the 2003 PRVI NSHM. We evaluate how changes in model selection, weighting, aleatory variability, and epistemic uncertainty influence seismic hazard estimates. Trends with distance, magnitude, and spectral period are analyzed to evaluate how the scaling behavior of the newer GMMs differs from that of earlier models. Relative to the GMMs used in the 2003 NSHM for this region, the 2025 models generally predict lower ground motions. Comparisons with additional GMMs indicate that the adjustments applied for PRVI are consistent with regional-specific modifications developed elsewhere globally. The increase in aleatory variability and epistemic uncertainty in the 2025 update results in a notable increase in hazard levels from these wider uncertainty bounds. These changes can result in as much as a 10%–20% variation in probabilistic ground motion at the 2% in 50 years exceedance level for hazard maps computed across the region for representative site classes and periods.

Puerto Rico, Virgin Islands↗

USGS2023

Don’t Let Negatives Hold You Back: Accounting for Underlying Physics and Natural Distributions of Hydrothermal Systems When Selecting Negative Training Sites Leads to Better Machine Learning Predictions

Selecting negative training sites is an important challenge to resolve when utilizing machine learning (ML) for predicting hydrothermal resource favorability because ideal models would discriminate between hydrothermal systems (positives) and all types of locations without hydrothermal systems (negatives). The Nevada Machine Learning project (NVML) fit an artificial neural network to identify areas favorable for hydrothermal systems by selecting 62 negative sites where the research team had confidence that no hydrothermal resource exists. Herein, we compare the implications of the expert selection of negatives (i.e., the NVML strategy) with a random sample strategy, where it is assumed that areas outside the favorable structural ellipses defined by NVML are negative. Because hydrothermal systems are sparse, it is highly probable that, in the absence of a favorable geological structure, hydrothermal favorability is low. We compare three training strategies: 1) the positive and negative labeled examples from NVML; 2) the positive examples from NVML with randomly selected negatives in equal frequency as NVML; and 3) the positive examples from NVML with randomly selected negatives reflecting the expected natural distribution of hydrothermal systems relative to the total area. We apply these training strategies to the NVML feature data (input data) using two ML algorithms (XGBoost and logistic regression) to create six favorability maps for hydrothermal resources. When accounting for the expected natural distribution of hydrothermal systems, we find that XGBoost performs better than the NVML neural network and its negatives. Model validation was less reliable using F1 scores, a common performance metric, than comparing probability estimates at known positives, likely because of the extreme natural class imbalance and the lack of negatively labeled sites. This work demonstrates that expert selection of negatives for training in NVML likely imparted modeling bias. Accounting for the sparsity of hydrothermal systems and all the types of locations without hydrothermal systems allows us to create better models for predicting hydrothermal resource favorability.

Geothermal Resources Council Transactions↗

USGS2025

Capturing the uncertainty of seismicity observations in earthquake rate estimates: Implications for probabilistic seismic hazard analysis and the USGS National Seismic Hazard Model

The rate of earthquakes in a region is a fundamental input to Probabilistic Seismic Hazard Analysis. We present a Monte Carlo method for computing that rate from seismicity catalogs while including a range of data and analysis uncertainties. This method is applied to regions for which the b value is assumed to be spatially invariant. Each region is broken down into epochs for which each epoch is estimated to have a uniform magnitude of completeness ( ⁠⁠ M c ). The distribution of earthquake rates for M ≥ M c is determined for each epoch by considering the Poisson likelihood of rates given the number of observed earthquakes with M ≥ M c ⁠ . We use a Monte Carlo process to include the uncertainty in b , ⁠ M c , and individual event magnitudes. The result for each epoch is the joint distribution of the Poisson rate of earthquakes with magnitudes larger than the minimum value used to calculate hazard ( ⁠⁠ M 1 ) and the Gutenberg–Richter b values, which control the extrapolation to other magnitudes. The rate for each region is either the duration‐weighted average over the epochs or, to better capture temporal variations, we also consider mixture models. The mixture models also provide an avenue to allow temporal variations in b values. To implement this joint distribution in a logic tree, we use the mean and 95% confidence branches, each of which is parameterized with an M ≥ M 1 rate and b value. We explore different ways of defining those branches, as well as non‐Gutenberg–Richter branches, and their impact on hazard estimates. The mean hazard, but not the fractiles, is robust with respect to these choices. To illustrate these new methods, we use synthetic data and catalogs from recent U.S. Geological Survey National Seismic Hazard Models for the Central and Eastern United States and for Puerto Rico and the U.S. Virgin Islands.

central and eastern United States, Puerto Rico, U.↗

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↗

USGS2026

Navigating uncertainty and competing objectives: Spring chinook salmon recovery in the upper Willamette river

Globally, anadromous fish populations are threatened with extinction due to multiple factors, including river impoundments that block migration, widespread alteration of physical habitat, temperature, and flow regimes, commercial and recreational harvest, changes in biological communities, and long-term climate trends. In this synthesis, we describe how physical and biological factors, policy and law, and public resource allocation affect management of spring-run Chinook Salmon ( Oncorhynchus tshawytscha ) in the Upper Willamette River (UWR), Oregon, USA. Efforts to restore salmon populations often involve trade-offs, for example, water management decisions must account for flood control, agricultural and municipal water use, and management of other species. UWR Chinook Salmon recovery efforts are also governed by multiple laws, policies, and action agencies that in some cases have incongruent objectives. Competing and conflicting objectives, combined with uncertainty regarding the outcomes of proposed management actions, create significant challenges for decision makers. Resource limitations, regulatory constraints, and incongruent objectives have created challenges for managers in the UWR, but opportunities exist for improved adaptive management and collaboration, which could build trust and engagement among stakeholders and ultimately support the successful implementation of decision tools. A suite of decision support tools that have been used in the Willamette River basin and elsewhere across the range of Pacific salmon have the potential to enhance decision processes. These tools can be effective when decision makers and stakeholders commit to long-term, collaborative utilization of the tools for decision making, and when available resources and regulatory environments are conducive to implementation and adaptive refinement of both the decision-support tools and management plans. The general biological and institutional principles discussed for UWR Chinook Salmon in this work apply more generally to populations of anadromous fishes across their range, especially to populations in highly regulated river systems.

Oregon↗

USGS2024

A methodology to estimate CO2 and energy gas storage resources in depleted conventional gas reservoirs

Depleted hydrocarbon reservoirs are subsurface geological structures capable of sequestering vast quantities of carbon dioxide (CO 2 ) as well as storing other energy gases for later usage, such as natural gas, and potentially hydrogen (H 2 ). Here we outline a methodology to quantify multi-gas storage resources in depleted conventional gas reservoirs for usage in assessments by the United States Geological Survey (USGS) at the scale of sedimentary basins. The methodology consists first of quantifying accessible pore volume in a depleted reservoir for natural gas storage using up to three equations. Input data are derived from commonly reported or estimated reservoir parameters and natural gas production volumes, and equations may be combined in linear models to improve pore volume estimates. Storage estimates from these equations are tested and validated for 31 reservoirs in the Michigan Basin Province, USA that were previously converted to underground gas storage facilities and have known (federally reported) natural gas storage capacities. Secondly, natural gas storage capacities can be transformed via fluid substitution calculations to estimate the storage resources for non-native fluids, applied here for, CO 2 , H 2 , and methane-H 2 blends, accounting for molecule-specific deviations from ideal gas behavior at reservoir pressures and temperatures as well as differing storage efficiencies. Importantly, the storage of non-native fluids may not be appropriate in all depleted gas reservoir settings due to potential risks like leakage, in particular in the case of H 2 storage, requiring additional knowledge of caprock sealing capacity. Given this caveat, we demonstrate the fluid substitution method for natural gas reservoirs of the Northern Niagaran Reef and Southern Niagaran Reef USGS plays in the Michigan Basin Province, as these trends of Silurian pinnacle reefs are capped with tight-sealing evaporite facies. The deterministic equations outlined from this methodology can be incorporated into future probabilistic USGS gas storage assessments for CO 2 , H 2 , and natural gas resources in the United States.

Michigan↗

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

Mapping potential sensitivity to hydrogeomorphic change in the UMRS riverscape

In 2020 the U.S. Geological Survey (USGS), as part of the Upper Mississippi River Restoration (UMRR) Program, began a new project to characterize potential hydrogeomorphic change associated with hydrogeomorphic units (HGUs) and their catenae (units linked by their association with sediment sources and flow origins). The goal of the project was to develop a geographic information system (GIS) database of HGUs for the Upper Mississippi River System (UMRS) available to both scientists and river managers working on UMRR studies and HREP planning and design studies. The characterization was based on a hydrogeomorphic change hierarchical classification developed previously for the UMRS. The products were generated with automated techniques in a GIS using systemic datasets. Landforms were mapped from the 2015 UMRS topobathymetric dataset with geomorphon (shape-based) tools tailored for the large riverscape, valley bottom environments in the UMRS. A clustering analysis was applied to the resulting landforms to identify HGUs associated with zones of perennial low flows, bankfull flows, and overbank floodplains. Catenae were assembled based on the proximity of the units to the main channel, tributary mouths, and side channels from previously published aquatic areas (USACE, 2018) coupled with least-cost flowpath linkages between potential sediment origins and planform change units developed by Rogala, Fitzpatrick, and Henderson (2020). These GIS-based analyses were successful at identifying a range of HGUs using an automated technique with available data across the entire riverscape, with emphasis on those that have the potential for hydrogeomorphic change. Most of the resulting features are depositional, as expected in a large river system. However, this is the first attempt of linking tributary inputs, side channel erosion and levee breaches with their depositional counterparts. The approach was successfully piloted in Pools 8 and 10 in the Upper Impounded Reach and Pool 14 in the Lower Impounded Reach, with next steps for application in reaches of the unimpounded section and Illinois River. This report emphasizes results from Pool 10, which was the focus of most of our attention during the pilot phase.

Illinois, Indiana, Iowa, Minnesota, Missouri, Sout↗

USGS2024

Separating signals in elevation data improves supervised machine learning predictions for hydrothermal favorability

A recent study identified topography (land surface elevation above sea level) as an important input dataset (feature) for predicting the location of hydrothermal systems in the Great Basin in Nevada. Yet, topography is generally a result of more than one geological process and may consequently contain multiple distinct signals. For example, the geologic evolution of the Great Basin has produced both crustal thickening (i.e., regional-scale trends in elevation) and thinning via Basin and Range extensional faulting (i.e., valley-scale topographic relief). We postulate that these geologic processes may affect the occurrence of hydrothermal systems differently. Therefore, we separate the regional trend from the valley-scale signal in the Great Basin, and then use them separately to evaluate the importance of each as predictors for hydrothermal favorability. Our prior work applying supervised machine learning (ML) using the data from the Nevada Machine Learning Project demonstrated that employing a training strategy that randomly selects negative training sites produces better performing models for predicting hydrothermal favorability than a training strategy that uses expert-selected negatives. The models created using both training strategies exhibited a west-east geographic trend in the predictions for the favorability of hydrothermal resources. These models generally predicted higher favorability in western Nevada and lower favorability in eastern Nevada. This west-east trend in predicted favorability correlates with elevation across the Great Basin, which trends higher from west to east. By separating the original elevation feature into distinct features for elevation trend (i.e., regional-scale topography) and detrended elevation (i.e., valley-scale or local relative topography), we find that models using the separated topographic signals consistently outperform competing models that use the original elevation feature. Although western Nevada still exhibits higher favorability than eastern Nevada, using separated signals for regional elevation and local structure reduces the west-east prediction trend in the region and emphasizes structures associated with hydrothermal upflow. This work emphasizes how carefully engineering features to represent geological conditions relevant to hydrothermal systems allows ML algorithms to detect important patterns for predicting hydrothermal resource favorability and leads to better model performance.

Conference Paper↗

USGS2024

Plan to coordinate post-earthquake investigations supported by the National Earthquake Hazards Reduction Program (NEHRP)

Introduction This report presents a plan supported by the National Earthquake Hazards Reduction Program (NEHRP) to coordinate domestic and international post-earthquake investigations (herein called “the Plan”). Post-earthquake scientific and engineering investigations are undertaken to capture critical information to understand the causes and impacts of the event, lessons from which can substantially improve the Nation’s resilience after future earthquakes. NEHRP is the Federal Government’s coordinated nationwide program to reduce risks to life and property from earthquakes. The Plan describes the activation and coordination of the four designated NEHRP Agencies in the Federal Government: Federal Emergency Management Agency, National Institute of Standards and Technology, U.S. National Science Foundation, and U.S. Geological Survey (USGS). The Plan also describes coordination between NEHRP Agencies and other organizations that may participate in pre-event and post-earthquake investigations, including non-NEHRP Federal agencies; State, regional, local, Tribal, and territorial agencies; domestic nongovernmental organizations; academic institutions and affiliated organizations; private companies; foreign governmental agencies and nongovernmental organizations; and international organizations. The Plan delineates the coordination of NEHRP post-earthquake scientific and engineering investigations to document the direct, indirect, and cascading physical and societal impacts from fault rupture and ground shaking hazards and from secondary hazards such as landslides, liquefaction, and tsunamis. In addition, the Plan identifies pre-event activities necessary to ensure that post-earthquake investigations are executed effectively. The USGS is the lead NEHRP Agency for activating and coordinating NEHRP post-earthquake investigations and for implementing this Plan. The USGS also leads coordination of the NEHRP Agencies in completing the pre-event activities identified in the Plan. The Plan has new information and supersedes USGS Circular 1242, “The Plan to Coordinate NEHRP Post-Earthquake Investigations,” which was published in 2003. This second-generation Plan was developed with the assistance of the Applied Technology Council (ATC) of Redwood City, California, under USGS contract 140G0121P0309, ATC-155. A 17-member Project Review Panel provided guidance on plan development, and input was solicited from subject matter experts representing key stakeholder groups and from participants at a public workshop.

Circular↗

USGS2024

Application of non-stationary shear-wave velocity randomization approach to predict 1D seismic site response and its variability at two downhole array recordings

Accounting for uncertainties in seismic site response is crucial to improving the performance of one-dimensional (1D) ground response analyses (GRAs) at downhole array recording sites. In addition to site effects, uncertainties in 1D-GRAs can also be contributed from the seismic source and/or path. Though often representing not more than one percent of the distance (path) from the source, site conditions are known to have an enormous influence on ground shaking. In this study, we focus on the site shear-wave velocity ( V S ) structure, which is the main ingredient for estimating the variability of site response. As such, V S can manifest aleatory uncertainties related to the effects of small-scale spatial heterogeneities within the near surface, thus V S can substantially modify ground shaking during earthquakes. We apply a novel V S randomization approach to propagate the small-scale heterogeneities of V S to estimate seismic site response within a non-stationary probabilistic framework. The randomization approach generates samples of V S profiles that are used to perform several 1D-GRAs and obtain an averaged site response and related variability. The proposed method is implemented on data recorded at two downhole array sites with different subsurface soil conditions: a soft soil site on Treasure Island (California, United States of America) and a rock outcrop site in Cadarache (South-East France). We show that synthetic surface-to-borehole transfer functions from 1D-GRAs provide an acceptable fit to the empirical transfer functions from low-motion earthquake records and succeed in reproducing most of the site-specific seismic response variability. The remaining mismatch between transfer functions is likely due to insufficient precision on the seismic bedrock and the impedance contrast. The variability in site response is discussed with emphasis on the role of V S small-scale heterogeneities, attenuation, and input motion incidence angle in ground motion variability for the site and soil conditions at both locations.

California↗

USGS2023

Thickness of the Saudi Arabian crust

As part of a joint Saudi Geological Survey (SGS) and U.S. Geological Survey (USGS) project, we analyzed P-wave receiver functions from seismic stations covering most of the Kingdom of Saudi Arabia to map the thickness of the crust across the Arabia Plate. We present an update of crustal-thickness estimates and fill in gaps for the western Arabian Shield and the rifted margin at the Red Sea (the coastal plain), as well as the eastern Arabian Platform. We applied a conventional H-k stacking algorithm and included careful attention to stacking weights, two forms of sedimentary corrections for stations located on the Arabian Platform, and additional processing for noisy stations. We obtained useful results at 154 stations from 898 teleseismic events over a 2-year period from 1995–1997 (for non-SGS stations) and a 6-year period from 2008–2014 (for SGS stations). Average crustal thickness (that is, depth to the Mohorovičić discontinuity [Moho] below the surface) beneath the Red Sea coastal plain (the rift margin) is 29 kilometers (km), beneath the volcanic fields (known in Arabic as harra [plural] or harrat [singular]) is 35 km, beneath the Arabian Shield (excluding harrats) is 37 km, and beneath the Arabian Platform is 38 km. Crustal thinning appears not to extend east of the rift escarpment, suggesting uniform extension that is no broader at depth than at the surface. In contrast to some previous claims that the Arabian Platform crust is thicker than that of the Arabian Shield, we find no statistically significant difference between their whole crustal thicknesses. However, the average subsedimentary crustal thickness (that is, the crystalline crust) for stations on the Arabian Platform is 34 km, 3 km thinner than the crust of the Arabian Shield. Individual station P-wave (pressure) velocity and S-wave (shear) velocity ratios ( V P / V S ) are highly variable for the Arabia Plate, ranging from 1.60 to 1.97 and averaging 1.75, with a standard deviation of 0.07. There are no statistically significant differences between V P / V S ratios of the different geologic regions of Saudi Arabia. Similar V P / V S ratios, coupled with similar crustal thicknesses for harrats and the Arabian Shield, indicate that Cenozoic magmatism has contributed negligibly to crustal growth.

Professional Paper↗

USGS2000

Use of a ground-penetrating radar system to detect pre- and post-flood scour at selected bridge sites in New Hampshire, 1996-98

Ground-penetrating radar was used to measure the depth and extent of existing and infilled scour holes and previous scour surfaces at seven bridges in New Hampshire from April 1996 to November 1998. Ground-penetrating-radar survey techniques initially were used by the U.S. Geological Survey to study streambed scour at 30 bridges. Sixteen of the 30 bridges were re-surveyed where floods exceeded a 2-year recurrence interval. A 300-megahertz signal was used in the ground-penetrating radar system that penetrated through depths as great as 20 feet of water and as great as 32 feet of streambed materials. Existing scour-hole dimensions, infilled thickness, previous scour surfaces, and streambed materials were detected using ground-penetrating radar. Depths to riprap materials and pier footings were identified and verified with bridge plans. Post data-collection-processing techniques were applied to assist in the interpretation of the data, and the processed data were displayed and printed as line plots. Processing included distance normalization, migration, and filtering but processing was kept to a minimum and some interference from multiple reflections was left in the record. Of the 16 post-flood bridges, 22 ground-penetrating-radar cross sections at 7 bridges were compared and presented in this report. Existing scour holes were detected during 1996 (pre-flood) data collection in nine cross sections where scour depths ranged from 1 to 3 feet. New scour holes were detected during 1998 (post-flood) data collection in four cross sections where scour depths were as great as 4 feet deep. Infilled scour holes were detected in seven cross sections, where depths of infilling ranged from less than 1 to 4 feet. Depth of infilling by means of steel rod and hammer was difficult to verify in the field because of cobble and boulder streambeds or deep water. Previous scour surfaces in streambed materials were identified in 15 cross sections and the depths to these surfaces ranged from 1 to 10 feet below the streambed. Riprap materials or pier footings were identified in all cross sections. Calculated record depths generally agree with bridge plans. Pier footings were exposed at two bridges and steel pile was exposed at one bridge. Exposures were verified by field observations.

New Hampshire↗

USGS2024

Probabilistic assessment of postfire debris-flow inundation in response to forecast rainfall

Communities downstream of burned steep lands face increases in debris-flow hazards due to fire effects on soil and vegetation. Rapid postfire hazard assessments have traditionally focused on quantifying spatial variations in debris-flow likelihood and volume in response to design rainstorms. However, a methodology that provides estimates of debris-flow inundation downstream of burned areas based on forecast rainfall would provide decision-makers with information that directly addresses the potential for downstream impacts. We introduce a framework that integrates a 24 h lead-time ensemble precipitation forecast with debris-flow likelihood, volume, and runout models to produce probabilistic maps of debris-flow inundation. We applied this framework to simulate debris-flow inundation associated with the 9 January 2018 debris-flow event in Montecito, California, USA. When the observed debris-flow volumes were used to drive the probabilistic forecast model, analysis of the simulated inundation probabilities demonstrates that the model is both reliable and sharp. In the fully predictive model, however, in which debris-flow likelihood and volume were computed from the atmospheric model ensemble's predictions of peak 15 min rainfall intensity, I 15 , the model generally under-forecasted the inundation area. The observed peak I 15 lies in the upper tail of the atmospheric model ensemble spread; thus a large fraction of ensemble members forecast lower I 15 than observed. Using these I 15 values as input to the inundation model resulted in lower-than-observed flow volumes which translated into under-forecasting of the inundation area. Even so, approximately 94 % of the observed inundated area was forecast to have an inundation probability greater than 1 %, demonstrating that the observed extent of inundation was generally captured within the range of outcomes predicted by the model. Sensitivity analyses indicate that debris-flow volume and two parameters associated with debris-flow mobility exert significant influence on inundation predictions, but reducing uncertainty in postfire debris-flow volume predictions will have the largest impact on reducing inundation outcome uncertainty. This study represents a first step toward a near-real-time hazard assessment product that includes probabilistic estimates of debris-flow inundation and provides guidance for future improvements to this and similar model frameworks by identifying key sources of uncertainty.

California↗

USGS2026

Quantitative mineral resource assessment of lithium pegmatite deposits in the northern Appalachian orogen, USA

Lithium demand is projected to increase more than 48 times by 2040 due to electric vehicle production and other energy storage needs. Most lithium production is outside of the USA, thereby increasing supply chain vulnerability. The combined end use importance and heightened supply risk of lithium make this lightest metallic element a critical commodity to the USA. To mitigate this supply risk, the US Geological Survey is actively assessing lithium deposits in the USA. Herein, we detail an assessment for lithium-mineralized pegmatites in the US northern Appalachian Mountains. Permissive tracts were generated by cross-referencing tectonic and geologic maps and mineral occurrence data with mappable criteria derived from generalized and region-specific lithium pegmatite ore deposit models; tracts were then ranked as having high, medium, or low permissibility. Available geophysical and geochemical data were found to be of minimal utility for this deposit type at the scale of the assessment. The number of undiscovered deposits were estimated and integrated into probabilistic simulations, which included an expanded and updated global grade and tonnage model of pegmatite-hosted lithium ore. The estimated total amount of undiscovered resources for the northern Appalachian Orogen has a median value of 1,410,000 metric tons of Li 2 O when considering moderate correlation across sub-regions. At a confidence level of 90%, a resource of at least 90,000 metric tons of Li 2 O remains undiscovered, and at a 10% confidence level, a resource of as much as 7,380,000 metric tons Li 2 O remains undiscovered. After applying an up-to-date economic filter to convert median contained lithium to recoverable material, a correlated total of 900,000 metric tons of Li 2 O may be economically extractable, equating to enough Li 2 O to provide the current annual US lithium supply deficit (presently obtained through net imports) for 127 years at 2025 rates of apparent consumption. This period of provision will inevitably shorten with projected increasing consumption rates, emphasizing that further research could be completed to better delineate regions of high lithium resource potential and support exploration and domestic production.

Connecticut, Delaware, Maine, Massachusetts, New H↗

USGS1997

Analysis of minimum 7-day discharges and estimation of minimum 7-day, 2-year discharges for streamflow-gaging stations in the Brazos River basin, Texas

The minimum 7-day, 2-year discharge (7Q2) was estimated for each of 97 streamflow-gaging stations in the Brazos River Basin, Texas, with at least 10 years of unregulated or 10 years of regulated discharge data using the Pearson Type III distribution fitted by L-moment statistics of the annual minimum 7-day discharges. Forty-six stations have at least 10 years of unregulated discharge data only, 26 stations have at least 10 years of regulated discharge data only, and 25 stations have at least 10 years of unregulated and 10 years of regulated discharge data. Statistical tests were applied to the annual minimum 7-day discharges to determine if the unregulated and regulated data have trends; and to determine if the data before regulation are significantly different from the data after regulation. Results of the Mann-Kendall test indicate a significant trend at the 95-percent confidence level for 8 of the 46 unregulated-only stations, for 6 of the 26 regulated-only stations, and for 5 stations before regulation and for 10 stations after regulation of the 25 stations with both unregulated and regulated data. Results of the Wilcoxon rank-sum test indicate a significant difference between the unregulated and regulated annual minimum 7-day discharges at the 95-percent confidence level for 16 of the 25 stations with both types of data. The 7Q2 is zero for 31 of the 46 unregulated-only stations and for 9 of the 26 regulated-only stations. The 7Q2 also is zero for 9 stations before regulation and for 3 stations after regulation of the 25 stations with both unregulated and regulated data. The 7Q2 ranges from 33 to 631 cubic feet per second for the eight regulated-only stations on the mainstem Brazos River. For the seven stations on the mainstem Brazos River with at least 10 years of unregulated and 10 years of regulated discharge data, the 7Q2 ranges from 0 to 568 cubic feet per second before regulation and from 0.30 to 670 cubic feet per second after regulation. The 7Q2 increased during regulation for 19 of the 25 stations with both unregulated and regulated data. The effect of regulation for most stations generally was an increase in the magnitude of the 7Q2.

Texas↗

USGS2025

Update of the 2008 provisional Enhanced Geothermal Systems (EGS) assessment for the Great Basin, USA

In response to the Energy Act of 2020, the U.S. Geological Survey (USGS) is updating the Enhanced Geothermal Systems (EGS) resource assessment for the Great Basin, USA. The previous 2008 provisional assessment estimated how much electricity could be generated from EGS resources of the western United States using models of electric-grade heat, models of heat extraction over time, and estimates of how much rock might be stimulated to produce viable amounts of heat. Herein, a similar conceptual strategy is applied, using updated models of heat extraction as a function of fracture spacing and well distance. Previously used reservoir heat delivery models are updated to have a dependence on fracture and well spacing, potentially improving future estimates of EGS resources as ongoing research provides a better understanding about the success of reservoir stimulation as a function of geology and location. For a range of well distances (250-1000 m) and fracture spacings (1-50 m), heat extraction efficiency ranges from 25-62%, demonstrating the importance of accounting for the most likely results of proven viable fracturing technologies. Although fracturing is important, the biggest uncertainty by far in estimating the EGS resource for the Great Basin is estimating which geologic units at what depths can be stimulated sufficiently to produce geothermal energy economically and efficiently. Uncertainties in these factors yield estimates that range over two orders of magnitude with an upper limit of ~174 terawatts-thermal (TWth) produced for 30 years from the upper 7 km of the crust. This upper limit would require significant technological advances to access most of the electric-grade resource across the Great Basin. Assuming that 1% of this estimate will be accessible in the next few decades gives a resource estimate similar to that made in the 2008 provisional assessment. These estimated EGS heat extraction rates far exceed ( greater than 100x) the natural geothermal heat production rate, thus geothermal electricity production at these rates might not be sustainable unless heat is also recharged from other sources (e.g., excess solar energy when supply exceeds demand). In addition to assessment maps and cumulative estimates, the new models of fractured reservoirs developed herein can be used to estimate steady power production given a set of fractures and well spacing, and estimates can be made for setback distances to ensure no thermal interference with nearby powerplants.

Arizona, California, Idaho, Nevada, Oregon, Utah↗