Engineering Papers⌕ Search

Engineering topics

Donnellan, Andrea

Publications and source records attributed to Donnellan, Andrea.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Time Series Foundation Models and Deep Learning Architectures for Earthquake Temporal and Spatial Nowcasting

Advancing the capabilities of earthquake nowcasting, the real-time forecasting of seismic activities, remains crucial for reducing casualties. This multifaceted challenge has recently gained attention within the deep learning domain, facilitated by the availability of extensive earthquake datasets. Despite significant advancements, the existing literature on earthquake nowcasting lacks comprehensive evaluations of pre-trained foundation models and modern deep learning architectures; each focuses on a different aspect of data, such as spatial relationships, temporal patterns, and multi-scale dependencies. This paper addresses the mentioned gap by analyzing different architectures and introducing two innovative approaches called Multi Foundation Quake and GNNCoder. We formulate earthquake nowcasting as a time series forecasting problem for the next 14 days within 0.1-degree spatial bins in Southern California. Earthquake time series are generated using the logarithm energy released by quakes, spanning 1986 to 2024. Our comprehensive evaluations demonstrate that our introduced models outperform other custom architectures by effectively capturing temporal-spatial relationships inherent in seismic data. The performance of existing foundation models varies significantly based on the pre-training datasets, emphasizing the need for careful dataset selection. However, we introduce a novel method, Multi Foundation Quake, that achieves the best overall performance by combining a bespoke pattern with Foundation model results handled as auxiliary streams.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The NASA ISRO SAR (NISAR) Mission - Validation of Science Measurement Requirements

The NASA ISRO Synthetic Aperture Radar (NISAR) is scheduled for launch early in 2024 from the Satish Dhawan Space Centre (SDSC), at Sriharikota, near Chennai, India. This mission is the result of a collaboration between NASA and Indian Space Research Organization (ISRO), where NASA has contributed elements of the mission such as an L-band SAR, and ISRO has contributed other elements, such as an S-band SAR. After successful launch, the NISAR mission will collect left-looking L-band SAR data over most of the Earth’s land areas twice during every 12-day exact repeat orbit. (once while in an ascending orbit direction and once while in a descending orbit direction). NASA and ISRO have individual and joint requirements on the mission that include the performance of the imaging radars onboard the spacecraft. For example, NASA must demonstrate that this L-band SAR will achieve a set of identified science measurement accuracy requirements that span Ecosystem science, Solid Earth science, and Cryosphere science disciplines. Likewise, ISRO has several applications objectives on both the L-band and S-band data from NISAR that the ISRO science team and project will be developing and testing. Pre-launch and post-launch activities have been planned to validate that these requirements are met. Here, we will discuss how the NASA plans are being executed and will present any initial results at the conference.

Chapman, Bruce↗
DOE OSTI2022

Optimizing Earthquake Nowcasting With Machine Learning: The Role of Strain Hardening in the Earthquake Cycle

Abstract Nowcasting is a term originating from economics, finance, and meteorology. It refers to the process of determining the uncertain state of the economy, markets or the weather at the current time by indirect means. In this paper, we describe a simple two‐parameter data analysis that reveals hidden order in otherwise seemingly chaotic earthquake seismicity. One of these parameters relates to a mechanism of seismic quiescence arising from the physics of strain‐hardening of the crust prior to major events. We observe an earthquake cycle associated with major earthquakes in California, similar to what has long been postulated. An estimate of the earthquake hazard revealed by this state variable time series can be optimized by the use of machine learning in the form of the Receiver Operating Characteristic skill score. The ROC skill is used here as a loss function in a supervised learning mode. Our analysis is conducted in the region of 5° × 5° in latitude‐longitude centered on Los Angeles, a region which we used in previous papers to build similar time series using more involved methods (Rundle & Donnellan, 2020, https://doi.org/10.1029/2020EA001097 ; Rundle, Donnellan et al., 2021, https://doi.org/10.1029/2021EA001757 ; Rundle, Stein et al., 2021, https://doi.org/10.1088/1361-6633/abf893 ). Here we show that not only does the state variable time series have forecast skill, the associated spatial probability densities have skill as well. In addition, use of the standard ROC and Precision (PPV) metrics allow probabilities of current earthquake hazard to be defined in a simple, straightforward, and rigorous way.

58 GEOSCIENCES↗
DOE OSTI2022

Earthquake Nowcasting with Deep Learning

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.

Fox, Geoffrey Charles (ORCID:0000000310171391)↗
DOE OSTI2021

Nowcasting Earthquakes: Imaging the Earthquake Cycle in California With Machine Learning

We propose a new machine learning-based method for nowcasting earthquakes to image the time-dependent earthquake cycle. The result is a timeseries that may correspond to the process of stress accumulation and release. The timeseries are constructed by using principal component analysis of regional seismicity. The patterns are found as eigenvectors of the cross-correlation matrix of a collection of seismicity timeseries in a coarse grained regional spatial grid (pattern recognition via unsupervised machine learning). The eigenvalues of this matrix represent the relative importance of the various eigenpatterns. Using the eigenvectors and eigenvalues, we compute the weighted correlation timeseries of the regional seismicity. This timeseries has the property that the weighted correlation generally decreases prior to major earthquakes in the region, and increases suddenly just after a major earthquake occurs. As in a previous paper, we find that this method produces a nowcasting timeseries that resembles the hypothesized regional stress accumulation and release process characterizing the earthquake cycle. We then address the problem of whether the timeseries contain information regarding future large earthquakes. For this, we compute a receiver operating characteristic and determine the decision thresholds for several future time periods of interest (optimization via supervised machine learning). We find that signals can be detected that can be used to characterize the information content of the timeseries. These signals may be useful in assessing present and near-future seismic hazards.

58 GEOSCIENCES↗
DOE OSTI2020

Nowcasting Earthquakes in Southern California With Machine Learning: Bursts, Swarms, and Aftershocks May Be Related to Levels of Regional Tectonic Stress

Abstract Seismic bursts in Southern California are sequences of small earthquakes strongly clustered in space and time and include seismic swarms and aftershock sequences. A readily observable property of these events, the radius of gyration ( R G ), allows us to connect the bursts to the temporal occurrence of the largest M ≥ 7 earthquakes in California since 1984. In the Southern California earthquake catalog, we identify hundreds of these potentially coherent space‐time structures in a region defined by a circle of radius 600 km around Los Angeles. We compute R G for each cluster then filter them to identify those bursts with large numbers of events closely clustered in space, which we call “compact” bursts. Our basic assumption is that these compact bursts reflect the dynamics associated with large earthquakes. Once we have filtered the burst catalog, we apply an exponential moving average to construct a time series for the Southern California region. We observe that the R G of these bursts systematically decreases prior to large earthquakes, in a process that we might term “radial localization.” The R G then rapidly increases during an aftershock sequence, and a new cycle of “radial localization” then begins. These time series display cycles of recharge and discharge reminiscent of seismic stress accumulation and release in the elastic rebound process. The complex burst dynamics we observe are evidently a property of the region as a whole, rather than being associated with individual faults. This new method allows us to improve earthquake nowcasting, which is a technique to evaluate the current state of hazard in a seismically active region.

58 GEOSCIENCES↗
DOE OSTI2020

Characterizing the uppermost 100 m structure of the San Jacinto fault zone southeast of Anza, California, through joint analysis of geological, topographic, seismic and resistivity data

We present results from complementary geological, topographic, seismic and electrical resistivity surveys at the Sagebrush Flat (SGB) site along the Clark fault (CF) strand of the San Jacinto fault zone trifurcation area southeast of Anza, California. Joint interpretation of these data sets, each with unique spatiotemporal sensitivities, allow us to better characterize the shallow (<100 m) fault zone at this structurally complex site. Geological mapping at the surface shows the CF has three main subparallel strands within a <100 m zone with varying degrees of rock damage. These strands intersect units of banded gneiss and tonalite, and various sedimentary units. Near the surface, the weathered but more intact tonalite and gneiss to the southwest have relatively high VP. The low-lying flat sedimentary basins around the two southwestern-most CF strands and elevated damaged gneiss to the northeast have lowest VP <500 m s–1. The high relief of the northeast gneiss unit may in part be explained by its extensive damage and inferred increased relative rock uplift. Resistivity imaging shows the unconsolidated dry basin sediments (maximum >1300 Ohm.m) contrasted against the compacted fine-grained (potentially wet) materials within the CF core and the Bautista Formation (minimum <40 Ohm.m), which is slightly elevated above the flat basins. The inverse relationship between VP (increases) and resistivity (decreases) in the uppermost ~15 m can be characterized as log–log linear with slopes of –2.6 to –4. At depths >30 m, the velocity heterogeneity near the surface merges into larger-scale structures that are generally slower on the northeast side of the CF core compared to the southwest side (as much as ~40 per cent reduction in average VP). A previous study revealed a 20–37 per cent variability in peak ground velocities across the SGB site from local earthquakes. The upper end of that range is associated with the near-surface unconsolidated sedimentary basins and northeast damaged gneiss unit. Preliminary analysis of time-dependent topography mostly shows effects of changing vegetation and anthropogenic activity.

Geochemistry & Geophysics↗