The Sentinel-1 azimuth and range offsets came from the COMET-LiCS portal (https://comet.nerc.ac.uk/comet-lics-portal/) and are processed using Gamma package from Level 1 (L1) data from the European Space Agency (https://scihub.copernicus.eu/dhus/#/home). 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.