SEARCH · Engineering Papers
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.
Goals for the project included: Testing the limits of RSTT resolving capability; Increasing accuracy of travel time predictions (used for earthquake location) through data-dictated grid refinement; and, refining subsurface images of Israel and the Eastern Mediterranean and tectonic structures in the area.
We review previous approaches to nowcasting earthquakes and introduce new approaches based on deep learning using three distinct models based on recurrent neural networks and transformers. We discuss different choices for observables and measures presenting promising initial results for a region of Southern California from 1950–2020. Earthquake activity is predicted as a function of 0.1-degree spatial bins for time periods varying from two weeks to four years. The overall quality is measured by the Nash Sutcliffe efficiency comparing the deviation of nowcast and observation with the variance over time in each spatial region. The software is available as open source together with the preprocessed data from the USGS.
A plan was developed during 1971 to determine gross tectonic plate motions along the San Andreas Fault System in California. Knowledge of the gross motion along the total fault system is an essential component in the construction of realistic deformation models of fault regions. Such mathematical models will be used in the future for studies which will eventually lead to prediction of major earthquakes. The main purpose of the experiment described is the determination of the relative velocity of the North American and the Pacific Plates. This motion being so extremely small, cannot be measured directly but can be deduced from distance measurements between points on opposite sites of the plate boundary taken over a number of years.
Papers are presented dealing with prediction and forecasting earthquakes, the Haicheng earthquake of 1975, earthquake resistivity of factory structures, and earthquake studies and actions in China.
NASA's aeronautical research not only encompasses advancing aircraft technology and aviation safety for today and the future, but it also has a branch that deals with Airborne Science. This Airborne Science branch possesses unique airborne flight vehicles (manned and unmanned) whose sole purpose is to take scientists and their experiments anywhere in the world, into nearly any climatic condition, to gather the data they need for their research. This research spans such fields as; global warming, weather analysis, air and water pollution, mineral and archaeology search, earthquake damage and prediction, and yes, hurricanes. The data we obtain will hopefully, one day, allow us to better understand our planet, its natural forces and man's interaction with them.
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.
Accurate understanding and quantification of the risk to critical infrastructure posed by future large earthquakes continues to be a very challenging problem. Earthquake phenomena are quite complex and traditional approaches to predicting ground motions for future earthquake events have historically been empirically based whereby measured ground motion data from historical earthquakes are homogenized into a common data set and the ground motions for future postulated earthquakes are probabilistically derived based on the historical observations. This procedure has recognized significant limitations, principally due to the fact that earthquake ground motions tend to be dictated by the particular earthquake fault rupture and geologic conditions at a given site and are thus very site-specific. Historical earthquakes recorded at different locations are often only marginally representative. There has been strong and increasing interest in utilizing large-scale, physics-based regional simulations to advance the ability to accurately predict ground motions and associated infrastructure response. However, the computational requirements for simulations at frequencies of engineering interest have proven a major barrier to employing regional scale simulations. In a U.S. Department of Energy Exascale Computing Initiative project, the EQSIM application development is underway to create a framework for fault-to-structure simulations. This framework is being prepared to exploit emerging exascale platforms in order to overcome computational limitations. This article presents the essential methodology and computational workflow employed in EQSIM to couple regional-scale geophysics models with local soil-structure models to achieve a fully integrated, complete fault-to-structure simulation framework. Here, the computational workflow, accuracy and performance of the coupling methodology are illustrated through example fault-to-structure simulations.
Earthquake strong-motion predictions using kinematic source modeling require knowledge of the slip-rate functions (SRFs) along the rupture and their distinct characteristics in asperities, background (off-asperity) areas and near the surface. Here, in this study we analyzed SRFs from well-validated, self-consistent, and fully dynamic rupture models from earthquake cycles obeying a rate-and-state friction law, from our companion study (Galvez et al., 2021). The shapes of SRFs in asperities are well described by the regularized Yoffe function (RYF), which has only two parameters: rise time T r and smoothing time T s , which control the generation of long- and short-period ground motions, respectively. In background areas, we demonstrate that, in addition to the primary rupture, multiple secondary ruptures may also nucleate from rupture heterogeneities related to asperities, resulting in SRFs with multiple peaks. Because it is impossible to fit a multiple-peak SRF by the single-peak RYF, we describe SRFs in background areas in an effective way by fitting their amplitude spectra with the RYF spectra. Such spectrally effective RYFs capture salient aspects of seismic-wave generation and can be used in rupture generators for strong motion prediction. We found that small T s values correlate with small characteristic weakening distances, large peak slip rates (PSRs), and large rupture velocities. T r values are larger in background areas and smaller in asperities. Within the shallow aseismic zone, T s values approximately quadruple whereas T r values approximately double. Because of this dominant T s increase, PSR values decrease in the near-surface zone. These features indicate that the generation of strong motions by the near-surface portions of the rupture is negligible in the studied scenarios.
To study the response of ex-core detectors during an earthquake, an approach using Monte Carlo code MCNP was developed. As a first step to predict the response of the ex-core detectors a detailed MCNP model of the Krsko reactor core was created. A detailed MCNP core model was used to calculate 3D power distributions inside the core. Calculated power distributions were verified by comparison to the CORD-2 calculations, which is currently used for core design calculation verification of the Krsko nuclear power plant. Due to the outside the core position of the neutron detectors, the hybrid code ADVANTG is used to generate weight windows to speed up neutron transport outside the reactor core. To be able to use ADVANTG, fixed neutron source had to be reconstructed from the criticality core calculation. It was found that the power range ex-core detector movement has a negligible effect on the value of the thermal neutron flux at the active part of the detector. However, the 5 cm radial movement of the intermediate-range detector leads to 7 % - 8 % change in thermal neutron flux within the active intermediate-range detector region. The analysis continued with the evaluation of the effect of core barrel movement on the ex-core detector response. It was determined that the 2 mm core barrel radial oscillation can lead to ∼ 4% change in thermal neutron flux within the active detector region. It was evaluated that reactor core movement inside the baffle can lead up to ∼ 1.9% and ∼ 1.3% deviation in power-range and intermediate-range detector signal respectively. Analysis showed that the mechanical movement of ex-core neutron detectors could not explain the fluctuations in the ex-core detector signal. However, combined core barrel and reactor core inside baffle oscillations could be a probable reason for the observed fluctuations in the ex-core detector signal during an earthquake. (authors)
A verified prediction of the theory of plate tectonics is that the focal mechanisms of earthquakes on the mid-ocean ridge system indicate either normal faulting on ridge segments or strike-slip faulting on transform faults. A broad region in the southern Indian Ocean which differs from typical ridge and intraplate regimes has been identified. In this region a number of large off-ridge earthquakes have occurred in the last 20 years. The mechanisms for these events, where known, all involve normal faulting. Nine of these earthquakes have been studied in detail using a formal inversion technique based on matching synthetic body waves to observed seismograms.
The mechanics of the South Iceland Seismic Zone is examined by means of boundary element modeling. The differences between the observed geometry of the major tectonically active features and that of a geometrically orthogonal ridge-transform system (which assumes that the upper part of the earth's crust is an elastic plate containing vertical cuts) are examined. It is suggested that north-south faults can accommodate transform deformation only if the faults are longer or more numerous than those observed so far. This is considered reasonable because earthquake surface rupture lengths are commonly less than the inferred fault length at depth. The South Iceland Seismic Zone is subject to sequences of large earthquake every 45-112 years. In comparing the seismic moment release derived from earthquake magnitudes with that predicted by the models, it is confirmed that the system of north-south faults can act as a transform fault.
Computational simulations have become central to the seismic analysis and design of major infrastructure over the past several decades. Most major structures are now “proof tested” virtually through representative simulations of earthquake-induced response. More recently, with the advancement of high-performance computing (HPC) platforms and the associated massively parallel computational ecosystems, simulation is beginning to play a role in increased understanding and prediction of ground motions for earthquake hazard assessments. However, the computational requirements for regional-scale geophysics-based ground motion simulations are extreme, which has restricted the frequency resolution of direct simulations and limited the ability to perform the large number of simulations required to numerically explore the problem parametric space. In this article, recent developments toward an integrated, multidisciplinary earth science-engineering computational framework for the regional-scale simulation of both ground motions and resulting structural response are described with a particular emphasis on advancing simulations to frequencies relevant to engineered systems. This multidisciplinary computational development is being carried out as part of the US Department of Energy (DOE) Exascale Computing Project with the goal of achieving a computational framework poised to exploit emerging DOE exaflop computer platforms scheduled for the 2022–2023 timeframe.
Neural phase pickers—neural networks designed and trained to pick seismic phase arrivals—have proven to be a powerful tool for developing earthquake catalogs. However, these pickers suffer from prediction inconsistency in which the results they produce change, sometimes substantially, even under a small perturbation to the input waveform. This problem has not been addressed by the developers and users of these pickers. In this study, we show how prediction inconsistency can negatively affect the completeness of earthquake catalogs developed using neural phase pickers. Further, we show that simply using a small step size for the sliding window when processing continuous waveform data and aggregating the results significantly mitigates this problem. We also highlight the importance of training datasets for increasing the consistency and other performance metrics.
Earthquake stress drops are inferred to be independent of their source depth,contradicting standard linear scaling predictions for frictional stick-slip modelsof earthquakes, assuming increasing fault normal stress due to rockoverburden. Here, we examine the scaling between averaged stress drops andincreasing normal stress for simulated earthquakes sequences in continuumrate-and-state fault models. Our models exhibit a power-law-like scaling that isweaker than the linearity predicted by traditional friction models. This resultoccurs as the fault dimension becomes increasingly larger than the earthquakenucleation scale with increasing normal stress. Consequently, the averagedbehavior of ruptures becomes increasingly dominated by conditions for rupturepropagation, reflecting more heterogeneous shear stress conditions. As naturalfaults can be considerably larger than the smallest earthquakes they host, suchweaker scaling between averaged rupture conditions and normal stress maypartially explain the lack of an inferred depth dependence of earthquake stressdrops.
Measuring and enabling resiliency of electric distribution systems with increasing weather and cyber events are important. Some of the extreme events (e.g. Earthquakes, Hurricanes) and associated paths are predicted and monitored closely in advance and allow to take pre-event proactive control actions. The Distribution Phasor Measurement Units (D-PMUs) provide new opportunities and supporting such proactive actions. Here, a synchrophasor based resiliency driven pre-event reconfiguration can ensure minimizing impact of the expected event on the power distribution system and associated performance. However, the D-PMUs will also face challenges in terms of data quality similar to the transmission PMUs. The focus of this paper is to provide data mining approaches for anomaly detection in D-PMUs and proposing resiliency-driven pre-event reconfiguration with islanding as a proactive mechanisms to minimize the impact of adverse events on system using processed synchrophasors data. Results are validated for real industrial feeders and test cases with satisfactory response
Recent advancements in high performance computing platforms and computational workflow for regional-scale simulations are enabling unprecedented modeling of fault-to-structure earthquake processes. Regional simulations resolving ground motions at frequencies relevant to engineered systems are becoming computationally viable and provide a new capability to improve understanding of the geographical distribution and intensity of risk to buildings and critical infrastructure. As computational capabilities advance, it is essential to move beyond illustrative single rupture realizations for scenario earthquake events towards the development of a full suite of rupture realizations that appropriately characterize the range of risk to building systems. The work described in this article investigates the application of a suite of fault rupture realizations with the objective of assessing near-fault, site-specific seismic demand variability for building structures. A representative high-performance regional-scale computational model is utilized to execute ground motion and building response simulations based on 18 kinematic rupture realizations of an M7 strike-slip scenario earthquake. The fault rupture models for the scenario earthquake are created by systematically perturbing the hypocenter location and stochastically generating rupture parameters (slip, rise time, rake angle) to represent a breadth of ground motion intensities resulting from the spatial and temporal variabilities of an earthquake rupture process. The resulting seismic demand variability for three-story (short period) and forty-story (long period) steel moment-resisting frame buildings is characterized in terms of the median and distribution of peak inter-story drift ratio for a range of near-fault sites. The full suite of 18 fault rupture realizations and approximately 280,000 nonlinear dynamic building simulations indicate that the three-story building undergoes higher median seismic demand and significantly greater variability of demand at a given site than the forty-story building, which has important implications for the level of certainty in predicting building performance during an earthquake. The simulations performed provide deeper insight into the relationship between fault rupture parameterization and building response, which is essential information for developing a representative suite of rupture realizations for specific earthquake scenarios.
This Web-based tsunami simulation system allows users to remotely run a model on JPL s supercomputers for a given undersea earthquake. At the time of this reporting, predicting tsunamis on the Internet has never happened before. This new code directly couples the earthquake model and the ocean model on parallel computers and improves simulation speed. Seismometers can only detect information from earthquakes; they cannot detect whether or not a tsunami may occur as a result of the earthquake. When earthquake-tsunami models are coupled with the improved computational speed of modern, high-performance computers and constrained by remotely sensed data, they are able to provide early warnings for those coastal regions at risk. The software is capable of testing NASA s satellite observations of tsunamis. It has been successfully tested for several historical tsunamis, has passed all alpha and beta testing, and is well documented for users.
Exploring the spatiotemporal distribution of earthquake activity, especially earthquake migration of fault systems, can greatly to understand the basic mechanics of earthquakes and the assessment of earthquake risk. By establishing a three-dimensional strike-slip fault model, to derive the stress response and fault slip along the fault under regional stress conditions. Our study helps to create a long-term, complete earthquake catalog. We modelled Long-Short Term Memory (LSTM) networks for pattern recognition of the synthetical earthquake catalog. The performance of the models was compared using the mean-square error (MSE). Our results showed clearly the application of LSTM showed a meaningful result of 0.08% in the MSE values. Our best model can predict the time and magnitude of the earthquakes with a magnitude greater than Mw = 6.5 with a similar clustering period. These results showed conclusively that applying LSTM in a spatiotemporal series prediction provides a potential application in the study of earthquake mechanics and forecasting of major earthquake events.