Differential Synthetic Aperture Radar Interferometry for Land Movements Mapping Case of El Jebha (North of Morocco)
Total Page：16
File Type：pdf, Size：1020Kb
International Journal of Civil Engineering and Technology (IJCIET) Volume 9, Issue 5, May 2018, pp. 105–113, Article ID: IJCIET_09_05_013 Available online at http://iaeme.com/Home/issue/IJCIET?Volume=9&Issue=5 ISSN Print: 0976-6308 and ISSN Online: 0976-6316 © IAEME Publication Scopus Indexed DIFFERENTIAL SYNTHETIC APERTURE RADAR INTERFEROMETRY FOR LAND MOVEMENTS MAPPING CASE OF EL JEBHA (NORTH OF MOROCCO) Soukaina MEZIANE, Lahcen BAHI, Latifa OUADIF L3GIE Laboratory, Mohammadia Engineering School, Mohammed V University, Rabat, Morocco ABSTRACT An analysis of land movements is presented using remote sensing data through applying Differential Synthetic Aperture Radar Interferometry Processing (DInsAR).The method used is based on the DinSAR detection of ground subsidence. In this paper, the study focuses El Jebha area (North of Morocco) where occurs land instability because of the conjunction of several factors triggering landslides such as heavy rainfall, saturated soil, heterogeneous lithology with marl and clay soil especially in overlapping fractures and along feeder roads. The DInSAR method computes the displacement between two date’s single look complex IW acquisition of 05 November 2016 and 05 March 2017 acquired from Sentinel-1A satellite. The produced ground deformation maps allow detecting land subsidence and measure vertical displacements reaching up to 50 cm and uplift to -54cm in the observed time interval. The process applications are Sentinels Application Platform (SNAP) software and SNAPHU for unwrapping processing. DInSAR results with rain, geotechnics and geologic data explain the impact of natural and human activities on land movements. Key words: Mapping, Land Movements, Aperture Radar, El Jebha, Morocco Cite this Article: Soukaina MEZIANE, Lahcen BAHI, Latifa OUADIF, Differential Synthetic Aperture Radar Interferometry for Land Movements Mapping Case of El Jebha (North of Morocco), International Journal of Civil Engineering and Technology, 9(5), 2018, pp. 105–113. http://iaeme.com/Home/issue/IJCIET?Volume=9&Issue=5 1. INTRODUCTION El Jebha is located in northern Morocco on the Rif’s chain on its Mediterranean facade. It is a fragile area and vulnerable to land movements. Several studies show that this zone does not miss problems of instabilities and this because of several triggers [1-12]. There are a number of methods for mapping with GIS methods for vulnerable area to landslides and land movements in Mediterranean basins [13-15]. In parallel, several studies focused on DInSAR methods from the algorithm development to its direct application on determining http://iaeme.com/Home/journal/IJCIET 105 [email protected] Differential Synthetic Aperture Radar Interferometry for Land Movements Mapping Case of El Jebha (North of Morocco) displacements of the Earth ground. Among DInSAR algorithms figures Persistent Scatterers Interferometry (PSI) [16-19] and Small Baseline techniques [20]. This article discusses the ground deformation in El Jebha using Synthetic Aperture Radar (SAR). This Differential interferometric SAR is based on processing two single look complex interferometric wide images (SLC IW). Each pixel of SAR data contains a complex number (amplitude and phase) of the microwave field backscattered, which corresponds to resolution cell, projected on the ground [21]. In fact, through TOPS coregistration process and phase unwrapping, interferometric analysis allows obtaining land subsidence. This paper presents ground deformations mapping in El Jebha, in order to monitor these areas and to help strengthen them especially if it presents risks of human and material damage. 2. DESCRIPTION OF AREA STUDY El Jebha is located in the Northern Internal Rif. This zone is composed of three geological layers. Indeed, the layer of ghomarids rests in abnormal contact on the metamorphic substratum of the sheet of the Sebtides in the West. This substratum is composed mainly of peridotites surrounded by gneiss, mechachists and sericitoschists associated with sericitous quartzites and cipolins. These series end locally with metamorphic dolomites of the Triassic. While in the East, it is limited by the domain of flysches to the tablecloth Tisirene (According to the geological map of El Jebha). Figure 1 Extract of trust map of Rif chain El Jebha has a Mediterranean climate characterized by heavy rainfall. Maximum rainfall is in November and December, while July and August are dry months, with thunderstorms on the top. These rains can be continuous, with torrential rains [8]. They feed the highly fractured and weakened soil, causing infiltration of water to marl, marl-calcareous and clay formations. Due to the importance of precipitation [22], the imperviousness of the ground cover and the mountainous character makes that the runoff is relatively important in this zone and that the streams have irregular regimes characterized by a torrential flow in periods of floods. In addition, there is discordance in the stratigraphy marked by steep slopes and heterogeneous and contrasted lithology, especially in faults (figure1). There are N40-N50 oriented fractures with the most dominant to the crossover faults associated with the N60-N80 direction of the http://iaeme.com/Home/journal/IJCIET 106 [email protected] Soukaina MEZIANE, Lahcen BAHI, Latifa OUADIF Jebha-Chrafate fault [22]. Road improvements and quarrying on hillsides often lead to a break in the natural slope. They can also constitute corridors of circulation and infiltration of water. They modify the surface runoff and destabilize the slopes. 3. DATA AND METHODS 3.1. Satellite Data Sentinel-1 is a satellite for Earth observation developed by the European Space Agency [23] as part of the Copernicus program. The first satellite Sentinel-1A is into orbit in 03 April 2014 and the Sentinel-1B in 25 April 2016 in order to take over the ENVISAT and ERS missions. Each satellite carries a payload constituted by the synthetic aperture radar operating in C-band (wavelength 5.6 cm) which provides views with a resolution of up to 5 meters [23]. Compared to other SAR missions, Sentinel-1 provides short revisit time (Sentinel-1A with twelve days and six days with Sentinel-1B). The InSAR uses complex single look image, which is a digital SAR image, composed by columns and rows of small elements pixels forming a two-dimensional array [23]. Each pixel corresponds to cell resolution associated with a small area of the Earth’s surface. Each image has different rows and columns associated respectively with different azimuth locations and different slant range locations (also named Line of Sight direction LOS) [23]. The resolution cell coordinates in azimuth and slant-range depend only on the SAR system characteristics. In addition, the satellite supports four different operational imaging modes (Interferometric Wide swath (IW) mode, Extra Wide swath (EW) mode, Strip-Map (SM) mode, and Wave (WV) mode [23]. In this study, two single look complex Sentinel-1A TOPS mode Interferometric data in Single polarization IW mode for November 05th 2016 and March 5th 2017 are downloaded from Sentinels Scientific Data Hub (Table 1) and are processed with SNAP –S1TBX (Sentinel Application Platform –Sentinel 1 Toolbox) software [24] and with SNAPHU processing for unwrapping [25]. Figure 2 and 3 show the localization of El Jebha in footprint images and in SLC 20170503 Sentinel-1A image. Table 1 Data sets summary Acquisition Processing Orbit Sensor Stack Polarisation Date type Direction 05 nov2016 Sentinel1-A C Master VV Ascending 05march 2017 Sentinel1-A C Slave VV Ascending Figure 2 Footprints of 03March 2017 and 05November 2016 SLC IW Sentinel-1A images http://iaeme.com/Home/journal/IJCIET 107 [email protected] Differential Synthetic Aperture Radar Interferometry for Land Movements Mapping Case of El Jebha (North of Morocco) Figure 3 SLC Image of 2017-05-03 Sentinel-1A 3.2. Sentinel-1 TOPS Interferometric Processing DinSAR process includes several steps since the TOPS coregistration of images into a stack and the interferogram formation to phase filtering (earth flatting, topographic phase removal, Goldstein filter). Once unwrapping phase is computed, the phase is converted to displacement, which is geocoded and converted to (.kmz) format for visualizing map deformation on Google Earth. DInSAR method [26] proposes to observe deformation in El Jebha. Interferometric SAR derives the phase difference based on two images of Sentinel -1 observations taken between 05 November 2016 and 05 March 2017. According to the processing of interferometric SAR, those images are coregistered into a stack. Before coregistration, InSAR stack overview determines 05 November 2016 as master image and 05 march 2017 as slave image. Coregistration is essential for accurate determination of difference phases. Phase φ has the following linear dependence on the slant range coordinate r where λ is the SAR wavelength [26] (1) Considering that the phase of the transmitted signal is zero [26], the received signal that covers the distance 2R travelling from the satellite to the target and back is: radians (2) The SAR interferogram is complex multiplication of coregistrated SAR images. It is generated by cross-multiplying the master SAR image pixel with the complex conjugate pixel of the slave [26]. Thus, the interferogram amplitude is the amplitude of the first image multiplied by that of the second one, whereas its phase is the difference between the acquisitions. The interferometric phase variation ΔΦ is then proportional to ΔR