Geology Reports⌕ Search

Geology topics

William D. Barnhart

Publications and source records attributed to William D. Barnhart.

At least 19 recordsLinked to original sources
USGS2026

Surface rupture and slip distribution of the 2025 Mw7.7 Mandalay earthquake and updated length scaling of supershear earthquakes

The 2025 M w 7.7 Mandalay, Burma (Myanmar), earthquake ruptured 475 km of the central Sagaing fault and is the longest continental strike-slip rupture on record. The observed rupture length is 1.6–4.7 times the value expected (100–300 km) from existing length-magnitude scaling relations for strike-slip earthquakes. The earthquake resulted from shallow dextral faulting and ruptured bilaterally with supershear speeds south of the epicenter, rupturing close to three major cities in Myanmar and exposing over six million people to violent or extreme shaking. We report on the surface rupture character, length, and slip distribution based on sub-pixel correlation of Sentinel-2 (10 m) and Planet Dove (3 m) optical images and visual analysis of SkySat and WorldView (0.3–0.5 m) optical images. The earthquake had moderate surface slip (average = 3.3 m, maximum = 5.6 m, 25–75% range = 3.0–4.0 m), narrow deformation zone width (1–10 pixels in sub-pixel correlation and up to 190 meters for the detailed surface rupture mapping), and simple fault geometry (no stepovers or large changes in strike, 87% of the rupture that was mapped in detail is single-stranded). We attribute the extreme length of the Mandalay earthquake to supershear rupture speed, simple fault geometry, narrow down-dip width, and moderate surface slip. Based on a compilation of 25 supershear strike-slip earthquakes ( M w 6.5–8.6; 1979–2025), we find that the rupture length of supershear earthquakes does not fit empirical scaling relationships for strike-slip earthquakes that predict length from magnitude. A length-magnitude scaling relationship based on supershear earthquakes has a best fit of log 10 (surface rupture length) = 0.89 M w – 4.44, indicating that supershear earthquakes tend to be longer than their subshear counterparts for any given magnitude and thus may expose a greater population to shaking.

Mandalay↗
USGS2025

Surface rupture from an aftershock: Remote observations from the January 2024 Wushi (Aykol), China, earthquakes

On 22 January 2024, a moment magnitude ( M w ⁠⁠ ) 7.0 earthquake occurred in the southern Tian Shan region of western China, followed by an M w 5.7 aftershock on 29 January. The M w 7.0 mainshock occurred on a northwest‐dipping fault plane, at a hypocentral depth of 13 km. Both Interferometric Synthetic Aperture Radar (InSAR) and optical satellite imagery observations indicate a lack of surface rupture for the M w 7.0 earthquake. However, mapping and subpixel correlation of high‐resolution optical satellite imagery and InSAR observations spanning the 29 January M w 5.7 thrust‐faulting aftershock support the interpretation of fault surface rupture. We observed ∼4.8 km of discontinuous surface rupture associated with the M w 5.7 aftershock. The surface rupture trace trends subparallel to the northeast–southwest bedding orientation within the limb of a syncline. Subpixel image correlation and digital elevation model (DEM) differencing indicate mean on‐fault offsets of 1 m, with a peak of ∼1.4 m, along the fault. Given these observations, we interpret the M w 5.7 aftershock produced surface rupture on an unmapped flexural‐slip fault within a syncline in the hanging wall of the Maidan fault or on a back thrust associated with the Maidan fault. Pre‐existing fault features along the surface rupture trace in the landscape from pre‐earthquake DEMs indicate that the fault that ruptured in the M w 5.7 aftershock may have ruptured in the past. Surface rupture only from an aftershock in a seismic sequence is a critical observation that highlights the need to understand secondary or distributed faulting in larger earthquakes for seismic and fault displacement hazard analysis.

Tian Shan region↗
USGS2025

Estimating earthquake source depth using teleseismic broadband waveform modeling at the USGS National Earthquake Information Center

The U.S. Geologic Survey National Earthquake Information Center (NEIC) monitors global seismicity, producing a catalog of earthquake source parameters in near-real-time to provide information that can help mitigate the societal impact of earthquakes. The NEIC commonly relies on teleseismic observations to constrain earthquake source parameters (e.g., location, depth, magnitude, and mechanism) due to a lack of local and regional observations. For these ‘teleseismic’ events, depth phase (i.e., pP , sP ) arrival time observations provide the best estimate on source depth. However, depth phases are often difficult to accurately identify and/or pick. Therefore, NEIC relies on waveform modeling, such as those determined from W-phase ( M ww ), body wave ( M wb ), and regional ( M wr ) moment tensor estimations, to provide constraints on source depth. While depth estimates from these approaches are informative, higher frequency observations provide more precise estimates because depth phases are more prominently observed at higher frequencies. Here, we present NEIC’s relatively high-frequency (~0.04 to 1 Hz) teleseismic waveform modeling approach, termed Synthetic Depth Phase Modeling (SynDepth), for determining source depth. SynDepth was developed to provide NEIC with a tool that enables rapid, accurate, and quantifiable estimates of earthquake source depth in cases where locator depths are not reliable. This relatively simple and fast procedure searches over 1 km-incremented source depths and an expanding triangular source-time function to find the best-fitting solution. We compare automatic SynDepth solutions for a dataset of 1,216 earthquakes (M5.5-M7.6) between 2017 and 2021 to NEIC-derived depth estimates from other methods. Our approach provides a robust depth estimate for earthquakes lacking local arrival time data, and it minimizes the need for analyst review of depth-phase picks ( pP , sP ) or using predefined ‘fixed’ depths.

Seismological Research Letters↗
USGS2025

Ultralong, supershear rupture of the 2025 Mw 7.7 Mandalay earthquake reveals unaccounted risk

The 28 March 2025 moment magnitude ( M w ) 7.7 earthquake in Mandalay, Burma (Myanmar), ruptured 475 kilometers of the Sagaing Fault, which was more than twice the length predicted by magnitude scaling relationships. Kinematic slip models and observation of a Rayleigh Mach wave that passed through parts of Thailand confirmed that rupture occurred at supershear velocities of greater than 5 kilometers per second. The anomalous length exposed a vast population to violent near-fault shaking. The Mandalay earthquake is a modern analog for the M w 7.9 1906 San Francisco earthquake, another atypically long and fast rupture. Probabilistic seismic hazard analyses use scaling relations that do not account for such long ruptures at moderate magnitudes. This limitation, in conjunction with a likely increased population and infrastructure exposure for atypically long ruptures, contributes to a potential mischaracterization of seismic risk.

Mandalay↗
USGS2024

U.S. Geological Survey Earthquake Hazards Program decadal science strategy, 2024–33

Executive Summary Earthquakes represent one of our Nation’s most significant and costly natural hazards, with estimated annual loses from earthquakes close to $15 billion in 2023. Over the past two centuries, 37 U.S. States have experienced an earthquake exceeding a magnitude of 5, and 50 percent of States have a significant potential for future damaging shaking; these statistics speak to the need for nationwide interest and investment in earthquake hazard characterization and risk reduction. Authorized under the Earthquake Hazards Reduction Authorization Act, the U.S. Geological Survey (USGS) Earthquake Hazards Program (EHP) provides the scientific information, situational awareness, and knowledge necessary to reduce deaths, injuries, and economic losses from earthquakes and earthquake-induced tsunamis, landslides, and soil liquefaction. The EHP supports activities in three focused topical areas: (1) earthquake monitoring, (2) hazard assessment, and (3) applied research, using the results of each—and the coordination among them—to further support risk translation and communication in regions at risk nationwide. For earthquake monitoring, the Advanced National Seismic System (ANSS), a cooperative effort of USGS networks, university partner regional seismic networks, and real-time geodetic networks, collects and analyzes data on earthquakes; issues timely, reliable notifications of their occurrence and impacts; and provides data for earthquake research, hazard, and risk assessment as a foundation for building an earthquake-resilient Nation. The USGS-operated ShakeAlert Earthquake Early Warning system is a recent addition to EHP’s ANSS infrastructure. In the realm of earthquake hazard assessment, the EHP contributes to earthquake risk mitigation strategies by developing the National Seismic Hazard Model and maps, and other related products, that describe the likelihood and potential effects of earthquakes nationwide, especially in the urban areas of highest risk. The EHP also conducts research on the causes, characteristics, and effects of earthquakes and prioritizes work that directly increases the accuracy and precision of earthquake hazards assessments, earthquake forecasts, and earthquake monitoring and situational-awareness products and that supports the Nation’s earthquake mitigation practices. Bridging the EHP’s efforts across research, hazard assessments, and earthquake monitoring is a broad and comprehensive collection of earthquake information products, including the National Seismic Hazard Model, ShakeAlert, and other products describing impact, such as ShakeMap and PAGER (Prompt Assessment of Global Earthquakes for Response), which have been developed and integrated into EHP’s real-time monitoring systems. EHP funds external partners to carry out many important collaborative activities through an active external grants program—one of the largest in the USGS—and through cooperative agreements with other partners such as the university-operated regional seismic networks, funded as part of the ANSS. To continue its support of earthquake hazard characterization and risk reduction, the EHP aims to strengthen its foundational products and practices while positioning itself to respond to the evolving needs of the Nation and follow best practices of the scientific community. This document describes a strategy for the program to ensure it can meet these demands. The foundational priorities outlined in this strategy represent those activities that remain critical to the core functionality of the program and those that can be supported under current fiscal year 2024-level appropriations. Priorities described as aspirational are important for future growth, and to maintain the program’s position as a leading global resource in earthquake science, but would require increases in appropriated funding to be fully realized. Across the program’s portfolio of activities, several major themes have been identified as the most critical activities to advance EHP science over the coming decade. Together, these activities provide the framework necessary to integrate critical hazard characterization and risk reduction activities across the program. They provide the structure for research to advance the understanding of where, when, and why earthquakes occur and how we can use improved knowledge to drive short-term and actionable forecasts of seismic activity. They expand the usefulness of critical earthquake products and advance the sophistication of those products to keep pace with the rapidly evolving needs of an ever-expanding user base while maintaining the position of the USGS as a global leader in earthquake science. Focus on system-level science. Establish an automated earthquake-processing pipeline. Enhance the accuracy and reliability of the ShakeAlert Earthquake Early Warning system and plan for extension to other regions. Implement time-dependent earthquake forecasting. Develop physically realistic models. Expand computational capacity. This science strategy is organized into three primary sections. The first section provides an overview of the EHP and its budget, governance, and program council. Readers familiar with the program may wish to focus on the second section, which describes the core of the science strategy, including priorities across each of the EHP’s major program activities in monitoring, hazard assessment, and targeted research. The third section outlines science priorities that cut across program activities, including those involving collaborations external to the EHP.

American Samoa, Guam, Northern Marianas Islands, U↗
USGS2024

Preliminary observations of the April 5th, 2024, Mw4.8 New Jersey earthquake

On 5 April 2024, 10:23 a.m. local time, a moment magnitude 4.8 earthquake struck Tewksbury Township, New Jersey, about 65 km west of New York City. Millions of people from Virginia to Maine and beyond felt the ground shaking, resulting in the largest number (>180,000) of U.S. Geological Survey (USGS) “Did You Feel It?” reports of any earthquake. A team deployed by the Geotechnical Extreme Events Reconnaissance Association and the National Institute of Standards and Technology documented structural and nonstructural damage, including substantial damage to a historic masonry building in Lebanon, New Jersey. The USGS National Earthquake Information Center reported a focal depth of about 5 km, consistent with a lack of signal in Interferometric Synthetic Aperture Radar data. The focal mechanism solution is strike slip with a substantial thrust component. Neither mechanism’s nodal plane is parallel to the primary northeast trend of geologic discontinuities and mapped faults in the region, including the Ramapo fault. However, many of the relocated aftershocks, for which locations were augmented by temporary seismic deployments, form a cluster that parallels the general northeast trend of the faults. The aftershocks lie near the Tewksbury fault, north of the Ramapo fault.

New Jersey↗
USGS2024

Constraints from GPS measurements on plate coupling within the Makran subduction zone and tsunami scenarios in the western Indian Ocean

Plate-coupling estimates and previous seismicity indicate that portions of the Makran megathrust of southern Pakistan and Iran are partially coupled and have the potential to produce future magnitude 7+ earthquakes. However, the GPS observations needed to constrain coupling models are sparse and lead to an incomplete understanding of regional earthquake and tsunami hazard. In this study, we assess GPS velocities for plate coupling of the Makran subduction zone with specific attention to model resolution and the accretionary prism rheology. We use finite element model-derived Green's functions to invert for the interseismic slip deficit under both elastic and viscoelastic Earth assumptions. We use the model resolution matrix to characterize plate-coupling scenarios that are consistent with the limited spatial resolution afforded by GPS observations. We then forward model the corresponding tsunami responses at major coastal cities within the western Indian Ocean basin. Our plate-coupling results show potential segmentation of the megathrust with varying coupling from west to east, but do not rule out a scenario where the entire length of the megathrust could rupture in a single earthquake. The full subduction zone rupture scenarios suggest that the Makran may be able to produce earthquakes up to M w 9.2. The corresponding tsunami model from the largest earthquake event ( M w 9.2) estimates maximum wave heights reaching 2–5 m at major port cities in the northern Arabian Sea region. Cities on the west coast of India are less affected (1–2 m). Coastlines bounding eastern Africa, and the Strait of Hormuz, are the least affected (<1 m).

Geophysical Journal International↗
USGS2023

Rapid Source Characterization of the 2023 Mw 6.8 Al Haouz, Morocco, Earthquake

The U.S. Geological Survey (USGS) National Earthquake Information Center (NEIC) estimates source characteristics of significant damaging earthquakes, aiming to place events within their seismotectonic framework. Contextualizing the 8 September 2023, M w 6.8 Al Haouz, Morocco, earthquake is challenging, because it occurred in an enigmatic region of active surface faulting, and low seismicity yet produced significant damage and loss of life. Here, we present the rapid earthquake source products produced by the USGS NEIC, describing how the source model was derived using both seismic and geodetic observations. Our analysis indicates that the earthquake was the result of oblique‐reverse faulting in the lower crust on either a steeply north‐dipping fault or a moderately south‐dipping fault. Finite‐slip models using seismic and geodetic data reveal a compact source, with slip occurring at depths of 15–35 km. The causative fault is not apparent, because the rupture did not break the surface, and it is not possible to definitively attribute the earthquake to a known structure. The earthquake centroid depth of 25 km is noteworthy, because it shows slip extending beyond common estimates of seismogenic depth. This earthquake highlights that the seismogenic processes associated with mountain building in this wide plate boundary region are poorly understood.

The Seismological Record↗
USGS2023

Rapid surface rupture mapping from satellite data: The 2023 Kahramanmaraş, Turkey (Türkiye), earthquake sequence

The 6 February 2023 Kahramanmaraş, Turkey (Türkiye), earthquake sequence produced > 500 km of surface rupture primarily on the left‐lateral East Anatolian (~345 km) and Çardak (~175 km) faults. Constraining the length and magnitude of surface displacement on the causative faults is critical for loss estimates, recovery efforts, rapid identification of impacted infrastructure, and fault displacement hazard analysis. To support these efforts, we rapidly mapped the surface rupture from satellite data with support from remote sensing and field teams, and released the results to the public in near‐real time. Detailed surface rupture mapping commenced on 7 February and continued as high‐resolution (< 1.0 m/pixel) optical images from WorldView satellites (2023 Maxar) became available. We interpreted the initial simplified rupture trace from subpixel offset fields derived from Advanced Land Observation Satellite2 and Sentinel‐1A synthetic aperture radar image pairs available on 8 and 10 February, respectively. The mapping was released publicly on 10 February, with frequent updates, and published in final form four months postearthquake ( Reitman, Briggs, et al. , 2023 ). This publicly available, rapid mapping helped guide fieldwork and constrained U.S. Geological Survey finite‐fault and loss estimate models, as well as stress change estimates and dynamic rupture models.

Kahramanmaraş↗
USGS2023

Rapid characterization of the February 2023 Kahramanmaraş, Turkey, earthquake sequence

The 6 February 2023 M w "> M w 7.8 Pazarcık and subsequent M w "> M w 7.5 Elbistan earthquakes generated strong ground shaking that resulted in catastrophic human and economic loss across south‐central Türkiye and northwest Syria. The rapid characterization of the earthquakes, including their location, size, fault geometries, and slip kinematics, is critical to estimate the impact of significant seismic events. The U.S. Geological Survey National Earthquake Information Center (NEIC) provides real‐time monitoring of earthquakes globally, including rapid source characterization and impact estimates. Here, we describe the seismic characterization products generated and made available by the NEIC over the two weeks following the start of the earthquake sequence in southeast Türkiye, their evolution, and how they inform our understanding of regional seismotectonics and hazards. The kinematics of rupture for the two earthquakes was complex, involving multiple fault segments. Therefore, incorporating observations from rupture mapping was critical for characterizing these events. Dense local datasets facilitated robust source characterization and impact assessment once these observations were obtained and converted to NEIC product input formats. We discuss how we may improve the timeliness of NEIC products for rapid assessment of future seismic hazards, particularly in the case of complex ruptures.

The Seismic Record↗
USGS2023

Rapid shallow megathrust afterslip from the 2021 M8.2 Chignik, Alaska earthquake revealed by seafloor geodesy

The shallower portions of subduction zone megathrust faults host Earth’s most hazardous tsunamigenic earthquakes, yet understanding how and when they slip remains elusive because of challenges making seafloor observations. We performed Global Navigation Satellite System Acoustic seafloor geodetic surveys before and ~2.5 months after the 29 July 2021 M w (moment magnitude) 8.2 Chignik, Alaska, earthquake and determine ~1.4 meters cumulative co- and post-seismic horizontal displacement ~60 kilometers from the megathrust front. Only for the 2011 M w 9 Tohoku event have closer subduction zone earthquake displacements been observed. We estimate ~2 to 3 meters of megathrust afterslip shallower than 20 kilometers, a portion of the megathrust on which both inter- and co-seismic slip likely had occurred previously. Our analysis demonstrates that by 2.5 months, shallower and deeper moment had effectively equilibrated on the megathrust, suggesting that its tsunamigenic potential remains no more elevated than before the earthquake.

Alaska↗
USGS2022

Power-law viscoelastic flow of the lower accretionary prism in the Makran subduction zone following the 2013 Baluchistan Earthquake

Subduction zone accretionary prisms are commonly modeled as elastic structures where permanent deformation is accommodated by faulting and folding of otherwise elastic materials, yet accretionary prisms may exhibit other deformation styles over relatively short time scales. In this study, we use 6.5-year (2014–2021) Sentinel-1 interferometric synthetic aperture radar (InSAR) time-series of post-seismic deformation in the Makran accretionary prism of southeast Pakistan to characterize non-linear viscoelastic deformation within an active accretionary prism on short timescales (months to years). We constructed a series of 3-D finite-element models of the Makran subduction zone, including an accretionary prism, and constrained the elastic thickness of the upper wedge and the flow-law parameters (power-law exponent, activation enthalpy, and pre-exponential constant) of the lower wedge through forward model fits to the InSAR time-series. Our results show that the prism is elastically thin (8–12 km) and the non-linear viscoelastic relaxation of the deep portions of the prism alone can sufficiently explain the post-seismic surface deformation. Our best fitting flow-law parameters ( n = 3.76 ± 0.39, Q = 82.2 ± 37.73 kJ mol −1 , and A = 10 −3.36±4.69 ) are consistent with triggering of low temperature dislocation creep within fluid-saturated siliciclastic rocks. We believe that the fluids necessary for this weakening originate from sedimentary underplating and/or the presence the hydrocarbons. The presence of power-law rheology within the lower wedge impacts the estimated plate coupling and the stress state in the subduction system, with respect to the conventional elastic wedge model, and hence should to be considered in future earthquake cycle models.

Makran accretionary prism↗
USGS2022

The Geodetic Centroid (gCent) Catalog: Global earthquake monitoring with satellite imaging geodesy

Remote sensing geodetic observations (Interferometric Synthetic Aperture Radar [InSAR] and optical correlation [“pixel tracking”]) serve an increasingly diverse and important role in earthquake monitoring and response. This study introduces the Geodetic Centroid (gCent) catalog—an earthquake catalog derived solely from space‐based geodetic observations—and analysis of 74 earthquakes ( ⁠ M W 4.3–7.4) imaged from 1 August 2019 to 01 February 2022. For gCent, we use InSAR and optical correlation observations derived from the Sentinel‐1 satellites and various publicly available optical satellites to systematically image all global earthquakes M W 5.5 or larger and shallower than 25 km, M W 7.0 or larger at any depth, and other high‐impact earthquakes or seismic events of special interest. We invert surface displacements from successfully imaged earthquakes for the location, orientation, and dimensions of a single slipping fault patch that describes the centroid characteristics of the earthquake. These centroid models, in turn, are compiled into a catalog and used in U.S. Geological Survey/Advanced National Seismic System (ANSS) operational earthquake response products such as ShakeMaps and finite‐fault models. We provide a comparison of the gCent catalog to the ANSS Comprehensive Catalog and Global Centroid Moment Tensor (Global CMT) catalog to compare reported locations, depths, and magnitudes. We find that global earthquake catalogs not only generally provide reasonably comparable locations (within 10 km on average), but also they systematically overestimate depth that may have implications for earthquake shaking predictions based solely on earthquake origin information. Geodetic magnitudes are comparable to seismically inferred magnitudes, indicating that gCent models are unlikely to be systematically biased by the presence of postseismic deformation. We additionally highlight limitations of the gCent catalog induced by both the limitations of remote sensing imaging of earthquakes and our imposition of a simplified earthquake source description that does not include spatially distributed slip.

Bulletin of the Seismological Society of America↗
USGS2021

Coseismic surface displacement in the 2019 ridgecrest earthquakes: Comparison of field measurements and optical image correlation results

A fundamental topic in earthquake studies is understanding the extent to which fault rupture at the surface is localized on primary fault strands as opposed to distributed tens to hundreds of meters away from primary ruptures through off‐fault deformation (OFD) via a combination of discrete secondary faulting and bulk deformation. The 2019 Ridgecrest, CA M w 6.4 and M w 7.1 earthquakes provide an opportunity to explore this problem via comparison of published field‐based and mostly on‐fault offset measurements, with new lateral displacement measurements made over length scales of hundreds of meters across the primary rupture using WorldView satellite images collected before and after the earthquakes. A mean fit to the field observations underestimates net slip (pixel‐correlation results) by an average of 0.4 m along the M w 6.4 rupture (∼41% of net left‐lateral displacement) and 0.7 m along the M w 7.1 rupture (∼65% of net right‐lateral displacement). We attribute these differences to substantial OFD along the lengths of the M w 6.4 (∼59%) and M w 7.1 (∼35%) ruptures. A maximum fit to the field observations provides an improved match to the pixel correlation results (difference of 0.1–0.2 m or ∼76%–98% of net displacement). OFD may in part depend on displacement magnitude, where locations along the rupture with smaller displacements (<2 m) are characterized by a greater percentage of OFD. Furthermore, OFD increases with rupture zone width, especially at the scale of individual rupture strands (hundreds of meters). These findings contribute to a growing understanding of the important role of OFD in accommodating deformation near the surface.

California↗
USGS2020

Localized fault-zone dilatancy and surface inelasticity of the 2019 Ridgecrest earthquakes

Earthquakes produce a spectrum of elastic and inelastic deformation processes that are reflected across various length and time scales. While elasticity has long dominated research assumptions in active tectonics, increasing interest has focused on the inelastic characteristics of earthquakes, particularly those of the surface fault rupture zone itself, and how they relate to ground rupture hazard and the mechanics of damage zones. Here we present detailed co-seismic surface-strain analysis of the 2019 Ridgecrest, California, earthquakes. We derive three-dimensional high-resolution surface displacements from satellite optical imagery, which we then invert for the co-seismic surface-strain tensors. We show that fault-zone dilation is pervasive throughout these earthquakes and that inelastic failure is present but relatively localized (median width of 31 m). The width of the inelastic failure zone is not correlated to off-fault deformation, surface geology or displacement magnitude. Instead, the extent and kinematics of inelastic failure reflect active, mylonitic deformation of the fault damage zone that is influenced by rupture velocity and fault maturity. These results highlight how a single earthquake contributes to the long-term, permanent geologic record of faulting.

California↗
USGS2019

Surface rupture and distributed deformation revealed by optical satellite imagery: The intraplate 2016 Mw 6.0 Petermann Ranges earthquake, Australia

High-resolution optical satellite imagery is used to quantify vertical surface deformation associated with the intraplate 20 May 2016 Mw 6.0 Petermann Ranges earthquake, Northern Territory, Australia. The 21 ╓ 1 km long NW-trending rupture resulted from reverse motion on a northeast-dipping fault. Vertical surface offsets of up to 0.7 ╓ 0.1 m distributed across a 0.5-to-1 km wide deformation zone are measured using the Iterative Closest Point (ICP) algorithm to compare pre- and post-earthquake digital elevation models (DEMs) derived from Worldview imagery. The results are validated by comparison with field-based observations and interferometric synthetic aperture radar (InSAR). The pattern of surface uplift is consistent with distributed shear above the propagating tip of a reverse fault, leading to both an emergent fault and folding proximal to the rupture. This study demonstrates the potential for quantifying modest (<1 m) vertical deformation on a reverse fault using optical satellite imagery.

Geophysical Research Letters↗
USGS2019

Vertical coseismic offsets from differential high-resolution stereogrammetric DSMs: The 2013 Baluchistan, Pakistan earthquake

The recent proliferation of high-resolution (< 3-m spatial resolution) digital topography datasets opens a spectrum of geodetic applications in differential topography, including the quantification of coseismic vertical displacement fields. Most investigations of coseismic vertical displacements to date rely, in part, on pre- or post-event lidar surveys that are intractable or non-existent in many locales. Stereogrammetric digital surface models (DSMs) derived from high-resolution satellite optical imagery provide a new avenue for the retrieval of spatially-dense vertical coseismic displacements on a global scale. In this study, we generated 2-m resolution pre- and post-seismic DSMs from satellite optical imagery spanning the 2013 Mw7.7 Baluchistan strike-slip earthquake that occurred on the Hoshab fault in southern Pakistan. We applied the Iterative Closest Point algorithm to the DSMs to quantify the coseismic vertical displacement field at a spatial resolution of 10-30 m and to generate 3D coseismic strain tensors. We found that across-fault vertical offsets alternated between uplift and subsidence and varied between ~1-3 m in a non-systematic manner along the Hoshab fault. We show that the pre-existing topography and near-fault geomorphology are variably consistent and inconsistent with the displacement kinematics of the 2013 earthquake, and we argue that these relationships highlight varied slip sense history along the Hoshab fault. Notably, topography along the southern extents of the Hoshab fault requires different surface displacement kinematics than occurred in the 2013 earthquake, suggesting that the Hoshab fault accommodates varying senses of slip (bimodal slip) through time.

Baluchistan Province↗