ENVISAT Synthetic Aperture Radar Research Articles

Neural Network Pattern Recognition Experiments Toward a Fully Automatic Detection of Anomalies in InSAR Time Series of Surface Deformation

We present a neural network-based method to detect anomalies in time-dependent surface deformation fields given a set of geodetic images of displacements collected from multiple viewing geometries. The presented methodology is based on a supervised classification approach using combinations of line of sight multitemporal, multi-geometry interferometric synthetic aperture radar (InSAR) time series of displacements. We demonstrate this method with a set of 170 million time series of surface deformation generated for the entire Italian territory and derived from ERS, ENVISAT, and COSMO-SkyMed Synthetic Aperture Radar satellite constellations. We create a training dataset that has been compared with independently validated data and current state-of-the-art classification techniques. Compared to state-of-the-art algorithms, the presented framework provides increased detection accuracy, precision, recall, and reduced processing times for critical infrastructure and landslide monitoring. This study highlights how the proposed approach can accelerate the anomalous points identification step by up to 147 times compared to analytical and other artificial intelligence methods and can be theoretically extended to other geodetic measurements such as GPS, leveling data, or extensometers. Our results indicate that the proposed approach would make the anomaly identification post-processing times negligible when compared to the InSAR time-series processing.

Indian Ocean Crossing Swells: New Insights from “Fireworks” Perspective Using Envisat Advanced Synthetic Aperture Radar

Synthetic Aperture Radar (SAR) in wave mode is a powerful sensor for monitoring the swells propagating across ocean basins. Here, we investigate crossing swells in the Indian Ocean using 10-years Envisat SAR wave mode archive spanning from December 2003 to April 2012. Taking the benefit of the unique “fireworks” analysis on SAR observations, we reconstruct the origins and propagating routes that are associated with crossing swell pools in the Indian Ocean. Besides, three different crossing swell mechanisms are discriminated from space by the comparative analysis between results from “fireworks” and original SAR data: (1) in the mid-ocean basin of the Indian Ocean, two remote southern swells form the crossing swell; (2) wave-current interaction; and, (3) co-existence of remote Southern swell and shamal swell contribute to the crossing swells in the Agulhas Current region and the Arabian Sea.

The pattern of surface deformations in the Dąbrowa Basin (Poland) based on PSInSAR analysis

Slow, long-term ground deformations in the Dąbrowa Basin (southern Poland) were identified based on ERS-1 , ERS-2 and ENVISAT Synthetic Aperture Radar (SAR) images that were processed by means of Permanent Scatterer SAR Interferometry (PSInSAR). The Dąbrowa Basin is a region where two major factors can affect surface stability: intensive coal exploitation and neotectonic processes. In this study, in order to clarify the origin of surface deformations, the authors propose applying a newly developed algorithm of spatio-temporal PSInSAR data analysis. This analysis revealed that subsidence is a characteristic feature of the Dąbrowa Basin. A significant correlation exists between slow, long-term ground deformations and the location of the main tectonic structure of this region. The proposed spatiotemporal analysis of the PSInSAR data additionally showed some degree of correlation between mining activity and the studied deformations. This interconnection is a significant achievement of this study since the deformation values determined by means of PSInSAR were identified in earlier works solely on the basis of Dąbrowa Basin neotectonics. Normal 0 21 false false false PL X-NONE X-NONE /* Style Definitions */ table.MsoNormalTable {mso-style-name:Standardowy; mso-tstyle-rowband-size:0; mso-tstyle-colband-size:0; mso-style-noshow:yes; mso-style-priority:99; mso-style-parent:""; mso-padding-alt:0cm 5.4pt 0cm 5.4pt; mso-para-margin-top:0cm; mso-para-margin-right:0cm; mso-para-margin-bottom:10.0pt; mso-para-margin-left:0cm; line-height:115%; mso-pagination:widow-orphan; font-size:11.0pt; font-family:"Calibri","sans-serif"; mso-ascii-font-family:Calibri; mso-ascii-theme-font:minor-latin; mso-hansi-font-family:Calibri; mso-hansi-theme-font:minor-latin; mso-bidi-font-family:"Times New Roman"; mso-bidi-theme-font:minor-bidi; mso-fareast-language:EN-US;}

Sinkhole occurrence monitoring over shallow abandoned coal mines with satellite-based persistent scatterer interferometry

Unexpectedly occurring sinkholes caused by shallow submerged voids are one of the key problems of intensively built-up post-mining areas. Five percent of Poland’s territory is still under the influence of post-mining deformation. Impending sinkholes, which tend to develop randomly over very large areas, cannot be detected by using traditional geophysical methods. Moreover, geodetic measurement methods like Global Positioning System (GPS) analyses, tachymetry, laser scanning, or photogrammetry may be useful only for the registration of the dimensions and locations of sinkholes that have already occurred. Here, we investigate an area in Upper Silesia, Poland, where 345 sinkholes were recorded over a period of more than 20 years (1992–2013). Most of the events occurred in intensively built-up areas, and thus, these sinkholes posed a direct threat to the population. In the test area, 11 sinkholes were detected with a maximum depth of 16 m and maximum dimension of 25 m. The root cause of sinkhole formation was the collapse of post-mining shallow voids. We used satellite radar technology to detect and monitor ground movements potentially associated with the impending sinkholes. The findings showed that the application of Persistent Scatterer Interferometry can support the identification of zones where sinkholes will occur, given adequate spatio-temporal sampling. A stack of Envisat Synthetic Aperture Radar (SAR) images acquired between March 2003 and August 2010 was used in the analysis, and the results confirmed that precursory ground movements were detectable at an early stage of sinkhole development. The land subsidence rate observed prior to sinkhole collapse was not constant in time. Accelerated ground movements within 100 m of an observed sinkhole were detected. We conclude that satellite measurements may provide significant support in the early identification of areas prone to sinkhole occurrence.

Rapid detection to long ship wake in synthetic aperture radar satellite imagery

The maritime administrative department employs synthetic aperture radar (SAR) satellite remote sensing technology to obtain evidence of illegal discharge of ships. If the ship is discharged during navigation, it forms a long dark wake on the SAR image due to the suppression of the Bragg wave by the oil film. This study investigates key techniques for rapid detection of long ship wakes, thereby providing law enforcement agencies with candidate ships for possible discharge. This paper presents a rapid long ship wake detection method that uses satellite imaging parameters and the axial direction of the ship in images to determine the potential detection area of the wake. Then, the threshold of long ship wake detection is determined using statistical analysis, the area is binarized, and isolated points are removed using a morphological filter operator. The method was tested with ENVISAT Synthetic Aperture Radar and GF-3 SAR data, and results showed that the method was effective, and the overall accuracy of the decision reaches 71%. We present two innovations; one is a method that draws a Doppler shift curve, and uses the SAR imaging parameters to determine the detection area of the long wake to achieve rapid detection and reduce the image detection area. The other is where a classical linear fitting method is used to quickly and accurately determine whether the detected dark area is a long ship wake and realizes the twisted long ship wake detection caused by the sea surface flow field, which is otherwise difficult to detect by the traditional Radon and Hough transform methods. This method has good suppression performance for the dark spot false alarm formed by low speed wind region or upward flow. The method is developed for maritime ship monitoring system and will promote the operational application of maritime ship monitoring system.

Approaches for land cover monitoring over Europe based on backscatter and coherence properties of Envisat and ERS SAR data

The purpose of this study is the exploitation of Synthetic Aperture Radar (SAR) satellite data for land cover monitoring. The properties of the recorded radar data heavily depend on the type of land cover, season and weather conditions. Therefore, it is possible to utilize the variability of these factors, in order to develop various techniques and methodologies that can be used for classifying land surfaces. In this context, this paper proposes approaches, which include the application of mathematical expressions and application of thresholds on multi-temporal data, for recognizing and classifying various types of land cover, on the basis of ERS and Envisat C-band SAR backscatter and coherence properties. These can be useful for any kind of contemporary SAR data, such as those of the current two Sentinel-1 satellites. Although this study focuses on four main land cover types (urban, mountainous, agricultural-low vegetation and forested areas) over specific areas in Europe, the same principles can be extended worldwide, leading to useful insights for designing future SAR satellite missions or for establishing guidelines for in-depth studies of specific land cover types.

THE MONITORING AND SPATIAL-TEMPORAL EVOLUTION CHARACTERISTIC ANALYSIS FOR LAND SUBSIDENCE IN BEIJING

Abstract. At present the land subsidence has been the main geological disaster in the plain area of China, and became one of the most serious disaster that restrict the social and economic sustainable development, it also is an important content in the project of national geographic conditions monitoring. With the development of economy and society, Beijing as the capital of China has experienced significant population growth in the last few decades which led to over-exploitation of the ground water to meet the water demand of more than 20 million inhabitants, especially in the urban region with high population density. However, the rainfall and surface runoff can’t satisfy the need of aquifer recharge that product the land subsidence. As China’s political center and a metropolis, there are a lot of large constructions, underground traffic projects and complicated municipal pipeline network, and Beijing is also an important traffic hub for national railway and highway network, all of them would be threatened by the land subsidence disaster. In this article the author used twenty ENVISAT Synthetic Aperture Radar (SAR) images acquired in 2008 June–2010 August and ten TerraSAR images acquired in 2011 June–2012 September were processed with Small Baseline Subset SAR Interferometry (SBAS-InSAR) techniques, to investigate spatial and temporal patterns of land subsidence in the urban area of Beijing.

Strain Partitioning and Present‐Day Fault Kinematics in NW Tibet From Envisat SAR Interferometry

AbstractAn 8 year archive of Envisat synthetic aperture radar (SAR) data over a 300 × 500 km2 wide area in northwestern Tibet is analyzed to construct a line‐of‐sight map of the current surface velocity field. The resulting velocity map reveals (1) a velocity gradient across the Altyn Tagh fault, (2) a sharp velocity change along a structure following the base of the alluvial fans in southern Tarim, and (3) a broad velocity gradient, following the Jinsha suture. The interferometric synthetic aperture radar velocity field is combined with published GPS data to constrain the geometry and slip rates of a fault model consisting of a vertical fault plane under the Altyn Tagh fault and a shallow flat décollement ending in a steeper ramp on the Tarim side. The solutions converge toward 0.7 mm/yr of pure thrusting on the décollement‐ramp system and 10.5 mm/yr of left‐lateral strike‐slip movement on the Altyn Tagh fault, below a 17 km locking depth. A simple elastic dislocation model across the Jinsha suture shows that data are consistent with 4–8 mm/yr of left‐lateral shear across this structure. Interferometric synthetic aperture radar processing steps include implementing a stepwise unwrapping method starting with high‐quality interferograms to assist in unwrapping noisier interferograms, iteratively estimating long‐wavelength spatial ramps, and referencing all interferograms to bedrock pixels surrounding sedimentary basins. A specific focus on atmospheric delay estimation using the ERA‐Interim model decreases the uncertainty on the velocity across the Tibet border by a factor of 2.

Subsidence related to groundwater pumping for breweries in Merchtem area (Belgium), highlighted by Persistent Scaterrer Interferometry

ERS, ENVISAT and TerraSAR-X Synthetic Aperture Radar scenes covering the time span 1992–2014 were processed using a Persistent Scatterer technique to study the ground movements in Merchtem (25km NW of Brussels, Belgium). The processed datasets, covering three consecutive time intervals, reveal that the investigated area is affected by a global subsidence trend related to the extraction of groundwater in the deeper Cambro-Silurian aquifer. Through time the subsidence pattern is reduced and replaced by an uplift related to the rising water table attested by piezometers located in this aquifer. The subsidence is finally reduced to a zone where currently three breweries are very active and pump groundwater in the Ledo-Paniselian aquifer and in the Cambro-Silurian for process water for the production.

Large areas surface deformation analysis through a cloud computing P-SBAS approach for massive processing of DInSAR time series

We present in this work a methodology for computing surface deformation time series and mean velocity maps of large areas. Our approach relies on the availability of a multi-temporal set of synthetic aperture radar (SAR) data collected from ascending and descending orbits over an area of interest, and also permits us to estimate the vertical and horizontal (East-West) components of the Earth's surface deformation. The adopted methodology is based on an advanced cloud computing implementation of the differential SAR interferometry (DInSAR) Parallel Small Baseline Subset (P-SBAS) processing chain which allows the unsupervised processing of large SAR data volumes, from the raw data (level-0) imagery up to the generation of the corresponding DInSAR time series and maps. The solution presented, which is highly scalable, has been tested on ascending and descending ENVISAT SAR archives comprising approximately 400GB of data, which have been acquired over a large area of southern California (US) that extends over about 90,000km2. Such an input dataset has been processed in parallel by exploiting 280 computing nodes of the Amazon Web Services Cloud environment. The overall processing lasted about 8h and cost approximately $1900 USD. Moreover, to produce the final mean deformation velocity maps of the vertical and horizontal (East-West) displacement components of the whole investigated area, we also took advantage of the information available from external GPS measurements that permit us to account for possible regional trends not easily detectable by DInSAR and to refer the P-SBAS measurements to an external geodetic datum.The results presented clearly demonstrate the effectiveness of the proposed approach that paves the way to the extensive use of the available ERS-1/2 and ENVISAT SAR data archives. Furthermore, the proposed methodology can be particularly suitable to deal with the very large data flow provided by the Sentinel-1 constellation, thus permitting the extension of the DInSAR analyses at a nearly global scale.

Measurement of small co-seismic deformation field from multi-temporal SAR interferometry: application to the 19 September 2004 Huntoon Valley earthquake

ABSTRACTInterferometric synthetic aperture (InSAR) has been widely applied to natural disaster monitoring. However, it has limitations due to the influence of noise sources such as atmospheric and topographic artefacts, data processing errors, etc. In particular, atmospheric effect is one of the most prominent noise sources in InSAR for the monitoring of small magnitude deformations. In this paper, we proposed an efficient multi-temporal InSAR (MTInSAR) approach to measure small co-seismic deformations by minimizing atmospheric anomalies. This approach was applied to investigate the 18 September 2004 earthquake over Huntoon Valley, California, using 13 ascending and 22 descending ENVISAT synthetic aperture radar (SAR) images. The results showed that the co-seismic deformation was ±1.5 and ±1.0 cm in the horizontal and vertical directions, respectively. The earthquake source parameters were estimated using an elastic dislocation source from the ascending and descending acquisitions. The root mean square errors between the observed and modelled deformations were improved by the proposed MTInSAR approach to about 3.8 and 1.8 mm from about 4.0 and 5.2 mm in the ascending and descending orbits, respectively. It means that the MTInSAR approach presented herein remarkably improved the measurement performance of a small co-seismic deformation.

Multi-Stack Persistent Scatterer Interferometry Analysis in Wider Athens, Greece

The wider Athens metropolitan area serves as an interesting setting for conducting geodetic studies. On the one hand, it has a complex regional geotectonic characteristic with several active and blind faults, one of which gave the deadly M w 5.9 Athens earthquake on September 1999. On the other hand, the Greek capital is heavily urbanized, and construction activities have been taking place in the last few decades to address the city’s needs for advanced infrastructures. This work focuses on estimating ground velocities for the wider Athens area in a period spanning two decades, with an extended spatial coverage, increased spatial sampling of the measurements and at high precision. The aim is to deliver to the community a reference geodetic database containing consistent and robust velocity estimates to support further studies for modeling and multi-hazard assessment. The analysis employs advanced persistent scatterer interferometry methods, covering Athens with both ascending and descending ERS-1, ERS-2 and Envisat Synthetic Aperture Radar data, forming six independent interferometric stacks. A methodology is developed and applied to exploit track diversity for decomposing the actual surface velocity field to its vertical and horizontal components and coping with the post-processing of the multi-track big data. Results of the time series analysis reveal that a large area containing the Kifisia municipality experienced non-linear motion; while it had been subsiding in the period 1992–1995 (−12 mm/year), the same area has been uplifting since 2005 (+4 mm/year). This behavior is speculated to have its origin on the regional water extraction activities, which when halted, led to a physical restoration phase of the municipality. In addition, a zoom in the area inflicted by the 1999 earthquake shows that there were zones of counter-force horizontal movement prior to the event. Further analysis is suggested to investigate the source and tectonic implications of this observation.

PSInSAR Analysis in the Pisa Urban Area (Italy): A Case Study of Subsidence Related to Stratigraphical Factors and Urbanization

Permanent Scatterer Interferometry (PSI) has been used to detect and characterize the subsidence of the Pisa urban area, which extends for 33 km2 within the Arno coastal plain (Tuscany, Italy). Two SAR (Synthetic Aperture Radar) datasets, covering the time period from 1992 to 2010, were used to quantify the ground subsidence and its temporal evolution. A geotechnical borehole database was also used to make a correspondence with the detected displacements. Finally, the results of the SAR data analysis were contrasted with the urban development of the eastern part of the city in the time period from 1978 to 2013. ERS 1/2 (European Remote-Sensing Satellite) and Envisat SAR data, processed with the PSInSAR (Permanent Scatterer InSAR) algorithm, show that the investigated area is divided in two main sectors: the southwestern part, with null or very small subsidence rates (<2 mm/year), and the eastern portion which shows a general lowering with maximum deformation rates of 5 mm/year. This second area includes deformation rates higher than 15 mm/year, corresponding to small groups of buildings. The case studies in the eastern sector of the urban area have demonstrated the direct correlation between the age of construction of buildings and the registered subsidence rates, showing the importance of urbanization as an accelerating factor for the ground consolidation process.

An On-Demand Web Tool for the Unsupervised Retrieval of Earth’s Surface Deformation from SAR Data: The P-SBAS Service within the ESA G-POD Environment

This paper presents a web tool for the unsupervised retrieval of Earth’s surface deformation from Synthetic Aperture Radar (SAR) satellite data. The system is based on the implementation of the Differential SAR Interferometry (DInSAR) algorithm referred to as Parallel Small BAseline Subset (P-SBAS) approach, within the Grid Processing on Demand (G-POD) environment that is a part of the ESA’s Geohazards Exploitation Platform (GEP). The developed on-demand web tool, which is specifically addressed to scientists that are non-expert in DInSAR data processing, permits to set up an efficient on-line P-SBAS processing service to produce surface deformation mean velocity maps and time series in an unsupervised manner. Such results are obtained by exploiting the available huge ERS and ENVISAT SAR data archives; moreover, the implementation of the Sentinel-1 P-SBAS processing chain is in a rather advanced status and first results are already available. Thanks to the adopted strategy to co-locate both DInSAR algorithms and computational resources close to the SAR data archives, as well as the provided capability to easily generate the DInSAR results, the presented web tool may contribute to drastically expand the user community exploiting the DInSAR products and methodologies.

A First Assessment of the P-SBAS DInSAR Algorithm Performances Within a Cloud Computing Environment

We present in this work a first performance assessment of the Parallel Small BAseline Subset (P-SBAS) algorithm, for the generation of Differential Synthetic Aperture Radar (SAR) Interferometry (DInSAR) deformation maps and time series, which has been migrated to a Cloud Computing (CC) environment. In particular, we investigate the scalable performances of the P-SBAS algorithm by processing a selected ENVISAT ASAR image time series, which we use as a benchmark, and by exploiting the Amazon Web Services (AWS) CC platform. The presented analysis shows a very good match between the theoretical and experimental P-SBAS performances achieved within the CC environment. Moreover, the obtained results demonstrate that the implemented P-SBAS Cloud migration is able to process ENVISAT SAR image time series in short times (less than 7 h) and at low costs (about USD 200). The P-SBAS Cloud scalable performances are also compared to those achieved by exploiting an in-house High Performance Computing (HPC) cluster, showing that nearly no overhead is introduced by the presented Cloud solution. As a further outcome, the performed analysis allows us to identify the major bottlenecks that can hamper the P-SBAS performances within a CC environment, in the perspective of processing very huge SAR data flows such as those coming from the existing COSMO-SkyMed or the upcoming SENTINEL-1 constellation. This work represents a relevant step toward the challenging Earth Observation scenario focused on the joint exploitation of advanced DInSAR techniques and CC environments for the massive processing of Big SAR Data.

Application of a new algorithm using Doppler information to retrieve complex wind fields over the Black Sea from ENVISAT SAR images

Several algorithms have been proposed to retrieve near-surface wind fields from C-band synthetic aperture radar (SAR) images acquired over the ocean. They mainly differ in the way they retrieve the wind direction. Conventionally, the wind direction is taken from atmospheric models or is extracted from the linear features sometimes visible in SAR images. Recently, a new wind retrieval algorithm has been proposed, which also includes the Doppler shift induced by motions of the sea surface. In this article, we apply three wind retrieval algorithms, including the one using Doppler information, to three complex wind events encountered over the Black Sea and compare the SAR-derived wind fields with model wind fields calculated using the high-resolution weather research and forecasting (WRF) model. It is shown that the new algorithm is very efficient in resolving the 180° ambiguity in the wind direction, which is often a problem in the streak-based wind retrieval algorithms. However, the Doppler-based algorithm only yields good results for wind directions that have a significant component in the look direction of the SAR antenna. Furthermore, it is dependent on good separation of the contributions to the Doppler shift induced by surface currents and wind-related effects (wind drift and wind-sea components of the ocean wave spectrum). We conclude that an optimum wind retrieval algorithm should consist of a combination of the algorithms based on linear features and Doppler information.

A Novel Variable Index and Excision CFAR Based Ship Detection Method on SAR Imagery

When applying the constant false alarm rate (CFAR) detector to ship detection on synthetic aperture radar (SAR) imagery, multiple interferers such as upwelling, breaking waves, ambiguities, and neighboring ships in a dense traffic area will degrade the probability of detection. In this paper, we propose a novel variable index and excision CFAR (VIE-CFAR) based ship detection method to alleviate the masking effect of multiple interferers. Firstly, we improve the variable index (VI) CFAR with an excision procedure, which censors the multiple interferers from the reference cells. And then, the paper integrates the novel CFAR concept into a ship detection scheme on SAR imagery, which adopts the VIE-CFAR to screen reference cells and the distribution to derive detection threshold. Finally, we analyze the performances of the VIE-CFAR under different environments and validate the proposed method on both ENVISAT and TerraSAR-X SAR data. The results demonstrate that the proposed method outperforms other existing detectors, especially in the presence of multiple interferers.

Offshore wind climatology based on synergetic use of Envisat ASAR, ASCAT and QuikSCAT

The offshore wind climatology in the Northern European seas is analysed from ten years of Envisat synthetic aperture radar (SAR) images using a total of 9256 scenes, ten years of QuikSCAT and two years of ASCAT gridded ocean surface vector wind products and high-quality wind observations from four meteorological masts in the North Sea. The traditional method for assessment of the wind resource for wind energy application is through analysis of wind speed and wind direction observed during one or more years at a meteorological mast equipped with well-calibrated anemometers at several levels. The cost of such measurements is very high and therefore they are only sparsely available. An alternative method is the application of satellite remote sensing. Comparison of wind resource statistics from satellite products is presented and discussed including the uncertainty on the wind resource. The diurnal wind variability is found to be negligible at some location but up to 0.5ms−1 at two sites. Synergetic use of observations from multiple satellites in different orbits provides wind observations at six times in the diurnal cycle and increases the number of observations. At Horns Rev M2, FINO1 and Greater Gabbard satellite and in situ collocated samples show differences in mean wind speed of −2%, −1% and 3%, respectively. At Egmond aan Zee the difference is 10%. It is most likely due to scatterometer data sampled further offshore than at the meteorological mast. Comparing energy density with all samples at Horns Rev M2 shows overestimation 7–19% and at FINO1 underestimation 2–5% but no clear conclusion can be drawn as the comparison data are not collocated. At eight new offshore wind farm areas in Denmark, the variability in mean energy density observed by SAR ranges from 347Wm−2 in Sejerøbugten to 514Wm−2 at Horns Rev 3. The spatial variability in the near-shore areas is much higher than at areas located further offshore.

Large‐scale inflation of Tungurahua volcano (Ecuador) revealed by Persistent Scatterers SAR interferometry

AbstractThe Tungurahua volcano, in Ecuador, has been experiencing a substantial activity period since 1999, with several eruptions, including those of 2006 and 2008. We use a persistent scatterers approach to analyze a time series of Envisat synthetic aperture radar (SAR) data over the period 2003–2009, to investigate surface deformation in the region of the volcano. We measure a continuous large‐scale uplift with a maximum line of sight displacement rate of about 8 mm/yr, which is the first evidence of a sustained inflation in the Andes for an active volcano encompassing several eruptions. We model this signal as magma emplacement in a permanent storage zone at 11.5 km below sea level, with a net inflow rate of 7 million m3/yr. The paroxysmal eruptions in 2006 and 2008 did not seem to disrupt this long‐term signal. However, we observe significant deformation during the 2006 eruption consistent with an additional intrusion of 4.5 million m3 of magma.

Shoreline movement monitoring based on SAR images in Shanghai, China

In this article, a divide-and-combine approach and multiscale normalized cut techniques are combined to extract water–land boundaries from synthetic aperture radar (SAR) images. A validation study shows that waterline extraction from SAR using the proposed algorithm is consistent with that derived from Landsat optical remote-sensing imagery. Analysis of six European Remote Sensing-1/2 (ERS-1/2) satellite and Envisat SAR images acquired at spring tides in Lingang New City, a new land-reclamation city in the suburb of Shanghai, China, between 1993 and 2005 shows the shoreline movement of the study region. We find that the shoreline had minimum seaward movement in the region due to natural forces between 1993 and 1996. Between 1996 and 1999, the shoreline moved more in the eastern and slower in the western part due to artificial impelling siltation and deposition of silt in the Yangtze River. Later, between 1999 and 2005, the shoreline significantly changed following commencement of a tidal flat reclamation project.