Search NASA⌕ Search

SEARCH · Search NASA

Results for “Postseismic”

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 55 records · Page 3

Crustal deformation and earthquakes

The manner in which the Earth's surface deforms during the cycle of stress accumulation and release along major faults is investigated. In an investigation of the crustal deformation associated with a thin channel asthenosphere displacements are reduced from those computed for a half space asthenosphere. A previous finding by other workers that displacements are enhanced when flow is confined to a thin channel is based on several invalid approximations. The major predictions of the finite element model are that the near field postseismic displacements and strain rates are less than those for a half space asthenosphere and that the postseismic strain rates at intermediate distances are greater (in magnitude). The finite width of the asthenosphere ceases to have a significant impact on the crustal deformation pattern when its magnitude exceeds about three lithosphere thicknesses.

Cohen, S. C.↗

Earthquake deformation cycle in the Andean Back Arc, western Argentina

It is pointed out that the 1977 Caucete earthquake was an east-west compressive event which occurred at a depth of about 17 km in the Andean back arc approximately 450 km east of the Chile trench. The present study is primarily concerned with the mechanics of fault behavior prior to and after the 1977 earthquake, as evidenced by preseismic and postseismic vertical deformation. Simple fault models for each phase of deformation are presented. The models of preseismic slip, coseismic slip, and postseismic slip consist of uniform dislocations in an elastic half-space following the method of Chinnery (1961). The basically elastic models presented demonstrate that deformations in the vicinity of the 1977 Caucete, Argentina, earthquake are consistent with activity on some portion of the coseismic fault. This is evidence that these deformations are directly related to the processes of strain accumulation and release associated with the event.

Reilinger, R.↗

Interseismic deformation for normal fault earthquakes

Comparison between the coseismic vertical deformation associated with normal fault earthquakes and permanent deformation indicated by geologic structure indicates that large deformation must occur during the period between earthquakes. When sufficient geodetic and geologic information is available it is possible to estimate the spatial character of this interseismic deformation. A case in point is the 1983 Borah Peak, Idaho earthquake. This normal fault earthquake produced roughly 5 times as much basin subsidence as it did uplift of the adjacent mountain ranges. In contrast, geophysical and geological observations show that the basin is roughly as deep as the bounding range is high. A similar normal fault event in the Basin and Range, the 1959 Hebgen Lake earthquake, was also accompanied by substantially larger subsidence than uplift. In this case, postseismic geodetic measurements show broad regional uplift with a spatial pattern which is roughly consistent with that proposed for the Borah Peak event. Modeling suggests that postseismic viscoelastic relaxation and strain accumulation in an elastic lithosphere overlying a viscoelastic asthenosphere are possible physical mechanisms to generate interseismic uplift. These mechanisms may be the contemporary expression of those processes responsible for the high elevation of the Basin and Range Province.

Reilinger, Robert↗

Transient Creep of a Composite Lower Crust: Constitutive Theory - 1

A composite model is proposed to describe the time-dependent response of the Earth's lower crust. The motivation for such it model is twofold: First, new observations of widespread postseismic deformation indicate that the deep continental crust responds viscoelastically, having both long-and short-term decay times. Second, by any number of observationally based rationales, the lower crust is compositionally and structurally heterogeneous over many length scales. For heterogeneities that have much smaller characteristic lengths than the minimum deformation wavelength of interest, the aggregate rheology can be described by composite media theory. For wavelengths of the order of the thickness of the lower crust (approx. = 25-40 km) and larger, composite theory may be applied to heterogeneities that are smaller than about several hundred meters, or equivalent to the vertical extent of a thick lower crustal mylonitic shear zone. The composite media theory developed here is constructed using both Eshelhy-Mori-Tanaka theory for aligned generalized spheroidal inclusions and a generalized self-consistent method. The inclusions and matrix are considered to be Maxwellian viscoelastic: a rheology that is consistent with past homogeneous models of postseismic stress relaxation. The composite theory presented here introduces a transient response to a suddenly imposed stress field which does not appear in homogeneous Maxwell models. Analytic expressions for the amplitude and duration of the transient and for the effective long-and short-term viscosities of the composite are given which describe the sensitivity to inclusion concentration (phi), to shape, and to ratio of inclusion-to-matrix viscosity (R).

Ivins, Erik R.↗

Kinematics of the New Zealand Plate Boundary: Relative Motion by GPS Across Networks of 1000 km and 50 km Spacing

From the 1995-97 data we are able to estimate interseismic velocities - the original aim of the project - though corrections must first be made for the coseismic effects of the Cass earthquake and postseismic effects of the Arthur's Pass event. In the following, we present four sets of results: (1) Interseismic velocities, (2) Arthur's Pass earthquake - coseismic, (3) Cass earthquake, and (4) Arthur's Pass earthquake postseismic. Work on all these fronts is still incomplete, but is being actively continued with funding from the NZ Foundation for Research, Science and Technology (FRST).

Scholz, Chris↗

Satellite Geodetic Constraints On Earthquake Processes: Implications of the 1999 Turkish Earthquakes for Fault Mechanics and Seismic Hazards on the San Andreas Fault

Our principal activities during the initial phase of this project include: 1) Continued monitoring of postseismic deformation for the 1999 Izmit and Duzce, Turkey earthquakes from repeated GPS survey measurements and expansion of the Marmara Continuous GPS Network (MAGNET), 2) Establishing three North Anatolian fault crossing profiles (10 sitedprofile) at locations that experienced major surface-fault earthquakes at different times in the past to examine strain accumulation as a function of time in the earthquake cycle (2004), 3) Repeat observations of selected sites in the fault-crossing profiles (2005), 4) Repeat surveys of the Marmara GPS network to continue to monitor postseismic deformation, 5) Refining block models for the Marmara Sea seismic gap area to better understand earthquake hazards in the Greater Istanbul area, 6) Continuing development of models for afterslip and distributed viscoelastic deformation for the earthquake cycle. We are keeping close contact with MIT colleagues (Brad Hager, and Eric Hetland) who are developing models for S. California and for the earthquake cycle in general (Hetland, 2006). In addition, our Turkish partners at the Marmara Research Center have undertaken repeat, micro-gravity measurements at the MAGNET sites and have provided us estimates of gravity change during the period 2003 - 2005.

Reilinger, Robert↗

Uplift and Seismicity Driven by Magmatic Inflation at Sierra Negra Volcano, Galápagos Islands

Abstract Although episodes of surface uplift and elevated seismicity precede many volcanic eruptions, their temporal evolution is often complex, and apparently in contradiction to simple trends predicted by mechanical deformation models. Here, we use continuous global positioning system and seismic data recorded at Sierra Negra volcano, Galápagos Islands, to show how the edifice responded to stress changes driven by magma accumulation in a shallow sill. The rate of uplift varied during the 13 years and 6.5 m of inflation before the 2018 eruption. The number of earthquakes per unit of uplift increased exponentially with total uplift as the differential stress increased. Accordingly, the temporal seismicity rate varied in time as a function of both the total uplift and the uplift rate. The Gutenberg‐Richter b ‐value decreased as a function of total uplift. In the final six months before the eruption, a sequence of large ( M > 4) earthquakes regulated the state of stress on the fault, each being followed by 2–3 days of postseismic quiescence, and retarding the increase in seismicity rate. These earthquakes did not affect the overall uplift rate. Subsidence of 8.5 m accompanied the 2‐month eruption. On resumption of uplift, the number of earthquakes per unit of uplift was very low, and the b ‐value high, reflecting the relaxed stress state of the fault system. These observations show that crustal deformation becomes increasingly brittle at higher stress states, and supports theoretical models based on elastic‐brittle mechanics. They suggest that joint interpretation of deformation and seismicity is key for forecasting future eruptions in similar volcanic settings.

Bell, Andrew Forbes↗

The Role of Dislocations in the Anelasticity of the Upper Mantle

Dislocation‐based dissipation mechanisms potentially control the viscoelastic response of Earth's upper mantle across a variety of geodynamic contexts, including glacial isostatic adjustment, postseismic creep, and seismic‐wave attenuation. However, there is no consensus on which dislocation‐based, microphysical process controls the viscoelastic behavior of the upper mantle. Although both intergranular (plastic anisotropy) and intragranular (backstress) mechanisms have been proposed, there is currently insufficient laboratory data to discriminate between those mechanisms. Here, we present the results of forced‐oscillation experiments in a deformation‐DIA apparatus at confining pressures of 3–7 GPa and temperatures of 298–1370 K. Our experiments tested the viscoelastic response of polycrystalline olivine—the main constituent of the upper mantle—at stress amplitudes from 70 to 2,800 MPa. Mechanical data are complemented by microstructural analyses of grain size, crystallographic preferred orientation, and dislocation density. We observe amplitude‐ and frequency‐dependent attenuation and modulus relaxation and find that numerical solutions of the backstress model match our results well. Therefore, we argue that interactions among dislocations, rather than intergranular processes (e.g., plastic anisotropy or grain boundary sliding), control the viscoelastic behavior of polycrystalline olivine in our experiments. In addition, we present a linearized version of the constitutive equations of the backstress model and extrapolate it to conditions typical of seismic‐wave propagation in the upper mantle. Our extrapolation demonstrates that the backstress model can explain the magnitude of seismic‐wave attenuation in the upper mantle, although some modification is required to explain the weak frequency dependence of attenuation observed in nature and in previous experimental work.

Hein, Diede [Univ. of Minnesota, Minneapolis, MN (↗

Probing the evolution of fault properties during the seismic cycle with deep learning

We use seismic waves that pass through the hypocentral region of the 2016 M6.5 Norcia earthquake together with Deep Learning (DL) to distinguish between foreshocks, aftershocks and time-to-failure (TTF). Binary and N-class models defined by TTF correctly identify seismograms in test with > 90% accuracy. We use raw seismic records as input to a 7 layer CNN model to perform the classification. Here we show that DL models successfully distinguish seismic waves pre/post mainshock in accord with lab and theoretical expectations of progressive changes in crack density prior to abrupt change at failure and gradual postseismic recovery. Performance is lower for band-pass filtered seismograms (below 10 Hz) suggesting that DL models learn from the evolution of subtle changes in elastic wave attenuation. Tests to verify that our results indeed provide a proxy for fault properties included DL models trained with the wrong mainshock time and those using seismic waves far from the Norcia mainshock; both show degraded performance. Our results demonstrate that DL models have the potential to track the evolution of fault zone properties during the seismic cycle. If this result is generalizable it could improve earthquake early warning and seismic hazard analysis.

58 GEOSCIENCES↗

Dislocation theory of steady and transient creep of crystalline solids: Predictions for olivine

In applications critical to the geological, materials, and engineering sciences, deformation occurs at strain rates too small to be accessible experimentally. Instead, extrapolations of empirical relationships are used, leading to epistemic uncertainties in predictions. To address these problems, we construct a theory of the fundamental processes affecting dislocations: storage and recovery. We then validate our theory for olivine deformation. This model explains the empirical relationships among strain rate, applied stress, and dislocation density in disparate laboratory regimes. It predicts the previously unexplained dependence of dislocation density on applied stress in olivine. The predictions of our model for Earth conditions differ from extrapolated empirical relationships. For example, it predicts rapid, transient deformation in the upper mantle, consistent with recent measurements of postseismic creep.

36 MATERIALS SCIENCE↗

Crustal Strain Rates in the Western United States and Their Relationship with Earthquake Rates

Abstract We present a suite of strain rate models for the western United States based on geologic and geodetic data. The geologic data consist of Quaternary fault-slip rates and the geodetic data consists of a new compilation of Global Positioning System (GPS) velocities derived from continuous, semicontinuous, and campaign measurements. We remove postseismic deformation from the GPS time series in order for our geodetic strain rate model to best capture the interseismic strain accumulation rate. We present models based on either geologic or geodetic data, but also create a hybrid model. Although there are some differences between the models, the large-scale features are the same, with the noticeable exception for the Pacific Northwest where interseismic strain is naturally more distributed than the long-term strain release. We also present a map of earthquake rate densities based on mainshocks, and the result has similar spatial features similar to the strain rate models (at least in the southwestern United States). We perform a general correlation analysis between strain rate and seismicity rate (south of Cascadia) and find a change in linearity between seismicity and strain rates from slow to faster deforming areas with seismicity rates relatively lower for the latter. The extent of that change depends a bit on assumptions made on the declustering and completeness of the catalog, but the finding of a change in slope is robust across the different strain rate models. Linearity for all areas is only expected when Gutenberg–Richter parameters and parameters involved in the conversion from strain to moment rate are uniform across the study area. We discuss these qualifications, but find no single satisfactory explanation for our observation. Moreover, when considering a rather short time and space, theoretical considerations of sampling from a power-law distribution actually predict there to be a power law instead of a linear relationship, generally consistent with our observation.

Geochemistry & Geophysics↗

Mechanism of underthrusting in southwest Japan: A model of convergent plate interactions

An elastic rebound mechanism consistent with underthrusting at the time of the magnitude 8.2 Nankaido earthquake of 21 December, 1946 accounts for a reversal in sense between seismic and preseismic changes in elevation throughout a large portion of southwest Japan. It is shown that the seismic movements during the earthquake were generated by slip on a complex thrust fault that is inferred to intersect the surface near the base of the continental slope. The preseismic deformation is explained by strain accumulation equivalent to a virtual dislocation with the same orientation as the fault surface but a sense opposite to that of the real dislocation that occurred at the time of the earthquake. Adjustments by slip along the fault surface and extensions of that surface can account for postseismic movements.

Fitch, T. J.↗

The behavior of a convergent plate boundary - Crustal deformation in the South Kanto district, Japan

The northwesternmost part of the Sagami trough, a part of the Philippine Sea-Eurasian plate boundary, was ruptured during the great South Kanto earthquake in 1923. Very extensive and frequent geodetic measurements of crustal deformation have been made in the South Kanto district since the 1890's, and these constitute the most complete data set on crustal movements in the world. These data were reanalyzed and interpreted and according to our interpretation indicate the following sequence of events. The coseismic movements were due to oblique thrust and right lateral slip of about 8 m on a fault outcropping at the base of the Sagami trough. This was followed by postseismic deformation resulting from reversed afterslip of 20-60 cm that occurred at an exponentially decaying rate in time. The interseismic deformation is produced by steady subduction at a rate of about 1.8 cm/yr. During subduction the top 10-15 km of the plate boundary is apparently locked, while deeper parts slip aseismically at an irregular rate. No significant precursory deformation was observed. The recurrence time for 1923 type earthquakes is 200-300 years. The Boso and Miura peninsulas are broken into a series of fault-bound blocks that move semi-independently of the surrounding region. The subduction zone itself, where it is exposed on land, is shown to be a wide zone encompassing several faults that are active at different times.

Scholz, C. H.↗

The viscoelastic stiffness model of seismicity

A viscoelastic stiffness model of seismicity is developed by introducing a viscoelastic element into the stiffness model for fault dynamics. The introduction of this element permits modeling of transient anelastic deformations in response to stress loading and relaxation and provides a mechanism for partial stress recovery following an earthquake. As a consequence, several phenomena not present in elastic stiffness theory emerge. These include postseismic creep, foreshocks, and aftershocks. Numerical simulations of fault motion also reveal episodes of stable sliding, tertiary creep preceeding earthquakes, and long-term aseismic creep.

Cohen, S. C.↗

Theory of time-dependent rupture in the Earth

Fracture mechanics is used to develop a theory of earthquake mechanism which includes the phenomenon of subcritical crack growth. The following phenomena are predicted: slow earthquakes, multiple events, delayed multiple events (doublets), postseismic rupture growth and afterslip, foreshocks, and aftershocks. The theory predicts a nucleation stage prior to an earthquake, and suggests a physical mechanism by which one earthquake may 'trigger' another.

Das, S.↗

Stress corrosion cracking in the earth

Two fundamental concepts of fracture mechanics are used to develop a theory of the earthquake mechanism which specifically predicts observed time-dependent rupture phenomena such as slow earthquakes, postseismic rupture growth and afterslip, multiple events, foreshocks, and aftershocks. The theory also predicts that there must be a nucleation stage prior to an earthquake, and suggests a physical mechanism by which one earthquake may trigger another. Investigations show that all earthquakes must be preceded by a quasi-static slip over a portion of the rupture surfaces, although it may be difficult to detect in practice, and a study of delayed multiple events characterizes the strength of some barriers in the earth as having a stress corrosion index of about 24.

Das, S.↗

Vertical movements following a dip-slip earthquake

The results of finite element (FEM) computations based on elastic dislocation theory are presented. It is shown that the ratio of fault depth to the lithospheric thickness is a significant factor in the sign of the rebound. Subsidence occurs if the fault penetrates less than approximately 0.68 of the way through the lithosphere. An uplift occurs if the fault cuts deeper. The FEM employed consisted of a 546 node, 500 element grid and a quasi-plane strain algorithm. The situation modeled consisted of a 30 deg dip fault in an initially planar elastic lithosphere overlying a Maxwell viscoelastic asthenosphere. It is suggested that postseismic rebound may be measurable over a several year interval, thus providing data on lithosphere thicknesses.

Melosh, H. J.↗

Strain accumulation and surface deformation along the San Andreas, California

Stressing and rupture of a locked zone adjacent to a creeping fault segment was studied with special reference to strength heterogeneity depthwise and along-strike. The resulting precursory temporal and spatial variations of surface strain rate profiles were compared to geodetic measurements on the San Andreas fault in central California. Crustal deformation in great California earthquake cycles was also studied with special reference to the temporal decay of strain rate observed since the 1957 and 1906 great earthquakes, and comtemporary surface strain rate and velocity profiles at several locations along the San Andreas. The effect of viscoelastic response in the deep aseismic shear zone on the surface deformation behavior was examined. Work was begun on a fundamental reformulation of the crustal deformation problem focusing on the crustal deformation process affected by deep aseismic slip as the slip zone progresses toward an instability and as deep seismic slip continues postseismically, the 3-D nature of the problem due to geometry and material heterogeneity, and the time-dependent source coming from the lithosphere/astenospheric coupling process.

Li, Victor C.↗