Next Article in Journal
SAR Target Classification Based on Sample Spectral Regularization
Next Article in Special Issue
Improvement of Earthquake Risk Awareness and Seismic Literacy of Korean Citizens through Earthquake Vulnerability Map from the 2017 Pohang Earthquake, South Korea
Previous Article in Journal
High-Quality UAV-Based Orthophotos for Cadastral Mapping: Guidance for Optimal Flight Configurations
Previous Article in Special Issue
Integration of InSAR Time-Series Data and GIS to Assess Land Subsidence along Subway Lines in the Seoul Metropolitan Area, South Korea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Land Subsidence Susceptibility Mapping in Jakarta Using Functional and Meta-Ensemble Machine Learning Algorithm Based on Time-Series InSAR Data

by
Wahyu Luqmanul Hakim
,
Arief Rizqiyanto Achmad
and
Chang-Wook Lee
*
Division of Science Education, Kangwon National University, Gangwon-do, Chuncheon-si 24341, Korea
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(21), 3627; https://0-doi-org.brum.beds.ac.uk/10.3390/rs12213627
Submission received: 10 September 2020 / Revised: 23 October 2020 / Accepted: 3 November 2020 / Published: 4 November 2020

Abstract

:
Areas at risk of land subsidence in Jakarta can be identified using a land subsidence susceptibility map. This study evaluates the quality of a susceptibility map made using functional (logistic regression and multilayer perceptron) and meta-ensemble (AdaBoost and LogitBoost) machine learning algorithms based on a land subsidence inventory map generated using the Sentinel-1 synthetic aperture radar (SAR) dataset from 2017 to 2020. The land subsidence locations were assessed using the time-series interferometry synthetic aperture radar (InSAR) method based on the Stanford Method for Persistent Scatterers (StaMPS) algorithm. The mean vertical deformation maps from ascending and descending tracks were compared and showed a good correlation between displacement patterns. Persistent scatterer points with mean vertical deformation value were randomly divided into two datasets: 50% for training the susceptibility model and 50% for validating the model in terms of accuracy and reliability. Additionally, 14 land subsidence conditioning factors correlated with subsidence occurrence were used to generate land subsidence susceptibility maps from the four algorithms. The receiver operating characteristic (ROC) curve analysis showed that the AdaBoost algorithm has higher subsidence susceptibility prediction accuracy (81.1%) than the multilayer perceptron (80%), logistic regression (79.4%), and LogitBoost (79.1%) algorithms. The land subsidence susceptibility map can be used to mitigate disasters caused by land subsidence in Jakarta, and our method can be applied to other study areas.

Graphical Abstract

1. Introduction

Several cities in Indonesia suffer from degradation at the ground level of buildings, known as land subsidence [1,2]. In Jakarta, this process has had a severe impact on urban infrastructure, leading to cracks in buildings, roads and damage to drainage systems [3]. These conditions are problematic because land subsidence may expand coastal flood areas due to sea-level rise [4]. Heavy monsoon rainfall [5] has caused frequent river flooding; if this occurs again, Jakarta could be submerged entirely underwater [1,4].
Recent studies of land subsidence in Jakarta have used various geodetic measurement methods, such as leveling surveys [3] and a global positioning system (GPS) surveys [6,7]. These studies indicated that excessive groundwater extraction is the leading cause of land subsidence and compaction to the vulnerable aquifer system [8]. This compaction may be exacerbated by natural consolidation since Jakarta’s landform mostly comprises young alluvial soils that cannot support the weight of human-made structures [9,10]. Therefore, it is essential to monitor land subsidence in Jakarta to predict further possible occurrences and mitigate damage [11,12].
Over the last decade, land subsidence susceptibility maps have been generated using geological, geomorphological, topographical, and hydrological data; these are considered the main factors influencing land subsidence [11,13,14]. Various methods are used to generate land subsidence susceptibility maps, including frequency ratios (FR) [12,15], weight of evidence (WOE) [16,17], logistic regression (LR) [18], evidential belief functions [11], analytical hierarchy processes (AHP) [19], support vector machines (SVM) [14,20], decision trees [21,22], fuzzy logic [12,23], adaptive neuro-fuzzy inference systems (ANFIS) [24,25], and artificial neural networks (ANN) [26,27]. In general, using a single modeling method leads to lower predictive accuracy than an ensemble method that uses a combination of models and machine learning algorithms [13,28]. Machine learning algorithms have the advantage of finding unpredictable relationships in datasets at multiple scales and have, thus, been recommended to obtain accurate land subsidence susceptibility maps [22].
Another challenge in generating land subsidence susceptibility maps is the low availability of land subsidence inventory maps. In this study, a land subsidence inventory map was generated via time-series interferometry synthetic aperture radar (InSAR) analysis. This technique can be applied to measure the displacement of the earth’s surface with an accuracy of up to millimeters per year by improving the selection of coherent pixels and reducing atmospheric noise [29,30,31,32,33]. It has been widely used to monitor land subsidence and generate land subsidence maps, for example, in large cities in Mexico [33,34], Kurdistan, and Iran [35], in Yuncheng Basin, China [36], and in coastal cities and areas such as Venice, Italy [37], New Orleans, United States [38], and Shanghai, China [39]. The InSAR algorithm has been successfully applied to earthquakes [40], volcanic activities [41,42], crustal deformation [31], landslides [43], manmade deformations [44], excessive groundwater extraction [3], mining activities [45], and natural consolidation of young alluvial soil [9].
The recent studies of monitoring land subsidence in Jakarta using interferometry synthetic aperture radar (InSAR) techniques was conducted using the Small Baseline Subset (SBAS) algorithm from 2007 to 2009 [1] and using Geodesy and Earth Observing Systems-Persistent Scatterer Interferometry (GEOS-PSI) algorithm from 2007 to 2010 [46]. Both studies utilized Advanced Land Observing Satellite (ALOS) phased array type-L synthetic aperture radar (PALSAR) data to produce a land subsidence map, and both compared their results with GPS survey data. However, Jakarta’s land subsidence map requires updating as it has remained unchanged for the past 10 years, and research on land subsidence susceptibility maps was not found in Jakarta.
Therefore, this study’s objective was to generate the updated land subsidence map in Jakarta using the Stanford Method for Persistent Scatterers (StaMPS) with the Sentinel-1 synthetic aperture radar (SAR) dataset from March 2017 to May 2020 in ascending track and March 2017 to April 2020 in a descending track. The mean vertical deformation maps in ascending and descending tracks were compared to validate Jakarta’s land subsidence location. After that, the land subsidence map obtained with this method was used as an inventory map to generate a land subsidence susceptibility map in Jakarta that predicted the areas at risk of land subsidence in the future. Two meta-ensemble machine learning algorithms (adaptive boosting (AdaBoost) and Logit Boost), and two functional machine learning algorithms (logistic regression and multilayer perceptron) were used. The result of the land subsidence susceptibility map produced by these algorithms was evaluated to compare all models’ accuracy and reliability using receiver operating characteristic (ROC) curve analysis.

2. Study Area

Jakarta is Indonesia’s capital city (Figure 1a), located at 6°17’ south (S) and 106°82’ east (E), on the northern coast of western Java. It is considered a lowland area, with an average altitude of ±7 m above sea level [3,47]. In 2019, its population was 10.5 million, with a population growth rate of about 1.19% per year. The population density was 15,900 people per km2, within a total land area of 662.33 km2 [47]. Historically, the population of Jakarta gradually migrated to the other districts and municipalities that are part of the Jakarta Metropolitan Region (JMR), which covers a total area of 5897 km2 and includes Bogor, Depok, Tangerang, and Bekasi [48] (Figure 1b). This migration to Jakarta’s outskirts led to increased urban development that could cause land subsidence in these areas [3]. Our study area covers the JMR; however, for simplicity, we refer to it as Jakarta. The three administrative areas were chosen as study area due to land subsidence reported on those areas from the recent study using GPS surveys, leveling surveys, and InSAR techniques. Training and testing data in Figure 1b used in this study to generate a land subsidence susceptibility map were assessed from the land subsidence inventory map.
The geological and geomorphological area in Jakarta lithologically was dominated by alluvium landform (50.20%) and alluvium fans (19.66%), followed by Serpong form (7.03%) that dominated by fragmented pumice sandstones, some limestone, and andesite. A volcanic pyroclastic flow formed tuff Banten (6.77%) and upper Banten tuff (4.88%), Sandstone unit (4.05%), swamp deposit (1.75%), Cihoe form (1.61%), beach ridge deposit (1.53%), and old alluvium (0.74%). Subang form (0.64%) was dominated by layered claystone lithology with limestone and marl found locally, marine deposit (0.52%), and young volcanic rocks (0.22%); Bojongmanik form (0.22%) was dominated by alternating sandstones and claystone inserted by limestones and coastal deposit (0.13%); Parigi form (0.05%) was dominated by medium limestone, lake (0.01%), and sandstone tuff (0.01%). The domination of alluvium landform has a risk of land subsidence due to the compaction of natural consolidation worsened by a human-made structure [3,9] covering the Jakarta area’s land use. Land use in Jakarta was dominated by settlement area with 45.72% and followed by 39.79% of rice field, 5.94% of dryland agriculture, 4.75% of fish pond, 1.98% of shrub-mixed dryland farms, 0.84% of airport, 0.51% of swamp, 0.35% of estate crop plantation, 0.09% of barren land, 0.03% of secondary mangrove forest, and 0.004% of swamp shrub.
Research on land subsidence in Jakarta has been conducted since 1982 using leveling surveys. There were two main monitoring periods: (1) 1982–1991 and (2) 1991–1997. Land subsidence was found in the period 1982–1991 in three regions with the highest accumulated subsidence value compared to other regions, namely, two regions in the northwestern part (Cengkareng and Kalideres districts) and the third in the northeastern part of Jakarta (Kemayoran-Sunter district); accumulated subsidence was found in Kalideres district with up to 68.5 cm between 1982 and 1991 with an annual rate of subsidence around 8 cm/year. Accumulated land subsidence was found in the Cengkareng district, with up to 60 cm between 1982 and 1991, with an annual subsidence rate of around 7 cm/year. Accumulated land subsidence was found in the Kemayoran-Sunter district, with up to 70 cm with an annual subsidence rate of around 6 cm/year. The second period of monitoring land subsidence using a leveling survey between 1991 and 1997 highlighted one region with the highest accumulated subsidence value than other regions in the Kalideres district, with up to 154.1 cm with an annual subsidence rate of around 23 cm/year [7].
Following this research, land subsidence in Jakarta was monitored using GPS surveys from 1997 to 2005, and two regions affected by land subsidence were reported in two stations: (1) Kwitang district and (2) Pantai Mutiara district. The accumulated land subsidence was found in the Kwitang district, with up to 48 cm with an annual rate of subsidence rate of 5 cm/year; accumulated land subsidence in Pantai Mutiara district was around 50 cm with an annual rate of subsidence around 4.6 cm/year. [6]. Land subsidence in Jakarta measured using leveling and GPS surveys positively correlated with excessive groundwater extraction and sea-level rise [6,7,49].
Land subsidence was also reported using InSAR techniques based on ALOS PALSAR satellite data in 2007–2009 using the SBAS method, with land subsidence found in Pluit district with an annual subsidence rate of 21.6 cm/year. Land subsidence was found in Cengkareng district with an annual subsidence rate around 21.8 cm/year, in Bekasi district with an annual rate of subsidence of around 10.6 cm/year, and in Karawang district with an annual rate of around 16.4 cm/year [1]. Research was also conducted with ALOS PALSAR satellite data within the period 2007–2010 using the GEOS-PSI method, finding land subsidence in the coastal area and lowland area in northwestern Jakarta within Penjaringan and Cengkareng districts with an annual subsidence rate of up to 26 cm/year, with accumulated subsidence rates of up to 86.5 cm between 2007 and 2010. Land subsidence was observed in the Bekasi district with an annual subsidence rate of up to 11.5 cm/year [46].

3. Material and Methods

3.1. SAR Datasets

A land subsidence inventory map for generating the land subsidence susceptibility map in Jakarta was generated using Sentinel-1 SAR C-band data (5.5 cm wavelength) provided by the European Space Agency (ESA). The SAR data were acquired from March 2017 to May 2020 (91 datasets in the ascending track with path number 98 and frame number 1160, with vertical–vertical (VV) polarization) and in the period of March 2017 to April 2020 (89 datasets in the descending track with path number 47 and frame number 614, with vertical–vertical (VV) polarization). The ascending and descending datasets are listed in Table 1 and Table 2, and the reference dates with zero delta day and zero perpendicular baselines from the ascending track are shown on 15 October 2018, whereas those from the descending track are shown on 16 November 2018; both reference dates are shown in bold text in table number 45. A perpendicular baseline graph (generated from Table 1 and Table 2 and shown in Figure 2) was used to visualize the temporal baseline from the reference date.

3.2. Land Subsidence Conditioning Factors

In total, 14 land subsidence conditioning factors consisting of geological, geomorphological, topographical, and hydrological data were chosen as the conditioning factors (Table 3) in this study as they are considered the main factors influencing land subsidence [11,13,14]. Each factor was classified using the quantile method, and the factors with different cell size resolutions were resampled into raster datasets with 30 m cell size from each conditioning factor to standardize each factor’s resolution.
A recent study in Jakarta found that the leading causes of land subsidence in Jakarta are groundwater extraction, the load of buildings and construction, natural consolidation of alluvium soil, and tectonic activities [3]. In this study, groundwater drawdown data were collected from the Groundwater Conservation Center, The Ministry of Energy and Mineral Resources of Indonesia. The observation of groundwater drawdown data is carried out periodically through an automatic water level record (AWLR) system. The annual change in groundwater level from 2019 to 2020 was calculated from 15 monitoring wells, and the map was constructed (Figure 3a) using the inverse distance weighting (IDW) interpolation method from the annual change in groundwater level data. The obtained groundwater drawdown data in this study were insufficient compared to the study area due to the limited monitoring wells over the study area. Although the available data are few, the use of groundwater drawdown data is essential to determine the relationship between land subsidence and groundwater extraction.
The groundwater level can increase due to conditional factors indirectly associated with land subsidence, such as rainfall, distance from a river, and river density, which can recharge the groundwater level [13,50]. Four years of daily rainfall data (from 2017 to 2020) from seven weather stations were acquired from the Indonesian Meteorology, Climatology, and Geophysical Agency of Indonesia, and the annual rainfall intensity was calculated and interpolated using the IDW tool in geographic information system (GIS) software (ArcGIS; ESRI, Redlands, CA, USA) (Figure 3b). However, the availability of rainfall data was limited compared to the study area due to the cost and constraints of acquiring data. Nevertheless, the utilization of the rainfall data in this study is important to correlate the topographical factor related to the infiltration flow affecting the soil’s strength.
Road, river, and fault networks were acquired from the Geospatial Information Agency of Indonesia from the main road map, main river map, and fault map at a scale of 1:250,000 of polygonal-shaped data using the Atlas map of Indonesia. The information on the location of roads, rivers, and faults was used to construct a distance map (Figure 3c, Figure 3d, and Figure 3e, respectively) using the Euclidean distance tool, and the maps were classified using the quantile method to provide suitable classes within a 30 m cell size. The drainage density or river map density was estimated using the kernel density tool (Figure 3f). Distance to the road and land use are related to urban development in Jakarta, which can affect land subsidence [3]. The land-use map (Figure 3g) was acquired from the Ministry of Environment and Forestry of Indonesia that used Landsat data to generate a land cover map. The geological parameters describe the spatial correlation of lithological landform and land subsidence in Jakarta as being caused by the compaction of the alluvium soil landform [3]. The lithology map (Figure 3h) was acquired from the Geological Atlas map of Indonesia from the Ministry of Energy and Mineral Resources, and the polygonal-shaped map was converted into a raster map with 30 m cell size using GIS tools.
The topographical map, which includes elevation, slope, aspect, profile curvature, plan curvature, and topographic wetness index (TWI) (Figure 3i, Figure 3j, Figure 3k, Figure 3l, Figure 3m, and Figure 3n, respectively) data extracted from the digital elevation model (DEM) of the Shuttle Radar Topography Mission (SRTM) [51], was constructed using the basic terrain analysis tools. Elevation influences the hydrological properties and soil moisture, whereby a lower-elevation area possibly gains more precipitation than a higher-elevation area [13,14]. The slope is associated with land subsidence because it affects the infiltration of rainfall (a steeper surface slope decreases infiltration) [13,21], and the aspect is the second derivative of the slope that has a relationship with land subsidence because the slope aspect affects the strength of the soil due to the moisture preservation and the amount of vegetation [11]. Profile curvature is associated with flow speed, sediment, and erosion quantity, while plan curvature is perpendicular to the slope and indirectly affects land subsidence by influencing convergence and divergence of flow across the surface [13]. The TWI defines the degree of deposition of water at a specific site [52]. The topographical factors extracted from the SRTM DEM are widely used as conditioning factors in land subsidence susceptibility mapping [11,12,13,28].
The relationship of land subsidence occurrence with the conditioning factors in each class is described in Table 4. A ratio greater than 1 denotes that the class in the conditioning factor is more correlated with the land subsidence occurrence [11]. The calculation of frequency ratio shows that the land subsidence occurred in areas with groundwater level data between 20.27 and 21.00 and between 23.31 and 34.55 m below ground level, while land subsidence also occurred in areas with more annual rainfall intensity due to the recharge of groundwater level and the usage of groundwater. Areas between 0 and 126 m with roads were more correlated with land subsidence occurrence. The land subsidence areas correlated with fault distance were between 15,944 and 68,975 m from the fault location. There were three drainage density classes correlated with land subsidence occurrence and settlement areas correlated with land subsidence. In terms of lithological factors, alluvium, alluvium fans, beach ridge deposits, and sandstone landforms were considered more correlated with the land subsidence occurrences. There were three classes in the elevation map, four classes in the aspect map, two classes in the slope map, one class in the plan curvature, one class in the profile curvature, and three classes in the topographic wetness index map correlated with land subsidence occurrences.

3.3. Illustration of Methodology

The methods to generate land subsidence susceptibility maps using functional and meta-ensemble algorithms are described below and illustrated in Figure 4.
  • Land subsidence occurrences were identified by exploiting Sentinel-1 SAR datasets from 2017 to 2020 from both ascending and descending tracks using time-series InSAR techniques based on the StaMPS algorithm. The persistent scatterer points from co-registered single master images showing a deformation value were used as the land subsidence inventory map.
  • Preparation of training and testing datasets was conducted by randomly dividing the persistent scatterer (PS) points of time-series InSAR showing a vertical deformation into 50% training data to generate land subsidence susceptibility models and 50% testing data to validate the land subsidence susceptibility map, as done in other studies finding optimal results [28,53]. The distribution of training and test data is shown in Figure 1b.
  • Preparation of land subsidence conditioning factors for spatial correlation analysis was done using the frequency ratio method to find the correlation between each factor and land subsidence occurrence [53]. We used each model’s ratio value and then used as the conditioning factors related to land subsidence occurrences. First, the conditioning factors were classified using quantile methods in GIS tools with a similar environment of 30 m cell size for each factor; then, the number of subsidence occurrences in each class was calculated using the cross-tabulation tool in GIS. Next, we calculated the ratio between the percentage of pixels of each conditioning factor class and the percentage of subsidence occurrence pixels to obtain the FR value as follows:
Frequency   Ratio = %   of   class   conditioning   factor %   of   land   subsidence   occurrence
4.
The conditioning factors consisting of frequency ratio values were used to generate land subsidence susceptibility models using two functional algorithms (logistic regression and multilayer perceptron) and two meta-ensemble algorithms (AdaBoost and LogitBoost).
5.
After all land subsidence susceptibility maps were generated, all maps were validated using the test data prepared before and analyzed using ROC curve analysis.

3.4. StaMPS Processing

StaMPS (Stanford Method for Persistent Scatterers) is an analysis method used to facilitate the generation of time-series deformation images of all terrains, including nonurban areas. The StaMPS algorithm uses the spatial correlation of phase measurements rather than a functional temporal model to identify PS pixels. StaMPS processing does not use a model to describe how the displacement rate varies with time. To identify PS pixels in a single master stack of interferograms, StaMPS uses the phase characteristics from the dominant point scatterer in each area and creates interferograms from SAR images. It also reduces decorrelation [54]. Thus, the StaMPS algorithm can identify and extract the deformation signal from stable pixels in all terrains.
The slave images in the acquired SAR datasets were resampled to perfectly match the master image through a co-registration process before generating an interferogram. The co-registration process was applied to two different images to produce refined SAR images, which were then cropped to focus only on the area of interest. Next, differential InSAR (DInSAR) images were generated by subtracting the topographic InSAR images generated using the interferograms from the topographic phases of the SRTM DEM [55,56]. We used the PS method to measure the displacement of the earth’s surface [57] to generate the time-series deformation map. The main processes generated a single master stack of interferograms and topographic phase removal [33].
The StaMPS algorithm is shown in Figure 5. Multiple images were co-registered to generate a single master image for the ascending track on 15 October 2018 and the descending track on 16 November 2018; co-registered images with the topographic phases removed were subjected to amplitude and phase noise analysis to derive a subset comprising all PS pixels, with weeding performed using a threshold value of 0.4. The wrapped phase of the selected pixels was corrected for the spatially uncorrelated look-angle error. The phase was then unwrapped, and PS outputs were generated. The deformation map from the line of sight (LOS) displacement could be converted into vertical deformation data [58,59] by assuming the horizontal deformation as very small compared to the vertical deformation caused by land subsidence [60,61,62]. Recent studies using GPS and leveling surveys reported that the land subsidence in Jakarta shows a vertical deformation [7], and a vertical deformation pattern was also found in research using InSAR to monitor land subsidence in Jakarta [1,46]; hence, the deformation from the line of sight (LOS) in this study could be assumed as negligible and could be converted directly into vertical deformation value using Equation (2) by dividing the displacement or deformation from the line of sight ( d LOS ) by the cosine of incident angle (θ) from the radar signal. The results of vertical deformation were assigned a negative value from the initial ground-level observation point, indicating that the land subsidence occurred vertically at that point [9].
V =   d LOS cos θ

3.5. AdaBoost

AdaBoost is a machine learning algorithm introduced by Freund and Schapire (1997) [63]. AdaBoost’s classifier uses an adaptive resampling technique that produces a series of individual classifiers to classify training samples accurately. The frequency of variables selected by a weak learner was examined, and the relative importance of the variables could be determined [64]. AdaBoost combines multiple weak learners to derive a single strong learner by repeatedly calling a weak classifier and adjusting the attributed weight to the sample. The new classifier focuses more on the misclassified sample as the misclassified sample is increased compared to the correct sample. The final weight is obtained by adding or subtracting the updated weight in each iteration. The final model is obtained by dividing each final weight by the total adjusted weight [65]. The general form of the AdaBoost algorithm is as follows [66]:
  • Start with weights   w i = 1 N for i   =   1 ,   ,   N ;
  • Repeat this step for  m   =   1 ,   ,   M :
    • Fit the classifier  f m x     1 , 1   using weights  w i with the training data;
    • Compute  err m =   E w 1 ( y f m x ) ,   c m = log 1     err m err m ;
    • Set  w i w i   exp c m 1 ( y f m x ) ,   i   =   1 ,   2 ,   ,   N , and renormalize so that i w i =   1 ;
  • Output the classifier:  sign m = 1 M c m f m ( x ) .

3.6. LogitBoost

LogitBoost is a modified version of the AdaBoost algorithm, introduced by Friedman, Hastie, and Tibshirani [66], where the exponential loss function is replaced with the log-likelihood loss function. This method reduces classification error and is less sensitive to noise [28,67]. LogitBoost can handle multiple class problems and uses a regression scheme as the base learner for classification [64]. The general form of the LogitBoost algorithm is as follows [66]:
  • Start with weights   w i = 1 N   for i   =   1 ,   2 ,   ,   N ,   F x   =   0 , and probability estimates p x i   =   1 2 ;
  • Repeat this step for m   =   1 ,   ,   M :
    • Compute the working response and weights:
      l i   = y i p x i p x i   1     p x i ,
      w i   =   p x i   1     p x i ;
    • Fit the function by weighted least-squares regression of  l i  to  x i using weight  w i ;
    • Update the function as follows:
      f ( x )     f ( x )   +   1 2 f m   ( x ) ,
      p ( x )     e f ( x ) e f ( x ) + e f ( x ) ;
  • Output the classifier:  sign F ( x )   =   sign m = 1 M f m ( x ) .

3.7. Logistic Regression

Logistic regression is a statistical method used to find the best model to describe the correlation between a dependent variable and several independent variables. This method’s advantage is that the variables do not need to be normally distributed [68]. It also offers several ways of selecting the best predictor for use in the P probability model [69,70]. The equations describing the logistic regression are as follows [28,69,70,71]:
f x = logit P = ln P 1 P = c 0 + c 1 x 1 + + c n x n ,
P   =   1 1   +   e f ( x )   =   1 1   +   e ( c 0 + c 1 x 1 + + c n x n ) ,
where f x is a linear combination function called logit P , P is the probability of subsidence occurrence, 1 P is the probability that subsidence will not occur, x 1 ,   x 2 , ,   x n are input variables, c 0 is the model intercept, and c 1 , ,   c n are the approximated coefficients of regression.

3.8. Multilayer Perceptron

The multilayer perceptron is a machine learning algorithm based on the ANN technique that consists of input layers, hidden layers, and output layers [72]. The advantages of the multilayer perceptron algorithm are as follows: the distribution of the training data does not require any assumptions, most inputs are selected during the training process based on the weight adjustment, and the relative importance of the various input measures does not need to be determined [73]. The equation representing the multilayer perceptron for land subsidence classification is as follows [28,72,73]:
m   =   f ( c ) ,
where f(c) is a hidden function that is optimized during the training process for a given network architecture via the adjustable weights, and c   =   c i for i = 1 , , 14 , which is a vector containing 14 land subsidence conditioning factors.

4. Results

4.1. Land Subsidence Inventory Map

The mean vertical deformation maps for Jakarta shown in Figure 6a,b were derived by time-series InSAR analysis based on the StaMPS algorithm. We overlaid a true color red/green/blue (RGB) composite image from Sentinel-2 taken on 28 August 2020. Six areas from both ascending and descending tracks showed high deformation; thus, we plotted a deformation trend for these areas.
The vertical deformation rate in the ascending track at Point P1 (Figure 6c) represents the Kosambi area, which subsided 189.48 mm from 2017 to 2020. Point P2 (Figure 6c) represents the Cengkareng area, which subsided 184.02 mm from 2017 to 2020. Point P3 (Figure 6e) represents the Ciledug area, which subsided up to 155.22 mm from 2017 to 2020. Point P4 (Figure 6e) represents the Penjaringan area, which subsided 148.36 mm from 2017 to 2020. Point P5 (Figure 6g) represents the Bekasi area, which subsided 128.17 mm from 2017 to 2020. Point P6 (Figure 6g) represents the Cikarang area, which subsided 271.84 mm from 2017 to 2020.
The vertical deformation rate in the descending track at Point P1 (Figure 6d) represents the Kosambi area, which subsided 210.07 mm from 2017 to 2020. Point P2 (Figure 6d) represents the Cengkareng area, which subsided 216.19 mm from 2017 to 2020. Point P3 (Figure 6f) represents the Ciledug area, which subsided up to 155.92 mm from 2017 to 2020. Point P4 (Figure 6f) represents the Penjaringan area, which subsided 148.31 mm from 2017 to 2020. Point P5 (Figure 6h) represents the Bekasi area, which subsided 107.19 mm from 2017 to 2020. Point P6 (Figure 6h) represents the Cikarang area, which subsided 257.94 mm from 2017 to 2020.
Figure 6g,h show that the deformation pattern was mostly linear. However, Figure 6c–f, representing the Kosambi, Cengkareng, Ciledug, and Penjaringan areas, show quite periodic subsidence with the standard deviation of the vertical deformation rate being higher than the vertical deformation rate in Figure 6g,h. These results occurred due to the seasonal effect of groundwater extraction, and those areas were surrounded by a residential area that mostly used groundwater as the water source. Meanwhile, P6 is one of the most significant industrial areas in Indonesia. The mean vertical deformation maps from ascending and descending tracks were compared to validate the accuracy of the land subsidence inventory map using the StaMPS algorithm, and the result showed a good correlation in Figure 6i with a coefficient of correlation (R2) up to 0.9584 between ascending and descending tracks.
The persistent scatterer density of both tracks in the east and west areas was relatively low due to them being wetland areas and more vegetated than other areas. The SAR dataset used in this study from Sentinel-1 SAR C-band data with 5.5 cm wavelength could not deeply penetrate beneath the trees. Thus, to overcome that limitation, the interpolation of the persistent scatterer points from the descending track was constructed using kriging interpolation in GIS tools to provide land subsidence information over the study area as shown in Figure 6j [74].

4.2. Land Subsidence Susceptibility Map

The land subsidence susceptibility model’s performance depended on the calculated parameters (Table 5) for optimization used in this study.
A land subsidence susceptibility map was generated using 14 land subsidence conditioning factors, training data from our land subsidence inventory map, and four different algorithms: LogitBoost (Figure 7a), AdaBoost (Figure 7b), logistic regression (Figure 7c), and multilayer perceptron (Figure 7d). Land subsidence susceptibility indices were generated for all unique pixels in the study area. The susceptibility indices were reclassified using the quantile method to identify feature pairs in five susceptibility classes: very low, low, moderate, high, and very high [28,64].
The proportion of very high susceptibility land was quite similar for all four algorithms. However, the map generated by AdaBoost differed from those of the other three methods because the AdaBoost model could only classify susceptibility into four classes due to the limits of the probability range. Therefore, the very low susceptibility class was excluded. The susceptibility class proportions (pixel distributions) are shown in Figure 8. For the Adaboost algorithm, the proportions were 0%, 57.88%, 2.54%, 22.02%, and 17.56% for the very low, low, moderate, high, and very high classes, respectively; the respective values for the LogitBoost algorithm were 33.33%, 27.09%, 7.03%, 16.77%, and 15.78%, and those for the logistic regression algorithm were 32.62%, 13.86%, 16.44%, 18.51%, and 18.57%. Finally, the multilayer perceptron algorithm’s respective values were 40.64%, 18.70%, 14.02%, 13.42%, and 13.23%.
The distribution of pixels in the very high class and high class in Figure 8 showed similar results for each algorithm, resulting in the maps of these classes being quite similar. The very high class was considered to be the land subsidence areas shown in the mean vertical deformation map in Figure 6a,b, with similarity seen because of the training data used in this study being acquired from the land subsidence inventory map with a large spatial resolution, which allowed more effectively defining the land subsidence area in the susceptibility map. The moderate and high classes were considered land subsidence areas in the future, and the very low and low classes were areas with the lowest probability of land subsidence in the future. The difference between the moderate class from AdaBoost and that from other algorithms is that AdaBoost did not consider other areas far from the land subsidence occurrences. Other algorithms showed that residential areas located in alluvium landforms had a reasonable possibility of land subsidence.

4.3. Model Validation

The accuracy of all four algorithms in this study was evaluated by ROC curve analysis [12,17]. ROC curve analysis is a standard way of validating the probability models used to generate land subsidence susceptibility maps, according to the area under the curve (AUC) [22,28]. Higher values indicate more accurate and reliable models. If the AUC, which ranges from 0 to 1, is lower than 0.5, the model is considered unacceptably inaccurate [75].
Land subsidence susceptibility maps produced using functional (AdaBoost, LogitBoost) and meta-ensemble (logistic regression and multilayer perceptron) algorithms were compared. The ROC curves for the four algorithms are shown in Figure 9. The largest AUC of 0.811 was from the AdaBoost algorithm (blue line in Figure 9). The multilayer perceptron algorithm had the next largest AUC (0.800; purple line in Figure 9), followed by the logistic regression (0.794; green line in Figure 9) and LogitBoost algorithms (0.791; red line in Figure 9).
Variables having at least one tie between the positive actual state group and the negative actual state group can bias the results. Since all of the AUC values were higher than 0.5, the land subsidence susceptibility maps produced by all algorithms used in this study are acceptable for predicting high-risk subsidence areas in Jakarta [75]. However, the AdaBoost algorithm map had the best performance.

5. Discussion

5.1. Land Subsidence Inventory Map

A land subsidence inventory map was successfully created through time-series InSAR analysis of Sentinel-1 datasets, using the StaMPS algorithm. The map is displayed as a vertical deformation map in Figure 6a,b for descending and ascending tracks. The comparison was made to validate the mean vertical deformation’s accuracy used as a land subsidence inventory map. A coefficient of correlation of 0.9584 was found. As stated above, there were six high-deformation areas (Figure 6c–h) within ascending and descending tracks, which were equally affected (in terms of subsidence) by increased groundwater usage and existing buildings [3].
Figure 6g,h show areas with largely linear deformation rates, the standard deviations of which were smaller than those of the area shown in Figure 6c–f. The periodic subsidence shown was affected due to the variable climate in Jakarta, which is generally characterized by high rainfall but a dry summer. In the rainy season, the groundwater level of all aquifers beneath Jakarta increases, which manifests as deformation, i.e., an uplift in the ground level due to a change in the underlying material’s thickness [34,76,77,78,79].
As shown by the land subsidence occurrence (Figure 6a,b) and groundwater drawdown (Figure 3a) data, land subsidence for P1–P5, which are between 20.58 and 34.55 m below ground level, correlated with the groundwater level. The correlation between land subsidence and groundwater level might have been stronger if groundwater level data were available for the whole study area. The groundwater level is deepest for P1 and P2, ranging from 22.47 to 34.55 m. P6 has more significant deformation because it contains the largest industrial estate in Indonesia.
Jakarta is mostly situated on young alluvial soil (Figure 3h), which cannot tolerate the maximum compression force of many buildings [9]. Thus, compaction of the unconsolidated alluvial soil occurs and is exacerbated by groundwater extraction (because pore pressure is reduced, leading to further clay compaction) [80,81].

5.2. Land Subsidence Susceptibility Map

The method used to create land subsidence susceptibility maps strongly affects the quality of the mapping. Machine learning techniques are effective [12,13,28,64]. In particular, the method used to generate training and testing data is important. Accurate land subsidence inventory maps can be obtained using InSAR; we combined InSAR and GIS spatial data to produce an accurate land subsidence susceptibility map [12]. Nevertheless, the distribution data of groundwater drawdown and rainfall intensity data from the interpolation process might not represent the whole study area. They could have affected the spatial distribution of the raster map. Access to the monitoring wells and rainfall station outside the Jakarta area could provide more accurate analysis to determine the relationship between those factors and land subsidence.
The land subsidence susceptibility maps of all four algorithms used in this study could be ordered according to the accuracy and time taken to build the model due to the similarities in the conditioning factors, training and testing data, and study area [28]. We used ROC curve analysis to assess the accuracy of the maps. The AUC data showed that the AdaBoost algorithm had the highest susceptibility class predictive accuracy of 81.1%, which was 1.1% higher than the multilayer perceptron algorithm, 1.7% higher than the logistic regression algorithm, and 2% higher than the LogitBoost algorithm. Since the accuracy of the algorithms used was so closely related, we analyzed the time consumption of data preparation needed to build the model. The AdaBoost algorithm needed 218.63 s to build the model, the LogitBoost algorithm needed 174.94 s to build the model, logistic regression needed 16667.96 s or 4.63 h to build the model, and multilayer perceptron needed 47,538.65 s or 13.20 h to build the model, thus identifying the LogitBoost algorithm as the fastest.
However, the LogitBoost algorithm had the lowest prediction accuracy (79.10%). The Adaboost algorithm needed 218.63 s or 3.64 min to build the model, along with having the highest prediction accuracy of 81.10% in the ROC curve analysis. Logistic regression needed 4.63 h to build the model, along with a prediction accuracy of about 79.40%, and the multilayer perceptron needed 13.20 h, along with a prediction accuracy of 80%. The maps generated using LogitBoost, logistic regression, and multilayer perceptron (Figure 7b–d, respectively) showed similar weights and indicated that the settlement area was at a high risk of land subsidence. The land-use map generated in this study indicated a relatively strong correlation between land subsidence and urban development, likely due to excessive groundwater extraction in urban areas [3,6,7]. Our map could be used by local environmental authorities and those in charge of urban development to identify areas with high subsidence risk. The method of combining SAR and GIS spatial data to generate land subsidence susceptibility maps employed in this study could be applied to other regions.

6. Conclusions

We generated a land subsidence inventory map using a time-series InSAR method based on the StaMPS algorithm from Sentinel-1 SAR datasets in ascending and descending tracks. The comparison of both tracks was conducted, finding a coefficient of correlation between the two tracks of 0.9548. The inventory map could be used as a training and testing dataset to create a land subsidence susceptibility map for Jakarta using meta-ensemble (AdaBoost and LogitBoost) and functional (logistic regression and multilayer perceptron) machine learning algorithms. We created a land subsidence inventory map through time-series InSAR analysis of Sentinel-1 SAR datasets using the StaMPS algorithm. The land subsidence susceptibility map produced by the AdaBoost machine learning algorithm had higher accuracy (AUC = 0.811) and only need 3.64 min to build the model compared with the maps created using the other algorithms (multilayer perceptron, AUC = 0.800; logistic regression, AUC = 0.794; LogitBoost, AUC = 0.791). LogitBoost was the fastest algorithm in building the model but had the lowest predictive accuracy. Logistic regression and multilayer perceptron needed 4.63 and 13.20 h, respectively, to build the model. All of the maps showed acceptable accuracy, as the AUC values were all higher than 0.5; thus, they can all be used for analyzing land subsidence susceptibility in Jakarta. Our approach based on time-series InSAR analysis, machine learning, and GIS spatial data yielded reasonable predictions of areas with high risk of land subsidence. Further research could use alternative algorithms and conditioning factors to generate land subsidence susceptibility maps in other regions.

Author Contributions

Conceptualization, C.-W.L., and W.L.H.; methodology, C.-W.L., A.R.A., and W.L.H.; software, A.R.A. and W.L.H.; validation, C.-W.L. and W.L.H.; formal analysis, W.L.H.; investigation, C.-W.L.; resources, C.-W.L.; data curation, C.-W.L. and W.L.H.; writing—original draft preparation, W.L.H.; writing—review and editing, C.-W.L., A.R.A., and W.L.H.; visualization, C.-W.L., A.R.A., and W.L.H.; supervision, C.-W.L.; project administration, C.-W.L.; funding acquisition, C.-W.L. All authors read and agreed to the published version of the manuscript.

Funding

This research was supported by a grant from the National Research Foundation of Korea provided by the government of Korea (No. 2019R1A2C1085686).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Chaussard, E.; Amelung, F.; Abidin, H.; Hong, S.-H. Sinking cities in Indonesia: ALOS PALSAR detects rapid subsidence due to groundwater and gas extraction. Remote Sens. Environ. 2013, 128, 150–161. [Google Scholar] [CrossRef]
  2. Machowski, R.; Rzetala, M.A.; Rzetala, M.; Solarski, M. Geomorphological and Hydrological Effects of Subsidence and Land use Change in Industrial and Urban Areas. Land Degrad. Dev. 2016, 27, 1740–1752. [Google Scholar] [CrossRef]
  3. Abidin, H.Z.; Andreas, H.; Gumilar, I.; Fukuda, Y.; Pohan, Y.E.; Deguchi, T. Land subsidence of Jakarta (Indonesia) and its relation with urban development. Nat. Hazards 2011, 59, 1753–1771. [Google Scholar] [CrossRef]
  4. Takagi, H.; Esteban, M.; Mikami, T.; Fujii, D. Projection of coastal floods in 2050 Jakarta. Urban Clim. 2016, 17, 135–145. [Google Scholar] [CrossRef]
  5. Budiyono, Y.; Aerts, J.C.J.H.; Tollenaar, D.; Ward, P.J. River flood risk in Jakarta under scenarios of future change. Nat. Hazards Earth Syst. Sci. 2016, 16, 757–774. [Google Scholar] [CrossRef] [Green Version]
  6. Abidin, H.Z.; Andreas, H.; Djaja, R.; Darmawan, D.; Gamal, M. Land subsidence characteristics of Jakarta between 1997 and 2005, as estimated using GPS surveys. GPS Solut. 2008, 12, 23–32. [Google Scholar] [CrossRef]
  7. Abidin, H.Z.; Djaja, R.; Darmawan, D.; Hadi, S.; Akbar, A.; Rajiyowiryono, H.; Sudibyo, Y.; Meilano, I.; Kasuma, M.A.; Kahar, J.; et al. Land subsidence of Jakarta (Indonesia) and its geodetic monitoring system. Nat. Hazards 2001, 23, 365–387. [Google Scholar] [CrossRef]
  8. Notti, D.; Mateos, R.M.; Monserrat, O.; Devanthéry, N.; Peinado, T.; Roldán, F.J.; Fernández-Chacón, F.; Galve, J.P.; Lamas, F.; Azañón, J.M. Lithological control of land subsidence induced by groundwater withdrawal in new urban AREAS (Granada Basin, SE Spain). Multiband DInSAR monitoring. Hydrol. Process. 2016, 30, 2317–2331. [Google Scholar] [CrossRef]
  9. Yastika, P.E.; Shimizu, N.; Abidin, H.Z. Monitoring of long-term land subsidence from 2003 to 2017 in coastal area of Semarang, Indonesia by SBAS DInSAR analyses using Envisat-ASAR, ALOS-PALSAR, and Sentinel-1A SAR data. Adv. Space Res. 2019, 63, 1719–1736. [Google Scholar] [CrossRef]
  10. Fiaschi, S.; Tessitore, S.; Bonì, R.; Di Martire, D.; Achilli, V.; Borgstrom, S.; Ibrahim, A.; Floris, M.; Meisina, C.; Ramondini, M.; et al. From ERS-1/2 to Sentinel-1: Two decades of subsidence monitored through A-DInSAR techniques in the Ravenna area (Italy). GIScience Remote Sens. 2017, 54, 305–328. [Google Scholar] [CrossRef]
  11. Pradhan, B.; Abokharima, M.H.; Jebur, M.N.; Tehrany, M.S. Land subsidence susceptibility mapping at Kinta Valley (Malaysia) using the evidential belief function model in GIS. Nat. Hazards 2014, 73, 1019–1042. [Google Scholar] [CrossRef]
  12. Bianchini, S.; Solari, L.; Del Soldato, M.; Raspini, F.; Montalti, R.; Ciampalini, A.; Casagli, N. Ground Subsidence Susceptibility (GSS) mapping in Grosseto plain (Tuscany, Italy) based on satellite InSAR data using frequency ratio and fuzzy logic. Remote Sens. 2019, 11, 2015. [Google Scholar] [CrossRef] [Green Version]
  13. Arabameri, A.; Saha, S.; Roy, J.; Tiefenbacher, J.P.; Cerda, A.; Biggs, T.; Pradhan, B.; Thi Ngo, P.T.; Collins, A.L. A novel ensemble computational intelligence approach for the spatial prediction of land subsidence susceptibility. Sci. Total Environ. 2020, 726, 138595. [Google Scholar] [CrossRef]
  14. Abdollahi, S.; Pourghasemi, H.R.; Ghanbarian, G.A.; Safaeian, R. Prioritization of effective factors in the occurrence of land subsidence and its susceptibility mapping using an SVM model and their different kernel functions. Bull. Eng. Geol. Environ. 2019, 78, 4017–4034. [Google Scholar] [CrossRef]
  15. Oh, H.J.; Lee, S. Integration of ground subsidence hazard maps of abandoned coal mines in Samcheok, Korea. Int. J. Coal Geol. 2011, 86, 58–72. [Google Scholar] [CrossRef]
  16. Oh, H.J.; Lee, S. Assessment of ground subsidence using GIS and the weights-of-evidence model. Eng. Geol. 2010, 115, 36–48. [Google Scholar] [CrossRef]
  17. Pradhan, B.; Oh, H.-J.; Buchroithner, M. Weights-of-evidence model applied to landslide susceptibility mapping in a tropical hilly area. Geomat. Nat. Hazards Risk 2010, 1, 199–223. [Google Scholar] [CrossRef]
  18. Kim, K.D.; Lee, S.; Oh, H.J.; Choi, J.K.; Won, J.S. Assessment of ground subsidence hazard near an abandoned underground coal mine using GIS. Environ. Geol. 2006, 50, 1183–1191. [Google Scholar] [CrossRef]
  19. Hu, B.; Zhou, J.; Wang, J.; Chen, Z.; Wang, D.; Xu, S. Risk assessment of land subsidence at Tianjin coastal area in China. Environ. Earth Sci. 2009, 59, 269–276. [Google Scholar] [CrossRef]
  20. Zhi-xiang, T.; Pei-xian, L.; Li-li, Y.; Ka-zhong, D. Study of the method to calculate subsidence coefficient based on SVM. Procedia Earth Planet. Sci. 2009, 1, 970–976. [Google Scholar] [CrossRef] [Green Version]
  21. Rahmati, O.; Falah, F.; Naghibi, S.A.; Biggs, T.; Soltani, M.; Deo, R.C.; Cerdà, A.; Mohammadi, F.; Tien Bui, D. Land subsidence modelling using tree-based machine learning algorithms. Sci. Total Environ. 2019, 672, 239–252. [Google Scholar] [CrossRef]
  22. Tien Bui, D.; Shahabi, H.; Shirzadi, A.; Chapi, K.; Pradhan, B.; Chen, W.; Khosravi, K.; Panahi, M.; Bin Ahmad, B.; Saro, L. Land Subsidence Susceptibility Mapping in South Korea Using Machine Learning Algorithms. Sensors 2018, 18, 2464. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Choi, J.K.; Kim, K.D.; Lee, S.; Won, J.S. Application of a fuzzy operator to susceptibility estimations of coal mine subsidence in Taebaek City, Korea. Environ. Earth Sci. 2010, 59, 1009–1022. [Google Scholar] [CrossRef]
  24. Jaafari, A.; Panahi, M.; Pham, B.T.; Shahabi, H.; Bui, D.T.; Rezaie, F.; Lee, S. Meta optimization of an adaptive neuro-fuzzy inference system with grey wolf optimizer and biogeography-based optimization algorithms for spatial prediction of landslide susceptibility. Catena 2019, 175, 430–445. [Google Scholar] [CrossRef]
  25. Park, I.; Choi, J.; Jin Lee, M.; Lee, S. Application of an adaptive neuro-fuzzy inference system to ground subsidence hazard mapping. Comput. Geosci. 2012, 48, 228–238. [Google Scholar] [CrossRef]
  26. Lee, S.; Park, I.; Choi, J.K. Spatial prediction of ground subsidence susceptibility using an artificial neural network. Environ. Manag. 2012, 49, 347–358. [Google Scholar] [CrossRef]
  27. Ambrožič, T.; Turk, G. Prediction of subsidence due to underground mining by artificial neural networks. Comput. Geosci. 2003, 29, 627–637. [Google Scholar] [CrossRef] [Green Version]
  28. Oh, H.-J.; Syifa, M.; Lee, C.-W.; Lee, S. Land Subsidence Susceptibility Mapping Using Bayesian, Functional, and Meta-Ensemble Machine Learning Models. Appl. Sci. 2019, 9, 1248. [Google Scholar] [CrossRef] [Green Version]
  29. Ferretti, A.; Prati, C.; Rocca, F. Permanent scatterers in SAR interferometry. IEEE Trans. Geosci. Remote Sens. 2001, 39, 8–20. [Google Scholar] [CrossRef]
  30. Ferretti, A.; Savio, G.; Barzaghi, R.; Borghi, A.; Musazzi, S.; Novali, F.; Prati, C.; Rocca, F. Submillimeter accuracy of InSAR time series: Experimental validation. IEEE Trans. Geosci. Remote Sens. 2007, 45, 1142–1153. [Google Scholar] [CrossRef]
  31. Hooper, A.; Bekaert, D.; Spaans, K.; Arikan, M. Recent advances in SAR interferometry time series analysis for measuring crustal deformation. Tectonophysics 2012, 514–517, 1–13. [Google Scholar] [CrossRef]
  32. Hooper, A. A multi-temporal InSAR method incorporating both persistent scatterer and small baseline approaches. Geophys. Res. Lett. 2008, 35. [Google Scholar] [CrossRef] [Green Version]
  33. Osmanoğlu, B.; Sunar, F.; Wdowinski, S.; Cabral-Cano, E. Time series analysis of InSAR data: Methods and trends. ISPRS J. Photogramm. Remote Sens. 2016, 115, 90–102. [Google Scholar] [CrossRef]
  34. Osmanoǧlu, B.; Dixon, T.H.; Wdowinski, S.; Cabral-Cano, E.; Jiang, Y. Mexico City subsidence observed with persistent scatterer InSAR. Int. J. Appl. Earth Obs. Geoinf. 2011, 13, 1–12. [Google Scholar] [CrossRef]
  35. Haji-Aghajany, S.; Amerian, Y. Atmospheric phase screen estimation for land subsidence evaluation by InSAR time series analysis in Kurdistan, Iran. J. Atmos. Solar Terr. Phys. 2020, 205, 105314. [Google Scholar] [CrossRef]
  36. Yang, C.; Zhang, F.; Liu, R.; Hou, J.; Zhang, Q.; Zhao, C. Ground deformation and fissure activity of the Yuncheng Basin (China) revealed by multiband time series InSAR. Adv. Space Res. 2020, 66. [Google Scholar] [CrossRef]
  37. Tosi, L.; Da Lio, C.; Teatini, P.; Strozzi, T. Land subsidence in coastal environments: Knowledge advance in the Venice coastland by TerraSAR-X PSI. Remote Sens. 2018, 10, 1191. [Google Scholar] [CrossRef] [Green Version]
  38. Kim, S.W.; Dixon, T.; Amelung, F.; Wdowinski, S. A time-series deformation analysis from TerraSAR-X SAR data over New Orleans, USA. In Proceedings of the 2011 3rd International Asia-Pacific Conference on Synthetic Aperture Radar, APSAR 2011, Seoul, Korea, 26–30 September 2011; pp. 832–833. [Google Scholar]
  39. Zhao, Q.; Ma, G.; Wang, Q.; Yang, T.; Liu, M.; Gao, W.; Falabella, F.; Mastro, P.; Pepe, A. Generation of long-term InSAR ground displacement time-series through a novel multi-sensor data merging technique: The case study of the Shanghai coastal area. ISPRS J. Photogramm. Remote Sens. 2019, 154, 10–27. [Google Scholar] [CrossRef]
  40. Imakiire, T.; Koarai, M. Wide-area land subsidence caused by “the 2011 off the Pacific Coast of Tohoku Earthquake”. Soils Found. 2012, 52, 842–855. [Google Scholar] [CrossRef] [Green Version]
  41. Hong, S.-J.; Baek, W.-K.; Jung, H.-S. Mapping Precise Two-dimensional Surface Deformation on Kilauea Volcano, Hawaii using ALOS2 PALSAR2 Spotlight SAR Interferometry. Korean J. Remote Sens. 2019, 35, 1235–1249. [Google Scholar] [CrossRef]
  42. Jo, M.; Osmanoglu, B.; Jung, H.-S. Detecting Surface Changes Triggered by Recent Volcanic Activities at Kīlauea, Hawai’i, by using the SAR Interferometric Technique: Preliminary Report. Korean J. Remote Sens. 2018, 34, 1545–1553. [Google Scholar] [CrossRef]
  43. Liu, P.; Li, Z.; Hoey, T.; Kincal, C.; Zhang, J.; Zeng, Q.; Muller, J.P. Using advanced inSAR time series techniques to monitor landslide movements in Badong of the Three Gorges region, China. Int. J. Appl. Earth Obs. Geoinf. 2012, 21, 253–264. [Google Scholar] [CrossRef]
  44. Samsonov, S.; d’Oreye, N.; Smets, B. Ground deformation associated with post-mining activity at the French-German border revealed by novel InSAR time series method. Int. J. Appl. Earth Obs. Geoinf. 2013, 23, 142–154. [Google Scholar] [CrossRef]
  45. Jung, H.C.; Kim, S.W.; Jung, H.S.; Min, K.D.; Won, J.S. Satellite observation of coal mining subsidence by persistent scatterer analysis. Eng. Geol. 2007, 92, 1–13. [Google Scholar] [CrossRef]
  46. Ng, A.H.M.; Ge, L.; Li, X.; Abidin, H.Z.; Andreas, H.; Zhang, K. Mapping land subsidence in Jakarta, Indonesia using persistent scatterer interferometry (PSI) technique with ALOS PALSAR. Int. J. Appl. Earth Obs. Geoinf. 2012, 18, 232–242. [Google Scholar] [CrossRef]
  47. BPS-Statistics of DKI Jakarta Province. DKI Jakarta Province in Figures; Division of Integration Processing and Statistics Dissemination BPS-Statistics of DKI Jakarta Province, Ed.; BPS-Statistics of DKI Jakarta Province: Jakarta, Indonesia, 2020; ISBN 978-602-0922-38-6. [Google Scholar]
  48. Hudalah, D.; Firman, T. Beyond property: Industrial estates and post-suburban transformation in Jakarta Metropolitan Region. Cities 2012, 29, 40–48. [Google Scholar] [CrossRef]
  49. Abidin, H.Z.; Andreas, H.; Gamal, M.; Gumilar, I.; Napitupulu, M.; Fukuda, Y.; Deguchi, T.; Maruyama, Y. Edi Riawan Land subsidence characteristics of the Jakarta basin (Indonesia) and its relation with groundwater extraction and sea level rise. Groundw. Response Chang. Clim. 2010, 113–130. [Google Scholar] [CrossRef]
  50. Arabameri, A.; Lee, S.; Tiefenbacher, J.P.; Ngo, P.T.T. Novel Ensemble of MCDM-Artificial Intelligence Techniques for Groundwater-Potential Mapping in Arid and Semi-Arid Regions (Iran). Remote Sens. 2020, 12, 490. [Google Scholar] [CrossRef] [Green Version]
  51. Farr, T.G.; Rosen, P.A.; Caro, E.; Crippen, R.; Duren, R.; Hensley, S.; Kobrick, M.; Paller, M.; Rodriguez, E.; Roth, L.; et al. The Shuttle Radar Topography Mission. Rev. Geophys. 2007, 45, RG2004. [Google Scholar] [CrossRef] [Green Version]
  52. Pourghasemi, H.R.; Beheshtirad, M. Assessment of a data-driven evidential belief function model and GIS for groundwater potential mapping in the Koohrang Watershed, Iran. Geocarto Int. 2015, 30, 662–685. [Google Scholar] [CrossRef]
  53. Lee, S.; Park, I. Application of decision tree model for the ground subsidence hazard mapping near abandoned underground coal mines. J. Environ. Manag. 2013, 127, 166–176. [Google Scholar] [CrossRef]
  54. Hooper, A. Persistent scatter radar interferometry for crustal deformation studies and modeling of volcanic deformation. Ph.D. Thesis, Stanford University, Stanford, CA, USA, 2006. [Google Scholar]
  55. Wegnüller, U.; Werner, C.; Strozzi, T.; Wiesmann, A.; Frey, O.; Santoro, M. Sentinel-1 Support in the GAMMA Software. Procedia Comput. Sci. 2016, 100, 1305–1312. [Google Scholar] [CrossRef] [Green Version]
  56. Werner, C.; Wegmüller, U.; Strozzi, T.; Wiesmann, A. GAMMA SAR and interferometric processing software. In Proceedings of the ERS—ENVISAT Symposium, Gothenburg, Sweden, 16–20 October 2000. [Google Scholar]
  57. Crosetto, M.; Monserrat, O.; Cuevas-González, M.; Devanthéry, N.; Crippa, B. Persistent Scatterer Interferometry: A review. ISPRS J. Photogramm. Remote Sens. 2016, 115, 78–89. [Google Scholar] [CrossRef] [Green Version]
  58. Hooper, A.; Segall, P.; Zebker, H. Persistent scatterer interferometric synthetic aperture radar for crustal deformation analysis, with application to Volcán Alcedo, Galápagos. J. Geophys. Res. Solid Earth 2007, 112, 1–21. [Google Scholar] [CrossRef] [Green Version]
  59. Sousa, J.J.; Hooper, A.J.; Hanssen, R.F.; Bastos, L.C.; Ruiz, A.M. Persistent Scatterer InSAR: A comparison of methodologies based on a model of temporal deformation vs. spatial correlation selection criteria. Remote Sens. Environ. 2011, 115, 2652–2663. [Google Scholar] [CrossRef]
  60. Pepe, A.; Bonano, M.; Zhao, Q.; Yang, T.; Wang, H. The Use of C-/X-Band Time-Gapped SAR Data and Geotechnical Models for the Study of Shanghai’s Ocean-Reclaimed Lands through the SBAS-DInSAR Technique. Remote Sens. 2016, 8, 911. [Google Scholar] [CrossRef] [Green Version]
  61. Floris, M.; Fontana, A.; Tessari, G.; Mulè, M. Subsidence Zonation Through Satellite Interferometry in Coastal Plain Environments of NE Italy: A Possible Tool for Geological and Geomorphological Mapping in Urban Areas. Remote Sens. 2019, 11, 165. [Google Scholar] [CrossRef] [Green Version]
  62. Ren, H.; Feng, X. Calculating vertical deformation using a single InSAR pair based on singular value decomposition in mining areas. Int. J. Appl. Earth Obs. Geoinf. 2020, 92, 102115. [Google Scholar] [CrossRef]
  63. Freund, Y.; Schapire, R.E. A Decision-Theoretic Generalization of On-Line Learning and an Application to Boosting. J. Comput. Syst. Sci. 1997, 55, 119–139. [Google Scholar] [CrossRef] [Green Version]
  64. Kadavi, P.; Lee, C.-W.; Lee, S. Application of Ensemble-Based Machine Learning Models to Landslide Susceptibility Mapping. Remote Sens. 2018, 10, 1252. [Google Scholar] [CrossRef] [Green Version]
  65. Sprenger, M.; Schemm, S.; Oechslin, R.; Jenkner, J. Nowcasting foehn wind events using the AdaBoost machine learning algorithm. Weather Forecast. 2017, 32, 1079–1099. [Google Scholar] [CrossRef]
  66. Friedman, J.; Hastie, T.; Tibshirani, R. Additive logistic regression: A statistical view of boosting. Ann. Stat. 2000, 28, 337–407. [Google Scholar] [CrossRef]
  67. Zhang, G.; Fang, B. LogitBoost classifier for discriminating thermophilic and mesophilic proteins. J. Biotechnol. 2007, 127, 417–424. [Google Scholar] [CrossRef]
  68. Lee, S. Application of logistic regression model and its validation for landslide susceptibility mapping using GIS and remote sensing data. Int. J. Remote Sens. 2005, 26, 1477–1491. [Google Scholar] [CrossRef]
  69. Erener, A.; Mutlu, A.; Sebnem Düzgün, H. A comparative study for landslide susceptibility mapping using GIS-based multi-criteria decision analysis (MCDA), logistic regression (LR) and association rule mining (ARM). Eng. Geol. 2016, 203, 45–55. [Google Scholar] [CrossRef]
  70. Ozdemir, A.; Altural, T. A comparative study of frequency ratio, weights of evidence and logistic regression methods for landslide susceptibility mapping: Sultan mountains, SW Turkey. J. Asian Earth Sci. 2013, 64, 180–197. [Google Scholar] [CrossRef]
  71. David, W.; Hosmer, J.; Lemeshow, S.; Sturdivant, R.X. Applied Logistic Regression, 3rd ed.; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2013. [Google Scholar]
  72. Pham, B.T.; Tien Bui, D.; Prakash, I.; Dholakia, M.B. Hybrid integration of Multilayer Perceptron Neural Networks and machine learning ensembles for landslide susceptibility assessment at Himalayan area (India) using GIS. Catena 2017, 149, 52–63. [Google Scholar] [CrossRef]
  73. Gardner, M.W.; Dorling, S.R. Artificial neural networks (the multilayer perceptron)—A review of applications in the atmospheric sciences. Atmos. Environ. 1998, 32, 2627–2636. [Google Scholar] [CrossRef]
  74. Kim, S.-W.; Wdowinski, S.; Dixon, T.H.; Amelung, F.; Won, J.-S.; Kim, J.W. InSAR-based mapping of surface subsidence in Mokpo City, Korea, using JERS-1 and ENVISAT SAR data. Earth Planets Space 2008, 60, 453–461. [Google Scholar] [CrossRef] [Green Version]
  75. Fawcett, T. An introduction to ROC analysis. Pattern Recognit. Lett. 2006, 27, 861–874. [Google Scholar] [CrossRef]
  76. Van Leijen, F.J. Persistent Scatterer Interferometry based on geodetic estimation theory. Ph.D. Thesis, Technische Universiteit Delft, Delft, The Netherlands, 2014. [Google Scholar]
  77. Colesanti, C.; Ferretti, A.; Novali, F.; Prati, C.; Rocca, F. SAR monitoring of progressive and seasonal ground deformation using the permanent scatterers technique. IEEE Trans. Geosci. Remote Sens. 2003, 41, 1685–1701. [Google Scholar] [CrossRef] [Green Version]
  78. Smith, R.G.; Knight, R.; Chen, J.; Reeves, J.A.; Zebker, H.A.; Farr, T.; Liu, Z. Estimating the permanent loss of groundwater storage in the southern San Joaquin Valley, California. Water Resour. Res. 2017, 53, 2133–2148. [Google Scholar] [CrossRef]
  79. Reeves, J.A.; Knight, R.; Zebker, H.A.; Kitanidis, P.K.; Schreüder, W.A. Estimating temporal changes in hydraulic head using InSAR data in the San Luis Valley, Colorado. Water Resour. Res. 2014, 50, 4459–4473. [Google Scholar] [CrossRef]
  80. Erban, L.E.; Gorelick, S.M.; Zebker, H.A. Groundwater extraction, land subsidence, and sea-level rise in the Mekong Delta, Vietnam. Environ. Res. Lett. 2014, 9. [Google Scholar] [CrossRef]
  81. Erban, L.E.; Gorelick, S.M.; Zebker, H.A.; Fendorf, S. Release of arsenic to deep groundwater in the Mekong Delta, Vietnam, linked to pumping-induced land subsidence. Proc. Natl. Acad. Sci. USA 2013, 110, 13751–13756. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. (a) Study area location in Indonesia and (b) training and testing data sets from land subsidence location within Jakarta Metropolitan Region (JMR) depicted from Sentinel-2 satellite imagery taken on 28 August 2020, divided into three administrative areas: Jakarta, Tangerang, and Bekasi.
Figure 1. (a) Study area location in Indonesia and (b) training and testing data sets from land subsidence location within Jakarta Metropolitan Region (JMR) depicted from Sentinel-2 satellite imagery taken on 28 August 2020, divided into three administrative areas: Jakarta, Tangerang, and Bekasi.
Remotesensing 12 03627 g001
Figure 2. Perpendicular baseline graph from (a) ascending track and (b) descending track.
Figure 2. Perpendicular baseline graph from (a) ascending track and (b) descending track.
Remotesensing 12 03627 g002
Figure 3. Land subsidence conditioning factors: (a) groundwater drawdown; (b) rainfall intensity; (c) distance to roads; (d) distance to rivers; (e) distance to faults; (f) drainage density; (g) land use; (h) lithology; (i) elevation; (j) slope; (k) aspect; (l) profile curvature; (m) plan curvature; (n) topographic wetness index.
Figure 3. Land subsidence conditioning factors: (a) groundwater drawdown; (b) rainfall intensity; (c) distance to roads; (d) distance to rivers; (e) distance to faults; (f) drainage density; (g) land use; (h) lithology; (i) elevation; (j) slope; (k) aspect; (l) profile curvature; (m) plan curvature; (n) topographic wetness index.
Remotesensing 12 03627 g003aRemotesensing 12 03627 g003b
Figure 4. The study workflow.
Figure 4. The study workflow.
Remotesensing 12 03627 g004
Figure 5. Flowchart of Stanford Method for Persistent Scatterers (StaMPS) Processing.
Figure 5. Flowchart of Stanford Method for Persistent Scatterers (StaMPS) Processing.
Remotesensing 12 03627 g005
Figure 6. Mean vertical deformation map for Jakarta depicted using Sentinel-2 image in (a) ascending tracks and (b) descending tracks; the vertical deformation rate at P1 (Kosambi) and P2 (Cengkareng) in (c) ascending and (d) descending tracks; the vertical deformation rate at P3 (Ciledug) and P4 (Penjaringan) in ascending (e) and descending (f) tracks; the vertical deformation at P5 (Bekasi) and P4 (Cikarang) in ascending (g) and descending tracks (h). (i) The comparison of mean vertical deformation between ascending and descending tracks. (j) Kriging interpolation of time-series deformation map from mean vertical deformation map of descending track, resulting in the extension of the land subsidence inventory map; the blue area is considered as the nonoccurrence area of land subsidence.
Figure 6. Mean vertical deformation map for Jakarta depicted using Sentinel-2 image in (a) ascending tracks and (b) descending tracks; the vertical deformation rate at P1 (Kosambi) and P2 (Cengkareng) in (c) ascending and (d) descending tracks; the vertical deformation rate at P3 (Ciledug) and P4 (Penjaringan) in ascending (e) and descending (f) tracks; the vertical deformation at P5 (Bekasi) and P4 (Cikarang) in ascending (g) and descending tracks (h). (i) The comparison of mean vertical deformation between ascending and descending tracks. (j) Kriging interpolation of time-series deformation map from mean vertical deformation map of descending track, resulting in the extension of the land subsidence inventory map; the blue area is considered as the nonoccurrence area of land subsidence.
Remotesensing 12 03627 g006aRemotesensing 12 03627 g006b
Figure 7. Land subsidence susceptibility map for Jakarta, generated using four different algorithms: (a) AdaBoost, (b) LogitBoost, (c) logistic regression, and (d) multilayer perceptron.
Figure 7. Land subsidence susceptibility map for Jakarta, generated using four different algorithms: (a) AdaBoost, (b) LogitBoost, (c) logistic regression, and (d) multilayer perceptron.
Remotesensing 12 03627 g007
Figure 8. Proportions of susceptibility classes for land subsidence susceptibility maps generated using four different machine learning algorithms.
Figure 8. Proportions of susceptibility classes for land subsidence susceptibility maps generated using four different machine learning algorithms.
Remotesensing 12 03627 g008
Figure 9. Receiver operating characteristic (ROC) curves for the land subsidence susceptibility maps produced by functional (logistic regression and multilayer perceptron) and meta-ensemble (AdaBoost and LogitBoost) algorithms.
Figure 9. Receiver operating characteristic (ROC) curves for the land subsidence susceptibility maps produced by functional (logistic regression and multilayer perceptron) and meta-ensemble (AdaBoost and LogitBoost) algorithms.
Remotesensing 12 03627 g009
Table 1. Acquisition dates (format ddmmyyyy) of the Sentinel-1 synthetic aperture radar (SAR) datasets in ascending track. Delta days = number of days between each acquisition date. B⊥ = perpendicular baseline. The reference dates in ascending tracks are shown in bold text.
Table 1. Acquisition dates (format ddmmyyyy) of the Sentinel-1 synthetic aperture radar (SAR) datasets in ascending track. Delta days = number of days between each acquisition date. B⊥ = perpendicular baseline. The reference dates in ascending tracks are shown in bold text.
No.Acquisition Date
(ddmmyyyy)
DeltaDaysB
(m)
No.AcquisitionDate
(ddmmyyyy)
DeltaDaysB
(m)
No.Acquisition Date
(ddmmyyyy)
DeltaDaysB
(m)
118032017−576773230042018−16864620705201920493
230032017−564653312052018−156−2631905201921616
311042017−55233424052018−14428643105201922820
423042017−540203505062018−1321276512062019240166
505052017−528−253617062018−1201036606072019264104
617052017−516173711072018−9695671807201927644
729052017−5041103823072018−8462683007201928890
810062017−492213904082018−72986911082019300−9
922062017−480214016082018−607570230820193121
1004072017−4681124128082018−4861710409201932446
1109082017−432544209092018−36607216092019336106
1221082017−420914321092018−24557328092019348−14
1302092017−408504403102018−121157410102019360−110
1414092017−396224515102018007522102019372−125
1526092017−3844346271020181253760311201938419
1608102017−37248470811201824857715112019396−2
1720102017−3607248201120183685782711201940838
1801112017−3484349021220184885790912201942098
1913112017−33691501412201860142802112201943287
2025112017−324305126122018726810201202044492
2107122017−31214052070120198474821401202045640
2219122017−3006053190120199646832601202046822
2331122017−288135543101201910840840702202048036
2412012018−276815512022019120103851902202049274
2524012018−26478562402201913229860203202050493
2605022018−252375708032019144−26871403202051631
2717022018−24015582003201915626882603202052841
2801032018−2283059010420191688890704202054040
2913032018−216101601304201918075901904202055214
3006042018−1921396125042019192−31910105202056467
3118042018−180126
Table 2. Acquisition dates (format ddmmyyy) of the Sentinel-1 Sentinel-1 synthetic aperture radar SAR datasets in descending track. Delta days = number of days between each acquisition date. B⊥ = perpendicular baseline. The reference dates in descending tracks are shown in bold text.
Table 2. Acquisition dates (format ddmmyyy) of the Sentinel-1 Sentinel-1 synthetic aperture radar SAR datasets in descending track. Delta days = number of days between each acquisition date. B⊥ = perpendicular baseline. The reference dates in descending tracks are shown in bold text.
No.Acquisition Date
(ddmmyyyy)
DeltaDaysB
(m)
No.AcquisitionDate
(ddmmyyyy)
DeltaDaysB
(m)
No.Acquisition Date
(ddmmyyyy)
DeltaDaysB
(m)
126032017−600223102042018−228−20612705201919224
207042017−588763214042018−2164062080620192040
319042017−576833326042018−20406320062019216−7
401052017−564603408052018−19226640207201922816
513052017−55203520052018−1800651407201924029
606062017−52803613062018−15626660708201926487
718062017−516483725062018−14426671908201927684
830062017−50403831072018−108−426831082019288−5
912072017−49203924082018−84−169120920193000
1024072017−480−104017092018−6007024092019312117
1105082017−468414129092018−4830710610201932494
1217082017−45604211102018−3651721810201933616
1329082017−444574323102018−240733010201934819
1410092017−43204404112018−120741111201936045
1522092017−420−445161120180075231120193720
1604102017−4080462811201812−1876051220193845
1716102017−396164710122018243771712201939654
1828102017−3849648221220183611782912201940863
1909112017−37204903012019480791001202042035
2021112017−3606250150120196024802201202043225
2103122017−348−25127012019720810302202044411
2215122017−336−235208022019841821502202045624
2327122017−324−50532002201996348327022020468−8
2408012018−31205404032019108101841003202048091
2520012018−3000551603201912023852203202049272
2601022018−2881456280320191320860304202050433
2713022018−276395709042019144−1287150420205160
2825022018−2640582104201915613388270420205280
2909032018−252−3659030520191688289090520205400
3021032018−240−68601505201918042
Table 3. Land subsidence conditioning factors in the study area. DEM, digital elevation model; SRTM, Shuttle Radar Topography Mission.
Table 3. Land subsidence conditioning factors in the study area. DEM, digital elevation model; SRTM, Shuttle Radar Topography Mission.
CategoryFactorSource
Hydrological factorsGroundwater drawdownGroundwater Conservation Center of Indonesia, The Ministry of Energy and Mineral Resources
Hydrological factorsRainfall intensityMeteorology, Climatology, and Geophysical Agency of Indonesia
Land cover factorsRoad networkGeospatial Information Agency of Indonesia
Hydrological factorsRiver networkGeospatial Information Agency of Indonesia
Geological factorsFaultsGeospatial Information Agency of Indonesia
Land cover factorsLand useThe Ministry of Environment and Forestry of Indonesia
Geological factorsLithologyThe Ministry of Energy and Mineral Resources
Topographical factorsElevationDEM SRTM 1 Arc-Second Global
Topographical factorsSlopeDEM SRTM 1 Arc-Second Global
Topographical factorsAspectDEM SRTM 1 Arc-Second Global
Geomorphological factorsProfile curvatureDEM SRTM 1 Arc-Second Global
Geomorphological factorsPlan curvatureDEM SRTM 1 Arc-Second Global
Hydrological factorsTopographic wetness indexDEM SRTM 1 Arc-Second Global
Table 4. Relationship between land subsidence occurrence and conditioning factors using frequency ratio (FR) model.
Table 4. Relationship between land subsidence occurrence and conditioning factors using frequency ratio (FR) model.
No.Conditioning FactorClass/CategoryRatio each ClassRatio of OccurrenceFR
1Groundwater drawdown
(m below ground level)
7.77–20.270.18310.12010.6561
20.27–21.000.20210.29531.4616
21.00–21.840.23610.21270.9009
21.84–23.310.19140.12820.6697
23.31–34.550.18740.24371.3003
2Rainfall intensity map (mm/year)1,549–1,7810.19990.10090.5046
1,781–1,8740.19450.16350.8404
1,874–1,9080.19770.22411.1338
1,908–1,9750.20900.30691.4680
1,975–2,1240.19890.20471.0290
3Distance to road map (m)0–1260.21140.22011.0412
126–3280.19780.19640.9931
328–6320.19720.19430.9853
632–1,1630.19680.19640.9979
1,163–6,4510.19680.19280.9795
4Distance to river map (m)0–3400.20090.22621.1262
340–1,0200.19990.23691.1848
1,020–2,2540.19980.20301.0159
2,254–4,4650.19970.14220.7120
4,465–10,8450.19970.19170.9601
5Distance to fault map (m)0–15,9440.20000.03860.1929
15,944–29,1150.20000.20291.0145
29,155–53,3780.20000.33751.6874
53,378–68,9750.20000.28451.4223
68,975–88,3860.20000.13660.6829
6Drainage density (km/km2)00.23310.24081.0331
0–60.19170.16550.8632
6–160.19170.20711.0801
16–480.19170.21011.0960
48–1570.19170.17650.9206
7Land-use mapAirport0.00840.00220.2644
Barren land0.00090.00040.4666
Dryland agriculture0.05940.00760.1275
Estate crop plantation0.00350.00030.0739
Fish pond0.04750.00120.0250
Rice field0.39790.06500.1634
Secondary mangrove forest0.00030.00010.3295
Settlement area0.45720.91942.0110
Shrub-mixed dryland farms0.01980.00280.1440
Swamp0.00510.00090.1831
Swamp shrub0.00000.00000.1441
8Lithology mapAlluvium0.50200.55771.1111
Alluvium fans0.19660.32001.6274
Beach ridge deposits0.01530.03412.2385
Bojongmanik form0.00220.00000.0000
Cihoe form0.01610.00020.0149
Coastal deposit0.00130.00000.0000
Lake0.00010.00000.1178
Marine deposits0.00520.00000.0024
Old alluvium0.00740.00020.0239
Parigi form0.00050.00000.0000
Sandstone tuff0.00010.00000.0000
Sandstone unit0.04050.04401.0873
Serpong form0.07030.00140.0198
Subang form0.00640.00000.0030
Swamp deposits0.01750.00050.0300
Tuff banten0.06770.03910.5779
Upper banten tuff0.04880.00260.0543
Young volcanic rocks0.00220.00000.0000
9Elevation map (m)0–30.22270.10880.4883
3–90.20290.21701.0699
9–200.19760.37671.9061
20–370.19680.26371.3397
37–1560.18000.03380.1878
10Slope (degree)00.20110.17270.8587
0–1.360.19970.16570.8294
1.36–3.190.19970.23991.2009
3.19–5.470.19970.24701.2367
>5.470.19970.17480.8753
AspectFlat0.12140.13391.1030
North0.14900.17781.1937
Northeast0.10910.13231.2128
East0.11770.12321.0468
11Southeast0.10060.08880.8829
South0.12120.09240.7626
Southwest0.09370.07110.7592
West0.09370.08800.9386
Northwest0.09370.09250.9869
12Profile curvatureConcave0.33320.32530.9764
Flat0.33360.30170.9045
Convex0.33320.37291.1192
13Plan curvatureConcave0.33320.32530.9764
Flat0.33360.30170.9045
Convex0.33320.37291.1192
Topographic wetness index2.52–6.810.14300.15331.0722
6.81–8.000.19390.22521.1614
148.00–10.140.21400.22791.0647
10.14–11.960.22580.20980.9293
11.96–22.900.22330.18390.8232
Table 5. Calculated parameters for the algorithms used in this study.
Table 5. Calculated parameters for the algorithms used in this study.
AlgorithmParameters
AdaBoostThe number of iterations: 10; seed: 1; weight threshold: 100.
LogitBoostNumber of iterations: 10; Seed: 1; weight threshold: 100; likelihood threshold: -1.7976E308; shrinkage: 1.0; max threshold: 3; thread pool: 1; thread to batch prediction: 1.
Logistic RegressionRidge: 1.0E-8; max iterations: -1; number of decimal places: 4.
Multilayer PerceptronHidden layers: a; learning rate: 0.3; momentum: 0.2; number of decimal places: 2; seed: 0; training time: 500; validation set size: 0; validation threshold: 20.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hakim, W.L.; Achmad, A.R.; Lee, C.-W. Land Subsidence Susceptibility Mapping in Jakarta Using Functional and Meta-Ensemble Machine Learning Algorithm Based on Time-Series InSAR Data. Remote Sens. 2020, 12, 3627. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12213627

AMA Style

Hakim WL, Achmad AR, Lee C-W. Land Subsidence Susceptibility Mapping in Jakarta Using Functional and Meta-Ensemble Machine Learning Algorithm Based on Time-Series InSAR Data. Remote Sensing. 2020; 12(21):3627. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12213627

Chicago/Turabian Style

Hakim, Wahyu Luqmanul, Arief Rizqiyanto Achmad, and Chang-Wook Lee. 2020. "Land Subsidence Susceptibility Mapping in Jakarta Using Functional and Meta-Ensemble Machine Learning Algorithm Based on Time-Series InSAR Data" Remote Sensing 12, no. 21: 3627. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12213627

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop