Source parameters for historical earthquakes worldwide are compiled to develop a series of empirical relationships among moment magnitude ( M ), surface rupture length, subsurface rupture length, downdip rupture width, rupture area, and maximum and average displacement per event. The resulting data base is a significant update of previous compilations and includes the additional source parameters of seismic moment, moment magnitude, subsurface rupture length, downdip rupture width, and average surface displacement. Each source parameter is classified as reliable or unreliable, based on our evaluation of the accuracy of individual values. Only the reliable source parameters are used in the final analyses. In comparing source parameters, we note the following trends: (1) Generally, the length of rupture at the surface is equal to 75% of the subsurface rupture length; however, the ratio of surface rupture length to subsurface rupture length increases with magnitude; (2) the average surface displacement per event is about one-half the maximum surface displacement per event; and (3) the average subsurface displacement on the fault plane is less than the maximum surface displacement but more than the average surface displacement. Thus, for most earthquakes in this data base, slip on the fault plane at seismogenic depths is manifested by similar displacements at the surface. Log-linear regressions between earthquake magnitude and surface rupture length, subsurface rupture length, and rupture area are especially well correlated, showing standard deviations of 0.25 to 0.35 magnitude units. Most relationships are not statistically different (at a 95% significance level) as a function of the style of faulting: thus, we consider the regressions for all slip types to be appropriate for most applications. Regressions between magnitude and displacement, magnitude and rupture width, and between displacement and rupture length are less well correlated and have larger standard deviation than regressions between magnitude and length or area. The large number of data points in most of these regressions and their statistical stability suggest that they are unlikely to change significantly in response to additional data. Separating the data according to extensional and compressional tectonic environments neither provides statistically different results nor improves the statistical significance of the regressions. Regressions for cases in which earthquake magnitude is either the independent or the dependent parameter can be used to estimate maximum earthquake magnitudes both for surface faults and for subsurface seismic sources such as blind faults, and to estimate the expected surface displacement along a fault for a given size earthquake.

New Empirical Relationships among Magnitude, Rupture Length, Rupture Width, Rupture Area, and Surface Displacement

A global plate motion model, named NUVEL-1, which describes current plate motions between 12 rigid plates is described, with special attention given to the method, data, and assumptions used Tectonic implications of the patterns that emerged from the results are discussed It is shown that wide plate boundary zones can form not only within the continental lithosphere but also within the oceanic lithosphere; eg, between the Indian and Australian plates and between the North American and South American plates Results of the model also suggest small but significant diffuse deformation of the oceanic lithosphere, which may be confined to small awkwardly shaped salients of major plates

Current plate motions

SUMMARY 
We determine best-fitting Euler vectors, closure-fitting Euler vectors, and a new global model (NUVEL-1) describing the geologically current motion between 12 assumed-rigid plates by inverting plate motion data we have compiled, critically analysed, and tested for self-consistency. We treat Arabia, India and Australia, and North America and South America as distinct plates, but combine Nubia and Somalia into a single African plate because motion between them could not be reliably resolved. The 1122 data from 22 plate boundaries inverted to obtain NUVEL-1 consist of 277 spreading rates, 121 transform fault azimuths, and 724 earthquake slip vectors. We determined all rates over a uniform time interval of 3.0m.y., corresponding to the centre of the anomaly 2A sequence, by comparing synthetic magnetic anomalies with observed profiles. The model fits the data well. Unlike prior global plate motion models, which systematically misfit some spreading rates in the Indian Ocean by 8–12mm yr−1, the systematic misfits by NUVEL-1 nowhere exceed ∼3 mm yr−1. The model differs significantly from prior global plate motion models. For the 30 pairs of plates sharing a common boundary, 29 of 30 P071, and 25 of 30 RM2 Euler vectors lie outside the 99 per cent confidence limits of NUVEL-1. Differences are large in the Indian Ocean where NUVEL-1 plate motion data and plate geometry differ from those used in prior studies and in the Pacific Ocean where NUVEL-1 rates are systematically 5–20 mm yr−1 slower than those of prior models. The strikes of transform faults mapped with GLORIA and Seabeam along the Mid-Atlantic Ridge greatly improve the accuracy of estimates of the direction of plate motion. These data give Euler vectors differing significantly from those of prior studies, show that motion about the Azores triple junction is consistent with plate circuit closure, and better resolve motion between North America and South America. Motion of the Caribbean plate relative to North or South America is about 7 mm yr−1 slower than in prior global models. Trench slip vectors tend to be systematically misfit wherever convergence is oblique, and best-fitting poles determined only from trench slip vectors differ significantly from their corresponding closure-fitting Euler vectors. The direction of slip in trench earthquakes tends to be between the direction of plate motion and the normal to the trench strike. Part of this bias may be due to the neglect of lateral heterogeneities of seismic velocities caused by cold subducting slabs, but the larger part is likely caused by independent motion of fore-arc crust and lithosphere relative to the overriding plate.

Empirical relations involving seismic moment M_o, magnitude M_S, energy E_S and fault dimension L (or area S) are discussed on the basis of an extensive set of earthquake data (M_S ≧ 6) and simple crack and dynamic dislocation models. The relation between log S and log M_o is remarkably linear (slope ∼ 2/3) indicating a constant stress drop Δσ; Δσ = 30, 100 and 60 bars are obtained for inter-plate, intra-plate and “average” earthquakes, respectively. Except for very large earthquakes, the relation M_S ∼ (2/3) log M_o ∼ 2 log L is established by the data. This is consistent with the dynamic dislocation model for point dislocation rise times and rupture times of most earthquakes. For very large earthquakes M_S ∼ (1/3) log M_o ∼ log L ∼ (1/3) log E_S. For very small earthquakes M_S ∼ log M_o ∼ 3 log L ∼ log E_S. Scaling rules are assumed and justified. This model predicts log E_S ∼ 1.5 M_S ∼ 3 log L which is consistent with the Gutenberg-Richter relation. Since the static energy is proportional to σL^3, where σ is the average stress, this relation suggests a constant apparent stress ησ where η is the efficiency. The earthquake data suggest ησ ~ 1/2 Δσ. These relations lead to log S ∼ M_S consistent with the empirical relation. This relation together with a simple geometrical argument explains the magnitude-frequency relation log N ∼ − M_S.

/pdf/theoretical-basis-of-some-empirical-relations-in-seismology-1e9arnhgn6.pdf

Theoretical basis of some empirical relations in seismology

Seismic techniques provide the highest-resolution measurements of the structure of the crust and have been conducted on a worldwide basis. We summarize the structure of the continental crust based on the results of seismic refraction profiles and infer crustal composition as a function of depth by comparing these results with high-pressure laboratory measurements of seismic velocity for a wide range of rocks that are commonly found in the crust. The thickness and velocity structure of the crust are well correlated with tectonic province, with extended crust showing an average thickness of 30.5 km and orogens an average of 46.3 km. Shields and platforms have an average crustal thickness nearly equal to the global average. We have corrected for the nonuniform geographical distribution of seismic refraction profiles by estimating the global area of each major crustal type. The weighted average crustal thickness based on these values is 41.1 km. This value is 10% to 20% greater than previous estimates which underrepresented shields, platforms, and orogens. The average compressional wave velocity of the crust is 6.45 km/s, and the average velocity of the uppermost mantle (Pn velocity) is 8.09 km/s. We summarize the velocity structure of the crust at 5-km depth intervals, both in the form of histograms and as an average velocity-depth curve, and compare these determinations with new measurements of compressional wave velocities and densities of over 3000 igneous and metamorphic rock cores made to confining pressures of 1 GPa. On the basis of petrographic studies and chemical analyses, the rocks have been classified into 29 groups. Average velocities, densities, and standard deviations are presented for each group at 5-km depth intervals to crustal depths of 50 km along three different geotherms. This allows us to develop a model for the composition of the continental crust. Velocities in the upper continental crust are matched by velocities of a large number of lithologies, including many low-grade metamorphic rocks and relatively silicic gneisses of amphibolite facies grade. In midcrustal regions, velocity gradients appear to originate from an increase in metamorphic grade, as well as a decrease in silica content. Tonalitic gneiss, granitic gneiss, and amphibolite are abundant midcrustal lithologies. Anisotropy due to preferred mineral orientation is likely to be significant in upper and midcrustal regions. The bulk of the lower continental crust is chemically equivalent to gabbro, with velocities in agreement with laboratory measurements of mafic granulite. Garnet becomes increasingly abundant with depth, and mafic garnet granulite is the dominant rock type immediately above the Mohorovicic discontinuity. Average compressional wave velocities of common crustal rock types show excellent correlations with density. The mean crustal density calculated from our model is 2830 kg/m3, and the average SiO2 content is 61.8%.

Seismic velocity structure and composition of the continental crust: A global view

The Chilean earthquake sequence of May 21–22, 1960, was accompanied by linear zones of tectonic warping, including both uplift and subsidence relative to sea level. The region involved is more than 200 km wide and about 1000 km long, and lies along the continental margin between latitude 37° and 48° S. Significant horizontal strains accompanied the vertical movements in parts of the subsided zone for which triangulation data are available. Displacements were initiated near the northern end of the deformed region during the opening earthquake of the sequence (M s ≅ 7.5) on May 21 at 10h 02m 50s GMT and were extended over the remainder of the region during the culminating shock (M s ≅ 8.5) on May 22 at 19h llm 17s GMT. During the latter event, sudden uplift of adjacent portions of the continental shelf and much or all of the continental slope apparently generated the destructive tsunami that immediately followed the main shock. Available data suggest that the primary fault or zone of faulting along which displacement occurred probably is a complex thrust fault roughly 1000 km long and at least 60 km wide; it dips eastward at a moderate angle beneath the continental margin and intersects the surface on the continental slope. Dip slip required to satisfy the surface displacements is at least 20 m and perhaps as large as 40 m. There is some evidence that there was a minor component of right-lateral slip on the fault plane.

Mechanism of the Chilean Earthquakes of May 21 and 22, 1960

The 1964 Alaskan earthquake (Ms ≈ 8.4) involved a segment of the eastern Aleutian arc 800 gm long; the 1960 Chilean earthquake sequence (Ms ≈ 8.5) affected roughly 100 km of the southern Peru-Chile arc. These two major events are strikingly similar in that (1) seismicity was shallow (<70 km), the earthquake focal regions and most of the associated tectonic deformation being between the oceanic trenches and volcanic chains of the two arcs; (2) regional vertical displacements were characterized by broad asymmetric downwarps elongate parallel to the arcs with flanking zones of marked uplift on the seaward sides and minor, possibly local, uplift on the landward sides; and (3) horizontal displacements, where determined by retriangulation, involved systematic shifts in a generally seaward direction and transverse tensile strains across the zones of subsidence. Surface displacements and seismicity for both events are compatible with dislocation models involving predominantly dip-slip movement of 20 meters or more on major complex thrust faults (megathrusts) inclined at average angles of about 9° beneath the eastern Aleutian arc and perhaps 20° beneath the Peru-Chile arc. The thrust-fault mechanism deduced for both the Alaskan and Chilean earthquakes is broadly consistent with the concept that the sectors of the Pacific rim in which they occurred are major zones of convergence along which the oceanic plates progressively underthrust the less mobile America plate. Directions of convergence between lithospheric plates at these arcs as deduced primarily from paleomagnetic data are in reasonably good agreement with the observed earthquake-related deformation; the deduced rates of convergence, however, appear to be too high in the eastern Aleutian arc and too low in the southern Peru-Chile arc. Despite gross similarities in tectonic setting and the present style of earthquake-related deformation, the geologies of the continental margins in the eastern Aleutian arc and southern Peru-Chile arc differ significantly. This difference suggests that Mesozoic and Cenozoic sediments and volcanic rocks conveyed into the eastern Aleutian trench have progressively accreted to the Alaskan continental margin, whereas most or all of the material carried into the southern Peru-Chile trench has disappeared beneath the Chilean continental margin.

Alaskan Earthquake of 1964 and Chilean Earthquake of 1960: Implications for Arc Tectonics

The MW (moment magnitude) 7.9 Denali fault earthquake on 3 November 2002 was associated with 340 kilometers of surface rupture and was the largest strike-slip earthquake in North America in almost 150 years. It illuminates earthquake mechanics and hazards of large strike-slip faults. It began with thrusting on the previously unrecognized Susitna Glacier fault, continued with right-slip on the Denali fault, then took a right step and continued with right-slip on the Totschunda fault. There is good correlation between geologically observed and geophysically inferred moment release. The earthquake produced unusually strong distal effects in the rupture propagation direction, including triggered seismicity.

https://science.sciencemag.org/content/sci/300/5622/1113.full.pdf

The 2002 Denali fault earthquake, Alaska: a large magnitude, slip-partitioned event.

The Trans-Alaskan Crustal Transect in the southern Copper River Basin and Chugach Mountains traverses the margins of the Peninsular and Wrangellia terranes, and the adjacent accretionary oceanic units of the Chugach terrane to the south The southern Wrangellia terrane margin consists of a polymetamorphosed magmatic arc complex at least in part of Pennsylvanian age (Strelna Metamorphics and metagranodiorite) and tonalitic metaplutonic rocks of the Late Jurassic Chitina magmatic arc The southern Peninsular terrane margin is underlain by rocks of the Late Triassic (?) and Early Jurassic Talkeetna magmatic arc (Talkeetna Formation and Border Ranges ultra-mafic-mafic assemblage) on Permian or older basement rocks The Peninsular and Wrangellia terranes are parts of a dominantly oceanic superterrane (composite Terrane II) that was amalgamated by Late Triassic time and was accreted to terranes of continental affinity north of the Denali fault system in the mid- to Late Cretaceous The Chugach terrane in the transect area consists of three successively accreted units: (1) minor greenschist and intercalated blueschist, the schist of Liberty Creek, of unknown protolith age that was metamorphosed and probably accreted during the Early Jurassic, (2) the McHugh Complex (Late Triassic to mid-Cretaceous protolith age), a melange of mixed oceanic, volcaniclastic, and olistostromal rocks that is metamorphosed to prehnite-pumpellyite and lower greenschist facies that was accreted by middle Cretaceous time, and (3) the Upper Cretaceous Valdez Group, mainly magmatic arc-derived flysch and lesser oceanic volcanic rocks of greenschist facies that was accreted by early Paleocene time A regional thermal event that culminated in early middle Eocene time (48–52 Ma) resulted in widespread greenschist facies metamorphism and plutonism

Bedrock geology and tectonic evolution of the Wrangellia, Peninsular, and Chugach Terranes along the Trans‐Alaska Crustal Transect in the Chugach Mountains and Southern Copper River Basin, Alaska

George Plafker

Papers

Mechanism of the Chilean Earthquakes of May 21 and 22, 1960

Alaskan Earthquake of 1964 and Chilean Earthquake of 1960: Implications for Arc Tectonics

The 2002 Denali fault earthquake, Alaska: a large magnitude, slip-partitioned event.

Bedrock geology and tectonic evolution of the Wrangellia, Peninsular, and Chugach Terranes along the Trans‐Alaska Crustal Transect in the Chugach Mountains and Southern Copper River Basin, Alaska

Geology of the southern Alaska margin