Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ground Motion Models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 127 records · Page 7

Chrono-Validation of Near-Real-Time Landslide Susceptibility Models via Plugin Statistical Simulations

The idea behind any validation scheme in landslide susceptibility studies is to test whether a model calibrated on a certain data can predict an unknown dataset of the same nature (landslide presences/absences and covariates). Almost the entirety of landslide susceptibility studies are validated by subsetting a single dataset into a training and test sets. This dataset usually corresponds either to event-specific or to historical inventories. Very rarely, a multi-temporal inventory is available and, in the few cases where this condition is met, the validation practices involve training a model on a specific landslide inventory, deriving a single predictive equation and validating it on a subsequent landslide inventory. This commonly leads landslide predictive studies, even those with a strong statistical rigor, to neglect the uncertainty estimation in their modeling scheme. In statistics, validation can also be performed via statistical simulations. This means that after fitting a given model, one can generate any number of predictive functions and test their predictive skills on any type and number of unknown datasets. In this work, we take a similar direction and we apply it to model and validate three separate co-seismic inventories, including an uncertainty estimation phase. We mapped these inventories within the same area in Indonesia, for three earthquakes occurred in 2012, 2017 and 2018. Specifically, we build three event-specific Bayesian Generalize Additive Models of the binomial family. From each model we then simulate 1000 predictive realizations over the remaining two inventories, by using a plug-in scheme where all the morphometric covariates are kept fixed and only the ground motion is replaced according to the prediction target. By doing so, we introduce a new analytical tool for near-real-time landslide predictive purposes, which is able to produce a probabilistic model which stands in between the definitions of susceptibility and hazard. In fact, our model is able to accurately estimate “where” and “when” - although not “how frequently” - landslide have occurred by featuring the multitemporal information of the trigger. In our findings, the simulations are quite similar to the fitted models; and the nine combinations we analyse produce excellent performance. This result confirms the assumption that “the past is the key to the future”, as we show that the relative contribution of each variable and their interactions in each probabilistic model remains practically the same across temporal replicates. This information is not trivial because it supports the routines implemented in global near-real-time applications.

Temporal validation↗

Establishment of terrestrial reference frames by new observational techniques

It is anticipated that terrestrial reference systems for geodynamics studies which include adopted plate motion models will be introduced for the analysis of both LAGEOS satellite and very long baseline interferometry ranging data. One of the possible approaches involves adjustment of ground station coordinates in conjunction with solutions for Universal Time 1 (UT1) and polar motions as functions of time. Another method uses principal value decomposition to reduce the number of degrees of freedom being solved for by three. In a third alternative, the values of UT1 and polar motion derived from the available data by means of the initial set of coordinates are kept fixed, and an appropriate block of data is reanalyzed using the previously determined values of UT1 and polar motion so that a new set of station coordinates can be derived.

Bender, P. L.↗

Multi‐Scale Seismic Imaging of the Ridgecrest, CA, Region With Waveform Inversion of Regional and Dense Array Data

Abstract We develop an inversion procedure for deriving multi‐scale velocity models with waveform inversions of earthquake and ambient noise data at multi‐frequency bands recorded by regional and dense sensor configurations. The method is applied for the area around the 2019 Ridgecrest earthquake rupture zones, utilizing data recorded by regional stations and dense 2D and 1D arrays with station spacings of ∼5 km and ∼100 m, respectively. Starting with regional Vp, Vs models and locations of Ridgecrest aftershocks, the velocity models and event locations are improved iteratively by inversions of waveforms recorded by regional stations and the 2D array, using a minimum spectral element size of ∼600 m. Waveforms from local events recorded by dense 1D arrays across the M7.1 rupture zone with frequencies of up to 10 Hz are used to resolve small‐scale features of the rupture zone and shallow crust with a local spectral element size of 80 m. The refined models provide self‐consistent descriptions of the rupture zone and the shallow crust embedded in the regional structures. The results reveal pronounced low Vs and high Vp/Vs in the M6.4 and M7.1 rupture zones coinciding with concentrations of seismicity, and also around the Garlock fault and in several local basins. We also observe clear velocity contrasts across the Garlock fault with polarity reversals along strike and with depth. The obtained multi‐scale velocity models can be used to improve derivations of earthquake source properties, simulations of dynamic ruptures and ground motions, and the understanding of fault and tectonic processes in the region.

Geochemistry & Geophysics↗

Transformational Regional-Scale Earthquake Simulations with the DOE EarthQuake SIMulation Exascale Framework

Earthquakes present worldwide risk to economic and human safety. The 2023 earthquakes in Turkiye provided a reminder of the potential for catastrophic consequences with 50,700 deaths and 15.7 million people affected. The ability to predict ground motions and infrastructure damage for earthquakes continues to be a challenging problem for scientists and engineers. Until now, estimates of ground motions have been performed empirically by looking at sparse data from past earthquakes. This approach can provide statistical information on intensity amplitudes but cannot inform site-specific ground motions essential to developing the most effective resilience. Interest has grown in large-scale computational models to simulate earthquakes at regional scale. The U.S. Department of Energy EarthQuake SIMulation (EQSIM) framework was developed for regional-scale earthquake simulations at unprecedented fidelity, taking advantage of emerging GPU-accelerated systems. This article describes the EQSIM workflow and demonstrates regional-scale simulations with the new computational capability available to scientists in their quest to mitigate future disasters.

58 GEOSCIENCES↗

Feasibility of Earthbound Motion in Lunar Gravity

BACKGROUND: Marginal stability of astronaut movement while performing lunar surface tasks has been well documented, and is clearly demonstrated in videos of falls, and near falls, during Apollo Lunar Extravehicular Activities. Referencing mission reports from Apollo 15 and 16 [1, 2], suspected causes for falls include: surface conditions, visibility, and gravitational effects (hypogravity). In this preliminary test, we employ the open-source biomechanical tool OpenSim [3, 4] to analyze the impact of lunar gravity (Lg) on two object-pickup motions performed by a single shirt-sleeved subject. Specifically, we attempt to answer the following questions based on an estimation of the Center of Mass Projection (CoMP) and ground reaction force Center of Pressure (COP) as it relates to the astronaut Base of Support (BOS) for 1g and Lg conditions: 1. Is the task motion, as performed in 1g, dynamically feasible in 1g and Lg? 2. Can we make the motion dynamically feasible in Lg by slowing it down? 3. Is the Lg COP equal to the 1g COP at a theoretically predicted reduction in motion speed? METHODS AND RESULTS: To answer the first question, the gravitational acceleration in the OpenSim model is modified from a nominal 9.81 m/s2 to 1.64 m/s2, and the 1g joint trajectory is input to an OpenSim based method [5] for estimating ground reaction forces and moments. From this method, the position of the CoMP and COP can be estimated and checked to see if they remain within a simulated BOS formed from the footprint of the OpenSim model to determine whether the motion is dynamically feasible. As expected, both of the motions were estimated to be feasible in a 1g environment, however, both motions had periods of infeasibility in Lg. It is well known that crew members make adjustments to motion trajectories in altered gravity fields to maintain balance. As a first step, we considered the simple adjustment of slowing the motion in Lg by a constant factor. This was accomplished by scaling the time stamps in the motion trajectory file by that factor. For the two motions considered, it was found that scale factors of 1.3 and 1.4 kept the COP just within the subject BOS. The CoMP is unchanged by the gravity field. Simple analysis of an inverted pendulum in the Lg environment, which generalizes to a general multibody system, leads to a theoretical prediction that a reduction in speed factor of √1g/Lg, or 2.445, will make the COP trajectory in Lg equivalent to that in 1g. When the above procedure was performed with a factor of 2.445, the estimated COP in Lg, was observed to be very close to that in 1g. In summary, we have developed a method for estimating the CoMP and COP in Lg, for subject motion collected in 1g. We believe this method can prove to be a valuable check and balance for simulated Lg training and testing by exposing potential simulator-induced artifacts that make the simulated task motion seem possible, when in fact, it would violate the above criteria. We also note that a reduction in task speed should tend the task motion towards stability, with a theoretical slowdown factor of √1g/Lg making the motion stability equal to that in 1g according to the CoMP and COP criteria.

R K Huffman↗

Modeling and Simulation for Exercise Vibration Isolation and Stabilization System Design

The microgravity environment that crew members experience on orbit presents a well-known health challenge, particularly when it comes to loss of muscle and bone mass. To counteract these negative effects, exercise countermeasures play a critical role in the daily routine of the crew on the International Space Station (ISS). To help inform requirements for upcoming exploration missions such as the Gateway Program, a new device, called the European Enhanced Exploration Exercise Device (E4D), is being built by the European Space Agency through their contractor, the Danish Aerospace Company. The E4D is being demonstrated on the ISS and is unique from other current exercise devices in that it provides four separate modalities in a single device: resistive, cycle ergometry, seated aerobic rowing, and rope pulling. To support the integration of E4D on ISS, a passive Vibration Isolation and Stabilization (VIS) system was required by NASA, and this responsibility was given to the Johnson Space Center. This paper describes the end-to-end process of modeling, simulation, and analysis used to inform the mechanical design of the VIS system. The process begins with the collection of representative exerciser motion capture (MoCap) through ground-based testing with the developmental E4D in both the Prototype Immersive Technology (PIT) Laboratory and the Active Response Gravity Offload System (ARGOS) facility at the Johnson Space Center. These collected MoCap data are processed through human biomechanics modeling to create forcing functions as input to a multibody dynamics simulation of the combined E4D/VIS system, with numerous resulting outputs. These outputs include microgravity accelerations, overall system displacements, internal and transmitted loads, as well as collisions. Microgravity accelerations are compared for compliance against ISS requirements while displacements are used for sway space determinations on the design and volumetric constraints within the targeted module. Internal loads are supplied to the supporting stress analysis teams and external loads for structural loads and dynamics teams, both at NASA and ESA. Finally, contact and clearance analysis is performed using the simulation to eliminate potential design issues. To ensure that the elements of the multibody simulation were verified and validated, correlation against multiple VIS related ground hardware testbeds was performed and characterized. In addition to the isolation part of the VIS problem, stabilization is also key to the integrated performance and evaluation. Due to the difficulty in defining ISS requirements in this area, the stability of the exerciser was inspected via analysis. Loss of balance was defined analytically as occurring when the resultant force vector acting on the exerciser lies outside the base of support of the feet. Both the VIS and E4D teams have recently gone through their Critical Design Reviews (CDRs) and the iterative model-based approach has been integral to inform mechanical design, particularly in the case of the VIS. This same end-to-end approach is now being applied for the Gateway Program, where an Exploration Exercise Device (EED) derived from the E4D and notional VIS for the device are under concept development.

Countermeasures↗

IMERG V06: Changes to the Morphing Algorithm

As the US Science Team’s globally gridded precipitation product from the NASA/JAXA Global Precipitation Measurement (GPM) mission, the Integrated Multi-satellitE Retrievals for GPM (IMERG) estimates the surface precipitation rates at 0.1° every half-hour using spaceborne sensors for various scientific and societal applications. One key component of IMERG is the morphing algorithm, which uses motion vectors to perform quasi-Lagrangian interpolation to fill in gaps in the passive microwave precipitation field using motion vectors. Up to IMERG V05, the motion vectors were derived from the large-scale motions of infrared observations of cloud tops. This study details the changes introduced in IMERG V06 to derive motion vectors from large-scale motions of selected atmospheric variables in numerical models, which allow IMERG estimates to be extended from the 60°N/S latitude band to the entire globe. Evaluation against both instantaneous passive microwave retrievals and ground measurements demonstrates the general improvement in the precipitation field of the new approach. Most of the model variables tested exhibited similar performance, but total precipitable water vapor was chosen as the source of the motion vectors for IMERG V06 due to its competitive performance and global completeness. Continuing assessments will provide further insights into possible refinements of this revised morphing scheme in future versions of IMERG.

Jackson Tan↗

Surface elevation, snow depth, vegetation height, and color imagery from multiple UAV surveys from 2018 to 2023 across a watershed near Teller road mile marker 27, Seward Peninsula, Alaska

Multiple photogrammetric surveys using an unoccupied aerial system (UAS) were performed across a watershed along Teller road on the Seward Peninsula (Alaska) to investigate terrain, vegetation and snowpack properties. This work was led by the environmental geophysics team from Lawrence Berkeley National Laboratory as part of the Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic). The multiple photogrammetric surveys used mostly a DJI Matrice 210 UAS with a DJI X5S camera, Ground Control Points (GCPs) surveyed with a Real Time Kinematic (RTK)-GPS, and a structure from motion (SfM) technique for photogrammetric reconstruction. Digital Surface elevation Models (DSMs) were inferred at various times of year, including near peak in plant and leaf density (July 19 2017), close to peak in snow depth (April 1 2019 and April 11 2022), and at low plant and leaf density of tall after the first bare-ground date (June 9 2019). Snow and canopy height were obtained by subtracting DSMs from a Digital Terrain Model (DTM) proxy inferred from the June 9 2019 survey. The derived products include ortho-mosaiced RGB map, DSMs, a DTM proxy, snow depth (April 2019 and 2022) and canopy height (July 2017). Unprocessed and processed data products are included in this package (processing levels 0-2). Data and metadata are provided as text (*.csv), image (*.jpg), and GeoTiff (*.tif) formats. This metadata document contains a description of the survey and processing steps and the inferred products.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic) was a 15-year research effort (2012-2027) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy’s Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy’s Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Dynamic Rupture Simulations of the Mw7.2 1992 Landers,California, Earthquake

An important part of the simulation technique implemented in the Recipe is the earthquake rupture characterization based on established empirical scaling relationships between source parameters and seismic moment. The scaling relationships used in the rupture characterization are derived from recorded data and are tested against dynamic and kinematic rupture models of crustal earthquakes. Currently, a revised three-stage source scaling relationship between the source rupture area and seismic moment of inland crustal earthquakes, initially proposed by Irikura and Miyake, is being used. Recent analysis of results obtained with the Recipe for the 2016 M7.1 Kumamoto earthquake performed by several researchers, including the PI, Arben Pitarka, indicate that in order to reproduce the near-fault long-period ground motions and the permanent displacements, observed during surface-rupturing earthquakes, the recipe needs modifications to the kinematic rupture parametrization in the shallow crustal weak zone. In particular, improvements are required in characterizing the slip rate function at shallow depths, above the sesimogenic zone. The main objective of this research project was the development of a dynamic rupture model for the Mw7.2 1992 Landers, California earthquake. The dynamic rupture model was used to derive kinematic rupture parameters, such as slip, slip-rate, rupture time distributions, and slip rate functions in the weak (0-4km) and seismogenic crustal zones. In addition to the knowledge gained from the dynamic rupture modeling of Chi-chi, Ridgecrest, Kumamoto and other crustal earthquakes, studied by the PI, the rupture characterization of the Landers earthquake can be helpful in updating the conventional recipe for broadband ground motion simulations of crustal earthquakes with shallow rupture. The project emphasizes the need for performing physics-based earthquake rupture modeling which ultimately will inform further improvements of rupture characterization. In particular, the modeling of rupture dynamics for recorded earthquakes, and numerical simulations for different rupture scenarios, provides physical constrains on slip evolution in time and space, used in adjustments of empirical source scaling relationships for large crustal earthquakes.

58 GEOSCIENCES↗

What Happens Before a Southward IMF Turning Reaches the Magnetopause?

Previous observations have shown an approximately 10-15 minute time delay in the ionospheric response to solar wind directional discontinuities marked by either southward or northward interplanetary magnetic field (IMF) turnings. We have studied one southward IMF turning observed by Time History of Events and Macroscale Interactions during Substorms (THEMIS) and GOES in the dayside magnetosphere. Using a global MHD (MagnetoHydroDynamics) model, we have reproduced the magnetopause motion in this event. We find that the observed delay in the ground response can be completely explained by deceleration of the directional discontinuity in the subsolar magnetosheath. We show that the speed of the discontinuity significantly decreases in the vicinity of the magnetopause where the magnetic barrier formed during the previous northward IMF interval. The southward turning can reach the magnetopause only after complete disruption of the magnetic barrier. The disruption or dissipation occurs via magnetosheath reconnection, as confirmed by high-speed jets in the magnetosheath. The magnetopause moves sunward as the directional discontinuity transits the magnetosheath. This sunward motion is followed by the earthward motion when the discontinuity strikes the magnetopause and magnetopause reconnection begins.

Samsonov, A. A.↗

Identifying Different Classes of Seismic Noise Signals Using Unsupervised Learning

Abstract Proper classification of nontectonic seismic signals is critical for detecting microearthquakes and developing an improved understanding of ongoing weak ground motions. We use unsupervised machine learning to label five classes of nonstationary seismic noise common in continuous waveforms. Temporal and spectral features describing the data are clustered to identify separable types of emergent and impulsive waveforms. The trained clustering model is used to classify every 1 s of continuous seismic records from a dense seismic array with 10–30 m station spacing. We show that dominate noise signals can be highly localized and vary on length scales of hundreds of meters. The methodology demonstrates the complexity of weak ground motions and improves the standard of analyzing seismic waveforms with a low signal‐to‐noise ratio. Application of this technique will improve the ability to detect genuine microseismic events in noisy environments where seismic sensors record earthquake‐like signals originating from nontectonic sources.

Johnson, Christopher W.↗

Improved lithospheric attenuation structure of the Arabian Peninsula through the use of national network data

We characterize the attenuation structure of the Arabian Peninsula through the measurement of regional phase amplitudes. High-resolution is achieved by combining stations from global networks with national network data through the cooperative effort of several countries in the region, including Saudi Arabia, Oman, Iraq, and Kuwait. The result is an improved attenuation model of the crust and upper mantle for a broad frequency band that extends from 0.5 to 10 Hz. The observed attenuation is in accordance with various elements of earth structure, including plate boundary type, style of tectonism, thermo-tectonic age, and temperature. Emerging features from the model include details in the structure along the Red Sea, and improved imaging of the southern Arabian Peninsula extending north from the Gulf of Aden. Finally, the resulting attenuation model can be employed for better magnitude estimates, in isolating tectonic and structural features, and in characterizing strong ground motion in the Arabian Peninsula.

58 GEOSCIENCES↗

A flux rope followed by recurring encounters with traveling compression regions: GEOTAIL observations

This paper reports a peculiar series of events observed with Geotail on October 24, 1992 at X approximately equals -80 R(sub E). On this day a combination of a bipolar perturbation in B(sub y) and a double-peaked perturbation in B(sub z), signature of a flux rope, was observed. Following this flux rope, bipolar B(sub z) events, identifiable as Traveling Compression Regions (TCRs), were recurrently observed. The axis of the flux rope was deduced to be field-aligned. Energetic ions (electrons) streamed earthward (tailward) in the flux rope. During the TCRs, the E x B flow showed a bipolar north-then-south perturbation, consistent with expected lobe plasma motion. The TCRs corresponded to weak ground magnetic disturbances and Pi2 waves. The observations can be explained by a model in which quasi-stagnant plasmoids are recurrently generated.

Kawano, H.↗

Spinning Rocket Simulator Turntable Design

Contained herein is the research and data acquired from the Turntable Design portion of the Spinning Rocket Simulator (SRS) project. The SRS Project studies and eliminates the effect of coning on thrust-propelled spacecraft. This design and construction of the turntable adds a structural support for the SRS model and two degrees of freedom. The two degrees of freedom, radial and circumferential, will help develop a simulated thrust force perpendicular to the plane of the spacecraft model while undergoing an unstable coning motion. The Turntable consists of a ten-foot linear track mounted to a sprocket and press-fit to a thrust bearing. A two-inch high column grounded by a Triangular Baseplate supports this bearing and houses the slip rings and pressurized, air-line swivel. The thrust bearing allows the entire system to rotate under the moment applied through the chain-driven sprocket producing a circumferential degree of freedom. The radial degree of freedom is given to the model through the helically threaded linear track. This track allows the Model Support and Counter Balance to simultaneously reposition according to the coning motion of the Model. Two design factors that hinder the linear track are bending and twist due to torsion. A Standard Aluminum "C" channel significantly reduces these two deflections. Safety considerations dictate the design of all the components involved in this project.

Miles, Robert W.↗

Design Considerations for LTV HITL Testing of Pressurized Suited Crew on a Motion-Based Platform

The upcoming NASA Lunar Terrain Vehicle (LTV) will succeed the Apollo Lunar Roving Vehicle, performing both surface exploration and logistics transfer in NASA’s return to the Moon. Human-in-the-Loop (HITL) testing will play a key role in refining the design of the LTV, ensuring its usability and ability to accommodate the astronaut population. A motion-based platform and a lunar terrain and lighting model can simulate the conditions of the lunar south pole region. It can be used by NASA to conduct HITL testing in concert with HITL testing of the drivable Ground Reference Unit (GTU) concept vehicle. The dynamic motion of the platform combined with the mobility restrictions of pressurized suits could help improve NASA understanding of vehicle-suit-astronaut interfaces. NASA human factors practitioners have proposed multiple HITL test series to use the motion-based platform in conjunction with testing of the GTU in hopes that lessons learned can be applied to help select a commercial partner to develop the Artemis LTV.

LTV↗

Study to eliminate ground resonance using active controls

The effectiveness of active control blade feathering in increasing rotor body damping and the possibility to eliminate ground resonance instabilities were investigated. An analytical model representing rotor flapping and lead-lag degrees of freedom and body pitch, roll, longitudinal and lateral motion is developed. Active control blade feathering is implemented as state variable feedback through a conventional swashplate. The influence of various feedback states, feedback gain, and weighting between the cyclic controls is studied through stability and response analyses. It is shown that blade cyclic inplane motion, roll rate and roll acceleration feedback can add considerable damping to the system and eliminate ground resonance instabilities, which the feedback phase is also a powerful parameter, if chosen properly, it maximizes augmentation of the inherent regressing lag mode damping. It is shown that rotor configuration parameters, like blade root hinge offset, flapping stiffness, and precone considerably influence the control effectiveness. It is found that active control is particularly powerful for hingeless and bearingless rotor systems.

Straub, F. K.↗

Investigation of Reduced-Order Modeling for Aircraft Stability and Control Prediction

High-fidelity computational fluid dynamics tools offer the potential to approximate increments for ground-to-flight scaling effects, as well as to augment the dynamic damping derivative data for motion-based flight simulators. Unfortunately, the computational expense is currently prohibitive for populating a complete simulator database. This work investigates an existing surrogate-based, indicial response reduced-order model methodology as a means to efficiently augment a flight simulator database with high-fidelity nonlinear aerodynamic damping derivatives. Creation of the reduced-order model is based on the superposition integrals of the step response with the derivative of its corresponding input signal. Step responses are calculated using a computational grid motion approach that separates the effects of angle of attack and sideslip angle from angular rates, and rates from angle of attack and sideslip. It is demonstrated that the transients produced during the start of a forced-oscillation motion are captured by the reduced-order model to the level of fidelity of a comparable computational solution. Aerodynamic coefficients computed within minutes by the reduced-order model for an aircraft undergoing an 18-second half Lazy-8 maneuver and a 25-second Immelmann turn maneuver are compared with those from full computational flight solutions that required days to complete. Finally, a cost-benefit assessment is included that demonstrates a compelling advantage for this approach. d for maneuvering, flexible vehicles.

Frink, Neal T.↗

Real-time reconstruction of ground motion during small magnitude earthquakes: A pilot study

This study presents a pilot investigation into a novel method for reconstructing real-time ground motion during small magnitude earthquakes (M < 4.5), removing the need for computationally expensive source characterization and simulation processes to assess ground shaking. Small magnitude earthquakes, which occur frequently and can be modeled as point sources, provide ideal conditions for evaluating real-time reconstruction methods. Utilizing sparse observation data, the method applies the Gappy Auto-Encoder (Gappy AE) algorithm for efficient field data reconstruction. This is the first study to apply the Gappy AE algorithm to earthquake ground motion reconstruction. Numerical experiments conducted with SW4 simulations demonstrate the method’s accuracy and speed across varying seismic scenarios. The reconstruction performance is further validated using real seismic data from the Berkeley area in California, USA, demonstrating the potential for practical application of real-time earthquake data reconstruction using Gappy AE. As a pilot investigation, it lays the groundwork for future applications to larger and more complex seismic events.

58 GEOSCIENCES↗