Rupture process of large earthquakes in the northern Mexico subduction zone

The Cocos plate subducts beneath North America at the Mexico trench. The northernmost segment of this trench, between the Orozco and Rivera fracture zones, has ruptured in a sequence of five large earthquakes from 1973 to 1985; the Jan. 30, 1973 Colima event (M _s 7.5) at the northern end of the segment near Rivera fracture zone; the Mar. 14, 1979 Petatlan event (M _s 7.6) at the southern end of the segment on the Orozco fracture zone; the Oct. 25, 1981 Playa Azul event (M _s 7.3) in the middle of the Michoacan gap; the Sept. 19, 1985 Michoacan mainshock (M _s 8.1); and the Sept. 21, 1985 Michoacan aftershock (M _s 7.6) that reruptured part of the Petatlan zone. Body wave inversion for the rupture process of these earthquakes finds the best: earthquake depth; focal mechanism; overall source time function; and seismic moment, for each earthquake. In addition, we have determined spatial concentrations of seismic moment release for the Colima earthquake, and the Michoacan mainshock and aftershock. These spatial concentrations of slip are interpreted as asperities; and the resultant asperity distribution for Mexico is compared to other subduction zones. The body wave inversion technique also determines theMoment Tensor Rate Functions; but there is no evidence for statistically significant changes in the moment tensor during rupture for any of the five earthquakes. An appendix describes theMoment Tensor Rate Functions methodology in detail.The systematic bias between global and regional determinations of epicentral locations in Mexico must be resolved to enable plotting of asperities with aftershocks and geographic features. We have spatially shifted all of our results to regional determinations of epicenters. The best point source depths for the five earthquakes are all above 30 km, consistent with the idea that the down-dip edge of the seismogenic plate interface in Mexico is shallow compared to other subduction zones. Consideration of uncertainties in the focal mechanisms allows us to state that all five earthquakes occurred on fault planes with the same strike (N65°W to N70°W) and dip (15±3°), except for the smaller Playa Azul event at the down-dip edge which has a steeper dip angle of 20 to 25°. However, the Petatlan earthquake does prefer a fault plane that is rotated to a more east-west orientation—one explanation may be that this earthquake is located near the crest of the subducting Orozco fracture zone. The slip vectors of all five earthquakes are similar and generally consistent with the NUVEL-predicted Cocos-North America convergence direction of N33°E for this segment. The most important deviation is the more northerly slip direction for the Petatlan earthquake. Also, the slip vectors from the Harvard CMT solutions for large and small events in this segment prefer an overall convergence direction of about N20°E to N25°E.All five earthquakes share a common feature in the rupture process: each earthquake has a small initial precursory arrival followed by a large pulse of moment release with a distinct onset. The delay time varies from 4 s for the Playa Azul event to 8 s for the Colima event. While there is some evidence of spatial concentration of moment release for each event, our overall asperity distribution for the northern Mexico segment consists of one clear asperity, in the epicentral region of the 1973 Colima earthquake, and then a scattering of diffuse and overlapping regions of high moment release for the remainder of the segment. This character is directly displayed in the overlapping of rupture zones between the 1979 Petatlan event and the 1985 Michoacan aftershock. This character of the asperity distribution is in contrast to the widely spaced distinct asperities in the northern Japan-Kuriles Islands subduction zone, but is somewhat similar to the asperity distributions found in the central Peru and Santa Cruz Islands subduction zones. Subduction of the Orozco fracture zone may strongly affect the seismogenic character as the overlapping rupture zones are located on the crest of the subducted fracture zone. There is also a distinct change in the physiography of the upper plate that coincides with the subducting fracture zone, and the Guerrero seismic gap to the south of the Petatlan earthquake is in the wake of the Orozco fracture zone. At the northern end, the Rivera fracture zone in the subducting plate and the Colima graben in the upper plate coincide with the northernmost extent of the Colima rupture zone.     

Seismic coupling and uncoupling at subduction zones

Seismic coupling has been used as a qualitative measure of the “interaction” between the two plates at subduction zones. Kanamori (1971) introduced seismic coupling after noting that the characteristic size of earthquakes varies systematically for the northern Pacific subduction zones. A quantitative global comparison of many subduction zones reveals a strong correlation of earthquake size with two other variables: age of the subducting lithosphere and convergence rate. The largest earthquakes occur in zones with young lithosphere and fast convergence rates, while zones with old lithosphere and slow rates are relatively aseismic for large earthquakes. Results from a study of the rupture process of three great earthquakes indicate that maximum earthquake size is directly related to the asperity distribution on the fault plane (asperities are strong regions that resist the motion between the two plates). The zones with the largest earthquakes have very large asperities, while the zones with smaller earthquakes have small scattered asperities. This observation can be translated into a simple model of seismic coupling, where the horizontal compressive stress between the two plates is proportional to the ratio of the summed asperity area to the total area of the contact surface. While the variation in asperity size is used to establish a connection between earthquake size and tectonic stress, it also implies that plate age and rate affect the asperity distribution. Plate age and rate can control asperity distribution directly by use of the horizontal compressive stress associated with the “preferred trajectory” (i.e. the vertical and horizontal velocities of subducting slabs are determined by the plate age and convergence velocity). Indirect influences are many, including oceanic plate topography and the amount of subducted sediments.All subduction zones are apparently uncoupled below a depth of about 40 km, and we propose that the basalt to eclogite phase change in the down-going oceanic crust may be largely responsible. This phase change should start at a depth of 30–35 km, and could at least partially uncouple the plates by superplastic deformation throughout the oceanic crust during the phase change.     

Seismicity and the subduction process

There is considerable variation between subduction zones in the largest characteristic earthquake within each zone. Assuming that coupling between downgoing and upper plates is directly related to characteristic earthquake size, we have tested for correlations between variation in coupling and other physical features of subduction zones: the lateral extent and penetration depth of Benioff zones, age of subducting lithosphere, convergence rate, and back-arc spreading. Using linear multivariate regression, coupling is correlated with two variables: convergence rate and lithosphere age. Secondary correlations within the data set are penetration depth versus lithosphere age, and lateral extent versus convergence rate. An important additional correlation is that back-arc spreading is found to be associated with subduction zones where coupling is low (those characterised by small earthquakes). Taken together, the observed correlations suggest a simple qualitative model where convergence rate and lithosphere age determine the horizontal and sinking rates, respectively, of slabs: these parameters influence the seismic coupling in the subduction zone. In the limit of a fast sinking rate and slow convergence rate, back-arc spreading occurs and thereby appears to be a passive process.     

Distribution and characteristics of Adirondack-class basalt as observed by Mini-TES in Gusev crater,Mars and its possible volcanic source

We identified 104 unique rock targets belonging to the olivine-rich Adirondack class using Mini-TES data. Rare rocks on the West Spur of the Columbia Hills and on the plains east of the Hills also belong to this class. We present evidence that Adirondack-class basaltic lavas may have had their origin at Apollinaris Tholus. Linear modeling of Adirondack-class rock spectra shows only minor variations in mineralogy and the primary phases identified are consistent with olivine basalt having an average olivine composition of ～Fo₄₅. We used factor analysis and target transformation to identify variability within a single class of rocks for the first time, and we find that olivine abundance varies independently of a basaltic matrix. A spectral component previously attributed to downwelling radiance is proposed here to be attributable to optically thin dust on rocks, which has a greater effect on Mini-TES spectra than previously recognized, but can explain the sloping continuum observed in Adirondack-class rock spectra as well as the apparent distortion of an olivine absorption. Spectral mixing is not completely linear, leading to the overestimation of sulfate and olivine fractions and the underestimation of plagioclase feldspar, although linear mixing appears to successfully replicate the majority of the observed signal. Published TES spectra of low albedo, low dust cover, olivine-bearing materials do not exhibit the spectral telltales of non-linear behavior, probably because orbital spectra represent mostly mobile regolith that has not accumulated a sufficient dust cover over ～3 × 6 km areas.     

Landslide susceptibility analysis with a heuristic approach in the Eastern Alps (Vorarlberg, Austria)

The aim of the study is landslide hazard assessment carried out on a working scale of 1:25 000. The study area within the Northern Calcareous Alps was geologically and geotechnically mapped in order to identify causes and mechanisms of active mass movements. The field surveys were digitised by a Geographical Information System and divided into data layers. The geological units were classified according to their geotechnical properties. All layers were converted into grids and spatially analysed together with a Digital Elevation Model. Comparing the layers with the inventory of active landslides, the prevailing factors leading to sliding movements were identified. Because of the complex tectonic setting and the small number of active landslides, a statistical method of hazard assessment was not applicable. Using the heuristic approach of an index method, the data layers of geotechnical class, bedding conditions, tectonic layouts, slope angles, slope orientations, vegetation and erosion were analysed. The susceptibility of each layer has been evaluated with help of bivariate statistics. The layers have been weighted with indices due to their importance iteratively and were combined into a landslide susceptibility map.     

The rupture process and asperity distribution of three great earthquakes from long-period diffracted P-waves

Maximum earthquake size varies considerably amongst the subduction zones. This has been interpreted as a variation in the seismic coupling, which is presumably related to the mechanical conditions of the fault zone. The rupture process of a great earthquake indicates the distribution of strong (asperities) and weak regions of the fault. The rupture process of three great earthquakes (1963 Kurile Islands, M_W = 8.5; 1965 Rat Islands, M_W = 8.7; 1964 Alaska, M_W = 9.2) are studied by using WWSSN stations in the core shadow zone. Diffraction around the core attenuates the P-wave amplitudes such that on-scale long-period P-waves are recorded. There are striking differences between the seismograms of the great earthquakes; the Alaskan earthquake has the largest amplitude and a very long-period nature, while the Kurile Islands earthquake appears to be a sequence of magnitude 7.5 events.The source time functions are deconvolved from the observed records. The Kurile Islands rupture process is characterized by the breaking of asperities with a length scale of 40–60 km, and for the Alaskan earthquake the dominant length scale in the epicentral region is 140–200 km. The variation of length scale and M_W suggests that larger asperities cause larger earthquakes. The source time function of the 1979 Colombia earthquake (M_W = 8.3) is also deconvolved. This earthquake is characterized by a single asperity of length scale 100–120 km, which is consistent with the above pattern, as the Colombia subduction zone was previously ruptured by a great (M_W = 8.8) earthquake in 1906.The main result is that maximum earthquake size is related to the asperity distribution on the fault. The subduction zones with the largest earthquakes have very large asperities (e.g. the Alaskan earthquake), while the zones with the smaller great earthquakes (e.g. Kurile Islands) have smaller scattered asperities.     

What controls the lateral variation of large earthquake occurrence along the Japan Trench?

Abstract The lateral (along trench axis) variation in the mode of large earthquake occurrence near the northern Japan Trench is explained by the variation in surface roughness of the subducting plate. The surface roughness of the ocean bottom near the trench is well correlated with the large-earthquake occurrence. The region where the ocean bottom is smooth is correlated with'typical'large underthrust earthquakes (e.g. the 1968 Tokachioki event) in the deeper part of the seismogenic plate interface, and there are no earthquakes in the shallow part (aseismic zone). The region where the ocean bottom is rough (well-developed horst and graben structure) is correlated with large normal faulting earthquakes (e.g. the 1933 Sanriku event) in the outer-rise region, and large tsunami earthquakes (e.g. the 1896 Sanriku event) in the shallow region of the plate interface zone. In the smooth surface region, the coherent metamorphosed sediments form a homogeneous, large and strong contact zone between the plates. The rupture of this large strong contact causes great under-thrust earthquakes. In the rough surface region, large outer-rise earthquakes enhance the well-developed horst and grabens. As these structure are subducted with sediments in the graben part, the horsts create enough contact with the overriding block to cause an earthquake in the shallow part of the interface zone, and this earthquake is likely to be a tsunami earthquake. When the horst and graben structure is further subducted, many small strong contacts between the plates are formed, and they can cause only small underthrust earthquakes.     

Large earthquakes in the macquarie ridge complex: Transitional tectonics and subduction initiation

While most aspects of subduction have been extensively studied, the process of subduction initiation lacks an observational foundation. The Macquarie Ridge complex (MRC) forms the Pacific-Australia plate boundary between New Zealand to the north and the Pacific-Australia-Antarctica triple junction to the south. The MRC consists of alternating troughs and rises and is characterized by a transitional tectonic environment in which subduction initiation presently occurs. There is a high seismicity level with 15 large earthquakes (M>7) in this century. Our seismological investigation is centered on the largest event since 1943: the 25 MAY 1981 earthquake. Love, Rayleigh, andP waves are inverted to find: a faulting geometry of right-lateral strike-slip along the local trend of the Macquarie Ridge (N30°E); a seismic moment of 5×10²⁷ dyn cm (M _w=7.7) a double event rupture process with a fault length of less than 100km to the southwest of the epicenter and a fault depth of less than 20km. Three smaller thrust earthquakes occurred previous to the 1981 event along the 1981 rupture zone; their shallow-dipping thrust planes are virtually adjacent to the 1981 vertical fault plane. Oblique convergence in this region is thus accommodated by a dual rupture mode of several small thrust events and a large strike-slip event. Our study of other large MRC earthquakes, plus those of other investigators, produces focal mechanisms for 15 earthquakes distributed along the entire MRC; thrust and right-lateral strike-slip events are scattered throughout the MRC. Thus, all of the MRC is characterized by oblique convergence and the dual rupture mode. The true best-fit rotation pole for the Pacific-Australia motion is close to the Minster & Jordan RM2 pole for the Pacific-India motion. Southward migration of the rotation pole has caused the recent transition to oblique convergence in the northern MRC. We propose a subduction initiation process that is akin to crack propagation; the 1981 earthquake rupture area is identified as the crack-tip region that separates a disconnected mosaic of small thrust faults to the south from a horizontally continuous thrust interface to the north along the Puysegur trench. A different mechanism of subduction initiation occurs in the southernmost Hjort trench region at the triple junction. newly created oceanic lithosphere has been subducted just to the north of the triple junction. The entire MRC is a soft plate boundary that must accommodate the plate motion mismatch between two major spreading centers (Antarctica-Australia and Pacific-Antarctica). The persistence of spreading motion at the two major spreading centers and the consequent evolution of the three-plate system cause the present-day oblique convergence and subduction initiation in the Macquarie Ridge complex.     

Spectral evidence for zeolite in the dust on Mars

Spectral features observed in Mars Global Surveyor Thermal Emission Spectrometer data (∼1670-220 cm⁻¹) of martian surface dust provide clues to its mineralogy. An emissivity peak at ∼1630 cm⁻¹ is consistent with the presence of an H₂O-bearing mineral. This spectral feature can be mapped globally and shows a distribution related to the classical bright regions on Mars that are known to be dust covered. An important spectral feature at ∼830 cm⁻¹ present in a newly derived average spectrum of surface dust likely is a transparency feature arising from the fine particulate nature of the dust. Its shape and location are consistent with plagioclase feldspars and also zeolites, which essentially are the hydrous form of feldspar. The generally favored visible/near-infrared spectral analog for martian dust, JSC Mars-1 altered tephra, does not display the ∼830 cm⁻¹ feature. Zeolites commonly form from the interaction of low temperature aqueous fluids and volcanic glass in a variety of geologic settings. The combination of spectral features that are consistent with zeolites and the likelihood that Mars has (or had) geologic conditions necessary to produce them makes a strong case for recognizing zeolite minerals as likely components of the martian regolith.