Preliminary documentation of coseismic ground failure triggered by the February 6, 2023 Türkiye earthquake sequence

Tolga Gorum; Hakan Tanyas; Furkan Karabacak; Abdussamet Yilmaz; Serkan Girgin; Kate E. Allstadt; M. Lutfi Suzen; Paula Madeline Burgi

doi:10.1016/j.enggeo.2023.107315

USGS · 70256619

Preliminary documentation of coseismic ground failure triggered by the February 6, 2023 Türkiye earthquake sequence

Tolga Gorum·Hakan Tanyas·Furkan Karabacak·Abdussamet Yilmaz·Serkan Girgin·Kate E. Allstadt·M. Lutfi Suzen·Paula Madeline Burgi

View original source ↗Export metadata ↓

Abstract

The devastating Kahramanmaraş earthquake sequence occurred on February 6, 2023. Two main events, M w 7.8 and M w 7.5 occurred 9 hours apart, affected 11 cities in Turkey, and subjected an area of ∼90,000 km 2 to shaking levels known to trigger landslides (peak ground acceleration > 0.08 g). Extensive landsliding was expected given the hilly terrain affected by this significant ground shaking—about 15% of the topography is steeper than 20°—but was not initially apparent in early satellite imagery, mostly because of obscuring snow that fell just after the earthquakes. However, after a more detailed investigation using high-resolution satellite images, aerial photos, and a field survey, we confirmed that this earthquake sequence did, indeed, trigger numerous landslides. In this study, we present those findings and provide a preliminary characterization of the spatial distribution, general characteristics, and dominant types of landslides and hillslope deformation triggered by the earthquake sequence. We mapped 3673 coseismic landslides, mostly concentrated in the northern half of the impacted area. Rock falls are the most abundant landslide type, but bedrock rotational landslides, translational slides and lateral spreads are also numerous. Surface rupture through mountainous terrain caused several large, and in some cases fatal, landslides. Incipient landslides and ground cracks are also widespread, especially in the north. Lithology, spatial variability of ground shaking, and topographic relief appear to be the main variables controlling the spatial distribution of coseismic landslides. There are few detailed studies of earthquake-triggered landslides in arid and semi-arid regions such as this one, nor for such complex earthquake sequences. Therefore, this contribution provides valuable information for future hazard and modeling efforts in arid and semi-arid regions.

Explore related subjects

Engineering Geology Article

Source-reported bounding extent: 35.973045074631955° to 38.71377557689037° latitude; 35.469730632443316° to 39.69683217133675° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Tolga Gorum, Hakan Tanyas, Furkan Karabacak, Abdussamet Yilmaz, Serkan Girgin, Kate E. Allstadt, M. Lutfi Suzen, Paula Madeline Burgi. 2023. Preliminary documentation of coseismic ground failure triggered by the February 6, 2023 Türkiye earthquake sequence. https://doi.org/10.1016/j.enggeo.2023.107315

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

KEEP EXPLORING

Related USGS reports

USGS2026

Sackung at Bald Eagle ridge, central Colorado: An updated interpretation of ridge-spreading movement, structures, and mechanisms from 50 years (1975–2025) of U.S. Geological Survey research

Slow gravitational failures of mountain peaks and ridges are poorly understood. Herein, we report on 50 years of studies at a slowly spreading castellate ridge in the Sawatch Range in central Colorado. The orientations of geomorphic-structural features indicate that the fractured Precambrian granitic rock underlying the ridge has extended and spread northwestward toward the formerly glacier-covered Busk Creek valley. Results from surveying, field-based geomorphic-structural mapping using lidar, rock mass quality measurements, a passive seismic survey, and satellite radar provide a major update to research started by U.S. Geological Survey researchers in the 1970s and 1980s. New insights include a recognition that the entire ridge has slowly moved by concurrent sliding along an inferred northwest dipping, compound basal-slip surface (or zone), and through the formation of multiple grabens by normal faulting and flexural toppling along sets of pre-existing fractures that dip moderately (∼45°) to the southeast. We were unable to distinguish the presence of a sudden and strong contrast in seismic velocity across the inferred slip surface. Movement during the 50-year study period has been episodic and gradually decreasing, in correspondence with decreasing cumulative annual precipitation and increasing mean annual air temperatures. The fastest moving area, just upslope from the glacier trimline, had an average horizontal velocity of 3 – 4 mm/yr. Evidence suggests that movement started as a paraglacial response mechanism, but because of the site’s proximity to the Rio Grande Rift, we cannot exclude earthquake shaking as a mechanism for initiation or enhancement of slope movement. An estimate of longer-term horizontal movement from the exposed basal-slip surface at the uphill side of the ridgetop graben is ∼1.1 mm/yr for the 13 – 14 ky post-glacial period. Broad implications of our work are that: (1) long-term measurements (decades or longer) of slope movement can add insights into how sackungen form and evolve through time; (2) the identification of thrust faults and toes in zones of compression near valley bottoms can be crucial for interpreting sackung failure mechanisms, and (3) the use of passive-seismic techniques to identify the depth to a slip surface may not be successful in granitic terrain dominated by planar fractures and subtle changes in rock-mass characteristics.

Colorado↗

USGS2026

From landslide susceptibility to risk assessment in the conterminous U.S.

Understanding the spatial distribution of landslide prone-areas and what consequences they may have is important for risk management and land-use planning. In the United States, although landslides occur in every state, a comprehensive landslide risk assessment is still missing. Existing efforts, such as the Federal Emergency Management Agency (FEMA)’s National Risk Index, rely on aggregated products and coarse cartographic units, limiting their geomorphological and practical accuracy. In this study, we present a methodological advance for landslide risk assessment across large areas with incomplete and sparse data. We apply our procedures to the conterminous United States by integrating geomorphologically meaningful partitions and spatial and temporal probability data-driven models. Landslide susceptibility is estimated using a Generalized Additive Mixed Model incorporating a bias capture/correction scheme to account for inventory inaccuracies (reference Area Under the Curve = 0.75). The exceedance probabilities of landslide occurrence are defined for three temporal scenarios (2, 5, and 10 year). Then, we explore the associated potential economic consequences for human settlements and agricultural areas. The findings indicate that the spatial variability of risk is primarily controlled by exposure rather than by susceptibility/hazard alone. The mean risk increases by ∼170% from the 2-year to the 10-year scenario. Beyond its quantitative outcomes, this study offers a blueprint for continental or sub-continental scale landslide risk assessments, demonstrating both the opportunities and current limitations.

Engineering Geology↗

USGS2025

Global patterns of coseismic landslide runout mobility differ from aseismic landslide trends

Coseismic landslides significantly contribute to human and economic losses during and immediately following earthquakes, yet very little data on the runout of such landslides exist. While well-established behavior of aseismic (e.g., hydrologically triggered) landslide runout mobility suggests strong correlation between landslide size and mobility, limited studies of coseismic landslide runout find conflicting mobility trends. We present a global dataset of runout lengths produced from a new automated method for estimating landslide runout, developed and validated using 1726 manually mapped landslides from five unique earthquakes. We then apply the automated runout tool to 23 global earthquake-induced landslide inventories, producing a compiled database of 73,665 measured and estimated runout lengths of coseismic landslides to assess mobility trends. We find a significant divergence between well-established aseismic mobility trends and that of coseismic landslides, with far greater scatter and more complex mobility patterns in earthquake-triggered landslides. As a function of landslide size, we observe global coseismic landslide mobility patterns are bilinear, becoming increasingly less mobile with increasing size above some threshold. This discordance between aseismic and coseismic landslide mobility may be a function of landslide type, kinematics, hydrology, and or setting that systematically differ between triggering mechanisms and should be explored in more depth to develop predictive models of these unique runout patterns. These results suggest hazard and risk models for coseismic landslides may significantly under-predict or over-predict impacts, depending on the size of triggered landslides.

Engineering Geology↗