Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Computational Physics”

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,661 records · Page 2Linked to original sources
USGS2026

HyFlood: A surrogate-model-based framework for compound coastal flooding

Compound coastal flooding is a major threat to low-lying coastal regions and is expected to intensify under future climate change projections. However, modeling the joint interaction of waves, storm surge, tides, and rainfall remains computationally demanding, limiting the development of fast and reliable forecast tools. Here we present HyFlood, a hybrid statistical-numerical downscaling framework capable of computing and mapping high-resolution compound flood hazards while substantially reducing the computational cost compared with fully process-based hydrodynamic modeling. HyFlood combines statistical sampling and selection algorithms with a cascade of reduced-complexity surrogate models that emulate nearshore wave transformation, surf-zone hydrodynamics, and coastal, fluvial, and pluvial flooding. The surrogate models employ machine-learning and regression algorithms applied to a low-dimensional representation of the flooding outputs, obtained through statistical dimensionality reduction. The framework is demonstrated in southern O'ahu, Hawai'i, a region exposed to elevated sea levels driven by tides, waves, and storm surge along with frequent precipitation-driven flash flooding. Validation of the surrogates against the physics-based model outputs demonstrates that HyFlood accurately reproduces daily maxima of spatially distributed flooding depths. This hybrid approach offers a scalable and efficient tool to better quantify how changes in flooding drivers translate into hazard and impact assessments, and to support compound-flood risk assessments and climate-change adaptation planning.

Hawaii↗
USGS2025

Advancing current understanding of Martian impact-generated hydrothermal systems through novel coupled modeling: Insights from Gale, Jezero, and other craters

Impact-generated hydrothermal systems are considered potentially habitable environments on Mars, Earth, and other planetary bodies for microbial life. However, there is an ongoing debate regarding what geological features on Mars provide definitive evidence for such systems. Although earlier studies have modeled hydrothermal processes in Martian craters, they often lacked integration with shock physics hydrocodes to constrain initial impact conditions. The importance of this two-code coupling was demonstrated by successfully replicating alteration signatures in the Earth's Haughton impact structure. In this study, we use a similar two-code approach, combining the iSALE hydrocode with the HYDROTHERM hydrothermal model to simulate the full evolution of impact-generated hydrothermal systems. We apply this method to craters the size of Jezero (∼50 km) and Gale (∼154 km) in diameter. Although Jezero's interior is largely buried, our results align with hypothesized hydrothermal vents and alteration minerals near central uplifts in similarly sized exposed craters, such as Toro and Auki. Furthermore, our models correspond to alteration patterns observed by the Curiosity in the lower layers of Mount Sharp, which may represent remnants of impact-driven hydrothermal activity. A key finding is that these systems may persist much longer than previously estimated. Our simulations suggest that a Jezero-sized system could remain habitable for thermophiles for approximately 720,000 years, whereas a Gale-sized system could persist for nearly 2 million years. Additionally, simulations under unsaturated crustal conditions reveal that air-dominated near-surface layers can suppress vertical fluid flow, enabling deep subsurface alteration without producing detectable mineral signatures at the surface.

Journal of Geophysical Research - Planets↗
USGS2024

Framework for mapping liquefaction hazard–Targeted design ground motions

Liquefaction-induced ground failure poses substantial challenges to geotechnical earthquake engineering design. Current approaches for designing against liquefaction hazards, as specified in most seismic provisions, focus on estimating a liquefaction factor of safety ( 𝐹⁢𝑆𝐿 ) and typically characterize earthquake loading using design parameters based on probabilistic or deterministic ground motion levels. Because 𝐹⁢𝑆𝐿 is estimated deterministically, this basis of design neglects considerable uncertainties for estimating liquefaction triggering and its consequences and results in a lack of liquefaction-specific design criteria, particularly as structural design has advanced toward risk-targeted performance objectives. This study presents a framework for developing liquefaction-targeted design criteria based on a minimum acceptable return period of liquefaction, informed by probabilistic liquefaction hazard analysis (PLHA). PLHA quantifies annualized rates of liquefaction by considering contributions from (1) the full ground-motion probability space, and (2) uncertainties in liquefaction triggering using probabilistic models. PLHA is used in this study to characterize the current, effective return periods of 𝐹⁢𝑆𝐿 ( 𝑇𝑅,𝐹⁢𝑆 ) obtained from conventional liquefaction hazard analysis (CLHA) using uniform-hazard ground motions. 𝑇𝑅,𝐹⁢𝑆 is evaluated in a parametric study of nearly 100 sites throughout the conterminous United States. The results indicate large geographic variations in acceptable liquefaction hazard levels, with implied 𝑇𝑅,𝐹⁢𝑆 ranging between approximately 1,000 to 3,000 years. To address these inconsistencies without the computational demands of full PLHA, a framework is proposed for developing a liquefaction-targeted design peak ground acceleration, 𝑃⁢𝐺⁢𝐴𝐿 , for use in liquefaction models that result in consistent liquefaction design levels across all geographic locations. The mapped 𝑃⁢𝐺⁢𝐴𝐿 is shown to be somewhat sensitive to site-specific properties, and adjustment factors are developed and presented. The proposed 𝑃⁢𝐺⁢𝐴𝐿 mapping procedure produces 𝐹⁢𝑆𝐿 estimates that are consistent with those obtained from full PLHA at a target 𝑇𝑅,𝐹⁢𝑆 , providing a promising roadmap to incorporating PLHA concepts into current liquefaction design methods.

Journal of Geotechnical and Geoenvironmental Engin↗
USGS2026

Groundwater drought in the United States: Spatial and temporal variability

Many communities and ecosystems in the United States that are dependent on groundwater are potentially adversely affected by groundwater drought. We computed yearly groundwater-drought metrics and mean groundwater levels at well locations across the conterminous United States (CONUS), using data from wells and remotely sensed and modeled Gravity Recovery and Climate Experiment Drought Monitor Data Assimilation (GRACE-DADM). We also modeled the probability of low or high human impact at each well location. The spatial distribution of groundwater-drought duration and severity from 2001 to 2020 for 1,510 wells shows longer maximum duration and higher maximum severity events in drier regions like the Southwest than in wetter regions like the Northeast. Based on 613 wells in CONUS from 1981 to 2020, there are many significant decreases in drought duration and severity in the Northeast and many significant increases in annual-mean groundwater levels. In contrast, there are many significant increases in drought metrics and decreases in mean water levels in parts of the Southeast. There are major differences in trends from 2001 to 2020 between well-based and GRACE-DADM-based groundwater metrics in some CONUS regions and a very low correlation between trends at individual locations across CONUS. A potential reason for this disparity is the low GRACE-DADM resolution (∼12 km) and the potential for a large amount of groundwater variation at the local scale. Also, GRACE-DADM represents shallow, unconfined aquifers which may not match the screened interval of the monitoring wells we evaluated. Large spatial gaps in long-term, high frequency, and quality-assured groundwater-well monitoring data present a challenge for understanding groundwater-drought variability across CONUS. Remote sensing tools such as GRACE can help but cannot fully replace well monitoring, as highlighted by our study results. Substantially more long-term monitoring wells would more accurately represent groundwater-drought trends and spatial variability across CONUS, particularly in western regions.

conterminous United States↗
USGS2024

Near-real-time earthquake-induced fatality estimation using crowdsourced data and few-shot large-language models

When a damaging earthquake occurs, immediate information about casualties (e.g., fatalities and injuries) is critical for time-sensitive decision-making by emergency response and aid agencies in the first hours and days. Systems such as the Prompt Assessment of Global Earthquakes for Response (PAGER) by the U.S. Geological Survey (USGS) were developed to provide a forecast of such impacts within about 30 min of any significant earthquake globally. However, existing disaster-induced human loss estimation systems often rely on early casualty reports manually retrieved from global traditional media, which are labor-intensive, time-consuming, and have significant time latencies. Recent approaches use keyword matching and topic modeling to identify human casualty-relevant information from social media but tend to be error-prone when dealing with complex semantics in multi-lingual text data and parsing dynamically changing and conflicting human death and injury numbers shared by various unvetted sources in social media platforms. In this work, we introduce an end-to-end framework to significantly improve the timeliness and accuracy of global earthquake-induced human loss forecasting using multi-lingual, crowdsourced social media. Our framework integrates (i) a hierarchical casualty extraction model built upon large language models, prompt design, and few-shot learning to retrieve quantitative human loss claims from social media, (ii) a physical constraint-aware, dynamic-truth discovery model that discovers the truthful human loss from massive noisy and potentially conflicting human loss claims, and (iii) a Bayesian updating loss projection model that dynamically updates the final loss estimation using discovered truths. We test the framework in real-time on a series of global earthquake events in 2021 and 2022 and show that our framework effectively automates the retrieval of casualty information faster but with comparable accuracy to those now retrieved manually by the USGS. The code associated with this work is made available at: https://github.com/SusuXu-s-Lab/Hierarchical-Earthquake-Casualty-Information-Retrieval

International Journal of Disaster Risk Reduction↗
USGS2026

CRESCENT earthquake dynamic rupture, earthquake cycle, and tsunami code verification platform

Physics-based simulations are critical for understanding natural hazards. The increasing complexity of numerical codes requires benchmark exercises to verify that different computational methods yield consistent results when solving the same governing equations. Here, we present an open-access web platform designed for the verification of earthquake dynamic rupture, seismic cycle, and tsunami simulations. The platform architecture utilizes a modular, serverless backend on Amazon Web Services (AWS) to provide scalable file processing and visualization. A lightweight static web application provides a secure interface for uploading and managing results, while the browser-based data visualization enables interactive analysis of time series and surface grid data. By using structured JavaScript Object Notation (JSON) text files to define benchmark structures, the system remains fully extensible, allowing the addition of new scenarios without modifying the underlying software logic. The platform hosts the "The Tsunami Problem Versions" (TTPV) 1 & 2, two benchmarks for 3D fully coupled earthquake dynamic rupture and tsunami generation, and provides a framework for earthquake cycle models. This community resource aims to build trust in numerical simulations and facilitate long-term collaborative code verification as modeling software continues to evolve.

Seismica↗
USGS2024

Turbidite correlation for paleoseismology

Marine turbidite paleoseismology relies on the assumption of synchronous triggering of turbidity currents by earthquake shaking to infer rupture extent and recurrence. Such inference commonly depends on age dating and correlation of the physical stratigraphy of deposits carried by turbidity currents (i.e., turbidites) across great distances. Along the Cascadia subduction zone, which lies offshore the Pacific Northwest, USA, turbidite facies in core photographs, X-ray computed tomography images, and magnetic susceptibility (MS) data exhibit differences in character over relatively short distances, which implies that not all deposits can be correlated with confidence. Thus, subjective correlation based on expected similarity over great distances and weak age constraints does not independently support paleoseismic models. We present a new method for correlating turbidites along the Cascadia margin that can yield a more objective and repeatable stratigraphic framework to underpin earthquake recurrence. We use dynamic time warping to correlate MS logs and measure correlation coefficients of core pairs to evaluate correlation strength. We then compare these measures to a distribution of correlation coefficients of randomly generated turbidite sequences and find that only a small number of core pairs can be correlated more confidently than randomly stacked turbidites. This methodology promises a more robust correlation strategy for future stratigraphic studies.

Oregon, Washington↗
USGS1998

Methods and guidelines for effective model calibration; with application to UCODE, a computer code for universal inverse modeling, and MODFLOWP, a computer code for inverse modeling with MODFLOW

This report documents methods and guidelines for model calibration using inverse modeling. The inverse modeling and statistical methods discussed are broadly applicable, but are presented as implemented in the computer programs UCODE, a universal inverse code that can be used with any application model, and MODFLOWP, an inverse code limited to one application model. UCODE and MODFLOWP perform inverse modeling, posed as a parameter-estimation problem, by calculating parameter values that minimize a weighted least-squares objective function using nonlinear regression. Minimization is accomplished using a modified Gauss-Newton method, and prior, or direct, information on estimated parameters can be included in the regression. Inverse modeling in many fields is plagued by problems of instability and nonuniqueness, and obtaining useful results depends on (1) defining a tractable inverse problem using simplifications appropriate to the system under investigation and (2) wise use of statistics generated using calculated sensitivities and the match between observed and simulated values, and associated graphical analyses. Fourteen guidelines presented in this work suggest ways of constructing and calibrating models of complex systems such that the resulting model is as accurate and useful as possible.

Water-Resources Investigations Report↗
USGS1990

Preconditioned Conjugate-Gradient 2 (PCG2), a computer program for solving ground-water flow equations

This report documents PCG2: a numerical code to be used with the U.S. Geological Survey modular three-dimensional, finite-difference, ground-water flow model. PCG2 uses the preconditioned conjugate-gradient method to solve the equations produced by the model for hydraulic head. Linear or nonlinear flow conditions may be simulated. PCG2 includes two reconditioning options: modified incomplete Cholesky preconditioning, which is efficient on scalar computers; and polynomial preconditioning, which requires less computer storage and, with modifications that depend on the computer used, is most efficient on vector computers. Convergence of the solver is determined using both head-change and residual criteria. Nonlinear problems are solved using Picard iterations. This documentation provides a description of the preconditioned conjugate gradient method and the two preconditioners, detailed instructions for linking PCG2 to the modular model, sample data inputs, a brief description of PCG2, and a FORTRAN listing.

Water-Resources Investigations Report↗
USGS2024

Recent applications of the USGS National Crustal Model for Seismic Hazard Studies

The U.S. Geological Survey is developing the National Crustal Model (NCM) for seismic hazard studies to facilitate modeling site, path, and source components of seismic hazard across the conterminous United States. The NCM is composed of a 1km grid of geophysical profiles, extending from the Earth’s surface into the upper mantle. It is constructed from a threedimensional (3D) geologic framework and geophysical rules that use (1) a petrologic and mineral physics database; (2) a 3D temperature model; and (3) a calibrated rock type- and age-dependent porosity model. Parameters needed to estimate site response for existing ground motion models (GMMs), including the time-averaged velocity in the upper 30 meters (VS30), the depths to 1.0 and 2.5 km/s shear-wave velocity (Z1.0 and Z2.5), and sediment thickness, can be computed from the NCM. As GMMs continue to improve in the future, other metrics could also be extracted or derived from the NCM, such as fundamental period, site attenuation (ko), a fully frequency-dependent site response function, or 3D geophysical volumes for wavefield simulations. Application of the NCM may also benefit other aspects of seismic hazard analysis, including better accounting for path-dependent attenuation and geometric spreading, more accurate estimation of earthquake source properties such as hypocentral location and stress drop, and calculation of crustal strength profiles that inform estimates of the base of seismicity.

conterminous United States↗
USGS2026

pySATSI: A Python package for computing focal mechanism stress inversions

We introduce pySATSI, a Python package for computing earthquake focal mechanism stress inversions. This algorithm can handle a wide variety of types of stress inversion problems with a single script and can duplicate many capabilities of preceding methodologies. We also add new capabilities that include spatiotemporally variable inversion grids, damped stress estimates for clusters with few or no focal mechanisms, and variable fault‐plane ambiguities that the user can assign to individual events. In addition, we added the ability to use damped stress inversions with fault‐plane ambiguity probabilities that are weighted by fault instabilities. Our algorithm is computationally efficient with faster runtimes than previous algorithms, scales well for large datasets, and can be easily parallelized.

Seismological Research Letters↗
USGS2024

Methods to evaluate and improve the modeling of rupture directivity in assessment of seismic hazard

In recent years, there have been several advancements related to the modelling of near-source effects of earthquake rupture on strong ground shaking, leading to an improved characterization of ground motions and resulting seismic hazard. Some of these modifications have stemmed from physics-based numerical modelling of the earthquake rupture process, using physics-based dynamic rupture simulations. These contributions have led to a better understanding of how fault rupture characteristics, geometry, and the style of faulting can interact with the hypocenter-dependence on the path from source to site that may ultimately guide the development of seismic directivity models. Moving forward, the application of modern techniques can be used to incorporate these source characteristics and near-fault ground motion behavior that contribute to the azimuthally varying effects that result in rupture directivity. One example is the application of machine learning methods to support more automated integration of new predictor variables in model development and open more evaluation opportunities to access residuals. Here, we utilize several techniques to take advantage of the plethora of synthetic data and its ability to supplement preexisting trends observed in data. We showcase two examples of how models can be either developed, expanded upon, or constrained using artificial neural network model (ANNs). We evaluate the performance of the ANN with existing methods, comparing misfit, potential limitations, and ability to continue to improve upon these methods in the future. One approach uses a set of simulations with corresponding synthetic ground motions from the Southern California Earthquake Center (SCEC) CyberShake study to develop a ground motion model adapted to incorporate seismic directivity information using an ANN. This large database (TBs) enables us to train the model to capture magnitude, period, and distance variations and how these parameters relate to amplification from hypocenters located along finite-faults. In some cases, there is reduced misfit from better representing source features that aren’t included in base ground motion models that neglect hypocenter location (e.g. azimuthal variation, source-to-site terms). Another ANN method uses a shallow-layered neural network model to better fit a hypocenter-independent model. This method adjusts the median and aleatory variability to account for the averaged impact of various hypocenter distributions to fit the underlying directivity adjustment model. This method serves as a template to apply to other directivity models, improving computational efficiency and more readily enabling integration in hazard codes.

California↗
USGS2026

Water-withdrawal and groundwater-level fluctuations, physical and infiltration characteristics of Ohio River sediment, and isotope-based estimates of water sources relative to a riverbank filtration system in an alluvial aquifer, Louisville, Kentucky, 2000–20

The U.S. Geological Survey and the Louisville Water Company, from 2000 to 2020, cooperatively investigated physical and hydrologic conditions in the Ohio River alluvial aquifer and adjacent riverbed sediments in the northeast part of Louisville, Kentucky, to increase knowledge of processes affecting water sources to collector wells during operation of a riverbank filtration system (RBF). The RBF, in 2020, included an initial collector well (CW1) completed in 1999 and a tunnel and collector well system (TCWS) composed of four collector wells (CW2–CW5) completed in late 2010 that were connected by a tunnel that routed water to a pumping station and a treatment plant. Daily mean water-withdrawal rates typically ranged from about 10.4 to 20.7 million gallons per day (Mgal/d) from 2005 to 2007, then were cyclic from 2008 to 2010 and ranged from as low as about 9 to 13 Mgal/d in winter-spring months up to about 19 to 20 Mgal/d in summer-autumn months. From 2011 through 2017, peak daily mean water-withdrawal rates from CW2–CW5 fluctuated between winter-spring lows of about 20–30 Mgal/d to summer-autumn peak rates of about 41.6 Mgal/d in 2013 to about 48.2 Mgal/d in 2012. From 2017 to 2020, TCWS peak daily mean water-withdrawal rates decreased from about 45.7 Mgal/d in 2017 to about 37.5 Mgal/d in 2020 and peak daily mean rates from CW1 ranged from about 20.7 Mgal/d in 2020 to about 21.6 Mgal/d in 2018. Periods in 2019–20 with tunnel pumping station daily mean water-withdrawal rates smaller than about 30 Mgal/d had generally stable water levels in CW2 and indicated the sustainable yield of the TCWS collector wells under the prevailing hydrologic conditions. Drawdown during July–August of 2008 and 2009 was associated with CW1, the sole operating collector well. Groundwater-flow directions outside the CW1 cone of depression during July–August of 2008 and 2009 were from areas near the alluvial valley margin toward the Ohio River and in 2009 toward Harrods Creek. Collector wells CW2 through CW5 derived water during 2011–19 surveys from surface-water infiltration, as indicated by drawdown that extended beneath the Ohio River and Harrods Creek. Increases in drawdown toward the alluvial valley margins in July–August of 2013 and 2015, compared with July 2011, indicated less available recharge from those flow directions and decreased yield of groundwater from the aquifer after progressively longer TCWS water withdrawals. Drawdown at CW2–CW5 and at observation wells between them lessened in October 2019 and indicated that onshore groundwater levels can recover when collector well withdrawals decrease. Overlapping cones of depression in 2011–15 potentiometric surfaces centered around CW2, CW3, CW4, and CW5, and drawdown along divides between those wells indicated well interference between collector wells. Well-interference effects increased during 2013 and 2015 survey periods, despite decreased water-withdrawal rates compared with July 2011, and expanded to include CW1 when its use resumed. Groundwater-level fluctuations in observation wells near the Ohio River were consistent with gradual depletion of groundwater yield during tunnel pumping station operation at river elevations of about 420 feet above the National Geodetic Vertical Datum of 1929 (NGVD 29) and enhanced recharge of the aquifer when river elevations were above flood stage of about 430 feet above NGVD 29. Groundwater-level declines in observation wells near the Ohio River were most apparent in summer and some autumn months between 2012 and 2017 during periods when tunnel pumping station mean daily water-withdrawal rates were about 35 Mgal/d or more. Periods of relatively stable groundwater levels during 2019–20 in observation wells between CW2–CW5, when tunnel pumping station water-withdrawal rates were about 30–33 Mgal/d or less, also indicated that the prevailing sustainable yield was about 30 Mgal/d or less. Wells associated with CW3 and CW4 had the largest increases in slopes of best-fit lines relating water-level difference to water-withdrawal rate between 2012 and 2020, an indication of decreased efficiencies of those wells to supply water or of changes in groundwater yield around the well. Increased slopes of best-fit lines relating water-level difference to water-withdrawal rate from 2019 to 2020, when groundwater levels were otherwise relatively stable, also indicated that the specific capacity of one or more TCWS collector wells may have declined during that period. Riverbed surface elevation, roughness, and hardness changes between hydroacoustic surveys indicated active accretion and erosion of bed sediment offshore from the riverbank filtration system in 2002–4. Riverbed hydraulic conductivity computed from 2013 to 2015 in tests at 11 sites ranged from 0.2 to 230 feet per day. Ten temporary piezometer sites with adjacent operating collector wells had drawdowns that ranged from about 4 to 11.7 feet and hammer blows greater than or equal to 49 counts, which indicated a relatively great density of alluvial sediments and the potential for diminished riverbed infiltration rates. Contributions of Ohio River infiltration to collector wells estimated from deuterium and oxygen-18 isotopic ratios in representative Ohio River and groundwater samples were considerably larger in warmer months than cooler months and ranged from 80–90 percent in CW1–CW5 in May 2020 to 100 percent in May 2019 at CW2, CW4, and CW5. Groundwater contributions to water withdrawals by collector wells increased in cooler autumn months and ranged from 40 percent at CW3 and CW5 in November 2020 to 70 percent at CW2, CW4, and CW5 in November–December 2019. November 2019 and November 2020 isotopic compositions of water from the collector wells had enriched deuterium and oxygen-18 isotopic ratios consistent with mixtures of evaporation-modified water from the Ohio River and meteoric-origin groundwater.

Indiana, Kentucky↗
USGS2025

Machine learning provides reconnaissance-type estimates of carbon dioxide storage resources in oil and gas reservoirs

Oil and gas reservoirs represent suitable containers to sequester carbon dioxide (CO 2 ) in a supercritical state because they are accessible, reservoir properties are known, and they previously contained stored buoyant fluids. However, planners must quantify the relative magnitude of the CO 2 storage resource in these reservoirs to formulate a comprehensive strategy for CO 2 mitigation. Even reconnaissance-type estimates of CO 2 storage resources of known oil and gas reservoirs may require complicated calculations involving 1) estimates of recoverable oil and gas, 2) reservoir properties (depth, temperature, pressure, etc.), and 3) the physical qualities of the retained fluids. We demonstrate the application of machine learning (ML) algorithms to bypass these computations to yield more rapid estimates of CO 2 storage resources in reservoirs capable of hosting CO 2 in a supercritical state. ML algorithms are computationally efficient because they do not impose the strong assumptions on the data-generating process that standard statistical or engineering procedures require. Further, ML algorithms can capture highly complex, particularly nonlinear, relationships among predictor variables. We demonstrate the application of four different ML algorithms using data from onshore and offshore oil and gas reservoirs in Europe, and show they perform well when predictions are compared to engineering estimates. The proposed methods and models provide an effective and novel way to more rapidly and directly determine the subsurface CO 2 storage capacity of oil and gas reservoirs around the world, information that operators, researchers, and policymakers alike require to meet energy transition and decarbonization goals.

Frontiers in Enviornmental Science↗
USGS2024

Estimating groundwater level records using MOVE.1 and computing monthly percentiles from estimated groundwater records in Massachusetts

The U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, performed record extensions on groundwater levels at select wells using the Maintenance of Variance Extension type 1 (MOVE.1) method. The groundwater levels estimated from these record extensions were used to compute monthly percentiles to improve future determinations of a groundwater index. In Massachusetts, 27 of 29 short-record study wells with continuous groundwater levels between 0.8 and 8.1 years were suitable for record extensions; 37 long-record index wells were used to extend the groundwater level records at the study wells. The index well selected to pair with a study well was chosen based on Pearson correlation coefficient values; cross-correlation between the two wells; geologic and topographic similarity; and smallest distance spanning the wells. Each study well and its corresponding index well have 1 or more years of concurrent, overlapping data; a Pearson correlation coefficient that exceeded a threshold value of 0.8; and a similar aquifer type and hydrologic characteristics. Of the 29 study wells, 2 showed poor correlations with all index wells and were not considered for record extensions. Performance metrics used to assess the accuracy of the MOVE.1 models indicated that most models provided reasonable estimates of groundwater levels. Root mean square error values ranged from 0.097 to 2.292 feet, with a median of 0.536 foot. Nash-Sutcliffe efficiency coefficient values ranged from 0.623 to 0.996, with a median value of 0.759. Generally, study wells in close geographical proximity to their index well resulted in stronger model performance. The average length of groundwater level records was extended by 14.1 years to a new average of 18.1 years. The estimated groundwater level records from the MOVE.1 models resulted in an increase in the range of highest and lowest groundwater levels at 23 of 27 wells. The increase in range of groundwater levels was between 0.08 to 7.95 feet. Monthly percentiles for State drought indices were computed from the estimated MOVE.1 records and observed records through December 31, 2021. Percentiles computed from estimated records show an average groundwater level about 1.0 foot lower than observed data at the 2d percentile and 0.1 foot lower at the 30th percentile.

Massachusetts↗
USGS2025

An exploration of the relative influence of physical models for Omori’s law

Omori’s law states that the rate of aftershocks decays as a function of inverse time. There are multiple physical explanations that we reduce into a nonlinear mixed effects relation of three terms: (1) a Rate/State expression that can account for static/dynamic and viscoelastic triggering caused directly by the mainshock, (2) a fluid diffusion triggering term, and (3) a randomized secondary triggering (cascade) term. We fit free physical-model parameters to an observed aftershock sequence through two nonlinear regression methods to find the relative contributions of physics-based models in an observed aftershock sequence. Results from both methods show that Rate/State models overpredict aftershock rates by ∼0–30%. Secondary aftershocks cause a net negative contribution (seismicity rate reduction that corrects overprediction by other terms) ranging between ∼0 and 30%. All regression solutions yield negative secondary triggering contributions without being guided to do so. A physical explanation for this is that aftershock occurrence relieves stress from the crust, ultimately causing the sequence to extinguish itself. Fluid diffusion triggering contributions range from ∼0 to 20%. Diffusion processes are observed to be shorter in time than the full duration of an aftershock sequence and they are also spatially limited, diminishing their influence. Our results apply to an aftershock decay curve from the 2016 Central Apennines earthquake sequence, meaning that our specific results may not be general. Our primary conclusion is that any one physical model cannot alone fit the observed sequence as well as the combination of three we investigated.

Frontiers in Earth Science↗
USGS2025

Modeling seawater intrusion along the Alabama coastline using physical and machine learning models to evaluate the effects of multiscale natural and anthropogenic stresses

Seawater intrusion threatens groundwater resources in coastal regions, including southern Baldwin County, Alabama, where the freshwater-saltwater interface dynamics remain poorly understood. To address this gap, this study uses combined physics-based and machine-learning models to quantify seawater intrusion caused by natural (storm surges) and anthropogenic (human activities) perturbations. The long short-term memory network and wavelet analysis were used to assess vertical aquifer vulnerabilities, revealing that the shallow part of the Coastal lowlands aquifer system (CL1) in the southern Baldwin County region is more susceptible to sea level rise and groundwater extraction than deeper aquifers. Based on these findings, a cross-sectional numerical model (physics approach) for the CL1 aquifer was developed to evaluate tidal and storm surge effects, using Tropical Storm Claudette (June 2021) as a case study. Results showed that tidal fluctuations had a minimal impact on the saltwater-freshwater interface location, whereas storm surges caused substantial inland movement, with effects lasting for nine months. The steady-state version of the three-dimensional (3D) physical model predicted seawater intrusion across the entire area, and convolutional neural network-based modeling further validated the model results. The 3D physical model was also applied to a smaller area to assess human impact on the saltwater interface due to two groundwater pumping scenarios (± 50% of the baseline pumping rate). Results revealed that a 50% increase in groundwater withdrawals caused seawater to advance ~ 320 m inland, whereas a 50% reduction led to a ~ 270-meter retreat. This study highlights the vulnerability of Alabama’s shallow coastal aquifers to seawater intrusion due to storm surges and human activities, and demonstrates that combining physics-based models with machine learning approaches can improve groundwater predictions, though its accuracy depends on the availability of site-specific data.

Alabama↗
USGS2025

Deterministic physics-based earthquake sequence simulators match empirical ground-motion models and enable extrapolation to data poor regimes: Application to multifault multimechanism ruptures

We use the deterministic earthquake simulator RSQSim to generate complex sequences of ruptures on fault systems used for hazard assessment. We show that the source motions combined with a wave propagation code create surface ground motions that fall within the range of epistemic uncertainties for the Next Generation Attenuation‐West2 set of empirical models. We show the model is well calibrated where there are good data constraints, and has good correspondence in regions with fewer data constraints. We show magnitude, distance, and mechanism dependence all arising naturally from the same underlying friction. The deterministic physics‐based approach provides an opportunity for better understanding the physical origins of ground motions. For example, we find that reduced stress drops in shallow layers relative to constant stress drop with depth lead to peak ground velocities in the near field that better match empirical models. The simulators may also provide better extrapolations into regimes that are poorly empirically constrained by data because physics, rather than surface shaking data parameterizations, is underlying the extrapolations. Having shown the model is credible, we apply it to a problem where observations are lacking. We examine the case of crustal faults above a shallow subduction interface seen to break coseismically in simulations of the New Zealand fault system. These types of events were left out of consideration in the most recent New Zealand national seismic hazard model due to the modeling complexity and lack of observational data to constrain ground‐motion models (GMMs). Here, we show that in the model, by breaking up the coseismic crustal and interface rupturing fault motions into two separate subevents, and then recombining the resulting ground‐motion measures in a square‐root‐of‐sum‐of‐squares incoherent manner, we reproduce well the ground‐motion measures from the full event rupture. This provides a new method for extrapolating GMMs to more complex multifault ruptures.

Seismological Research Letters↗