skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


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 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.  more » « less
Award ID(s):
2026099
PAR ID:
10269145
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Bulletin of the Seismological Society of America
Volume:
111
Issue:
2
ISSN:
0037-1106
Page Range / eLocation ID:
1032 to 1042
Format(s):
Medium: X
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 of 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. 
    more » « less
  2. 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. 
    more » « less
  3. Geometric moment is accumulated and released throughout the earthquake cycle and can be imaged over both co- and inter-seismic time intervals with geodetic data. Recent dynamic earthquake cycle models have shown that substantial coseismic slip may occur in regions that exhibit low levels of moment accumulation during decadal-scale pre-earthquake epochs. We estimate that <1% of the region that slipped coseismically during the 2024 MW = 7.1 Hyuganada Sea (offshore southeastern Japan) earthquake was co-located with a region of pre-event moment accumulation. This stands in contrast with the behavior of the 2011 MW = 9.1 Tohoku-Oki(offshore northeastern Japan) earthquake, where ∼98% of the coseismic rupture area was co-located with pre-seismic moment accumulation. This co-location of coseismic slip and pre-seismic moment release proximal to the 2024 Hyuganada Sea earthquake serves as a real-world observation consistent with a novel theoretical prediction and is consistent with a perspective that decadal-scale pre-earthquake coupling estimates may provide snapshots of a highly time-variable fault system behavior, and may not necessarily serve as diagnostic identifiers of regions facing short-term earthquake hazard. 
    more » « less
  4. ABSTRACT Afterslip could help to reveal seismogenic fault structure. The 2020 Mw 6.3 Nima earthquake happened in a pull-apart basin within the Qiangtang block, central Tibetan plateau. Previous studies have explained the coseismic and early (<6 mo) postseismic deformation by rupture and afterslip on a normal fault bounding the western side of the basin. Here, we resolved the 19-month Interferometric Synthetic Aperture Radar-measured sequences of postseismic displacements that revealed a second postseismic displacement center ~12 km to the east of the main fault. Fitting the postseismic displacement requires afterslip on both the main fault and an antithetic fault that probably forms a y-shaped pair of conjugate faults in a negative flower structure. Stress-driven afterslip models suggest that the required afterslip on the antithetic fault could be triggered by coseismic rupture of the main fault or by a simultaneous rupture on the antithetic fault. The afterslip on both faults occurred mainly up-dip to the coseismic slip and has released moment ~15%–19% of that by the coseismic rupture. These results provide insights into active extension in the central Tibetan plateau and highlight the complex nature of fault rupture and afterslip. 
    more » « less
  5. On February 6, 2023, two large earthquakes occurred near the Turkish town of Kahramanmaraş. The moment magnitude (Mw) 7.8 mainshock ruptured a 310 km-long segment of the left-lateral East Anatolian Fault, propagating through multiple releasing step-overs. The Mw 7.6 aftershock involved nearby left-lateral strike-slip faults of the East Anatolian Fault Zone, causing a 150 km-long rupture. We use remote-sensing observations to constrain the spatial distribution of coseismic slip for these two events and the February 20 Mw 6.4 aftershock near Antakya. Pixel tracking of optical and synthetic aperture radar data of the Sentinel-2 and Sentinel-1 satellites, respectively, provide near-field surface displacements. High-rate Global Navigation Satellite System data constrain each event separately. Coseismic slip extends from the surface to about 15 km depth with a shallow slip deficit. Most aftershocks cluster at major fault bends, surround the regions of high coseismic slip, or extend outward of the ruptured faults. For the mainshock, rupture propagation stopped southward at the diffuse termination of the East Anatolian fault and tapered off northward into the Pütürge segment, some 20 km south of the 2020 Mw 6.8 Elaziğ earthquake, highlighting a potential seismic gap. These events underscore the high seismic potential of immature fault systems. 
    more » « less