USGS Science⌕ Search

USGS · 70283506

Rapid building damage estimates from the M7.8 Turkey earthquake sequence via causality-informed Bayesian inference from satellite imagery

Abstract

On February 6, 2023, a major earthquake of 7.8 magnitude and its aftershocks caused widespread destruction in Turkey and Syria, causing more than 55,000 deaths, displacing 3 million people in Turkey and 2.9 million in Syria, and destroying or damaging at least 230,000 buildings. Our research presents detailed city-scale maps of landslides, liquefaction, and building damage from this earthquake, utilizing a novel variational causal Bayesian network. This network integrates InSAR-derived change detection with new empirical ground failure models and building footprints, enabling us to (1) rapidly estimate large-scale building damage, landslides, and liquefaction from remote sensing data, (2) jointly attribute building damage to landslides, liquefaction, and shaking, (3) improve regional landslide and liquefaction predictions impacting infrastructure, and (4) simultaneously identify damage degrees in thousands of buildings. For city-scale, building-by-building damage assessments, we use building footprints and satellite imagery with a spatial resolution of approximately 30 meters. This allows us to achieve a high resolution in damage assessment, both in timeliness and scale, enabling damage classification at the individual building level within days of the earthquake. Our findings detail the extent of building damage, including collapses, in Hatay, Osmaniye, Adıyaman, Gaziantep, and Kahramanmaras. We classified building damages into five categories: no damage, slight, moderate, partial collapse, and collapse. We evaluated damage estimates against preliminary ground-truth data reported by the civil authorities. Our results demonstrate the accuracy of our classification system, as evidenced by the area under the curve (AUC) scores on the receiver operating characteristic (ROC) curve, which ranged from 0.9588 to 0.9931 across different damage categories and regions. Specifically, our model achieved an AUC of 0.9931 for collapsed buildings in the Hatay/Osmaniye area, indicating a 99.31% probability that the model will rank a randomly chosen collapsed building higher than a randomly chosen non-collapsed building. These accurate, building-specific damage estimates, with greater than 95% classification accuracy across all categories, are crucial for disaster response and can aid agencies in effectively allocating resources and coordinating efforts during disaster recovery.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 35.6328999° to 42.4361248° latitude; 25.0394104° to 44.5715137° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Xuechun Li, Xiao Yu, Paula Madeline Burgi, David J. Wald, Xie Hu, Susu Xu. 2025. Rapid building damage estimates from the M7.8 Turkey earthquake sequence via causality-informed Bayesian inference from satellite imagery. https://doi.org/10.1177/87552930241290501

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

An examination of the reasons for, and variability of, the differences between response spectra and pseudo response spectra

For many decades, it has been common practice to overlook the differences between response spectra and pseudo response spectra. This practice stems from analyses indicating that the average difference between response spectra and pseudo response spectra is negligible for a viscous damping ratio of 5%, which is near-universally adopted. As others have shown, the differences depend strongly on the damping ratio, earthquake magnitude, and oscillator period, with larger differences for larger damping, smaller magnitudes, and longer periods. To quantify the full statistical distribution of the ratio of the absolute acceleration response spectrum (SA) to the pseudo absolute acceleration response spectrum (PSA), we fit empirical ratios with the generalized gamma distribution. Data for magnitudes from 4.0 to 5.5 give 95% prediction intervals for the ratio of 1.00–1.02 at 1 s and 1.01–2.24 at 10 s for 5% damping. For 30% damping, these intervals increase to 1.14–2.01 at 1 s and 1.62–11.8 at 10 s, respectively. Further, we demonstrate how the shape of the ground-motion Fourier amplitude spectrum interacts with the natural period of the oscillator to control the differences in the response spectra type. Routinely treating the two measures of acceleration response spectra as interchangeable may lead to unacceptable inconsistencies when the response spectra of individual records are important or when larger damping ratios are used, such as in direct displacement-based design. If one must convert between them, the uncertainty of the conversion should be accounted for during analysis.

Earthquake Spectra↗

Accounting for biases in the analysis of building damage data for the 2023 M7.8 Türkiye/Syria earthquake sequence

In the aftermath of a large earthquake, several practical constraints affect the systematic collection of data for assessing building damage. The available data are therefore typically incomplete, frequently missing key information such as the type of construction, age and geographic location. Moreover, the observations often favour certain locations and damage grades. In the absence of robust methods, using such data directly for characterising the building stock or conducting fragility assessments can introduce systematic bias. Using the 2023 M7.8 Türkiye/Syria earthquake sequence as a case study, we propose a statistical model that enables joint inference over the inventory and fragility of the building stock while explicitly accounting for sources of data error and missing information. The results show higher-than-expected vulnerability amongst several building types, particularly mid- and high-rise reinforced concrete buildings, which performed extremely poorly during the 2023 earthquake sequence. The proposed Bayesian framework provides rigorous uncertainty quantification that can be propagated to future applications.

Earthquake Spectra↗

Suite of regional seismic velocity models for the U.S. Atlantic and Gulf Coastal Plains (CPVMv2.0): Incorporating epistemic uncertainty in velocity models

The U.S. Atlantic and Gulf Coastal Plain (CP) sediments are kilometers-thick, low-velocity deposits underlain by high-velocity bedrock that create regional patterns of amplification and attenuation of ground motions. Capturing these effects is challenging because of the limited subsurface information available, particularly the shear wave velocity ( V S ) structure down to bedrock (or assumed reference rock condition) and the associated uncertainties. This paper introduces a new suite of regional velocity models for the CPs, called CPVMv2.0, which builds upon previous efforts integrating geology, V S measurements, and sediment thickness in the region. CPVMv2.0 uses 750 measured V S profiles and identifies nine surficial-geologic-age-based and geography-based groups (with alternative physiographic-province-based subgroupings) to create surface-based median V S profiles. This suite of new models also includes six subsurface-based median profiles developed using geologic age boundaries from depth-dependent stratigraphic sequences extracted from the U.S. Geological Survey National Crustal Model. In addition to quantifying aleatory variability in sediment velocity, weathered-rock velocity, bedrock velocity, and the depth to the top of weathered rock and bedrock, CPVMv2.0 provides alternative profile models that allow for the characterization of the epistemic uncertainty. This is a unique and relevant feature of CPVMv2.0. Our suite of models was validated using measured and estimated fundamental resonant frequencies as well as other V S models in the region. CPVMv2.0 enables the incorporation of crustal velocities (and their associated variability and uncertainty) into physics-based ground-motion simulations in the eastern United States and into the development of new ground-motion models aiming to characterize systematic site effects via V S30 , depths to V S isosurfaces, and/or sediment thickness in the region.

U.S. Atlantic and Gulf Coastal Plains↗