Hector Mine Research Articles

Evaluation of the b Maps on the Faults of the Major (M > 7) South California Earthquakes

AbstractWe use the Godano et al. (2022, https://doi.org/10.1029/2021ea002205) method for evaluating the b maps of the faults associated with the largest earthquakes M ≥ 7.0 that occurred in California. The method allows an independent evaluation of the b parameter, avoiding the overlap of the cells and the omission of some earthquakes, while keeping all the available information in the catalog. We analyzed four large earthquakes: Landers, Hector Mine, Baja California, and Searles Valley. The maps obtained confirm that the b value can be considered as a strain meter and allow us to elucidate the presence of barriers, such as obstacles to the propagation of the fracture, on the fault of the analyzed earthquakes. A further estimated parameter is the time window during which aftershocks occur in the cell, Δt. This quantity is very useful for a better definition of the aftershock generation mechanism. It reveals where the stress is released in a short time interval and how the complexity of the faulting process controls the occurrence of aftershocks on the fault, and also the duration of the entire sequence.

Evidence for a prehistoric multifault rupture along the southern Calico fault system, Eastern California Shear Zone, USA

Abstract Geomorphic mapping and paleoseismologic data reveal evidence for a late Holocene multifault surface rupture along the Calico-Hidalgo fault system of the southern Eastern California Shear Zone (ECSZ). We have identified ~18 km of continuous surface rupture along the combined Calico and Hidalgo faults in the vicinity of Hidalgo Mountain in the southern Mojave Desert. Based on the freshness of geomorphic fault features and continuity of surface expression, we interpret this feature to reflect a simultaneous paleorupture of both faults. Displacement along the paleorupture is defined by 39 field measurements to be generally pure right-slip with a mean offset of 2.3 m. Scaling relationships for this offset amount imply that the original surface rupture length may have been ~82 km (corresponding to a M7.4 earthquake) and that much of the rupture trace was erased by subsequent erosion of sandy and unconsolidated valley alluvium. Eight luminescence ages from a paleoseismic trench across the paleorupture on the Hidalgo fault bracket the timing of the most recent rupture to 0.9–1.7 ka and a possible penultimate event at 5.5–6.6 ka. This timing is generally consistent with the known earthquake clusters in the southern ECSZ based on previous paleoseismic investigations. The ages of these earthquakes also overlap with the age brackets of the most recent events on the Calico fault 42 km to the north and the Mesquite Lake fault 40 km to the south from earlier work. Based on these age constraints and the expected surface rupture length, we propose that the Calico fault system experienced a major, multifault rupture that spanned the entire length of the fault system between the historical Landers and Hector Mine ruptures but preceded these events by ~1–2 k.y. Coulomb stress change modeling shows that the Calico paleorupture may have delayed the occurrence of the Landers-Hector Mine cluster by placing their respective faults in stress shadows and may have also prevented a triggered event from occurring on the Calico fault following the historic events. This work implies that closely spaced ruptures in complex shear zones may repel each other and thereby stretch out the duration of major earthquake clusters. These results also suggest that complex multifault ruptures in the ECSZ may not follow simple, repeatable patterns.

2019 M7.1 Ridgecrest earthquake slip distribution controlled by fault geometry inherited from Independence dike swarm

Faults often form through reactivation of pre-existing structures, developing geometries and mechanical properties specific to the system’s geologic inheritance. Competition between fault geometry and other factors (e.g., lithology) to control slip at Earth’s surface is an open question that is central to our knowledge of fault processes and seismic hazards. Here we use remote sensing data and field observations to investigate the origin of the 2019 M7.1 Ridgecrest, California, earthquake rupture geometry and test its impact on the slip distribution observed at Earth’s surface. Common geometries suggest the fault system evolved through reactivation of structures within the surrounding Independence dike swarm (IDS). Mechanical models testing a range of fault geometries and stress fields indicate that the inherited rupture geometry strongly controlled the M7.1 earthquake slip distribution. These results motivate revisiting the development of other large-magnitude earthquake ruptures (1992 M7.3 Landers, 1999 M7.1 Hector Mine) and tectonic provinces within the IDS.

Physical Properties of the Crust Influence Aftershock Locations

AbstractAftershocks do not uniformly surround a mainshock, and instead occur in spatial clusters. Spatially variable physical properties of the crust may influence the spatial distribution of aftershocks. I study four aftershock sequences in Southern California (1992 Landers, 1999 Hector Mine, 2010 El Mayor—Cucapah, and 2019 Ridgecrest) to investigate which physical properties are spatially correlated with aftershock occurrence. I find that aftershocks correlate with several properties, including measures of stress and stress change from the mainshock, fault structure, kinematics, seismic velocity, and heat flow. Aftershock spatial density exhibits an order of magnitude or more variation as a function of these properties. I determine simple empirical relations between each of the properties and the aftershock spatial density, and use these relations to construct new spatial models that describe aftershock locations. The new spatial models are a significant improvement over a simple base model, but do not fully capture the dense spatial clustering of aftershocks. Numerous spatially varying physical properties exhibit no (or poor) correlation with aftershock spatial density, including temperature, rock composition, and rheological properties that might be expected to control aftershock occurrence. These results suggest that while spatially variable physical properties appear to influence aftershock locations, more work is necessary in order to establish the connections between aftershock occurrence and the causative physical properties.

Effects of Parameters on Non-Linear Seismic Analysis of Fixed Offshore Platform

Fixed offshore structures in seismically active zones are likely to have deformations large enough to exceed the elastic limits of materials under strong ground motions. Therefore, investigation of nonlinear response of offshore structures under dynamic loadings is very interesting to structural engineers. Nonlinear behavior of structures is highly related with the soil properties and soil behavior. The main purpose of this study is to investigate and show the methodology for inelastic dynamic analyses of fixed-offshore structures under the 10,000-year earthquake loadings. The investigation in this study includes the effect of soil radiation damping, soil-pile gapping, and soil degradation on fixed offshore structures. Additionally, the effect of the direction of earthquake is included in this study. An offshore structure model ‘Platform A’ is modeled in CAPFOS. 10,000-year earthquake events are considered and for this purpose, two different earthquakes are used: Hector Mine Earthquake (1999) and Landers Earthquake (1992). These ground motion records are scaled to have a resulting earthquake that has a return period of 10,000-year. CAPFOS software is used to complete time-history analyses with these earthquakes. Results indicate that: (i) Considering the soil radiation damping negligible is not correct, because without the radiation damping, the total damping ratio is estimated as % 3, while including radiation damping raises the ratio to % 4.5, (ii) soil-pile gapping and soil radiation damping generate adverse effects on analyses outcome, and (iii) using high values of soil degradation factors is not always conservative.

Potential Fault Displacement Hazard Assessment Using Stochastic Source Models: A Retrospective Evaluation for the 1999 Hector Mine Earthquake

Surface fault displacement due to an earthquake affects buildings and infrastructure in the near-fault area significantly. Although approaches for probabilistic fault displacement hazard analysis have been developed and applied in practice, there are several limitations that prevent fault displacement hazard assessments for multiple locations simultaneously in a physically consistent manner. This study proposes an alternative approach that is based on stochastic source modelling and fault displacement analysis using Okada equations. The proposed method evaluates the fault displacement hazard potential due to a fault rupture. The developed method is applied to the 1999 Hector Mine earthquake from a retrospective perspective. The stochastic-source-based fault displacement hazard analysis method successfully identifies multiple source models that predict fault displacements in close agreement with observed GPS displacement vectors and displacement offsets along the fault trace. The case study for the 1999 Hector Mine earthquake demonstrates that the proposed stochastic-source-based method is a viable option in conducting probabilistic fault displacement hazard analysis.

From Hector Mine M7.1 to Ridgecrest M7.1 Earthquake. A Look from a 20-Year Perspective

The paper provides a comparative analysis of precursory phenomena in the ionosphere and atmosphere for two strong earthquakes of the same magnitude M7.1 that happened in the same region (North-East from Los Angeles) within a time span of 20 years, the Hector Mine and Ridgecrest earthquakes. Regardless of the similarity of their location (South-Eastern California, near 160 km one from another), there was one essential difference: the Hector Mine earthquake happened during geomagnetically disturbed conditions (essential in the sense of ionospheric precursors identification). In contrast, the quiet geomagnetic conditions characterized the period around the time of the Ridgecrest earthquake. The Hector mine earthquake happened in the middle of the rising phase of the 23-rd solar cycle characterized by high solar activity, while the Ridgecrest earthquake happened by the very end of the 24th cycle under very low solar activity conditions. We provide a comprehensive multi-factor analysis, determine the precursory period for both earthquakes and demonstrate the close similarity of ionospheric precursors. Unlike the majority of papers dealing with earthquake precursor identification based on the “abnormality” of observed time-series mainly determined by amplitude difference between “normal” (usually climatic) behavior and “abnormal” behavior with amplitudes exceeding some pre-established threshold, we used the technique of cognitive recognition of the precursors based on the physical mechanisms of their generation and the morphology of their behavior during the precursory period. These permits to uniquely identify precursors even in conditions of disturbed environment as it was around the time of the Hector Mine earthquake. We demonstrate the close similarity of precursors’ development for both events. The leading time of precursor appearance for the same region and similar magnitude was identical. For the Hector Mine it was 11 October 1999—5 days in advance—and for 2019 Ridgecrest it was 28 June—7 days before the mainshock and five days before the strongest foreshock.

Rupture Process of the Mainshock of the 2019 Ridgecrest Earthquake Sequence from Waveform Inversion with Empirical Green’s Functions

ABSTRACTThe 2019 Mw 7.1 mainshock of the Ridgecrest earthquake sequence, which was the first event exceeding Mw 7.0 in California since the 1999 Hector Mine earthquake, caused near-fault ground motions exceeding 0.5g and 70 cm/s. In this study, the rupture process and the generation mechanism of strong ground motions of the mainshock were investigated through waveform inversions of strong-motion data in the frequency range of 0.2–2.0 Hz using empirical Green’s functions (EGFs). The results suggest that the mainshock involved two large slip regions: the primary one with a maximum slip of approximately 4.4 m was centered ∼3 km northwest of the hypocenter, which was slightly shallower than the hypocenter, and the secondary one was centered ∼25 km southeast of the hypocenter. Outside these regions, the slip was rather small and restricted to deeper parts of the fault. A relatively small rupture velocity of 2.1 km/s was identified. The robustness of the slip model was examined by conducting additional inversion analyses with different combinations of EGF events and near-fault stations. In addition, using the preferred slip model, we synthesized strong motions at stations that were not used in the inversion analyses. The synthetic waveforms captured the timing of the main phases of observed waveforms, indicating the validity of the major spatiotemporal characteristics of the slip model. Our large slip regions are also generally visible in the models proposed by other researchers based on different datasets and focusing on lower frequency ranges (generally lower than 0.5 Hz). In particular, two large slip regions in our model are very consistent with two of the four subevents identified by Ross et al. (2019), which may indicate that part of the large slip regions that generated low-frequency ground motions also generated high-frequency ground motions up to 2.0 Hz during the Ridgecrest mainshock.

The Lavic Lake Fault: A Long-Term Cumulative Slip Analysis via Combined Field Work and Thermal Infrared Hyperspectral Airborne Remote Sensing

The 1999 Hector Mine earthquake ruptured to the surface in eastern California, with >5 m peak right-lateral slip on the Lavic Lake fault. The cumulative offset and geologic slip rate of this fault are not well defined, which inhibits tectonic reconstructions and risk assessment of the Eastern California Shear Zone (ECSZ). With thermal infrared hyperspectral airborne imagery, field data, and auxiliary information from legacy geologic maps, we created lithologic maps of the area using supervised and unsupervised classifications of the remote sensing imagery. We optimized a data processing sequence for supervised classifications, resulting in lithologic maps over a test area with an overall accuracy of 71 ± 1% with respect to ground-truth geologic mapping. Using all of the data and maps, we identified offset bedrock features that yield piercing points along the main Lavic Lake fault and indicate a 1036 +27/−26 m net slip, with 1008 +14/−17 m horizontal and 241 +51/−47 m vertical components. For the contribution from distributed shear, modern off-fault deformation values from another study imply a larger horizontal slip component of 1276 +18/−22 m. Within the constraints, we estimate a geologic slip rate of <4 mm/yr, which does not increase the sum geologic Mojave ECSZ rate to current geodetic values. Our result supports previous suggestions that transient tectonic activity in this area may be responsible for the discrepancy between long-term geologic and present-day geodetic rates.

Thin crème brûlée rheological structure for the Eastern California Shear Zone

AbstractSince the occurrence of the 1992 CE Mw 7.3 Landers and 1999 Mw 7.1 Hector Mine earthquakes in the Mojave Desert (California, USA), postseismic deformation following both earthquakes has been intensively studied, and models with a strong crust overlying a low-viscosity mantle asthenosphere have been favored. However, we recently found that the near-field postseismic transients after the two earthquakes have lasted longer than previously thought, which requires a revision of the postseismic modeling. Our new modeling results based on the revised postseismic transients show that: (1) the effective viscosity of the lower crust beneath the Mojave region at the decadal time scale is ∼2 × 1020 Pa·s (transient viscosity ∼2 × 1019 Pa·s), i.e., only ∼5 times that of the underlying mantle asthenosphere, and (2) the transient viscosity of the upper mantle exhibits a time-dependent increase, providing fresh geodetic evidence for frequency-dependent rheology (e.g., Andrade or extended Burgers rheology). The inferred transient rheology for the first year agrees well with that obtained for the July 2019 Mw 6.4 and Mw 7.1 Ridgecrest earthquakes ∼180 km north of the two Mojave events. Our modeling results support a thin crème brûlée model for the Eastern California Shear Zone (part of the Pacific-North America plate boundary) in which both the lower crust and the upper mantle exhibit ductility at decadal time scales.

Simulation of Broadband Ground Motion by Superposing High-Frequency Empirical Green’s Function Synthetics on Low-Frequency Spectral-Element Synthetics

Seismic wave-propagation simulations are limited in their frequency content by two main factors: (1) the resolution of the seismic wave-speed structure of the region in which the seismic waves are propagated through; and (2) the extent of our understanding of the rupture process, mainly on the short length scales. For this reason, high-frequency content in the ground motion must be simulated through other means. Toward this end, we adopt a variant of the classical empirical Green’s function (EGF) approach of summing, with suitable time shift, recorded seismograms from small earthquakes in the past to generate high-frequency seismograms (0.5–5.0 Hz) for engineering applications. We superimpose these seismograms on low-frequency seismograms, computed from kinematic source models using the spectral element method, to produce broadband seismograms. The non-uniform time- shift scheme used in this work alleviates the over-estimation of high-frequency content of the ground motions observed. We validate the methodology by simulating broadband motions from the 1999 Hector Mine and the 2006 Parkfield earthquakes and comparing them against recorded seismograms.

Introduction to the Special Section on the 2019 Ridgecrest, California, Earthquake Sequence

Following the 16 October 1999 M_w 7.1 Hector Mine earthquake, no earthquake larger than magnitude 6.0 occurred in southern California for nearly 20 yr. That lull came to an end with the occurrence of the 4 July 2019 M_w 6.4 earthquake at 22:40 UTC (10:33 a.m. local time). Located in a relatively sparsely populated part of the state, with most of its rupture within the China Lake Naval Air Weapons Station (NAWSCL) adjacent to the civilian town of Ridgecrest (population: 26,000), the earthquake was felt widely across most of southern California. Immediate response efforts were mounted to investigate the rupture and assess damage (Stewart et al., 2019) and to calculate aftershock probabilities (see Hardebeck, 2020; also see Data and Resources). The earthquake occurred in a region historically characterized by higher than average aftershock productivity (Hardebeck et al., 2019), leading to a higher‐than‐average probabilities for both large aftershocks and subsequent larger events. The latter odds, although still small, were borne out with the occurrence of 6 July M_w 7.1 earthquake at 03:19:53 UTC (8:19:53 p.m. local time on 5 July) ∼34 hr after the initial M_w 6.4. Arguably, the M_w 7.1 could be considered an aftershock, but, by convention, this earthquake is regarded as the mainshock, whereas the initial M_w 6.4 is considered a foreshock.

Coulomb stress models for the 2019 Ridgecrest, California earthquake sequence

We present a Coulomb stress analysis for the July 4th Mw 6.4 and July 5th Mw 7.1 earthquakes of the 2019 Ridgecrest, California earthquake sequence. We find that the Coulomb stress on the Mw 6.4 rupture plane was increased by the 1872 Owens Valley earthquake and slightly increased by the 1952 Kern County earthquake. Like- wise, the rupture plane of the July 5th Mw 7.1 earthquake was promoted by both the Owens Valley and Kern County earthquakes, as well as being promoted slightly by the 1992 Landers earthquake and the 1999 Hector Mine earthquake. To investigate the regional Coulomb stress change from the Ridgecrest earthquake sequence, we calculate Coulomb stress change on nearby faults from both the Mw 6.4 and Mw 7.1 earthquakes. We find that the Mw 6.4 earthquake promoted failure at the epicen- ter of the Mw 7.1 earthquake, which struck 34 h later. Furthermore, the Mw 6.4 and Mw 7.1 Ridgecrest earthquakes increased Coulomb stress on the Garlock fault to the west of the rupture and decreased Coulomb stress to the east. The Panamint Valley fault was relaxed by the Ridgecrest earthquakes. The most significant stress change is an increase of Coulomb stress to the south of the Ridgecrest Sequence on the nearby Blackwater fault. The Blackwater fault, situated between the Lan- ders and Ridgecrest earthquakes, was loaded significantly by the Mw 7.1 Ridgecrest earthquake (up to 0.12 MPa) and by lesser amounts by the Landers and Hector Mine earthquakes. Of the known active faults in the region, including the Garlock fault, the Blackwater fault appears to have been loaded the most by the Ridgecrest earth- quakes.

Rupture Process of the 2019 Ridgecrest, California Mw 6.4 Foreshock and Mw 7.1 Earthquake Constrained by Seismic and Geodetic Data

ABSTRACTThe 2019 Ridgecrest earthquake sequence culminated in the largest seismic event in California since the 1999 Mw 7.1 Hector Mine earthquake. Here, we combine geodetic and seismic data to study the rupture process of both the 4 July Mw 6.4 foreshock and the 6 July Mw 7.1 mainshock. The results show that the Mw 6.4 foreshock rupture started on a northwest-striking right-lateral fault, and then continued on a southwest-striking fault with mainly left-lateral slip. Although most moment release during the Mw 6.4 foreshock was along the southwest-striking fault, slip on the northwest-striking fault seems to have played a more important role in triggering the Mw 7.1 mainshock that happened ∼34 hr later. Rupture of the Mw 7.1 mainshock was characterized by dominantly right-lateral slip on a series of overall northwest-striking fault strands, including the one that had already been activated during the nucleation of the Mw 6.4 foreshock. The maximum slip of the 2019 Ridgecrest earthquake was ∼5 m, located at a depth range of 3–8 km near the Mw 7.1 epicenter, corresponding to a shallow slip deficit of ∼20%–30%. Both the foreshock and mainshock had a relatively low-rupture velocity of ∼2 km/s, which is possibly related to the geometric complexity and immaturity of the eastern California shear zone faults. The 2019 Ridgecrest earthquake produced significant stress perturbations on nearby fault networks, especially along the Garlock fault segment immediately southwest of the 2019 Ridgecrest rupture, in which the coulomb stress increase was up to ∼0.5 MPa. Despite the good coverage of both geodetic and seismic observations, published coseismic slip models of the 2019 Ridgecrest earthquake sequence show large variations, which highlight the uncertainty of routinely performed earthquake rupture inversions and their interpretation for underlying rupture processes.

Localization and coalescence of seismicity before large earthquakes

SUMMARY We examine localization processes of low magnitude seismicity in relation to the occurrence of large earthquakes using three complementary analyses: (i) estimated production of rock damage by background events, (ii) evolving occupied fractional area of background seismicity and (iii) progressive coalescence of individual earthquakes into clusters. The different techniques provide information on different time scales and on the spatial extent of weakened damaged regions. Techniques (i) and (ii) use declustered catalogues to avoid the occasional strong fluctuations associated with aftershock sequences, while technique (iii) examines developing clusters in entire catalogue data. We analyse primarily earthquakes around large faults that are locked in the interseismic periods, and examine also as a contrasting example seismicity from the creeping Parkfield section of the San Andreas fault. Results of analysis (i) show that the M > 7 Landers 1992, Hector Mine 1999, El Mayor-Cucapah 2010 and Ridgecrest 2019 main shocks in Southern and Baja California were preceded in the previous decades by generation of rock damage around the eventual rupture zones. Analysis (ii) reveals localization (reduced fractional area) 2–3 yr before these main shocks and before the M > 7 Düzce 1999 earthquake in Turkey. Results with technique (iii) indicate that individual events tend to coalesce rapidly to clusters in the final 1–2 yr before the main shocks. Corresponding analyses of data from the Parkfield region show opposite delocalization patterns and decreasing clustering before the 2004 M6 earthquake. Continuing studies with these techniques, combined with analysis of geodetic data and insights from laboratory experiments and model simulations, might improve the ability to track preparation processes leading to large earthquakes.

Surface Displacement Distributions for the July 2019 Ridgecrest, California, Earthquake Ruptures

ABSTRACTSurface rupture in the 2019 Ridgecrest, California, earthquake sequence occurred along two orthogonal cross faults and includes dominantly left-lateral and northeast-striking rupture in the Mw 6.4 foreshock and dominantly right-lateral and northwest-striking rupture in the Mw 7.1 mainshock. We present >650 field-based, surface-displacement observations for these ruptures and synthesize our results into cumulative along-strike displacement distributions. Using these data, we calculate displacement gradients and compare our results with historical strike-slip ruptures in the eastern California shear zone. For the Mw 6.4 rupture, we report 96 displacements measured along 18 km of northeast-striking rupture. Cumulative displacement curves for the rupture yield a mean left-lateral displacement of 0.3–0.5 m and maximum of 0.7–1.6 m. Net mean vertical displacement based on the difference of down-to-the-west (DTW) and down-to-the-east (DTE) displacement curves is close to zero (0.02 m DTW). The Mw 6.4 displacement distribution shows that the majority of displacement occurred southwest of the intersection with the Mw 7.1 rupture. The Mw 7.1 rupture is northwest-striking and 50 km long based on 576 field measurements. Displacement curves indicate a mean right-lateral displacement of 1.2–1.7 m and a maximum of 4.3–7.0 m. Net vertical displacement in the rupture averages 0.3 m DTW. The Mw 7.1 displacement distributions demonstrate that maximum displacement occurred along a 12-km-long portion of the fault near the Mw 7.1 epicenter, releasing 66% of the geologically based seismic moment along 24% of the total rupture length. Using our displacement distributions, we calculate kilometer-scale displacement gradients for the Mw 7.1 rupture. The steepest gradients (∼1–3 m/km) flank the 12-km-long region of maximum displacement. In contrast, gradients for the 1992 Mw 7.3 Landers and 1999 Mw 7.1 Hector Mine earthquakes are <0.6 m/km. Our displacement distributions are important for understanding the influence of cross-fault rupture on Mw 6.4 and 7.1 rupture length and displacement and will facilitate comparisons with distributions generated remotely and at broader scales.

What Can Surface-Slip Distributions Tell Us about Fault Connectivity at Depth?

ABSTRACT Fault systems with stepovers and gaps along strike are ubiquitous in nature, and many modern earthquakes (e.g., 1992 Landers, 1999 Hector Mine, 2016 Kaikōura, and 2019 Ridgecrest) have shown that ruptures can readily propagate across some disconnections, while being halted by others. It is quite possible, however, that many faults that appear discontinuous at the surface are in fact connected at depth, facilitating throughgoing rupture, and potentially increasing earthquake size. The present work explores whether the mapped surface slip in an earthquake is indicative of the connectivity of the fault system at depth. If there is a signal of subsurface connectivity in the surface-slip pattern, then the connectivity of the system could potentially be inferred. Through 3D dynamic rupture modeling of faults with along-strike gaps of various depths, I explore whether the amplitude or the spatial distribution of slip after an earthquake could be used to diagnose fault connectivity at depth. I find that, in general, fault segments that are connected up to shallow depths of 1–2 km and are relatively long along strike compared to the seismogenic depth tend to have higher slip gradients at their edges than faults that are connected at greater depth, or that are disconnected to the bottom of the seismogenic zone. Systematic slip gradient differences at fault segment edges have been recorded in past earthquakes, giving hope that the modeled effect can be detected in many cases, even though mapped surface slip is affected by a number of different sources of heterogeneity. The results provide an alternative explanation for observations that stepovers that allow throughgoing rupture tend to have higher slip gradients than those at which rupture terminates: perhaps many such stepovers are connected at depth, which could persistently favor throughgoing rupture. There may be implications for interpretation of apparent fault discontinuities worldwide.

Coseismic Displacements and Surface Fractures from Sentinel-1 InSAR: 2019 Ridgecrest Earthquakes

Abstract Interferometric Synthetic Aperture Radar is an important tool for imaging surface deformation from large continental earthquakes. Here, we present maps of coseismic displacement and strain from the 2019 Ridgecrest earthquakes using multiple Sentinel-1 images. We provide three types of interferometric products. (1) Standard interferograms from two look directions provide an overview of the deformation and can be used for modeling coseismic slip. (2) Phase gradient maps from stacks of coseismic interferograms provide high-resolution (∼30 m) images of strain concentration and surface fracturing that can be used to guide field surveys. (3) High-pass filtered, stacked, unwrapped phase is decomposed into east–west and up–down, south–north components and is used to determine the sense of fault slip. The resulting phase gradient maps reveal over 300 surface fractures, including triggered slip on the Garlock fault. The east–west component of high-pass filtered phase reveals the polarity of the strike-slip offset (right lateral or left lateral) for many of the fractures. We find a small number of fractures that have slip polarity that is retrograde to the background tectonic stress. This is similar to observations of retrograde slip observed near the 1999 Mw 7.1 Hector Mine rupture, but the Ridgecrest observations are more completely imaged by the frequent and high-quality acquisitions from the twin Sentinel-1 spacecrafts. Determining whether the retrograde features are triggered slip on existing faults, or compliant fault deformation in response to stress changes from the Ridgecrest earthquakes, or new Coulomb-style failures, will require additional field work, modeling, and analysis.

Transient Deformation in California From Two Decades of GPS Displacements: Implications for a Three-Dimensional Kinematic Reference Frame.

Our understanding of plate boundary deformation has been enhanced by transient signals observed against the backdrop of time‐independent secular motions. We make use of a new analysis of displacement time series from about 1,000 continuous Global Positioning System (GPS) stations in California from 1999 to 2018 to distinguish tectonic and nontectonic transients from secular motion. A primary objective is to define a high‐resolution three‐dimensional reference frame (datum) for California that can be rapidly maintained with geodetic data to accommodate both secular and time‐dependent motions. To this end, we compare the displacements to those predicted by a horizontal secular fault slip model for the region and construct displacement and strain rate fields. Over the past 19 years, California has experienced 19 geodetically detectable earthquakes and widespread postseismic deformation. We observe postseismic strain rate variations as large as 1,000 nstrain/year with moment releases equivalent up to an Mw6.8 earthquake. We find significant secular differences up to 10 mm/year with the fault slip model, from the Mendocino Triple Junction to the southern Cascadia subduction zone, the northern Basin and Range, and the Santa Barbara channel. Secular vertical uplift is observed across the Transverse Ranges, Coastal Ranges, Sierra Nevada, as well as large‐scale postseismic uplift after the 1999 Mw7.1 Hector Mine and 2010 Mw7.2 El Mayor‐Cucapah earthquakes. We also identify areas of vertical land motions due to anthropogenic, natural, and magmatic processes. Finally, we demonstrate the utility of the kinematic datum by improving the accuracy of high‐spatial‐resolution 12‐day repeat‐cycle Sentinel‐1 Interferometric Synthetic Aperture Radar displacement and velocity maps.

Fractal behaviour of distances between consecutive aftershocks: the examples of Landers, Northridge and Hector Mine Southern California mainshocks

Distances between consecutive aftershocks are analysed by means of mono- and multifractal theory with the aim of quantifying the complexity of the physical mechanism governing them, as well as their predictability and predictive instability. Hausdorff, Ha, and Hurst, H, exponents are determined by semivariograms and rescaled analysis, respectively. The exponent β of the power law describing power spectral contents is also quantified. These three parameters permit a generation of fractional Gaussian noise, fGn, simulating distances. The complexity and predictive instability of physical mechanism generating the series of distances is quantified by means of the correlation dimension, μ*, the Kolmogorov entropy, κ, and the Lyapunov exponents, λi, which are based on the reconstruction theorem formulation. Additionally, the multifractal detrended fluctuation analysis, MF-DFA, contributes with a different point of view to quantify the complexity of the series, in terms of fractal spectral width, W, spectral asymmetry, B, and the critical Holder exponent, a0. By one hand, the MF-DFA is applied to the complete set of distances characterising the whole aftershock process. By the other hand, the MF-DFA is applied to segments of the series of distances with the aim of determining the evolution of the complexity since the mainshock up to the end of the stress relaxation process. Finally, an ARIMA multilinear regression process is applied to obtain some improvements, in comparison with fGn simulations, on the prediction of distances. The database for this analysis is obtained from the Southern California Seismic Network (SCSN) catalogue. Three series of aftershocks equalling to or exceeding magnitudes of 2.0, assuring seismic catalogue completeness, and associated with Landers (06/28/1992), Northridge (01/17/1994) and Hector Mine (10/16/1999) mainshocks are obtained. It is worth mentioning that common mono-multifractal behaviour for the three aftershocks series is not detected, whatever aftershock periods or segments of them are considered.