Jakarta Land Subsidence and Inundation Vulnerability Based on SAR Data

Transcription

1 Jakarta Land Subsidence and Inundation Vulnerability Based on SAR Data Agustan 1, Hartanto Sanjaya 1, Takeo Ito 2 1 Center of Technology for Natural Resources Inventory (PTISDA) Agency for the Assessment and Application of Technology (BPPT) BPPT 1 st Bulding, 20 th Floor, Jl. M.H. Thamrin No. 8 Jakarta, Indonesia agustan@bppt.go.id hartantosanjaya@gmail.com 2 Graduate School of Environmental Studies, Nagoya University Furo-cho, Chikusa-ku, Nagoya Japan takeo_ito@nagoya-u.jp Abstract: Heavy rains in mid-january 2013 caused terrible floods in Jakarta, the capital city of Indonesia. The floods inundated about 41 square kilometers, or eight percent of Jakarta. It is less than the 2007 event where floods inundated 231 square kilometer, however, the total material losses of the 2013 event was much higher, reaching 32 trillion rupiahs in 2013 compared to 4.3 trillion rupiahs in Apart from regular flooding, Jakarta is facing a condition of high rate land subsidence. Land subsidence is a hazard resulting in negative impacts and could lead to serious problems, such as increasing risk of flooding, cracking the buildings and infrastructures, destructing local groundwater systems, and generating tension cracks on land and reactivating faults. This research describes the relation of inundation vulnerability with land subsidence in Jakarta based on Synthetic Aperture Radar (SAR) data. Land subsidence is identified by implementing interferometry (InSAR) technique of ALOS- PALSAR data using GAMMA SAR software, whereas the inundated area of January event is revealed from RADARSAT data using NEST software. The principle behind using radar for flood mapping is that the strength of the radar backscatter depends on the roughness of the surface it interacts with. Water bodies generally create a very smooth and homogenous surface so the radar image over that area will be black or dark. To distinguish flooded areas from permanent water bodies, an image acquired during or after the flood is compared with another image acquired before the flood. Both images need to have been acquired from the same sensor geometry or orbit direction. By processing four years time span ALOS-PALSAR data, it is found that long-term observation gives clear information of land subsidence in Jakarta and tends to increase in term of area width and number of locations. By comparing RADARSAT data before and after the peak of heavy rain, it is found that areas affected by land subsidence were also inundated in longer time. Keywords: SAR, land subsidence, inundation, Jakarta 1. Introduction As a capital city with more then ten millions population, Jakarta faces a serious problem related to disaster i.e. flood during rainy season. Flood history in Jakarta is recorded on 1979, 1996, 1999, 2002, 2007, and During the flood event, inundated area usually occurred on the low elevated area that was estimated to cover 40% of Jakarta s total area or approximately 780 km2. In 2007, the flood cause 4.3 trillion rupiahs with 57 casualties (Sagala et al., 2010), whereas the 2013 event cause 41 casualties with 32 trillion rupiahs as summarized in Table 1.! "!

2 Table 1. Summary of Jakarta Flood Events (Modified from Sagala et al., 2013, BNPB and OCHA, 2013) Variables 1996 event 2002 event 2007 event 2013 event Causes The capacity of river is smaller than the incoming water s runoff. Low stage of river capacity and major canal are caused by the high conversion of area around these rivers and canals into settlement function, sedimentation, and illegal waste disposal Land use in urban areas which lots of buildings and settlements has led to the decreasing of land absorption ability as well as narrowing the river canal in downstream area. Beside of poor drainage system, flood was preceded by heavy rain from afternoon on Feb 1 to the next day on Feb 2. It was worsened by the high volume of water in 13 rivers in Jakarta which originated from Bogor Puncak Cianjur and the tide of Jakarta s sea water. Flood in Jan 2013 was less intense compared to 2007 but the rainfall was widely distributed upstream and downstream. 17th Jan was considered the worst but it is predicted to be continued and the city should be put on alert level 1 till 27 th Jan. Inundation Point (Google crisis map) Rainfall Intensity mm mm mm (geographically concentrated) mm (tentative) (geographically widely distributed) Evacuee (people) Dead (people) 30 thousand 380 thousand 398 thousand >100 thousands (tentative) Losses > IDR 1 Billion IDR 1.8 Trillion IDR 4.3 Trillion >IDR 32 Trillion The inundated and vulnerable area due to the increasing rainfall rate is identified in the same location from one event to another, mainly in the western region of Jakarta as illustrated in Figure 1. During the 2013 event, some locations are identified as suffering from longer inundated time. One factor that increases the inundated area is land subsidence factor (Murakami et al., 2012). Figure 1. Inundated Area due to the 2013 Jakarta flood event (Source: BNPB Map)! #!

3 Land subsidence is generally related to geological subsidence i.e. sediment consolidation due to its own weight and tectonic movements; or related to human activities such as withdrawal of ground water and geothermal fluid, oil and gas extraction from underground reservoirs, and collapse of underground mines (Raucoules et al., 2007; Yuill et al., 2009). The amount of subsidence or uplift can be estimated from the number of concentric fringes that appear in the interferogram as illustrated in Figure 2. Land subsidence in Jakarta had been observed by utilizing geodetic observation i.e. by Global Positioning System (GPS) observation technique (Abidin et al., 2008) and by interferometric synthethic aperture radar (InSAR) technique (Abidin et al., 2011; Ng et al., 2012; Chaussard et al., 2013). All of these previous studies concluded that Jakarta suffer a sinking phenomena due to its rapid subsidence rate, approximately 260 mm/year in northern part of Jakarta. This research describes the relation of inundation vulnerability with land subsidence in Jakarta based on Synthetic Aperture Radar (SAR) data. 2. Data and Methods To understand the expansion of subsided region in Jakarta, InSAR technique is applied to identify the location in regional scale. Eight (8) ALOS-PALSAR (Advanced Land Observing Satellite with Phased-Array Synthetic-Aperture Radar, L-band with wavelength ~ 23.6 cm) data in fine beam mode, ranging from 2007 to 2011 are processed with GAMMA SAR Software (Wegmuller and Werner, 1997). These data are obtained from Japan Aerospace Exploration Agency Earth Observation Research Center (JAXA-EORC). However, since the ALOS satellite had stopped the mission on 2011, the inundated area of the 2013 Jakarta flood event is assessed by processing RADARSAT-2 data (C-band with wavelength ~ 5 cm), obtained from Canadian Space Agency (CSA) with NEST SAR toolbox software (Engdahl et al., 2012). Figure 2. Phase difference of L-band satellite based InSAR technique due to land subsidence The ALOS-PALSAR data were observed on January 31 st 2007, August 3 rd 2007, September 18 th 2007, February 3 rd 2008, May 5 th 2008, June 20 th 2008, August 8 th 2009, August 11 th 2009, November 11 th 2010 and February 11 th The first epoch is set as the master, so the time difference for each interferograms are: 184, 230, 460, 506, 920, 1380 and 1472 days respectively. The RADARSAT data were observed on December 26 th 2012, February 2 nd 2013, and February 5 th The data is observed in fine beam mode and in SAR Georeferenced Fine (SGF) format. SGF format indicates that the data is converted to ground range, is multilooked processed, and is oriented in the direction of orbit path. As stated on Table 1 before, the Jakarta 2013 flood event is mainly caused by two days high rainfall on January 17 th to 18 th 2013.! $!

4 The InSAR technique produce an interferogram which is obtained by crossmultiplying, pixel by pixel, of two SAR images. The SAR image can be seen as a mosaic of small picture elements (pixels) that corresponds to a small area of the Earth s surface that can be defined as a resolution cell. Each pixel contains a complex number that carries amplitude and phase information about the microwave field backscattered by all objects in corresponding resolution cell projected on the ground. These kinds of information are stored in complex format by adapting IQ (In-phase and Quadrature) data format. Therefore, SAR image is also known as single look complex (SLC) that is composed of a regular grid with complex values or phasors (Hanssen, 2001) and can be decomposed into amplitude A or real R and phase " or imaginer I components. y = Ae i"...(1) where: y is the SLC data that represents the electric field of a plane electromagnetic wave, A is amplitude of the electromagnetic pulse, and " is phase angle. The amplitude represents the quantity of electromagnetic field scattered back grouped in each SAR image-sampling cell or pixel, whereas the phase represents an ambiguous measure of distance between sensor and each area on the ground corresponding to an image pixel (Raucoules et al., 2007). The first SAR image is called master and the second one slave. As a result, the interferogram amplitude is the amplitude of master image multiplied by that of slave image, whereas the interferometric phase is the phase difference between the images. Interferogram = y(master) " y(slave)...(2) Interferogram = A M A S e i(# M $# S )...(3) % I ( # = # M $ # S = tan $1 ' * = $ 4+(R $ R ) M S & R),...(4) # = # curv + # topo + # orb + # defo + # atm + # noise...(5) # = 4+, Bsin(- $.) + 4+HB /,R M sin- + 4+, B // + 4+, D + # atm + # noise...(6) where, " is the interferometric phase or phase difference derived from master (M) and slave (S) images in one point, R M and R S are the geometric distances of satellite to target both acquisition times, and " is the SAR system wavelength. The distance between the two satellite positions is known as the baseline B, and can be decomposed into parallel baseline B //, the component along the radar s line of sight; and perpendicular baseline B ", the component which is perpendicular to the line of sight. Look angle in one point " and " is angle of satellite baseline, whereas D is possible deformation. The ALOS PALSAR image has range pixel spacing of approximately 4.6 m, and 3.2 m for the azimuth pixel spacing. Further, co-registration and resampling of slave image to master image to conform the geometry were performed. Then phase co-registration by the range spectral shift and the azimuth common bandwidth filtering. It was found that the coregistration of the two images was 0.1. This indicated the geometry of both images was conformed in sub-pixel accuracy. This information ensured the high level of interferometric correlation. 3. Results and Discussion The interferogram was derived after the ALOS-PALSAR data was looked down into 4 and 6 in range and azimuth direction respectively. Hence, the interferogram pixel spacing was approximately 19 m in range and 19 m in azimuth. The Earth s curvature interferogram was removed based on the first term of Equation (6). Since the perpendicular baseline for all interferograms was small (<200 m), the flattened interferograms thus contained the deformation, atmospheric artifacts, and noise quantity in wrapped phase. The result can be seen in Figure 3.! %!

5 Figure 3. Series of the interferogram from ALOS-PALSAR based on January 31 st Interferogram (a), (b), (c), (d), (e), (f), and (g) has time difference of 184, 230, 460, 506, 920, 1380 and 1472 days respectively; whereas (h) illustrates the backscatter image of Jakarta. Figure 3 illustrates the interferogram of Jakarta and based on previous studies that stated the main subsided area is located on the northern and western part of Jakarta, so the region of interest region can be select to this area and illustrated in Figure 4.! &!

6 Figure 4. Interferogram of the northern and the western part of Jakarta; whereas (h) illustrates the optical image of Jakarta. One fringe in the interferogram means half of the wavelength and therefore it represents 11.8 cm. Figure 4 shows the expansion of the subsided region in the northern and the western part of Jakarta. In 4 years, at least 6 fringes are identified, and represent the target moving far from the satellite approximately 70.8 cm. This region then is validated with RADARSAT image and inundation map. Figure 5. RADARSAT (C-band) image of Jakarta, observed on December 26 th 2012, February 2 nd 2013, and February 5 th 2013.! '!

7 Figure 6. Inundation map of the 2013 Jakarta flood event issued by BNPB Figure 5 shows that the area indicated with red circle and green square is still inundated on February 2nd after the heavy rain on 17th and 18th January (after 14 days), After 17 days the area with red circle has dried up, in constrast to area in the green square, which is still inundated. Inundation map in Figure 6 also shows the inundated area on January 24th, 6 days after the peak of the rainfall, and it confirms that the subsided area still inundated. 4. Conclusion The intereferometric technique of satellite based SAR data to identify the deformed area in large scale is useful. Based on this technique, it is found that there are some locations in Jakarta that suffer subsidence phenomena. The subsided area in the northern and western part of Jakarta is inundated in longer time, more than 14 days after the peak of rainfall during the 2013 Jakarta flood event. Acknowledgement The authors thank the Canadian Space Agency for providing RADARSAT data; Japan Aerospace Exploration Agency-Earth Observation Research Center (JAXA-EORC) for providing ALOS-PALSAR data. References Abidin, H.Z., Andreas, H., Djaja, R., Darmawa, D., Gamal, M. (2008). Land subsidence characteristics of Jakarta between 1997 and 2005, as estimated using GPS surveys. GPS Solutions 12 (1), ! (!

8 Abidin, H., Andreas, H., Gumilar, I., Fukuda, Y., Pohan, Y. E., & Deguchi, T. (2011). Land subsidence of Jakarta (Indonesia) and its relation with urban development. Natural Hazards, 59(3), , Chaussard, E., Amelung, F., Abidin, H., & Hong, S.-H. (2013). Sinking cities in Indonesia: ALOS PALSAR detects rapid subsidence due to groundwater and gas extraction. Remote Sensing of Environment, 128(0), doi: Engdahl, M., Minchella, A., Marinkovic, P., Veci, L., and Lu, J. (2012). NEST: An esa open source Toolbox for scientific exploitation of SAR data, Geoscience and Remote Sensing Symposium (IGARSS), 2012 IEEE International, pp. 5322,5324, doi: /IGARSS Hanssen, R.F. (2001). Radar Interferometry: Data Interpretation and Error Analysis. Kluwer Academic Publishers, Dordrecht, 328 pp. Murakami, S., Kawase, M. and Komine, H. (2012). Land subsidence in Mekong delta by using InSAR and future development for vulnerability assessment in consideration of global climate change, International Workshop on Geo-engineering for responding to climate change and sustainable development of infrastructure (HUEGEO 2012), Dec Kusumaputra R, Rainayati M, Yendra M. (2010). Banjir Kanal Timur Karya Anak Bangsa. Grasindo, Jakarta (in Indonesian) Ng, A. H.-M., Ge, L., Li, X., Abidin, H. Z., Andreas, H., & Zhang, K. (2012). Mapping land subsidence in Jakarta, Indonesia using persistent scatterer interferometry (PSI) technique with ALOS PALSAR. International Journal of Applied Earth Observation and Geoinformation, 18(0), doi: Sagala, S., Lassa, J.A., Yasaditama, H. and Hudalah, D. (2013). The evolution of risk and vulnerability in Greater Jakarta: contesting government policy in dealing with a megacity s exposure to flooding, IRGSC Working Paper No OCHA. (2013). December January 2013 Monthly Humanitarian Bulletin, United Nations Office for the Coordination of Humanitarian Affairs. p_dec12-jan13_eng.pdf Raucoules, D., Colesanti, C., and Carnec, C. (2007). Use of SAR interferometry for detecting and assessing ground subsidence. Comptes Rendus Geosciences, 339(5), , Wegmuller, U. and Werner, C., GAMMA SAR Processor and Interferometry Software. Third ERS Symposium on Space at the service of our Environment, held in Florence, Italy, March, Compiled by T.-D. Guyenne and D. Danesy. European Space Agency, p Yuill, B., Lavoie, D., & Reed, D. (2009). Understanding subsidence processes in coastal Louisiana. Journal of Coastal Research, ! )!