Direction Of Plate Convergence Research Articles

Oblique plate convergence along arcuate trenches on a spherical Earth. An example from the Western Sunda Arc

We propose a new definition of oblique plate convergence. Our model considers arc-trench curvature on a spherical Earth, and defines obliquity as the angle between the direction of plate convergence and the azimuth from the central point along the arc to the center of curvature of the arcuate trench. We also show how this model may be extended to multiple-arc trenches.We apply this definition to the Western Sunda Arc, which has been traditionally considered the prime example of oblique plate convergence, particularly the segment corresponding to Sumatra. We define eight arc-like segments along this arc-trench system. Obliquity angles on these arc-like segments vary from 92° in the northernmost (Burmese) arc to 7° offshore the Sunda strait (between Sumatra and Java). The angle of obliquity offshore Sumatra is 31°, almost within the definition of slightly oblique convergence. Given the low angle of obliquity, partitioning along the Sumatran fault, in the southern portion of the Western Sunda Arc, is difficult to explain by oblique plate convergence alone. We suggest that the entire segment (or forearc) inboard of the Western Sunda Arc be dragged from the Burmese arc, where overriding and subducting plates are completely coupled and oblique plate convergence is high.

Laramide crustal detachment in the Rockies: Cordilleran shortening of fluid-weakened foreland crust

ABSTRACTWhat causes previously stable continental crust in the forelands of Cordilleran orogenic systems to shorten during low-angle subduction? The National Science Foundation/EarthScope Bighorn Project combined seismic imaging of the crust and Moho with kinematic modeling of Laramide (Late Cretaceous–Paleogene) basement-involved deformation to address this question. In north-central Wyoming, asymmetrical ENE-verging upper-crustal folds are highly discordant with broader, N-trending warps in the Moho, indicating crustal detachment. Restorable cross sections of ENE-directed detachment at a depth of ~30 km, combined a smaller component of NNW–SSE shortening due to the east-narrowing shape of the crustal allochthon, can explain the anastomosing network of Laramide basement-cored arches without major deformation of the underlying mantle lithosphere.Thrust-related fold geometries and west-to-east initiation of deformation in the Laramide and Sevier thrust belts point to Cordilleran end-loading from the west. Differences between Laramide (~N65E) and plate (~N25E) convergence directions, along with the fanning of Laramide shortening directions from nearly E–W to the south to NE–SW to the north, indicate slip partitioning during end-loading west of the Rockies.Sub-horizontal detachment with a near-zero critical taper within cratonic crust suggests an extremely weak Laramide detachment zone during deformation. Analogous lower-crustal deformation in subduction forearcs is associated with slow earthquakes and slab dehydration. We hypothesize that low-angle subduction of the Farallon Plate suppressed fluid-consuming melting and corner-flow processes that characterize higher-angle subduction. This allowed subduction-generated fluids to escape upward into the overlying continental lithosphere, causing retrograde metamorphism and increased fluid pressure that facilitated crustal detachment. This hydration-based hypothesis predicts that crustal detachment will accompany major earthquakes in active analog orogens.

Strain Partitioning in the Southern Ryukyu Margin Revealed by Seafloor Geodetic and Seismological Observations

AbstractThe southern Ryukyu subduction zone is one of the potential sources for tsunamigenic earthquakes. Despite a great seismic risk, the deformation pattern remains poorly known, primarily due to the absence of seafloor constraints. With GNSS‐acoustic measurements over years, we characterize the convergence rate across this margin growing from 92 mm/yr offshore eastern Taiwan to 123 mm/yr near the Gagua Ridge. The new data suggest the subduction interface is capable of hosting Mw 7.5–8.4 earthquakes. The orientations of seafloor movement and P‐axes in the Nanao Basin are both subnormal to the trench, notably deviate from the direction of plate convergence. By considering the combined effect of plate convergence and backarc rifting, different trends between the forearc convergence, P‐axes, and seafloor movement may indicate some degree of slip‐partitioning. The trench‐parallel component is likely accommodated in part by earthquakes near Taiwan, lower plate deformation, and strike‐slip faults within the accretionary wedge.

Structure and Stress Field of the Lithosphere Between Pamir and Tarim

AbstractThe Pamir plateau protrudes ∼300 km between the Tajik‐ and Tarim‐basin lithosphere of Central Asia. Whether its salient location and shape are caused by forceful indentation of a promontory of Indian mantle lithosphere is debated. We present a new local‐seismicity and focal‐mechanism catalog, and a P‐wave velocity model of the eastern part of the collision system. The data suggest a south‐dipping Asian slab that overturns in its easternmost segment. The largest principal stress at depth acts normal on the slab and is orientated parallel to the plate convergence direction. In front (south) of the Asian slab, a volume of mantle with elevated velocities and lined by weak seismicity constitutes the postulated Indian mantle indenter. We propose that the indenter delaminates and overturns the Asian slab, underthrusts the Tarim lithosphere along a compressive transform boundary, and controls the location and shape of the Pamir plateau.

Earthquake Shocks Around Delhi-NCR and the Adjoining Himalayan Front: A Seismotectonic Perspective

An increase in the number of earthquakes and subsequent clustering in northwest India, particularly around the Delhi-National Capital Region (NCR) and adjacent NW Himalayan front, provides a good opportunity to understand the underpinning tectonic controls and the likelihood of any large earthquake in the future. The 2001 Mw7.7 Bhuj, 2011 Mw6.9 Sikkim and 2015 Mw7.8 and 7.3 Nepal earthquakes (and 2004 Mw9.2 Sumatra event) are important in this context. We analyzed the seismicity around the Delhi-NCR and the adjoining Himalayan front, including event clustering and the spatio-temporal distribution of b-values, in the context of kinematics and the regional geodynamics. The overall moderate-to-low b-values, both in time and space, since 2016, provide information regarding an increase and subsequent stabilization of the stress field in the study area. The analysis led to the identification of (1) a structurally guided stress field in the region between the Kachchh and the NW Himalaya that coincides with the direction of Indian plate convergence and (2) frequent occurrences of earthquakes particularly in the Delhi, Kangra and Uttarkashi areas. We propose that faults in western Peninsular India, which pass through the margins of the Aravalli Range, the Marwar basin, and the isostatically over-compensated Indo-Gangetic Plains beneath the under-plated Indian lithosphere, act as stress guides; concentrating and increasing stress in regions of lithospheric flexure. This enhanced stress may trigger a large earthquake.

Modeling In-situ tectonic stress state and maximum horizontal stress azimuth in the Central Algerian Sahara – A geomechanical study from El Agreb, El Gassi and Hassi Messaoud fields

Central Algerian Sahara hosts many prolific hydrocarbon accumulations in the Paleozoic successions. In this work a contemporary stress field of the Saharan platform has been evaluated using the dataset from recently drilled wells in El Agreb, El Gassi and Hassi Messaoud fields. A pore fluid pressure gradient of 0.56 PSI/feet is interpreted from the in-situ measurements in the Paleozoic reservoir units. Vertical stress (Sv) modeled from the bulk-density data indicates an average of 1.02 PSI/feet gradient. Rock elastic property-based approach is employed to model the magnitudes of minimum (Shmin) and maximum horizontal stress (SHmax) components, which were calibrated with leak off test/minifrac and breakout widths, respectively. Paleozoic stress profiles reveal Shmin/Sv range of 0.74–0.84, while SHmax/Sv varies between 1.1 and 1.33. Subsurface stress distribution indicates that the present-day stress field in the Saharan platform is principally strike-slip faulting (SHmax > Sv > Shmin). A cumulative 1490 m of B-D quality wellbore breakouts, inferred from the acoustic image logs, suggest a NW-SE/WNW-ESE SHmax orientation, which is parallel to the absolute African plate motion and Africa-Eurasia plate convergence direction, implying ridge push force to be the dominant contributor to the tectonic stress field. Mean SHmax orientation shows slightly anticlockwise rotation (126°N to 144°N) from south (El Agreb) to north (Hassi Messaoud field). Inferences are discussed regarding the fault slip potential and hydrocarbon reservoir development.

Three‐Dimensional P Wave Velocity Structure of the Northern Hikurangi Margin From the NZ3D Experiment: Evidence for Fault‐Bound Anisotropy

AbstractWe present a high‐resolution three‐dimensional (3‐D) anisotropic P wave velocity (Vp) model in the northern Hikurangi margin offshore Gisborne, New Zealand, constructed by tomographic inversion of over 430,000 first arrivals recorded by a dense grid of ocean bottom seismometers. Since the study area covers a region where shallow slow slip events (SSEs) occur repeatedly and the subduction of a seamount is proposed, it offers an ideal location to link our understanding of structural and hydrogeologic properties at megathrust faults to slip behavior. The Vp model reveals an ~30‐km‐wide, low‐velocity accretionary wedge at the toe of the overriding plate, where previous seismic reflection studies show a series of active thrust faults branching from the plate interface. We find some locations with significant Vp azimuthal anisotropy >5% near the branching faults and the deformation front. This finding suggests that the anisotropy is not ubiquitous and homogeneous within the overriding plate, but more localized in the vicinity of active thrust faults. The fast axes of Vp within the accretionary wedge are mostly oriented to the plate convergence direction, which is interpreted as preferentially oriented cracks in a compressional stress regime associated with plate subduction. We find that the magnitudes of anisotropy are roughly equivalent to values found at oceanic spreading centers, where the extensional stress regime is dominant and the crack density is expected to be higher than subduction zones. This consideration may indicate that additional effects such as fault foliation and clay mineral alignment also contribute to upper plate anisotropy along subduction margins.

Cenozoic increase in subduction erosion during plate convergence variability along the convergent margin off Trujillo, Peru

Subduction erosion is a long-term, large-scale geological process that dominates the structural inventory of a large number of convergent margins. Along the Peruvian convergent margin, this process has controlled the structural deformation of the forearc as a consequence of spatio-temporal variations during the Cenozoic. However, the pace at which tectonic erosion occurred has not been linked to regional, subduction-related events. The aim of this study is to show that subduction erosion along the Central Peruvian convergent margin off Trujillo (7∘S–9∘S) was influenced by plate convergence variability since the Late Cretaceous. Drilling data, multi-channel seismic (MCS) data and subsidence analyses are used to investigate the long-term effect of subduction erosion on the tectonic evolution of the Peru margin. Our study shows that subduction erosion (1) did not occur between the Cretaceous and Early Eocene (110–56 Ma), (2) increased and fluctuated between the Eocene and Pliocene (56–5 Ma) with a maximum average trench retreat of 3 km.Myr−1 between the Middle Miocene and Late Miocene (20–5 Ma), and (3) has decreased since the Pliocene (5 Ma) (1 km.Myr−1 of average trench retreat). A remarkable finding is that shorter periods of faster subduction erosion coincide with or follow major plate reorganization events that were marked by changes in the velocity and direction of relative plate convergence as well as with well-known phases of Andean orogeny. Based on our analyses, we therefore conclude that convergence variability off Peru is a primary mechanism controlling the pace at which tectonic erosion occurred along the Trujillo Basin during the Cenozoic.

Evidence for a large strike-slip component during the 1960 Chilean earthquake

SUMMARYThe strainmeter record observed at Isabella (ISA), California, for the 1960 Chilean earthquake (Mw = 9.5) is one of the most important historical records in seismology because it was one of the three records that provided the opportunity for the first definitive observations of free oscillations of the Earth. Because of the orientation of the strainmeter rod with respect to the back azimuth to Chile, the ISA strainmeter is relatively insensitive to G (Love) waves and higher order (order ≥ 6) toroidal modes, yet long-period G waves and toroidal modes were recorded with large amplitude on this record. This observation cannot be explained with the conventional low-angle thrust mechanism typical of great subduction-zone earthquakes and requires an oblique mechanism with half strike-slip and half thrust. The strain record at Ogdenburg, New Jersey, the Press–Ewing seismograms at Berkeley, California, and the ultra-long period displacement record at Pasadena, California, also support the oblique mechanism. We tested the performance of the ISA strainmeter using other events including the 1964 Alaskan earthquake and found no instrumental problems. Thus, the ISA observation of large G/R and toroidal/spheroidal ratios most likely reflects the real characteristics of the 1960 Chilean earthquake, rather than an observational artefact. The interpretation of the large strike-slip component is not unique, but it may represent release of the strike-slip strain that has accumulated along the plate boundary as a result of oblique convergence at the Nazca–South American plate boundary. The slip direction of the 2010 Chilean (Maule) earthquake ( Mw = 8.8) is rotated by about 10° clockwise from the plate convergence direction suggesting that right-lateral strain comparable to that of an Mw = 8.3 earthquake remained unreleased and accumulates near the plate boundary. One possible scenario is that the strike-slip strain accumulated over several great earthquakes like the 2010 Maule earthquake was released during the 1960 Chilean earthquake. If this is the case, we cannot always expect a similar behaviour for all the great earthquakes occurring in the same subduction zone and such variability needs to be considered in long-term hazard assessment of subduction-zone earthquakes.

Variations in Cascadia Low‐Frequency Earthquake Behavior With Downdip Distance

AbstractWe identify and characterize four different low‐frequency earthquake (LFE) families (LFEs 1–4, ordered updip to downdip) that span the width of the transition zone in the Cascadia Subduction Zone beneath western Washington State. We find that LFE swarm duration, recurrence interval, and event size decrease systematically with increasing depth. LFE moments are observed to follow an exponential distribution rather than a power law distribution, allowing us to determine a characteristic event size for each family. Absolute locations of these LFE families place them very near the inferred Juan de Fuca‐North America plate interface. Relative relocation of individual LFEs within each family reveal elongation parallel to the plate convergence direction and narrow depth distributions with dips corresponding to plate interface models. We interpret the behavior of our LFE families as indicative of slow slip at those locations on the plate interface. The tidal response of LFE activity at LFEs 1, 2, and 3 during a large episodic tremor and slip (ETS) event follows a previously observed evolution, in which tidal response during the initial 1 to 2 days of slip at an LFE patch is much weaker than that during the next few days of slip. The brief bursts of activity at LFE4, although little affected by large episodic tremor and slip events, respond strongly to tidal stresses. Our observed along‐dip contrasts in recurrence behavior and tidal response provide evidence for stress transfer and increased frictional strength with position updip through the tremor zone.

Crustal earthquakes in the Cook Inlet and Susitna region of southern Alaska

Several large (M ≥ 6) earthquakes have occurred in the vicinity of Anchorage, Alaska, within the past century. The presence of the underlying subducting Pacific plate makes it difficult to determine the origin of these older earthquakes as either crustal, slab, or the subduction plate interface. We perform a seismological study of historical and modern earthquakes within the Cook Inlet and Susitna region, west of Anchorage. We first estimate hypocenters for historical large earthquakes in order to assess their likelihood of origin as crustal, slab, or plate interface. We then examine modern crustal seismicity to better understand the style of faulting and the location of active structures, including within (and beneath) the Cook Inlet and Susitna basins. We perform double-couple moment tensor inversions using high frequency body waves (1–10 Hz) for small to moderate (M ≥ 2.5) crustal earthquakes (depth ≤ 30 km) occurring from 2007 to 2017. Our misfit function combines both waveforms differences as well as first-motion polarities in order to obtain reliable moment tensor solutions. The three focus regions—Beluga, upper Cook Inlet, and Susitna—exhibit predominantly thrust mechanisms for crustal earthquakes, indicating an overall compressive regime within the crust that is approximately consistent with the direction of plate convergence. Mechanisms within upper Cook Inlet have strike directions aligned with active anticlines previously identified in Cook Inlet from active-source seismic data. Our catalog of moment tensors is helpful for identifying and characterizing subsurface faults from seismic lineaments and from faults inferred from subsurface images from active-source seismic data.

Arcuate stress state in accretionary prisms from real-scale numerical sandbox experiments

The stress states in accretionary prisms are important for understanding the building and releasing of seismic energy. Numerous researchers have conducted sandbox experiments as a scaled physical analog model to understand the formation of accretionary prisms. However, measuring stress states in laboratory sandbox experiments is still practically infeasible. Here we performed real-scale numerical sandbox experiments using the discrete element method to understand the 3D stress state in the accretionary prism. Despite the nearly uniform initial conditions, macro-scale undulations of faults, which are similar to those observed in the trenches of an accretionary prism, appear. We reveal that these undulations are caused by the formation of stress arches. We show that the mechanism behind the arch formation is the discontinuous change in the stress orientation during the rearrangement of the stress chain. Furthermore, analyses demonstrate that the in-situ stress orientation from borehole data can be a signal of either the regional direction of plate convergence or the local stress orientation associated with the stress arch. The results may greatly enhance the outcome of long term monitoring in areas, such as the Nankai Trough.

Crustal stress field perturbations in the continental margin around the Korean Peninsula and Japanese islands

Seismic activity and focal mechanisms are governed by the effective stress field that is a combined result of regional tectonic processes and local stress perturbation. This study investigates the regional variation in the stress field in the eastern continental margin of the Eurasian plate around the Korean Peninsula and Japanese islands using a damped stress inversion technique based on the focal mechanism solutions of regional earthquakes. The dominant compressional stress is directed ENE–WSW around the Korean Peninsula and eastern China, E–W at the central East Sea and northern and southern Japan, NW–SE at central Japan, and N–S around the northern Nankai trough. The dominant compression direction changes rapidly in the East Sea and Japanese islands, which may be due to the combined effects of tectonic loading in the subduction zones off the Japanese islands and the India-Eurasia plate boundary. The crustal stress fields around the subduction zones off the Japanese islands present characteristic depth-dependent orientations. The orientations of the largest horizontal stress components, σH, in the subduction zones are subparallel with the plate convergence directions at shallow depths. The σH orientations are observed to rotate clockwise with the depth owing to slab subduction and lithospheric deformation. The regional stress field around the Japanese islands was perturbed temporally by the 2011 M9.0 Tohoku-Oki megathrust earthquake. The regional stress field was recovered in a couple of years. The stress field and tectonic structures are mutually affected, causing stress field distortion and a localized mixture of earthquakes in different faulting types.

High‐Precision Analysis of an Aftershock Sequence Using Matched‐Filter Detection: The 4 May 2015ML 6 Wanaka Earthquake, Southern Alps, New Zealand

ABSTRACT The 4 May 2015 M L 6.0 Wanaka earthquake was the largest earthquake in the Central Otago region of New Zealand’s Southern Alps for 70 years and the largest shallow event in the region since national records began. The earthquake occurred within 50 km of the alpine resort towns of Wanaka and Queenstown (combined population 35,000), at 9.1 km depth, and was felt widely across South Island, although no significant damage was reported. We utilize waveforms recorded on a temporary broadband network and the national network to conduct matched‐filter analysis of the aftershock sequence. We use the mainshock and 99 aftershocks distributed through the sequence as template events for the detection of further aftershocks. We detect 2544 aftershocks over 26 days, 27 times more than reported by the national network alone. Magnitudes of detected aftershocks range from M L 5.8 (34 s after the mainshock) to 0.3 and follow a Gutenberg–Richter exponential relationship with a b ‐value of 0.807±0.02 and a magnitude of completeness of M L 1.4. We compute high‐precision relocations of detected events using cross‐correlation‐derived pick corrections. Our relocated aftershocks highlight a steeply northwest dipping (∼70°) fault striking at 250°, which allows us to identify the structure responsible for the earthquake in a region where no active fault traces are mapped. This orientation is consistent with the mainshock focal mechanism for a dextral strike‐slip earthquake on a subvertical plane subparallel to the plate convergence direction (249°). Analysis of focal mechanism slip vectors shows that aftershocks principally occurred on secondary synthetic Riedel shears oriented 25° clockwise from the principal slip zone, possibly controlled by local structures. Aftershock locations migrate in the down‐dip direction away from the mainshock, slowing with a 1/log(time) relationship, consistent with previous numerical models of afterslip within the fault plane.

New insights on the tectonics of the Lampedusa Plateau from the integration of offshore, onland and space geodetic data

Joint analysis of multichannel seismic reflection profiles calibrated with well-logs across the northern part of the Lampedusa Plateau (central sector of the Pelagian Block, Sicily Channel), of structural data collected on Lampedusa island, and of GNSS geodetic velocities of sites on the islands and on the northern shore of the Channel, suggests that this part of the plateau forms an anticlinorium (Lampedusa Plateau Anticlinorium, LPA). The LPA developed during Paleogene to Early Miocene intraplate contraction followed by Miocene to current strike-slip deformation. It is formed by WNW-ESE striking highs and lows, which have an ~20 km average wavelength and culminate at the Lampione-Lampedusa High. These broad folds are bounded by high-angle faults with a reverse component of displacement, which cut Eocene to Lower Pliocene strata offshore, and Late Miocene strata on Lampedusa. Extensional faults, that have a bathymetric expression and are responsible for marked stratal tilting due to their listric geometry, are only found to the NE of the island and are associated to the rifting that affected the central part of the Sicily Channel in the Pliocene-Quaternary. Seismic reflection profiles show that normal fault activity peaked during the middle part of the Pliocene and strongly diminished afterward. Appraisal of recent plate motion reconstructions and of published and new structural data offshore and on-land suggest that the main growth phase of the LPA occurred during (Late Cretaceous?) Paleocene-Early Miocene ~N-S convergence between Nubia and Eurasia and associated intraplate shortening. Starting from Early Miocene, likely in response to a CCW rotation of the plate convergence direction, strike-slip deformation occurred with a ~NW-SE shortening axis and ~NE-SW extension axis. During this time span the previous contractional structures were locally reactivated in transpression. The two different strain regimes, extensional and transpressional that established since Miocene NE and W to NW of Lampedusa, respectively, still persist today as documented by geodetic velocities.

A comprehensive analysis of the Illapel 2015 Mw8.3 earthquake from GPS and InSAR data

The September 16, 2015 Mw8.3 Illapel Earthquake occurred on a locked segment of the South American subduction in Chile. This segment ruptured during comparable size earthquakes in the past, in 1880 and 1943, suggesting a somehow regular pattern of characteristic Mw8+ earthquakes occurring every 60 to 80 yr. This recurrence is in agreement with the accumulation of elastic deformation in the upper plate due to the Nazca–South America subduction at a constant rate of 6.5 cm/yr, leading to a deficit of ∼4.5 m of slip to be released every 70 yr. Previous studies consistently imaged the distribution of co-seismic slip along the fault based on geodetic, seismological and far field tsunami data and all described a significant amount of shallow slip resulting in a large tsunami. In addition, some models highlighted an apparent mismatch between the modeled rake of slip and the direction of plate convergence, suggesting the buildup of large strike-slip deficit. Some of these important questions remain open. Is shallow slip really well resolved and substantiated? Is the apparent principal direction of slip during the earthquake really required by the geodetic data?Here, using a comprehensive analysis of continuous GPS sites (including high rate and static displacements) and new survey data from acquired over more than 50 pre-existing sites, complemented with InSAR data, we show that the 2015 rupture overlaps very well the 1943 rupture, with the absence of significant slip south of 32°S and north of 30.2°S (peninsula Lingua de Vaca). Despite the wealth of geodetic data, the shallowest part of the subduction interface remains poorly resolved. We also show that the rake of the earthquake is fully compatible with the oblique plate convergence direction (rather than perpendicular to the trench), meaning that no subsequent trench-parallel motion is required by the data. We propose that the large Low Coupling Zone (LCZ) at the latitude of La Serena revealed by present day coupling distribution is stable over at least two seismic cycles. Inside the coupled area, peak coseismic slip is located precisely offshore the highest coastal topography and elevated terraces, adding weight to a potential correlation between the seismic cycle and long term permanent deformation. Finally, we show that early post-seismic after-slip occurs mostly down-dip of co-seismic asperity(ies), extending north and south of the 2015 rupture area.

The 16 April 2016, MW 7.8 (MS 7.5) Ecuador earthquake: A quasi-repeat of the 1942 MS 7.5 earthquake and partial re-rupture of the 1906 MS 8.6 Colombia–Ecuador earthquake

The 2016 Ecuador MW7.8 earthquake ruptured the subduction zone boundary between the Nazca and South American plates. Joint modeling of seismic and tsunami observations indicates an ∼120 km long rupture area beneath the coastline north of the 1998 MW7.2 rupture. The slip distribution reveals two discrete asperities near the hypocenter and around the equator. Their locations and the patchy pattern are consistent with the prior interseismic geodetic strain, which showed highly locked patches also beneath the coastline. Aftershocks cluster along two streaks, one aligned nearly parallel to the plate convergence direction up-dip of the main slip patches, and the other on a trench-perpendicular lineation south of the 1958 rupture zone. Comparisons of seismic waveforms and magnitudes show that the 2016 event and 1942 earthquakes have similar surface wave magnitude (MS7.5), overlapping rupture areas, and similar main pulses of moment rate. The same area ruptured as the southernmost portion of the larger earthquake of 1906 (MW8.6, MS8.6). The seismic behavior reflects persistent heterogeneous frictional properties of the Colombia–Ecuador megathrust.

Comparative study of tectonic tremor locations: Characterization of slow earthquakes in Guerrero, Mexico

AbstractDeep tectonic tremor in Guerrero, Mexico, has been observed using dense temporal seismic networks (i.e., the Meso‐American Subduction Experiment and Guerrero Gap Experiment (G‐GAP) arrays) during two different time periods. We apply a set of seismic waveform analysis methods to these data sets to constrain the locations of tremors and determine the associated moment tensors. First we detect and locate the tremors. Next, very low frequency (VLF) signals are identified by stacking waveform data during tremor bursts, and their moment tensors are determined. Finally, to better investigate the link between tremors and VLF earthquakes, we detect VLF events using a matched filtering algorithm to search continuous seismic records. None of the 11 VLF events detected by this method occurred in the absence of tremor bursts suggesting they are indeed part of the same phenomena. Unlike previous investigations, our results for the G‐GAP period reveal that downdip tremor activity (i.e., in the so‐called “sweet spot”) is segmented into two patches separated by 40 km in the along‐trench direction, indicating possible variations in the geometry of the plate interface and/or slab effective pressure. Moment tensors of VLF signals are consistent with shear slip on the near‐horizontal plate interface, but source depths are about 5 km deeper than the established plate interface. The slip directions of the VLF events are slightly (~10°) counterclockwise of the plate convergence direction, indicating that strain energy promoting left‐lateral strike‐slip motion may accumulate in the continental crust during the interseismic period.

The M w6.7 12 October 2013 western Hellenic Arc main shock and its aftershock sequence: implications for the slab properties

The 12 October 2013 M w6.7 earthquake offshore Crete Island is one of the few strong earthquakes to have occurred in the last few decades in the southwestern part of the Hellenic subduction zone (HSZ), providing the opportunity to evaluate characteristics of the descending slab. The HSZ has experienced several strong (M ≥ 7.0) earthquakes in historical times with the largest one being the 365 AD, M w = 8.4 earthquake, the largest known ever occurred in the Mediterranean region. The 2013 main shock occurred in close proximity with the 365 event, on an interplate thrust fault at a depth of 26 km, onto the coupled part of the overriding and descending plates. GCMT solution shows a slightly oblique (rake = 130°) thrust faulting with downdip compression on a nearly horizontal (dip = 3°) northeast-dipping fault plane with strike (340°) parallel to the subduction front, with the compression axis being oriented in the direction of plate convergence. The subduction interface can be more clearly resolved with the integration of aftershock locations and CMT solution. For this scope, the aftershocks were relocated after obtaining a v p/v s ratio equal to 1.76, a one-dimensional velocity model and time delays that approximate the velocity structure of the study area, and the employment of double-difference technique for both phase pick data and cross-correlation differential times. The first-day relocated seismicity, alike aftershocks in the first 2 months, shows activation of an area at the upper part of the descending slab, with most activity being concentrated between 13 and 27 km, where the main shock is also encompassed. Aftershocks are rare near to the main shock, implying homogeneous slip on a large patch of the rupture plane. Based on the aftershock distribution, the size of the activated area estimated is about 24 km long and 17 km wide. Coulomb stress changes resolved for transpressive motion reveal negligible off-fault aftershock triggering, evidencing a comparatively stable regime in the downdip part of the slab or different fault mechanism.

Far-field effects of Appalachian orogenesis: A view from the craton

The sedimentary cover of the North American craton preserved little evidence of the Paleozoic tectonic events that shaped the Appalachian orogen on its eastern side. A notable exception is the NNW-trending Hudson Bay central high, which corresponds to a normal-fault array extending for a minimum length of 500 km. A working hypothesis is proposed in which stresses applied to the continental margin during the Silurian–earliest Devonian Salinian orogeny were transmitted over a distance of >1400 km in the continental interior, where they induced the normal-fault reactivation of older structural discontinuities. The shutdown of tectonic activity along the Hudson Bay central high during the latest Early Devonian to earliest Middle Devonian is interpreted as resulting from a change in the direction of plate convergence during the Acadian orogeny.