Search NASA⌕ Search

SEARCH · Search NASA

Results for “Earthquake Prediction”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Use of QuakeSim and UAVSAR for Earthquake Damage Mitigation and Response

Spaceborne, airborne, and modeling and simulation techniques are being applied to earthquake risk assessment and response for mitigation from this natural disaster. QuakeSim is a web-based portal for modeling interseismic strain accumulation using paleoseismic and crustal deformation data. The models are used for understanding strain accumulation and release from earthquakes as well as stress transfer to neighboring faults. Simulations of the fault system can be used for understanding the likelihood and patterns of earthquakes as well as the likelihood of large aftershocks from events. UAVSAR is an airborne L-band InSAR system for collecting crustal deformation data. QuakeSim, UAVSAR, and DESDynI (following launch) can be used for monitoring earthquakes, the associated rupture and damage, and postseismic motions for prediction of aftershock locations.

fault systems↗
NASA NTRS1993

Earthquake slip vectors and estimates of present-day plate motions

Two alternative models for present-day global plate motions are derived from subsets of the NUVEL-1 data in order to investigate the degree to which earthquake slip vectors affect the NUVEL-1 model and to provide estimates of present-day plate velocities that are independent of earthquake slip vectors. The data set used to derive the first model excludes subduction zone slip vectors. The primary purpose of this model is to demonstrate that the 240 subduction zone slip vectors in the NUVEL-1 data set do not greatly affect the plate velocities predicted by NUVEL-1. A data set that excludes all of the 724 earthquake slip vectors used to derive NUVEL-1 is used to derive the second model. This model is suitable as a reference model for kinematic studies that require plate velocity estimates unaffected by earthquake slip vectors. The slip-dependent slip vector bias along transform faults is investigated using the second model, and evidence is sought for biases in slip directions along spreading centers.

Demets, Charles↗
NASA NTRS1980

Proposed method for determining the thickness of glass in solar collector panels

An analytical method was developed for determining the minimum thickness for simply supported, rectangular glass plates subjected to uniform normal pressure environmental loads such as wind, earthquake, snow, and deadweight. The method consists of comparing an analytical prediction of the stress in the glass panel to a glass breakage stress determined from fracture mechanics considerations. Based on extensive analysis using the nonlinear finite element structural analysis program ARGUS, design curves for the structural analysis of simply supported rectangular plates were developed. These curves yield the center deflection, center stress and corner stress as a function of a dimensionless parameter describing the load intensity. A method of estimating the glass breakage stress as a function of a specified failure rate, degree of glass temper, design life, load duration time, and panel size is also presented.

Moore, D. M.↗
NASA NTRS1984

The distribution of earthquakes with depth and stress in subducting slabs

The global variation of Benioff zone seismicity with depth and the orientation of stress axes of deep and intermediate earthquakes is explained using numerical models of subducting slabs. Models that match the seismicity and stress require a barrier to flow at the 670 km seismic discontinuity. The barrier may be a viscosity increase of at least an order of magnitude or a chemical discontinuity. Instantaneous flow is subparallel to the slabs for models with a viscosity increase but contorted for models with a chemical barrier. Log N (number of earthquakes) decreases linearly to 250-300 km depth and increases thereafter. Stress magnitude in the models shows the same pattern, in accord with experiments showing N proportional to e(k-sigma), with k a constant and sigma the stress magnitude. The models predict downdip compression in the slabs at depths below 300-400 km, as observed for earthquake stress axes.

Vassiliou, M. S.↗
NASA NTRS1977

Satellite Relaying of Geophysical Data

Data Collection Platforms (DCPs) for transmitting surface data to an orbiting satellite for relaying to a central data distribution center are being used in a number of geophysical applications. "Off-the-shelf" DCP's, transmitting through Landsat or GOES satellites, are fully capable of relaying data from low-data-rate instruments, such as tiltmeters or tide gauges. In cooperation with the Lamont-Doherty Geological Observatory, Goddard has successfully installed DCP systems on a tide gauge and tiltmeter array on Anegada, British Virgin Islands. Because of the high-data-rate requirements, a practical relay system capable of handling seismic information is not yet available. Such a system could become the basis of an operational hazard prediction system for reducing losses due to major natural catastrophies such as earthquakes, volcanic eruptions, landslides or tsunamis.

Allenby, R. J.↗
NASA NTRS1991

A rheologically layered three-dimensional model of the San Andreas fault in central and southern California

The effects of rheological parameters and the fault slip distribution on the horizontal and vertical deformation in the vicinity of the fault are investigated using 3D kinematic finite element models of the San Andreas fault in central and southern California. It is shown that fault models with different rheological stratification schemes and slip distributions predict characteristic deformation patterns. Models that do not include aseismic slip below the fault locking depth predict deformation patterns that are strongly dependent on time since the last earthquake, while models that incorporate the aseismic slip below the locking depth depend on time to a significantly lesser degree.

Williams, Charles A.↗
NASA NTRS1991

Models of recurrent strike-slip earthquake cycles and the state of crustal stress

Numerical models of the strike-slip earthquake cycle, assuming a viscoelastic asthenosphere coupling model, are examined. The time-dependent simulations incorporate a stress-driven fault, which leads to tectonic stress fields and earthquake recurrence histories that are mutually consistent. Single-fault simulations with constant far-field plate motion lead to a nearly periodic earthquake cycle and a distinctive spatial distribution of crustal shear stress. The predicted stress distribution includes a local minimum in stress at depths less than typical seismogenic depths. The width of this stress 'trough' depends on the magnitude of crustal stress relative to asthenospheric drag stresses. The models further predict a local near-fault stress maximum at greater depths, sustained by the cyclic transfer of strain from the elastic crust to the ductile asthenosphere. Models incorporating both low-stress and high-stress fault strength assumptions are examined, under Newtonian and non-Newtonian rheology assumptions. Model results suggest a preference for low-stress (a shear stress level of about 10 MPa) fault models, in agreement with previous estimates based on heat flow measurements and other stress indicators.

Lyzenga, Gregory A.↗
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↗
NASA NTRS1995

Coseismic fault slip associated with the 1992 M(sub w) 6.1 Joshua Tree, California, earthquake: Implications for the Joshua Tree-Landers earthquake sequence

Coseismic surface deformation associated with the M(sub w) 6.1, April 23, 1992, Joshua Tree earthquake is well represented by estimates of geodetic monument displacements at 20 locations independently derived from Global Positioning System and trilateration measurements. The rms signal to noise ratio for these inferred displacements is 1.8 with near-fault displacement estimates exceeding 40 mm. In order to determine the long-wavelength distribution of slip over the plane of rupture, a Tikhonov regularization operator is applied to these estimates which minimizes stress variability subject to purely right-lateral slip and zero surface slip constraints. The resulting slip distribution yields a geodetic moment estimate of 1.7 x 10(exp 18) N m with corresponding maximum slip around 0.8 m and compares well with independent and complementary information including seismic moment and source time function estimates and main shock and aftershock locations. From empirical Green's functions analyses, a rupture duration of 5 s is obtained which implies a rupture radius of 6-8 km. Most of the inferred slip lies to the north of the hypocenter, consistent with northward rupture propagation. Stress drop estimates are in the range of 2-4 MPa. In addition, predicted Coulomb stress increases correlate remarkably well with the distribution of aftershock hypocenters; most of the aftershocks occur in areas for which the mainshock rupture produced stress increases larger than about 0.1 MPa. In contrast, predicted stress changes are near zero at the hypocenter of the M(sub w) 7.3, June 28, 1992, Landers earthquake which nucleated about 20 km beyond the northernmost edge of the Joshua Tree rupture. Based on aftershock migrations and the predicted static stress field, we speculate that redistribution of Joshua Tree-induced stress perturbations played a role in the spatio-temporal development of the earth sequence culminating in the Landers event.

Bennett, Richard A.↗
NASA NTRS1993

A model for the motion of the Philippine Sea plate consistent with NUVEL-1 and geological data

We investigate angular velocity vectors of the Philippine Sea (PH) plate relative to the adjacent major plates, Eurasia (EU) and Pacific (PA), and the smaller Caroline (CR) plate. Earthquake slip vector data along the Philippine Sea plate are inverted, subject to the constraint that EU-PA motion equals that predicted by the global relative plate model NUVEL-1. The resulting solution fails to satisfy geological constraints along the Caroline-Pacific boundary: convergence along the Mussau Trench and divergence along the Sorol Trough. We then seek solutions satisfying both the CR-PA boundary conditions and the Philippine Sea slip vector data, by adjusting the PA-PH and EU-PH best fitting poles within their error ellipses. We also consider northern Honshu to be part of the North American plate and impose the constraint that the Philippine Sea plate subducts beneath northern Honshu along the Sagmi Trough in a NNW-NW direction. Of the solutions satisfying these conditions, we select the best EU-PH as 48.2 deg N, 157.0 deg E, 1.09 deg/my, corresponding to a pole far from Japan and south of Kamchatka, and PA-PH, 1.2 deg N, 134.2 deg E, 1.00 deg/my. Predicted NA-PH and EU-PH convergence rates in central Honshu are consistent with estimated seismic slip rates. Previous estimates of the EU-PH pole close to central Honshu are inconsistent with extension within the Bonin backarc implied by earthquake slip vectors and NNW-NW convergence of the Bonin forearc at the Sagami Trough.

Seno, Tetsuzo↗
NASA NTRS2013

PHyTIR - A Prototype Thermal Infrared Radiometer

This paper describes the PHyTIR (Prototype HyspIRI Thermal Infrared Radiometer) instrument, which is the engineering model for the proposed HyspIRI (Hyperspectral Infrared Imager) earth observing instrument. The HyspIRI mission would be comprised of the HyspIRI TIR (Thermal Infrared Imager), and a VSWIR (Visible Short-Wave Infra-Red Imaging Spectrometer). Both instruments would be used to address key science questions related to the earth's carbon cycle, ecosystems, climate, and solid earth properties. Data gathering of volcanic activities, earthquakes, wildfires, water use and availability, urbanization, and land surface compositions and changes, would aid the predictions and evaluations of such events and the impact they create. Even though the proposed technology for the HyspIRI imager is mature, the PHyTIR prototype is needed to advance the technology levels for several of the instrument's key components, and to reduce risks, in particular to validate 1) the higher sensitivity, spatial resolution, and higher throughput required for this focal plane array, 2) the pointing accuracy, 2) the characteristics of several spectral channels, and 4) the use of ambient temperature optics. The PHyTIR telescope consists of the focal plane assembly that is housed within a cold housing located inside a vacuum enclosure; all mounted to a bulkhead, and an optical train that consists of 3 powered mirrors; extending to both sides of the bulkhead. A yoke connects the telescope to a scan mirror. The rotating mirror enables to scan- a large track on the ground. This structure is supported by kinematic mounts, linking the telescope assembly to a base plate that would also become the spacecraft interface for HyspIRI. The focal plane's cooling units are also mounted to the base plate, as is an overall enclosure that has two viewing ports with large exterior baffles, shielding the focal plane from incoming stray light. PHyTIR's electronics is distributed inside and near the vacuum enclosure, and in a nearby rack. The data acquisition technique would be to take measurements over a 51deg wide swath in the cross spacecraft velocity direction, which is brought into view through the rotating scan mirror. A landscape mosaic thus can be assembled by overlaying rows of measurements. The paper briefly outlines the proposed HyspIRI mission and its data acquisition technique; it then describes the prototype PHyTIR instrument.

Space Instrument↗
NASA NTRS2020

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↗
NASA NTRS1987

Crustal deformation in Great California Earthquake cycles

A model in which coupling is described approximately through a generalized Elsasser model is proposed for computation of the periodic crustal deformation associated with repeated strike-slip earthquakes. The model is found to provide a more realistic physical description of tectonic loading than do simpler kinematic models. Parameters are chosen to model the 1857 and 1906 San Andreas ruptures, and predictions are found to be consistent with data on variations of contemporary surface strain and displacement rates as a function of distance from the 1857 and 1906 rupture traces. Results indicate that the asthenosphere appropriate to describe crustal deformation on the earthquake cycle time scale lies in the lower crust and perhaps the crust-mantle transition zone.

Li, Victor C.↗
NASA NTRS1982

The earthquake cycle in subduction zones

A simplified model of a subduction zone is presented, which incorporates the mechanical asymmetry induced by the subducted slab to anchor the subducting plate during post-seismic rebound and thus throw most of the coseismic stream release into the overthrust plate. The model predicts that the trench moves with respect to the deep mantle toward the subducting plate at a velocity equal to one-half of the convergence rate. A strong extensional pulse is propagated into the overthrust plate shortly after the earthquake, and although this extension changes into compression before the next earthquake in the cycle, the period of strong extension following the earthquake may be responsible for extensional tectonic features in the back-arc region.

Melosh, H. J.↗
NASA NTRS2019

Completeness Index for Earthquake-Induced Landslide Inventories

Understanding the global relation between Earthquake-Induced Landslides (EQILs) and the factors contributing to their initiation is still an open topic within the geomorphological community. Accessing EQIL inventories and analyzing them concerning their potential causes is the key to explore such relation. However, each of the existing EQIL inventories has its level of completeness and associated uncertainty which makes any unified relationship which is challenging to obtain. So far, the completeness of EQIL inventories has never been clearly defined. As a result, it has never been accounted for in global EQIL predictive models. In this technical note, we propose a simple definition for the completeness of EQIL inventories. We analyze 30 digital EQIL for 21 earthquakes and develop a semi-quantitative method to estimate the completeness level, which we refer to as Completeness Index (CI). The CI results from a combination of topographic factors, ground shaking parameters, and a measure derived from landslide size statistics. The proposed CI consists of Low, Moderate, and High completeness classes which can be used to evaluate any landslide inventory. We made the whole procedure to compute the CI accessible in an ArcGIS toolbox together with a test dataset.

Hakan Tanyas↗
NASA NTRS1991

A synthetic seismicity model for the Middle America Trench

A novel iterative technique, based on the concept of fault segmentation and computed using 2D static dislocation theory, for building models of seismicity and fault interaction which are physically acceptable and geometrically and kinematically correct, is presented. The technique is applied in two steps to seismicity observed at the Middle America Trench. The first constructs generic models which randomly draw segment strengths and lengths from a 2D probability distribution. The second constructs predictive models in which segment lengths and strengths are adjusted to mimic the actual geography and timing of large historical earthquakes. Both types of models reproduce the statistics of seismicity over five units of magnitude and duplicate other aspects including foreshock and aftershock sequences, migration of foci, and the capacity to produce both characteristic and noncharacteristic earthquakes. Over a period of about 150 yr the complex interaction of fault segments and the nonlinear failure conditions conspire to transform an apparently deterministic model into a chaotic one.

Ward, Steven N.↗
NASA NTRS1986

Viscoelastic deformation near active plate boundaries

Model deformations near the active plate boundaries of Western North America using space-based geodetic measurements as constraints are discussed. The first six months of this project were spent gaining familarity with space-based measurements, accessing the Crustal Dynamics Data Information Computer, and building time independent deformation models. The initial goal was to see how well the simplest elastic models can reproduce very long base interferometry (VLBI) baseline data. From the Crustal Dynamics Data Information Service, a total of 18 VLBI baselines are available which have been surveyed on four or more occasions. These data were fed into weighted and unweighted inversions to obtain baseline closure rates. Four of the better quality lines are illustrated. The deformation model assumes that the observed baseline rates result from a combination of rigid plate tectonic motions plus a component resulting from elastic strain build up due to a failure of the plate boundary to slip at the full plate tectonic rate. The elastic deformation resulting from the locked plate boundary is meant to portray interseismic strain accumulation. During and shortly after a large interplate earthquake, these strains are largely released, and points near the fault which were previously retarded suddenly catch up to the positions predicted by rigid plate models. Researchers judge the quality of fit by the sum squares of weighted residuals, termed total variance. The observed baseline closures have a total variance of 99 (cm/y)squared. When the RM2 velocities are assumed to model the data, the total variance increases to 154 (cm/y)squared.

Ward, S. N.↗
NASA NTRS1993

An analysis of nearfield normal mode amplitude anomalies of the Landers earthquake

The 1992 Landers earthquake (M(sub w) = 7.3) occurred in the middle of the TERRAscope network. Long-period Rayleigh waves recorded at the TERRAscope stations (delta less than or = 3 deg) after traveling around the Earth show large amplitude anomalies, one order of magnitude larger than spherical Earth predictions up to a period of about 600 s. The ground motions over the epicentral region at and after the arrival of R4-5 are in phase at all stations. These observations are inconsistent with the nearly vertical strike slip mechanism of the Landers earthquake. Synthetic seismograms for a rotating, elliptic, and laterally heterogeneous Earth model calculated by the variational method agree well with the observed waveforms. Calculations for various 3D Earth models demonstrate that the amplitudes are very sensitive to the large scale aspherical structure in the crust and the mantle. The anomalies for modes shorter than 300 s period can be explained by lateral heterogeneity shallower than the upper mantle. Rotation of the Earth and lower mantle heterogeneity are required to explain mode amplitudes at longer periods. Current whole mantle seismic tomographic models can fully explain the observed amplitudes longer than 300 s. To assess the effect of the high order lateral heterogeneity in the mantle more precise estimate of the crustal correction is required.

Watada, Shingo↗