Spatiotemporal Inversion and Mechanism Analysis of Surface Subsidence in Shanghai Area Based on Time-Series InSAR

: Abnormal surface subsidence has become a widespread geological problem being faced by cities. As a ﬁrst-tier city in the world, there are problems such as unclear analysis of subsidence mechanism. Therefore, the monitoring of the Shanghai surface is particularly important. In this paper, the Sentinel-1A satellite SAR image data of 36 scenes covering Shanghai area from January 2018 to March 2020 were processed on the basis of time-series interferometry synthetic aperture radar technique. The subsidence rate ﬁeld and accumulated surface subsidence in Shanghai area during the study period were obtained, and the spatial–temporal distribution characteristics of subsidence in the study area were discussed and analyzed from many different aspects. Moreover, the correlation between regional subsidence and geological structure, precipitation, urbanization, and other inﬂuencing factors were analyzed and established. Results show that the non-uniform subsidence in Shanghai area is clear, and those of the districts of eastern Songjiang, southern Jinshan, and Fengxian are more serious with a maximum subsidence rate of − 26.2 mm/year. After analyzing the causes of subsidence, the special foundation of soft soil in Shanghai area is determined as the main reason for the subsidence, and the uneven subsidence is mainly caused by the over-exploitation of groundwater, human activities, and the subsidence of soil layer. Through comparative analysis, factors, such as rainfall, groundwater, and urbanization process, have high correlation with surface subsidence.


Introduction
As a serious geological disaster, surface subsidence will directly lead to major safety accidents, such as the collapsing of grounds and houses, which seriously threatens the safety of human life and property [1,2]. With the rapid development of urbanization in China, the surface non-tectonic deformation caused by human activities (e.g., natural resource exploitation and underground space engineering exploitation) has become a major geological disaster being faced by China [3]. According to the data, abnormal subsidence, which has caused great losses to people's life and property safety, has occurred in many areas in China (e.g., Beijing [4], Foshan [5]). Therefore, surface subsidence monitoring is particularly important. Traditional subsidence monitoring is mainly based on leveling and GNSS measurement [6,7], and these technologies are mainly based on the target of monitoring points and lack certain timeliness in the observation process [8,9]. For the monitoring of large-scale subsidence in some Chinese cities, the aforementioned technologies have the disadvantages of low spatial resolution, long period, and high cost [10,11].
In the 1960s, interferometry synthetic aperture radar (InSAR) technique achieved rapid development. This technique has been widely used in mine subsidence [12], volcano monitoring [13], landslide monitoring [14], earthquake [15], dam [16], and other aspects by virtue of its all-day, all-weather, and wide coverage. With the development of InSAR technique, differential interferometry synthetic aperture radar (DInSAR) technique was introduced. In 1989, Gabriel et al. [17] obtained the surface displacement of the Impirale Valley, California, USA by using DInSAR technique. However, this technique still has shortcomings (e.g., incoherence and atmospheric delay) that will directly lead to large errors in monitoring results. To solve the shortcomings of DInSAR, several time-series InSAR techniques have been proposed. In 2000, Ferretti et al. [18] proposed a persistent scatterer InSAR (PSInSAR) technique that works by calculating some persistent scatterer (PS) points during the study period. Then, the influence of atmospheric effect and incoherent factors is reduced. Unlike PSInSAR, the technique proposed by Berardino [19] in 2002, which was called small baseline subset InSAR (SBASInSAR), changed the point target modeling method into distributed target modeling analysis. In recent years, with the development of time-series InSAR technique, various time-series InSAR technologies (e.g., MTInSAR [20], TCPInSAR [21]) have emerged, and these time-series InSAR technologies all have good application prospect in surface subsidence monitoring. For example, Calof et al. [22] analyzed the driving factors of surface deformation in Kongya area by using SBASInSAR technique and combining the data of climate, strata, and land cover change. Zhou et al. [23] used SBASInSAR to investigate wide-area surface subsidence characteristics in Wuhan and analyzed and established the correlations between surface subsidence and urban construction, precipitation, industrial development, carbonate karstification, and water level changes in Yangtze River in Wuhan. Brunori et al. [24] used ENVISAT ASAR and RADARSAT-2 data to detect obvious subsidence in the northwest of Ciudad Guzmán. After field investigation, the feasibility of multi-temporal InSAR technique in surface subsidence was verified.
The subsidence monitoring of infrastructure in Shanghai has always been a research hotspot of scholars at home and abroad. Yang et al. [25] monitored Shanghai's reclamation area (Lingang New Area) based on time series InSAR by using various radar data and showed that the ground deformation after reclamation in Shanghai experienced three obvious stages: primary consolidation, slight rebound, and platform stage. Wang et al. [26] processed 49 TerraSAR-X images with TSInSAR technique, obtained the subsidence time series of viaduct in Shanghai, and analyzed it in combination with construction time, load, and other factors. Zhang et al. [27] extracted the average subsidence rate of coherent targets in line of sight by using the ALOSPALSAR data that covered the Shanghai area from 2007 to 2010 with the small baseline time series analysis method and verified the feasibility of SBAS technique in Shanghai by comparing with the leveling data. The obvious subsidence funnel in Shanghai area has been detected. Antonio et al. [28] used SBASInSAR technique to study the temporal evolution characteristics of surface deformation in Shanghai land reclamation area from 2007 to 2016. Combined with C-band and X-band data, the average cumulative deformation of Shanghai land reclamation from 2007 to 2016 reached 25 cm. Lei et al. [29] obtained the deformation situation of Shanghai reclamation area in 2016 by using COSMS-SkyMed and Sentinel-1A data and obtained that the subsidence rate of Shanghai reclamation area was relatively high with a maximum subsidence rate reaching 30 mm/year. Qin et al. [30] used PSInSAR technique to monitor the main traffic infrastructure in Shanghai area and compared with the level data. They found that some traffic lines in Shanghai had serious subsidence and verified the feasibility of PSInSAR in subsidence monitoring in Shghai area. Zhao et al. [31] used InSAR technique to monitor Shanghai Lupu Bridge and verified the feasibility of PSInSAR technique in deformation monitoring of a single building. The subsidence monitoring in Shanghai is mainly focused on the local monitoring of infrastructure and other aspects. Studies on the large-scale surface subsidence monitoring in Shanghai are relatively few, and discussions on the subsidence mechanism are lacking. Moreover, the research data are mostly from 5 years ago, thereby lacking certain timeliness.
In this paper, we first processed and analyzed the Sentinel-1A image data of 36 scenes based on PSInSAR technique, and then obtained the deformation rate field from January 2018 to March 2020 in Shanghai. Second, we investigated the surface subsidence spatialtemporal characteristics in Shanghai area in detail. Finally, the driving factors of subsidence were analyzed in combination with urbanization construction, groundwater level change, rainfall, and other factors in Shanghai. These factors provide scientific data support for the surface subsidence prevention and control in Shanghai area.

Methodology
Aiming at the problems of interference incoherence and atmospheric delay error that exists in DInSAR technique, PSInSAR technique can reduce the influence of interference incoherence and atmospheric delay error on the calculation results effectively. Unlike DInSAR technique, PSInSAR focuses on the use of SAR image sequences acquired from the same area during a certain period. First, according to the scattering characteristics of ground objects and statistical analysis method, the permanent scatterers with high correlation in time series in the study area are detected. Second, based on the phase time series of these specific targets, modeling and analysis are carried out to obtain high-precision deformation measurement results. The main steps are detailed as follows [32][33][34]: (1) Assume the N+1 SAR images covering the same area. Arrange the SAR images according to the time series. Calculate the temporal baseline, spatial baseline, and doppler centroid frequency difference between the images. Establish the relevant model according to the three aforementioned factors to select the reference image.
where ρ m is the comprehensive correlation coefficient of an interference pair; B k,m ⊥ , T k,m , and B k,m DC are the perpendicular baseline, temporal baseline, and Doppler centroid frequency difference between interference images k and m, respectively. The critical conditions of the aforementioned parameters are B c , T c , and f c , respectively. The image with the maximum value is the reference image to select the optimal reference image.
(2) After calculation, the image on 29 June 2019 was selected as the reference image in this paper. All the secondary images were registered to the reference image, and N interference pairs were generated after differential interference processing. The spatialtemporal baseline distribution is shown in Figure 1, and the specific spatial-temporal baseline distribution parameters of Sentinel-1A data are listed in Table 1.
(3) The amplitude deviation index method was used to select the high-amplitude value as the candidate PS point according to the strong reflection characteristics of PS points. The algorithm of amplitude deviation index is shown in (2). The interference phase information contained in the candidate PS points can be expressed as Equation (3) [35,36]: where D is the amplitude deviation index, σ i is the amplitude standard deviation of the images with time series, and m i is the mean amplitude of the time series image. In Equation (3), φ is the interference phase; φ re f represents the flat-earth phase; φ de f refers to the surface deformation phase; φ topo corresponds to the phase caused by topographic relief; φ atm denotes atmospheric phase error; and φ noi depicts the rando//m noise phases (e.g., orbital errors, thermal noise, and spatial decorrelation). Among them, φ re f can be eliminated by precise orbit data, φ topo can be eliminated by external DEM. The red lines depict interferometric pairs. The yellow and green diamonds correspond to the reference image and secondary images, respectively. (4) After separating the linear deformation and elevation error from using the correlation characteristics of atmospheric phase with low correlation in time and high correlation in space, the high-pass filtering in time and low-pass filtering in space are used to obtain the estimated values of atmospheric phase and nonlinear deformation in each scene of the differential interferogram.
where ∆ represents the interferometric phase after removing the flat-earth and topographic phases; and and denote the linear and nonlinear surface deformation phases, respectively. refers to the deformation phase. The red lines depict interferometric pairs. The yellow and green diamonds correspond to the reference image and secondary images, respectively. (4) After separating the linear deformation and elevation error from φ using the correlation characteristics of atmospheric phase with low correlation in time and high correlation in space, the high-pass filtering in time and low-pass filtering in space are used to obtain the estimated values of atmospheric phase and nonlinear deformation in each scene of the differential interferogram.
where ∆φ di f f represents the interferometric phase after removing the flat-earth and topographic phases; and φ lin and φ non denote the linear and nonlinear surface deformation phases, respectively. φ de f refers to the deformation phase.
The deformation phase can be obtained by adding the nonlinear and linear deformation, and then the deformation information in the study area can be solved.

Study Area
Shanghai (120 • 52 E-122 • 12 E, 30 • 40 N-31 • 53 N) is located in the easternmost part of the Yangtze River Delta in China. The average altitude is approximately 4 m above sea level. It is at the confluence of the Yangtze River and the Huangpu River. The place borders the Yangtze River to the north and the East China Sea to the east. Shanghai is a typical soft soil area with complex geological structure and large pore ratio. A large number of infrastructure construction can easily lead to abnormal subsidence in the area. Many faults exist in the Shanghai area, but these faults are mainly the early and middle Pleistocene faults and are classified as inactive. Therefore, the geological structure of this area is relatively stable, and it is not easy for major structural changes to occur.
Shanghai is one of the largest industrial bases in China. Many scientific research and production bases of light and heavy industrial products are located in Shanghai. Shanghai has a population of 24.28 million, thereby making it the second most populous city in China by 2019. The geographical location of the study area in this paper is represented by the area in the red rectangle in Figure 2a, which contains most of Shanghai. Figure 2b shows the backscattering intensity image of Sentinel-1A SAR image that covers the study area, and the measures in east-west and south-north directions are approximately 79.46 km and 78.43 km, respectively. The background shows the Google image that covers Shanghai area.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 5 of 22 (5) The deformation phase can be obtained by adding the nonlinear and linear deformation, and then the deformation information in the study area can be solved.

Study Area
Shanghai (120°52′ E-122°12′ E, 30°40′ N-31°53′ N) is located in the easternmost part of the Yangtze River Delta in China. The average altitude is approximately 4 m above sea level. It is at the confluence of the Yangtze River and the Huangpu River. The place borders the Yangtze River to the north and the East China Sea to the east. Shanghai is a typical soft soil area with complex geological structure and large pore ratio. A large number of infrastructure construction can easily lead to abnormal subsidence in the area. Many faults exist in the Shanghai area, but these faults are mainly the early and middle Pleistocene faults and are classified as inactive. Therefore, the geological structure of this area is relatively stable, and it is not easy for major structural changes to occur.
Shanghai is one of the largest industrial bases in China. Many scientific research and production bases of light and heavy industrial products are located in Shanghai. Shanghai has a population of 24.28 million, thereby making it the second most populous city in China by 2019. The geographical location of the study area in this paper is represented by the area in the red rectangle in Figure 2a, which contains most of Shanghai. Figure 2b shows the backscattering intensity image of Sentinel-1A SAR image that covers the study area, and the measures in east-west and south-north directions are approximately 79.46 km and 78.43 km, respectively. The background shows the Google image that covers Shanghai area. Due to its geographical location, Shanghai enjoys convenient transportation, advanced culture, and rapid industrial development. In recent years, many infrastructure and rail transits have been constructed in Shanghai, thereby signifying an uneven subsidence phenomenon. Shanghai has a subtropical monsoon climate with abundant rainfall, and the typical rainy season occurs in September of each year, with approximately 60% of the annual precipitation. The rainfall change has a great impact on regional surface subsidence, and the seasonal variation of subsidence is evident. To alleviate the excessive groundwater in the region, the Shanghai municipal government adopts tap water Due to its geographical location, Shanghai enjoys convenient transportation, advanced culture, and rapid industrial development. In recent years, many infrastructure and rail transits have been constructed in Shanghai, thereby signifying an uneven subsidence phenomenon. Shanghai has a subtropical monsoon climate with abundant rainfall, and the typical rainy season occurs in September of each year, with approximately 60% of the annual precipitation. The rainfall change has a great impact on regional surface subsidence, and the seasonal variation of subsidence is evident. To alleviate the excessive groundwater in the region, the Shanghai municipal government adopts tap water recharge to maintain the basic balance of groundwater. It has also become one of the measures to alleviate the surface subsidence.

Data
The 36 ascending SAR images used in this paper are from Sentinel-1A launched by the European Space Agency (ESA) (https://scihub.copernicus.eu/dhus/#/home, 20 May 2020), the satellite carries C-band synthetic aperture radar antenna with a 12-day revisit cycle. It can monitor the ground micro-deformation effectively. The time span of the 36 SAR images are from January 2018 to March 2020. The incidence angle of the image center is 36.6 • . All of the SAR images are in VH polarization with a resolution of 1.2 m in range and 13.9 m in azimuth. Specific image parameters are illustrated in Table 2. The three arc-second Shuttle Radar Topography Mission (SRTM) DEM provided by the National Aeronautics and Space Administration (NASA) (http://srtm.csi.cgiar.org/srtmdata, 2 June 2020) was used to remove the topographic phase. We adopted the ESA-released precise orbit data (https://qc.sentinel1.eo.esa.int/aux_poeorb, 22 June 2020) to execute orbital refinement and phase re-flattening for all interferometric pairs. The rainfall data were provided by China National Meteorological Science Data Center, the groundwater level and the position of water level monitoring wells were provided by Shanghai Geological Data Information Sharing Platform, and the night-light remote sensing data were the monthly composite product of NPP-VIIRS. The above data were used to analyze the mechanisms of surface subsidence in the study area. The spatial distribution of data is shown in Figure 3. recharge to maintain the basic balance of groundwater. It has also become one of the measures to alleviate the surface subsidence.

Data
The 36 ascending SAR images used in this paper are from Sentinel-1A launched by the European Space Agency (ESA) (https://scihub.copernicus.eu/dhus/#/home, 20 May 2020), the satellite carries C-band synthetic aperture radar antenna with a 12-day revisit cycle. It can monitor the ground micro-deformation effectively. The time span of the 36 SAR images are from January 2018 to March 2020. The incidence angle of the image center is 36.6°. All of the SAR images are in VH polarization with a resolution of 1.2 m in range and 13.9 m in azimuth. Specific image parameters are illustrated in Table 2. The three arcsecond Shuttle Radar Topography Mission (SRTM) DEM provided by the National Aeronautics and Space Administration (NASA) (http://srtm.csi.cgiar.org/srtmdata, 2 June 2020) was used to remove the topographic phase. We adopted the ESA-released precise orbit data (https://qc.sentinel1.eo.esa.int/aux_poeorb, 22 June 2020) to execute orbital refinement and phase re-flattening for all interferometric pairs. The rainfall data were provided by China National Meteorological Science Data Center, the groundwater level and the position of water level monitoring wells were provided by Shanghai Geological Data Information Sharing Platform, and the night-light remote sensing data were the monthly composite product of NPP-VIIRS. The above data were used to analyze the mechanisms of surface subsidence in the study area. The spatial distribution of data is shown in Figure  3.

Results and Analysis of Subsidence Rate
After processing the 36 Sentinel-1A SAR images, the subsidence rate field and accumulated subsidence in Shanghai from January 2018 to March 2020 were obtained. Figure 4 shows the surface subsidence rate map in Shanghai area derived by PSInSAR (the positive values indicate surface uplift, and the negative values represents surface subsidence), and the map is superimposed on a Google Earth image of Shanghai acquired in 2020.   Figure 5 shows the cumulative subsidence time series in Shanghai from 2018 to 2020. During the study period, the maximum cumulative subsidence in Shanghai is −68 mm, and it is located in Region c in the south of Fengxian District. As a whole, no large-scale subsidence is observed in Shanghai area during 2018-2019, and the subsidence has been slow. However, after 2019, the subsidence in Shanghai has shown an accelerating trend, which may be affected by the regional rainfall and the rainy season in the lower reaches of the Yangtze River.

Results and Analysis of Subsidence Time Series
The subsidence in the northern part of Shanghai is relatively stable, and the accumulated subsidence range from −10 mm to 20 mm. Figure 5 presents that the subsidence in the northern part of Shanghai has been showing a slow subsidence trend without serious large-scale subsidence. In the middle part of Shanghai, the subsidence near Songjiang District is more serious, with the maximum accumulated subsidence reaching −44 mm. In the cumulative time series, the Songjiang District has been showing a downward trend, with a relatively stable settlement range but without any obvious diffusion trend. Jinshan District and Fengxian District in the south of Shanghai have the most obvious subsidence phenomenon, the subsidence amount is gradually increasing, and the subsidence area is gradually expanding. The subsidence is relatively stable in most areas of Shanghai. However, the uneven subsidence phenomenon is obvious. Region a is located in the urban area of Shanghai. The subsidence in this region is relatively stable and without obvious subsidence. The main subsidence areas in Shanghai include Region b (SJ), Region c (FX), and Region d (JS). Large-scale subsidence occurs in the three regions. The subsidence range of Region b presents an expanding trend, and the subsidence rates range from −15.0 mm/year to 5mm/year. The overall subsidence rate in Region d is relatively high, which are between −25 mm/year and 6 mm/year. Many large industrial areas, such as Jinfengyuan Paper (JFYP), Yadong Petrochemical (YDP), Shanghai Paimilei Industrial Park (SPIP), are found in Region c. Moreover, this region is the most serious subsidence area in Shanghai with multiple significant subsidence funnels.
Region a is the main urban area of Shanghai. As shown in the subsidence rate map, the subsidence rate in this area is around −3 mm/year, and the subsidence is relatively stable without significant subsidence. According to the "Shanghai City Master Plan (1999-2020)" released by the Shanghai Municipal Urban Planning Administration, Shanghai has made strict restrictions on the population and urban construction land in the main urban area. Moreover, the central urban area has been built earlier, in which the subsidence has tended to stabilize.
Region b is located in Songjiang District, Shanghai. According to the "Shanghai Songjiang District Master Plan and Land Use Master Plan (2017-2035)" published by the People's Government of Songjiang District of Shanghai, Songjiang District will adhere to the development trend of combining multiple towns and Songjiang New City as the center. Under the condition of ensuring the steady growth of the regional ecological area, the construction of urban infrastructure and manufacturing factories should be strengthened.
These conditions have led to the rapid development of the industry in Songjiang District from 2017 to 2020. Population data show that by 2020, the region's permanent population reached 1582,400. Industrial and domestic water mainly come from the exploitation of underground water. In Section 5.3, the relationship between the surface subsidence and the change of the water table will be further discussed. Figure 5 shows the cumulative subsidence time series in Shanghai from 2018 to 2020. During the study period, the maximum cumulative subsidence in Shanghai is −68 mm, and it is located in Region c in the south of Fengxian District. As a whole, no large-scale subsidence is observed in Shanghai area during 2018-2019, and the subsidence has been slow. However, after 2019, the subsidence in Shanghai has shown an accelerating trend, which may be affected by the regional rainfall and the rainy season in the lower reaches of the Yangtze River. To assess the internal precision of the subsidence rate field extracted from Sentinel-281 1A SAR data processed by PSInSAR technique in Shanghai area, the standard deviation 282 of the obtained subsidence rate is statistically analyzed in this section. Figure 6 shows the 283 distributions of the subsidence rate standard deviations of all PS points in the study area. 284 According to the statistical results, 637213 PS points were detected in the study area, 98% 285 of PS points standard deviation are less than 1 mm/yr, and the maximum standard devi-286 ation reaches 2.28mm/yr. In conclusion, PSInSAR technique exhibits high reliability and 287 stability in monitoring surface subsidence in Shanghai area. The subsidence in the northern part of Shanghai is relatively stable, and the accumulated subsidence range from −10 mm to 20 mm. Figure 5 presents that the subsidence in the northern part of Shanghai has been showing a slow subsidence trend without serious large-scale subsidence. In the middle part of Shanghai, the subsidence near Songjiang District is more serious, with the maximum accumulated subsidence reaching −44 mm. In the cumulative time series, the Songjiang District has been showing a downward trend, with a relatively stable settlement range but without any obvious diffusion trend. Jinshan District and Fengxian District in the south of Shanghai have the most obvious subsidence phenomenon, the subsidence amount is gradually increasing, and the subsidence area is gradually expanding.

Accuracy Evaluation of Time Series InSAR Results
To assess the internal precision of the subsidence rate field extracted from Sentinel-1A SAR data processed by PSInSAR technique in Shanghai area, the standard deviation of the obtained subsidence rate is statistically analyzed in this section. According to the statistical results, 637,213 PS points were detected in the study area, 98% of PS points standard deviation are less than 1 mm/year, and the maximum standard deviation reaches 2.28 mm/year. In conclusion, PSInSAR technique exhibits high reliability and stability in monitoring surface subsidence in Shanghai area.

Time Series Analysis of Feature Points
To analyze the temporal and spatial characteristics of subsidence distribution in Shanghai better, this section randomly selects six feature points (i.e., PS1-PS6) scattered in different areas of Shanghai for specific analysis. The specific distribution of feature points is shown in Figure 3, and the accumulated subsidence amount of six time series points is shown in Figure 6.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 10 of 22 PS4 is located near Fei Shangwangjia in Fengxian District. There are a large number of large factories such as aluminum products and machinery parts production plants. Industrial production and industrial water using become one of the main factors for subsidence. The maximum accumulated subsidence in this area is −44.6 mm, and the subsidence phenomenon is more obvious.
PS5 is located near Shanghai Keen Machinery Manufacturing, Fengxian District. The surface of this area shows a continuous subsidence trend, and many machinery part production factories are found in this area, which becomes one of the main reasons for the subsidence in this area.
PS6 is located near Kangqiao Peninsula in Pudong, and the subsidence in this area is basically stable without obvious subsidence.
In general, the temporal and spatial inhomogeneity of subsidence in Shanghai is obvious, which shows uplift in the east and decline in the west and uplift in the north and decline in the south. The surface subsidence in Shanghai is nonlinear because of the overexploitation of groundwater in the region.

Influence of Regional Stratigraphic Structure on Surface Subsidence
Shanghai, which is located in the southeast front of the estuary of the Yangtze River Delta, is one of the world's famous mega cities located on soft soil. It belongs to the coastal plain landform type, and its geological age is Holocene Q4 and Upper Pleistocene Q3. As the longest river in Asia, the main stream of the Yangtze River flows through 11 provinces and cities before reaching the Huangpu River in Shanghai and entering the sea. The size of soil particles carried by the Yangtze River changed significantly in the process of flowing through various regions. Shanghai has suffered from the long-term impact of the Yangtze River estuary, and a huge amount of land was lost every year. According to the "Shanghai Water Resources Bulletin in 2019," the total soil erosion area of Shanghai in 2019 is 2.89 square kilometers, which accounts for 0.05% of the total land area of the city. Moreover, the sediment carried by the Yangtze River to Shanghai is mainly concentrated in particles of 0.005 mm to 0.075 mm. The particles, which are typical clay particles and accumulate to form silty soil, are relatively small.
To analyze the influence of geological structure in Shanghai on regional subsidence better, we obtained the geological data of two Quaternary Geological Boreholes (the green pentagram in Figure 3 In terms of Quaternary geology, the soil layer in Shanghai is evenly distributed. The upper layer is silty powdery clay, and the lower layer is silty clay, which together PS1 is located in the vicinity of Sichuan North Road Park in Hongkou District, which was constructed earlier. The area is the location of the old city of Shanghai, and the subsidence here has basically stabilized. Therefore, the subsidence of PS1 is relatively stable during the study period, and no excessive subsidence phenomenon occurred.
PS2 is located near Xinglian Central Village, Jiuting Town, Songjiang District. In the subsidence area, the subsidence fluctuation is relatively small, and no obvious subsidence phenomenon occurred.
PS3 is located near the Pailoutou Building in Minhang District. The point shows a slow subsidence trend with an accumulated subsidence of −14.4 mm. Many industrial parks and factories are found in this area. Moreover, the industrial water is extracted from groundwater, which is the main cause of subsidence.
PS4 is located near Fei Shangwangjia in Fengxian District. There are a large number of large factories such as aluminum products and machinery parts production plants. Industrial production and industrial water using become one of the main factors for subsidence. The maximum accumulated subsidence in this area is −44.6 mm, and the subsidence phenomenon is more obvious.
PS5 is located near Shanghai Keen Machinery Manufacturing, Fengxian District. The surface of this area shows a continuous subsidence trend, and many machinery part production factories are found in this area, which becomes one of the main reasons for the subsidence in this area.
PS6 is located near Kangqiao Peninsula in Pudong, and the subsidence in this area is basically stable without obvious subsidence.
In general, the temporal and spatial inhomogeneity of subsidence in Shanghai is obvious, which shows uplift in the east and decline in the west and uplift in the north and decline in the south. The surface subsidence in Shanghai is nonlinear because of the over-exploitation of groundwater in the region.

Influence of Regional Stratigraphic Structure on Surface Subsidence
Shanghai, which is located in the southeast front of the estuary of the Yangtze River Delta, is one of the world's famous mega cities located on soft soil. It belongs to the coastal plain landform type, and its geological age is Holocene Q4 and Upper Pleistocene Q3. As the longest river in Asia, the main stream of the Yangtze River flows through 11 provinces and cities before reaching the Huangpu River in Shanghai and entering the sea. The size of soil particles carried by the Yangtze River changed significantly in the process of flowing through various regions. Shanghai has suffered from the long-term impact of the Yangtze River estuary, and a huge amount of land was lost every year. According to the "Shanghai Water Resources Bulletin in 2019," the total soil erosion area of Shanghai in 2019 is 2.89 square kilometers, which accounts for 0.05% of the total land area of the city. Moreover, the sediment carried by the Yangtze River to Shanghai is mainly concentrated in particles of 0.005 mm to 0.075 mm. The particles, which are typical clay particles and accumulate to form silty soil, are relatively small.
To analyze the influence of geological structure in Shanghai on regional subsidence better, we obtained the geological data of two Quaternary Geological Boreholes (the green pentagram in Figure 3 In terms of Quaternary geology, the soil layer in Shanghai is evenly distributed. The upper layer is silty powdery clay, and the lower layer is silty clay, which together constitute the typical soft soil in Shanghai. In the soft soil layer, the soil moisture content is higher, and the hardness is relatively soft. Geological layers in Shanghai are mainly divided into artificial filled soil, cohesive soil, silty clay, and powder silt soil. The specific soil layer distribution of boreholes a and b are shown in Figure 7. The soil layer in Shanghai is mostly composed of soft clay soil, which has poor engineering geological properties, such as high water content, large void, and high compressibility. According to existing research, Shanghai is dominated by clay. The natural pore ratio of the soil is 0.7-1.4, the natural water content is between 26.9-49.7%, and the compression modulus of the soil is between 2.56-6.43/Mpa(Es 0.1-0.2) [37,38]. There exists a certain linear relationship among the three, which reflects the characteristics of high-water content, large voids and high compressibility of the soils in Shanghai. With the development of surface urbanization and excessive development of urban underground space, the soil pressure in Shanghai is increasing. In addition, due to the characteristics of the soil, the bearing capacity of the area is low, and the overall subsidence and uneven subsidence are prone to occur [39,40]. Under the influence of geology and underlying structure, excessive subsidence occurred at the construction site of Shanghai Rail Transit Line 16 in 2011, which resulted in collapse. In 2012, the construction site of Metro Line 12 collapsed and have caused serious harm to people's property and safety. In Shanghai, the soft soil layer composed of soft clay soil is relatively thick. In terms of engineering construction, subway tunnels in Shanghai area generally pass through the soft soil layer, and viaducts and highways are all affected by the soft soil layer, thereby resulting in subsidence. Therefore, the influence of soft soil layer on regional engineering construction and daily maintenance cannot be ignored.
forces are mostly on a large time scale range, and the study span in this paper is relatively small. The subsidence in Shanghai is dominantly affected by groundwater and urban development, so the natural subsidence caused by tectonic forces in this paper is extremely small, and it is not the focus of this paper and can be ignored.
Overall, the geological structure of Shanghai area is mainly dominated by the soft soil layer, and the soft soil layer has a relatively large impact on engineering construction, which is also a main reason that the Shanghai area is prone to subsidence. In terms of deep bedrock, the bedrock borehole data used in this paper is about 300 m below ground. As Shanghai is a deltaic impact plain, the deep bedrock is dominated by mudstone and mud-intercalated sand, and the composition is dominated by feldspar, quartz and small gravels. The bedrock core of this type is harder, while the water content in the mud is still larger. The bedrock is hard plastic and the mud with high water content is still much more, which also reflects, to some extent, the poor geological structure of Shanghai and the importance of continuous monitoring of the subsidence phenomenon in Shanghai.
In terms of tectonic forces, there are no concentrated and large geological faults in Shanghai, and the regional geological structure is relatively stable. The effects of tectonic forces are mostly on a large time scale range, and the study span in this paper is relatively small. The subsidence in Shanghai is dominantly affected by groundwater and urban development, so the natural subsidence caused by tectonic forces in this paper is extremely small, and it is not the focus of this paper and can be ignored.
Overall, the geological structure of Shanghai area is mainly dominated by the soft soil layer, and the soft soil layer has a relatively large impact on engineering construction, which is also a main reason that the Shanghai area is prone to subsidence.

Correlation Analysis between Surface Subsidence Time Series and Rainfall
According to previous studies [41,42], urban rainfall and urban surface subsidence have a certain relationship. The greater the rainfall the more stable or somewhat uplift the land will be. The surface subsidence caused by groundwater exploitation is effectively alleviated because the increase in rainfall can replenish the groundwater content.
Shanghai is close to the ocean and is greatly affected by the ocean. Meanwhile, it is located in the subtropical monsoon climate. Thus, the precipitation times in Shanghai is longer, and the influence of rainfall on the change in groundwater level is greater.
In this paper, the monthly precipitation information of two meteorological stations QX1 and QX2, which are located in Baoshan and Longhua of Shanghai in 2019, respectively. They are selected from the National Meteorological Data Network. In order to better analyze the correlation between the time series of surface subsidence and rainfall, the subsidence data used in this section are the average values of the subsidence of PSs within 200 m around the weather station. In this experiment, a buffer zone of 200 m was established, with two weather stations as the center. The PSs within 200 m around the two weather stations were extracted separately, and the mean subsidence of these PSs were obtained as the basis for comparison.
The results of comparison between the two meteorological stations and the subsidence are shown in Figures 8 and 9. In the months with relatively more precipitation, the surface subsidence is significantly smaller or uplift, whereas in the months with relatively less precipitation, the surface subsidence is easy to occur significantly. thereby confirming the above inference.
The result shows that the seasonal variation of land subsidence is related to rainfall. We can combine the conclusion of linear relationship between subsidence and groundwater exploitation, which has been verified by previous studies. The over exploitation of groundwater will lead to surface subsidence, but the timely supplement of rainfall to groundwater will reduce the amount of subsidence effectively and can reduce the occurrence of excessive subsidence [43].

Surface Subsidence and Groundwater Level Change
The extraction/injection of subsurface resources (e.g., groundwater, natural gas, etc.) in the ground can directly affect the surface subsidence, while the surface subsidence can also reflect the subsurface resource content. Scholars have discovered that in the areas where the extraction of underground resources is more serious, the surface subsidence is more obvious, indicating that there exists a certain relationship between the exploitation of underground resources and the subsidence of the surface [44][45][46]. As a key city in China, Shanghai always has a serious problem of water exploitation. According to the "2018 Shanghai Water Resources Bulletin" and "2019 Shanghai Water Resources Bulletin" published by Shanghai Municipal Water Bureau, Shanghai's water resources in 2018 and 2019 are shown in Table 3. For nine consecutive years, the amount of artificial recharge in Shanghai is far greater than the amount of mining, which is also the fundamental reason for the stability of subsidence in most areas of Shanghai.
To further study the groundwater level changes and the relationship between the surface subsidence, this paper obtained the distribution in Songjiang district, Jinshan district, and Fengxian in four water level monitoring well data. Compared with the average subsidence in the surrounding 200-m area, the distribution of groundwater level monitoring wells is shown in Figure 3, and the comparison result is shown in Figure 10. The specific analysis is detailed as follows: monitoring wells A and B are located near Songjiang Erzhong, East Zhongshan Road, Songjiang District, and near Tinglin Vocational School, Songjin Road, Jinshan District; monitoring wells C and D are located near Caojing Chemical Plant, Huhang Road, Jinshan District, and near Zhujia Village, Qianqiao Town, Fengxian District.
The gray correlation analysis method was used to compare the P1-P4 accumulated The 2018 rainy season in the middle and lower reaches of the Yangtze River begins in June. The rainfall data of Shanghai area show that the rainfall in Shanghai has increased significantly since June. As shown in Figures 8 and 9, the subsidence amount of the subsidence feature points located in the QX1 region significantly decreased, whereas the subsidence feature points located in the QX2 region showed a small amplitude of uplift, thereby confirming the above inference.
The result shows that the seasonal variation of land subsidence is related to rainfall. We can combine the conclusion of linear relationship between subsidence and groundwater exploitation, which has been verified by previous studies. The over exploitation of groundwater will lead to surface subsidence, but the timely supplement of rainfall to groundwater will reduce the amount of subsidence effectively and can reduce the occurrence of excessive subsidence [43].

Surface Subsidence and Groundwater Level Change
The extraction/injection of subsurface resources (e.g., groundwater, natural gas, etc.) in the ground can directly affect the surface subsidence, while the surface subsidence can also reflect the subsurface resource content. Scholars have discovered that in the areas where the extraction of underground resources is more serious, the surface subsidence is more obvious, indicating that there exists a certain relationship between the exploitation of underground resources and the subsidence of the surface [44][45][46]. As a key city in China, Shanghai always has a serious problem of water exploitation. According to the  Table 3. For nine consecutive years, the amount of artificial recharge in Shanghai is far greater than the amount of mining, which is also the fundamental reason for the stability of subsidence in most areas of Shanghai. To further study the groundwater level changes and the relationship between the surface subsidence, this paper obtained the distribution in Songjiang district, Jinshan district, and Fengxian in four water level monitoring well data. Compared with the average subsidence in the surrounding 200-m area, the distribution of groundwater level monitoring wells is shown in Figure 3, and the comparison result is shown in Figure 10. The specific analysis is detailed as follows: monitoring wells A and B are located near Songjiang Erzhong, East Zhongshan Road, Songjiang District, and near Tinglin Vocational School, Songjin Road, Jinshan District; monitoring wells C and D are located near Caojing Chemical Plant, Huhang Road, Jinshan District, and near Zhujia Village, Qianqiao Town, Fengxian District.

Surface Subsidence and Underground Rail Transit Construction
Subway has now become an important means of transportation in people's daily lives. As of August 2020, Shanghai Metro has 17 operating lines (including maglev lines), 416 stations (including 2 maglev lines), and a total mileage of 705 km (including 29 km maglev lines). At the same time, four rail transit lines are under construction in Shanghai. Given the factors of geological soft soil in Shanghai, the construction of underground rail lines are greatly affected by geology, and abnormal subsidence occurs from time to time. Therefore, it is particularly necessary to monitor the subsidence of the rail transit lines that are under construction and their surrounding areas [54][55][56].
The construction of Shanghai Metro Line 15 was started in August 2016 and the whole line was opened in December 2019, with underground stations. In the time span of this experimental study, the line is in the period of accelerated construction and natural consolidation of the ground, so it is of great engineering significance to monitor its subsidence, and the study on the correlation between surface subsidence and underground The gray correlation analysis method was used to compare the P1-P4 accumulated subsidence curve with the groundwater level change curve. The grey correlation degree was obtained as: 0.65, 0.92, 0.59, 0.64, thereby showing high consistency [47,48]. When the groundwater level rises, the subsidence of the surface detection points tends to be stable or rises to some extent. Meanwhile, when the local water level decreases obviously, the subsidence of the surface detection points also occurs obviously. In July 2019, the rainy season gradually receded in the lower reaches of the Yangtze River, and the rainfall decreased significantly, thereby decreasing the groundwater level at the monitoring wells significantly. At the same time, P1-P4 experienced significant subsidence to varying degrees, thereby verifying the correlation between the groundwater level and the surface subsidence [49][50][51][52][53].

Surface Subsidence and Underground Rail Transit Construction
Subway has now become an important means of transportation in people's daily lives. As of August 2020, Shanghai Metro has 17 operating lines (including maglev lines), 416 stations (including 2 maglev lines), and a total mileage of 705 km (including 29 km maglev lines). At the same time, four rail transit lines are under construction in Shanghai. Given the factors of geological soft soil in Shanghai, the construction of underground rail lines are greatly affected by geology, and abnormal subsidence occurs from time to time. Therefore, it is particularly necessary to monitor the subsidence of the rail transit lines that are under construction and their surrounding areas [54][55][56].
The construction of Shanghai Metro Line 15 was started in August 2016 and the whole line was opened in December 2019, with underground stations. In the time span of this experimental study, the line is in the period of accelerated construction and natural consolidation of the ground, so it is of great engineering significance to monitor its subsidence, and the study on the correlation between surface subsidence and underground rail transit construction is the representative of Shanghai Metro Line 15.
In this section, Shanghai Metro Line 15 is selected as the research object, and the subsidence field is extracted from the line and the surrounding area within 600 m. In 2016, the construction of Shanghai Metro Line 15 began, and in December 2019, the whole line was successfully completed. The total length of Line 15 is 42.3 km with 30 stations, all of which are underground stations.
According to the monitoring of this line and its surrounding areas (the subsidence rate figure is shown in Figure 11), two obvious subsidence areas of A and B are found in Metro Line 15 (Respectively shown in a and b in Figure 11). Area a is located in Minhang District, in the south of Route 15, including Zizhu High-tech Zone Station to Jingxi Road Station. The accumulated subsidence of Area a is shown in Figure 12. Two obvious subsidence funnel areas are found on the subsidence trend line. By comparing Figure 11, the two places are the area of Yongde Road Station and the area between Yuanjiang Road Station and Shuangbai Road Station. The maximum accumulated subsidence of the two places is −40.2 mm and −37.8 mm respectively. Area b is located in Changning District, Shanghai, and includes Yaohong Road Station to Guilin Road Station. The subsidence of different degrees appeared in this area, and the subsidence near Wuzhong Road Station was the most obvious, with the maximum accumulated subsidence of −28.9 mm. In December 2019, the two underground construction lines of Metro Line 15 converged on the Guilin Station. The intersection of the two construction lines greatly increases the bearing capacity of the ground, which is the main reason for the obvious subsidence of the area [57,58]. After the completion of the subway, the natural consolidation subsidence becomes the main reason for the subsidence of Areas a and b [59,60].
In Zizhu high-tech area standing nearby, the subsidence phenomenon is more evident. The area of Minxing District of Shanghai Zizhu science park situates many universities and scientific research bases to gather hereat. Residents' living water all come from the exploitation of groundwater, and the consolidation subsidence caused by subway construction becomes the settling of the main reasons. [57,58]. After the completion of the subway, the natural consolidation subsidence becomes the main reason for the subsidence of Areas a and b [59,60].
In Zizhu high-tech area standing nearby, the subsidence phenomenon is more evident. The area of Minxing District of Shanghai Zizhu science park situates many universities and scientific research bases to gather hereat. Residents' living water all come from the exploitation of groundwater, and the consolidation subsidence caused by subway construction becomes the settling of the main reasons.   the main reason for the subsidence of Areas a and b [59,60].
In Zizhu high-tech area standing nearby, the subsidence phenomenon is more evident. The area of Minxing District of Shanghai Zizhu science park situates many universities and scientific research bases to gather hereat. Residents' living water all come from the exploitation of groundwater, and the consolidation subsidence caused by subway construction becomes the settling of the main reasons.

Influence of Industrialization Development on Land Surface Subsidence
In this section, feature region C in Figure 4 is selected for specific analysis. Area C is located in the south of Fengxian District, Shanghai. It is the gathering place of Fengxian District Industrial Area, where many large industrial plants are located. Figure 13 shows the deformation rate diagram in the characteristic area, and the three-dimensional accumulated subsidence diagram is shown in Figure 14. As shown in Figure 13, obvious subsidence occurs in area C, and the subsidence rate in many areas exceeds −20 mm/year. Among them, GY-a area is located near the intersection of Fengzhe Road and Zhenqian Road in Fengxian District with the maximum accumulated subsidence of −54 mm. A number of large industrial plants, such as Rubber & Plastic Products Co., Ltd., Thermal Energy Equipment Co., Ltd., are clustered in this area. By comparison with Figure 14, an obvious subsidence funnel in this area is observed. GY-b area is located in Xinghuo Development Industrial Zone, Fengxian District, where the subsidence phenomenon is more serious, with the maximum accumulated subsidence of −75 mm/year. near Dongfeng Automobile Special Parts Co., Ltd. SPIP area is Shanghai Paimilei Industrial Park, which has a number of heavy industrial bases. And Figure 15 shows the change of water level at the water level monitoring well D near the feature area in 2019. Compared with other areas of Shanghai, the groundwater level in this area is significantly lower because of the over-exploitation of groundwater by industrial production in this area. Overall, the overexploitation of groundwater caused by industrial development is the main cause of surface subsidence in the three industrial areas.    This area is the spark industrial area, gathered Jinfengyuan Paper Co., Ltd., Yadong Petrochemical and other large factories and heavy and light industrial bases.
The area with the most serious subsidence, as shown in the right part of Figure 13, is near Dongfeng Automobile Special Parts Co., Ltd. SPIP area is Shanghai Paimilei Industrial Park, which has a number of heavy industrial bases. And Figure 15 shows the change of water level at the water level monitoring well D near the feature area in 2019. Compared with other areas of Shanghai, the groundwater level in this area is significantly lower because of the over-exploitation of groundwater by industrial production in this area. Overall, the overexploitation of groundwater caused by industrial development is the main cause of surface subsidence in the three industrial areas.

Surface Subsidence and Urbanization Development
In recent years, with the rapid development of thermal infrared remote many countries have launched noctilucent remote sensing satellites.
Nighttime light remote sensing refers to the use of remote sensing techniqu ture the light distribution on the earth at night, which can effectively reflect th distribution of human activities, so it is often used for remote sensing inversion development [61][62][63].
In this paper, VIIRS noctilucent remote sensing satellite sensor was used t the noctilucent remote sensing data of the main urban area of Shanghai, and the d processing of the nighttime light imageries was carried out to obtain the urbaniza velopment changes in Shanghai from April 2018 to April 2020. Due to the poor effect of the monthly composite data from January 2018 to March 2018 published b this paper selected the monthly composite products from April 2018 and April 2 lished by VIIRS, as shown in Figures 16 and 17. After the difference processing, th in nighttime light imageries in Shanghai from April 2018 to April 2020 is obtaine paper (Figure 18). The overall trend of urbanization in Shanghai is developing

Surface Subsidence and Urbanization Development
In recent years, with the rapid development of thermal infrared remote sensing, many countries have launched noctilucent remote sensing satellites.
Nighttime light remote sensing refers to the use of remote sensing technique to capture the light distribution on the earth at night, which can effectively reflect the spatial distribution of human activities, so it is often used for remote sensing inversion of urban development [61][62][63].
In this paper, VIIRS noctilucent remote sensing satellite sensor was used to obtain the noctilucent remote sensing data of the main urban area of Shanghai, and the difference processing of the nighttime light imageries was carried out to obtain the urbanization development changes in Shanghai from April 2018 to April 2020. Due to the poor imaging effect of the monthly composite data from January 2018 to March 2018 published by VIIRS, this paper selected the monthly composite products from April 2018 and April 2020 published by VIIRS, as shown in Figures 16 and 17. After the difference processing, the change in nighttime light imageries in Shanghai from April 2018 to April 2020 is obtained in this paper ( Figure 18). The overall trend of urbanization in Shanghai is developing rapidly. Among them, it is located in the southeast of Shanghai Jinshan District and the southern cross of light brightness rising big scope, which reflect regional urbanization and the development of infrastructure construction. At the same time, compared with Figure 4, the subsidence of this area has a high correlation with the change in light brightness. This finding also proves that the process of urbanization and infrastructure construction will further increase the pressure on the surface, thereby resulting in uneven subsidence. The results show that urbanization and infrastructure construction are positively correlated with surface subsidence to a certain extent. velopment of infrastructure construction. At the same time, compared with Figu subsidence of this area has a high correlation with the change in light brightne finding also proves that the process of urbanization and infrastructure construct further increase the pressure on the surface, thereby resulting in uneven subsiden results show that urbanization and infrastructure construction are positively co with surface subsidence to a certain extent.   velopment of infrastructure construction. At the same time, compared with Figu subsidence of this area has a high correlation with the change in light brightn finding also proves that the process of urbanization and infrastructure construc further increase the pressure on the surface, thereby resulting in uneven subside results show that urbanization and infrastructure construction are positively co with surface subsidence to a certain extent.

Conclusions
In view of the lack of large-scale subsidence research, insufficient analysis o ence mechanism and low timeliness of experimental data in Shanghai. In this pap series InSAR technique was used to process the 36 Sentinel-1A SAR images from 2018 to March 2020. Combined with the changes in precipitation, groundwater l urbanization process in Shanghai, the spatial-temporal characteristics of surfac ence in Shanghai are analyzed in detail. The main conclusions are presented as f (1) The internal coincidence accuracy of PSInSAR technique is analyzed. Most standard deviations are less than 1 mm/year, which is relatively small. Ther accuracy and stability of PSInSAR technique for surface subsidence moni Shanghai area are relatively high. (2) The uneven subsidence in Shanghai is obvious. The surface subsidence in and east of Shanghai is stable, and the surface subsidence in some areas is However, the subsidence in Songjiang District in the east, Jinshan District in and Fengxian District in Shanghai are more serious, with large area subsid outward expansion trend. The subsidence near Xinghuocun in Fengxian D the most serious, with the maximum accumulated subsidence of −75 mm, related to the exploitation of groundwater caused by industrial water in thi (3) The study area is mainly a soft soil geological area, and the water content o is relatively high, which has a relatively large impact on engineering cons This reason is the main cause that Shanghai area is prone to subsidence. In seasonal changes in surface subsidence may be closely related to rainfall, an replenishment of groundwater by rainfall can reduce the subsidence effecti slow down the occurrence of excessive subsidence. The surface subsiden study area is highly correlated with the change in groundwater level. That the groundwater level rises, the surface will be uplifted or stable, and when t level drops, the surface will obviously sink. (4) Three subsidence areas were identified along the Shanghai Metro Line 15, w

Conclusions
In view of the lack of large-scale subsidence research, insufficient analysis of subsidence mechanism and low timeliness of experimental data in Shanghai. In this paper, time-series InSAR technique was used to process the 36 Sentinel-1A SAR images from January 2018 to March 2020. Combined with the changes in precipitation, groundwater level, and urbanization process in Shanghai, the spatial-temporal characteristics of surface subsidence in Shanghai are analyzed in detail. The main conclusions are presented as follows: (1) The internal coincidence accuracy of PSInSAR technique is analyzed. Most PS point standard deviations are less than 1 mm/year, which is relatively small. Therefore, the accuracy and stability of PSInSAR technique for surface subsidence monitoring in Shanghai area are relatively high. (2) The uneven subsidence in Shanghai is obvious. The surface subsidence in the north and east of Shanghai is stable, and the surface subsidence in some areas is uplifted. However, the subsidence in Songjiang District in the east, Jinshan District in the south and Fengxian District in Shanghai are more serious, with large area subsidence and outward expansion trend. The subsidence near Xinghuocun in Fengxian District is the most serious, with the maximum accumulated subsidence of −75 mm, which is related to the exploitation of groundwater caused by industrial water in this area. (3) The study area is mainly a soft soil geological area, and the water content of soft soil is relatively high, which has a relatively large impact on engineering construction. This reason is the main cause that Shanghai area is prone to subsidence. In addition, seasonal changes in surface subsidence may be closely related to rainfall, and timely replenishment of groundwater by rainfall can reduce the subsidence effectively and slow down the occurrence of excessive subsidence. The surface subsidence in the study area is highly correlated with the change in groundwater level. That is, when the groundwater level rises, the surface will be uplifted or stable, and when the water level drops, the surface will obviously sink. (4) Three subsidence areas were identified along the Shanghai Metro Line 15, which are the area of Yongde Road Station, the area between Yuanjiang Road Station and Shuangbai Road Station, and the area between Yaohong Road Station and Guilin Road Station. Among them, the subsidence of Yongde Road Station is the most serious with a maximum accumulated subsidence that reaches −40.2 mm. The main reason for the subsidence is the natural consolidation subsidence of the surface after the completion of the subway construction. Additionally, the industrial water in large industrial zones in Shanghai area increases the amount of underground water extraction, which is the main reason for the subsidence of industrial areas. Urbanization and infrastructure construction are also highly correlated with surface subsidence to a certain extent. Institutional Review Board Statement: Not applicable.