Project description:We provide here a first-hand description of the coseismic surface effects caused by the Mw 6.4 Petrinja earthquake that hit central Croatia on 29 December 2020. This was one of the strongest seismic events that occurred in Croatia in the last two centuries. Field surveys in the epicentral area allowed us to observe and map primary coseismic effects, including geometry and kinematics of surface faulting, as well as secondary effects, such as liquefaction, sinkholes and landslides. The resulting dataset consists of homogeneous georeferenced records identifying 222 observation points, each of which contains a minimum of 5 to a maximum of 14 numeric and string fields of relevant information. The earthquake caused surface faulting defining a typical 'conjugate' fault pattern characterized by Y and X shears, tension cracks (T fractures), and compression structures (P shears) within a ca. 10 km wide (across strike), NW-SE striking right-lateral strike-slip shear zone (i.e., the Petrinja Fault Zone, PFZ). We believe that the results of the field survey provide fundamental information to improve the interpretation of seismological, GPS and InSAR data of this earthquake. Moreover, the data related to the surface faulting may impact future studies focused on earthquake processes in active strike-slip settings, integrating the estimates of slip amount and distribution in assessing the hazard associated with capable transcurrent faults.
Project description:We predict, with a model (earthquake stress model) that inverts the displacements documented at 163 GNSS onshore stations of the GEONET, the change of shear and normal stresses on the megathrust near the Japan Trench over the seven years before the 2011 Mw 9.0 Tohoku-Oki earthquake. We find three areas on the megathrust with greater accumulations of shear and normal stresses before the earthquake, which match the ruptured areas of the mainshock and two largest aftershocks (Mw 7.8 and 7.4) that occurred within half an hour after the mainshock. We also find that the change of normal stress on the fault before the earthquake is not uniform but increases in the up-dip portion (shallower depth) of the fault from the hypocenter and decreases in the down-dip portion. We infer that the occurrence of the giant earthquake at the shallow portion of the megathrust may be attributed to the increase of the normal stress there, which leads to an increase of fault shear strength and allows more elastic strain energy to accumulate to prepare for the next big earthquake. Based on these results we propose a new concept of the seismogenic asperity as the area of greater accumulations of shear and normal stresses. The method presented here may be useful for predicting the rupture zone of future large earthquakes.
Project description:Several recent studies have presented evidence that significant induced earthquakes occurred in a number of oil-producing regions during the early and mid-twentieth century related to either production or wastewater injection. We consider whether the 21 July 1952 Mw 7.5 Kern County earthquake might have been induced by production in the Wheeler Ridge oil field. The mainshock, which was not preceded by any significant foreshocks, occurred 98 days after the initial production of oil in Eocene strata at depths reaching 3 km, within ~1 km of the White Wolf fault (WWF). Based on this spatial and temporal proximity, we explore a potential causal relationship between the earthquake and oil production. While production would have normally be expected to have reduced pore pressure, inhibiting failure on the WWF, we present an analytical model based on industry stratigraphic data and best estimates of parameters whereby an impermeable splay fault adjacent to the main WWF could plausibly have blocked direct pore pressure effects, allowing the poroelastic stress change associated with production to destabilize the WWF, promoting initial failure. This proof-of-concept model can also account for the 98-day delay between the onset of production and the earthquake. While the earthquake clearly released stored tectonic stress, any initial perturbation on or near a major fault system can trigger a larger rupture. Our proposed mechanism provides an explanation for why significant earthquakes are not commonly induced by production in proximity to major faults.
Project description:On 17-18 September 2022, an earthquake sequence with a moment magnitude of 6.6 foreshock and a 7.0 mainshock occurred in southeast Taiwan along the Longitudinal Valley. Several surface breaks and collapsed buildings were observed after the event and one person died. The focal mechanisms of the foreshock and mainshock both had a west-dipping fault plane, which is different from the known active east-dipping boundary fault between the Eurasian Plate and the Philippine Sea Plate. Joint source inversions were performed to better understand the rupture mechanism of this earthquake sequence. The results show that the ruptures mainly occurred on a west-dipping fault. In the mainshock, the slip originated from the hypocenter and propagated toward the north with a rupture velocity of approximately 2.5 km/s. The east-dipping Longitudinal Valley Fault also ruptured, which could be passive and dynamically triggered by the significant rupture on the west-dipping fault. Most importantly, this source rupture model together with the occurrence of large local earthquakes over the past decade strongly supports the existence of the Central Range Fault, which is a west-dipping boundary fault that lies along the north to south ends of the Longitudinal Valley suture.
Project description:The Mexican subduction zone is an ideal location for studying subduction processes due to the short trench-to-coast distances that bring broad portions of the seismogenic and transition zones of the plate interface inland. Using a recently generated seismicity catalog from a local network in Oaxaca, we identified 20 swarms of earthquakes (M < 5) from 2006 to 2012. Swarms outline what appears to be a steeply dipping structure in the overriding plate, indicative of an origin other than the plate interface. This steeply dipping structure corresponds to the northern boundary of the Xolapa terrane. In addition, we observed an interesting characteristic of slow slip events (SSEs) where they showed a shift from trenchward motion toward an along-strike direction at coastal GPS sites. A majority of the swarms were found to correspond in time to the along-strike shift. We propose that swarms and SSEs are occurring on a sliver fault that allows the oblique convergence to be partitioned into trench-perpendicular motion on the subduction interface and trench-parallel motion on the sliver fault. The resistivity structure surrounding the sliver fault suggests that SSEs and swarms of earthquakes occur due to high fluid content in the fault zone. We propose that the sliver fault provides a natural pathway for buoyant fluids attempting to migrate upward after being released from the downgoing plate. Thus, sliver faults could be responsible for the downdip end of the seismogenic zone by creating drier conditions on the subduction interface trenchward of the sliver fault, promoting fast-slip seismogenic rupture behavior.
Project description:The 2010 MW 7.2 El Mayor-Cucapah, Mexico, earthquake ruptured multiple faults with different faulting mechanisms. Resolving the earthquake rupture process and its relation to the geometric fault complexities is critical to our understanding of the earthquake source physics, but doing so by conventional finite-fault inversion is challenging because modelling errors due to inappropriate assumptions about the fault geometry distort the solution and make robust interpretation difficult. Here, using a potency density tensor approach to finite-fault inversion, we inverted the observed teleseismic P waveforms of the 2010 El Mayor-Cucapah earthquake to simultaneously estimate the rupture process and the fault geometry. We found that the earthquake consisted of an initial normal faulting rupture, which was followed by a strike-slip bilateral rupture towards the southeast and northwest that originated on the northwest side of the epicentre. The southeastern rupture propagated back through the initial rupture area, but with strike-slip faulting. Although the northwestern rupture propagated across the left step in the Puerta fault-accommodation zone, the rupture was temporarily stalled by the associated change of the fault geometry. These results highlight the irregular rupture process, which involved a back-propagating rupture and fluctuating rupture propagation controlled the complexity of the fault system.
Project description:Seismic rupture in carbonate rocks influences fault friction behavior through thermal evolution and mineral reactions. Focusing on the 1959 Mw 7.2 Hebgen Lake event in western Yellowstone, Montana, the largest earthquake on a normal fault in the United States, we analyze fault rock microstructures and mineralogical changes to constrain frictional heating on the fault plane. We combine thermal maturity of organic matter, magnetic fabric, and thermomagnetic methods with scanning electron microscopy to unravel variations in peak frictional temperature along the fault slip surface. The mineral changes caused by coseismic heating (e.g., nanocalcite formation or goethite to hematite reaction) occur in patches along the fault mirror, hence reflecting considerable differences in frictional heat. While coseismic thermal heterogeneities have been reported in other rock types, this is the first time they are documented and quantified specifically in carbonates. Furthermore, these results provide new mineralogical criteria to quantify coseismic frictional heat in natural faults at temperatures lower than that of decarbonation and highlight the need to consider coseismic friction processes at a scale larger than most deformation experiments. For example, we document the critical role played by fault plane attitude (dip) at the scale of a few tens of centimeters in production of frictional heat. Our results emphasize that while coseismic decarbonation dynamically weakens carbonate-hosted faults, it may generally not occur along an entire fault plane.
Project description:Soil behavior is studied during the Tohoku earthquake, where abnormally high accelerations > 1 g were recorded. Based on vertical array records, models of soil behavior are constructed at 28 sites in northern Honshu (Tohoku region). They are compared with previously studied models of soil behavior in southern Tohoku and Kanto regions, where shock waves were identified as possible causes of the recorded high accelerations. Shear moduli did not reduce during strong motion at many sites, and the behavior of softer and denser soils was similar to a large extent. The nonlinearity of soil response during the Tohoku earthquake was weaker than that observed earlier during the 1995 Kobe and 2000 Tottori earthquakes (Mw ~6.7-6.8). Instead, a widespread soil hardening was found, most expressed at stations recorded the highest PGAs. To explain the observed features of soil behavior, two possible mechanisms are suggested, such as, 1) shock wave fronts generated by rupture propagation along the fault plane induce soil hardening and high PGAs; 2) soil compaction and hardening is a soil response to long-lasting dynamic loadings during the earthquake. Most likely we may expect similar effects of soil hardening and generation of high PGAs during other mega-thrust earthquakes in future.
Project description:The 28th December 1908 Messina earthquake (Mw 7.1), Italy, caused >80,000 deaths and transformed earthquake science by triggering the study of earthquake environmental effects worldwide, yet its source is still a matter of debate. To constrain the geometry and kinematics of the earthquake we use elastic half-space modelling on non-planar faults, constrained by the geology and geomorphology of the Messina Strait, to replicate levelling data from 1907-1909. The novelty of our approach is that we (a) recognise the similarity between the pattern of vertical motions and that of other normal faulting earthquakes, and (b) for the first time model the levelling data using the location and geometry of a well-known offshore capable fault. Our results indicate slip on the capable fault with a dip to the east of 70° and 5 m dip-slip at depth, with slip propagating to the surface on the sea bed. Our work emphasises that geological and geomorphological observations supporting maps of capable non-planar faults should not be ignored when attempting to identify the sources of major earthquakes.
Project description:Surface faulting earthquakes are known to cluster in time from historical and palaeoseismic studies, but the mechanism(s) responsible for clustering, such as fault interaction, strain-storage, and evolving dynamic topography, are poorly quantified, and hence not well understood. We present a quantified replication of observed earthquake clustering in central Italy. Six active normal faults are studied using 36Cl cosmogenic dating, revealing out-of-phase periods of high or low surface slip-rate on neighboring structures that we interpret as earthquake clusters and anticlusters. Our calculations link stress transfer caused by slip averaged over clusters and anti-clusters on coupled fault/shear-zone structures to viscous flow laws. We show that (1) differential stress fluctuates during fault/shear-zone interactions, and (2) these fluctuations are of sufficient magnitude to produce changes in strain-rate on viscous shear zones that explain slip-rate changes on their overlying brittle faults. These results suggest that fault/shear-zone interactions are a plausible explanation for clustering, opening the path towards process-led seismic hazard assessments.