Search NASA⌕ Search

DOE OSTI · 3375649

Learning earthquake ground motions via conditional generative modeling

Ren, Pu [Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)] (ORCID:000000026354385X)·Nakata, Rie [Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States); Univ. of Tokyo (Japan); University of California, Berkeley, CA (United States). International Computer Science Inst.] (ORCID:0000000291889030)·Lacour, Maxime [University of California, Berkeley, CA (United States). International Computer Science Inst.] (ORCID:0000000283164435)·Naiman, Ilan [Ben Gurion Univ. of the Negev, Beer Sheva (Israel)]·Nakata, Nori [Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States); University of California, Berkeley, CA (United States). International Computer Science Inst.; Massachusetts Inst. of Technology (MIT), Cambridge, MA (United States)] (ORCID:0000000292959416)·Song, Jialin [University of California, Berkeley, CA (United States). International Computer Science Inst.; Simon Fraser Univ., Burnaby, BC (Canada)] (ORCID:0009000313816754)·Bi, Zhengfa [Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)]·Malik, Osman Asif [Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)] (ORCID:000000034477481X)·Morozov, Dmitriy [Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)]·Azencot, Omri [University of California, Berkeley, CA (United States). International Computer Science Inst.; Ben Gurion Univ. of the Negev, Beer Sheva (Israel)]·Erichson, N. Benjamin [Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States); University of California, Berkeley, CA (United States). International Computer Science Inst.] (ORCID:0000000306673516)·Mahoney, Michael W. [Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States); University of California, Berkeley, CA (United States). International Computer Science Inst.] (ORCID:0000000179204652)

Abstract

Predicting high-fidelity ground motions for future earthquakes is crucial for seismic hazard assessment and infrastructure resilience. Conventional empirical simulations suffer from sparse sensor distribution and geographically localized earthquake locations, while physics-based methods are computationally intensive and require accurate representations of Earth structures and earthquake sources. We propose an artificial intelligence (AI) spectrogram generator, Conditional Generative Modeling for Ground Motion (CGM-GM). CGM-GM leverages earthquake magnitudes and geographic coordinates of earthquakes and sensors as inputs, when postprocessed with phase information, capturing spatially continuous Fourier amplitude spectra (FAS) as well as properties such as P and S arrivals, and waveform durations, without explicit physics constraints. This is achieved through a probabilistic autoencoder that extracts latent distributions in the time-frequency domain and variational sequential models for prior and posterior distributions. We evaluate the performance of CGM-GM using small-magnitude earthquake records from the San Francisco Bay Area, a region with high seismic risks. Here, we report that CGM-GM demonstrates potential for complementing physics-based simulations and non-ergodic empirical ground motion models, as well as shows promise in seismology and beyond.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Ren, Pu [Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)] (ORCID:000000026354385X), Nakata, Rie [Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States); Univ. of Tokyo (Japan); University of California, Berkeley, CA (United States). International Computer Science Inst.] (ORCID:0000000291889030), Lacour, Maxime [University of California, Berkeley, CA (United States). International Computer Science Inst.] (ORCID:0000000283164435), Naiman, Ilan [Ben Gurion Univ. of the Negev, Beer Sheva (Israel)], Nakata, Nori [Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States); University of California, Berkeley, CA (United States). International Computer Science Inst.; Massachusetts Inst. of Technology (MIT), Cambridge, MA (United States)] (ORCID:0000000292959416), Song, Jialin [University of California, Berkeley, CA (United States). International Computer Science Inst.; Simon Fraser Univ., Burnaby, BC (Canada)] (ORCID:0009000313816754), Bi, Zhengfa [Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)], Malik, Osman Asif [Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)] (ORCID:000000034477481X), Morozov, Dmitriy [Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)], Azencot, Omri [University of California, Berkeley, CA (United States). International Computer Science Inst.; Ben Gurion Univ. of the Negev, Beer Sheva (Israel)], Erichson, N. Benjamin [Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States); University of California, Berkeley, CA (United States). International Computer Science Inst.] (ORCID:0000000306673516), Mahoney, Michael W. [Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States); University of California, Berkeley, CA (United States). International Computer Science Inst.] (ORCID:0000000179204652). 2026-03-16. Learning earthquake ground motions via conditional generative modeling. https://doi.org/10.1038/s41467-026-70719-2

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

KEEP EXPLORING

Related reports

Advances in geophysical forensic event monitoring

Forensic analysis of man-made, non-nuclear events (such as industrial accidents, explosion experiments and mine collapses) has become more frequent and detailed owing to advancements in geophysical monitoring. Here, in this Technical Review, we demonstrate how geophysical forensic monitoring using seismic, infrasound and hydroacoustic recordings provides insights on events in the solid earth, atmosphere and underwater. Advanced techniques, including machine-learning-based models, have been developed to detect, identify and investigate these events, providing information on location, subevents, sources and explosive yield. The increase in data availability, application of advanced methods and computation and the growth of multitechnology approaches have increased the accuracy of forensic event analysis and enabled more realistic characterization of uncertainties. For example, the 2020 Beirut explosion in Lebanon demonstrated that various seismic, acoustic and other methods could be used to estimate explosive yield (and yield uncertainties) of about 1 ktonne, providing confidence in the application of these methods to smaller events where data are available. However, forensic investigations remain largely limited to known events with identified sources. Increased access to data, sophisticated analysis methods and high-resolution earth models will improve forensic event analysis further, enabling civil and scientific applications, such as localization in the search for the lost ARA San Juan submarine.

geophysics↗

Exploring for Superhot Geothermal Targets in Magmatic Settings: 2022 Field Campaign at Newberry Volcano

This paper presents preliminary results from a subset of work carried out as part of a multinational research project entitled DErisking Exploration for multiple geothermal Plays in magmatic ENvironments (DEEPEN), supported by the U.S. Department of Energy (DOE) and Geothermica, a joint effort by EU member states and associated countries. The DEEPEN project will develop a customized approach to exploration for supercritical and superhot geothermal plays in magmatic systems, which will be applied to two demonstration sites. This paper summarizes field activities carried out at the U.S. demonstration site, Newberry Volcano in central Oregon. The objective of this work effort is to refine the subsurface model of Newberry Volcano, with special focus on deeper zones including the magmatic plumbing system and other key geologic elements. New data collection included gravity and wideband magnetotelluric (MT) surveys, as well as reinstallation of a seismic network. The National Renewable Energy Laboratory (NREL) and Enthalpion Energy LLC (Enthalpion) worked with the Deschutes National Forest Fort Rock District to use a low ground disturbance method of MT deployment to collect MT data inside the caldera and other restricted areas inside the National Volcanic Monument. This opened these areas to geophysical exploration for the first time in decades. Sites along and adjacent to the south rim of the volcano constituted the primary survey objectives. A team from Lawrence Berkeley National Laboratory (LBNL), the U.S. Geological Survey (USGS), and AltaRock also began the process of reinstalling the seismic network from the AltaRock enhanced geothermal system (EGS) demonstration in anticipation of further development activities at the site. The data ingestion, reduction, and analysis phase of the project is ongoing. We are currently processing the MT and gravity data and are developing a new, highly GPU-accelerated, 3D joint MT and gravity inversion to better localize the south rim/south flank conductive target and better understand its relationship to deep heat, fluid sources, and surface extrusive features. Joint inversions, which have not yet been undertaken at Newberry, will allow us to obtain constraints on the geologic model that cannot be determined from each method in isolation, improving our ability to image key geologic features at depth.

geophysics↗

Iapetus' Geophysics: Rotation Rate, Shape, and Equatorial Ridge

Iapetus has preserved evidence that constrains the modeling of its geophysical history from the time of its accretion until now. The evidence is (a) its present 79.33-day rotation or spin rate, (b) its shape that corresponds to the equilibrium figure for a hydrostatic body rotating with a period of approximately 16 h, and (c) its high, equatorial ridge, which is unique in the Solar System. This paper reports the results of an investigation into the coupling between Iapetus' thermal and orbital evolution for a wide range of conditions including the spatial distributions with time of composition, porosity, short-lived radioactive isotopes (SLRI), and temperature. The thermal model uses conductive heat transfer with temperature-dependent conductivity. Only models with a thick lithosphere and an interior viscosity in the range of about the water ice melting point can explain the observed shape. Short-lived radioactive isotopes provide the heat needed to decrease porosity in Iapetus? early history. This increases thermal conductivity and allows the development of the strong lithosphere that is required to preserve the 16-h rotational shape and the high vertical relief of the topography. Long-lived radioactive isotopes and SLRI raise internal temperatures high enough that significant tidal dissipation can start, and despin Iapetus to synchronous rotation. This occurred several hundred million years after Iapetus formed. The models also constrain the time when Iapetus formed because the successful models are critically dependent upon having just the right amount of heat added by SLRI decay in this early period. The amount of heat available from short-lived radioactivity is not a free parameter but is fixed by the time when Iapetus accreted, by the canonical concentration of Al-26, and, to a lesser extent, by the concentration of Fe-60. The needed amount of heat is available only if Iapetus accreted between 2.5 and 5.0Myr after the formation of the calcium aluminum inclusions as found in meteorites. Models with these features allow us to explain Iapetus? present synchronous rotation, its fossil 16-h shape, and the context within which the equatorial ridge arose.

geophysics↗