remote sensing Three-Dimensional Surface Displacements of the 8 January 2022 Mw6.7 Menyuan Earthquake, China from Sentinel-1 and ALOS-2 SAR Observations

: The 8 January 2022 Mw6.7 Menyuan earthquake was generated in the transition zone between the western Lenglongling fault and the eastern Tuolaishan fault, both being part of the Qilian–Haiyuan fault system with an important role in the adjustment of the regional tectonic regime. In this study, four pairs of SAR (synthetic aperture radar) data from Sentinel-1 and ALOS-2 (Advanced Land Observation Satellite-2) satellites were used to derive the surface displacement observations along the satellite line-of-sight (LOS) and azimuth directions using the differential interferometric SAR (InSAR, DInSAR), pixel offset-tracking (POT), multiple aperture InSAR (MAI), and burst overlap InSAR (BOI) methods. An SM-VCE method (i.e., a method for measuring three-dimensional (3D) surface displacements with InSAR based on a strain model and variance component estimation) was employed to combine these derived SAR displacement observations to calculate the 3D co-seismic displacements. Results indicate that the 2022 Menyuan earthquake was dominated by left-lateral slip, and the maximum horizontal and vertical displacements were 1.9 m and 0.6 m, respectively. The relative horizontal surface displacement across the fault was as large as 2–3 m, and the fault-parallel displacement magnitude was larger on the southern side of the fault compared with the northern side. Furthermore, three co-seismic strain invariants were also investigated, revealing that the near-fault area suffered severe deformation, and two obviously expanding and compressed zones were identiﬁed. We provide displacements/strains derived in this study in the prevailing geotiff format, which will be useful for the broad community studying this earthquake; in addition, the SM-VCE code used in this study is open to the public so that readers can better understand the method. Conceptualization, J.L. and J.H.; methodology, J.L. and J.S.; software, J.L.; validation, J.L. and Z.M.; formal analysis, J.L. and Z.M.; investigation, J.S.; resources, Z.L. and W.X.; writing—original draft preparation, J.L., Z.M. and J.S.; writing—review and editing, J.H., Z.M. Q.S.; visualization, J.L. and J.S.; supervision, J.H. and Q.S.; project administration, J.L.;


Introduction
At 1:45:30 a.m. on 8 January 2022 (local time), an Mw = 6.7 earthquake occurred about 54 km northwest of Menyuan county, Qinghai province, China with a focal depth of about 10 km [1,2]. This earthquake was generated in the transition zone between the western Lenglongling fault (LEF) and the eastern Tuolaishan fault (TF), both being part of the Qilian-Haiyuan fault system (QHFS) (Figure 1). The QHFS comprises from west to east the following structures: the TF, LEF, Jinqianghe fault (JF), Maomaoshan fault (MF), Laohushan fault (LAF), and Hauyuan fault (HF) [3]. This fault system is generally characterized by left-lateral slip and is one of the most important fault systems in the northeastern boundary of the Qinghai-Tibet Plateau block [4]. In this fault system, the northeastward movements of the Qinghai-Tibet Plateau block are resisted by the Gebi-Alashan block and the Huabeikelatong block (see Figure 1), and the tectonic movement direction changes clockwise slowly [5][6][7]. The left-lateral QHFS can just compensate for these complicated tectonics. Within this fault system, the segment between the LEF, JF, MF, and LAF is considered as a seismic gap and requires more attention [8,9]. Although two Mw > 6.0 events occurred near the LEF (see Figure 1), their fault slip mechanisms were obviously inconsistent with the overall left-lateral slip of the QHFS [5]. Investigations show that these two events may have been related to the northern Lenglongling fault (NLEF), which is an associated thrust fault located in the northern part of the western LEF [5,10]. According to the field investigation and geodesy observations [2,11], this 2022 Menyuan earthquake was a typical left-lateral seismic event, which is consistent with the overall pattern of the QHFS. Therefore, it is significant to swiftly obtain the ground surface displacement of this earthquake, thus providing reliable observations for the fault mechanism analysis and the overall seismic risk assessment of the QHFS.
Remote Sens. 2022, 14, x FOR PEER REVIEW 3 of 15 In this study, we focus on estimating a complete and accurate 3D co-seismic displacement field related to the 2022 Menyuan earthquake on the basis of the ascending/descending Sentinel-1 SAR data and descending ALOS-2 SAR data. The DInSAR, POT, MAI, and BOI [15] methods are all used to measure the ground surface displacements along the lineof-sight (LOS) or azimuth directions. The SM-VCE method was applied here to combine the DInSAR-, MAI-, POT-, and BOI-derived displacement observations to obtain 3D coseismic displacements of the 2022 Menyuan earthquake. Since more SAR displacement observations are used and an advanced method (i.e., SM-VCE) is applied to calculate 3D displacements, the obtained 3D displacements can accurately reveal the co-seismic deformation of this earthquake and provide a reliable dataset for further model interpretation. Moreover, the co-seismic strains were also calculated on the basis of the 3D displacements. In particularly, to benefit more researchers interested in co-seismic surface deformation assessment and the SM-VCE method, the obtained displacement data in geotiff format and the used SM-VCE code can be freely downloaded (see Data Availability Statement). Qinghai-Tibet Plateau block [6,7]. Blue squares are Qilian county and Haiyuan county.

The Used SAR Data
Error! Reference source not found. presents the basic information of the SAR data used in this paper. Their coverages can be found in Figure 1. All of these SAR data were acquired at several days after the earthquake; therefore, it is inevitable to include the postseismic deformations in the SAR-derived displacement observations, which is normal for SAR earthquake studies [22,23]. We also investigated the post-seismic displacement in the descending Sentinel-1 pair 20220110-20220122, and we found that the post-seismic displacement of this earthquake was negligible compared with the co-seismic displacement; therefore, it was practical to combine these SAR-data-derived observations to obtain a 3D displacement field even though the SAR acquisition dates were inconsistent.  [6,7]. Blue squares are Qilian county and Haiyuan county.
Benefiting from the open Sentinel-1 synthetic aperture radar (SAR) data with regular acquisitions, scientists soon obtained the ascending/descending Sentinel-1 SAR image pairs on 10 January 2022, 2 days after the earthquake, providing a data basis for revealing the surface displacement of this earthquake. In particular, Li et al. [2] quickly obtained the co-seismic displacements along the satellite line-of-sight (LOS) direction from the ascending/descending Sentinel-1 SAR images using the differential interferometric SAR (InSAR, DInSAR) method, and then investigated the source parameters and slip distributions of the 2016 and 2022 Menyuan earthquakes. Moreover, Yang et al. [11] used the same SAR data and displacement observations to construct a finite fault model of this event, indicating that multiple fault segments ruptured during this earthquake. However, the used ascending SAR image pair only covers part of the deforming zone and cannot reflect the complete deformation pattern of this earthquake. Furthermore, since the DInSAR LOS displacements can only capture the one-dimensional projection of the real three-dimensional (3D) dis-placements, ascending/descending DInSAR LOS displacements cannot directly reveal the real deformation characteristics and cannot adequately constrain the fault slip model.
After several days of the earthquake, multiple SAR data with complete coverage (e.g., advanced land observation satellite-2, ALOS-2, and Sentinel-1) were available; therefore, it is feasible to obtain different SAR observations with various imaging geometries using various data processing methods (e.g., DInSAR [12], pixel offset-tracking, POT [13], multiple aperture InSAR, MAI [14], and burst overlap InSAR, BOI [15]), as well as to calculate the 3D displacements of this earthquake. For calculating 3D displacements, the traditional weighted least square (WLS) method (e.g., [16,17]) calculates 3D displacements of a target point using only the SAR observations at this point, and it determines the variances/weights of SAR observations on the basis of some a priori information, e.g., the displacements within a limited range are constant. However, the real deformations at adjacent points are generally correlated under the action of continuous stress, and the a priori information cannot adequately reflect the real complex noise distributions of SAR observations. Therefore, the WLS-obtained 3D displacements generally have a lower accuracy. Compared with the WLS method, a method for measuring 3D surface displacements with InSAR based on a strain model (SM) and variance component estimation (VCE) (SM-VCE method) [18], which was developed by our team, has proven to be superior in estimating 3D displacements [19][20][21]. The superiority lies in the fact that (1) the SM is applied to describe the correlation of 3D displacements between adjacent points and to establish the relationship between the 3D displacements at a target point and SAR observations at surrounding points, and (2) the well-known VCE algorithm is used to accurately determine the weights of SAR observations in an a posteriori manner.
In this study, we focus on estimating a complete and accurate 3D co-seismic displacement field related to the 2022 Menyuan earthquake on the basis of the ascending/descending Sentinel-1 SAR data and descending ALOS-2 SAR data. The DInSAR, POT, MAI, and BOI [15] methods are all used to measure the ground surface displacements along the line-of-sight (LOS) or azimuth directions. The SM-VCE method was applied here to combine the DInSAR-, MAI-, POT-, and BOI-derived displacement observations to obtain 3D co-seismic displacements of the 2022 Menyuan earthquake. Since more SAR displacement observations are used and an advanced method (i.e., SM-VCE) is applied to calculate 3D displacements, the obtained 3D displacements can accurately reveal the co-seismic deformation of this earthquake and provide a reliable dataset for further model interpretation. Moreover, the co-seismic strains were also calculated on the basis of the 3D displacements. In particularly, to benefit more researchers interested in co-seismic surface deformation assessment and the SM-VCE method, the obtained displacement data in geotiff format and the used SM-VCE code can be freely downloaded (see Data Availability Statement). Table 1 presents the basic information of the SAR data used in this paper. Their coverages can be found in Figure 1. All of these SAR data were acquired at several days after the earthquake; therefore, it is inevitable to include the post-seismic deformations in the SARderived displacement observations, which is normal for SAR earthquake studies [22,23]. We also investigated the post-seismic displacement in the descending Sentinel-1 pair 20220110-20220122, and we found that the post-seismic displacement of this earthquake was negligible compared with the co-seismic displacement; therefore, it was practical to combine these SAR-data-derived observations to obtain a 3D displacement field even though the SAR acquisition dates were inconsistent.

The Used SAR Data
The DInSAR, POT, MAI, and/or BOI methods were applied here to obtain the displacement observations from the aforementioned SAR data, and the details of the data processing procedure can be found in [21]. As for the ALOS-2 MAI displacement observation, we conducted a directional filtering and interpolation procedure to mitigate the possible ionospheric noise [24]. Figure 2 shows the SAR displacement observations. As can be seen, the spatial patterns of the SAR displacements with similar imaging geometry were consistent with each other, e.g., the LOS displacements from DInSAR and POT method (i.e., the first and second columns in Figure 2). Nevertheless, different method-derived displacements have their own superiorities. For example, the DInSAR LOS displacement has a much higher accuracy compared with the POT LOS displacement, but the POT method can measure larger displacement gradient compared with the DInSAR method. Although the POT/MAI/BOI methods have a lower displacement accuracy, they can obtain the displacements along the azimuth direction (near the north-south direction), which is significant since the DInSAR LOS displacement is almost blind to the north-south displacement component. Since the accuracy of the POT method depends highly on the spatial resolution of the used SAR images, the POT-derived azimuth displacements of Sentinel-1 data and the LOS displacements of ALOS-2 data were seriously contaminated by noise. Hence, these four SAR displacement observations (i.e., Figure 2c,g,k,n) were not used to calculate the 3D co-seismic displacements of this earthquake.

The SM-VCE Method
To obtain the 3D co-seismic displacements of the 2022 Menyuan earthquake, the SM-VCE method [18] was employed to combine the SAR displacement observations. Compared with the traditional WLS method (e.g., [16,17]) that calculates the 3D displacements for a target point using only the SAR observations at this point, the SM-VCE method incorporates more observations at the surrounding points for calculating the target 3D displacements. This increment relies on the fact that the strain model [25,26] represents the geophysical relationship between 3D displacements of adjacent points and can be used to link the SAR displacement observations at surrounding points with the target 3D displacements.
For example, when calculating the 3D displacements it is necessary to predefine a window around P 0 . According to the strain model, for a point P k with 3D position and displacement components being T , respectively, we can obtain the following relationship [26]: where , superscripts e, n, v denote east-west, northsouth, and vertical components, respectively, and H is the displacement gradient matrix, which can be written as

The SM-VCE Method
To obtain the 3D co-seismic displacements of the 2022 Menyuan earthquake, the SM-VCE method [18] was employed to combine the SAR displacement observations. Compared with the traditional WLS method (e.g., [16,17]) that calculates the 3D displacements Generally, H is assumed to be constant in the predefined window. Equation (1) can be rewritten as Assuming that there are J types of SAR displacement observations at P k , i.e., , it is easy to establish the relationship between the 3D displacement components d k and the observation vector L k . where is a unit vector that projects the 3D displacements at point P k to the corresponding j-th observations, and when the L k j is the displacement observation along the LOS or AZI direction, respectively. α k j and θ k j are the satellite heading angle (clockwise from the north) and the radar incidence angle at P k , respectively [18]. By combining Equations (3) and (6), we can establish the relationship between the unknown vector l and the SAR observations L k .
For K j surrounding points of P 0 for the j-th type of observation, the overall observation functions can be constructed on the basis of Equation (9).
where the size of observation vector L and the coefficient matrix B are ∑ J 1 K j × 1 and ∑ J 1 K j × 12, respectively.
Benefiting from the increasing SAR displacement observations used for calculating the target 3D displacements (i.e., Equation (12)), the famous VCE algorithm [27,28] is used to determine the accurate weights of these J types of SAR displacement observations in a posteriori way.
In the SM-VCE method, one critical parameter is the window size when establishing Equation (11) according to the strain model. Since the displacement gradients are assumed to be constant in the predefined window, too large a window size will violate this assumption, and too small a window size will significantly decrease the number of SAR displacement observations (i.e., K j ), as well as the performance of the VCE algorithm.
In this paper, we empirically determined a window size of 40 pixels × 40 pixels, which corresponds to about 2 km × 2 km with the observations' spatial resolution of about 50 m × 50 m. In such a spatial scale of this case study, the displacement gradients can be regarded as constants, and the spatially high-frequency noises in POT observations can also be well suppressed [21]. Figure 3 shows the 3D co-seismic displacements of the Menyuan earthquake derived using the SM-VCE method as a function of the SAR displacement observations in Figure 2. It is easy to infer that this earthquake was dominated by left-lateral strike slip, which is consistent with the overall tectonic movement of the QHFS. The maximum vertical displacement occurred near the ruptured faults with the magnitude of 0.6 m. Since the direction of the ruptured faults was nearly east-west, the east-west displacement component showed much larger magnitude compared with the north-south and vertical components. Figure 3e is the projection of the horizontal displacements to the fault orientation (about 114 • clockwise from the north), where the negative and positive values represent northwestward and southeastward horizontal displacements, respectively. The spatial patterns of horizontal displacements on both sides of the fault are generally symmetrical. Figure 3f presents the across-fault displacements of six selected profiles along the fault, illuminating that the maximum horizontal displacement could reach about 1.9 m, and the relative horizontal displacement across both sides of the fault was as large as 2-3 m, which is consistent with the field investigations [11]. In addition, Figure 3f indicates that the absolute displacement magnitude was larger on the southern side of the fault compared with the displacement magnitude on the northern side. This may be attributed to the fact that, across the QHFS, the Qinghai-Tibet Plateau block on the southern side of this fault system is relatively active and continuously pushing northward, but being resisted by the Gebi-Alashan block and the Huabeikelatong block, which are relatively stable compared with the Qinghai-Tibet Plateau block [4,6,29]. In this case, this earthquake is of great significance for studying the tectonic stress across the QHFS. Figure 4 presents the corresponding variance of the 3D co-seismic displacements estimated using the SM-VCE method. The variances of the east-west and vertical components were generally at a centimeter to millimeter level, while the north-south component was at a decimeter level, which is reasonable since the north-south component was mainly contributed by the low-accuracy SAR azimuth displacements, while the relatively highaccuracy SAR displacements mainly contributed to the east-west and vertical components. In the north-south variance map, there were obvious striped zones with relatively lower values since, at these zones, the BOI observations were also used to calculate 3D displacements. Since the correlation between the north-south and vertical components was high when calculating 3D displacements with the current SAR observations [30], the striped zones also appeared in the vertical variance map. Moreover, a far-field zone (the red dashed rectangles in Figure 3a-c) was selected to assess the accuracy of the obtained 3D displacements. In this zone, we assumed that no co-seismic displacement occurred; hence, the obtained displacement signals in this zone were regarded as noise and could be used to estimate the variance of the corresponding component. The results show that the variance values were 8.6 mm, 27.5 mm, and 5.9 mm for the east-west, north-south, and vertical components, respectively, which are generally consistent with the results in Figure 4. Error! Reference source not found. presents the corresponding variance of the 3D co-seismic displacements estimated using the SM-VCE method. The variances of the eastwest and vertical components were generally at a centimeter to millimeter level, while the north-south component was at a decimeter level, which is reasonable since the northsouth component was mainly contributed by the low-accuracy SAR azimuth displacements, while the relatively high-accuracy SAR displacements mainly contributed to the east-west and vertical components. In the north-south variance map, there were obvious striped zones with relatively lower values since, at these zones, the BOI observations were also used to calculate 3D displacements. Since the correlation between the north-south and vertical components was high when calculating 3D displacements with the current SAR observations [30], the striped zones also appeared in the vertical variance map. Moreover, a far-field zone (the red dashed rectangles in Error! Reference source not found.ac) was selected to assess the accuracy of the obtained 3D displacements. In this zone, we assumed that no co-seismic displacement occurred; hence, the obtained displacement signals in this zone were regarded as noise and could be used to estimate the variance of the corresponding component. The results show that the variance values were 8.6 mm, 27.5 mm, and 5.9 mm for the east-west, north-south, and vertical components, respectively, which are generally consistent with the results in Error! Reference source not found..

Three Strain Invariants of the 2022 Menyuan Earthquake
The three strain invariants [21] comprise dilatation, differential rotation, and maximum shear, which are calculated from the 3D displacement gradients, and are independent of any shift and rotation of the coordinate systems [25]. Compared with the displacement fields, the strain fields can provide some unique view about the surface deformations, and they have been widely used for analyzing the inter-seismic deformations on

Three Strain Invariants of the 2022 Menyuan Earthquake
The three strain invariants [21] comprise dilatation, differential rotation, and maximum shear, which are calculated from the 3D displacement gradients, and are independent of any shift and rotation of the coordinate systems [25]. Compared with the displacement fields, the strain fields can provide some unique view about the surface deformations, and they have been widely used for analyzing the inter-seismic deformations on the basis of sparse GNSS (global navigation satellite system) data [31][32][33]. For example, with respect to the inter-seismic deformations, zones with higher magnitude of the strain invariants indicate a greater potential for stress accumulation and seismic risk [33]. For the high-spatialresolution co-seismic displacements, it is also practical to calculate the strain invariants and to assist the analysis of co-seismic behavior [34][35][36][37] (details about the calculation can be found in [21]). One point that should be noted is that, for the inter-seismic deformations, the deformations are generally continuously occurring; thus, the strain analysis can benefit the assessment of seismic risk. However, for the co-seismic events, since the displacements have occurred, the strain analysis can be used for the disaster assessment and even some post-event analysis.
As shown in Figure 5, the strain invariants were more obvious near the ruptured fault zones compared with the far-field area. This is reasonable since the strains were calculated on the basis of displacement gradients, which were of high magnitude near the fault compared with the far-field area. For the dilatation invariant, the positive and negative values indicate that the volume around this point was expanding and compressed, respectively. As shown in Figure 5a, there were two obvious positive (expanding) zones and two obvious negative (compressed) zones near the fault. By comparing with the horizontal displacements (Figure 3d), it can be found that the expanding or compressed zones corresponded to the areas with increasing or decreasing magnitude of horizontal displacements, respectively. Both the expanding and the compressed areas can be considered to be deformed due to the external force. This force may not only be related to the inter-seismic deformation but also affect the post-event behavior of the fault, and further investigations can be conducted with respect to the dilatation invariant. For the differential rotation invariant, negative or positive values represent the clockwise or anticlockwise rotation of the horizontal displacements around a target point. It is easy to infer that both sides of the fault experienced a similar clockwise rotation, which is consistent with the horizontal displacement pattern, whereby the horizontal displacement direction changed gradually to the right of the current direction (Figure 3d). For the maximum shear invariant, the value was positive, generated in the area where relative horizontal displacements occurred. For example, given two points with a similar horizontal direction, but with different magnitudes, there will be maximum shear strain. As can be seen, the maximum shear strain signal only concentrated near the fault area, indicating that significant relative deformations occurred in these zones.

The Contributions of Different SAR Observations to the 3D Displacements
It is necessary to combine SAR displacement observations from at least three imaging geometries with significant difference for calculating the real 3D displacements. For the current SAR satellites, two imaging geometries can be easily obtained, i.e., ascending and descending DInSAR LOS directions with the right-looking mode. The LOS observations are sensitive to the east-west and vertical displacement components, but almost blind to the north-south component. Therefore, by combing the easily acquired ascending/descending DInSAR LOS observations, many studies and applications have been conducted to only estimate the east-west and vertical displacements by ignoring the north-south component [30,[38][39][40][41][42]. In order to obtain the real 3D displacements, one must find at least one another SAR displacement observation that can complement the ascending/descending DInSAR LOS observations. Considering the current SAR data and the data processing methods, there are mainly two ways to realize the calculation of 3D displacements. One way is to incorporate the left-looking mode SAR data, whose imaging geometry is quite different from the right-looking mode SAR data. By combining ascending/descending left-/rightlooking mode SAR data, it is practical to calculate 3D displacements with high accuracy using only the DInSAR method [19,43]. The main shortcoming lies in that left-looking SAR data are rarely available. Another way is to incorporate the azimuth SAR displacement observations. In addition to the DInSAR method, the POT, MAI, and BOI methods can derive the displacement observations along the satellite azimuth direction, which is near the north-south direction. In this case, by combining the ascending/descending SAR data with the DInSAR and POT/MAI/BOI methods, it is also feasible to calculate reliable 3D displacements [44][45][46][47][48][49][50]. Since the accuracy of POT/MAI/BOI methods is not as high as that of the DInSAR method, the measurable north-south displacement magnitude is generally higher than a dozen centimeters. gradually to the right of the current direction (Error! Reference source not found.d). For the maximum shear invariant, the value was positive, generated in the area where relative horizontal displacements occurred. For example, given two points with a similar horizontal direction, but with different magnitudes, there will be maximum shear strain. As can be seen, the maximum shear strain signal only concentrated near the fault area, indicating that significant relative deformations occurred in these zones.

The Contributions of Different SAR Observations to the 3D Displacements
It is necessary to combine SAR displacement observations from at least three imaging geometries with significant difference for calculating the real 3D displacements. For the current SAR satellites, two imaging geometries can be easily obtained, i.e., ascending and descending DInSAR LOS directions with the right-looking mode. The LOS observations are sensitive to the east-west and vertical displacement components, but almost blind to the north-south component. Therefore, by combing the easily acquired ascending/descending DInSAR LOS observations, many studies and applications have been conducted to only estimate the east-west and vertical displacements by ignoring the north-south component [30,[38][39][40][41][42]. In order to obtain the real 3D displacements, one must find at least one another SAR displacement observation that can complement the ascending/descending DInSAR LOS observations. Considering the current SAR data and the data processing methods, there are mainly two ways to realize the calculation of 3D displacements. One way is to incorporate the left-looking mode SAR data, whose imaging geometry is quite different from the right-looking mode SAR data. By combining ascending/descending left-/right-looking mode SAR data, it is practical to calculate 3D displacements with high accuracy using only the DInSAR method [19,43]. The main shortcoming lies in that left-looking SAR data are rarely available. Another way is to incorporate the azimuth SAR displacement observations. In addition to the DInSAR method, the POT, MAI, and BOI methods can derive the displacement observations along the satellite azimuth direction, which is near the north-south direction. In this case, by combining the ascending/descending SAR data with the DInSAR and POT/MAI/BOI methods, it is also feasible to calculate reliable 3D displacements [44][45][46][47][48][49][50]. Since the accuracy of POT/MAI/BOI methods is not as high as that of the DInSAR method, the measurable north-south displacement magnitude is generally higher than a dozen centimeters.
For the 2022 Menyuan earthquake, it was easy to obtain the SAR displacement observations along two geometries, i.e., ascending and descending LOS direction, using only on the free Sentinel-1 SAR data (e.g., Error! Reference source not found.e,i). Due to the For the 2022 Menyuan earthquake, it was easy to obtain the SAR displacement observations along two geometries, i.e., ascending and descending LOS direction, using only on the free Sentinel-1 SAR data (e.g., Figure 2e,i). Due to the low azimuth pixel's spatial resolution of Sentinel-1 SAR images, the signal-to-noise ratio (SNR) of the POT-derived azimuth displacement observations (e.g., Figure 2g,k) was too low to identify any valuable displacement signal. Although the displacement accuracy was relatively high for the BOI observations compared with the POT azimuth observations, the BOI coverage was too limited to reveal the detailed spatial variety of surface displacements. In this case, if considering only the free Sentinel-1 SAR data, one can hardly obtain complete and accurate 3D displacements associated with this event. Unlike the Sentinel-1 SAR images, the azimuth pixel's spatial resolution of the ALOS-2 SAR images is as high as 2-3 m; therefore, the POT method is sufficient to measure the azimuth displacement with considerable accuracy. Simultaneously, the MAI method can also be applied to the ALOS-2 SAR data to derive the azimuth displacement observation, which shares the same imaging geometry with the POT azimuth displacement observation. In this sense, the ALOS-2 azimuth displacement observations (e.g., Figure 2o,p) play a very important role in calculating the 3D displacements of this event, especially by contributing to the north-south displacement component, since other SAR displacement observations in Figure 2 cannot adequately reveal the displacement pattern in the satellite azimuth geometry.

Calculating the Displacements/Strains near the Ruptured Fault Zone
When calculating the displacements or strains of a target point, it is inevitable to incorporate the surrounding points. For example, when calculating the 3D displacements using the SM-VCE method, one should establish the observation functions (i.e., Equation (11)) on the basis of the observations within a predefined window. When calculating the strain invariants, it is also necessary to use the 3D displacements at surrounding points. One should pay attention that the precondition for reliable calculation of the displacements/strains lies such that the displacement gradients among these surrounding points are constant. However, when the target point is located near the ruptured fault and these surrounding points are located on both sides of the fault, it is obvious that the displacement gradients are different among these points. In this case, if the surrounding points on the other side of the fault are not distinguished and not discarded in the calculation, the obtained displacements/strains will be contaminated and cannot appropriately reflect the real deformation pattern.
For the 2022 Menyuan earthquake, the fault ruptured the ground surface, and there were displacement jumps across the ruptured fault. In order to accurate calculate the displacements/strains near the ruptured fault zone, we first manually mapped the fault traces on the basis of POT observations and excluded those points on the opposite side of the fault in the calculation. One can image that if those points on the opposite side of the fault were not excluded, the obtained 3D displacements would be very smooth across the fault. Moreover, since the strain invariants were calculated on the basis of the displacement gradients, the magnitude of the strain invariants would be very large near the displacement jump zone compared with the far-field zone. These possible results are wrong and do not correspond with the real deformation of earthquakes.

Conclusions
This paper calculated the 3D surface displacements of the 2022 Menyuan earthquake using ascending/descending Sentinel-1 SAR data and descending ALOS-2 SAR data. To derive SAR displacement observations along different directions (i.e., the satellite line-ofsight and azimuth directions), DInSAR, POT, MAI, and/or BOI methods were applied to process these SAR data. With such abundant SAR displacement observations from various observing directions, it is clear that the co-seismic deformation of this earthquake could be well constrained. Afterward, an SM-VCE method was used to combine these SAR displacement observations to obtain a high-accuracy 3D displacement field. Since the SM was used to establish the relationship between SAR observations and 3D displacements, and since the VCE was used to determine the accurate weights of SAR observations in the SM-VCE method, the obtained 3D co-seismic displacement was more accurate and reliable compared with the traditional WLS method, and the precision was generally at the level of centimeters to millimeters in this study. Results indicate that the 2022 Menyuan earthquake was dominated by left-lateral slip, and the maximum horizontal and vertical displacements were~1.9 m and 0.6 m, respectively. The relative horizontal surface displacement across the fault was as large as 2-3 m, and the fault-parallel displacement magnitude was larger on the southern side of the fault compared with the northern side. Furthermore, three strain invariants were also calculated and revealed that the near-fault zones suffered the most severe strains, and two expanding and compressed zones were also identified.
To calculate reliable 3D displacements near the fault zone, the manually mapped fault traces from the SAR POT observations were used to assist the SM-VCE method. However, the manually mapped fault traces are not always available; thus, it is necessary to apply some automatic algorithms to increase the robustness of SM-VCE method near the displacement jump zone. Additionally, since the SM-VCE method uses more observations to calculate 3D displacements, the computational efficiency still has much room for improvement. Given that this study aimed to provide a quick and reliable 3D displacement field of the 2022 Menyuan earthquake, we did not investigate the fault model of this event. Although some fault models have been published on the basis of DInSAR LOS observations (e.g., [2,11]) or seismic wave data (e.g., [1,51]), it is worthy to conduct a slip model inversion using the obtained 3D displacements in the future.