USGS Science⌕ Search
Home/Archive/USGS

USGS · 70256116

Computationally efficient emulation of spheroidal elastic deformation sources using machine learning models: a Gaussian-process-based approach

Kyle Anderson·Mengyang Gu
View original source ↗Export metadata ↓

Abstract

Elastic continuum mechanical models are widely used to compute deformations due to pressure changes in buried cavities, such as magma reservoirs. In general, analytical models are fast but can be inaccurate as they do not correctly satisfy boundary conditions for many geometries, while numerical models are slow and may require specialized expertise and software. To overcome these limitations, we trained supervised machine learning emulators (model surrogates) based on parallel partial Gaussian processes which predict the output of a finite element numerical model with high fidelity but >1,000× greater computational efficiency. The emulators are based on generalized nondimensional forms of governing equations for finite non‐dipping spheroidal cavities in elastic halfspaces. Either cavity volume change or uniform pressure change boundary conditions can be specified, and the models predict both surface displacements and cavity (pore) compressibility. Because of their computational efficiency, using the emulators as numerical model surrogates can greatly accelerate data inversion algorithms such as those employing Bayesian Markov chain Monte Carlo sampling. The emulators also permit a comprehensive evaluation of how displacements and cavity compressibility vary with geometry and material properties, revealing the limitations of analytical models. Our open‐source emulator code can be utilized without finite element software, is suitable for a wide range of cavity geometries and depths, includes an estimate of uncertainties associated with emulation, and can be used to train new emulators for different source geometries.

Explore related subjects

JGR Machine Learning and ComputationArticle

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Kyle Anderson, Mengyang Gu. 2024-07-20. Computationally efficient emulation of spheroidal elastic deformation sources using machine learning models: a Gaussian-process-based approach. https://doi.org/10.1029/2024jh000161

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Dimensionality reduction techniques for analyzing tsunami simulations in hazard assessment and forecasting applications

A wide selection of linear and non-linear dimensionality reduction techniques is evaluated on synthetic tsunami data that is the basis for tsunami hazard assessment and computational forecasting. The data were computed from earthquake rupture forecasts (ERFs) supplying initial generation conditions and a linear long-wave Green's function approach to efficiently calculate tsunami propagation. The dimensionality reduction techniques include various forms of principal component analysis, manifold learning, and neural networks. Both external (e.g., mutual information) and internal (silhouette) metrics are used to evaluate clustering in the reduced dimension, focusing on peak-nearshore tsunami amplitudes (PNTA) along a nearshore isobath and time-series (marigrams) at specific nearshore locations. This evaluation is intentionally scoped to cluster-based downsampling to produce probabilistic inundation maps. For our test case along the Nankai subduction zone, manifold learning methods produced the highest clustering metric scores for PNTA profiles among the examined dimensionality reduction techniques. The marigram results are less consistent, owing primarily to uncertainty in defining appropriate clusters in the original high-dimensional (ambient) space to establish a “ground truth.” A novel application of using the integrative ERF-tsunami approach is that features extracted in the reduced domain, particularly from manifold learning, can be mapped back to rupture zones along the fault. This yields targeted fault information associated with unique PNTA profiles and marigram characteristics, such as late-arriving waves. Further testing would be needed for application-specific workflows and more complex fault systems, but results from this single subduction zone case study support the consideration of non-linear techniques (e.g., manifold learning) as alternatives to linear methods.

Nankai megathrust↗
USGS2026

SlideDetect: Spatio-temporal landslide detection using a three-dimensional convolutional neural network

Landslides pose a serious and ongoing threat to both human lives and infrastructure worldwide; therefore, it is of interest to predict where and when landslides are likely to occur. Advances in machine learning techniques have spurred numerous studies aimed at estimating relative landslide propensity, but are limited to spatial (as opposed to temporal) prediction due to the sparsity of landslide timing data. We address this data gap by training SlideDetect, a 3-dimensional convolutional neural network (3D CNN), to identify landslides based on their spatial and temporal occurrence within multitemporal image stacks. We use an inventory of landsides triggered by the 2018 Hokkaido earthquake and two years of monthly composite optical imagery spanning this event. The model can identify not only landslide location but also landslide date with an area under the precision-recall curve (PR-AUC) of 0.84. We further present a new standard for presenting PR curve results that explicitly compares model performance at different confidence thresholds, allowing for clearer model evaluation and comparison. Our new approach to constraining landslide timing paired with this more consistent and objective method for evaluating model performance shows considerable promise, and with further application and testing, SlideDetect could enhance the data availability and tools needed to advance landslide hazard and risk assessments.

JGR Machine Learning and Computation↗
USGS2026

Localization of spatiotemporally heterogeneous subsurface flows using autoencoder-based deep learning framework for time-lapse self-potential tomography

Self-potential (SP) monitoring has emerged as a valuable method for characterizing subsurface hydrogeological features and processes due to its sensitivity to fluid-induced electrokinetic effects. Despite advancements in SP inversion, challenges remain in imaging groundwater dynamics from SP activities due to complex hydrological settings and transient noise. In this study, a deep learning autoencoder (AE)-based framework is proposed for the spatiotemporal localization of subsurface fluid movement from time-lapse SP tomography. Temporal segments of time-lapse numerical inversions were first derived from long-term SP monitoring conducted from a floodplain site in Oak Ridge, Tennessee, known for active hyporheic exchange. Subsequently, AE models based on vision transformer (ViT), convolutional long short-term memory (ConvLSTM), convolutional neural network, and temporal convolutional network were individually trained and compared on the SP tomography segments for reconstruction performance. Finally, the reconstruction error over time serves as an anomaly score to identify moments of active SP variation, whereas spatial distributions of errors within these moments are analyzed to image and localize regions associated with anomalous subsurface fluid movement. The results demonstrate that ConvLSTM- and ViT-AE are most capable for the localization task with contrasting error distributions and consistent delineation of anomalies. Applying the method to both SP arrays parallel and perpendicular to the stream produced consistent anomaly zones near a fault or karst feature, validating the robustness and generalization of the approach. These results demonstrate the potential of the proposed framework as a scalable and interpretable tool for spatiotemporal analysis of subsurface flow dynamics in complex hydrogeological systems.

Tennessee↗