Normal Slip Component Research Articles

Variations in the present-day stress field and its implications for hydrocarbon development in the Hassi Terfa field, Hassi Messaoud, Algeria

Understanding the contemporary stress field and its variations is important in hydrocarbon production. These fluctuations are significant for borehole stability and trajectory, hydraulic stimulation, fault reactivation, and optimal well completion. In this study, we provide the first geomechanical model to determine the orientation and magnitude of the present-day stress of the Hassi Terfa field using geophysical logs to understand the origin of the spatial variation of the stress state. The breakouts from acoustic image logs provide an average maximum horizontal stress (SHmax) orientation of N115°, consistent with the Africa–Eurasia plate motion, which imposes the largest first-order stress source. However, our analysis demonstrates local stress perturbations near faults. This may mean that the stress field in the Hassi Terfa is superimposed by a third-order source of stress (<100 km). Our model demonstrates variations in stress magnitude with depth interpreted to be lithology-dependent, which act as a natural barrier against fracture propagation. We propose a new approach to quantify these perturbations via the stress regime index R′, using the magnitude of the principal stress components. The results demonstrate different layer-to-layer faulting regimes. These perturbations will limit the vertical propagation of hydraulic fracturing to the overlying formation, which enables it to access a larger reservoir volume and improve production. In general, a strike–slip regime with a slight normal slip component is interpreted for the study area. We develop the first-ever quantitative stress map using the relative principal stress magnitude parameter (Aϕ) to present the lateral variation of the stress regime throughout the Saharan platform. Results demonstrate a consistent strike–slip to normal/strike–slip regime with a uniform SHmax orientation. This map will help operators achieve optimal well direction and optimize hydraulic stimulation. Our results have considerable implications for hydrocarbon production in the central and southeastern Saharan platform.

Geodetic constraints of the 2015Mw6.5 Alor, East Indonesia earthquake: a strike-slip faulting in the convergent boundary

SUMMARYOn 2015 November 4, an Mw 6.5 earthquake struck the east region of Alor Island, eastern Indonesia. Here, we use Sentinel-1 and ALOS-2 data to explore the coseismic surface displacement of this earthquake. Based on the ALOS-2 coseismic interferograms, the first fault model and coseismic slip distribution of the 2015 Alor earthquake are presented in this study. The preferred slip model links to a blind south-southeast striking, south-southwest dipping strike-slip fault with amount of normal slip component and a peak slip of 2.09 m at 2.34 km of depth. Considering the fact that the ascending and descending InSAR prediction of the distributed slip model does not fit the InSAR observations at the northwest and southeast tips with the simple planar fault model. We tried to construct a strike-variable fault model to further improve data fit. The results of the calculated Coulomb stress change imply that the regions with positive CFS changes mainly located at the northwest and southeast extremities of the rupture of the Alor earthquake, and in the lobes north and south of the rupture. The most striking discovery from the InSAR observations of the Alor earthquake is that most of the displacements occurred on a fault whose existence was unknown before the earthquake.

Field observations of surface ruptures accompanying a tsunami and supershear earthquake along a plate boundary strike-slip fault

AbstractStrike-slip earthquakes near major subduction zones have received less attention than thrust or reverse earthquakes in subduction zone areas. The occurrence of the 2018 Palu Mw 7.5 earthquake in eastern Indonesia provides an unprecedented opportunity to investigate the characteristics of one of these events. The Palu earthquake occurred on the left-lateral, north–south-striking Palu–Koro fault, which is the main plate boundary structure accommodating the convergence between blocks in a triple junction area. It excited a significant tsunami, which unusually is associated with strike-slip earthquakes, and also ruptured at a supershear speed, which is mostly observed on strike-slip faults in continents. Based on our fieldwork, we speculate that the normal slip component of the offshore rupture section in Palu bay on the middle segment probably favours tsunami genesis. Our field investigation has revealed evidence of a simple geometry as well as slip partitioning of dip-slip and strike-slip motion on two subparallel strands on the main segment, both of which may have contributed to the supershear of the rupture propagation. Instead of only a transtensive behaviour of the middle segment, our results also illustrate the transpressional property of the northern and southern rupture segments, which shows more complex behaviour than that of a common continental strike-slip fault.

Refining slip distribution in moderate earthquakes using Sentinel-1 burst overlap interferometry: a case study over 2020 May 15Mw 6.5 Monte Cristo Range Earthquake

SUMMARYAzimuth and range measurements from Synthetic Aperture Radar (SAR) images are commonly used to depict the coseismic deformation of large earthquakes. Azimuth measurements associated with moderate-sized earthquakes are hardly applicable due to their limited accuracy. In this paper, we first explored the capability of Sentinel-1 azimuth measurements to map the coseismic deformation of a moderate earthquake. We used both range and azimuth offset measurements to map the coseismic deformation of the 2020 Mw 6.5 Monte Cristo Range earthquake in Nevada. Optimal dip angles of the two main faults and the slip model were obtained. By adding azimuth displacements as constraints, the inverted slip model reveals shallower and more refined slip than models only constrained by InSAR and GPS data, highlighting the importance of introducing the azimuth measurements to the moderate earthquake. The preferred fault model shows a mixture of left-lateral and normal faulting on the western segment and a left-lateral slip on the eastern segment. We found that the inferred faults might verify the clockwise rotation block model, which has been proposed to accommodate fault slip across the Mina deflection. Moreover, a shallow alluvial basin and the young left-lateral fault within the left bend can potentially be responsible for the orientation and normal slip components in the western fault segment, respectively.

Fault Geometry and Slip Distribution of the 2021 Mw 7.4 Maduo, China, Earthquake Inferred from InSAR Measurements and Relocated Aftershocks

Abstract A nearly 70 yr hiatus of major seismic activity in the central eastern Bayan Har block (BKB) ended on 22 May 2021, when a multislip-peak sinistral strike-slip earthquake struck western Maduo County, Qinghai. This earthquake, which ruptured the nearly 170 km long Kunlun Pass–Jiangcuo fault, is a rather unique event and offers a rare opportunity to probe the mechanical properties of the intraplate lithosphere of the central eastern BKB. Here, we inferred the fault geometry associated with the Maduo earthquake using Interferometric Synthetic Aperture Radar (InSAR), and relocated aftershocks and inverted the slip distribution through InSAR radar phases and range offsets. Our analysis revealed that the geometry of the fault varies along the strike: the southeastern end of the fault dips steeply to the northeast, whereas the northwestern end dips southwestward. Using the combined datasets to constrain a coseismic slip, we found that the 2021 Maduo event was dominated by sinistral strike-slip movement, with a slight normal-slip component at a shallow depth, rupturing the steep-dipping fault for nearly 170 km in length. Five asperities were detected along the fault strike in the shallow crust (0–12 km) with a peak slip of ∼4.2 m corresponding mostly to simple structures, namely, continuous and straight rupture segments, suggesting that the rupture propagated across geometrical barriers in a multiasperity way. Based on an analysis of the strain field and the focal mechanisms of both the 2021 Maduo earthquake and historical earthquakes that have occurred in the BKB, we propose that the fault zones within the BKB can also generate large earthquakes and have the ability to accommodate the ongoing eastward and northeastward penetration of the Indian plate into the Eurasian plate.

Revisiting the Source Rupture Process of the Mainshock of the 2016 Kumamoto Earthquake and Implications for the Generation of Near-Fault Ground Motions and Forward-Directivity Pulse

ABSTRACT Strong near-fault ground motions associated with the MJMA 7.3 mainshock of the 2016 Kumamoto, central Kyushu, Japan, earthquake sequence have received attention by seismological and engineering communities. In this study, the kinematic source rupture process was reanalyzed based on an improved approach for the representation of source faults. The slips at densely distributed point sources were defined via the bilinear interpolation of those at surrounding control points. The result shows that the rupture started on the Hinagu fault with a small initial rupture and propagated beyond the junction to the Futagawa fault. The rupture on the Futagawa fault mainly propagated up and northeastward. A large slip area with a peak slip of 4.9 m and peak slip velocity of 3.1 m/s was detected at depths ranging from 3 to 15 km in the central part of the Futagawa fault. This asperity spatially coincides with a body with moderate-seismic velocity (VP∼6 km/s) and low-seismic attenuation. The slips on the Futagawa fault have significant normal-slip components, whereas the slip vectors of the Hinagu fault represent almost pure right-lateral strike slip. The shallower part of the fault segments in the western Aso caldera is characterized by relatively large normal slips. The estimated slip-velocity functions at shallower depths (&amp;lt;3 km) are almost temporally symmetric and relatively long. The shallower portion of the source fault significantly contributes to the velocity and displacement waveforms at near-fault sites. On the contrary, the slip-velocity functions at deeper depths (&amp;gt;3 km) are temporally asymmetric and have a sharp peak. The simulation of the ground-motion evolution suggests that the lateral flow in the Aso Valley was primarily triggered by the strong forward-directivity pulse generated from the asperity on the Futagawa fault.

Geodetic Observation and Modeling of the Coseismic and Postseismic Deformations Associated With the 2020 Mw 6.5 Monte Cristo Earthquake

AbstractThe 2020 Mw 6.5 Monte Cristo earthquake occurred in the northeastern Mina deflection of the central Walker Lane Belt (WLB) in Nevada, USA. The Mina deflection represents a typical stepover zone that transfers the dextral slip in the northern Eastern California Shear Zone onto the dextral faults in the central WLB. The Monte Cristo earthquake provides a rare opportunity to investigate the strain accumulation and stress transition mechanisms in the WLB. In this study, ascending and descending Sentinel‐1 images were utilized to generate coseismic and early postseismic deformations associated with this earthquake. Combined with global positioning system measurements, these images were inverted for the coseismic slip and afterslip of the Monte Cristo earthquake. The preferred coseismic slip model suggests that the causative fault is characterized by two fault segments with southward dips of 64° and 79°. The coseismic rupture was dominated by sinistral slip with obvious normal slip components in the western segment of the source fault. The coseismic slip was mainly concentrated in the 3–12 km depth range and decreased at shallower depths, suggesting a moderate amount of shallow coseismic slip deficit. Rapid afterslip was mainly confined to the 0–3 km depth range, largely compensating for the shallow slip deficit caused by the mainshock. The widely distributed normal slips in the northeastern Mina deflection revealed by the Monte Cristo earthquake suggest the transtensional model is more applicable due to its ability to account for the slip transition kinematics in the Mina deflection.

Multi-Segment Rupture Model of the 2016 Kumamoto Earthquake Revealed by InSAR and GPS Data

The 2016 Kumamoto earthquake, including two large (Mw ≥ 6.0) foreshocks and an Mw 7.0 mainshock, occurred in the Hinagu and Futagawa fault zones in the middle of Kyushu island, Japan. Here, we obtain the complex coseismic deformation field associated with this earthquake from Advanced Land Observation Satellite-2 (ALOS-2) and Sentinel-1A Interferometric Synthetic Aperture Radar (InSAR) data. These InSAR data, in combination with available Global Positioning System (GPS) data, are then used to determine an optimal four-segment fault geometry with the jRi method, which considers both data misfit and the perturbation error from data noise. Our preferred slip distribution model indicates that the rupture is dominated by right-lateral strike-slip, with a significant normal slip component. The largest asperity is located on the northern segment of the Futagawa fault, with a maximum slip of 5.6 m at a 5–6 km depth. The estimated shallow slips along the Futagawa fault and northern Hinagu fault are consistent with the displacements of surface ruptures from the field investigation, suggesting a shallow slip deficit. The total geodetic moment release is estimated to be 4.89 × 1019 Nm (Mw 7.09), which is slightly larger than seismological estimates. The calculated static Coulomb stress changes induced by the preferred slip distribution model cannot completely explain the spatial distribution of aftershocks. Sensitivity analysis of Coulomb stress change implies that aftershocks in the stress shadow area may be driven by aseismic creep or triggered by dynamic stress transfer, requiring further investigation.

Geometric Variation in the Surface Rupture of the 2018 Mw7.5 Palu Earthquake from Subpixel Optical Image Correlation

We obtained high-resolution (10 m) horizontal displacement fields from pre- and post-seismic Sentinel-2 optical images of the 2018 Mw7.5 Palu earthquake using subpixel image correlation. From these, we calculated the curl, divergence, and shear strain fields from the north-south (NS) and east-west (EW) displacement fields. Our results show that the surface rupture produced by the event was distributed within the Sulawesi neck (0.0974–0.6632°S) and Palu basin (0.8835–1.4206°S), and had a variable strike of 313.0–355.2° and strike slip of 2.00–6.62 m. The NS and EW displacement fields within the Palu basin included fine-scale displacements in both the near- and far-fault, the deformation patterns included a small restraining bend (localized shortening), a distributed rupture zone, and a major releasing bend (net extension) from the curl, divergence, and shear strain. Surface rupture was dominated by left-lateral strike-slip from initiation to termination, with a localized normal slip component peaking at ~3.75 m. The characteristics and geometric variation of the ruptured fault controlled both the formation of these surface deformation patterns and sustained supershear rupture.

Seismic potential of the Dead Sea Fault in the northern Gulf of Aqaba-Elat: New evidence from liquefaction, seismic reflection, and paleoseismic data

The cities of Elat, Israel and neighboring Aqaba, Jordan are major economic, cultural, and seaport centers. They are located on the northern shore of the Gulf of Aqaba/Elat (GAE) directly on the Dead Sea Transform. Yet the precise location of the fault trace and its tectonic activity are lacking. The interpretation of seismic reflection profiles across the GAE beach and paleoseismic trench data located 2.2 km north of the shoreline provide evidence that the active offshore mapped Avrona Fault extends onland along the eastern side of the Elat Sabkha (mudflat), where three prominent fault strands crosscut the sedimentary fill. Mismatch of reflector geometry across the faults and flower structures indicate strike-slip faulting with a normal-slip component. Subsurface data from two trenching sites provide evidence for a minimum of two surface ruptures and two paleoliquefaction events. Faulting is constrained by radiocarbon dating for an Event 1 between 897 and 992 CE and Event 2 after 1294 CE. We suggest that the historically documented 1068 CE, and at least one later earthquake in 1458 or 1588 CE, ruptured the Elat Sabkha site. Based on fault mapping, we suggest a minimum value of M 6.6 for the 1068 CE earthquake. Whereas no surface rupture was observed for the 1212 CE historical earthquake, fluidized strata radiocarbon dated to before 1269–1389 CE identified as paleoliquefaction may be attributed to it. Two liquefaction sand-blows mapped in the trench likely formed after 1337 CE and before 1550 CE, which possibly occurred at the same time as in the second faulting event. Our data suggest that no large event occurred along the Avrona segment in the past ~430–550 years. Given a ~ 5 mm/yr slip rate, we conclude that a significant period of time passed since the last surface rupturing on the Avrona Fault, increasing its seismic potential.

A hybrid source mechanism of the 2017 Mw 6.5 Jiuzhaigou earthquake revealed by the joint inversion of strong-motion, teleseismic and InSAR data

As occurred during the tourist season, the 2017 Mw 6.5 Jiuzhaigou (China) earthquake led to destructive damages. The seismogenic fault of this event was merely speculated to be the northern extension of the Huya fault, while no apparent surface ruptures were discovered in the field investigation. Previous studies and released moment tensor solutions indicated the Jiuzhaigou earthquake was dominated by left-lateral strike slip with partial normal-slip component. It seems unintelligible that the normal slip occurred in this event as the epicenter was located near the boundary of convergent blocks. Hence, a reasonable and elaborate source rupture model is necessary to investigate such a blind fault and the role it plays in the complex fault system of this region. Combining space-based geodetic, teleseismic or regional seismic observations can provide detailed information about earthquake ruptures. We first attempted to determine the fault geometry using the Bayesian approach with synthetic aperture radar interferograms (InSAR). And then the fault geometry was refined based on the relocated aftershock distribution, and a two-segment fault model was constructed. Based on the two-segment model, we resolved the source rupture process of the 2017 Jiuzhaigou earthquake through the joint inversion of strong-motion, teleseismic body-wave and InSAR data. The inversion results reveal a hybrid source mechanism, in which normal and thrust slips coexist besides the strike-slip component. We suggest that the eastward motion with the extrusion of the lower crustal flow in the northeastern margin of the Bayan Har block is responsible for such a faulting behavior. The co-seismic Coulomb stress changes show a significant stress loading in the western segment of the Tazang fault, increasing its seismic hazard. Due to the lack of aftershocks to the southeast of the seismogenic fault, the probably enhanced seismicity in the northern segment of the Huya fault is also worth further attention.

Hydrodynamics of Palu Bay during the Event of 2018 Palu Earthquake and Tsunami

The devastation of coastal area in Palu Bay few minutes after the September 28th, 2018 Sulawesi earthquake showed high variation of tsunami arrival time as well as the tsunami run-up and inundation. Recent findings showed that both local submarine landslides and the normal-slip components inside the Palu Bay may contribute to the generation of tsunami. However, the fact that the event occurred during high tide, the hydrodynamic characteristics of this narrow bay and their role in the dynamics of the generated of tsunami were unknown. Hydrodynamics simulation (Mike21-flow model) using the latest available bathymetry field data (the 2018 deep water of the Indonesian navy data and 2015 shallow water of the BIG data) was conducted to investigate the variation of sea levels and tidal currents within the bay during the event of earthquake and tsunami or within the first 8 minutes timeframe. Results showed that significant increase of water elevation up to 6 cm and current velocity up to 1 cm/s directed towards the city of Palu were observed that may contribute to the dynamics of the tsunami e.g. the speed of tsunami arrival time and the transformation of tsunami. Therefore, considering that multiple tsunami arrivals were in few minutes after the earthquakes, the hydrodynamics of Palu Bay during the event should also be considered in future tsunami simulation scenarios.

Source parameters and slip distribution of the 2018 Mw 7.5 Palu, Indonesia earthquake estimated from space-based geodesy

The 2018 Mw 7.5 Palu earthquake is one of the most geometrically complex events ever recorded on Central Sulawesi Island. It caused devastating tsunamis and landslides after the mainshock, resulting in great casualties and damages to the infrastructures in Palu city and surrounding areas. In this study, we combined far-field InSAR, near-field SAR and optical sub-pixel correlation data to determine the fault geometry, source parameters and fault-slip distribution of this event along the Palu-Koro fault zone. Our best-fitting fault model contains three fault segments with variable strike angles, showing that the coseismic rupture is dominated by a left-lateral strike-slip motion accompanied by minor normal-slip and reverse-slip components. On the onshore segments, the maximum slip is up to ∼ 8.4m with a left-lateral strike-slip component of ∼ 6.8m at a depth of 4km on the central segment near Palu city, a normal-slip component of ∼ 4.3m in the left releasing step-over and a reverse-slip component of ∼ 3.2m in the offshore area near the Palu bay. On the offshore segment, the maximum slip is up to ∼ 4.2m with a normal-slip component of ∼ 1.8m. We argue that the offshore reverse-slip motion might have played a key role in the generation of the Palu tsunami. We also suggest that the linking normal segment contributed to the rupture propagation across the left releasing step-over from the central segment to the southern segment, while the linking reverse segment arrested the rupture propagation across the right restraining step-over to north. The Coulomb stress change reveals that there is possibility of future great earthquakes occurring on several major faults in the northeast of Palu city. Our results have great significance to qualitatively assess the risk of future seismic hazards in such a geometrically complex fault system.

Seismogenic fault and tectonic significance of 1996 Karakoram Pass earthquake (Ms 7.1) based on InSAR

Due to hard observation condition of the western Tibet region, the slip behaviors of the Ms7.1 Karakoram Pass earthquake occurred in Hetian, Xinjiang on November 19, 1996 remains unclear. Using ERS 1/2 SAR data and InSAR technique, we obtain the co-seismic deformation of the earthquake. The north and south deformation areas show asymmetric pattern, with the maximum LOS displacement of the southern part approximately 24.6 cm, and the maximum LOS displacement in the northern part approximately − 18.5 cm. Nonlinear and linear inversion algorithms are used to determine the geometric parameters and slip distribution of the earthquake fault. Our results show that the co-seismic displacement is dominated by deformation fields are clearly visible sinistral strike-slip accompanied by a small amount of normal slip component. The co-seismic slip occurred between 0 and 18 km at depth. The maximum slip is ~ 81 cm, occurring at a depth of 8.5 ± 0.5 km at (35.36°N 78.03°E), indicating a shallow event with a moment magnitude of Mw 6.5. The seismogenic fault is a secondary fault in the Karakoram fault zone with strike 96°, dip 84°, and rake – 24°. This earthquake shows that the Karakoram fault zone undergoes a complex tectonic deformation process, with central part of the fault zone showing minor tensional deformation behaviors.

Complete three-dimensional near-field surface displacements from imaging geodesy techniques applied to the 2016 Kumamoto earthquake

The recent development of imaging geodesy, an advanced technique with a high spatial resolution and large-scale coverage, has enabled researchers to obtain multiple high-quality surface displacement estimates at low labor-cost, thereby improving the capability to monitor and manage geological disasters. The different sources (e.g., radar, optical and LiDAR sensors) and analysis approaches (e.g., differential interferometric synthetic aperture radar, DInSAR; multiple-aperture InSAR; pixel offset tracking; and iterative closest point, ICP) in imaging geodesy used to derive displacement estimates have unique benefits and drawbacks. However, the inherent differences among these data sources and methods in the construction of three-dimensional (3D) deformation maps, particularly in the near field, remain poorly understood and require further discussion. In this study, we acquired three pairs of ALOS-2 stripmap mode images, two pairs of Sentinel-1 TOPS mode images and pre- and post-event LiDAR data for the 2016 Kumamoto earthquake to explore the 3D near-field displacements using various imaging geodesy techniques with different types of image information, i.e., SAR phase data, SAR amplitude data and LiDAR point cloud data. Our results show that each image type is independently capable of producing a high-quality 3D deformation map for the 2016 Kumamoto earthquake with an on-fault accuracy of <43 cm determined from field work measurement, and an off-fault accuracy of <14 cm determined from GNSS observations. The major contributors to the uncertainty in the 3D deformation estimates from the SAR phase, SAR amplitude, and LiDAR data methods are the effective Doppler bandwidth, pixel resolution and topographic roughness, respectively. In the near-field deformation region, more secondary fault ruptures were revealed by SAR amplitude and LiDAR information than by the SAR phase information, thereby overcoming the sensitivity of the SAR phase signal to incoherence. Finally, an integrated complete 3D map was generated to constrain the coseismic rupture behavior of the Kumamoto earthquake sequence associated with the Futagawa and Hinagu faults. Our slip model suggests the main slip rupture a length of 32 km and terminated near the rim of the Aso caldera; additionally, the slip pattern delineated by three asperities that was dominated by right-lateral strike-slip with a minor normal slip component at a depth of 7 km. Furthermore, the main rupture triggered two secondary faults with right-lateral strike-slip and normal slip. Our source model yields a geodetic moment of 5.77 × 1019 N·m, which corresponds to Mw 7.11.

A NE-Trending Oblique-Slip Fault Responsible for the 2016 Zaduo Earthquake (Qinghai, China) Revealed by InSAR Data

Twelve Sentinel-1 images are used to measure the surface displacements caused by the Zaduo earthquake on 17 October 2016 in Qinghai, China. Comprehensive analysis of the local geologic structure, focal mechanisms, relocated aftershocks, and fault geometrical parameters from nonlinear inversions of the InSAR data indicates that a NE-trending fault with a strike of ~ 61°, namely the Aduo fault, was responsible for this event. Further geodetic inversions for finite fault models show that the dip angle of this fault is 75° to the SE. The majority of the slip occurred at the depth from 5 km to 17 km, and the normal-slip component was comparable to the strike-slip component. The peak slip of ~ 31 cm occurred at a depth of ~ 10 km. The InSAR determined moment is 9.61 × 1017 N m with a rigidity of 30 GPa, equivalent to Mw 5.9. The geological map reveals that the Aduo fault was characterized by a left-lateral strike-slip component in the past, as evidenced by the offset of Paleogene sedimentary rocks. Our results for the Zaduo earthquake suggest that the Qiangtang block probably consists of two or more sub-blocks with different eastward extrusion velocities and that the deformation of the Qiangtang block is occurring not only along its boundaries but also in its interior.

Simultaneous estimation of the dip angles and slip distribution on the faults of the 2016 Kumamoto earthquake through a weak nonlinear inversion of InSAR data

At the 2016 Kumamoto earthquake, surface ruptures were observed not only along the Futagawa fault, where main ruptures occurred, but also along the Hinagu fault. To estimate the slip distribution on these faults, we extend a method of nonlinear inversion analysis (Fukahata and Wright in Geophys J Int 173:353-364, 2008) to a two-fault system. With the method of Fukahata and Wright (2008), we can simultaneously determine the optimal dip angle of a fault and the slip distribution on it, based on Akaike’s Bayesian information criterion by regarding the dip angle as an hyperparameter. By inverting the InSAR data with the developed method, we obtain the dip angles of the Futagawa and Hinagu faults as 61° ± 6° and 74° ± 12°, respectively. The slip on the Futagawa fault is mainly strike slip. The largest slip on it is over 5 m around the center of the model fault (130.9° in longitude) with a significant normal slip component. The slip on the Futagawa fault quickly decreases to zero beyond the intersection with the Hinagu fault. On the other hand, the slip has a local peak just inside Aso caldera, which would be a cause of severe damage in this area. A relatively larger reverse fault slip component on a deeper part around the intersection with Aso caldera suggests that something complicated happened there. The slip on the Hinagu fault is almost a pure strike slip with a peak of about 2.4 m. The developed method is useful in clarifying the slip distribution, when a complicated rupture like the Kumamoto earthquake happens in a remote area.

Deformation style in the damage zone of the Mondy fault: GPR evidence (Tunka basin, southern East Siberia)

The Mondy strike-slip fault connects the W-E Tunka and N-S Hovsgol basins on the southern flank of the Baikal rift system. Ground penetrating radar (GPR) surveys in its damage zone provide constraints on thicknesses, dips, and plunges of fault planes, as well as on the amount and sense of vertical slip. Strike-slip faulting in the southern segment of the Mondy fault within the territory of Russia bears a normal slip component of motion along the W-E and NW planes. These motions have produced negative flower structures in shallow crust appearing as grabens upon Pleistocene fluvioglacial terraces. The amount of normal slip estimated from the displacement of reflection events varies over the area and reaches its maximum of 3.4 m near Mondy Village. In the Kharadaban basin link, left-lateral strike slip displaces valleys of ephemeral streams to 22 m, while normal slip detected by GPR reaches 2.2 m; this normal-to-strike slip ratio corresponds to a direction of ~6° to the horizon. The angles of dips of faults are in the range 75°-79°; the thicknesses of fault planes marked by low- or high-frequency anomalies in GPR records vary from 2.5 to 17.0 m along strike and decrease with depth within a few meters below the surface, which is common to near-surface coseismic motions. Many ruptures fail to reach the surface but appear rather as sinkholes localized mainly in fault hanging walls. The deformation style in the damage zone of the Mondy fault bears impact of the NW Yaminshin fault lying between its two segments. According to photoelasticity, the stress field changes locally at the intersection of the two faults, under NE compression at 38°, till the inverse orientations of principal compression and extension stresses. This stress pattern leads to a combination of normal and left-lateral strike slip components.

Source rupture processes of the foreshock and mainshock in the 2016 Kumamoto earthquake sequence estimated from the kinematic waveform inversion of strong motion data

The 2016 Kumamoto earthquake sequence started with an MJMA 6.5 foreshock occurring along the northern part of the Hinagu fault, central Kyushu, Japan, and the MJMA 7.3 mainshock occurred just 28 h after the foreshock. We analyzed the source rupture processes of the foreshock and mainshock by using the kinematic waveform inversion technique on strong motion data. The foreshock was characterized by right-lateral strike-slip occurring on a nearly vertical fault plane along the northern part of the Hinagu fault, and it had two large-slip areas: one near the hypocenter and another at a shallow depth. The rupture of the mainshock started from the deep portion of a northwest-dipping fault plane along the northern part of the Hinagu fault, then continued to transfer to the Futagawa fault. Most of the significant slip occurred on the Futagawa fault, and the shallow portion of the Hinagu fault also had a relatively large slip. The slip amount on the shallowest subfaults along the Futagawa fault was approximately 1–4 m, which is consistent with the emergence of surface breaks associated with this earthquake. Right-lateral strike-slip dominated on the fault segment along the Hinagu fault, but normal-slip components were estimated to make a significant contribution on the fault segment along the Futagawa fault. The large fault-parallel displacements recorded at two near-fault strong motion stations coincided with the spatiotemporal pattern of the fault slip history during the mainshock. The spatial relationship between the rupture areas of the foreshock and mainshock implies a complex fault structure in this region.

Tectonic inversion in the Caribbean-South American plate boundary: GPS geodesy, seismology, and tectonics of theMw6.7 22 April 1997 Tobago earthquake

On 22 April 1997 the largest earthquake recorded in the Trinidad-Tobago segment of the Caribbean-South American plate boundary zone (Mw 6.7) ruptured a shallow (~9 km), ENE striking (~250° azimuth), shallowly dipping (~28°) dextral-normal fault ~10 km south of Tobago. In this study, we describe this earthquake and related foreshock and aftershock seismicity, derive coseismic offsets using GPS data, and model the fault plane and magnitude of slip for this earthquake. Coseismic slip estimated at our episodic GPS sites indicates movement of Tobago 135 ± 6 to 68 ± 6 mm NNE and subsidence of 7 ± 9 to 0 mm. This earthquake was anomalous and is of interest because (1) its large component of normal slip and ENE strike are unexpected given the active E-W dextral shearing across the Caribbean-South American plate boundary zone, (2) it ruptured a normal fault plane with a low (~28°) dip angle, and (3) it reactivated and inverted the preexisting Tobago terrrane-South America ocean-continent (thrust) boundary that formed during early Tertiary oblique plate convergence.