DEM Generation from GF-7 Satellite Stereo Imagery Assisted by Space-Borne LiDAR and Its Application to Active Tectonics

: China’s ﬁrst optical stereo mapping satellite with a sub-meter resolution, GaoFen-7 (GF-7), launched in November 2019, shows signiﬁcant potential for providing high-resolution topographic and geomorphic data for quantitative research on active tectonics


Introduction
With the rapid development of spatial measurement technology and dating methods, active tectonic research has gradually developed from the early qualitative to the quantitative research stage [1,2]. Quantitative research on active tectonics should acquire some quantitative parameters that represent the active tectonic characteristics, for which the acquisition depends on topographic and geomorphic data with high precision and a high resolution [3][4][5][6].
Various methods have been presented to obtain high-precision digital elevation models (DEMs), including traditional topographic field surveys, synthetic aperture interferometry radar (InSAR), light detection and ranging (LiDAR), unmanned aerial vehicle (UAV) surveys, and high-resolution satellite stereo mapping. Since traditional measurement methods are time-consuming and labor-intensive, it is challenging to obtain large-scale, high-precision topographic data along the entire fault [3]. Although airborne LiDAR can obtain sub-meter high-precision topographic data, it cannot be applied to large areas due to its high acquisition cost [2,7,8]. Recently, the combination of UAV technology and the structure from motion (SfM) photogrammetry method has provided high portability and flexibility, significantly reducing the cost of acquiring high-resolution topographic and geomorphic data [9]. However, applying this method to large areas with a harsh climate and complex terrain is challenging. Although InSAR can obtain topographic data for a large area, the accuracy of its data cannot be guaranteed because the data processing technique is not yet mature [10,11]. The emergence of DEM extraction based on high-resolution satellite stereo imagery provides a new solution for the rapid acquisition of large-scale topographic and geomorphic data in fault zones [6,12,13], significantly reducing the cost and improving the efficiency of acquiring topographic and geomorphic data.
Many previous studies have used commercial high-resolution stereo images such as QuickBird, Pleiades 1A/1B, and WorldView-1/2/3/4 to perform quantitative research on active tectonics [14][15][16][17]. However, these commercial images are expensive and have poor timeliness, restricting their application in the field of active tectonics. China's first optical stereo mapping satellite with a sub-meter resolution, Gaofen-7 (GF-7), was successfully launched in 2019. The satellite is equipped with a dual-line-array stereo camera comprising a front-view camera and a rear-view camera. It can effectively acquire panchromatic stereo images with a width of 20 km and a resolution of 0.8 m, and multi-spectral images with a resolution of 3.2 m [18]. Although GF-7 exhibits significant potential for providing largescale high-resolution topographic and geomorphic data, research has yet to explore the potential of applying the GF-7 data to the quantitative study of active tectonics.
Ground control points (GCPs) are necessary to realize accurate DEM generation based on GF-7 stereo images [19,20]. However, acquiring GCPs in the field is often timeconsuming and laborious, especially in deserts, woodlands, and areas with complex terrain. Many previous studies have utilized public geographic information data such as Google Earth images and SRTM DEM to improve the accuracy of block adjustments, thereby improving the accuracy of DEMs generated from stereo images [21][22][23]. Recently, with the rapid development of space-borne LiDAR, reliable GCPs have been provided for the block adjustment of satellite stereo images. Existing studies have extracted GCPs using the multi-beam micro-pulse photon-counting LiDAR Ice, Cloud, and land Elevation Satellite-2 (ICESat-2). Although it can effectively improve the accuracy of DEM extraction based on stereo images [24,25], no research has explored the effect of DEM extraction with or without GCPs on the measurement accuracy of active tectonic parameters.
This study aimed to validate the accuracy of the GF-7-extracted DEMs with or without ICESat-2-derived GCPs, and evaluate the potential of applying GF-7 DEMs to active tectonics. In particular, there were three specific objectives: (1) to verify the accuracy of GF-7 DEMs with or without GCPs; (2) to assess the ability of GF-7 DEMs to identify active faults, fault scarps, and horizontal offsets; and (3) to evaluate the potential application of GF-7 DEMs to horizontal and vertical offset measurements with or without GCPs. Ultimately, this study attempted to answer the following questions: (1) Is there any difference between the accuracies of DEMs extracted from GF-7 with and without ICESat-2 derived GCPs? (2) What horizontal and vertical offset levels can be identified on the basis of GF-7 DEMs? (3) Are there any differences in the measurement accuracies of horizontal and vertical offsets with or without the GCPs extracted from ICESat-2 data?

Study Site
The Haiyuan fault is one of the most critical faults in the northeastern margin of the Tibetan Plateau. It connects with the eastern Liupanshan pediment fault to the east and the Maomaoshan fault to the west (Figure 1a). It is one of the fault zones with the most significant tectonic deformation in the northeastern margin of the Tibetan Plateau [26]. Many large earthquakes have been documented in the history of the Haiyuan fault, such as the M8.5 Haiyuan earthquake in 1920 and the M8-8.3 Gulang earthquake in 1927. The Haiyuan fault is dominated by a left-lateral strike-slip with a small thrust component. Because of their long-term activity, numerous typical strike-slip dislocation landforms have developed along the Haiyuan fault, forming various offset levels. Moreover, the dry climatic environment of the Haiyuan fault makes the landforms well-preserved. In summary, the Haiyuan fault is ideal for quantitative research on active faults [4].
tical offsets with or without the GCPs extracted from ICESat-2 data?

Study Site
The Haiyuan fault is one of the most critical faults in the northeastern margin of the Tibetan Plateau. It connects with the eastern Liupanshan pediment fault to the east and the Maomaoshan fault to the west (Figure 1a). It is one of the fault zones with the most significant tectonic deformation in the northeastern margin of the Tibetan Plateau [26]. Many large earthquakes have been documented in the history of the Haiyuan fault, such as the M8.5 Haiyuan earthquake in 1920 and the M8-8.3 Gulang earthquake in 1927. The Haiyuan fault is dominated by a left-lateral strike-slip with a small thrust component. Because of their long-term activity, numerous typical strike-slip dislocation landforms have developed along the Haiyuan fault, forming various offset levels. Moreover, the dry climatic environment of the Haiyuan fault makes the landforms well-preserved. In summary, the Haiyuan fault is ideal for quantitative research on active faults [4].
The study area is located in the Hasishan section of the Haiyuan fault (Figure 1b), including various fault scarps, water system dislocations, and trough landforms. The Hasishan section has obvious topographic relief and less vegetation coverage, which is helpful for evaluating the application potential of DEMs extracted from GF-7 data on active tectonics.  The study area is located in the Hasishan section of the Haiyuan fault (Figure 1b), including various fault scarps, water system dislocations, and trough landforms. The Hasishan section has obvious topographic relief and less vegetation coverage, which is helpful for evaluating the application potential of DEMs extracted from GF-7 data on active tectonics.

GF-7 Data
The GF-7 is China's first civil optical stereo mapping satellite with a sub-meter resolution. It was successfully launched in November 2019 with an orbital altitude of about 506 km and a return cycle period of fewer than 60 days. GF-7 aims to generate 1:10,000 stereo mapping products and update large-scale basic geographic information products [27,28]. It can effectively acquire 20 km wide front-view panchromatic stereoscopic images with a 0.8 m resolution, rear-view panchromatic stereoscopic images with a 0.65 m resolution, and multi-spectral images with a 2.6 m resolution. In this study, four GF-7 stereo images from December 2019 to December 2020 with less cloud and snow, a clear texture, and geometric structures were provided by the High-resolution Remote Sensing Data Center, China Earthquake Administration. Table 1 lists the specific information of these four GF-7 stereo images.

ICESat-2/ATLAS Data
ICESat-2 data were utilized to obtain GCPs to improve the accuracy of block adjustments of the GF-7 stereo images. ICESat-2 was launched in 2018 and is equipped with a multi-beam, micro-pulse photon-counting LiDAR system named the Advanced Topographic Laser Altimeter System (ATLAS), which can generate overlapping footprints with a diameter of 12 m spaced at 0.7 m intervals, covering an areas from 88 • S to 88 • N [29]. ATLAS emits three pairs of beams, each containing strong and weak beams. The cross-track distances between and within the pairs are about 3.3 km and 90 m [30].
ICESat-2 data products are divided into Level 1, Level 2, Level 3A, and Level 3B. The Level-3 A product, ATL08, records the elevation of every 100 m segment. The latest ATL08 data (Version 5) from October 2018 to May 2022 were downloaded through the website https://search.earthdata.nasa.gov/ (accessed on 2 October 2022). In order to select high-precision GCPs, 15 parameters were extracted from the ATL08 data, including the longitude, latitude, elevation, and other parameters ( Table 2).

Airborne LiDAR Data
The airborne LiDAR data used in this study were collected in 2011 with a Leica ALS-60 system [3,4]. The swath width of the airborne LiDAR data was 1 km, with an average point density of 4-6 points/m 2 . The triangulated irregular network (TIN) interpolation algorithm was utilized to generate DEM products with a spatial resolution of 1 m based on the cloud points. The coordinate system of the produced DEM was the WGS84 universal transverse Mercator (UTM) projection system. The accuracy of the airborne LiDAR DEM's horizontal and vertical elevation was 0.3 m and 0.2 m, respectively [3,31]. On the basis of the airborne LiDAR-derived DEM, 20 checkpoints with apparent marks were selected to evaluate the horizontal and vertical accuracies of the DEMs extracted from the GF-7 stereo images ( Figure 2a). In addition, 8 positions were selected for measurements of the vertical offset and 13 for measurements of the horizontal offset ( Figure 2b). subset_te_flag: quality flag terrain_slope: the along-track terrain slope of each 100 m segment

Airborne LiDAR Data
The airborne LiDAR data used in this study were collected in 2011 with a Leica ALS-60 system [3,4]. The swath width of the airborne LiDAR data was 1 km, with an average point density of 4-6 points/m 2 . The triangulated irregular network (TIN) interpolation algorithm was utilized to generate DEM products with a spatial resolution of 1 m based on the cloud points. The coordinate system of the produced DEM was the WGS84 universal transverse Mercator (UTM) projection system. The accuracy of the airborne LiDAR DEM's horizontal and vertical elevation was 0.3 m and 0.2 m, respectively [3,31]. On the basis of the airborne LiDAR-derived DEM, 20 checkpoints with apparent marks were selected to evaluate the horizontal and vertical accuracies of the DEMs extracted from the GF-7 stereo images ( Figure 2a). In addition, 8 positions were selected for measurements of the vertical offset and 13 for measurements of the horizontal offset ( Figure 2b).

Ancillary Data
Google Earth images of the study area with a 1 m resolution were downloaded using Shuijingzhu software, and SRTM DEM data with a 30 m resolution were downloaded from the Google Earth Engine (https://code.earthengine.google.com/ (accessed on 11 September 2022)). The geographic coordinate system of the Google Earth images and SRTM DEM was the World Geodetic System 1984 (WGS84). The elevation datum of the SRTM DEM was the WGS84 ellipsoid.

Extraction of GCPs from ICESat-2/ATLAS Data
Previous studies have demonstrated that ICESat-2/ATLAS data can be used to extract GCPs. In this study, we extracted a total of 58,419 ATL08 points in the study area as the candidates of the GCPs. However, because of the influence of the atmospheric environment, the targets' characteristics, the hardware, and other factors, only some ATL08 points had high-accuracy elevation data. Thus, only the reliable ATL08 points were selected as the GCPs. This study utilized the GCP extraction strategy proposed by [25] to obtain highprecision GCPs from the ATL08 data. Previous studies have indicated that the strong beams had better elevation accuracy than weak beams, and the beams collected at nighttime also had better elevation accuracy than those collected during the daytime [32]. Additionally, considering that some ATL08 points are densely distributed, the densely

Ancillary Data
Google Earth images of the study area with a 1 m resolution were downloaded using Shuijingzhu software, and SRTM DEM data with a 30 m resolution were downloaded from the Google Earth Engine (https://code.earthengine.google.com/ (accessed on 11 September 2022)). The geographic coordinate system of the Google Earth images and SRTM DEM was the World Geodetic System 1984 (WGS84). The elevation datum of the SRTM DEM was the WGS84 ellipsoid.

Extraction of GCPs from ICESat-2/ATLAS Data
Previous studies have demonstrated that ICESat-2/ATLAS data can be used to extract GCPs. In this study, we extracted a total of 58,419 ATL08 points in the study area as the candidates of the GCPs. However, because of the influence of the atmospheric environment, the targets' characteristics, the hardware, and other factors, only some ATL08 points had high-accuracy elevation data. Thus, only the reliable ATL08 points were selected as the GCPs. This study utilized the GCP extraction strategy proposed by [25] to obtain highprecision GCPs from the ATL08 data. Previous studies have indicated that the strong beams had better elevation accuracy than weak beams, and the beams collected at nighttime also had better elevation accuracy than those collected during the daytime [32]. Additionally, considering that some ATL08 points are densely distributed, the densely distributed ATL08 points were further filtered to ensure the uniform spatial distribution of the GCPs. Therefore, three more filter criteria were added, in line with the method proposed in [25]. Table 3 shows all the filter criteria for the extraction of GCPs. This study only selected the ATL08 points that met all the filtering criteria as the GCPs. Five flags of subset_te_flag greater than −1 with the middle three flags equal to 1 10 cloud_flag_atm <= 2 11 segment_landcover = 60 12 The distances between ATL08 points should be larger than 500 m Figure 3 shows the number of points remaining after each data filtering step. Finally, 114 ATL08 points were selected as the GCPs by the 12 filter criteria ( Figure 4).
Remote Sens. 2023, 15, x FOR PEER REVIEW 6 of 20 distributed ATL08 points were further filtered to ensure the uniform spatial distribution of the GCPs. Therefore, three more filter criteria were added, in line with the method proposed in [25]. Table 3 shows all the filter criteria for the extraction of GCPs. This study only selected the ATL08 points that met all the filtering criteria as the GCPs. Five flags of subset_te_flag greater than −1 with the middle three flags equal to 1 10 cloud_flag_atm <= 2 11 segment_landcover = 60 12 The distances between ATL08 points should be larger than 500 m Figure 3 shows the number of points remaining after each data filtering step. Finally,114 ATL08 points were selected as the GCPs by the 12 filter criteria ( Figure 4).

Extraction of DEMs from GF-7 Stereo Images
For the GF-7 stereo images, the transformation relationship between the images' coordinates and the object space coordinates (Equation (1)) was expressed by a rational function model (RFM), which contained multiple rational polynomial coefficients (RPCs). These RPCs are often physically meaningless and inaccurate [33]. Therefore, the block ad-

Extraction of DEMs from GF-7 Stereo Images
For the GF-7 stereo images, the transformation relationship between the images' coordinates and the object space coordinates (Equation (1)) was expressed by a rational function model (RFM), which contained multiple rational polynomial coefficients (RPCs). These RPCs are often physically meaningless and inaccurate [33]. Therefore, the block adjustment method was used to add the affine transformation parameters of the image for compensation (Equation (2)) and improve the geometric imaging model of the stereo images and the DEM product's accuracy [27,34,35].
where (c, r) are the image coordinates; (X, Y, Z) is the object space coordinates; p 1 , p 2 , p 3 , and p 4 are the RPCs; (∆c, ∆r) are the compensation values of the image coordinates; and a i , b i (i = 0, 1, 2) are the affine transformation coefficients. Currently, there are three commonly used block adjustment methods: • Block adjustment without GCPs. This method was essential for selecting some tie points in the GF-7 images and calculating the object space coordinates for each pair of GF-7 stereo pixels, utilizing the space-forward intersection. The final object space coordinate of the tie point was the average coordinate of the tie points between each pair of GF-7 stereo pixels. These tie points were utilized as the virtual GCPs in the block adjustment process. • Block adjustment with the aid of geographic information system (GIS) data. This method utilized the existing GIS data (such as digital orthophoto maps (DOM) and DEMs) to assist the block adjustment of stereo images. First, a large number of cognominal points were obtained through the automatic registration of the GF-7 satellite images and the existing DOM. Second, the DOM and DEM were used to obtain these cognominal points' horizontal and elevation coordinates, respectively. Finally, these cognominal points were taken as the control points. In addition to the basic geographic information data (the DOM and DEM), public geographic information data can also be utilized. The most commonly used public geographic information data are Google Earth images and SRTM DEM [21][22][23]. • Block adjustment with GCPs. This method used high-precision GCPs to constrain the elevation values of the forward intersection of the GF-7 stereo images. The most commonly used GCPs can be obtained from space-borne LiDAR data, such as ICESat-1/GLAS and ICESat-2/ATLAS.
Space Data Processor (SDP), a large-scale remote sensing image processing software package, was utilized to perform the block adjustment of the GF-7 satellite stereo images. SDP is a comprehensive, highly automated system for processing remote sensing images independently developed by the Land Satellite Remote Sensing Application Center of the Ministry of Natural Resources of the People's Republic of China. It can realize the three abovementioned block adjustment methods of GF-7 stereo images. The DEMs obtained by block adjustment without the GCPs, with the aid of 1 m Google Earth maps and 30 m SRTM DEM, and with the GCPs extracted from the ICESat-2 data were named GF-7 DEM Method1 , GF-7 DEM Method2 , and GF-7 DEM Method3 , respectively. The coordinate system of these generated GF-7 DEMs was WGS 84.

Measurement of the Horizontal and Vertical Offsets
LaDiCaoz_v2.1 software was adopted to measure the horizontal offsets automatically [14,36]. First, the contour lines were generated based on the DEM to identify the similar geomorphic features of the hanging wall and footwall of the fault. Second, the locations of the fault were identified according to the topographical information and the dislocated geological bodies of the landmark on both sides of the fault were selected for measurement of the horizontal offset (Figure 5a). Third, the characteristics of similar geological bodies on both sides of the fault were compared, and the offsets were automatically matched through horizontal movement, vertical movement, and stretching adjustment (Figure 5c). Finally, the geological body was automatically matched to obtain the optimal horizontal offsets and the dislocation recovery results of the fault (Figure 5b). In order to extract the height of the fault scarp, global mapper software was utilized to manually extract the survey lines perpendicular to the fault and obtain information on the relative distance and elevation. The vertical offsets and errors were calculated on the basis of the fitting lines of the hanging wall and the footwall of fault according to Equations (3) and (4) [37]. However, undulations in the terrain may cause differences in the elevation between the fitting lines and the DEM. Therefore, the fitting errors of hanging In order to extract the height of the fault scarp, global mapper software was utilized to manually extract the survey lines perpendicular to the fault and obtain information on the relative distance and elevation. The vertical offsets and errors were calculated on the basis of the fitting lines of the hanging wall and the footwall of fault according to Equations (3) and (4) [37]. However, undulations in the terrain may cause differences in the elevation between the fitting lines and the DEM. Therefore, the fitting errors of hanging wall and footwall lines were considered ( Figure 6). h h h 1 where h1 and h2 are the fault scarp heights extracted from the fitting line of the hanging wall (L1) and the footwall (L2), respectively; zi is the DEM's elevation value at the location xi;

Validation of Accuracy
Airborne LiDAR data were utilized to assess the accuracy of the GF-7 DEMs. Specifically, the 20 checkpoints extracted from the airborne LiDAR data were used to quantitatively evaluate the accuracies of the GF-7 DEMs extracted by the three different block adjustment methods. To further verify the accuracy of the elevation of the GF-7 DEMs, a comparison of the difference between the elevation of the GF-7 DEMs and of the airborne LiDAR data was conducted. As the GF-7 DEMs and airborne LiDAR data were acquired by different methods at different times, they were co-registered to make a reliable comparison. Considering that the accuracy of registration of the LiDAR and GF-7 data based on a point cloud may be higher than that of the interpolated DEMs, we compared the

Validation of Accuracy
Airborne LiDAR data were utilized to assess the accuracy of the GF-7 DEMs. Specifically, the 20 checkpoints extracted from the airborne LiDAR data were used to quantitatively evaluate the accuracies of the GF-7 DEMs extracted by the three different block adjustment methods. To further verify the accuracy of the elevation of the GF-7 DEMs, a comparison of the difference between the elevation of the GF-7 DEMs and of the airborne LiDAR data was conducted. As the GF-7 DEMs and airborne LiDAR data were acquired by different methods at different times, they were co-registered to make a reliable comparison. Considering that the accuracy of registration of the LiDAR and GF-7 data based on a point cloud may be higher than that of the interpolated DEMs, we compared the differences in the elevation of the LiDAR and GF-7 point clouds rather than the interpolated DEMs utilizing the software CloudCompare (http://www.danielgm.net/cc/ (accessed on 21 October 2019)). First, the GF-7 DEMs were projected into the UTM system and transformed to point clouds. Second, from the positioning error between GF-7 and airborne LiDAR point clouds, the two point clouds were aligned using the well-known iterative closest point (ICP) algorithm [38]. The GF-7 point clouds were precisely registered with the LiDAR point clouds by iteratively optimizing the rotation matrix and the translation vector. Finally, the differences in the point-to-point elevation between the GF-7 and airborne point clouds were calculated, which only represented the differences in the internal elevation of the two without being affected by their spatial differences [39]. We also evaluated the performance of the GF-7 DEMs by comparing the terrain slopes obtained by the different GF-7 DEMs with terrain slopes calculated on the basis of the airborne LiDAR DEM.
To quantitatively assess the accuracy of the measured horizontal and vertical offsets based on the GF-7 DEMs, the horizontal and vertical offsets measured by the airborne LiDAR-derived DEM were taken as the standard values, and three statistical indicators were included, namely the coefficient of determination (R 2 , Equation (5)), bias (Equation (6)), and the root mean square error (RMSE, Equation (7)).
where y i is the i th horizontal or vertical offset value extracted from the airborne LiDARderived DEM, x i is the i th horizontal or vertical offset value extracted from the GF-7 DEM, y is the average value of the horizontal or vertical offsets extracted from the airborne LiDAR-derived DEM, and n is the number of horizontal or vertical offsets. Figure 7 shows the GF-7 DEMs generated by the three different block adjustment methods. The accuracy of the GF-7 DEMs evaluated by 20 checkpoints in the X, Y, and Z directions is shown in Table 4. As shown in Figure 7, the three GF-7 DEMs were not visually significantly different. Table 4 shows that the RMSE values of GF-7 DEM Method1 were extremely large, namely 178.24 m, 198.90 m, and 110.58 m in the X, Y, and Z directions. The ∆X values gradually increased while the ∆Z values decreased with an increase in the number of checkpoints along the fault. The GF-7 DEM generated by [28] using the block adjustment method without GCPs yielded a plane positioning accuracy of 4.26 m and an elevation accuracy of 8 m. However, the plane positioning and elevation errors of GF-7 DEM Method1 exceeded 100 m. The reason may be that the GF-7 stereo images used in this paper were mainly collected during the initial operation phase. At this stage, the geometric calibration of GF-7 in orbit was not yet mature, leading to relatively large errors in the horizontal and vertical positioning. The obtained accuracies of GF-7 DEM Method2 and GF-7 DEM Method3 were superior to that obtained with DEM Method1 . The RMSE values of GF-7 DEM Method3 in the X, Y, and Z directions were 1.38 m, 1.73 m, and 1.35 m, respectively. This indicates that the use of GCPs extracted from ICESat-2/ATLAS data significantly improved the accuracy of the generated GF-7 DEMs. Figure 8 shows the spatial distribution and histogram of the differences in elevation between the GF-7 and LiDAR point clouds. The differences in elevation between the airborne LiDAR and GF-7 DEM Method1 point clouds mainly ranged from −10 m to 10 m, with a bias of 1.2 m and a standard deviation (STD) of 3.23 m. In contrast to GF-7 DEM Method1 , GF-7 DEM Method2 and GF-7 DEM Method3 yielded lower biases of -1.03 m and 0.03 m, and lower STD values of 2.05 m and 1.54 m, respectively. The average difference in the elevation of 0.03 m indicated that GF-7 DEM Method3 was highly consistent with the airborne LiDAR data, demonstrating that the use of GCPs extracted from ICESat-2/ATLAS data can decrease the differences in elevation between GF-7 DEM and airborne LiDAR data, and improve the accuracy of the elevation of GF-7 DEMs.     Method2 and GF-7 DEMMethod3 yielded lower biases of -1.03 m and 0.03 m, and lower STD value of 2.05 m and 1.54 m, respectively. The average difference in the elevation of 0.03 m indicate that GF-7 DEMMethod3 was highly consistent with the airborne LiDAR data, demonstrating tha the use of GCPs extracted from ICESat-2/ATLAS data can decrease the differences in elevatio between GF-7 DEM and airborne LiDAR data, and improve the accuracy of the elevation o GF-7 DEMs. To explore the ability of GF-7 DEMs to reproduce the slope, the values of the terrain slope were calculated by the DEMs. Figure 9 shows the distribution of the slopes of th terrain calculated by the airborne LiDAR DEM and GF-7 DEMs. The results indicate tha the slopes generated by the GF-7 DEMs agreed with those extracted from the airborn LiDAR data in most areas. The differences in the slope between the airborne LiDAR an GF-7 DEMs were also calculated, as shown in Figure 10. The results indicated that GF-DEMMethod1 and GF-7 DEMMethod3 were superior to GF-7 DEMMethod2 for reproducing th slope of the terrain. To explore the ability of GF-7 DEMs to reproduce the slope, the values of the terrain's slope were calculated by the DEMs. Figure 9 shows the distribution of the slopes of the terrain calculated by the airborne LiDAR DEM and GF-7 DEMs. The results indicate that the slopes generated by the GF-7 DEMs agreed with those extracted from the airborne LiDAR data in most areas. The differences in the slope between the airborne LiDAR and GF-7 DEMs were also calculated, as shown in Figure 10. The results indicated that GF-7 DEM Method1 and GF-7 DEM Method3 were superior to GF-7 DEM Method2 for reproducing the slope of the terrain.     Figure 11 compares the results of active faults, fault scarps, and horizontal fault o sets identified by airborne LiDAR-derived DEM with a 1 m resolution and the GF-7 DE with a 2 m resolution. As shown in Figure 11, most faults, vertical scarps, and horizon offsets could be effectively identified by the GF-7 DEM. The GF-7 DEM can be used identify horizontal offsets larger than 5 m (Figure 11f), while it cannot identify horizon offsets smaller than 3 m because of the limited resolution (Figure 11d). In addition, t fault scarps, such as V6, with a vertical offset of 0.89 ± 0.04 m, and V8, with a vertical off of 1.12 ± 0.16 m, could not be effectively identified on the basis of the GF-7 DEM. In sho the GF-7 DEM is suitable for identifying apparent horizontal and vertical offsets over lar areas, while it had difficulty identifying fine fault landforms.  Figure 11 compares the results of active faults, fault scarps, and horizontal fault offsets identified by airborne LiDAR-derived DEM with a 1 m resolution and the GF-7 DEM with a 2 m resolution. As shown in Figure 11, most faults, vertical scarps, and horizontal offsets could be effectively identified by the GF-7 DEM. The GF-7 DEM can be used to identify horizontal offsets larger than 5 m (Figure 11f), while it cannot identify horizontal offsets smaller than 3 m because of the limited resolution ( Figure 11d). In addition, the fault scarps, such as V6, with a vertical offset of 0.89 ± 0.04 m, and V8, with a vertical offset of 1.12 ± 0.16 m, could not be effectively identified on the basis of the GF-7 DEM. In short, the GF-7 DEM is suitable for identifying apparent horizontal and vertical offsets over large areas, while it had difficulty identifying fine fault landforms.   Figure 12 shows the vertical offsets of V1~V8 measured by the airborne LiDAR rived DEM, GF-7 DEMMethod1, GF-7 DEMMethod2, and GF-7 DEMMethod3. The GF-7 D could reconstruct the fault profile shape with a scarp height larger than 2 m, but sig cant errors could be observed in the fault profile shape reconstructed by the GF-7 D for small vertical scarps, such as V6 and V8. Figure 13 represents the error scatterplo tween the vertical offsets extracted from the airborne LiDAR-derived DEM and the different GF-7 DEMs. The bias and RMSE values of the vertical offsets measured by DEMMethod1 were -0.32 m and 0.55 m. Compared with the significant error in the elev of GF-7 DEMMethod1 (RMSE = 110.58 m), the accuracy of the vertical offsets measure GF-7 DEMMethod1 improved significantly, indicating that the absolute elevation error o 7 DEMMethod1 was large while the relative elevation error was small. The RMSE and values of the vertical offsets measured by GF-7 DEMMethod2 were the same as or sm than those of GF-7 DEMMethod1. Among these three GF-7 DEMs, the accuracy of the ve offsets measured by GF-7 DEMMethod3 is the highest, with a bias value of -0.19 m an RMSE value of 0.41 m. This means that the GCPs extracted from ICESat-2/ATLAS can slightly improve the accuracy of the measurements of the vertical offsets.  Figure 12 shows the vertical offsets of V1~V8 measured by the airborne LiDAR-derived DEM, GF-7 DEM Method1 , GF-7 DEM Method2 , and GF-7 DEM Method3 . The GF-7 DEMs could reconstruct the fault profile shape with a scarp height larger than 2 m, but significant errors could be observed in the fault profile shape reconstructed by the GF-7 DEMs for small vertical scarps, such as V6 and V8. Figure 13 represents the error scatterplot between the vertical offsets extracted from the airborne LiDAR-derived DEM and the three different GF-7 DEMs. The bias and RMSE values of the vertical offsets measured by GF-7 DEM Method1 were -0.32 m and 0.55 m. Compared with the significant error in the elevation of GF-7 DEM Method1 (RMSE = 110.58 m), the accuracy of the vertical offsets measured by GF-7 DEM Method1 improved significantly, indicating that the absolute elevation error of GF-7 DEM Method1 was large while the relative elevation error was small. The RMSE and bias values of the vertical offsets measured by GF-7 DEM Method2 were the same as or smaller than those of GF-7 DEM Method1 . Among these three GF-7 DEMs, the accuracy of the vertical offsets measured by GF-7 DEM Method3 is the highest, with a bias value of -0.19 m and an RMSE value of 0.41 m. This means that the GCPs extracted from ICESat-2/ATLAS data can slightly improve the accuracy of the measurements of the vertical offsets.

Comparison of the Measurements of the Vertical Offsets
From Figure 13, we can also see most of the measured vertical offsets were generally below the 1:1 line for all three DEMs, suggesting that the vertical offsets may have been underestimated compared with the airborne LiDAR DEM, especially for smaller offset values. The possible reason may be that GF-7 DEMs had a lower resolution (2 m) than the LiDAR DEM (1 m) and cannot reflect the subtle changes in the surface elevation. In other words, the GF-7 DEMs can be regarded as the result obtained after smoothing the LiDAR DTM. In this case, the variation in the elevation of the GF-7 DEMs was smaller than that of the LiDAR DTM, resulting in the vertical offsets being underestimated, especially for smaller offset values.
LiDAR DEM (1 m) and cannot reflect the subtle changes in the surface elevation. In other words, the GF-7 DEMs can be regarded as the result obtained after smoothing the LiDAR DTM. In this case, the variation in the elevation of the GF-7 DEMs was smaller than that of the LiDAR DTM, resulting in the vertical offsets being underestimated, especially for smaller offset values.   Figure 14 shows the horizontal offsets and errors of H1-H13 measured by the airborne LiDAR-derived DEM, GF-7 DEMMethod1, GF-7 DEMMethod2, and GF-7 DEMMethod3. The horizontal offsets extracted by the GF-7 DEMs were compatible with those extracted by the airborne LiDAR-derived DEM. Table 5 lists the R 2 , bias, and RMSE values of the horizontal offsets extracted by the airborne LiDAR-derived DEM and the GF-7 DEMs. The errors in the horizontal offsets measured by GF-7 DEMMethod1 were the highest, with a bias value of -1.80 m and an RMSE value of 3.98 m. The possible reason is that GF-7 DEMMethod1 had a specific deformation (the horizontal positioning errors in different locations were not constant). Compared with GF-7 DEMMethod1 and GF-7 DEMMethod2, the measured horizontal offsets obtained with GF-7 DEMMethod3 yielded the highest accuracy, with bias and the RMSE values of -0.81 m and 1.37 m, respectively. In other words, the GCPs extracted from ICESat-2/ATLAS data can improve the accuracy of the measurements of horizontal offsets.    Figure 14 shows the horizontal offsets and errors of H1-H13 measured by the airborne LiDAR-derived DEM, GF-7 DEM Method1 , GF-7 DEM Method2 , and GF-7 DEM Method3 . The horizontal offsets extracted by the GF-7 DEMs were compatible with those extracted by the airborne LiDAR-derived DEM. Table 5 lists the R 2 , bias, and RMSE values of the horizontal offsets extracted by the airborne LiDAR-derived DEM and the GF-7 DEMs. The errors in the horizontal offsets measured by GF-7 DEM Method1 were the highest, with a bias value of -1.80 m and an RMSE value of 3.98 m. The possible reason is that GF-7 DEM Method1 had a specific deformation (the horizontal positioning errors in different locations were not constant). Compared with GF-7 DEM Method1 and GF-7 DEM Method2 , the measured horizontal offsets obtained with GF-7 DEM Method3 yielded the highest accuracy, with bias and the RMSE values of -0.81 m and 1.37 m, respectively. In other words, the GCPs extracted from ICESat-2/ATLAS data can improve the accuracy of the measurements of horizontal offsets.  Figure 14 shows the horizontal offsets and errors of H1-H13 measured by the airborne LiDAR-derived DEM, GF-7 DEMMethod1, GF-7 DEMMethod2, and GF-7 DEMMethod3. The horizontal offsets extracted by the GF-7 DEMs were compatible with those extracted by the airborne LiDAR-derived DEM. Table 5 lists the R 2 , bias, and RMSE values of the horizontal offsets extracted by the airborne LiDAR-derived DEM and the GF-7 DEMs. The errors in the horizontal offsets measured by GF-7 DEMMethod1 were the highest, with a bias value of -1.80 m and an RMSE value of 3.98 m. The possible reason is that GF-7 DEMMethod1 had a specific deformation (the horizontal positioning errors in different locations were not constant). Compared with GF-7 DEMMethod1 and GF-7 DEMMethod2, the measured horizontal offsets obtained with GF-7 DEMMethod3 yielded the highest accuracy, with bias and the RMSE values of -0.81 m and 1.37 m, respectively. In other words, the GCPs extracted from ICESat-2/ATLAS data can improve the accuracy of the measurements of horizontal offsets.

Conclusions
This study quantitatively evaluated the accuracy of DEMs extracted from GF-7 stereo images with a 2 m resolution with or without GCPs and their potential for application to the quantitative study of active tectonics. The conclusions are as follows: (1) The GF-7 can provide large-scale high-resolution topographical and geomorphological data for quantitatively studying active structures. (2) The GF-7 DEM extracted without GCPs yielded high RMSE values of 178.24 m, 198.90 m, and 110.58 m in the X, Y, and Z directions, respectively. The GCPs extracted from ICESat-2/ATLAS data were used to significantly improve the RMSE values of the GF-7 DEMs in the X, Y, and Z directions to 1.38 m, 1.73 m, and 1.35 m. (3) On the basis of the GF-7 DEMs, the apparent fault scarps and horizontal offsets could be effectively identified, but the fine fault landforms could not be easily identified. (4) Although all three GF-7 DEMs achieved satisfactory accuracy in the measurements of the horizontal and vertical offsets, the GCPs extracted from the ICESat-2/ATLAS data could still improve the accuracy of the measurements of the horizontal and vertical offsets.
In summary, this research indicated, for the first time, the significant potential of China's first sub-meter stereo mapping satellite, GF-7, for studying active tectonics. However, this study also has some deficiencies, as follows. (1) Only few reliable GCPs were obtained from the ICESat-2 data in the areas with a higher elevation and a rough, steep topography, which may have affected the accuracy of the elevation of the generated GF-7 DEMs. (2) This study only selected the Hasishan section of the Haiyuan fault for research. Therefore, in-depth research should be conducted to explore how to obtain high-precision GCPs from ATL08 data in the area with a rough, steep topography and to analyze the performance of different global DEMs (e.g., the Advanced Land Observing Satellite World 3D (AW3D30) DEM, the Copernicus DEM) in assisting block adjustments of GF-7 data. In addition, we can also select global GF-7 data and sites for measuring the offsets of different levels to improve the reliability of the results.