This website uses cookies. By continuing to use this website you are agreeing to our use of cookies. 



3D Displacements and Strain from the 2023 February Turkey Earthquakes, version 1

Latest Data Update: 2023-03-10
Status: Completed
Online Status: ONLINE
Publication State: Citable
Publication Date: 2023-03-03
DOI Publication Date: 2023-03-14
Download Stats: last 12 months
Dataset Size: 73 Files | 5GB


This data set contains strain, motion magnitude, 3D displacements, and surface slip distributions from the 2023 February Türkiye Earthquakes. COMET presents 0.001 degree (~100 m) resolution 3D displacement fields, with associated uncertainties, jointly inverted from 4 tracks of Sentinel-1 range and azimuth offsets and a set of north and east displacements from Sentinel-2 pixel tracking. From the 3D displacement fields, we calculate 0.001-degree resolution horizontal motion magnitudes, from which we extract surface slip distribution along the two faults ruptured during the Mw7.8 and Mw7.5 earthquakes. Further calculation of a 0.01-degree resolution strain field from the east and north displacement fields highlights the surface ruptures caused by the two earthquakes.

Citable as:  Ou, Q.; Lazecky, M.; Watson, C.S.; Maghsoudi, Y.; Wright, T. (2023): 3D Displacements and Strain from the 2023 February Turkey Earthquakes, version 1. NERC EDS Centre for Environmental Data Analysis, 14 March 2023. doi:10.5285/df93e92a3adc46b9a5c4bd3a547cd242.
Abbreviation: Not defined
Keywords: COMET, Sentinel 1, Sentinel 2


Previous Info:
No news update for this record
Previously used record identifiers:
No related previous identifiers.
Access rules:
Public data: access to these data is available to both registered and non-registered users.
Use of these data is covered by the following licence: When using these data you must cite them correctly using the citation given on the CEDA Data Catalogue record.
Data lineage:

Data were produced by the project team and supplied for archiving at the Centre for Environmental Data Analysis (CEDA).
Data quality: Quantitative offset and Cartesian motion uncertainties are provided with the data. The northern half of the scene is noisier than the south due to snow cover. The Cartesian motions have mode uncertainties of 0.03 m for the east and up components and 0.06 m for the north component. The median and mean uncertainties of the 3 components are ~0.3 m and ~1 m, respectively.

Data Quality:
Data are as given by the data provider, no quality control has been performed by the Centre for Environmental Data Analysis (CEDA).
File Format:
These data are provided in GeoTIFF (.tif) or png format.

Process overview

This dataset was generated by a combination of instruments deployed on platforms and computations as detailed below.

Instrument/Platform pairings

Sentinel 2 Multispectral Instrument (MSI) Deployed on: Sentinel 2B

Instrument/Platform pairings

Sentinel 1 Synthetic Aperture Radar (SAR) Deployed on: Sentinel 1A

Instrument/Platform pairings

Sentinel 2 Multispectral Instrument (MSI) Deployed on: Sentinel 2A

Computation Element: 1

Title Computation for the 3D displacements and strain from the 2023 February Turkey Earthquakes
Abstract The Sentinel-1 azimuth and range offsets came from the COMET-LiCS portal ( and are processed using Gamma package from Level 1 (L1) data from the European Space Agency ( The offsets are from 4 tracks of interferometric pairs which are provided in links in the docs tab. These are results of processing by GAMMA command offset_pwr_tracking, cross-correlating 128x64 pixel windows (range x azimuth) over 2x oversampled deramped low-pass filtered intensity data, selecting data of cross-corr. coeff. over 0.1. Upsampled by linear interpolation. The optical pixel tracking east and north offsets are derived from the Sentinel-2 data using the ENVI COSI-Corr plugin and are also provided in the docs tab. Images: L1C 25th Jan and 9th Feb. Parameters: window size of 64 down to 32 with step of 4. 2 iterations, 0.9 mask threshold, no resampling and gridded output. All the offset data are referenced to a distribution of dummy zero points away from the coseismic ruptures by removing a planar ramp. The uncertainties of the offset data are empirically estimated from the offset data themselves by evaluating the absolute deviation of each pixel value from its local mean averaged across a 4x4 pixel window, assuming nan-values are zeros. These uncertainties are used to weight the 3D motion inversion and are propagated to the uncertainties of the decomposed displacements through a model covariance matrix. The motion magnitude field is a vector combination of the east and north motion fields, each masked by respective uncertainties. The arrows represent average east and north motions from 30 km windows evaluated at 0.2 degree intervals. The second invariant of horizontal strain is calculated from the horizontal displacement gradients of east and north motions after median filter with 30 km windows at 0.01 degree intervals.
Input Description None
Output Description None
Software Reference None
Output Description


Output Description


Output Description


No variables found.

Temporal Range
Start time:
End time:
Geographic Extent