Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Earthquakes”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Simulation of Physics-Based 0-10Hz Strong Motion Using High Performance Computing Supporting Refinements to Regional Ground Motion Models for the Central Eastern US

In collaboration with the U.S. Nuclear Regulatory Commission (NRC) the LLNL has developed a computationally efficient simulation platform designed to perform physics-based ground motion simulations for crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using high-performance computing. The main objective of the earthquake simulations was to use synthetic ground motion to provide constrains to refinements of existing ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-10Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model (1994). In simulations performed during the first stage of this project the Saikia’s velocity model demonstrated better performance in modelling high frequency regional wave propagation for the CEUS region recorded during the Mw5.0 November 7, 2016, Cushing Oklahoma (Taylor et al., 2017), and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes. The proposed regional 3D model includes random perturbations to the 1D background model using the stochastic scheme of Pitarka and Mellors (2021). In addition, validation analysis of the rupture generator and regional wave propagation models, using comparisons with different GMMs for Mw6.5 and Mw7.0 scenario earthquakes in the CEUS region resulted in a very good match between the simulated and empirical ground motion models. For the purposes of seismic hazard assessment at the existing and planned nuclear power plants, NRC is interested in studies aimed at improving the current ground motion models (GMM) for both Stable Continental Regions (SCR) in the Central and Eastern US and Active Crustal Regions (ACR) in the Western US. Due to lack of recorded data, these improvements require synthetic data for short fault distances and large magnitude earthquakes for which the existing recorded data is not enough to uniquely constrain the GMMs. The need for simulations and strong motion data is especially critical for the CEUS region where we do not have recorded data from potentially large damaging earthquakes with moment magnitudes 6.0 and higher. In this the project, we focused on 10Hz simulations of Mw7.0 scenario earthquakes with strike slip and thrust faulting mechanisms. We used more than 50 Mw7.0 earthquake rupture scenarios to investigate the ground motion uncertainty due to unknown earthquake rupture parameters, in particular, the slip distribution, rupture velocity, and faulting mechanism, and their implication on ground motion amplification due to forward rupture directivity effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2012

Seismic Imager Space Telescope

A concept has been developed for a geostationary seismic imager (GSI), a space telescope in geostationary orbit above the Pacific coast of the Americas that would provide movies of many large earthquakes occurring in the area from Southern Chile to Southern Alaska. The GSI movies would cover a field of view as long as 300 km, at a spatial resolution of 3 to 15 m and a temporal resolution of 1 to 2 Hz, which is sufficient for accurate measurement of surface displacements and photometric changes induced by seismic waves. Computer processing of the movie images would exploit these dynamic changes to accurately measure the rapidly evolving surface waves and surface ruptures as they happen. These measurements would provide key information to advance the understanding of the mechanisms governing earthquake ruptures, and the propagation and arrest of damaging seismic waves. GSI operational strategy is to react to earthquakes detected by ground seismometers, slewing the satellite to point at the epicenters of earthquakes above a certain magnitude. Some of these earthquakes will be foreshocks of larger earthquakes; these will be observed, as the spacecraft would have been pointed in the right direction. This strategy was tested against the historical record for the Pacific coast of the Americas, from 1973 until the present. Based on the seismicity recorded during this time period, a GSI mission with a lifetime of 10 years could have been in position to observe at least 13 (22 on average) earthquakes of magnitude larger than 6, and at least one (2 on average) earthquake of magnitude larger than 7. A GSI would provide data unprecedented in its extent and temporal and spatial resolution. It would provide this data for some of the world's most seismically active regions, and do so better and at a lower cost than could be done with ground-based instrumentation. A GSI would revolutionize the understanding of earthquake dynamics, perhaps leading ultimately to effective warning capabilities, to improved management of earthquake risk, and to improved public safety policies. The position of the spacecraft, its high optical quality, large field of view, and large field of regard will make it an ideal platform for other scientific studies. The same data could be simply reused for other studies. If different data, such as multi-spectral data, is required, additional instruments could share the telescope.

Sidick, Erkin↗
NASA NTRS2020

The Influence of Frequency and Duration Of Seismic Ground Motion On The Size Of Triggered Landslides—A Regional View

Observation, theory, and intuition all suggest that larger earthquakes should trigger larger landslides. Many factors could contribute to this, including depth-dependent shear strength or non-linearity of ground motion in soils and rock, but we hypothesize that the key characteristics of large earthquakes causing this phenomenon are (in addition to magnitude) the frequency and duration of the strong ground motion. Because of the paucity of site-specific data for detailed analysis, we take a regional approach to this question by analyzing strong-motion records and earthquake-induced landslide (EQIL) inventories from six well-documented earthquakes. Ground motion is characterized using earthquake magnitude and the median durations and frequencies (mean periods) of subsets of strong-motion records relevant to landslide triggering. EQIL inventories are characterized using the median landslide area of the entire inventory as well as the median areas of the largest 1% of the landslides and the largest 10 landslides. We then compare ground-motion characteristics with landslide size statistics to determine possible correlations. Comparisons of all earthquake- and landslide-size statistics show strong positive correlations between landslide size and (1) magnitude, (2) ground-motion duration, and (3) mean period. Although all the ground-motion measures yield highly correlated regressions, mean period appears to be the best overall predictor of landslide size. Landslide modeling using Newmark's sliding-block method also shows that longer mean periods and durations and larger magnitudes correlate strongly with increases in modeled displacements. These results support our hypothesis that increasing period and duration of seismic ground motion are the physical factors driving increased landslide sizes for larger earthquakes. Additional studies including data from a much larger set of earthquakes is needed to confirm the results of this initial study.

Hakan Tanyas↗
NASA NTRSDate not supplied

Adding GPU Support to the Markov Chain Monte Carlo Code Catmip

In geophysics, we are confronted with many under-determined inverse problems. For example, all of our observations of earthquakes are made at the Earth’s surface. So, when we try to infer how slip during an earthquake evolves in space and time, we find that there are many potential slip histories that are consistent with our limited observations and our understanding of earthquake physics. One way to approach these problems is with Bayesian analysis which allows us to infer the ensemble of all potential slip models that satisfy the observations and our prior knowledge of earthquake physics. In Bayesian analysis, our prior knowledge is known as the prior probability density function or prior PDF, the fit to the data is known as the data likelihood, and the target PDF that satisfies both the prior PDF and data likelihood is known as the posterior PDF. However, simulating the posterior PDF typically requires using Markov Chain Monte Carlo (MCMC) to draw tens of billions of random realizations of earthquake slip models, which may not be computationally feasible. To make this and similar geophysical inversions computationally tractable, we developed the Cascading Adaptive Transitional Metropolis In Parallel (CATMIP) algorithm. CATMIP is an efficient parallel Markov Chain Monte Carlo (MCMC) sampler that is used for model fitting and uncertainty quantification in geophysics. Example use cases are earthquake rupture modeling, determining mineral composition on Mars, reconstructing the history of ocean salinity, and historical earthquake relocation. CATMIP employs many parallel instances of the Metropolis algorithm for sampling in a transitioning framework. Transitioning is a process in which a set of random samples at equilibrium with a known probability density function (PDF) are used as seeds for the Markov chains to sample successive target PDFs that incrementally move the distribution from the starting seeds to the final desired PDF that describes the relative plausibility of potential values for the model parameters. The algorithm is implemented as a Master-Worker model employing MPI for communication. The worker processes are loosely coupled with global parameters periodically optimized by the master process. This provides a very high amount of parallelism with little communication between updates. During the presentation we will discuss the history of the algorithm and elaborate the earthquake rupture modeling use case for the CATMIP package. Our first step toward GPU optimization was to optimize the code for the CPU. CPU profiling revealed that most of the compute time is spent in calls to level 2 BLAS routines and calls to GSL random number generators. We revised the algorithm to employ level 3 BLAS routines instead. In our presentation we will describe how this was accomplished. Adding GPU support to CATMIP consisted mostly of replacing the calls to GSL with calls to GPU vendor-provided library routines. A small number of loops were directly implemented in CUDA. In the presentation will provide implementation details. Finally, we will discuss methods for profiling and opportunities for further optimizing GPU execution. By creating a code with the flexibility to run on either a CPU or GPU architecture, CATMIP can be used on systems ranging from large CPU-based HPC environments to single servers with GPU acceleration and everything in between.

HECC↗
DOE OSTI2022

Limited Dynamic Triggering in the Utah Region, USA

The state of Utah, USA, experiences around 3800 catalogued earthquakes per year, highlighting that the region is seismically active and susceptible to earthquakes. Following the 2002 Denali Fault (M7.9) earthquake in Alaska, the region showed an elevated seismicity rate for 3 weeks following the passage of high amplitude surface waves, suggesting that the region may be particularly susceptible to dynamic triggering. With over 23,396 faults and each fault presenting a potential fault for triggering, we systematically search for dynamic triggering throughout the state of Utah caused by large, global earthquakes with M ≥ 7. Specifically, we analyse earthquake catalogues and all available waveform data to determine statistical increases of seismicity rate following the passage of seismic arrivals. While we find instances of dynamic triggering, our results show that these events occur sparsely in the region. In total, less than 20 per cent of the 273 main shocks that occur from 2000 to the end of 2017 show a statistical indication of dynamic triggering throughout the Utah region, highlighting that dynamic triggering is limited for stresses created by transient signals from global M ≥ 7 earthquakes, with the exception being the Denali Fault (M7.9), Alaska earthquake (i.e. an instance of significant triggering).

58 GEOSCIENCES↗
DOE OSTI2023

Quantifying rupture characteristics of microearthquakes in the Parkfield Area using a high-resolution borehole network

It is well known that large earthquakes often exhibit significant rupture complexity such as well separated subevents. With improved recording and data processing techniques, small earthquakes have been found to exhibit rupture complexity as well. Studying these small earthquakes offers the opportunity to better understand the possible causes of rupture complexities. Specifically, if they are random or are related to fault properties. We examine microearthquakes (M < 3) in the Parkfield, California, area that are recorded by a high-resolution borehole network. We quantify earthquake complexity by the deviation of source time functions and source spectra from simple circular (omega-square) source models. We establish thresholds to declare complexity, and find that it can be detected in earthquakes larger than magnitude 2, with the best resolution above M2.5. Comparison between the two approaches reveals good agreement (>90 per cent), implying both methods are characterizing the same source complexity. For the two methods, 60–80 per cent (M 2.6–3) of the resolved events are complex depending on the method. The complex events we observe tend to cluster in areas of previously identified structural complexity; a larger fraction of the earthquakes exhibit complexity in the days following the M w 6 2004 Parkfield earthquake. Ignoring the complexity of these small events can introduce artefacts or add uncertainty to stress drop measurements. Focusing only on simple events however could lead to systematic bias, scaling artefacts and the lack of measurements of stress in structurally complex regions.

58 GEOSCIENCES↗
NASA NTRS2021

Capturing the Footprints of Ground Motion in the Spatial Distribution of Rainfall-Induced Landslides

The coupled effect of earthquakes and rainfall is rarely investigated in landslide susceptibility assessments although it could be crucial to predict landslide occurrences. This is even more critical in the context of early warning systems and especially in cases of extreme precipitation regimes in post-seismic conditions, where the rock masses are already damaged due to the ground shaking. Here, we investigate this concept by accounting for the legacy of seismic ground shaking in rainfall-induced landslide (RFIL) scenarios. We do this to identify whether ground shaking plays a role in the susceptibility to post-seismic rainfall-induced landslides and to identify whether this legacy effect persists through time. With this motivation, we use binary logistic regression and examine time series of landslides associated with four earthquakes occurred in Indonesia: 2012 Sulawesi (M(sub w) = 6.3), 2016 Reuleut (M(sub w) = 6.5), 2017 Kasiguncu (M(sub w) = 6.6) and 2018 Palu (M(sub w) = 7.5) earthquakes. The dataset includes one co-seismic and three post-seismic landslide inventories for each earthquake. We use the peak ground acceleration map of the last strongest earthquake in each case as a predisposing factor of landslides representing the effect of ground shaking. We observe that, at least for the study areas under consideration and in a probabilistic context, the earthquake legacy contributes to increase the post-seismic RFIL susceptibility. This positive contribution decays through time. Specifically, we observe that ground motion is a significant predisposing factor controlling the spatial distribution of RFIL in the post-seismic period 110 days after an earthquake. We also show that this effect dissipates within 3 years at most.

Landslides↗
DOE OSTI2022

Small Seismic Events in Oklahoma Detected and Located by Machine Learning–Based Models

A complete earthquake catalog is essential to understand earthquake nucleation and fault stress. Following the Gutenberg–Richter law, smaller, unseen seismic events dominate the earthquake catalog and are invaluable for revealing the fault state. The published earthquake catalogs, however, typically miss a significant number of small earthquakes. Part of the reason is due to a limitation of conventional algorithms, which can hardly extract small signals from background noise in a reliable and efficient way. To address this challenge, we utilized a machine learning method and developed new models to detect and locate seismic events. These models are efficient in processing a large amount of seismic data and extracting small seismic events. We applied our method to seismic data in Oklahoma, United States, and detected ~14 times more earthquakes compared with the standard Oklahoma Geological Survey catalog. The rich information contained in the new catalog helps better understand the induced earthquakes in Oklahoma.

58 GEOSCIENCES↗
DOE OSTI2025

A Deep Learning-Aided Workflow for Decoding the Stress Regime of Southern Nevada

The Rock Valley fault zone in southern Nevada has a notable history of seismic activity and is the site of a future direct comparison experiment of explosion and earthquake sources. This study aims to gain insight into regional tectonic processes by leveraging recent advances in seismic monitoring capabilities to elucidate the local stress regime. A crucial step in this investigation is the accurate determination of P-wave first-motion polarities, which play a vital role in resolving earthquake focal mechanisms of small earthquakes. Here, we deploy a deep learning-based method for automatic determination of first-motion polarities to vastly expand the polarity dataset beyond what has been reviewed by human analysts. By the integrating P-wave polarities with new measurements of S/P amplitude ratios, we obtain robust focal mechanism estimates for 1306 earthquakes with a local magnitude of 1 and above occurring between 2010 and 2023 in southern Nevada. We then use the focal mechanism catalog to examine the regional stress orientation, confirming an overall trans-tensional stress regime with smaller scale complexities illuminated by individual earthquake sequences. These findings demonstrate how detailed analyses of small earthquakes can provide fundamental information for understanding earthquake processes in the region and inform future experiments at the Nevada National Security Site.

58 GEOSCIENCES↗
DOE OSTI2025

Machine Learning Approaches to Predicting Induced Seismicity and Imaging Geothermal Reservoir Properties

This project developed machine learning (ML) methods, lab data sets, and field data to advance geothermal exploration and geothermal energy production. The work had three focus areas. One involved the development of ML methods to use microearthquakes (MEQs) for imaging geothermal reservoir properties and improving subsurface characterization – most importantly the evolution of permeability within the evolving reservoir. This part of the work included development of ML approaches for automated MEQ location, focal mechanism determination and identification of earthquake precursors. The second area focused on using MEQ signals generated by geothermal exploration and production to predict the relationship between fluid injection and seismicity. Here, we extended to reservoir scale our success in using ML to predict laboratory earthquakes and fault zone stress state. The third focus area was on lab experiments. Here, we developed new ML models for lab earthquake prediction and identification of precursors to failure to improve earthquake forecasting and early warning in geothermal settings. Major outcomes of our work include ML models that learn from MEQ signals during geothermal exploration and production to predict induced seismicity. MEQs occur naturally in connection with drilling and energy production. We developed ML methods to use the seismic waves from these events to characterize the elastic, hydraulic and poromechanical properties of reservoirs. Our work illuminated fracture geometry and the evolution of fracture permeability by incorporating seismic coda wave analysis and ML methods to relate fluid injection and seismicity. We significantly expanded laboratory earthquake prediction to include methods that use both passive measurements of microearthquakes within the lab fault zones and also active source acoustic measurements of fault zone elastic properties. These methods can now predict fault zone stress state, time to failure and the magnitude of lab earthquakes. Our work showed that repetitive stick- slip failure events during frictional sliding (the lab equivalent of earthquakes) are preceded by a cascade of micro-failure events that radiate energy in a manner that foretells unstable failure – manifest as laboratory MEQs. We documented a mapping between fracture properties and statistical attributes of elastic radiation. We extended existing works to geothermal reservoir scale and developed ML methods to determine reservoir permeability, fracture properties, and their evolution during geothermal energy production. An attractive feature of ML algorithms is their ability to handle big datasets and reveal patterns and correlations that may remain invisible to conventional analyses. Our work connected data from field, laboratory and intermediate scales to study permeability, stress, strength, fracture stiffness and geometry. At the field scale we used data from the Newberry Volcano field site, UtahFORGE, EGS Collab, and also the Bedretto underground research lab in Switzerland. These data sets are bridging the gap between the lab scale, theory, and reservoir scale. Our work produced plain language summaries to improve public understanding of DOE research. We also developed openly distributed ML and seismicity datasets for use by all researchers and we published connections between induced seismicity in geothermal areas and reservoir properties including permeability, fracture properties, and stress state. Our models are designed for the large data sets of induced seismicity typically associated with geothermal sites. We produced labeled event catalogs and used them on geothermal data to assess how ML can facilitate geothermal production and exploration. All datasets are available on the GDR Productivity: The project produced 32 publications in peer reviewed journals (two are in review). It supported the work of 6 PhD students, 40 conference presentations, 6 keynote talks at national meetings, and mentoring and professional development for 4 postdoctoral fellows.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

A Comprehensive Stress Drop Map From Trench to Depth in the Northern Chilean Subduction Zone

We compute stress drops for earthquakes in Northern Chile recorded between 2007 and 2021. By applying two analysis techniques, (a) the spectral ratio (SR) method and (b) the spectral decomposition (SDC) method, a stress drop map for the subduction zone consisting of 51,510 stress drop values is produced. We build an extended set of empirical Green’s functions (EGF) for the SR method by systematic template matching. Outputs are used to compare with results from the SDC approach, where we apply cell-wise obtained global EGF's to compensate for the structural heterogeneity of the subduction zone. We find a good consistency of results of the two methods. The increased spatial coverage and quantity of stress drop estimates from the SDC method facilitate a consistent stress drop mapping of the different seismotectonic domains. Albeit only small differences of median stress drop, strike-perpendicular depth sections clearly reveal systematic variations, with earthquakes at different seismotectonic locations exhibiting distinct values. In particular, interface seismicity is characterized by the lowest observed median value, whereas upper plate earthquakes show noticeably higher stress drop values. Intermediate depth earthquakes show comparatively high average stress drop and a rather strong depth-dependent increase of median stress drop. Additionally, we observe spatio-temporal variability of stress drops related to the occurrence of the two megathrust earthquakes in the study region. The presented study is the first coherent large scale 3D stress drop mapping of the Northern Chilean subduction zone. It provides an important component for further detailed analysis of the physics of earthquake ruptures.

58 GEOSCIENCES↗
DOE OSTI2026

Activating a Natural Fault Zone in the Swiss Alps

One major hurdle for understanding earthquake mechanics are observational limitations. Important phenomena like strain localisation, fault dilation, and fault healing are readily studied in rock mechanical laboratory experiments and with numerical models. At the scale of natural earthquakes, however, these phenomena are often unresolvable, even by state-of-the-art observatories. To overcome this limitation, we are currently building the Earthquake Physics Testbed at the Bedretto Underground Laboratory for Geosciences and Geoenergies (BedrettoLab), an experimental testbed where we can activate an extensively instrumented natural fault zone via hydraulic stimulation. The goal of the Fault Activation and Earthquake Rupture (FEAR) project is to induce earthquakes of up to Mw~1.0 on this exceptionally well characterised and instrumented fault zone. Here we summarize the main scientific goals and current FEAR project status, and present first results from conducted experiments. We discuss how this large-scale experimental approach may allow us to tackle both fundamental science as well as practical questions on earthquake physics, induced seismicity and seismic hazard.

Meier, Men-Andrin↗
DOE OSTI2024

A strainmeter array as the fulcrum of novel observatory sites along the Alto Tiberina Near Fault Observatory

Fault slip is a complex natural phenomenon involving multiple spatiotemporal scales from seconds to days to weeks. To understand the physical and chemical processes responsible for the full fault slip spectrum, a multidisciplinary approach is highly recommended. The Near Fault Observatories (NFOs) aim at providing high-precision and spatiotemporally dense multidisciplinary near-fault data, enabling the generation of new original observations and innovative scientific products. The Alto Tiberina Near Fault Observatory is a permanent monitoring infrastructure established around the Alto Tiberina fault (ATF), a 60 km long low-angle normal fault (mean dip 20°), located along a sector of the Northern Apennines (central Italy) undergoing an extension at a rate of about 3 mm yr –1 . The presence of repeating earthquakes on the ATF and a steep gradient in crustal velocities measured across the ATF by GNSS stations suggest large and deep (5–12 km) portions of the ATF undergoing aseismic creep. Both laboratory and theoretical studies indicate that any given patch of a fault can creep, nucleate slow earthquakes, and host large earthquakes, as also documented in nature for certain ruptures (e.g., Iquique in 2014, Tōhoku in 2011, and Parkfield in 2004). Nonetheless, how a fault patch switches from one mode of slip to another, as well as the interaction between creep, slow slip, and regular earthquakes, is still poorly documented by near-field observation. With the strainmeter array along the Alto Tiberina fault system (STAR) project, we build a series of six geophysical observatory sites consisting of 80–160 m deep vertical boreholes instrumented with strainmeters and seismometers as well as meteorological and GNSS antennas and additional seismometers at the surface. By covering the portions of the ATF that exhibits repeated earthquakes at shallow depth (above 4 km) with these new observatory sites, we aim to collect unique open-access data to answer fundamental questions about the relationship between creep, slow slip, dynamic earthquake rupture, and tectonic faulting.

58 GEOSCIENCES↗
NASA NTRS1989

The interpretation of crustal dynamics data in terms of plate motions and regional deformation near plate boundaries

A particularly detailed set of observations in the vicinity of an intraplate, thrust earthquake (M 7.4) in Argentina, indicate a cyclic pattern of deformation very similar to that reported previously for interplate earthquakes. This deformation cycle, which may be characteristic of many seismically active areas, consists of: (1) steady strain accumulation, possibly punctuated by strain reversals; (2) coseismic strain release; (3) a period of continued strain release due to afterslip (persisting for perhaps a year or so); and (4) rapid postseismic strain accumulation which decreases exponentially and grades into steady strain accumulation. Deformation associated with three earthquakes in the U.S. (1940, M7.1 Imperial Valley California; 1964, M8.4 Alaska; 1959, M7.5 Hebgen Lake, Montana) are interpreted in light of this general earthquake cycle and are used to investigate the mechanics of strain release for these events. These examples indicate that large postseismic movements can occur for strike-slip, thrust, and normal fault events, and both viscoelastic relaxation and postseismic after-slip must be incorporated in models of earthquake related deformation. The mechanics of the strain release process revealed by these examples has implications for estimating earthquake repeat times from geodetic observations near active faults.

Solomon, Sean C.↗
NASA NTRS1998

Permanent GPS Geodetic Array in Southern California

The southern California Permanent GPS Geodetic Array (PGGA) was established in the spring of 1990 to evaluate continuous Global Positioning System (GPS) measurements as a new too] for monitoring crustal deformation. Southern California is an ideal location because of the relatively high rate of tectonic deformation, the high probability of intense seismicity, the long history of conventional and space geodetic measurements, and the availability of a well developed infrastructure to support continuous operations. Within several months of the start of regular operations, the PGGA recorded far-field coseismic displacements induced by the June 28, 1992 (M(sub w)=7.3), Landers earthquake, the largest magnitude earthquake in California in the past 40 years and the first one to be recorded by a continuous GPS array. Only nineteen months later, on 17 January 1994, the PGGA recorded coseismic displacements for the strongest earthquake to strike the Los Angeles basin in two decades, the (M(sub e)=6.7) Northridge earthquake. At the time of the Landers earthquake, only seven continuous GPS sites were operating in southern California; by the beginning of 1994, three more sites had been added to the array. However, only a pair of sites were situated in the Los Angeles basin. The destruction caused by the Northridge earthquake spurred a fourfold increase in the number of continuous GPS sites in southern California within 2 years of this event. The PGGA is now the regional component of the Southern California Integrated GPS Network (SCIGN), a major ongoing densification of continuous GPS sites, with a concentration in the Los Angeles metropolitan region. Continuous GPS provides temporally dense measurements of surface displacements induced by crustal deformation processes including interseismic, coseismic, postseismic, and aseismic deformation and the potential for detecting anomalous events such as preseismic deformation and interseismic strain variations. Although strain meters yield much higher short-term resolution to a period of about 1 year, a single continuous GPS site is significantly less expensive than a single strain meter and probably has better long-term stability beyond a 1-year period. Compared to less frequent field measurements, continuous GPS provides the means to better characterize the errors in GPS position measurements and thereby obtain more realistic estimates of derived parameters such as site velocities.

Green, Cecil H.↗
NASA NTRS2002

Transients in Pacific/North American Plate Boundary Deformation: Synthesis and Modeling of GPS and Borehole Strain Observations

This is the Final Technical Report on research conducted between 1 June 1997 and 14 September 2001 entitled "Transients in Pacific/North American plate boundary deformation: Synthesis and modeling of GPS and borehole strain observations." As the project title implies, our effort involved a geodetic study of strain transients, i.e., temporal variations in deformation rates, that occur within plate boundary zones and their relationship to earthquakes and plate motions. Important transients occur during and following large earthquakes, and there are also strain transients not apparently associated with earthquakes. A particularly intriguing class of transients, for which there is a modest but growing list of examples, are preseismic anomalies. Such earthquake precursors, if further documented and understood, would have obvious importance for earthquake hazard mitigation. Because the timescales for these diverse transients range over at least 6 orders of magnitude (minutes to years), no single geodetic technique is optimum. We therefore undertook a systematic synthesis of Global Positioning Satellite (GPS) and borehole strainmeter data in three areas in California where there are adequate numbers of both types of instruments (or their equivalent): the San Francisco Bay region (within the Bay Area Regional Deformation network), southern California (within the Southern California Integrated GPS Network), and Parkfield (where a two-color laser system provides a proxy for continuous GPS measurements). An integral component of our study was the elucidation of the physical mechanisms by which such transients occur and propagate. We therefore initiated the development of multiple forward models, using two independent approaches. In the first, we explored the response to specified earthquake slip in viscoelastic models that incorporated failure criteria and the geometry of major faults in California. In the second approach, we examined the dynamical response of a complex rheological medium to the application of a far-field stress imposed by plate motions. The forward models were used both to gain insight into the range of strain transients to be expected under different assumed mechanical conditions and to develop representations for strain fields that allow GPS, borehole, and other strain data to be combined in a self-consistent, yet well-determined, manner. The models also provided a basis for hypothesis testing, by which data from a strain transient well characterized by GPS and borehole observations were utilized to distinguish among competing candidates for the causative physical mechanism and the governing physical characteristics. During the three years of this project, continued to a fourth year through a no-cost extension of the grant, we published 14 papers and presented or co-authored 37 papers at national scientific meetings.

Solomon, Sean C.↗
NASA NTRS2004

Crustal Deformation and the Seismic Cycle across the Kodiak Islands, Alaska

The Kodiak Islands are located approximately 130 to 250 km from the Alaska-Aleutian Trench where the Pacific plate is underthrusting the North American plate at a rate of about 57 mm/yr. The southern extent of the 1964 Prince William Sound (${M-w}$ = 9.2) earthquake rupture occurred offshore and beneath the eastern portion of the Kodiak Islands. Here we report GPS results (1993-2001) from northern Kodiak Island that span the transition between the 1964 uplift region along the eastern coast and the region of coseismic subsidence further inland. The horizontal velocity vectors range from 22.9 $\pm$ 2.2 mm/yr at N26.3$\deg$W $\pm$ 2.5$\deg$, about 150 km from the trench, to 5.9 $\pm$ 1.3 mm/yr at N65.9$\deg$W $\pm$ 6.6$\deg$, about 190 km from the trench. Near the northeastern coast of Kodiak the velocity vector above the shallow, locked main thrust zone is between the orientation of PCFC-NOAM plate motion (N22$/deg$W) and the trench-normal (N3O$\deg$W). Further west, our geodetic results suggest the accumulation of shear strain that will be released eventually as left-lateral motion on upper plate faults such as the Kodiak Island fault. These results are consistent with the hypothesis that the difference between the Pacific-North American plate motion and the orientation of the down going slab would lead to 4-8 mm/yr of left-lateral slip. Short-term geodetic uplift rates range from 2 - 14 mm/yr, with the maximum uplift located near the axis of maximum subsidence during the 1964 earthquake. We evaluated alternate interseismic models for Kodiak to test the importance of various mechanisms responsible for crustal deformation rates. These models are based on the plate interface slip history inferred from earlier modeling of coseismic and post-seismic geodetic results. The horizontal (trench perpendicular) and vertical deformation rates across Kodiak are consistent with a model that includes the viscoelastic response to : (1) a downgoing Pacific plate interface that is locked at shallow depths,(2) coseismic slip in the 1964 and (3) interseismic creep below the seismogenic zone. The change in orientation of the horizontal velocity vector occurs down-dip from the locked main thrust zone. In southern Kodiak, the coseismic slip in the 1964 earthquake was smaller than in the northern Kodiak region; yet, the horizontal, interseismic velocities as a function of distance from the trench are comparable to those in northern Kodiak. Based on the earthquake history prior to, and following the 1964 earthquake, we hypothesize that the plate interface in southern Kodiak slips in more frequent large earthquakes than in northern Kodiak.

Sauber, Jeanne↗
NASA NTRS2005

Crustal Deformation and the Seismic Cycle Across the Kodiak Islands, Alaska

The Kodiak Islands are located approximately 120 to 250 km from the Alaska-Aleutian Trench - and are within the southern extent of the 1964 Prince William Sound (M(sub w) = 9.2) earthquake rupture zone. Here we report new campaign GPS results (1993-2001) from northern Kodiak Island. The rate and orientation of the horizontal velocities, relative to a fixed North America, range from 25.3 plus or minus 1.4 mm/yr at N32.9 deg. W plus or minus 2.5 to 8.5 plus or minus 1.0 mm/yr at N59.7 deg. W plus or minus 6.5 deg. In addition to the northern Kodiak data, we analyzed data from three southern Kodiak Island stations. The inland stations from both the northern and southern networks indicate a counterclockwise rotation of the velocity vectors. These results are consistent with the hypothesis that the difference between the Pacific-North American plate motion and the orientation of the down going slab would lead to 4-8 mm/yr of left-lateral slip above the unlocked, down-dip portion of the main thrust zone. The northern and southern Kodiak geodetic data are consistent with a model that includes the viscoelastic response to (1) a downgoing Pacific plate interface that is locked at shallow depths, (2) local coseismic slip in the 1964 earthquake, and (3) interseismic creep down dip from the seismogenic zone. Based on the pre-1964 and post-1944 earthquake history, as well as the pattern of interseismic earthquakes across the plate boundary zone, we hypothesize that in southern Kodiak some strain is released in moderate to large earthquakes between the occurrences of great earthquakes like the 1964 event.

Sauber, Jeanne↗