Report of the panel on plate motion and deformation, section 2

Given here is a panel report on the goals and objectives, requirements and recommendations for the investigation of plate motion and deformation. The goals are to refine our knowledge of plate motions, study regional and local deformation, and contribute to the solution of important societal problems. The requirements include basic space-positioning measurements, the use of global and regional data sets obtained with space-based techniques, topographic and geoid data to help characterize the internal processes that shape the planet, gravity data to study the density structure at depth and help determine the driving mechanisms for plate tectonics, and satellite images to map lithology, structure and morphology. The most important recommendation of the panel is for the implementation of a world-wide space-geodetic fiducial network to provide a systematic and uniform measure of global strain

Crustal deformation at the Sumatran Subduction Zone revealed by coral rings

Analyses of coral rings grown in the interval 1970–1997 reveal a geographically distinct pattern of interseismic uplift off Sumatra's western coast. At distances less than 110 km from the Sumatran trench, coral reefs are submerging as fast as 5 mm/y. At 130 and 180 km distance from the trench, they are emerging at similar rates. We suggest that a locked, or partially locked patch, located above 30 km depth on the upper surface of the subducting oceanic plate, generates this pattern

Active megadetachment beneath the western United States

Geodetic data, interpreted in light of seismic imaging, seismicity, xenolith studies, and the late Quaternary geologic history of the northern Great Basin, suggest that a subcontinental-scale extensional detachment is localized near the Moho. To first order, seismic yielding in the upper crust at any given latitude in this region occurs via an M7 earthquake every 100 years. Here we develop the hypothesis that since 1996, the region has undergone a cycle of strain accumulation and release similar to “slow slip events” observed on subduction megathrusts, but yielding occurred on a subhorizontal surface 5–10 times larger in the slip direction, and at temperatures >800°C. Net slip was variable, ranging from 5 to 10 mm over most of the region. Strain energy with moment magnitude equivalent to an M7 earthquake was released along this “megadetachment,” primarily between 2000.0 and 2005.5. Slip initiated in late 1998 to mid-1999 in northeastern Nevada and is best expressed in late 2003 during a magma injection event at Moho depth beneath the Sierra Nevada, accompanied by more rapid eastward relative displacement across the entire region. The event ended in the east at 2004.0 and in the remainder of the network at about 2005.5. Strain energy thus appears to have been transmitted from the Cordilleran interior toward the plate boundary, from high gravitational potential to low, via yielding on the megadetachment. The size and kinematic function of the proposed structure, in light of various proxies for lithospheric thickness, imply that the subcrustal lithosphere beneath Nevada is a strong, thin plate, even though it resides in a high heat flow tectonic regime. A strong lowermost crust and upper mantle is consistent with patterns of postseismic relaxation in the southern Great Basin, deformation microstructures and low water content in dunite xenoliths in young lavas in central Nevada, and high-temperature microstructures in analog surface exposures of deformed lower crust. Large-scale decoupling between crust and upper mantle is consistent with the broad distribution of strain in the upper crust versus the more localized distribution in the subcrustal lithosphere, as inferred by such proxies as low P wave velocity and mafic magmatism

Mode of Strain Release Associated with Major Earthquakes in Japan

It is now widely accepted that most shallow earthquakes are caused by a sudden release of tectonic stress, and therefore strain, in the form of faulting (56, Chap. 14). The elastic rebound theory explains remarkably well the overall nature of major earthquakes. However, it is now evident that the strain release by actual earthquakes takes place in a variety of ways. For example, some earthquakes arc reported to have been preceded by significant premonitory crustal deformations with or without foreshocks. Most major earthquakes are followed by aftershocks and, in some cases, by creep-like postseismic deformations. This variety obviously results from the difference in the nature of the crust and the stress field in the respective epicentral area. Thus, detailed study on the mode of the strain release for individual earthquakes is important for understanding the physics of earthquakes and for predicting the mode of occurrence of future major earthquakes. The latter problem is of course very important for designing a practical method of earthquake prediction. The crustal deformation represents the low frequency spectrum of the strain release, and seismic waves represent the high frequency spectrum. To understand better the mode of the strain release, it is important to study it over a wide frequency range, from geodetic to seismic. From this point of view, this paper discusses the mode of strain release in several major earthquakes in Japan for which both geodetic and seismic data are available, and we will focus our attention on the relation between the geodetic slip and the seismic slip. Figure 1 shows the locations of the earthquakes discussed in this paper. We do not intend to make an exhaustive review of the literature on the subject that can be referred to (19, 21, 36, 39, 50, 63, 71)

An elastic plate model for interseismic deformation in subduction zones

Geodetic observations of interseismic surface deformation in the vicinity of subduction zones are frequently interpreted using simple kinematic elastic dislocation models (EDM). In this theoretical study, we develop a kinematic EDM that simulates plate subduction over the interseismic period (the elastic subducting plate model (ESPM)) having only 2 more degrees of freedom than the well-established back slip model (BSM): an elastic plate thickness and the fraction of flexural stresses due to bending at the trench that are released continuously. Unlike the BSM, in which steady state deformation in both plates is assumed to be negligible, the ESPM includes deformation in the subducting and overriding plates (owing to plate thickness), while still preserving the correct sense of convergence velocity between the subducting and overriding plates, as well as zero net steady state vertical offset between the two plates when integrated over many seismic cycles. The ESPM links elastic plate flexure processes to interseismic deformation and helps clarify under what conditions the BSM is appropriate for fitting interseismic geodetic data at convergent margins. We show that the ESPM is identical to the BSM in the limiting case of zero plate thickness, thereby providing an alternative motivation for the BSM. The ESPM also provides a consistent convention for applying the BSM to any megathrust interface geometry. Even in the case of nonnegligible plate thickness, the deformation field predicted by the ESPM reduces to that of the BSM if stresses related to plate flexure at the trench are released either continuously and completely at shallow depths during the interseismic period or deep in the subduction zone (below ∼100 km). However, if at least a portion of these stresses are not continuously released in the shallow portion of the subduction zone (via seismic or aseismic events), then the predicted surface velocities of these two models can differ significantly at horizontal distances from the trench equivalent to a few times the effective interseismic locking depth

Inference of Multiple Earthquake-Cycle Relaxation Timescales from Irregular Geodetic Sampling of Interseismic Deformation

Characterizing surface deformation throughout a full earthquake cycle is a challenge due to the lack of high‐resolution geodetic observations of duration comparable to that of characteristic earthquake recurrence intervals (250–10,000 years). Here we approach this problem by comparing long‐term geologic slip rates with geodetically derived fault slip rates by sampling only a short fraction (0.001%–0.1%) of a complete earthquake cycle along 15 continental strike‐slip faults. Geodetic observations provide snapshots of surface deformation from different times through the earthquake cycle. The timing of the last earthquake on many of these faults is poorly known, and may vary greatly from fault to fault. Assuming that the underlying mechanics of the seismic cycle are similar for all faults, geodetic observations from different faults may be interpreted as samples over a significantly larger fraction of the earthquake cycle than could be obtained from the geodetic record along any one fault alone. As an ensemble, we find that geologically and geodetically inferred slip rates agree well with a linear relation of 0.94±0.09. To simultaneously explain both the ensemble agreement between geologic and geodetic slip‐rate estimates with observations of rapid postseismic deformation, we consider the predictions from simple two‐layer earthquake‐cycle models with both Maxwell and Burgers viscoelastic rheologies. We find that a two‐layer Burgers model, with two relaxation timescales, is consistent with observations of deformation throughout the earthquake cycle, whereas the widely used two‐layer Maxwell model with a single relaxation timescale, is not, suggesting that the earthquake cycle is effectively characterized by a largely stress‐recoverable rapid postseismic stage and a much more slowly varying interseismic stage.Earth and Planetary Science

Geodynamics Branch research report, 1982

The research program of the Geodynamics Branch is summarized. The research activities cover a broad spectrum of geoscience disciplines including space geodesy, geopotential field modeling, tectonophysics, and dynamic oceanography. The NASA programs which are supported by the work described include the Geodynamics and Ocean Programs, the Crustal Dynamics Project, the proposed Ocean Topography Experiment (TOPEX) and Geopotential Research Mission. The individual papers are grouped into chapters on Crustal Movements, Global Earth Dynamics, Gravity Field Model Development, Sea Surface Topography, and Advanced Studies

Interseismic strain accumulation: Spin-up, cycle invariance, and irregular rupture sequences

Using models of infinite length strike-slip faults in an elastic layer above linear viscoelastic regions, we investigate interseismic deformation. In the models we investigate, interseismic strain accumulation on mature faults is the result of the cumulative effects of all previous ruptures and is independent of the fault loading conditions. The time for a fault to spin-up to a mature state depends on the rheologies and the fault loading conditions. After the model has spun-up, the temporal variation of shear stresses is determined by the fault slip rate and model rheologies. The change in stress during spin-up depends on the slip rate, rheologies, and fault loading conditions but is independent of the magnitude of the initial stress. Over enough cycles such that the cumulative deformation is block-like, the average mature interseismic velocities are equal to the interseismic velocities of an elastic model with the same geometry and distribution of shear moduli. In a model that has spun-up with the fault rupturing periodically, the cumulative deformation is block-like at the end of each seismic cycle, and the interseismic deformation is cycle-invariant (i.e., the same in all cycles). When the fault ruptures nonperiodically, the fault spins up to a mature state that is the same as if the fault had ruptured periodically with the mean slip rate. When the fault slip rate within each cycle varies, the interseismic deformation evolves toward the cycle-invariant deformation determined by the most recent fault slip rate. Around a fault whose slip rate has been faster (slower) than average, interseismic velocities are larger (smaller) than the cycle-invariant velocities and increase (decrease) from cycle to cycle

The interpretation of crustal dynamics data in terms of plate interactions and active tectonics of the Anatolian plate and surrounding regions in the Middle East

A detailed study was made of the consequences of the Arabian plate convergence against Eurasia and its effects on the tectonics of Anatolia and surrounding regions of the eastern Mediterranean. A primary source of information is time rates of change of baseline lengths and relative heights determined by repeated SLR measurements. These SLR observations are augmented by a network of GPS stations in Anatolia, Aegea, and Greece, established and twice surveyed since 1988. The existing SLR and GPS networks provide the spatial resolution necessary to reveal the details of ongoing tectonic processes in this area of continental collision. The effort has involved examining the state of stress in the lithosphere and relative plate motions as revealed by these space based geodetic measurements, seismicity, and earthquake mechanisms as well as the aseismic deformations of the plates from conventional geodetic data and geological evidence. These observations are used to constrain theoretical calculations of the relative effects of: (1) the push of the Arabian plate; (2) high topography of Eastern Anatolia; (3) the geometry and properties of African-Eurasian plate boundary; (4) subduction under the Hellenic Arc and southwestern Turkey; and (5) internal deformation and rotation of the Anatolian plate

Viscoelastic deformation near active plate boundaries

Very Long Baseline Interferometry (VLBI) now has the capacity to monitor geodetic positions with precisions of a few 1 mm over continental baselines. For tectonic applications, one of the major products of the VLBI program is the determination of the rate of change of station locations. Vector site velocities are now routinely produced. One of the novel techniques, VLBI Euler poles, is discussed