skip to main content

Title: Refining the 2018 Mw 7.5 Papua New Guinea Earthquake Fault-Slip Model Using Subpixel Offset
ABSTRACT The Mw 7.5 earthquake that struck central Papua New Guinea in 2018 was the largest event ever recorded in the region with modern seismic instruments. The ground motions associated with this event also triggered widespread landslides and affected more than 500,000 people. However, due to the absence of a local seismic and Global Positioning System network in the vicinity, the fault location, system, and slip distribution of this earthquake are not well documented. In this study, we use the subpixel offset method on the Copernicus Sentinel-1 Synthetic Aperture Radar (SAR) images to calculate the 3D coseismic displacement of the 2018 Papua New Guinea earthquake. The results show clear fault traces that suggest coseismic slip on the Mubi fault and the Mananda fault that triggered landslides that spread out in a more than 260  km2 region. Finite-source inversions based on the subpixel offset measurements show up to 4.1 and 6.5 m coseismic slip on the Mubi and Mananda faults, respectively. Despite higher data uncertainty (∼0.4–0.8  m) of the subpixel offset data, synthetic resolution tests show resolvable slip above 8 km in depth. The lack of shallower slip on the west side of the Mananda fault could be due to an inflated geothermal gradient near the more » dormant volcano, Mount Sisa, as a slip barrier. The result of the coulomb stress change suggests possible southeastward slip propagation from the Mananda fault to the Mubi fault. Our work successfully resolves 3D coseismic displacement in highly vegetated terrains and demonstrates the feasibility of using the subpixel offset on SAR images to help our understanding of regional active tectonic systems. « less
Award ID(s):
Publication Date:
Journal Name:
Bulletin of the Seismological Society of America
Page Range or eLocation-ID:
1032 to 1042
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The 2018 Palu tsunami contributed significantly to the devastation caused by the associated7.5 earthquake. This began a debate about how the moderate size earthquake triggered such a large tsunami within Palu Bay, with runups of more than 10 m. The possibility of a large component of vertical coseismic deformation and submarine landslides have been considered as potential explanations. However, scarce instrumental data have made it difficult to resolve the potential contributions from either type of source. We use tsunami waveforms derived from social media videos in Palu Bay to model the possible sources of the tsunami. We invert InSAR data with different fault geometries and use the resulting seafloor displacements to simulate tsunamis. The coseismic sources alone cannot match both the video‐derived time histories and surveyed runups. Then we conduct a tsunami source inversion using the video‐derived time histories and a tide gauge record as inputs. We specify hypothetical landslide locations and solve for initial tsunami elevation. Our results, validated with surveyed runups, show that a limited number of landslides in southern Palu Bay are sufficient to explain the tsunami data. The Palu tsunami highlights the difficulty in accurately capturing with tide gauges the amplitude and timing ofmore »short period waves that can have large impacts at the coast. The proximity of landslides to locations of high fault slip also suggests that tsunami hazard assessment in strike‐slip environments should include triggered landslides, especially for locations where the coastline morphology is strongly linked to fault geometry.

    « less
  2. Abstract

    Large earthquakes can construct mountainous topography by inducing rock uplift but also erode mountains by causing landslides. Observations following the 2008 Wenchuan earthquake show that landslide volumes in some cases match seismically induced uplift, raising questions about how the actions of individual earthquakes accumulate to build topography. Here we model the two‐dimensional surface displacement field generated over a full earthquake cycle accounting for coseismic deformation, postseismic relaxation, landslide erosion, and erosion‐induced isostatic compensation. We explore the related volume balance across different seismotectonic and topographic conditions and revisit the Wenchuan case in this context. The ratio (Ω) between landslide erosion and uplift is most sensitive to parameters determining landslide volumes (particularly earthquake magnitudeMw, seismic energy source depth, and failure susceptibility, as well as the seismological factor responsible for triggering landslides), and is moderately sensitive to the effective elastic thickness of lithosphere,Te. For a specified magnitude, more erosive events (higher Ω) tend to occur at shallower depth, in thicker‐Telithosphere, and in steeper, more landslide‐prone landscapes. For given landscape and seismotectonic conditions, the volumes of both landslides and uplift to first order positively scale withMwand seismic momentMo. However, higherMwearthquakes generate lower landslide and uplift volumes per unitMo, suggesting lower efficiency in themore »use of seismic energy to drive topographic change. With our model, we calculate the long‐term average seismic volume balance for the eastern Tibetan region and find that the net topographic effect of earthquakes in this region tends to be constructive rather than erosive. Overall, destructive events are rare when considering processes over the full earthquake cycle, although they are more likely if only considering the coseismic volume budget (as was the case for the 2008 Wenchuan earthquake where landsliding substantially offset coseismic uplift). Irrespective of the net budget, our results suggest that the erosive power of earthquakes plays an important role in mountain belt evolution, including by influencing structures and spatial patterns of deformation, for example affecting the wavelength of topography.

    « less
  3. 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-qualitymore »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.« less
  4. Abstract Despite a lack of modern large earthquakes on shallowly dipping normal faults, Holocene M w  > 7 low-angle normal fault (LANF; dip<30°) ruptures are preserved paleoseismically and inferred from historical earthquake and tsunami accounts. Even in well-recorded megathrust earthquakes, the effects of non-linear off-fault plasticity and dynamically reactivated splay faults on shallow deformation and surface displacements, and thus hazard, remain elusive. We develop data-constrained 3D dynamic rupture models of the active Mai’iu LANF that highlight how multiple dynamic shallow deformation mechanisms compete during large LANF earthquakes. We show that shallowly-dipping synthetic splays host more coseismic slip and limit shallow LANF rupture more than steeper antithetic splays. Inelastic hanging-wall yielding localizes into subplanar shear bands indicative of newly initiated splay faults, most prominently above LANFs with thick sedimentary basins. Dynamic splay faulting and sediment failure limit shallow LANF rupture, modulating coseismic subsidence patterns, near-shore slip velocities, and the seismic and tsunami hazards posed by LANF earthquakes.
  5. Abstract

    We develop finite element models of the coseismic displacement field accounting for the 3D elastic structures surrounding the epicentral area of the 2019 Ridgecrest earthquake sequence containing two major events of Mw7.1 and Mw6.4. The coseismic slip distribution is inferred from the surface displacement field recorded by interferometric synthetic aperture radar. The rupture dip geometry is further optimized using a novel nonlinear‐crossover‐linear inversion approach. It is found that accounting for elastic heterogeneity and fault along‐strike curvilinearity improves the fit to the observed displacement field and yields a more accurate estimate of geodetic moment and Coulomb stress changes. We observe spatial correlations among the locations of aftershocks and patches of high slip, and rock anomalous elastic properties, suggesting that the shallow crust's elastic structures possibly controlled the Ridgecrest earthquake sequence. Most of the coseismic slip with a peak slip of 7.4 m at 3.6 km depth occurred above a zone of reducedS‐wave velocity and significant post‐Mw7.1 afterslip. This implies that viscous materials or fluid presence might have contributed to the low rupture velocity of the mainshock. Moreover, the zone of high slip on the northwest‐trending fault segment is laterally bounded by two aftershock clusters, whose location is characterized by intermediate rock rigidity.more »Notably, some minor orthogonal faults consistently end above a subsurface rigid body. Overall, these observations of structural controls improve our understandings of the seismogenesis within incipient fault systems.

    « less