An official website of the United States government
Abstract Numerical simulations of seismic wave propagation are crucial for investigating velocity structures and improving seismic hazard assessment. However, standard methods such as finite difference or finite element are computationally expensive. Recent studies have shown that a new class of machine learning models, called neural operators, can solve the elastodynamic wave equation orders of magnitude faster than conventional methods. Full waveform inversion is a prime beneficiary of the accelerated simulations. Neural operators, as end‐to‐end differentiable operators, combined with automatic differentiation, provide an alternative approach to the adjoint‐state method. State‐of‐the‐art optimization techniques built into PyTorch provide neural operators with greater flexibility to improve the optimization dynamics of full waveform inversion, thereby mitigating cycle‐skipping problems. In this study, we demonstrate the first application of neural operators for full waveform inversion on a real seismic data set, which consists of several nodal transects collected across the San Gabriel, Chino, and San Bernardino basins in the Los Angeles metropolitan area.
more »
« less
A Deep‐Learning Based Parameter Inversion Framework for Large‐Scale Groundwater Models
Abstract Hydrogeologic models generally require gridded subsurface properties, however these inputs are often difficult to obtain and highly uncertain. Parametrizing computationally expensive models where extensive calibration is computationally infeasible is a long standing challenge in hydrogeology. Here we present a machine learning framework to address this challenge. We train an inversion model to learn the relationship between water table depth and hydraulic conductivity using a small number of physical simulations. For a 31M grid cell model of the US we demonstrate that the inversion model can produce a reliable K field using only 30 simulations for training. Furthermore, we show that the inversion model captures physically realistic relationships between variables, even for relationships that were not directly trained on. While there are still limitations for out of sample parameters, the general framework presented here provides a promising approach for parametrizing expensive models.
more »
« less
- PAR ID:
- 10640848
- Publisher / Repository:
- American Geophysical Union
- Date Published:
- Journal Name:
- Geophysical Research Letters
- Volume:
- 52
- Issue:
- 8
- ISSN:
- 0094-8276
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Reconstructing the history of polar temperature from ice core water isotope (δ18O) calibration has remained a challenge in paleoclimate research, because of our incomplete understanding of various temperature–δ18O relationships. This paper resolves this classical problem in a new framework called the unified slope equations (USE), which illustrates the general relations among spatial and temporalδ18O–surface temperature slopes. The USE is applied to the Antarctica temperature change during the last deglaciation in model simulations and observations. It is shown that the comparable Antarctica-mean spatial slope with deglacial temporal slope inδ18O–surface temperature reconstruction is caused, accidentally, by the compensation responses between theδ18O–inversion layer temperature relation and the inversion layer temperature itself. Furthermore, in light of the USE, we propose that the present seasonal slope ofδ18O–inversion layer temperature is an optimal paleothermometer that is more accurate and robust than the spatial slope. This optimal slope suggests the possibility of reconstructing past Antarctic temperature changes using present and future instrumental observations. Significance StatementThis paper develops a new framework called the unified slope equations (USE) to provide, for the first time, a general relation among various spatial and temporal water isotope–temperature slopes. The application of the USE to Antarctic deglacial temperature change shows that the optimal paleothermometer is the seasonal slope of the inversion layer temperature.more » « less
-
Maini, Philip K (Ed.)The Cellular-Potts model is a powerful and ubiquitous framework for developing computational models for simulating complex multicellular biological systems. Cellular-Potts models (CPMs) are often computationally expensive due to the explicit modeling of interactions among large numbers of individual model agents and diffusive fields described by partial differential equations (PDEs). In this work, we develop a convolutional neural network (CNN) surrogate model using a U-Net architecture that accounts for periodic boundary conditions. We use this model to accelerate the evaluation of a mechanistic CPM previously used to investigatein vitrovasculogenesis. The surrogate model was trained to predict 100 computational steps ahead (Monte-Carlo steps, MCS), accelerating simulation evaluations by a factor of 562 times compared to single-core CPM code execution on CPU. Over short timescales of up to 3 recursive evaluations, or 300 MCS, our model captures the emergent behaviors demonstrated by the original Cellular-Potts model such as vessel sprouting, extension and anastomosis, and contraction of vascular lacunae. This approach demonstrates the potential for deep learning to serve as a step toward efficient surrogate models for CPM simulations, enabling faster evaluation of computationally expensive CPM simulations of biological processes.more » « less
-
Abstract PurposeTo develop and evaluate a novel method for computationally efficient reconstruction from noisy MR spectroscopic imaging (MRSI) data. MethodsThe proposed method features (a) a novel strategy that jointly learns a nonlinear low‐dimensional representation of high‐dimensional spectroscopic signals and a neural‐network‐based projector to recover the low‐dimensional embeddings from noisy/limited data; (b) a formulation that integrates the forward encoding model, a regularizer exploiting the learned representation, and a complementary spatial constraint; and (c) a highly efficient algorithm enabled by the learned projector within an alternating direction method of multipliers (ADMM) framework, circumventing the computationally expensive network inversion subproblem. ResultsThe proposed method has been evaluated using simulations as well as in vivo H and P MRSI data, demonstrating improved performance over state‐of‐the‐art methods, with about 6 fewer averages needed than standard Fourier reconstruction for similar metabolite estimation variances and up to 100 reduction in processing time compared to a prior neural network constrained reconstruction method. Computational and theoretical analyses were performed to offer further insights into the effectiveness of the proposed method. ConclusionA novel method was developed for fast, high‐SNR spatiospectral reconstruction from noisy MRSI data. We expect our method to be useful for enhancing the quality of MRSI or other high‐dimensional spatiospectral imaging data.more » « less
-
Abstract Simulating open quantum systems, which interact with external environments, presents significant challenges on noisy intermediate‐scale quantum (NISQ) devices due to limited qubit resources and noise. In this study, an efficient framework is proposed for simulating open quantum systems on NISQ hardware by leveraging a time‐perturbative Kraus operator representation of the system's dynamics. This approach avoids the computationally expensive Trotterization method and exploits the Lindblad master equation to represent time evolution in a compact form, particularly for systems satisfying specific commutation relations. The efficiency of this method is demonstrated by simulating quantum channels, such as the continuous‐time Pauli channel and damped harmonic oscillators, on NISQ trapped‐ion hardware, including IonQ Harmony and Quantinuum H1‐1. Additionally, hardware‐agnostic error mitigation techniques are introduced, including Pauli channel fitting and quantum depolarizing channel inversion, to enhance the fidelity of quantum simulations. These results show strong agreement between the simulations on real quantum hardware and exact solutions, highlighting the potential of Kraus‐based methods for scalable and accurate simulation of open quantum systems on NISQ devices. This framework opens pathways for simulating more complex systems under realistic conditions in the near term.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1029/2024GL114285
