Subduction Thrust Research Articles

Fluid upwelling across the Hikurangi subduction thrust during deep slow-slip earthquakes

AbstractSlow-slip events at global subduction zones relieve tectonic stress over days to years. Through slow-slip cycles, high fluid pressures observed at the top of subducting plates are thought to fluctuate, potentially due to the valving action of an impermeable layer near the plate interface. We model teleseismic scattering data at the Manawatu deep slow-slip patch at the Hikurangi margin in New Zealand and find high seismic P-to-S wave velocity ratios, VP/VS, in the upper ~5 km of the subducting Pacific Plate, reflecting sustained elevated fluid pressures that decrease during slow-slip and increase during inter-slow-slip periods. Within a ~ 3 km thick lower crustal layer of the overriding Australian Plate, decreasing VP/VS during inter-slow-slip periods reflects permeability reduction due to mineral precipitation. Increasing VP/VS during slow-slip reflects increasing permeability and crack density, facilitating upward fluid transfer through this layer. Our results suggest it acts as a valve to relieve high fluid pressures in the subducting slab.

Prediction of Fourier Amplitude Spectrum of Ground Motion in Mexico City from Subduction Thrust Earthquakes

Mexico City, one of the largest cities in the world, experiences frequent, devastating interplate earthquakes that originate on the subduction thrust along the Pacific coast of Mexico more than 250 km away. A notorious example is the 19 September 1985 Michoacán earthquake (Mw 8.0) which killed ~ 10,000 persons and caused wide-spread destruction in the city. The main cause of seismic damage in Mexico City is its subsoil that amplifies the ground motion. In view of the seismic hazard faced by the city, a reliable estimation of ground motion during future earthquakes is of vital importance. We present a new ground motion prediction equation (GMPE) for the Fourier amplitude acceleration spectrum (FAS) from subduction thrust earthquakes along the Pacific coast of Mexico at the reference station, CU, located in the firm zone of Mexico City. The GMPE is derived via Bayesian regression analysis and is based on an enlarged set of recordings (1965-2020; 40 earthquakes; 250 ≤ R≤ 500 km; 5 ≤ Mw ≤ 8.0). An important feature of the new GMPE is that it includes a term to model different attenuation along different ray paths. The inclusion of this term leads to a reduction in the aleatory variability of the GMPE, particularly at high frequencies. Since spectral amplification of seismic waves with respect to CU is known at many sites in the city, the FAS at these sites can be computed if it is known at CU. The new GMPE has been incorporated in a fully Fourier-based probabilistic seismic hazard analysis of Mexico City.

Repeating Large Earthquakes along the Mexican Subduction Zone

Abstract Repeating large earthquakes (M ≥ 7), waveforms for which are nearly identical, have been identified only on the Mexican subduction thrust near Acapulco. These earthquakes occurred on 1962 (Ms 7.0) and 2021 (Ms 7.0, Mw 7.0). Here, we report on two more sequences of three repeating large earthquakes each in eastern and western Oaxaca, Mexico. The repeating earthquakes in eastern Oaxaca occurred on 23 March 1928 (Ms 7.5), 1965 (Ms 7.6, Mw 7.5), and 2020 (Ms 7.4, Mw 7.4), and in western Oaxaca on 4 August 1928 (Ms 7.4), 1968 (Ms 7.2, Mw 7.3), and 2018 (Ms 7.2, Mw 7.2). Galitzin seismograms of the earthquakes in each sequence at DeBilt, The Netherlands or at Strasbourg, France are strikingly similar for at least 2600 s after the P-wave arrival. Similarity of waveforms of earthquakes in each sequence and tests with seismograms of events locations for which are accurately known suggest that their source areas were less than 10–20 km of each other. Moment-rate functions of these events are remarkably simple. We also document quasi-repeating earthquakes in central Oaxaca on 17 June 1928 (Ms 7.6) and 29 November 1978 (Ms 7.6, Mw 7.6). Such events have similar locations with large overlap in primary slip but are not identical. Recently, Michoacán–Colima earthquakes of 1973 (Ms 7.5, Mw 7.6) and 2022 (Ms 7.6, Mw 7.6) were reported as quasi-repeaters. Repeating or quasi-repeating large earthquakes imply that they are known for all the other events in the sequence if we know the location and gross source parameters of one of them. This permits the estimation of recurrence periods and the delineation of seismic gaps with greater confidence. Repeating and quasi-repeating large earthquakes in Oaxaca, an unique observation, shed new light on seismic hazard of the region, provide further support for the characteristic earthquake model, and reveal remarkably persistent behavior of ruptures through multiple earthquake cycles.

Structural map of Sulawesi derives from gravity data and its implications for geothermal systems

Gravity satellites are one of the alternatives to conducting preliminary surveys because the cost of exploring a large area using conventional gravity methods is very expensive and takes a lot of time. Therefore, this study on Sulawesi Island uses a gravity satellite to identify fault contact boundaries. The methodology in this study consists of several steps, starting with a literature study, data collection through ICGEM, and then processing and analysis by compiling the gravity anomaly with the Sulawesi regional geological structure map, further interpreted using the classification of rock density variations from a geological perspective. The regional gravity anomaly values for Sulawesi Island and its surroundings range from -270.9 mgal to 320.8 mgal, averaging 78.99 mgal. High anomalies are estimated to be between 180 mgal and 360 mgal, dominated by land with a lithology predominance of igneous rocks, while low anomaly areas are identified in the Batui thrust subduction zone. This study concludes that the convergence of high and low gravity anomalies characterizes the fault contact boundary in the study area. This has implications for the geothermal system on Sulawesi Island, which has a significant impact as a medium for the release of hot fluids from the subsurface.

Frictional Behavior Downdip Along the Subduction Megathrust: Insights From Laboratory Experiments on Exhumed Samples at In Situ Conditions

AbstractSubduction megathrusts exhibit a spectrum of slip modes, including catastrophic earthquakes. Although the mechanical and frictional properties of materials sampled from subduction zones have been studied extensively, few datasets have been collected for compositions and at pressure and temperature conditions representative of those in situ. The Nankai subduction zone in southwest Japan is a well‐studied margin, and abundant data provide an opportunity to advance our understanding of fault and earthquake physics. Here, we use samples exhumed in the Shimanto and Sanbagawa Belts on Shikoku Island of southwest Japan that represent analogs for materials along the present‐day megathrust at depths of ∼5–>25 km, and we shear these at their peak in situ pressure‐temperature (P‐T) conditions. Effective normal stresses range from 28 to 192 MPa, and temperatures from 105°C to 470°C. We used pore fluid pressures of 45–240 MPa, corresponding to fluid overpressure ratios λ of 0.65 and 0.90. Slip velocities of 0.1–100 μm/s were used, in order to focus on the nucleation of instability and earthquakes. We found predominantly velocity‐strengthening (inherently stable) behavior under all conditions for λ = 0.65. For λ = 0.90, velocity‐weakening behavior was observed at 350°C, with velocity‐strengthening behavior at lower and higher temperatures. The rate/state frictional stability parameter (a‐b) increases with slip velocity at temperatures up to ∼200°C and remains constant or decreases with slip velocity at higher temperatures. Overall, our results demonstrate the potentially important roles of both temperature and slip velocity in controlling the distribution of stress and frictional rheology along subduction thrusts.

A Source Study of the Mw 7.0 Acapulco, Mexico, Earthquake of 8 September 2021

AbstractThe Acapulco earthquake of 2021 broke a segment of the southeast Guerrero seismic gap along the Mexican subduction thrust. The rupture initiated offshore Acapulco (16.770° N, 99.942° W) and propagated down-dip toward northeast. This source directivity is confirmed from both (1) an analysis of local and regional recordings as a function of azimuth and (2) kinematic inversion of near-source, band-pass filtered (0.025–0.5 Hz) displacement seismograms and Global Positioning System static coseismic displacement vectors. The inversion reveals little slip near the hypocenter (<0.5 m) and significant slip distributed over an area of ∼184 km2, with the large slip patches in the northeast part of the fault. The estimated average slip and static stress drop are 260 cm and 18.6 MPa, respectively. Moment rate function reported by National Earthquake Information Center–U.S. Geological Survey from finite-fault modeling is simple, and it resembles other Mexican subduction earthquakes in the 7.0 ≤ M ≤ 7.5 range. Moment rate spectrum is well fit by the Brune ω−2 source model. Radiated seismic energy from teleseismic P waves is 7.5×1014 J, and ER/M0 is 2.1×10−5. Radiated energy enhancement factor—a measure of source complexity—is small, 5.8, similar to other Mexican subduction thrust earthquakes. Seismograms at DeBilt of the 2021 and the 11 May 1962 Acapulco earthquakes show an extraordinary similarity, seldom seen at M 7.0 level. The 2021 earthquake seems a repeat of the 1962 earthquake. The slip deficit since 1962 corresponding to a plate convergence rate of 6.2 cm/yr and perfect coupling is 366 cm. Thus, the seismic slip of 260 cm during that 2021 earthquake suggests a coupling ratio of 0.7, greater than 0.3 and 0.5 reported from geodetic measurements. Large moment release in the southeast seismic gap appears to have a periodicity of ~60 yr. Because 60 yr have elapsed since the last sequence earthquakes (1957 MS 7.5; 1962 MS 7.0 and 6.8), a renewal of large earthquakes in the region may be expected.

Embrittlement Within Viscous Shear Zones Across the Base of the Subduction Thrust Seismogenic Zone

AbstractGeophysical observations indicate that patches of localized fracturing occur within otherwise viscous regions of subduction plate boundaries. These observations place uncertainty on the possible down‐dip extent of the seismogenic zone, and as a result the maximum magnitude of subduction thrust earthquakes. However, the processes controlling where and how localized fracturing occurs within otherwise viscous shear zones are unclear. We examined three exposures of exhumed plate boundary on Kyushu, Japan, which contain subducted sediments and hydrated oceanic crust deformed at ∼300 to ∼500°C. These exposures preserve subduction‐related viscous deformation, which in two of the studied exposures has a mutually overprinting relationship with quartz veins, indicating localized cyclical embrittlement. Where observed, fractures are commonly near lithological contacts that form viscosity contrasts. Mineral equilibrium calculations for a metabasalt composition indicate that exposures showing cyclical embrittlement deformed at pressure‐temperature conditions near dehydration reactions that consume prehnite and chlorite. In contrast, dominantly viscous deformation occurred at intervening pressure‐temperature conditions. We infer that at conditions close to metamorphic dehydration reactions, only small stress perturbations are required for transient embrittlement, driven by localized dehydration reactions reducing effective stress, and/or locally increased shear stresses along rheological contrasts. Our results show that the protolith composition of the subducting oceanic lithosphere controls the locations and magnitudes of dehydration reactions, and the viscosity of metamorphosed oceanic crust. Therefore, compositional variations might drive substantial variations in slip style.

Aseismic slip and recent ruptures of persistent asperities along the Alaska-Aleutian subduction zone

The frictional properties and slip behaviors of subduction thrusts play a key role in seismic and tsunami hazard assessment, especially in weakly coupled “seismic gaps”. Here, we rely on GPS observations in the Shumagin Gap of the Aleutian subduction zone to derive the slip distribution of the 2020 Mw 7.8 Simeonof Island, Alaska earthquake and of the subsequent afterslip during the first 87-day period. Our modeling results show that the mainshock ruptured at depths of ∼30–40 km beneath Simeonof Island. Kinematic and stress-driven models indicate that the afterslip occurred both updip and downdip of the mainshock rupture. Physically plausible locking models derived from interseismic GPS velocities suggest that the 2020 Simeonof and 2021 Mw 8.2 Chignik earthquakes ruptured persistent asperities on the subduction thrust. We infer that there are several additional persistent asperities at depths of 20–50 km west ∼157°W. However, it is still uncertain whether there are additional locked asperities at shallow depths because of the current lack of geodetic observations close to the trench.

Spatial Variation of Shallow Stress Orientation Along the Hikurangi Subduction Margin: Insights From In‐Situ Borehole Image Logging

AbstractKnowledge of the contemporary in‐situ stress orientations in the Earth's crust can improve our understanding of active crustal deformation, geodynamic processes, and seismicity in tectonically active regions such as the Hikurangi Subduction Margin (HSM), New Zealand. The HSM subduction interface is characterized by varying slip behavior along strike, which may be a manifestation of combined variations in both stress state and the mechanical properties of faults and their hanging walls. Alternatively, these variations in subduction thrust slip behavior may drive heterogeneity in the stress state in space and time. In this study, we analyze borehole image and oriented four‐arm caliper logs acquired from 13 boreholes along the HSM to present a comprehensive stress orientation data set from borehole data within the shallow (<3 km) upper plate of the subduction thrust. Our results reveal a 065°/245° SHmax orientation within the central HSM (Hawke's Bay region) which rotates to 112°/292° and 140°/320° in the southern HSM. This rotation of SHmax orientation correlates spatially with along‐strike variations in subduction interface slip behavior, characterized by creep and/or shallow episodic slip events in the central HSM and interseismic locking in the south. The borehole SHmax orientations suggest that contemporary stress orientations may be caused in part by along‐strike variation in deformation style imposed by clockwise rotation of forearc. In the southern HSM, borehole‐derived SHmax orientations are inconsistent with SHmax orientations derived from focal mechanism solutions in the subducting plate, implying some degree of mechanical decoupling between the shallow hanging wall and subducting slab.

Heterogeneous material properties—as inferred from seismic attenuation—influenced multiple fault rupture and ductile creep of the KaikouraMw 7.8 earthquake, New Zealand

SUMMARYThe 2016 Mw 7.8 Kaikoura, New Zealand, earthquake occurred along the eastern margin of the transition region between active subduction in the North Island and oblique collision in the South Island. To infer crustal properties, we imaged Q (1/seismic attenuation) by combining selected M > 3.5 aftershocks with data from previous Q models. For 158 distributed aftershocks, we fit spectral decay on temporary stations and all Geonet stations, providing 6194 t*p and 19 497 t*s. Considering the varied rheology and faults, we also used 2.5-D numerical models to study ductile strain development. The complex earthquake ruptured an ∼180-km-long zone of multiple faults, which involved jumping around the complicated eastern end of the Hope fault, without significant slip on the Hope fault. The Qs and Qp results show features in the upper and lower crust which correlate to the distribution and types of fault rupture. This earthquake involved numerous faults over a region of greywacke crust, where the underlying high Q Cretaceous slab is about 30-km depth. It initiated with ∼5-m slip on the Humps fault in a region of background seismicity and low Q lower crust, adjacent to the Hope fault. The central region near Kaikoura shows a high Q crustal block, which appears to have inhibited rupture; as the rupture progressed over several small faults to jump offshore of the apparently strong block. Underlying the Kaikoura greywacke crust, below 20-km depth, there is a 40-km-long region of increased Vp, Qs and gravity, which likely represents an intraplate plutonic complex emplaced into the Hikurangi Plateau, forming an elevated section which influences deformation. In the northern section, in a region with relatively uniform moderate to low Q, the earthquake evolved into the relatively continuous ∼80-km-long major rupture along the Jordan, Kekerengu and Needles faults, with ∼6–20 m dextral slip at depth and surface displacements of ∼10 m dextral and ∼2 m vertical. The northern progression of the rupture stopped when it approached an abrupt change to high Q crust across Cook Strait. At 20–30-km depth northwest of the rupture, deeper zones with low Q are consistent with regions of distributed ductile shear and creep where the observed afterslip may have occurred, where the underlying slab is 25–40 km deep. The numerical model shows that ductile deformation localizes in this area of lower crust above the relatively strong slab, connecting outer faults (Kekerengu) to inland faults (Clarence, Awatere, Wairau), and demonstrates that no subduction thrust is required under the Marlborough region.

Pore Pressure Regime and Fluid Flow Processes in the Shallow Nankai Trough Subduction Zone Based on Experimental and Modeling Results from IODP Site C0023

AbstractPore fluid pressure is a critical parameter controlling the dynamics of subduction zones and related geological hazards, but quantitative constraints on pore pressure are sparse. Here we provide robust estimates of pore pressure in the shallow Nankai Trough subduction zone based on consolidation experiments and numerical models. Our experimental data reveal excess pore pressures with nearly constant ratios of pore pressure to lithostatic stress for the accreted and underthrust sediments, implying the décollement is not a fluid barrier as previously inferred. Also, our research suggests repeated lithostatic pore pressure transients in the last 3 kyrs, probably associated with the propagation of the frontal thrust. The coupled pore pressure dissipation and chemical transport model results let us conclude the updip migration fluid originates from an overpressured horizon nearby, confirming the view of a heterogeneous plate boundary with overpressured and permeable patches along the fault surface that shift in time and space. The resulting high excess pore pressure leads to the shear strength of Nankai subduction thrust as low as 0.4–0.8 MPa. This provides quantitative constraints on the weakness of subduction thrusts at the toe of accretionary prisms and may explain why seismic slip can propagate all the way to the trench in some earthquakes (e.g., 2011 Tohoku event).

Methoni Mw 6.8 rupture and aftershocks distribution from a dense array of OBS and land seismometers, offshore SW Hellenic subduction

Along the south-western offshore Hellenic subduction zone, the overriding Aegean upper plate above the Mediterranean oceanic lithosphere generates uncommon large earthquakes on the offshore megathrust fault. The largest subduction thrust event, for half a century, has been the 14 February 2008 Methoni earthquake (Mw = 6.8) that occurred offshore of the southwest coast of Peloponnesus. We conducted micro-seismicity experiments around the rupture area and forearc domain -between Peloponnesus and Crete- using ocean bottom seismometers (OBS) jointly with land-based seismological stations. Our first experiment in 2006, had revealed an association of the Matapan Trough, a 400-km-long forearc basin, with local seismicity clustering and a possible gap in activity over the later Methoni rupture area. Here we present new data of post-Methoni seismic activity, recorded during a time-span of 11 months, beginning in October 2008 within the period of proposed afterslip on the megathrust, by an extended and dense seismic array consisting of up to 33 OBS. A minimum 1D velocity model was constructed for the region to provide better constraints on absolute locations and double-difference relocation was applied to produce an enhanced image of the spatial distribution of hypocenters. The high resolution earthquake locations confirm correlation of the Matapan Trough with local seismicity as a regional feature, also filling up the previously observed gap. Over the Methoni rupture area, we constrain seismicity to be located mainly within the upper plate. Hypocenters are also resolved above the updip and downdip edges of the rupture area, respectively. Seismic activity provides hints of upper plate structures which were activated in response to post-seismic deformation spreading within the forearc crust. Our findings highlight the characteristics of a megathrust domain which is related with a highly deformable overriding plate and controlled by a segmented lower plate topography.

A Review of Tsunami Hazards in the Makran Subduction Zone

The uncertain tsunamigenic potential of the Makran Subduction Zone (MSZ) has made it an interesting natural laboratory for tsunami-related studies. This study aims to review the recent activities on tsunami hazard in the Makran subduction zone with a focus on deterministic and probabilistic tsunami hazard assessments. While almost all studies focused on tsunami hazard from the Makran subduction thrust, other local sources such as splay faults and landslides can be also real threats in the future. Far-field tsunami sources such as Sumatra-Andaman and Java subduction zones, commonly lumped as the Sunda subduction zone, do not seem to pose a serious risk to the Makran coastlines. The tsunamigenic potential of the western segment of the MSZ should not be underestimated considering the new evidence from geological studies and lessons from past tsunamis in the world. An overview of the results of tsunami hazard studies shows that the coastal area between Kereti to Ormara along the shoreline of Iran-Pakistan and the coastal segment between Muscat and Sur along Oman’s shoreline are the most hazardous areas. Uncertainties in studying tsunami hazard for the Makran region are large. We recommend that future studies mainly focus on the role of thick sediments, a better understanding of the plates interface geometry, the source mechanism and history of extreme-wave deposits, the contribution of other local tsunamigenic sources and vulnerability assessment for all coastlines of the whole Makran region.

Stress Orientations in the Nankai Trough Constrained Using Seismic and Aseismic Slip

AbstractThe strength of subduction thrust faults is key to understanding seismogenesis at the provenance of Earth's largest earthquakes. Earthquake focal mechanisms are routinely inverted to constrain the stress state at seismogenic depths. However, on some megathrusts, deformation is accommodated by both earthquakes and types of slow fault slip. We employ focal mechanisms of short‐term slow slip events (SSEs), a type of slow fault slip, and earthquakes in a regional stress inversion to investigate the stress state of the Nankai Trough megathrust and interpret the results in the context of regional tectonics Previous studies using earthquake‐only stress inversions found principal stress orientations in this region that are incompatible with thrust faulting on the megathrust. When both SSEs and earthquakes are considered, the stress state of the central and eastern Nankai Trough megathrust is well oriented for thrust faulting. Our results suggest that slow fault slip source regions may appear to have misoriented stress fields if slow fault slip constitutes a substantial proportion of fault slip and the stress field is not well constrained by earthquakes. In the SSE region, we find that faults are well oriented for failure, suggesting they have strengths similar to their surroundings. Combined with low Vp/Vs ratios and sensitivity to small stress changes, our results imply that the megathrust and surroundings operate at low deviatoric stresses in the SSE source region. Further, we show that the coefficient of friction for areas hosting SSEs is low (μ=0.19–0.50), implying frictionally weak materials in the SSE source region.

Stratigraphic architecture of Solander Basin records Southern Ocean currents and subduction initiation beneath southwest New Zealand

AbstractSolander Basin is characterized by subduction initiation at the Pacific‐Australia plate boundary, where high biological productivity is found at the northern edge of the Antarctic Circumpolar Current. Sedimentary architecture results from tectonic influences on accommodation space, sediment supply and ocean currents (via physiography); and climate influence on ocean currents and biological productivity. We present the first seismic‐stratigraphic analysis of Solander Basin based on high‐fold seismic‐reflection data (voyage MGL1803, SISIE). Solander Trough physiography formed by Eocene rifting, but basinal strata are mostly younger than ca. 17 Ma, when we infer Puysegur Ridge formed and sheltered Solander Basin from bottom currents, and mountain growth onshore increased sediment supply. Initial inversion on the Tauru Fault started at ca. 15 Ma, but reverse faulting from 12 to ca. 8 Ma on both the Tauru and Parara Faults was likely associated with reorganization and formation of the subduction thrust. The new seabed topography forced sediment pathways to become channelized at low points or antecedent gorges. Since 5 Ma, southern Puysegur Ridge and Fiordland mountains spread out towards the east and Solander Anticline grew in response to ongoing subduction and growth of a slab. Solander Basin had high sedimentation rates because (1) it is sheltered from bottom currents by Puysegur Ridge; and (2) it has a mountainous land area that supplies sediment to its northern end. Sedimentary architecture is asymmetric due to the Subtropical Front, which moves pelagic and hemi‐pelagic sediment, including dilute parts of gravity flows, eastward and accretes contourites to the shelf south of Stewart Island. Levees, scours, drifts and ridges of folded sediment characterize western Solander Basin, whereas hemi‐pelagic drape and secondary gravity flows are found east of the meandering axial Solander Channel. The high‐resolution record of climate and tectonics that Solander Basin contains may yield excellent sites for future scientific ocean drilling.

Strain partitioning across a subduction thrust fault near the deformation front of the Hikurangi subduction margin, New Zealand: A magnetic fabric study on IODP Expedition 375 Site U1518

Understanding the distribution of strain along thrust and splay faults in active accretionary systems is crucial to understand the mechanical properties of the sediments and the strength of the fault zone and its slip behavior. This paper investigates the distribution of strain through sediment compaction and texture development across the Pāpaku fault, a major splay fault near the deformation front of the Hikurangi subduction margin, New Zealand using the anisotropy of magnetic susceptibility technique (AMS). International Ocean Discovery Program Site U1518 penetrated hanging wall, fault zone and footwall sequences to a maximum depth of 484.9 meters below seafloor. A total of 330 discrete samples was subjected to AMS measurements and magnetic remanence data used to reconstruct the axial orientation of each sample in a geographic reference frame. The AMS display distinct fabric differences between hanging wall, through the fault zone and footwall domains, demonstrating that strain is partitioned across the fault zone. Hanging wall sequences show a strike-parallel northeasterly lineation of Kmax and weakly prolate shapes, typical for a component of northeast-southwest lateral shortening. In contrast, footwall sequences are more oblate and show a clustering of Kmax in northerly direction. This demonstrates that strain in the footwall is dominated by gravitational loading, however a component of sub-horizontal east-westerly strain, parallel to the convergence direction of the Pacific Plate exists. Strain decoupling between hanging- and footwall sequences occurs near the top of the Pāpaku fault zone. Differences in the degree of magnetic susceptibility between footwall sediments incorporated into the fault zone, and the underlying undeformed footwall sequences are indicative for the progressive dewatering of the underconsolidated footwall sequences.

Neotectonic structures and stress fields associated with oblique collision and forearc sliver formation in northern Hispaniola: Implications for the seismic hazard assessment

Active oblique collision between the Caribbean and North America plates has led to the formation of the Septentrional forearc sliver, which is a wedge-shaped crustal domain limited to the north by the North Hispaniola subduction thrust and to the south by the Septentrional Fault Zone (SFZ). In northern Hispaniola, oblique collision of the Bahamas Platform began in the uppermost Pliocene to Early Pleistocene and continues today. Collision-related tectono-stratigraphic processes include destruction of the shallow-water carbonate forearc platform, uplift and erosion of the Cordillera Septentrional and Samaná Peninsula, emergence of the Cibao Basin and syn-tectonic deposition of alluvial fan systems. The geometry and kinematics of neotectonics structures at all scales, together with stress field solutions obtained from late Cenozoic fault-slip data, show that the ENE-directed (far-) stress field associated with collision triggered the formation of the WNW-trending, left-lateral strike-slip SFZ. The SFZ thus represents the inland, concave master fault bounding the Septentrional forearc sliver. Subsidiary ENE to NE-trending left-lateral strike-slip faults (R-type) and NNE to NE-trending normal and oblique faults (T-type) controlled the internal deformation and subdivision into smaller tectonic blocks of the sliver. The highly oblique convergence and the increase in frictional plate coupling associated with the underthrusting of the Bahamas Platform are the likely causes of the formation of this forearc sliver. In the Western Cordillera Septentrional and Cibao Basin, an extensional regime is locally produced by the transfer of the left-lateral movement to a southern splay of the SFZ with the incorporation of a detached block to the Septentrional forearc sliver during the Middle Pleistocene-Holocene. The regional distribution of the modeled Peak Ground Acceleration values indicates that the SFZ is the greatest seismic hazard in northern Hispaniola.

Frictional Strengths of Subduction Thrust Rocks in the Region of Shallow Slow Earthquakes

AbstractEarthquakes are expected to nucleate within velocity‐weakening materials; however, at the updip limit of the subduction seismogenic zone, the principal lithologies exhibit velocity‐strengthening behavior. At an exhumed analogue for present‐day subduction at Nankai (the Mugi Mélange, Japan) two examples of paleoseismic features occur within altered basalts, suggesting that they may be velocity‐weakening. We shear altered basalt and shale matrix from the mélange in the triaxial saw cut configuration at in situ conditions of deformation (Pc = 120 MPa, T = 150 °C) for two pore fluid factors (λ = 0.36, 0.7). The shale matrix exhibits a coefficient of friction between 0.4 and 0.5, and velocity‐strengthening behavior. Altered basalt exhibits Byerlee friction and velocity‐weakening behavior. In most experiments deformation was partitioned into Riedel shears, which have a lower percentage of clasts, clast size distributions with higher fractal dimensions, and shape factors indicating rounder clasts compared to the matrix between Riedel shears. We hypothesize that earthquakes preferentially nucleate within altered basalt at the updip limit of the seismogenic zone. More complex forms of deformation (e.g., shallow very low frequency earthquakes) may occur by mixing velocity‐weakening altered basalt into the velocity‐strengthening shale matrix. In subduction zones where the matrix is composed of shale, this complex behavior is limited to shallow depths (T < ~250 °C), above the updip transition to velocity‐weakening behavior for the matrix. Altered basalt is a ubiquitous subducting material, and future studies on its behavior through a range of subduction zone conditions are required for a full understanding of the mechanical behavior of subduction zones.

Scaling of Peak Ground Displacement with Seismic Moment above the Mexican Subduction Thrust

AbstractWe use accelerograms, seismograms, and data from sparse continuous Global Positioning System (GPS) and campaign-mode GPS stations, deployed along the Pacific coast of Mexico, to study scaling of horizontal peak ground displacement (PGD) with seismic moment (M0) in the epicentral zone above the Mexican subduction thrust. The thrust interface is located at a depth of ∼25 km below the coast. We select recordings with (S–P) time ≤5.9 s (R≤46 km) and reduce the amplitudes to (S–P) time of 3.2 s (R = 25 km). The dataset consists of 58 events and covers a M0 range of 1013–1021 N·m. We find that the double integration of accelerograms, using piecewise linear detrending schemes, leads to sufficiently accurate estimation of PGD to study the scaling relation. The sparse data for great earthquakes are complemented with theoretical static displacement computed using the model of Okada (1992). For earthquakes with M0≤1.26×1018 N·m (Mw≤6.0) the point-source, far-field approximation holds, and the PGD data follows theoretically expected M02/3 scaling. For great earthquakes (M0>1.26×1021 N·m; Mw>8.0), static offset (which is approximately equivalent to PGD) scales as M01/3. About two-thirds of the observed PGD data fall within a factor of 0.67 and 1.5 of the relation given above. The relationship may be useful in earthquake engineering as well as for rapid estimation of magnitude for early tsunami alert.

Implications of basement rock alteration in the Nankai Trough, Japan for subduction megathrust slip behavior

The Nankai Trough is an exceptionally well-studied convergent margin known to host damaging megathrust earthquakes as well as various forms of slow fault slip. Outcrop studies of exhumed analogues for the modern subduction thrust, as well as seismic reflection images of the active subduction zone, suggest that the basaltic basement participates in plate-boundary shearing at depths of a few km and greater. We obtained altered basaltic upper basement from the modern Nankai Trough seaward of the trench, recovered by drilling on IODP Expedition 333. We performed laboratory friction experiments on bare surfaces and gouge powders of this material, in addition to and in combination with other materials for comparison. The altered basalt exhibits predominantly velocity-strengthening frictional behavior, indicating a tendency for stable slip. For bare surface experiments, few instances of velocity-weakening friction occur and are restricted to slip rates of <10−6 m/s. Thin sections and XRD analyses of the starting material indicate that the velocity-strengthening behavior is likely associated with the presence of clay minerals, mostly Mg-smectite, which coat the rims of larger, stronger grains. Based our friction data, we suggest that in addition to creep the altered Nankai basalts may also allow the possibility of slow slip events for two reasons: (1) velocity-weakening friction may be expected at low slip rates, but velocity-strengthening at higher rates will damp any potential slip instabilities; and (2) based on a critical stiffness criterion for accelerating (stable) slip, SSEs are capable of nucleating in the altered Nankai basalt despite velocity-strengthening friction. In light of these observations, and given the depth at which altered basement appears to be incorporated along the plate interface, we raise the possibility that some shallow slow slip events in the Nankai accretionary prism could originate from shearing in altered basalt.