Multitemporal Sentinel-1SAR Interferometry for Surface Deformation Monitoring near High Dam in Aswan, Egypt

The area in the downstream of Aswan High Dam is subjected to many factors that may cause a deformation on the earth surface. The existence of Lake Naser, the hesitated water table and its bad effect and the tectonic activities increase the possibility of surface elevation changes. In this study, Multitemporal Synthetic Aperture Radar (SAR) images and Differential SAR Interferometry (DinSAR) technique were used to derive terrain deformation measurements over the area near Aswan high dam. The data were processed using SNAP, Sentinel Application Platform and a two-pass interferometry technique with the help of the Shuttle Radar Topography Mission (SRTM) is performed. Based on the results, all interferograms performed have high coherence while the pixels represent water surface had a low coherence value. The surface elevation changes range between -6.5 cm and +2.4 cm in 2017, while in 2018 it ranges between -5.2 cm and +6.6 cm. it is recommended to use a longer observation time period with three-pass interferogram technique to give the independence of the external low-resolution DEM and map the annual rates of movement.


Introduction
Radar interferometry is an effective method for obtaining Digital Elevation Models (DEM) from SAR images by exploiting the phase information of the signal backscattered by the terrain. When using differential mode, it can also provide information on changes in terrain caused by various phenomena such as earthquakes, landslides, subsidence, variation in the ground water table, etc [1][2][3][4]. Differential Synthetic Aperture Radar Interferometry (DInSAR) can be used to measure ground deformation by relying on the phase difference between multi temporal SAR acquisitions, with accuracies in the range of mm to cm [5,6]. Although the principles have been stated long time ago, radar interferometry did not really used until the 1990s, with the launch of the ERS-1 satellite followed by the ERS-2 satellite (tandem mission). It must be noted, however, that this technology is not yet widely used in the field of remote sensing in general and radar in particular, due to the difficulty of obtaining data, which is not free of charge, or due to the sophisticated software and methodologies. In April 2014 the European Union's Copernicus programme lunched the first sentinel-1A followed by sentinel-1B, each carries a C-band Synthetic Aperture Radar. The main advancement of Sentinel-1 is the new imaging technique, i.e., the Terrain Observation by Progressive Scans (TOPS) [7]. This new scan technique allows to cover large footprints (about 250 km across the orbit track) by three overlapping sub-swaths. The data acquired by this satellite can be easily downloaded for free from The Copernicus Open Access Hub.
The Aswan High Dam is considered a unique structure among all the large irrigation and electric power projects in the world, see Figure 1. It was chosen to be a rock fill dam provided with an impervious core in its main part. The construction of the dam was completed in 1971 with the financial and technological assistance of the Soviet Union. The lake behind the dam, Lake Naser, is the major freshwater body supplying Egypt with water for different purposes. The lake now provides an important source of fish for Aswan and the rest of Egypt [8]. Lake Naser has a length of about 500 Km, 325 Km in Egypt and 175 Km in Sudan. It can store up to 130 billion m3 of water, water level reach 182-m above sea level in case of full storage scenario [9].
At the beginning of the water year (1 August), the water level is kept at 175 m above mean sea level to fulfil high and low flow requirement. When the water level upstream reaches an elevation between 178 m and 183 m, excess water is directed to the Toshka Depression [10]. Furthermore, the lake is located in a very hot, dry climate with an annual evaporation ranging from 2.1 to 2.6 m/y [11]. The annual rainfall over the lake is negligible [12]. According to El Shirbeny at el, 2018. The lake surface area has shrunk from 5637.53 km 2 in january 2015 to 4971.12 km 2 in January 2016 [13]. The operation rules of Nasser reservoir aim at ensuring adequate water supply and the safety of the Aswan High Dam. At the beginning of the water year (1 August), the water level is kept at 175 m AMSL to fulfil high and low flow requirement. When the water level upstream reaches an elevation between 178 m and 183 m, excess water is directed to the Toshka Depression and, if necessary, by means of the emergency spillways on the western bank of the Nile [10]. The maximum retention at 180 m AMSL is obtained in November and, subsequently, the reservoir levels decrease from January to July as water is released. In fact, as the water level increase in the summer, wet land area appeared everywhere at the areas have low terrain around the downstream of high dam.  Figure 2(a) in El Shalaal region near Kima factory, appeared in satellite images as dry soil, it was a mine for sand then the water table in this area raised up. The groundwater rising phenomena that happened can be directly refereed to the cessation of groundwater pumping from El-Shallal wells, and to the reduction of pumping from Kima factory wells [14]. Figure  2 (b) shows the pond accumulated there. Figure 3(a,b) shows the rising groundwater level that affecting many areas.
On November 14, 1981, an earthquake with magnitude 5.2 occurred in the Kalabsha area, which lies in the western side of Lake Nasser at a distance of about 60 km south of Aswan High Dam [15]. This earthquake occurred in depth of 25 km on the east-west Kalabsha fault beneath Gebel Marawa. From that date this area was considered to be seismically active. The seismic activity in this area might be related to both tectonic activities along these active faults and/or reservoir induced seismicity due to the Nasser Lake [16,17].
Fluctuations of the water table through the year, the huge amount of water body stored behind the dam, the water body accumulated in El Shalaal region and the seismicity of the site make studying terrain deformation and monitoring the surface an essential and crucial issue. The objective of this study is to use Multitemporal Synthetic Aperture Radar (SAR) images and Differential SAR Interferometry (DinSAR) technique to derive terrain deformation measurements over the area near Aswan high dam.

Interferometry and Differentia Interferometry
The radar image contains data of a complex nature, namely, the amplitude of the signal and also its phase. The phase is considered a valuable source of information on the change of relief in the image area. The measurement of the phase is at the base of the so-called interferometry Radar technique. This technique has been proposed and detailed for topographic applications by Zebker and Goldstein [18]. The principle of interferometry is the measurement of the phase difference between signals coming from the same object and picked up in two different positions by the radar. This phase difference is directly related to the altitude of the illuminated object.
The technique of differential interferometry was introduced by Gabriel et al, 1989 [19]. It is used to measure sub centimeter variations at altitude in the relief of a nominal terrain. The basic concept is to use radar images of different periods to detect any variation in altitude in the terrain. Technically, via the subtraction of two interferograms obtained before and after the manifestation of the phenomenon that caused the deformation. After the subtraction, the remaining represent the magnitude of the altimetric change that the terrain has subjected. According to Massonnet et al, 1993, in order to take advantage of the technique of interferometry, two conditions are decisive they are: the two images must have a high correlation coefficient between them and the target must have the same physical characteristics during acquisition times [20].

Interferogram Formation
A satellite SAR, when repeating orbits, can observe the same area at different times from slightly different look angles. The nominal repeat interval of the Sentinel1 acquisitions was 12 days. The distance between the two orbits in the plane perpendicular to the orbit is called the interferometer baseline (see Figure 3) and its projection perpendicular to the slant range is the perpendicular baseline (Bn). When some of the point on the ground slightly change their relative position in the time interval between two SAR observations (orbit 1 at time T1 and orbit 2 at time T2), the interferometric phase (∆ ), with respect to notations in figure 4, can be calculated as equation (1).
Where: ∆ is the relative scattered displacement projected on the slant range direction.
The interferogram is defined as the product of the complex values of a master image and complex conjugate of a slave image, i.e. the amplitudes of the corresponding pixels are averaged and the difference of the phase values for each point in the image is calculated. In interferometry, the perpendicular baseline plays a very important role in the application point of view. Particularly, when the application is related to the height or altitude of any target, the altitude of ambiguity (ha) is defined as the altitude difference that generates an interferometric phase changing of 2π after interferogram flattening. The altitude of ambiguity is inversely proportional to the perpendicular baseline, see equation 2. The larger the perpendicular baseline, the more is the sensitivity of the interferogram phase to topography. However, the smaller the perpendicular baseline, the more is sensitivity of the interferogram phase to the deformation. Where: is the perpendicular baseline, and is the incident angle.

Differential Interferometry Methods
There are four methods for measuring terrain elevation changes using differential interferometry. The first method used two images, by which one image represents the terrain before change, and the second is acquired after this change. However, both acquisitions must have exactly the same geometry, otherwise this technique will not produce any results. Areas that have no changes will have the same phase and a zero-phase difference. On the other hand, those with an altitude changes will backscatter waves with different phases. The rate of change can then be deduced geometrically. This technique remains difficult to achieve because even a small shift in the trajectory will certainly induce an error. The second method, named two-pass interferometry, can be implemented if we have two images beside a DEM for the studied area. The DEM represents the terrain before the change, one or both images are acquired after this change. An interferogram is simulated from the DEM and subtracted from the one generated based on the two images. The third method is the three-pass interferometry where tree images are used. In this method, the first two images of the scene must be taken before the altimetric variation of the relief, while third image is acquired after this change. Two interferograms are generated: one from the first two images, and the other from the second and third images, and then these two interferograms are subtracted. This technique results in a double difference interferogram. The fourth method is the fourth-pass interferometry that can be achieved with two separate couples of images. The first couple must be acquired before the change, and the second is taken after. By subtraction of the two interferograms obtained from the two pairs, the terrain elevation changes are then carried out.

Study Area
The study area is located in a downstream of Aswan High dam in the Eastern side of the Nile river, in Aswan, Egypt, bounded by latitudes 24°01'30" and 24°06'30"N and longitudes 32°52' and 32˚56'E. (Figure 2). This area can be classified as highlands where the elevations vary from 164 to 188 m above the main sea level. Few wadis incise this area and directed from east to west such as Wady El-Keimab and Wady Al-Haytah. The high land areas are underlain by basement rocks and or by Nubian sandstones. This region is subjected to an abrupt increase in the ground water levels and continuing up till now. Besides, it is considered to be seismically active place, based on the Egyptian National Seismic Network (ENSN) observations.

Data
This study was carried out by using the C-band SAR, wavelength λ = 5.6 cm. Data acquired in ascending orbit by Sentinel-1A (S1A) which was launched on 3 April 2014 and Sentinel-1B (S1B) on 25 April 2016. Their data are free downloadable from the European Space Agency (ESA) [21]. Dataset consists of seven SAR images acquired from October 2016 until January 2019 as shown in table 1. They are characterized by small spatial and temporal baselines. Repeating cycle between observations varies between 6 and 12 days. Level1 interferometric Wide (IW) Single Look Complex (SLC) data were used. Image pairs for interferogram conduction were selected thanks for the Interferometry Baseline Tool Published by Alaska Satellite Facility (ASF), University of Alaska Fairbanks (UAF) [22].

Methods
In this study we will use the two-pass interferometry with the help of the Shuttle Radar Topography Mission (SRTM) 3 second DEM for the study area, see figure 5. The datasets are required to generate a real interferogram in each time interval. The differential interferometric processing aims at separating the topographic and displacement term in an interferogram. The topographic phase has to be removed to identify the displacement component. Hence, two interferograms are formed. The master image was acquired on October 10, 2016, and the slave image depending on the time interval for formation the DEFO interferogram, where the TOPO interferogram is generated from the DEM. So, we have five deferential interferograms to estimate the displacement during the time intervals between the dates of these images. The data were processed using SNAP, Sentinel Application Platform. SNAP is an Eesa opensource software for scientific exploitation of earth observations. It is available to download from step web site [step]. The use of SNAP in tis study was due to the following technological innovations: Extensibility, Portability, a Graph Processing Framework, and Tiled Memory Management. First of all, images used in this study were resized to fit the study area using TOPS split tool (sub-swath IW2 and paths from 3-5 for (S1A) and from 4 -6 for (S1B) were chosen). The interferometric chain involves number of steps. The first step in the processing is applying orbit file to update the image with the precise orbit position and velocity information which will be available some days after acquisition date. The second step is co-registration of the tow sentinel-1 images.
Co-registration insures that each ground target contributes to the same pixels in range and azimuth in both master and slave image. This step includes both Back Geocoding and Enhanced Spatial Diversity. These steps insure that the images are precisely co-registered. The next step is to compute interferograms. Then Deburst is to be made to remove gaps between paths. Multilook is used to reduce speckles and to increase the image interpretability. In order to reduce the phase noise golden phase filter is performed. To apply the phase unwrapping the output in this stage must be exported to the Snaphu tool and the product will be imported to convert the phase into displacement. The final step is to make terrain correction for the displacement map using Range Dopplar Terrain Corection tool. A typical workflow of the proposed methodology is shown in figure 6.

Results and Discussion
The debursted interferograms contain interferometric phases band and coherence bands. In the phase bands represents the phase difference between the two images, after debursting, it is appeared as a spatially continuous image. The coherence map is an indirect measure of the quality of the interferogram. The coherence of all interferogram was very high on a scale from zero to one except the pixels represent water in the ground which appear to have low coherence. The next step is differential interferogram formation where the Topographic phase must be removed to obtain the deformation. Figure 7 shows the Topographic phase VV polarization removed from the differential interferogram created from orbits 13430 and 2534, i.e. between 10 and 16 October 2016. The maximum height removed was 28.755 m.
Snaphu was used to unwrap the phase, and the produced image represents the displacement in meter but in the lineofsight direction. The final step is to project the displacement using terrain correction tool. Table 2 presents the statistics of the surface elevation changes for all interferograms, while Figures 8, and 9 show a selection (one per year) of the accumulated displacement extracted from the differential interferograms during 2017 and 2018 respectively.

Conclusions
In this study the interferometric processing of the sen-tinel-1 TOPS data was presented, a straight forward ap-proach was followed to reach the goal of determining the surface elevation changes in the area near Aswan High Dam, Egypt. All interferograms performed have high coherence values, while the pixels represent water sur-face had low coherence values. High coherence permit-ting the retrieval of a time-series of displacement in a good manner. Based on the results, the surface elevation changes range between -6.5 cm and +2.4 cm in 2017, while in 2018 it ranges between -5.2 cm and +6.6 cm. The maximum height removed in differential interferogram due to topography was 28.755 m. However, a long-er observation time period is necessary because these data do not permit so far to estimate the annual rates of movement. The low spatial resolution of the DEM used prevent of reveal many spatial details and possibility to identify non-uniform motion became difficult. It is recommended that, a specific in-situ experiments are need-ed, and it should be well planned to validate the InSAR results. For further work, it will be better to use three-pass interferometry to give the independence of the external low-resolution DEM.