4.1. Testing Site
Gongming Reservoir (as depicted in
Figure 10) was built along the surrounding hills with six adjacent earth-rock-filled dams forming a fan-shaped water area. The reservoir is designed to maintain a water level of 59.7 m, boasting a volumetric storage capacity of 142 million cubic meters and covering an expansive 6 square kilometers. All those dams were built on hilly terrain with homogeneous earth-filled structures, except for the No. 3 dam, which is a clay core earth and rock-filled dam. The total length of the dam crest is 4.3 km. The earth filling reaches heights ranging from 30 to 50 m. The maximum elevation of the dam peaks at 50.7 m.
Gongming Reservoir’s geographical location and layout are illustrated in
Figure 10;
Figure 10a presents an optical image, while
Figure 10b displays the DEM of the six dams acquired by UAV in December 2016. The resolution of the DEM is 0.45 m and the height accuracy is approximately 20 cm. Optical images of the six dams are depicted in
Figure 10c–h.
The upstream slope of each dam is constructed with a 0.3-m-thick concrete panel, while the downstream slope is covered with soil ranging from 0.3 to 0.5 m in thickness and adorned with planted grass. On the dam crest, there is a concrete road with an 8 m width (refer to
Figure 11b). Additionally, a 0.7 m high waterproof wall has been constructed near the upstream slope on the crest, and a 0.3 m high protective wall is positioned on the opposite grass slope side.
It took 8 years to build up all six dams of Gongming Reservoir. The construction of No. 1–2 dams started in 2010 and was completed in 2013. No. 3 dam was completed by the end of 2016. And in the end of 2017, No. 4, 5, and 6 dams were completed. Since July 2017, the reservoir started to impound. To prevent the leakage of the dam basement, concrete leakage-proof walls were built at the bottom of the dam center, and curtain grouting was carried out in the weathered rock zone. The six dams are mainly homogenous earth dams with slight structural differences. Only the No. 3 dam has a clay core for waterproofing.
The earth-filling height of the six dams ranges approximately from 30 to 50 m, and their post-construction settlement period typically spans 1 to 3 years. Based on the projected construction schedule, the primary post-construction settlement period of the No. 1 and No. 2 dams occurred from 2013 to 2014, while that of the No. 3 dam transpired from 2016 to 2018. No. 4, No. 5, and No. 6 dams experienced their main settlement periods from 2017 to 2019. As the multi-satellite SAR data used in this study have a considerable overlap in the first half of 2017, the paper focuses on studying No. 4 and No. 5 with large deformation during this period.
4.2. Dam Surface Scatterers and Coherence in Different SAR Sensors
We gathered three types of ascending SAR data of the Gongming reservoir including TSX-Spotlight (TSX-SL), COSMO-SkyMed stripmode (CSK-SM), and Sentinel-1 TOPS (S1), all acquired between March to June 2017. The detailed imaging parameters are outlined in
Table 2. The range resolution of TSX-SL, CSK-SM, and S1 are 0.45 m, 1.2 m, and 2.3 m, respectively. In the processing of Sentinel-1 TOPS data, a 4:1 multi-looking was conducted to maintain a similar aspect ratio of the dam, considering the substantial difference between azimuth and range resolutions.
Table 2 reveals that the temporal baselines of multi-satellite SAR data are approximately identical. According to the dam post-construction mechanism, the dam surface deformation trend during three months can be regarded as partially linear. Consequently, a time interval difference of a few days is not expected to introduce significant deformation values to the DInSAR interferograms.
Figure 12a–f present the SAR images of six dams obtained from ascending TSX-SL data. The scattering intensity on the concrete panel of the upstream slope appears notably weak, which is similar to the intensity of the calm water surface This phenomenon can be attributed to the specular reflection of the smooth concrete panel and its large local incident angle. Consequently, the coherence of the concrete slope is considerably low, leading to the conclusion that the reliable interference phase cannot be extracted from the concrete panels in the interferograms. In contrast, the scattering intensity of the downstream grass slope becomes pronounced when the dam axis is parallel to the SAR heading direction. Meanwhile, when the intersection angle
between the dam axis and the SAR heading direction increases, the grass slope scattering intensity decreases correspondingly.
Furthermore, the waterproofing on the crest of No. 4, 5, and 6 dams exhibit a distinct and discernible linear texture in SAR images, due to dihedral angle reflection, resulting in high scattering intensity. However, for the remaining dams, the angle between the waterproofing and the satellite flight is so large that the scattering intensity is weak and difficult to identify.
The SAR image of the six dams in ascending CSK-SM data is shown in the
Appendix A (
Figure A1). Despite the lower resolution of CSK-SM, approximately 2.5 times less than that of TSX-SL, the radar scattering characteristics of the concrete panel and grass slope remain similar. However, for Sentinel-1 -TOP data, the dams are almost invisible (refer to
Figure A2 in the
Appendix A for details). A 4:1 multi-look processing has been completed for a similar aspect ratio, and the resolution of Sentinel data is around 10 m both in azimuth and in range direction. Notably, the concrete and grass slopes, spanning approximately 100 m in length, correspond to merely 7–9 pixels in the SAR image. This limitation leads to an unclear outline of these slopes in Sentinel-1 images. However, as the scattering intensity of the grass slope is much higher than that of the concrete panel, the boundary between the grass slope and the concrete surface of the dam body can be visually distinguished.
Taking the No. 4 dam of Gongming Reservoir as an example,
Figure 13 illustrates the geometric distortion characteristics of the dam in various SAR images. The TSX-SL data, with a resolution of 0.45 m, can reflect the structure of the dam well, covering approximately 250 × 300 pixels. In the CSK-SM data, with a resolution of 1.18 m, the linear features such as waterproofing on the dam crest can still be distinguished, and the dam coverage is about 80 × 110 pixels. For Sentinel-1 TOPS data with about 10 m resolution (after the 4:1 multi-look) in range direction, it is difficult to distinguish the dam body, and the dam covers about 15 × 20 pixels. The linear texture features produced by the dihedral corner reflection of the waterproofing on the crest can be observed in both the TSX and CSK data, but it is difficult to distinguish between the dam crest and the slopes in the Sentinel data. The scattering intensity of the concrete panel in the X-band and C-band Sar data is quite low, which is similar to the calm water surface.
4.3. The Dams’ Settlement in Multi-DInSAR Images
Figure 14 shows ascending differential interferograms of six dams in the spring of 2017 using different SAR sensors. In general, the differential interferometric phase of the grass slope is well maintained. To illustrate the resolution differences among various SAR data, the images are presented based on the pixel size in SAR coordinates. The low resolution of Sentinel data makes the original differential interferometric image nearly invisible (We enlarge it five times for better visualization). The differential interferometric phase of high-resolution TSX exhibits the best clarity with clear fringe visibility. The visibility of CSK fringes is slightly worse than that of TSX, while the S1 differential interferometric phase of the grass slope can hardly be recognized. According to the construction progress of Gongming Reservoir, the No. 1 and No. 2 dams tend to be stable and there is no post-construction settlement effect in 2017. The differential interference fringes on the surfaces of dams No. 1 and No. 2 have the same color, confirming their stability. Because of the shallow filling of the No. 6 dam, it is difficult to recognize the deformation signal from the interferometric phase. There is a slight deformation on the No. 3 dam during this time. However, there is significant deformation on the No. 4 and No. 5 dams as seen in dense fringes there.
Figure 15a–f illustrate the coherence map of No. 4–5 dams based on multi-satellite SAR data. The coherence map derived from the Sentinel data is enlarged five times for better visualization. In
Figure 15, it is evident that the coherence for the grass slope is high. The coherence map of TSX data exhibits smoothness with consistently high values. As for CSK data, the grass slope shows high coherence except for the location near the dam crest. Although the Sentinel data shows good coherence between the two slopes of the dam, it is difficult to distinguish the difference between the grass and concrete slopes there. Especially, the top of the No. 5 dam shows obvious decorrelation, resulting in unreliable phase values in the differential interferograms.
Figure 15 reveals that the waterproofing on the crest of No. 4 and 5 dams presents linear texture and good coherence in TSX-SL and CSK-SM differential interferograms (refer to the larger version of
Figure 15a). However, the Sentinel-1 TOPS data cannot distinguish waterproofing and its texture characteristics due to its lower resolution.
In
Figure 16, we have applied the minimum cost flow method for phase unwrapping to the differential interferograms of No. 4 and No. 5 dams from
Figure 14. Pixels with coherence values below 0.6 and intensity below a specific threshold are filtered out. The initial value for phase unwrapping, located in the stable area at the bottom of the dams, was set to 0. The unwrapped interferogram of the Sentinel data is enlarged five times.
For a better visualization in
Figure 16, we superimpose the unwrapped differential interferogram on the SAR intensity image. The TSX data, with its highest resolution and coherence, exhibits the most distinct deformation fringes. It is evident that the bottom of the dam remains stable, while the dam crest suffers the largest settlement. The deformation of linear features (waterproofing) at the crest of No. 4 and No. 5 dams can also be clearly observed (refer to the enlarged views in
Figure 16a,b).
Figure 16c,d show the unwrapped interferograms of CSK which suffer more noise compared with TSX data due to its lower resolution, but the deformation characteristics extracted from the No. 4 and No. 5 dams are almost consistent with TSX. For the Sentinel data, the unwrapped phase at the top of the No. 5 dam is filtered out due to the low coherence. It is observed that the deformation gradient near the dam top is substantial, leading to decorrelation issues.
4.4. Dam Settlement Analysis in Different SAR Sensors
To illustrate the post-construction settlement of No. 4 and No. 5 dams in various SAR sensors, we extract deformation curves from the multi-satellite SAR interferograms. In
Figure 16, for the No. 4 dam, a longitudinal section line K1–K2 at the top of the grass slope and waterproofing K3–K4 is drawn. The profile line of J1–J2 is selected at the middle of the dam where the largest deformation occurs. Similarly, for the No. 5 dam, the longitudinal profile lines are defined on the top of the grass slope Q1–Q2 and on the waterproofing on the crest Q3–Q4. Also, a profile P1–P2 is drawn in the area with the most significant deformation to further analyze the details of dam deformation.
Figure 17 illustrates the deformation curves of the longitudinal section line and the cross-section line of the area with the largest deformation on the No. 4 dam, which are extracted from the high-quality points (coherence higher than 0.6) in the unwrapped TSX, CSK, and S1 interferograms. The temporal baselines of different SAR sensors are between March and June 2017, lasting for nearly three months, with a slight mismatch (no more than 15 days). The specific time is indicated in
Figure 17b,d,f. Considering the small difference in temporal baselines, the post-construction settlement characteristics of the dam based on multi-DInSAR results should be consistent. In addition, the deformation by water pressure can be ignored because the impoundment started after the SAR acquisition time.
Figure 17a,c,e presents LOS deformation results of the longitudinal section line at the top of the grass slope K1–K2 and waterproofing on the crest K3–K4, extracted from TSX, CSK, and S1 unwrapped differential interferograms, respectively. The y-axis on the left of
Figure 17a,c,e represents the embankment filling height of the dam clay (in m), depicted by a purple line. The y-axis on the right represents the LOS deformation (in mm). The blue points show the LOS deformation curve extracted along the K1–K2 profile line at the top of the grass slope, and the red point shows the LOS deformation curve extracted from the K3–K4 profile line of the waterproofing on the crest. There are no valid measurement points extracted from the concrete slope, due to the low intensity and decorrelation.
Figure 17a illustrates the result of the TSX data, where high-quality points extracted from K1–K2 and K3–K4 profile lines allow a reliable assessment of the deformation. The three-month post-construction settlement of the grass slope and the waterproofing exhibits the same deformation trend on the top axis of the dam, and this trend has a good correlation with filling height. The thicker the filling, the larger the deformation. It should be noted that the deformation curve of the K1–K2 section line at the top of the grass slope is separated, to a certain extent, from the K3–K4 section line deformation curve of the waterproof in the deeply filled area. This phenomenon is associated with the deformation conversion parameter ε (refer to Equation (8)) proposed in
Section 2.2.
Figure 17c presents the result of CSK data; with the number of points on the K1–K2 and K3–K4 profile lines being nearly half of that in TSX, which is expected due to the difference in resolution of the two types of data. The CSK result exhibits some volatility due to its greater phase noise. In general, the two-month deformation trend of the grass slope and the waterproofing manifested in the CSK data is highly consistent with the three-month deformation trend of TSX, indicating that deformation in the area with deeper filling is higher.
Figure 17e illustrates the results of the K1–K2 profile line of the grass slope of the Sentinel data. Only 21 sparse points are obtained from the profile line, due to the low resolution. Meanwhile, the phase of the K3–K4 profile line of waterproofing could not be extracted. The deformation trend of the grass slope in the Sentinel data is still correlated with the depth of the dam body filling, and it is consistent with the trends in TSX and CSK. The amount of deformation in the middle area of the dam body has accumulated to about 9 mm in 3 months.
Figure 17b,d,f present the deformation curves of the J1–J2 cross-section lines in the middle of the No. 4 dam, obtained from the unwrapped interferograms of TSX, CSK, and Sentinel. To minimize the influence of phase noise in individual points and obtain more samples, we set the width of the J1–J2 cross profile line to 40 m, which is indicated by the white bands along J1–J2 in
Figure 17b. The blue point represents the LOS deformation of the grass slope, the averaged value is denoted by the yellow dashed line. The averaged value of the grass slope near the crest is marked with black points. The red point indicates the LOS deformation of the waterproof on the crest. The black dashed line represents the fitting averaged value of the red points and is marked with text.
Figure 17b depicts the results obtained from the TSX data. The post-construction settlement of the No. 4 dam at the lower part of the grass slope is minimal within three months. The post-construction deformation of the dam body gradually increases from the middle of the grass slope, reaching a maximum of 15 mm (equal to 2π phase cycle) at the top area of the grass slope within three months. However, the deformation of the waterproofing on the crest is approximately 11 mm, which is 3.8 mm less than that of the grass slope.
Figure 17d displays the deformation results of the CSK data along the cross-section line. Similar to the TSX results, the post-construction settlement of the dam on the bottom of the grass slope exhibits a small amount of deformation within two months. Starting from the middle of the grass slope, the deformation of the dam body gradually increases and reaches 10 mm at the top of the grass slope, while that of the waterproofing is about 6 mm. This indicates a 4 mm gradient between them. Referring to Equation (9), due to the different local incident angles of two planes, the LOS deformation of the dam crest and the top of the grass slope exhibits variations in both TSX and CSK data. It is necessary to correct the deformation with reference to the local incident angles of different satellites on different dam planes.
Figure 17f illustrates the results obtained from the Sentinel data, revealing a maximum deformation of approximately 9 mm along the cross-sectional line. It can be vaguely distinguished that there is no settlement deformation in the lower part of the grass slope while the deformation in the middle and upper parts gradually increases. The Sentinel results exhibit partial agreement with TSX and CSK, indicating that the low resolution of Sentinel-1 can capture dam deformation under suitable incident angles, baselines, and slope angles. However, compared to the deformation captured by TSX within three months, the deformation captured by Sentinel appears somehow underestimated. Moreover, due to the low resolution of Sentinel, it is hard to distinguish the specific location of the waterproofing on the crest and extract its deformation phase.
Figure 18 presents the deformation curves of the longitudinal section and the cross-section of the area with the largest deformation on the No. 5 dam, obtained from the high-quality points in multi-satellite SAR unwrapped interferograms.
Figure 18a,c,e depict LOS deformation results of the longitudinal section of the top of the grass slope Q1–Q2 and waterproofing on the crest Q3–Q4, extracted from TSX, CSK, and Sentinel, respectively.
Figure 18b,d,f illustrate the deformation curves of the P1–P2 cross-section lines in the middle of the No. 5 dam extracted from TSX, CSK, and Sentinel, respectively. The x-axis and y-axis settings, as well as the temporal baseline for SAR sensors, remain consistent with those in
Figure 17.
In
Figure 18b,d,f, the deformation curves of the P1–P2 cross-section lines in the middle of the No. 5 dam are displayed, extracted from the unwrapped interferograms of TSX, CSK, and Sentinel. The x-axis and y-axis settings, along with the color-coded points, maintain consistency with those in
Figure 17b. The width of the cross-section line is set at 50 m.
Regarding the TSX result presented in
Figure 18b, the deformation along the profile line of the grass slope is consistently depicted, indicating favorable phase quality for TSX. The settlement observed at the top of the grass slope measures approximately 20 mm (equivalent to 2.6π phase cycles). Notably, there exists a noticeable gradient difference of approximately 6.2 mm between the yellow dashed line, representing the averaged deformation at the top of the grass slope, and the black dashed line, signifying the averaged deformation at the crest.
In
Figure 18d, the cross-sectional deformation of the No. 5 dam obtained by CSK is presented. The deformation trend of the grass slope and the crest are consistent with TSX. Considering that the temporal baseline of CSK is less than TSX for about one month, the maximum deformation by CSK is reduced proportionally. Nonetheless, a deformation disparity of around 6.6 mm still exists between the top of the grass slope and the crest.
Figure 18f illustrates the deformation curve of the cross-section of the No. 5 dam obtained by the Sentinel data. This result significantly differs from the TSX and CSK findings. Notably, in
Figure 18f, the phase noise at the bottom of the grass slope is prominent, and points in the middle and upper sections of the grass slope are masked due to coherence loss due to large deformation. Furthermore, the low resolution of the Sentinel data hinders the extraction of deformation information at the waterproofing on the crest.
The LOS deformation results for the No. 4 and No. 5 dams highlight a noticeable correlation between post-construction settlement and earth-filling height within the first year after construction. Approaching the top of the dam, the deformation becomes greater. The result also reveals the difference between the deformation curve of the top of the grass slope and that of the waterproofing on the crest. With reference to the previous simulation and derivations, the difference could be attributed to the conversion parameter
. Therefore, utilizing Equation (9) and considering the deformation difference between the grass slope and crest from
Figure 17b,d and
Figure 18b,d, it is possible to perform an inverse calculation to determine the actual conversion parameter
. This derived value can then be validated against the theoretical one.
The theoretical conversion parameter
of slopes can be calculated by Equation (8), according to the dam slope angle and orientation from DEM.
Figure 19 presents the theoretical conversion parameter values of the six dam slopes. The conversion parameter values are labeled through a blue-red color bar from 0.5 to 1.0. Since the incident angles of the TSX and CSK data only differ by 1.4° in this case, the deformation conversion parameters of the two data are similar. Therefore, only the TSX results are displayed.
In
Figure 19, all the grass slopes for the dams have a slope angle of 20°, while the concrete slopes have a consistent slope angle of 18.5°. Considering that most grass slopes directly face the SAR LOS direction, the conversion parameter for the grass slope is about 0.9. However, as most concrete slopes deviate from the radar LOS direction, the conversion parameters for these slopes are mainly less than 0.6. Specifically, for the No. 1 dam, its axis direction is nearly perpendicular to the flight direction of the satellite, resulting in ϕ being close to 90°. The local incident angles and the conversion parameters of the grass slope and the concrete slope are almost identical. Regarding the No. 2 dam, the zigzag design of the dam results in four slopes with different orientations, each corresponding to a distinct conversion parameter. For the crest of all dams, the deformation conversion parameter remains constant regardless of dam orientation and is consistently 0.8.
Similar to the simulated dam deformation in
Section 3, we assume that the deformation on the top of the grass slope and the crest are consistent. Subsequently, the converted deformation of the crest is considered the truly vertical deformation of the top of the grass slope, and the practical conversion parameter
of the grass slope can be calculated by the observed LOS deformation and converted vertical deformation on the top of the grass slope.
Table 3 shows the results of the theoretical and practical conversion parameters
ε of the No. 4 and No. 5 dams. It also indicates the consistency between the theoretical and practical values, which is calculated by
. The theoretical deformation conversion parameters for the No. 4 and No. 5 dams are the same due to their identical slope angle and orientation. The conformity of the theoretical and practical conversion parameters
are higher than 80% for both the TSX and CSK data. The conformity of the CSK data is slightly higher than that of the TSX data in both dams, suggesting that the real dam post-construction settlement phenomenon may not be as ideal as estimated and simulated, and the discrepancy does show some highlights for us.