Changes of Chinese Coastal Regions Induced by Land Reclamation as Revealed through TanDEM-X DEM and InSAR Analyses

Chinese coastal topography has changed significantly over the last two decades due to human actions such as the development of extensive land reclamation projects. Newly-reclaimed lands typically have low elevations (<10 m) and often experience severe ground subsidence. These conditions, combined with the more frequent occurrence of extreme sea-level events amplified by global climate change, lead to an increased risk of flooding of coastal regions. This work focuses on twelve Chinese coastal areas that underwent significant changes from 2000 to 2015 in their environments, correlated to relevant land reclamation projects. First, the ground changes between 2000 and 2015 were roughly computed by comparing the TanDEM-X and the Shuttle Radar Topography Mission (SRTM) digital elevation models of the investigated areas. These results indicate that six of the analyzed coastal zones have reclaimed more than 200 km2 of new lands from 2000 to 2015, with five of them in northern China. Second, we focused specifically on the city of Shanghai, and characterized the risk of flood in this area. To this purpose, two independent sets of synthetic aperture radar (SAR) data collected at the X- and C-band through the COSMO-SkyMed (CSK) and the European Copernicus Sentinel-1 (S-1) sensors were exploited. We assumed that the still extreme seawater depth is chi-square distributed, and estimated the probability of waves overtopping the coast. We also evaluated the impact on the territory of potential extreme flood events by counting the number of very-coherent objects (at most anthropic, such as buildings and public infrastructures) that could be seriously affected by a flood. To forecast possible inundation patterns, we used the LISFLOOD-FP hydrodynamic model. Assuming that an extreme event destroyed a given sector of the coastline, we finally computed the extent of the flooded areas and quantified its impact in terms of coherent structures potentially damaged by the inundation. Experimental results showed that two coastline segments located in the southern districts of Shanghai, where the seawalls height is lower, had the highest probability of wave overtopping and the most significant density of coherent objects potentially subjected to severe flood impacts.


Introduction
Coastal zones are transitional regions connecting ocean and land systems containing highly concentrated human activities [1]. The rapid urbanization of coastal cities has yielded a scarcity of suitable space for development. Land reclamation is a common worldwide practice [2] as it can be an efficient way to mitigate land-use pressure. Around 52% of the world's megacities with populations greater than 10 million are located in coastal zones, and many of them have carried out some measure of coastal reclamation projects [3]. In China, coastal regions are essential for social-economic development. From 1991 to 2015, over 6418.9 km 2 of new lands were reclaimed from the sea [4], which have been utilized for agriculture (e.g., cultivation, aquaculture, and salt pan), industry (e.g., wharves and factories), and urban development (e.g., residences) [5]. Such reclamation projects have played a significant role in the socio-economic development of China [6]. Previous studies on the coastal evolution of China have mainly focused on calculating and analyzing the changes in coastline length, coastline types, reclamation areas, and reclamation intensity (i.e., the reclamation area per kilometer of coastline) [4,7,8]. Before the advent of remote sensing technologies [9,10], the length and location of coastlines were obtained from historical maps or field survey results, geomorphological maps [11], or historical documents that recorded land reclamation projects [12]. Currently, optical remote sensing images at multiple time scales are the main data source for coastal evolution research. Free and open access image archives such as Landsat, which have a long record and are geometrically and radiometrically consistent, are often employed [13][14][15][16]. Established methods of extracting coastlines include manual and automatic interpretation [17]. Manual interpretation is widely used, and its accuracy is determined by the spatial resolution of the images [18]. Automatic coastline delineation involves image segmentation [19], which detects the landsea boundary. Automatic interpretation methods include threshold segmentation [17,20], edge detection [21], object-oriented approaches [22], and active contour models [23]. Once the spatial and temporal changes of the coastline are determined, the area and rate of coastal erosion/accretion can also be obtained. However, the traditional methods applied for coastal evolution monitoring rely on identifying the land-sea boundary, extracting the coastline, and finally determining the evolution of coastlines. These conventional approaches usually cannot track with continuity the topographic changes of coastal environments and evaluate their consequent impacts on the environment and populations. Since the 1990s, the growing availability of radar satellites operating at microwave frequency bands has opened new possibilities to monitor the Earth's surface. In particular, the rapid development of synthetic aperture radar interferometry (InSAR) technology [24] has allowed for mapping the topography of the Earth on a global scale and generating digital elevation models (DEMs) of the terrain. InSAR technology relies on extracting the phase difference of couples of complex-valued SAR images of the terrain, acquired in different epochs and from slightly different orbital positions. In February 2000, the U.S. Shuttle Radar Topography Mission (SRTM) was launched and began acquiring data that, once adequately processed, has permitted the generation of a near-global DEM of the Earth [25], with a three-arc-second (and one-arc-second in some selected areas) resolution. More recently, in 2015, the TanDEM-X mission managed by DLR (Deutsches Zentrum für Luft-und Raumfahrt) [26] has yielded "updated" Earth DEMs. The comparison of InSAR-driven DEMs obtained from these different missions (SRTM and TanDEM-X), with a global scale coverage, allows for straightforward identification and study of the regions on Earth that have changed from 2000 to 2015. In particular, the DEMs comparison is useful to characterize the ground elevation of ocean-reclaimed lands. Although InSAR-driven DEMs have global coverage, their ground elevation accuracy is in the order of a few meters. Therefore, the comparison of the SRTM and TanDEM-X DEMs can only give a rough estimate of the occurred surface height changes. This is especially true for newly reclaimed lands located in low-elevation coastal zones (LECZ). On regional/continental scales, however, the results of such comparison could be sufficient to produce global measurements on the impact of land reclamation procedures on the coastal environment. Under climate change effects, potential by applying a multi-temporal SBAS algorithm [36] to two independent sets of SAR data collected from March 2018 to July 2021 by the COSMO-SkyMed (X-band) and the Sentinel-1 (C-band) sensors over complementary (ascending/descending) orbits. Subsequently, the vertical (subsidence) ground deformation rate and the relevant time-series were obtained by applying the multi-track InSAR combination method proposed in Pepe et al. [51]. Combined ground deformation time-series were produced on a georeferenced spatial grid with a sampling of 80 × 80 m. The impact of a potential flood for every selected segment was evaluated by counting the number of high-coherent objects (e.g., the point scatterers (PSs)), as obtained by analyzing the stack of available COSMO-SkyMed data at the single-look scale, located within the georeferenced grid.
The results of our investigation will further demonstrate the effectiveness of InSAR methodologies for the study of the ocean-reclaimed lands. Specifically, for the Shanghai case, we show that the coastline segments located in the south, where the height of the seawall is lower, are especially at risk in the case of extreme flood events because they are characterized by higher wave overtopping probabilities and by a high density of coherent objects that a severe inundation could impact.

Study Area
The coastal area of China comprises 14 provincial-level administrative regions between the latitudes of 18 • N to 43 • N and longitudes 104 • E to 125 • E. These regions only account for 13% of China's total territorial area but contain 43.5% of the population [12]. Figures from the National Bureau of Statistics of China [61] show that the gross domestic product (GDP) of the coastal provinces was 41.6 trillion yuan (approx. USD 6.2 trillion) in 2015, 56% of the total GDP of China. According to the China Marine Statistical Yearbook (2016) [62], along the 18,000 km of continental coastlines, there are 57 prefecture-level cities located along the 18,000 km coastline, with 23 cities having a population greater than five million (see Figure 1), which were, for our purposes, classified as large cities [63]. Among them, we selected twelve cities with high-intensity coastal reclamation projects as the primary subject of this paper that will be analyzed in detail.
Furthermore, we selected the Yangtze River Delta and the Shanghai region as case study areas. The YRD, an alluvial plain formed by sediment deposition from the Yangtze River [64], is one of the most crucial river delta regions of China and has the highest population density. In recent decades, intensive coastal reclamation has been conducted [7,64], yielding significant coastal changes. Shanghai is located in the Yangtze Estuary at the midpoint of the YRD and is one of the most densely populated and economically developed megacities globally. Most of Shanghai is at a low elevation (about 3-4 m above sea level) and has undergone high-intensity coastal reclamation (approx. 1192.5 km 2 from 1960 to 2015 [65]) as well as constructing flood-control seawalls [66].

Digital Elevation Models
Open-source medium-resolution global digital elevation model (DEM) datasets [67] have widely been used to extract geomorphological features, especially in regional-scale studies. The U.S. Shuttle Radar Topography Mission has provided near-global DEM products since 2003 [25]. In our study, NASA's SRTM-1 V3 DEM product of one arc-second resolution (approx. 30 m at the equator) was utilized (https://earthexplorer.usgs.gov/ (accessed on 3 July 2021)). The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) global DEM (Version 2) covers 99% of the Earth's surface at one arcsecond resolution, which is also available from the USGS Earth Explorer site. The Advanced Land Observing Satellite (ALOS) has multiple instruments including the Panchromatic Remote-sensing Instrument for Stereo Mapping (PRISM), which is designed to map digital elevation. The ALOS WORLD 3D (AW3D) DEM includes 5 m and 30 m resolution; the former is distributed commercially, and the latter is freely available from the Japan Aerospace Exploration Agency (JAXA) website. Finally, a newly released (September 2018) global DEM, which has a three-arc-second resolution (TanDEM-X 90 m; TDX90), has been developed by DLR (Deutsches Zentrum für Luft-und Raumfahrt) [26].
TDX90 DEM was derived from the standard TanDEM-X (TDX) DEM product that has a 0.4 arc second posting (approx. 12 m at the equator) and vertical accuracy of 2 m at a 90% confidence level [68]. TDX90 has a lower vertical accuracy than standard TDX, but is free of charge and more easily available. Globally, root mean square error (RMSE) of the elevation difference between the ground-based checkpoints and DEM extracted points for SRTM, ASTER, and AW3D30 DEMs were 6 m [69], 8.68 m [70], and 4.4 m [71], respectively. However, for the TDX90 DEM, there is not an overall global assessment of its vertical accuracy. Hawker [72] found that compared with SRTM, ASTER, and AW3D30 DEMs, TDX90 has the best vertical accuracy in 32 low-slope floodplain sites. More details about each of these global DEMs are shown in Table 1.

Digital Elevation Models
Open-source medium-resolution global digital elevation model (DEM) datasets [67] have widely been used to extract geomorphological features, especially in regional-scale studies. The U.S. Shuttle Radar Topography Mission has provided near-global DEM products since 2003 [25]. In our study, NASA's SRTM-1 V3 DEM product of one arc-second resolution (approx. 30 m at the equator) was utilized (https://earthexplorer.usgs.gov/ (accessed on 3 July 2021)). The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) global DEM (Version 2) covers 99% of the Earth's surface at one arcsecond resolution, which is also available from the USGS Earth Explorer site. The Advanced Land Observing Satellite (ALOS) has multiple instruments including the Panchromatic Remote-sensing Instrument for Stereo Mapping (PRISM), which is designed to map digital elevation. The ALOS WORLD 3D (AW3D) DEM includes 5 m and 30 m resolution; the former is distributed commercially, and the latter is freely available from the Japan Aerospace Exploration Agency (JAXA) website. Finally, a newly released (September 2018) global DEM, which has a three-arc-second resolution (TanDEM-X 90 m; TDX90), has been developed by DLR (Deutsches Zentrum für Luft-und Raumfahrt) [26].
TDX90 DEM was derived from the standard TanDEM-X (TDX) DEM product that has a 0.4 arc second posting (approx. 12 m at the equator) and vertical accuracy of 2 m at a 90% confidence level [68]. TDX90 has a lower vertical accuracy than standard TDX, but is free of charge and more easily available. Globally, root mean square error (RMSE) of the elevation difference between the ground-based checkpoints and DEM extracted points for SRTM, ASTER, and AW3D30 DEMs were 6 m [69], 8.68 m [70], and 4.4 m [71], respectively. However, for the TDX90 DEM, there is not an overall global assessment of its vertical accuracy. Hawker [72] found that compared with SRTM, ASTER, and AW3D30 DEMs, TDX90 has the best vertical accuracy in 32 low-slope floodplain sites. More details about each of these global DEMs are shown in Table 1. The vertical datum for TDX90 DEM is the World Geodetic System 1984 (WGS84) ellipsoid [68], whereas other DEMs are referenced vertically to the Earth Gravitational Model 1996 (EGM96) geoid [72]. Hence, pre-processing was conducted to convert the TDX90 DEM to the EGM96 height datum for data consistency. Figure 2 shows the final TDX90 DEM used in our study, with the location of major twelve coastal cities with highintensity land reclamation indicated by   [72]. Hence, pre-processing was conducted to convert the TDX90 DEM to the EGM96 height datum for data consistency. Figure 2 shows the final TDX90 DEM used in our study, with the location of major twelve coastal cities with highintensity land reclamation indicated by

Landsat Images
Landsat images acquired from 1983 to 2019 at 3-year intervals (78 total images) were selected to interpret the coastline of the YRD (Table 2). Six Landsat scenes with different paths/rows (Table 2) covered the entire YRD. The selected images had a criterion of zero or very little cloud coverage (<1%) and zero pixels near the coastline blocked by clouds. We performed geometric corrections, radiometric calibrations, atmospheric corrections, and image mosaics to obtain a series of images covering the Yangtze Delta. The coastline in 2019 was vectorized first, and the changed regions of the coastline were modified to acquire the coastline in 2016, which was then used as a reference for the coastline in 2013. We iterated this method in discrete steps back to 1983.

Landsat Images
Landsat images acquired from 1983 to 2019 at 3-year intervals (78 total images) were selected to interpret the coastline of the YRD (Table 2). Six Landsat scenes with different paths/rows (Table 2) covered the entire YRD. The selected images had a criterion of zero or very little cloud coverage (<1%) and zero pixels near the coastline blocked by clouds. We performed geometric corrections, radiometric calibrations, atmospheric corrections, and image mosaics to obtain a series of images covering the Yangtze Delta. The coastline in 2019 was vectorized first, and the changed regions of the coastline were modified to acquire the coastline in 2016, which was then used as a reference for the coastline in 2013. We iterated this method in discrete steps back to 1983.  1986,1989,1992,1995,1998,2001,2004,2007 23 August 2021 with ascending passes, have a revisit time of 12 days. The CSK and S-1A data were acquired in the StripMap (SM) and interferometric wide-swath (IW) modes, respectively.

Methodology
This section briefly introduces the main techniques and the methodologies applied to analyze the flood risk over the selected ocean-reclaimed lands.

Interferometric Synthetic Aperture Radar Techniques
Interferometric synthetic aperture radar (InSAR) techniques [24,73] rely on the extraction and analysis of the phase difference between couples of complex-valued SAR images collected over the same area on Earth and from slightly different orbital positions. They have been applied to generate elevation maps of the Earth's surface and produce DEMs of the terrain [74,75] (see Section 2.2.1). The differential SAR interferometry technology [24,76,77] is an evolution of InSAR, which was developed to obtain maps of the radar line-of-sight-projected ground displacements. Historically applied to analyze single deformation episodes such as the ground ruptures due to large earthquakes [78,79], DIn-SAR has then evolved to produce ground deformation time-series using multi-temporal InSAR methods [34,38,80,81]. Among these methods, a crucial role has been played by the SBAS technique [36]. It processes a group of short baseline interferograms characterized by low phase decorrelation effects [82]. Well-processed SAR pixels characterized by temporal coherence values [83] higher than a given threshold (in our work, this threshold value was 0.5) were identified and used for the subsequent analyses. To increase the density of well-processed coherent SAR pixels, we also applied a region-growing strategy [84]. SAR images collected from ascending and descending passes allow one to discriminate from the LOS-projected ground deformation time-series the vertical (subsidence) ground deformations from the lateral East-West movements by applying multi-track satellite combination methods such as the minimum acceleration (MinA) combination technique [51]. The ground deformation measurement accuracy in single interferograms depends on the used operational wavelength λ. If ∆φ 1,2 is the (unwrapped) phase difference between the SAR images collected at times t 1 and t 2 , the standard deviation σ de f of the ground displacement occurring between the observed epochs is equal to σ de f = λσ ∆φ 1,2 /4π where σ ∆φ 1,2 is the standard deviation of the unwrapped phases that depends on both the decorrelation noise and the phase unwrapping errors [82]. With X-band COSMO-SkyMed and Sentinel-1 data, the operational wavelength is about 3.1 cm and 5.4 cm. Accordingly, the measurement accuracy is on the order of a few millimeters, depending on the amount of decorrelation noise and phase unwrapping artifacts present. Quantitative analyses [54,[85][86][87], based on the comparison of SBAS-driven results with GPS/leveling data, have shown that the accuracy of the SBAS-based LOS-projected ground displacement and mean deformation rate is to the order of 3-5 mm and 1 mm/year, respectively. The accuracy of the 2-D (Up/Down and East-West) components of the ground deformations, as obtained by combining SAR data acquired from ascending and descending orbits of a sensor with approximately the same side illumination angle, can roughly be derived as a function of the accuracy of LOS-projected products σ LOS by σ U p−Down = 2σ LOS /cos ϑ and σ East−West = 2σ LOS /sin ϑ where ϑ is the sensor-to-target wave (local) incidence angle. Accordingly, Up/Down ground deformation products are more accurate than East-West products, whereas there is minimal sensitivity to ground movements along the North-South direction. Interested readers can find details on the SBAS and MinA techniques and additional information on the ground deformation measurements accuracy by Casu et al. and Pepe et al. [86,88].

Coastal Reclamation Area Detection Method
Coastal reclaimed lands were retrieved on the regional/continental scale by comparing the elevation height recorded in the oldest SRTM with the newest TDX90 global DEM (see Table 1). The land-sea border can be validated, since it is easy to distinguish between the land and sea areas of different elevations. We conducted superposition analyses of DEMs to detect coastal reclamation areas, called the DEM-based method. The specific steps performed were (i) produce the vector boundaries of SRTM and TDX90 DEM, denoted by B S and B T ; (ii) subtract B S from B T and obtain the vector boundary of reclamation land, denoted by B C ; (iii) using vector B C as a mask, clip the TDX90 DEM within it, and obtain the DEM of the newly-formed land, denoted by DEM C ; and (iv) clip the DEM C with the administrative boundary of coastal cities, and obtain the DEM for each prefecture-level city.

Quantitative Analysis Method of Coastline Change
We used the coastline interpretation method discussed in Section 2.2.2 to vectorize thirteen coastlines in different years of the YRD. We then employed the Digital Shoreline Analysis System (DSAS, version 5.0) [89] to compute and quantitatively analyze the rate of coastline changes. The coastline in 1983 was used as a baseline, and a series of transects perpendicular to the baseline were generated along the shoreline with equal spacing of 500 m. Note that the selected Chinese coastal regions have a smooth topography; hence, the considered spatial spacing is adequate to perform the analyses at the regional/continental scale. In the analysis of land area change, only the oldest (in 1983) and newest (in 2019) coastlines were considered, and the region between them was regarded as a land accretion/erosion area. We here define this coastline change analysis method the coastline-based method.

Coastline Segmentation and Coastal Flood Modeling
Under the combined effect of wave, subsidence, and the flood-control seawalls (with different heights), certain segments of seawalls are more vulnerable to burst or more likely to be overtopped by seawater. However, there is a low probability that all the vulnerable seawalls would simultaneously lose efficacy. The possibility of a single seawall segment failing is more realistic. Thus, we divided the seawall system of Shanghai into 40 segments, with the length of every section being about 2 km. We assumed that each section could fail individually, so there were 40 independent scenarios. For the flood modeling, we utilized the two-dimensional hydrodynamic model LISFLOOD-FP, a rasterbased flood inundation model, to simulate the extent and depth of inundation of the coastal floodplain [57,90]. LISFLOOD-FP has been widely used to simulate the evolution of coastal flooding and performed well in many coastal inundation studies [91][92][93][94]. We used the model's Acceleration Solver, which employs simplified shallow water equations by neglecting the convective acceleration term [58]. Three files need to be prepared to run the model. One is the elevation file (raster DEM). The second is a file with the information of boundary condition types used to simulate the initial state of flood along the coast. It contains coordinates of water points, water depth, and the form of water flow. Here, the fixed free surface elevation was selected, assuming seawater always maintains a specific water depth during inland flooding. The third is a parameter file, which states the solver, sets the storage path and name, etc. We used the Global Tide and Surge Reanalysis (GTSR) datasets [59] as the water depth. It simulates extreme sea levels with a 100-year return period and takes the effects of tides and surges into account. The sea levels in GTSR have been corrected to the same datum as SRTM (i.e., EGM96) [95].

Flood Risk Maps Evaluation
Here, we describe the strategy adopted to generate risk maps in coastal areas due to potential extreme weather flood events. We considered the presence of a network of seawalls along the coastline, and we applied the LISFLOOD-FP hydrodynamic model to estimate the extent of the regions inundated after a flood event. The model describes the dynamic propagation of waves on floodplains using continuity and momentum equations, discretized over a grid mesh [96,97]. During dynamic simulations, the wave propagation was controlled by friction (Manning's number) and gravity forces [94,98]. The computation of a risk map requires the knowledge of the probability that a hazard could occur and the evaluation of the impact of a disaster event on the population when it happens. To this aim, we studied the characteristics of the seawall networks, with a specific focus on those deployed along the Shanghai coastline. We remark that the design standard for seawall heights is based on the combination of extreme wave heights and extreme still sea levels. This standard is not universal for different coastal cities in China because of their various economies and natural environments [99]. Due to Shanghai's important place in the national economy, the protection standard of the seawalls in Shanghai is relatively high, and about 78% of these projects are at a standard that can withstand a 100-year storm tide plus waves caused by force-12 winds [100]. The dataset, which contains the crest heights of seawalls in 2010 along coastal Shanghai, shows that the minimal seawall height is around 7.6 m. The GTSR datasets show that the still extreme sea level of Shanghai is about 3.71 m and is much lower than the minimal seawall height. Thus, there is a small probability of overtopping. To estimate this probability, we considered that the extreme sea level is not a constant number and could be perturbed higher or lower under the influence of climatic factors (e.g., wind speed). Specifically, to calculate the probability of overtopping, we considered that the extreme still water level (i.e., from GTSR datasets) had a Pearson Type III distribution [99,101] that under proper conditions may be approximated with a chi-square distribution of an average known value. It is also worth remarking that another commonly used distribution for extreme water levels in hydraulics is also the generalized extreme value distribution [102][103][104].
The probability chi-square density function (PDF) is: where k (integer number) is the degree of freedom; x is the water level; and Γ is the gamma function. The expectation value of the chi-distribution equals k, which in the case of the Shanghai area is equal to 3.71 and has been rounded to 4. Once the PDF is determined, the probability of overtopping can be calculated as: To obtain some simple predictions on future risk conditions due to extreme flood events, we also corrected the seawall height for ten years from now, considering, as a worst-case, that the observed subsidence rate of the seawalls will remain constant. Previous works have demonstrated that the subsidence of ocean-reclaimed lands is not linear but has a time-decay behavior [105]. Nonetheless, for the scope of the present investigation, we focused exclusively on the ground (subsidence) deformation rate to identify the critical sectors of the Shanghai coast to provide some simple predictions on future evolution and impacts of extreme floods on the Shanghai coastal urban texture.
Accordingly, for pixels located in correspondence with the seawall location, in Equation (2), we used the corrected (possible) future height given by: where t 1 is the time when seawall leveling was measured and t 2 is a time in the future. V is the Up-Down mean velocity obtained from CSK/S-1A datasets, which is a negative value when subsidence occurs. H t1 is the seawall height at t 1 , and H t2 is the corrected seawall height at t 2 . For each segment of the seawalls, we also evaluated the impact of a potential flood by counting the number of high-coherent point scatterers (PSs) that fall within the flooded regions. The PSs points reflect, at most, the distribution of man-made structures such as buildings, roads, seawalls, and ports. We analyzed a set of CSK differential SAR interferograms computed at the single-look (full) spatial resolution. We applied an azimuth split-spectrum algorithm that allows for straightforward identification of those scatterers that preserve their coherence at different non-overlapped frequency bands. This method consists of splitting the azimuth spectra of couples of SAR images into two non-overlapping slices (the upper and lower Doppler bands) from which two distinctive interferograms (from the upper and lower bands) can be generated. The phase difference of the upper and lower interferograms leads to a figure on the phase stability of the scatterers. It permits the identification of objects that are coherent at both low and up frequencies. Specifically, the following coherence-like estimator: is calculated, where j = √ −1, and φ high , φ low are the relevant k-th DInSAR interferograms that consider the high and the low azimuthal (Doppler) frequencies, respectively. The detected high-coherent points (i.e., those related to PS objects) are identified by simply considering those pixels with a coherence higher than a given threshold, which in our case is equal to 0.35. The map of PS pixels can then be geo-referenced for the subsequent analyses. Details of the adopted azimuth split-spectrum method can be found in Falabella et al. [106].
Thus, the density of detected high-coherent points (denoted by D P ) in the inundation areas are used as an indication of the flood impacts on the coastal communities: where Num 1 is the number of high-coherent points in the given cell of the geocoded spatial grid, and Num 2 is the extent of every patch of the 90 × 90 m geocoded grid used to implement the risk flood analyses in coastal regions.

Landforms in Coastal Reclamation Areas of Twelve Large Coastal Cities of China
The SRTM DEM gives the location and elevation of land in 2000, and the TDX90 DEM provides the same information for 2015. Comparing both DEMs, we can estimate the coastal reclamation areas from 2000 to 2015. According to the area and shape of reclamation land (Figures 3 and 4), these twelve cities can be categorized into four types: (i) Cities that have expanded their coastlines significantly with a straight outline including Tangshan (b) and Tianjin (c), which have built a port and a harbor, respectively; (ii) Cities that have had relatively uniform reclamation projects including Dalian (a), Yancheng (e), and Nantong (f). In these three cities, many marine aquaculture ponds have been built in the reclaimed regions; (iii) Economically developed cities, with large reclamation areas and straight artificial seawalls built on the periphery including Shanghai (g) and Ningbo (h); and (iv) Cities where coastal reclamation projects have been small and intermittent including Zhanjiang (l), Qingdao (d), Wenzhou (i), Fuzhou (j), and Quanzhou (k). reclaimed regions; (iii) Economically developed cities, with large reclama straight artificial seawalls built on the periphery including Shanghai (g) a and (iv) Cities where coastal reclamation projects have been small and i cluding Zhanjiang (l), Qingdao (d), Wenzhou (i), Fuzhou (j), and Quanzho   From 2000 to 2015, the total reclamation area of the twelve large coasta ered here was 3167.60 km 2 (Table 3), accounting for 2.34% of the entire area. Specifically, the coastal reclamation area of the cities in the Bohai Bay in the YRD (e-h), and the cities in southern China (i-l) were 1290.28 km 2 , 13 409.25 km 2 , respectively. This means their proportions in the total reclama 40.73%, 43.44%, and 12.92%, respectively. The cities of the YRD had the la tion area, the largest inundation area of the whole city (11,608.18 km 2 , 69 largest reclaimed land (1136.94 km 2 , 59.97%). For the proportion of the re in the city, the maximum and minimum were 4.96% (h: Ningbo) and 0.62% From 2000 to 2015, the total reclamation area of the twelve large coastal cities considered here was 3167.60 km 2 (Table 3), accounting for 2.34% of the entire administrative area. Specifically, the coastal reclamation area of the cities in the Bohai Bay (a-c), the cities in the YRD (e-h), and the cities in southern China (i-l) were 1290.28 km 2 , 1375.97 km 2 , and 409.25 km 2 , respectively. This means their proportions in the total reclamation area were 40.73%, 43.44%, and 12.92%, respectively. The cities of the YRD had the largest reclamation area, the largest inundation area of the whole city (11,608.18 km 2 , 69.19%), and the largest reclaimed land (1136.94 km 2 , 59.97%). For the proportion of the reclamation area in the city, the maximum and minimum were 4.96% (h: Ningbo) and 0.62% (l: Zhanjiang). Although the administrative areas of the different cities are diverse, there was a positive correlation between reclamation areas and their proportions in the administrative regions (correlation coefficient was 0.91). The average elevations in the newly reclaimed land ranged from 1.47 m (l: Zhanjiang) to 7.15 m (j: Fuzhou), and as shown in Figures 3 and 4, there was almost no grid with an elevation of more than 10 m. Thus, the large majority of newly reclaimed land considered here was less than 10 m and belonged to the LECZ [107]. Table 3. Important statistics for coastal cities of China that have experienced significant land reclamation including population, reclamation area, average elevation, and projected inundation areas in total city area and reclaimed land. The last column shows the proportion of flood areas in reclaimed land to total flooded areas of each city.

Inundation Mapping of Twelve Coastal Cities
Overall, 12.42% of the total land in all cities considered will be flooded under extreme sea level, and this impact will be more intense in newly reclaimed land, of which 59.85% (1895.82 km 2 ) will be flooded Table 3. In particular, the inundation proportion of reclaimed land was highest in Yancheng (89.17%), followed by Wenzhou (87.72%), Nantong (86.92%), and Shanghai (80.45%). Yancheng borders Nantong to the south. The average elevation of these two cities is low (2.74 m and 4.05 m, respectively); in these locations, as shown in Table 3 and Figure 5 the inundation proportion of the entire city area is also very high (55.01% and 19.13%, respectively). In addition, Yancheng-Nantong includes a large area of aquaculture ponds and silty-sandy beach within reclaimed lands, and there are extensive deep-water (>4 m) flooded areas (Figure 6e,f). Conversely, for Shanghai and Wenzhou, the inundation proportion of the entire city is only about 6%, far less than reclaimed regions (approx. 85%). Table 3 and Figure 5 show that Ningbo and Quanzhou have a similar inundation proportion for reclaimed regions (approx. 70%). Ningbo's floodplain is concentrated in the north (Figure 7h), and its total inundation area (333.13 km 2 ) ranked second in our list. However, Quanzhou's inundation area (70.64 km 2 ) is relatively small. For Tangshan-Tianjin, most coastal wharves will be submerged (Figure 6b,c), and about 12% of the total city area will be flooded. As shown in Table 3, for Dalian and Qingdao, flood areas, flood percentage of the whole city area, and flood percentage of the reclaimed area were relatively low due to high average elevations (81.9 m and 55.1 m, respectively). The last column in Table 3 shows the effects of elevation changes in newly reclaimed land on flood inundation. Based on this index, we categorized these cities into four groups. For Dalain, >90% of the flood inundation occurred in newly reclaimed lands; followed by Quanzhou, Shanghai, and Ningbo (>35% and <50%); then, Qingdao, Nantong, Wenzhou, Fuzhou (>15% and <25%); and Tianjin, Tangshan, Zhanjiang, and Yancheng, where <10% of the flood inundation occurred in the newly reclaimed land.  Table 3. City names are displayed in the abbreviated format: Nantong is NT, Yancheng is YC, Wenzhou is WZ, Shanghai is SH, Ningbo is NB, Quanzhou is QZ, Zhanjiang is ZJ, Fuzhou is FZ, Dalian is DL, Qingdao is QD, Tangshan is TS, and Tianjin is TJ.   Table 3. City names are displayed in the abbreviated format: Nantong is NT, Yancheng is YC, Wenzhou is WZ, Shanghai is SH, Ningbo is NB, Quanzhou is QZ, Zhanjiang is ZJ, Fuzhou is FZ, Dalian is DL, Qingdao is QD, Tangshan is TS, and Tianjin is TJ.  Table 3. City names are displayed in the abbreviated format: Nantong is NT, Yancheng is YC, Wenzhou is WZ, Shanghai is SH, Ningbo is NB, Quanzhou is QZ, Zhanjiang is ZJ, Fuzhou is FZ, Dalian is DL, Qingdao is QD, Tangshan is TS, and Tianjin is TJ.

Coastal Evolution in the Yangtze River Delta
We compared the coastlines of the YRD in 1983 and 2019 using Landsat images a described above (Table 2) to obtain the land expansion/erosion over the past 36 years. W divided the Yangtze River Delta into six segments from north to south ( Figure 8) and quantitatively analyzed them for temporal and spatial change rates. Overall, most seg ments in the Yangtze Delta showed a seaward expansion (i.e., accretion). The average rat of shoreline expansion was 60 m/yr. Only a small section within segment one from th Guanhe River Estuary to Biandan Port in northern Jiangsu Province (labelled by (a) i Figure 8) and a section in southern Shanghai within segment 3 (labelled by (b) in Figur 8) experienced slight erosion (marked by red in Figure 8). The rate of coastline expansio in segment 2 (covering Yancheng and Nantong) and segment 4 (surrounding Ningbo) i much greater than in other locations. The time-series analysis shows that the rates of thes two segments fluctuate considerably, and the change rate was most significant from 200 to 2016, which indicates that most reclamation projects for these cities occurred durin this time. This is consistent with the results in Figures 3 and 4. Segment 3 covers the main

Coastal Evolution in the Yangtze River Delta
We compared the coastlines of the YRD in 1983 and 2019 using Landsat images as described above (Table 2) to obtain the land expansion/erosion over the past 36 years. We divided the Yangtze River Delta into six segments from north to south ( Figure 8) and quantitatively analyzed them for temporal and spatial change rates. Overall, most segments in the Yangtze Delta showed a seaward expansion (i.e., accretion). The average rate of shoreline expansion was 60 m/yr. Only a small section within segment one from the Guanhe River Estuary to Biandan Port in northern Jiangsu Province (labelled by (a) in Figure 8) and a section in southern Shanghai within segment 3 (labelled by (b) in Figure 8) experienced slight erosion (marked by red in Figure 8). The rate of coastline expansion in segment 2 (covering Yancheng and Nantong) and segment 4 (surrounding Ningbo) is much greater than in other locations. The time-series analysis shows that the rates of these two segments fluctuate considerably, and the change rate was most significant from 2001 to 2016, which indicates that most reclamation projects for these cities occurred during this time. This is consistent with the results in Figures 3 and 4. Segment 3 covers the mainland of Shanghai, and segment 6 covers three islands in the Yangtze Estuary (Chongming, Hengsha, and Changxing Island) under Shanghai's administration. Segments 3 and 6 jointly express the total reclamation situation of Shanghai. The overall expansion speed of Shanghai's shoreline was between 29.4 ± 7.5 and 44.8 ± 10.5 m/yr, with the expansion rate of the islands being larger than the mainland. Segment 1 is in the northernmost region of the Yangtze River Delta and has the smallest shoreline expansion rate. City-wise, Ningbo (segment 4) had the highest-intensity reclamation, followed by Yancheng and Nantong (segment 2), while Shanghai (segments 3 and 6) was relatively low. Hengsha, and Changxing Island) under Shanghai's administration. Segments 3 and 6 jointly express the total reclamation situation of Shanghai. The overall expansion speed of Shanghai's shoreline was between 29.4 ± 7.5 and 44.8 ± 10.5 m/yr, with the expansion rate of the islands being larger than the mainland. Segment 1 is in the northernmost region of the Yangtze River Delta and has the smallest shoreline expansion rate. City-wise, Ningbo (segment 4) had the highest-intensity reclamation, followed by Yancheng and Nantong (segment 2), while Shanghai (segments 3 and 6) was relatively low.  Figure 4 and Table 3 show that four southern cities (Wenzhou, Fuzhou, Quanzhou, and Zhanjiang) and Qingdao have carried out low-intensity and spatially intermittent reclamation projects. The coastal reclamation activities in Dalian, Yancheng, and Nantong were mainly for aquaculture, salt pan, and cultivation [2]. For Tangshan and Tianjin, the sea enclosure was primarily for constructing new ports and harbors. Both Shanghai and Ningbo have utilized a large reclamation land area for human habitation. Especially in Shanghai, a township with a population of 830,000 named Lingang New City was built here. Shanghai is unique for its continuous artificial seawalls, low-elevation topography, Figure 8. Area variation and rate of coastline change in the Yangtze Delta. The delta's shoreline was divided into six segments based on coastline change rates and geographical characteristics. The separation of segments 1 and 2 was due to the reclamation intensity. Segment 2 to segment 5 were divided by the Yangtze Estuary, the Qiantang River Estuary, and the Yongjiang River Estuary from north to south. Segment 6 consists of three islands administered by Shanghai. Sections (a,b) indicate locations where shoreline erosion has occurred. Figure 4 and Table 3 show that four southern cities (Wenzhou, Fuzhou, Quanzhou, and Zhanjiang) and Qingdao have carried out low-intensity and spatially intermittent reclamation projects. The coastal reclamation activities in Dalian, Yancheng, and Nantong were mainly for aquaculture, salt pan, and cultivation [2]. For Tangshan and Tianjin, the sea enclosure was primarily for constructing new ports and harbors. Both Shanghai and Ningbo have utilized a large reclamation land area for human habitation. Especially in Shanghai, a township with a population of 830,000 named Lingang New City was built here. Shanghai is unique for its continuous artificial seawalls, low-elevation topography, and large population. Additionally, subsidence has been widely confirmed in the reclamation area of Shanghai [45,46,108] from multi-platform InSAR measurements. In this work, we also show DInSAR-driven ground displacement maps that are useful for insights into the future evolution of coastline subsidence. The analyses presented here have been carried out at multi-look scales using different sets of SAR images, as detailed in the following Section 3.4.1.

Analysis of High-Coherent Objects in the Shanghai Coastal Area
Using the SBAS approach, we obtained separate deformation time series of two platforms in the same periods. For both, the LOS-projected deformation time series were calibrated to a common reference pixel (identified by the red star in Figure 9), whose deformation was nearly zero. Then, the LOS-projected deformation measurements were preliminarily geocoded to a reference (oversampled) latitude/longitude grid with a grid size of about 3 × 3 m. Finally, on the common geocoded area between the CSK and S-1A SAR data footprints, we applied the MinA combination approach [51] to decompose the LOS deformation time series into Up-Down and East-West directions. The mean deformation velocity map in the Up-Down direction is shown in Figure 9c. We can conclude that the Up-Down subsidence mainly occurred in the eastern region, especially in the triangular area near Pudong Airport. The areas that were subjected more to subsidence were the object of our subsequent investigations. We focused on the coastline sectors characterized by significant ground subsidence rates, and analyzed the deformations associated with the seawall instabilities.
To estimate the impact of a potential flood on the coastal regions, we also analyzed high-resolution CSK images using the azimuth split-spectrum method [106] and the strategy outlined in Section 2. The spatial distribution of the high-coherent objects in the investigated area is shown in Figure 10a.
As shown in Figure 10a, along the coast, high-coherent objects are significantly concentrated in the south-central part of Shanghai, reflected by high-density grids. In the southwest of Dishui Lake and the east-central part of Shanghai, the density of highcoherent objects is intermediate. In contrast, other regions have relatively low densities of coherent targets.

The Inundation Map for Different Segmentation Scenarios
Following Section 2, we split the coastline seawalls into forty non-overlapping segments and estimated the probability of wave overtopping and flood impacts in the case of an extreme weather event under the failure of each individual segment by considering the relevant seawall height and Equations (1) and (2). The impact is calculated by applying the LISFLOOD-FP hydrodynamic model [57,58], assuming that, for every selected sector, incoming waves overpass the barriers and the seawall sector has a severe failure. The worst-case where the seawall is destroyed was also considered in the simulations. The combined CSK/S-1 2018-2021 mean Up-Down ground deformation rates (see Figure 9) were also used to obtain rough estimates of the seawall height corrections due to subsidence for ten years from now (i.e., in 2030). However, the observed ground deformation (subsidence) rates along the Shanghai coastline from 2018 to 2021, which agreed with what was expected from previous work [105], were to the order of a few centimeters per year. Accordingly, the applied seawall height corrections did not significantly affect the estimates of the wave overtopping probabilities (see Equation (3)), and only marginally changed projections for the areas where the subsidence rate was higher ( Figure 9). Overall, the probability of wave overtopping remains very low in the northern sectors and more significant in the southern coastal sectors where the seawalls are lower. The simulation results that considered the failure of each segment are shown in Figures 11 and 12. As shown in Figure 10a, along the coast, high-coherent objects are significantly concentrated in the south-central part of Shanghai, reflected by high-density grids. In the southwest of Dishui Lake and the east-central part of Shanghai, the density of high-coherent objects is intermediate. In contrast, other regions have relatively low densities of coherent targets.

The Inundation Map for Different Segmentation Scenarios
Following Section 2, we split the coastline seawalls into forty non-overlapping segments and estimated the probability of wave overtopping and flood impacts in the case of an extreme weather event under the failure of each individual segment by considering the relevant seawall height and Equations (1) and (2). The impact is calculated by applying the LISFLOOD-FP hydrodynamic model [57,58], assuming that, for every selected sector, incoming waves overpass the barriers and the seawall sector has a severe failure. The worst-case where the seawall is destroyed was also considered in the simulations. The combined CSK/S-1 2018-2021 mean Up-Down ground deformation rates (see Figure 9) were also used to obtain rough estimates of the seawall height corrections due to subsidence for ten years from now (i.e., in 2030). However, the observed ground deformation (subsidence) rates along the Shanghai coastline from 2018 to 2021, which agreed with what was expected from previous work [105], were to the order of a few centimeters per year. Accordingly, the applied seawall height corrections did not significantly affect the estimates of the wave overtopping probabilities (see Equation (3)), and only marginally changed projections for the areas where the subsidence rate was higher ( Figure 9). Overall, the probability of wave overtopping remains very low in the northern sectors and more significant in the southern coastal sectors where the seawalls are lower. The simulation results that considered the failure of each segment are shown in Figures 11 and 12.
For segments 5 to 14 (Figures 9 and 10), the inundation areas were relatively low and distributed in a local region determined by the overtopping location. For segments 17 to 32, any failure would lead to extensive and continuous flooding along the coast from Pudong Airport to Dishui Lake. The deepest flooded regions (>3 m) were seen in a triangular-shaped area on the northeast side of Pudong Airport and around Dishui Lake. This is because these regions have a near-zero elevation, and thus will be quickly inundated. We can easily assess the damage extent from the inundation maps by calculating the flooded area. However, the more important factor is that the effects on inhabitants are still unknown due to the uneven distribution of the population and buildings. Overall, land in eastern Shanghai is more vulnerable to seawater intrusion. This region was reclaimed from the ocean after 2000, and its average elevation is lower than southern Shanghai. However, this region is mainly composed of farmland, wetland, and unused land. In contrast, more built-up areas are in the flooded area within southern Shanghai (Figure 10a). Furthermore, as mentioned in Section 2.3.5, we utilized the DP to measure the density of buildings in a flooded region.

Flood Risk Assessment of Wave Overtopping
We derived from the considered 40 scenarios described in Section 3.4.2 a risk map (Figure 13a) that relates the probability of wave overtopping (worst case) with the impact on the ground in terms of high-coherent objects located in the simulated inundated regions. The simulations considered the seawall heights and their potential displacements ten years from now (2030) using the obtained DInSAR-driven ground displacement time series. The probability of overtopping was relatively high (>7%) in segments 3 to 6. For segments 5 to 9, the density of high-coherent objects was more than twice that in other segments. The flooded domains were similar when overtopping occurred in segments 5 to 9, as shown in Figure 11. In this area, although the inundation areas were not large, the density of high-coherent objects was the highest when compared with other coastal areas. We focused on critical segments 5 and 6, which had high probability and density. There are many kinds of industrial plants found in the inundation areas of segments 5 and 6. Two ports (Luchaogang Port and Shanghai Southern Port) are also adjacent (see Figure  13b). The vertical deformations of the seawall were low for segments 5 and 6. However, because the seawalls were built before 2000 and are lower than the eastern seawalls, the projected seawall height in 2030 is the lowest, which leads to a high probability of overtopping. For segments 5 to 14 (Figures 9 and 10), the inundation areas were relatively low and distributed in a local region determined by the overtopping location. For segments 17 to 32, any failure would lead to extensive and continuous flooding along the coast from Pudong Airport to Dishui Lake. The deepest flooded regions (>3 m) were seen in a triangularshaped area on the northeast side of Pudong Airport and around Dishui Lake. This is because these regions have a near-zero elevation, and thus will be quickly inundated. We can easily assess the damage extent from the inundation maps by calculating the flooded area. However, the more important factor is that the effects on inhabitants are still unknown due to the uneven distribution of the population and buildings. Overall, land in eastern Shanghai is more vulnerable to seawater intrusion. This region was reclaimed from the ocean after 2000, and its average elevation is lower than southern Shanghai. However, this region is mainly composed of farmland, wetland, and unused land. In contrast, more built-up areas are in the flooded area within southern Shanghai (Figure 10a). Furthermore, as mentioned in Section 2.3.5, we utilized the D P to measure the density of buildings in a flooded region.

Flood Risk Assessment of Wave Overtopping
We derived from the considered 40 scenarios described in Section 3.4.2 a risk map (Figure 13a) that relates the probability of wave overtopping (worst case) with the impact on the ground in terms of high-coherent objects located in the simulated inundated regions. The simulations considered the seawall heights and their potential displacements ten years from now (2030) using the obtained DInSAR-driven ground displacement time series. The probability of overtopping was relatively high (>7%) in segments 3 to 6. For segments 5 to 9, the density of high-coherent objects was more than twice that in other segments. The flooded domains were similar when overtopping occurred in segments 5 to 9, as shown in Figure 11. In this area, although the inundation areas were not large, the density of high-coherent objects was the highest when compared with other coastal areas. We focused on critical segments 5 and 6, which had high probability and density. There are many kinds of industrial plants found in the inundation areas of segments 5 and 6. Two ports (Luchaogang Port and Shanghai Southern Port) are also adjacent (see Figure 13b). The vertical deformations of the seawall were low for segments 5 and 6. However, because the seawalls were built before 2000 and are lower than the eastern seawalls, the projected seawall height in 2030 is the lowest, which leads to a high probability of overtopping.

Discussion
The coastal reclamation areas of the twelve cities considered here were largest in Dalian (634.38 km 2 ), followed by Ningbo (451.13 km 2 ), Tangshan (409.65 km 2 ), Yancheng (407.93 km 2 ), Nantong (374.21 km 2 ), and Tianjin (246.25 km 2 ). Except for Ningbo, the other five cities were all located in the north of China. This result is consistent with Tian et al. [8] and Huang et al. [5], who found that from 1985 to 2010, more than 70% of coastal reclamation projects were implemented in northern China. In particular, two cities in Bohai Bay (Tangshan and Tianjin) have rapidly carried out reclamation projects in the past sev-

Discussion
The coastal reclamation areas of the twelve cities considered here were largest in Dalian (634.38 km 2 ), followed by Ningbo (451.13 km 2 ), Tangshan (409.65 km 2 ), Yancheng (407.93 km 2 ), Nantong (374.21 km 2 ), and Tianjin (246.25 km 2 ). Except for Ningbo, the other five cities were all located in the north of China. This result is consistent with Tian et al. [8] and Huang et al. [5], who found that from 1985 to 2010, more than 70% of coastal reclamation projects were implemented in northern China. In particular, two cities in Bohai Bay (Tangshan and Tianjin) have rapidly carried out reclamation projects in the past several decades up to 2015 [109]. We compared the delineated outlines of newly reclaimed land shown in Figure 3b,c with other studies for these two cites. The results over the same period, which were interpreted with Landsat images, showed highly similar patterns [6,12]. For the reclamation areas of the YRD, the high-intensity reclamation areas were mainly located on the coast of Yancheng and Nantong, the three islands in the Yangtze Estuary, southeastern Shanghai, and northern Ningbo [2,7,16,110], in line with our results.
We have demonstrated here that the simple cross-comparison of InSAR-driven DEMs produced at different epochs can give valuable (global) information on changes at the regional scale. We addressed the potential of InSAR technology to create worldwide topography maps with impacts on the study of ocean-reclaimed lands. DEM-assisted analyses were performed to quantify the extent of regions over China affected by floods. Flood risk increases in newly ocean-reclaimed lands were investigated by calculating the flood percentage of the reclaimed area and the ratio of reclaimed area flood zones to the total flood zones. On one hand, from this perspective, using newly-released DEMs (e.g., TDX90) is valuable to have global (at regional scale) estimates of the flood risk in coastal environments [107]. On the other hand, the vertical accuracy (to the order of a few meters) of the available DEMs is typically not accurate enough to implement detailed (localized) analyses, and the lack of data in areas close to water bodies both represent serious challenges that must be considered when studying the relations between extreme flood events and the state of consolidation of the ocean-reclaimed lands after the end of reclamation procedures. Seaward reclamation induces extreme events to occur earlier, and decreases the process of water-level and water-velocity reduction [30]. Before land reclamation, extensive continuous tidal wetlands (e.g., salt marshes and mangroves) can provide efficient wave damping [111,112]. In our work, we followed and improved upon previous approaches [29,94,113] in using global DEMs (i.e., TDX90) to map coastal inundation at middle-to-large spatial scales. The DEM-assisted differencing method can provide valuable information on coastal districts at medium-to-low resolutions. More accurate DEMs, derived from LiDAR (Light Detection and Ranging) missions can be determined to produce detailed spatially explicit inundation maps. However, LiDAR missions are costly. The generated DEMs have a limited spatial extent and result in time-consuming operations to run the hydrodynamic models and recover the inundated areas. The application of DInSAR technology allows us to easily estimate the extent of the areas to be potentially analyzed (e.g., EU Copernicus Sentinel-1 swaths are to the order of 100 × 200 km) and simultaneously perform detailed investigations over the critical areas previously identified by using simple DEM-assisted experiments. DInSAR-driven ground displacement time-series (with an accuracy of a few millimeters), in situ data (e.g., seawall heights), and other ancillary data (e.g., the seawater level) can jointly be used to evaluate the risk of flooding at the local scale.
Our work provides a general overview of the capabilities of InSAR technologies at both regional (using at most DEM-assisted methods) and local scales (using DInSAR-driven displacement maps) to obtain updated space-time information changes that have occurred over a set of coherent objects located in ocean-reclaimed lands. In relation to previous InSAR studies performed in the Shanghai area [45,55,[114][115][116][117], this work delineates a strategy to evaluate the risk maps of potential extreme flood events by considering the worst-case scenario that some of the barriers (seawall sectors) could fail due to wave action and be subsequently destroyed during an extreme event. To this aim, both the probabilities of wave overtopping and the impact of the potential inundations in the areas near the Shanghai coastline have been estimated.

Conclusions
This study analyzed flood risks in the coastal cities of China that have experienced significant land reclamation. First, we used a DEM-based method to analyze the coastal evolution and elevation changes in newly reclaimed regions on a regional scale. We provide statistics of the reclamation areas related to twelve high intensity reclaimed cities and the ratio of reclamation area to the total administrative area. Results show that, in total, 3167.6 km 2 of coastal land has been reclaimed in these cities. Among them, there are six cities with a reclamation area of more than 200 km 2 , mainly in northern China.
Moreover, we mapped the shape and hypsometry of newly reclaimed land and calculated the average elevation. Most of the ocean-reclaimed land belongs to the LECZ, which is expected to exacerbate flood risks. Quantitative analyses employ a simplified hydrodynamic model (LISFLOOD-FP) to map coastal inundation and calculate the flood zones in the whole city and reclaimed lands. A majority of the total flooded regions are in reclamation zones, and around 60% of the reclaimed land is likely to be flooded. These ratios for Yancheng, Nantong, Shanghai, and Wenzhou are more than 80%.
Subsequently, we focused on a local scale analysis relevant to the area of Shanghai. The multi-temporal SBAS approach and the minimum acceleration multi-track DInSAR-based combination method [36,51] were applied to update ground displacement maps of the coastal area of Shanghai from 2018 to 2021. The density of high-coherent objects within inundation areas was used to estimate the flood impacts of over 40 simulated flooding scenarios. These results demonstrated that two southern segments (segments 5 and 6) had the highest flood risk among the 40 segments. This result can be attributed to the fact that these sections were built before 2000 and at a lower height than those in northern Shanghai. Although the subsidence of seawalls in the eastern parts of Shanghai is much higher, the future seawall height on the east coast is still higher than that in the south. In addition, compared with the original seawall height, the observed ground subsidence and its actual rate (from 2018 to 2021) has minor impacts on the probability of wave overtopping over time. For segments 5 and 6, the inundation areas were not the largest among the 40 scenarios. However, within their inundation regions, there were numerous buildings, factories, and two artificial ports, leading to a high density of high-coherent object distribution and showing a high coastal flood impact. Shanghai is a unique situation compared to the other cities considered here, and is the only place we considered seawall subsidence and wave overtopping using the InSAR method.
This investigation shows the potential of using InSAR technology in evaluating flood risk in low-lying coastal regions and delineates a strategy to quantify it. However, it also reveals some significant challenges such as the need to have accurate ground-truth and ancillary data at the scale of the single structures affected by potential inundations to have more reliable quantitative results on the inundated areas in real scenarios, instead of considering worst-case scenarios by relaxing some of the hypotheses regarding the on-going ground displacements and observed phenomena. Future studies will focus on other coastal populated cities, especially those with continuous artificial seawalls.