Determining Optimal Location for Mangrove Planting Using Remote Sensing and Climate Model Projection in Southeast Asia

: The decreasing area of mangroves is an ongoing problem since, between 1980 and 2005, one-third of the world’s mangroves were lost. Rehabilitation and restoration strategies are required to address this situation. However, mangroves do not always respond well to these strategies and have high mortality due to several growth limiting parameters. This study developed a land suitability map for new mangrove plantations in di ﬀ erent Southeast Asian countries for both current and future climates at a 250-m resolution. Hydrodynamic, geomorphological, climatic, and socio-economic parameters and three representative concentration pathway (RCP) scenarios (RCP 2.6, 4.5, and 8.5) for 2050 and 2070 with two global climate model datasets (the Centre National de Recherches M é t é orologiques Climate model version 5 [CNRM-CM5.1] and the Model for Interdisciplinary Research on Climate [MIROC5]) were used to predict suitable areas for mangrove planting. An analytical hierarchy process (AHP) was used to determine the level of importance for each parameter. To test the accuracy of the results, the mangrove land suitability analysis were further compared using di ﬀ erent weights in every parameter. The sensitivity test using the Wilcoxon test was also carried out to test which variables had changed with the ﬁrst weight and the AHP weight. The land suitability products from this study were compared with those from previous studies. The di ﬀ erences in land suitability for each country in Southeast Asia in 2050 and 2070 to analyze the di ﬀ erences in each RCP scenario and their e ﬀ ects on the mangrove land suitability were also assessed. Currently, there is 398,000 ha of potentially suitable land for mangrove planting in Southeast Asia, and this study shows that it will increase between now and 2070. Indonesia account for 67.34% of the total land area in the “very suitable” and “suitable” class categories. The RCP 8.5 scenario in 2070, with both the MIROC5 and CNRM-CM5.1 models, resulted in the largest area of a “very suitable” class category for mangrove planting. This study provides information for the migration of mangrove forests to the land, alleviating many drawbacks, especially for ecosystems. each parameter and sub-parameter. After identifying the land suitable for planting mangroves, this study analyzed the socio-economic parameters of the area and their inﬂuences on the land’s suitability. The socio-economic

suitability, the results were compared with the results of the land selection using equal weights (i.e., without AHP). They were also compared with the land suitability products of other studies to analyze the uncertainties in this study. The suitable sites for mangrove growth in 2050 and 2070 were assessed using climate model data (precipitation and average air temperature parameters) from the Centre National de Recherches Météorologiques Climate model version 5 (CNRM-CM5.1) and the Model for Interdisciplinary Research on Climate (MIROC). Three RCP scenarios (RCP 2.6, 4.5, and 8.5) were used for each climate parameter. This study compared the results of the current land suitability with that in 2050 and 2070 for each RCP scenario with the two climate models to determine the differences in each scenario for mangrove planting.
The methodology used in this study was divided into stages. The first stage was to determine the parameters and build their hierarchy. The second stage was to create the base map. The third stage was to classify both the parameters and the sub-parameters into several classes. The fourth stage was to determine the weight of each parameter and sub-parameter using the AHP method. Finally, the fifth stage was to create scenarios to determine the land's suitability for mangroves. Figure 1 shows the data-processing schemes for the land suitability analysis for mangrove planting.

Determination of Parameters and Hierarchy Building
The determination of the appropriate parameters was the most important step for identifying land suitability. The parameters were first selected from a literature review and then the parameters were divided into three categories. The first category included the hydrodynamic, geomorphological, and climatic parameters that made the land suitable for planting mangroves. Moreover, the subparameters used in this study were those that had the highest weight when compared with the other sub-parameters, when weighted using the AHP method, and those with very low weights were not used in this study. The second category included parameters that make a land unsuitable for planting mangroves, owing to intensive social and economic anthropogenic activities. In this study, these parameters were classified as human pressure parameters. Finally, the third category included climate prediction model parameters. These parameters were used to predict the land suitability for

Determination of Parameters and Hierarchy Building
The determination of the appropriate parameters was the most important step for identifying land suitability. The parameters were first selected from a literature review and then the parameters were divided into three categories. The first category included the hydrodynamic, geomorphological, and climatic parameters that made the land suitable for planting mangroves. Moreover, the sub-parameters used in this study were those that had the highest weight when compared with the other sub-parameters, when weighted using the AHP method, and those with very low weights were not used in this study. The second category included parameters that make a land Remote Sens. 2020, 12, 3734 5 of 29 unsuitable for planting mangroves, owing to intensive social and economic anthropogenic activities. In this study, these parameters were classified as human pressure parameters. Finally, the third category included climate prediction model parameters. These parameters were used to predict the land suitability for the future growth of mangroves, while assuming that there would be no changes in the hydrodynamic, geomorphological, and human pressure parameters.

Hydrodynamic
The seawater parameter had one sub-parameter, which was tidal inundation. The tidal inundation data used in this study were obtained from tidal data models. The tidal data model was created from TPXO 9.0 data. TPXO is a tidal model based on bathymetry, which assimilates various data sources and is processed using OTIS software [80,81]. The tidal model was verified using tidal station data across Indonesia, which are available from the Indonesian Geospatial Information Agency. Subsequently, the tidal inundation was calculated by multiplying the amplitude value by two, which results in a tidal inundation map of 2018 that had a spatial resolution of 1 • .

Geomorphology
The geomorphology parameter had two sub-parameters, elevation, and slope. The elevation and slope sub-parameter data used in this study were obtained from the Multi Error Removed Improved Terrain Digital Elevation Model (MERIT DEM). MERIT DEM is the result of the development of DEM data from SRTM3 v2.1 and AW3D-30m v.1, which eliminated several error components, including absolute bias, noise stripe, noise speckle, and height of the tree bias. These data had a spatial resolution of 90 m [82].

Climate
The climatic parameters had two sub-parameters, namely, average air temperature and precipitation. The average air temperature data used in this study were the 2 m air temperature data obtained from the European Center for Medium-Range Weather Forecasts (ECMWF). These air temperature data pertained to 2 m above the land, sea, or water surface. The 2-m air temperature was calculated by interpolating between the lowest level model and the surface of the earth, while taking atmospheric conditions into account. The unit of measurement for this parameter was kelvin (K). The data used were monthly data derived from daily data accumulation with a spatial resolution of 0.125 • .
The precipitation data used in this study were obtained from the Climate Hazards Group InfraRed Precipitation with Stations version 2.0 (CHIRPS v.2.0). The CHIRPS data were a combination of pentadal rainfall data, global geostationary TIR satellite observations from Climate Prediction Center (CPC) and the National Climatic Data Center (NCDC), rainfall atmospheric models from the NOAA Climate Forecast System version 2 (CFSv2), and rainfall data collected from in situ observations at each station [83,84]. This CHIRPS product had a 0.05 • spatial resolution with pentadal, decadal, and monthly temporal resolutions [85].

Human Pressure
Human pressure was the parameter that made an area unsuitable for planting mangroves owing to over-population and development of urban areas. The human pressure parameter had four sub-parameters, land cover, population, gross domestic product (GDP), and night light. The land cover data used in this study were climate change land cover initiative (CCI-LC) data from the European Space Agency (ESA). These data describe the surface of the earth in 37 classes of the original land cover based on the United Nations Land Cover Classification System (UN-LCCS) with a spatial resolution of 300 m [86]. The population data used in this study were gridded populations sourced from the world version 4 (GPWv4) data. These data provide information on the population density (number of people per km 2 ) with a 30" × 30" grid size [87].
The GDP data used in this study were obtained from Kummu et al. [88]. They included USD GDP data, obtained from the GDP per capita multiplied by the History Database of the Global Environment (HYDE) population data 3.2, which had a grid size of 5 × 5 [88]. The night light data used in this study were the black marble nighttime light (VNP46) data of the National Aeronautics and Space Administration (NASA). These data had a spatial resolution of 500 m. NASA developed black marble as a series of daily products that were calibrated, correlated, and validated, so that the night light data could be used effectively for scientific observations [89].

Climate Prediction Model
The climate prediction data model had two sub-parameters, namely, average air temperature and precipitation. These sub-parameters were obtained from the CNRM-CM5.1 and the MIROC5 models. CNRM-CM5.1 and MIROC5 are the earth system models designed to facilitate simulations that contribute to phase five of the Coupled Model Intercomparison Project (CMIP5). CNRM-CM5.1 consists of several models, namely, the atmospheric model from ARPEGE-Climat version 5.2 and the marine model from The Nucleus for European Modelling of the Ocean version 3.2 (NEMO v3.2) and includes the Interaction Sol-Biosphère-Atmosphère (ISBA) land surface and the sea ice model from Global Experimental Leads and sea ice for Atmosphere and Ocean version 5 (GELATO v5). These models were designed independently and combined through OASIS v.3 software [90]. The CNRM-CM5.1 data had a spatial resolution of 1.4 • × 1.4 • and it is one of the best global climate models for temperature and rainfall parameters data in Southeast Asia [91]. MIROC5 is a climate prediction model created by the Japanese research community. It has 40 standard vertical resolution levels of up to 3 hPa with a spatial resolution of 1.4 • × 1.4 • . MIROC5 has better climatological features (precipitation, average zonal atmospheric fields, equatorial subsurface fields, and El Niño-Southern Oscillation simulations) than the previous version [92]. Both CNRM-CM5.1 and MIROC5 included three scenarios, namely, RCP 2.6, 4.5, and 8.5. In RCP 2.6, the greenhouse gas concentration was very low, while RCP 4.5 meant that the total radiative forcing was stabilized shortly. In RCP 8.5, the greenhouse gas emissions increase over time, resulting in higher concentrations. In this study, the results of the mangrove suitability site analyses were compared in the context of these RCPs to determine the scenarios that included the most suitable sites for planting mangroves.

Creation of the Basemap
Before initiating the data processing, a base map was constructed from coastline data and mangrove distributions and buffered 10 km inland. The mangrove distribution data used in this study were obtained from Giri et al. [13]. The primary data for the mangrove distribution resulted from a combination of the global land survey (GLS) data from 1997 to 2000, with a 30-m resolution, and the Landsat imagery available from the United States Geological Survey (USGS). The secondary data used were global mangrove data [14] and data from the national and local mangrove databases. The mangrove distribution data were buffered 10 km inland. This base map was then used for land selection and to analyze land suitability.

Classification of the Parameter
To select suitable land for planting mangroves, each parameter needed to be classified and scored. The classification and scoring in this study involved two stages. The first stage was the division of the classes of each land suitability parameter. Each parameter could be classified into four classes (very suitable, suitable, moderate, and unsuitable). The classes were devised based on the results of a literature review (Table 1). The second stage involved assigning a score to each class. Very suitable and suitable classes were assigned scores of four and three, respectively, whereas moderate and unsuitable classes were assigned scores of two and one, respectively. In addition, class divisions and scores were applied to each human pressure parameter. The human pressure parameters could be categorized into either of four classes (low, medium, high, and very high) except for the land cover parameter. The land cover parameter included two classes (urban and non-urban). Low, medium, high, and very high classes were assigned scores of one, two, three, and four, respectively.

Determination of the Parameters Weight
The AHP method was used to determine the parameter weights in this study. AHP is a decision-making technique developed by Saaty [79] and was carried out in three stages. First, a pairwise comparison matrix was created by creating a scale with values from one to nine for the hydrodynamic, geomorphological, and climate parameters ( Table 2). The scale of importance between one parameter and another was determined by expert judgment. Second, normalize the results of the pairwise comparison matrix were normalized and a column vector created that had n units with n components, so that the weight value and the total weight vector for each parameter could be obtained. Third, the consistency ratios were estimated. To calculate the consistency ratio, the lambda (λ) and consistency index (CI) parameters were required. The λ value was the average value of the consistency vector values of all parameters. The consistency vector value was obtained by dividing the number vector by the weight of each parameter with the following formula.
Vector of consistency VC n th parameter = Number o f vectors n th parameter Weight n th parameter Average consistency (λ) = VC 1 st parameter + VC 2 nd parameter + VC n th parameter The λ value was used to calculate the CI value, as described in Equation (3). The CI calculation was based on the observation that λ is always greater than or equal to the number of the criteria considered (n). If the paired comparison matrix is a consistent matrix, then it has a positive value, a reciprocal matrix, and λ = n. Therefore, if λ-n, then it can be considered as a measure of the inconsistency level. For more details, the CI formula can be written as follows.
The CI value in Equation (3) can be used to calculate the consistency ratio (CR), as described in Equation (4) below.
where RI (random index) is the consistency index of a paired comparison matrix that is randomly generated. If a CR value of less than 0.10 was obtained, it indicated a level of rational/reasonable consistency in pairwise comparisons. However, if the CR value obtained was greater than 0.10, the ratio values indicated inconsistent ratings. If the assessment was inconsistent, a correction needed to be made to the scoring when the comparison matrix was paired. Table 2 shows the pairwise comparison matrix of the main parameters and the value of the CR from this study. Moreover, pairwise comparison matrixes for the sub-parameters are shown in Tables A1-A7. The results of the weighting method using AHP for each parameter and sub-parameter are shown in Table 3.

Scenario Generation of Land Suitability
In this study, five scenarios were created to analyze the suitability of the land for mangrove planting in Southeast Asia. In the first scenario, a land suitability map for mangrove planting was created with the use of the AHP method. In the second scenario, the AHP method was utilized and the influences of the human pressure parameters were also considered. In the third scenario, all parameters were assumed to have the same level of importance, and the AHP method was not used. The influence of the human pressure parameters was not considered, as described in Equation (5). In the fourth scenario, all parameters were assumed to have the same level of importance, and the influence of the human pressure parameters was considered. Finally, in the fifth scenario, a land suitability map for mangrove planting in 2050 and 2070, using the data from the two models (CNRM-CM5.1 and MIROC5) was constructed with the use of the AHP method and influence of the human pressure parameters, as represented in Equations (6) and (7). The scenarios in which the human pressure parameters were considered aimed to assess their influence and how essential it was to consider them when selecting suitable areas for mangroves. The scenarios that did not utilize AHP were used to determine the accuracy of the weight results obtained from the AHP method based on the subjectivity of the experts' judgment. To assess the differences between the AHP effects, the Wilcoxon test was carried out as described previously by Chakraborty et al. [57]. The Wilcoxon test can also be used to compare the mean values of a variable from two paired sample data [100,101], whereas the Wilcoxon signed rank test is used only for interval or ratio type data that does not follow a normal distribution. In addition, Equations (5)-(7), respectively, represent the calculations for the replanting suitability of the mangroves' sites.
Replanting suitability score = Land suitability score − Human pressure score (5) where: m = number of sub-parameters, n = number of sub-parameters, wpi = weight of parameter I, wspij = weight of sub-parameter j in parameter I, wrij = ranking weight of the pixel's sub-parameter j class in parameter i.

Land Suitability for Mangrove Planting as Determined with the AHP Method
This study produced a land suitability map for mangrove planting in Southeast Asia, using the AHP technique, with a 250-m spatial resolution. The potential land suitability area were divided into three classes (very suitable, suitable, and moderate). The detailed results are shown for a selected region of Southeast Asia in Figure 2.
Remote Sens. 2020, 12, x FOR PEER REVIEW 10 of 28 2,725,000 ha (8.58%), 2,637,000 ha (8.31%), 1,679,000 ha (5.29%), and 1,530,000 ha (4.82%), respectively. The countries with the least land potential for mangrove planting were Cambodia (1.06%), Brunei Darussalam (0.29%), and Singapore (0.15%). The inclusion of the human pressure parameters resulted in a significant change in the potential land area that could be suitable for mangrove planting in Southeast Asia. Based on the land suitability map for mangrove planting that was generated using the AHP technique combined with the human pressure parameters (Figure 3), there were approximately 398,000 ha of land suitable for mangrove planting in Southeast Asia. The land that was classified as suitable and moderate for mangrove planting was 4,771,000 ha and 20,123,000 ha, respectively.
Additionally, the order of the countries with the highest land potential for mangrove planting changed. The country with the largest potential land area was still Indonesia, which accounted for 67.34% of the total land area, including land in the very suitable and suitable categories, which was approximately 2,926,000 ha. Indonesia was followed by Vietnam, the Philippines, Thailand, Malaysia, and Cambodia with areas of 457,000 ha (10.52%), 287,000 ha (6.62%), 234,000 ha (5.38%), 182,000 ha (4.18%), and 143,000 ha (3.29%), respectively. The countries with the smallest area of land for mangroves planting were Myanmar (2.57%), Brunei Darussalam (0.09%), and Singapore (0.02%). According to the land suitability map for mangrove planting that was generated using the AHP technique ( Figure 2), approximately 3,960,000 ha of land have the potential to be very suitable for mangrove planting in Southeast Asia. The land that was classified as suitable and moderate for mangrove planting extended over 27,791,000 ha and 16,357,000 ha, respectively. The country that had the greatest potential was Indonesia, as it accounted for 57.38% of the total land area, including the very suitable and suitable categories (approximately 18,220,000 ha). Indonesia was followed by the Philippines, Malaysia, Vietnam, Thailand, and Myanmar with areas of 4,483,000 ha (14.12%), 2,725,000 ha (8.58%), 2,637,000 ha (8.31%), 1,679,000 ha (5.29%), and 1,530,000 ha (4.82%), respectively. The countries with the least land potential for mangrove planting were Cambodia (1.06%), Brunei Darussalam (0.29%), and Singapore (0.15%).
The inclusion of the human pressure parameters resulted in a significant change in the potential land area that could be suitable for mangrove planting in Southeast Asia. Based on the land suitability map for mangrove planting that was generated using the AHP technique combined with the human pressure parameters (Figure 3), there were approximately 398,000 ha of land suitable for mangrove planting in Southeast Asia. The land that was classified as suitable and moderate for mangrove planting was 4,771,000 ha and 20,123,000 ha, respectively.
Additionally, the order of the countries with the highest land potential for mangrove planting changed. The country with the largest potential land area was still Indonesia, which accounted for 67.34% of the total land area, including land in the very suitable and suitable categories, which was approximately 2,926,000 ha. Indonesia was followed by Vietnam, the Philippines, Thailand, Malaysia, and Cambodia with areas of 457,000 ha (10.52%), 287,000 ha (6.62%), 234,000 ha (5.38%), 182,000 ha (4.18%), and 143,000 ha (3.29%), respectively. The countries with the smallest area of land for mangroves planting were Myanmar (2.57%), Brunei Darussalam (0.09%), and Singapore (0.02%).

Current Potential Land Suitability for Mangrove Planting without the AHP Method
The differences in the potentially suitable land for mangrove planting with and without the use of the AHP technique ( Figures 4 and 5), which was used to determine the weight of each parameter, were compared. This comparative analysis aimed to determine the amount of influence the weighting produced by the AHP technique exerted on the potential suitability of the land for mangrove planting.
Based on the map of the land suitability potential without the use of the AHP technique ( Figure 4), approximately 1,149,781 ha of land had the potential to be very suitable for planting mangroves in Southeast Asia. The land that was very suitable and suitable extended over 9,648,900 ha and 23,320,444 ha, respectively. The country with the greatest land potential for mangrove planting was Indonesia, as it accounted for 44.819% of the total land area, including very suitable and suitable categories of approximately 4,840,000 ha. Indonesia was followed by Myanmar, Vietnam, the Philippines, Thailand, and Malaysia with areas of 2,045,000 ha (18.940%), 1,331,000 ha (12.327%), 809,000 ha (7.494%), 781,000 ha (7.231%), and 753,000 ha (6.972%), respectively. Meanwhile, the countries with the least land potential for mangrove planting were Cambodia (1.956%), Brunei Darussalam (0.248%), and Singapore (0.014%).
In addition, this study mapped the potential land suitability for mangrove planting without using AHP, but considering the influence of the human pressure parameters ( Figure 5). The resultant map of the land suitability potential showed that there appeared to be approximately 131,756 ha of potentially very suitable land in Southeast Asia. The very suitable and suitable land areas for mangrove planting were 939,231 ha and 6,988,350 ha, respectively. The country with the greatest land potential for mangrove planting was Indonesia, as it accounted for 59.326% of the total land area, including approximately 635,375 ha of land in the very suitable and suitable categories. The country with the second greatest land area suitable for mangrove planting was the Philippines with 118,806 ha (11.093%). The Philippines was followed by Vietnam, Thailand, Cambodia, and Malaysia with 104,063 ha (9.716%), 92,600 ha (8.646%), 55,025 ha (5.138%), and 36,488 ha (3.407%), respectively. The countries with the least land potential for mangrove planting were Myanmar (2.652%), Singapore (0.015%), and Brunei Darussalam (0.007%).

Current Potential Land Suitability for Mangrove Planting without the AHP Method
The differences in the potentially suitable land for mangrove planting with and without the use of the AHP technique (Figures 4 and 5), which was used to determine the weight of each parameter, were compared. This comparative analysis aimed to determine the amount of influence the weighting produced by the AHP technique exerted on the potential suitability of the land for mangrove planting.     It is evident from the results that the analysis without the use of the AHP technique resulted in smaller areas of suitable and very suitable land (Table 4). This was because each parameter had a different weighed value, and these weight differences caused the higher-weighted parameters to be more important than the lower-weighted parameters. The results of the AHP technique indicated that the geomorphological parameters had the highest weight, accounting for approximately 38% of the total weight, which were followed by the hydrodynamic parameters and accounted for 32% of the total weight. The parameter with the smallest weight was the climate, which accounted for 30% of the total weight. These results showed that the geomorphological parameters were more important Figure 5. Map of the current potential land suitability for mangrove planting without the use of the analytic hierarchy process (AHP) method, but with the influence of human pressure.
Based on the map of the land suitability potential without the use of the AHP technique (Figure 4), approximately 1,149,781 ha of land had the potential to be very suitable for planting mangroves in Southeast Asia. The land that was very suitable and suitable extended over 9,648,900 ha and 23,320,444 ha, respectively. The country with the greatest land potential for mangrove planting was Indonesia, as it accounted for 44.819% of the total land area, including very suitable and suitable categories of approximately 4,840,000 ha. Indonesia was followed by Myanmar, Vietnam, the Philippines, Thailand, and Malaysia with areas of 2,045,000 ha (18.940%), 1,331,000 ha (12.327%), 809,000 ha (7.494%), 781,000 ha (7.231%), and 753,000 ha (6.972%), respectively. Meanwhile, the countries with the least land potential for mangrove planting were Cambodia (1.956%), Brunei Darussalam (0.248%), and Singapore (0.014%).
In addition, this study mapped the potential land suitability for mangrove planting without using AHP, but considering the influence of the human pressure parameters ( Figure 5). The resultant map of the land suitability potential showed that there appeared to be approximately 131,756 ha of potentially very suitable land in Southeast Asia. The very suitable and suitable land areas for mangrove planting were 939,231 ha and 6,988,350 ha, respectively. The country with the greatest land potential for mangrove planting was Indonesia, as it accounted for 59.326% of the total land area, including approximately 635,375 ha of land in the very suitable and suitable categories. The country with the second greatest land area suitable for mangrove planting was the Philippines with 118,806 ha (11.093%). The Philippines was followed by Vietnam, Thailand, Cambodia, and Malaysia with 104,063 ha (9.716%), 92,600 ha (8.646%), 55,025 ha (5.138%), and 36,488 ha (3.407%), respectively. The countries with the least land potential for mangrove planting were Myanmar (2.652%), Singapore (0.015%), and Brunei Darussalam (0.007%).
It is evident from the results that the analysis without the use of the AHP technique resulted in smaller areas of suitable and very suitable land (Table 4). This was because each parameter had a different weighed value, and these weight differences caused the higher-weighted parameters to be more important than the lower-weighted parameters. The results of the AHP technique indicated that the geomorphological parameters had the highest weight, accounting for approximately 38% of the total weight, which were followed by the hydrodynamic parameters and accounted for 32% of the total weight. The parameter with the smallest weight was the climate, which accounted for 30% of the total weight. These results showed that the geomorphological parameters were more important than the hydrodynamic and climate parameters. However, the weights of each of these parameters added up to the weights of each sub-parameter, which means that an area that has a large geomorphological parameter value may be no more suitable for mangrove planting than areas with high hydrodynamic values. This proves that the weight of each sub-parameter can be influential on the selection of potentially suitable lands for mangrove planting. Table 4. Areas of land classified as very suitable, suitable, and moderate in terms of planting mangroves in Southeast Asia both with and without the use of the analytic hierarchy process (AHP). The results of the AHP technique analysis showed that the tidal inundation sub-parameter had the greatest weight when compared with the other sub-parameters. This was approximately four times the value of the smallest weight. The greater weight of this sub-parameter suggests that an area with tidal inundation classified in the very suitable, suitable, or moderate categories for mangrove planting will be preferred over areas in which other sub-parameters may have greater weights, even if they are classified as very suitable. In addition to the tidal inundation sub-parameter, the air temperature sub-parameter played an important role in the selection of potentially suitable land for mangrove planting. The air temperature sub-parameter ranked second in terms of importance after the tidal inundation sub-parameter and its weight was approximately three times that of the smallest weight. This suggests that the air temperature sub-parameter should be prioritized over the other three sub-parameters when selecting land that is potentially suitable for mangrove planting. The elevation and precipitation sub-parameters had the same weight and were considered equally important. However, when compared with the tidal inundation and air temperature, they were less important. The sub-parameter with the smallest weight, according to the results of the AHP technique, was slope.

Weighting
The analysis of the potential land suitability for mangrove planting without the use of the AHP technique considered all parameters and sub-parameters as equally important in the selection of potentially suitable land for mangrove planting. The equalization of the parameter and sub-parameter weights gave a less accurate selection of the potentially suitable land than when the AHP method was used. Inaccurate land selection may result in the use of unsuitable land for mangrove planting or growth.

Land Suitability for Mangrove Planting Using the AHP Method and Climate Models
Two climate models were used to predict the potential land suitability for mangrove planting in Southeast Asia in 2050 and 2070. These were the CNRM-CM5.1 and the MIROC5. In the land suitability map for mangrove planting in 2050 constructed using the CNRM-CM5.1 model, the 2.6, 4.5, and 8.5 RCP scenarios were examined ( Figure 6). There was a slight difference in the areas of the suitable land for mangrove planting among these three scenarios. RCP 8.5 resulted in the largest potential land area (324,125 ha), which was followed by RCP 4.5 (307,313 ha). RCP 2.6 resulted in the smallest suitable land area (264,906 ha). The difference in the land area of the very suitable category between RCP 8.5 and 2.6 was 59,219 ha. In the land suitability map for mangrove planting in 2050 constructed using the CNRM-CM5.1 model, the 2.6, 4.5, and 8.5 RCP scenarios were examined ( Figure 6). There was a slight difference in the areas of the suitable land for mangrove planting among these three scenarios. RCP 8.5 resulted in the largest potential land area (324,125 ha), which was followed by RCP 4.5 (307,313 ha). RCP 2.6 resulted in the smallest suitable land area (264,906 ha). The difference in the land area of the very suitable category between RCP 8.5 and 2.6 was 59,219 ha. The results of the three RCPs were also different in terms of the area of the potential land area that was very suitable for mangrove planting in each country (Figure 7). Once again, Indonesia had the highest land suitability in all three RCPs. Additionally, Indonesia showed an increasing area of very suitable land from RCP 2.6 to 8.5 (52,181 ha). This increase was also observed for several other Southeast Asian countries, namely, the Philippines, Myanmar, and Malaysia, with increases of 6650 ha, 1331 ha, and 244 ha, respectively. However, Thailand, Vietnam, and Brunei Darussalam The results of the three RCPs were also different in terms of the area of the potential land area that was very suitable for mangrove planting in each country (Figure 7). Once again, Indonesia had the highest land suitability in all three RCPs. Additionally, Indonesia showed an increasing area of very suitable land from RCP 2.6 to 8.5 (52,181 ha). This increase was also observed for several other Southeast Asian countries, namely, the Philippines, Myanmar, and Malaysia, with increases of 6650 ha, 1331 ha, and 244 ha, respectively. However, Thailand, Vietnam, and Brunei Darussalam maintained the same area of very suitable potential land with all three RCPs. Singapore had the same potential land area with RCP 2.6 and 4.5, and it increased with RCP 8.5 by 13 ha. The suitable potential land area in Cambodia increased by 156 ha from RCP 2.6 to RCP 4.5. However, there was a drastic reduction of 1200 ha with RCP 8.5.  (Figures 8 and 9), the land area that was potentially very suitable for mangrove planting remained relatively similar for each RCP in most countries including Thailand, the Philippines, Vietnam, Malaysia, Brunei Darussalam, and Singapore. Indonesia experienced a decrease of 1219 ha in the potential land area from RCP 2.6 to RCP 4.5. However, the land area that was very suitable for mangrove planting in Indonesia increased by 3356 ha with RCP 8.5. Like Indonesia, the potential land area that was very suitable in Myanmar decreased by 281 ha from RCP 2.6 to RCP 4.5, and increased by 2519 ha with RCP 8.5. Conversely, in Cambodia, the very suitable potential land area decreased by 3813 ha from RCP 2.6 to RCP 8.5. In the land suitability map constructed using the MIROC5 model 2050, the 2.6, 4.5, and 8.5 RCP scenarios were also assessed ( Figure 8). There was a difference in the potential land area that was very suitable for mangrove planting among the three scenarios. The largest potential land area for mangrove planting in RCP 8.5 was 324,125 ha, whereas RCP 4.5 had a very suitable potential land area of 307,313 ha. RCP 2.6 had the least very suitable potential land area (264,906 ha). Therefore, the difference in the very suitable land area from RCP 8.5 to 2.6 was 82,838 ha.
The three RCPs also differed in terms of the potential land area that was very suitable in each country ( Figure 9). Indonesia had the largest potential land area suitable in all three RCPs, and it increased by 75,800 ha from RCP 2.6 to 8.5. An increase also occurred in several other Southeast Asian countries, namely, the Philippines, Myanmar, Malaysia, and Singapore. These experienced increases of 6413 ha, 1138 ha, 231 ha, and 13 ha, respectively. The area of the very suitable potential land in Vietnam and Brunei Darussalam remained the same in the three RCPs. The very suitable potential land area in Thailand decreased by 69 ha from RCP 2.6 to 4.5 and then increased by 3 ha with RCP 8.5. The suitable potential land area in Cambodia decreased by 719 ha from RCP 2.6 to RCP 8.5. In contrast to 2050, for the three RCPs in 2070 (Figures 8 and 9), the land area that was potentially very suitable for mangrove planting remained relatively similar for each RCP in most countries including Thailand, the Philippines, Vietnam, Malaysia, Brunei Darussalam, and Singapore. Indonesia experienced a decrease of 1219 ha in the potential land area from RCP 2.6 to RCP 4.5. However, the land area that was very suitable for mangrove planting in Indonesia increased by 3356 ha with RCP 8.5. Like Indonesia, the potential land area that was very suitable in Myanmar decreased by 281 ha from RCP 2.6 to RCP 4.5, and increased by 2519 ha with RCP 8.5. Conversely, in Cambodia, the very suitable potential land area decreased by 3813 ha from RCP 2.6 to RCP 8.5.

Potential Land Suitability for Mangrove Planting in 2050 and 2070 Using the MIROC5 Model and the AHP Method
In the land suitability map constructed using the MIROC5 model 2050, the 2.6, 4.5, and 8.5 RCP scenarios were also assessed ( Figure 8). There was a difference in the potential land area that was very suitable for mangrove planting among the three scenarios. The largest potential land area for mangrove planting in RCP 8.5 was 324,125 ha, whereas RCP 4.5 had a very suitable potential land area of 307,313 ha. RCP 2.6 had the least very suitable potential land area (264,906 ha). Therefore, the difference in the very suitable land area from RCP 8.5 to 2.6 was 82,838 ha.     The three RCPs also differed in terms of the potential land area that was very suitable in each country ( Figure 9). Indonesia had the largest potential land area suitable in all three RCPs, and it increased by 75,800 ha from RCP 2.6 to 8.5. An increase also occurred in several other Southeast Asian countries, namely, the Philippines, Myanmar, Malaysia, and Singapore. These experienced increases of 6413 ha, 1138 ha, 231 ha, and 13 ha, respectively. The area of the very suitable potential land in Vietnam and Brunei Darussalam remained the same in the three RCPs. The very suitable potential land area in Thailand decreased by 69 ha from RCP 2.6 to 4.5 and then increased by 3 ha with RCP 8.5. The suitable potential land area in Cambodia decreased by 719 ha from RCP 2.6 to RCP 8.5.
There were differences in the potentially suitable land area for mangrove planting (Figure 8). The largest area of potential very suitable land was found with RCP 2.6, which was 330,919 ha, whereas RCP 8.5 had the largest area of potential land area in the suitable category (329,381 ha). RCP 4.5 had the smallest land area in the very suitable category (321,506 ha) when compared with the other two RCPs. Therefore, there was a difference of 9413 ha in the very suitable land area from RCP 2.6 to RCP 4.5. Apart from the differences in the area of the very suitable potential land, the results of the three RCPs were relatively stable (Figure 8). The countries where the land area remained the same were Thailand, the Philippines, Vietnam, Malaysia, Brunei Darussalam, and Singapore. In contrast to the other countries, Indonesia experienced a drastic decrease of 5888 ha in a very suitable land area from RCP 2.6 to RCP 4.5. However, it experienced an increase of 5813 ha in the very suitable land area in RCP 8.5. Therefore, the decrease from RCP 2.6 to RCP 8.5 was 75 ha.

Comparison of Our Results with Those of Other Studies
An additional feature of this study was the comparison of the results of the potential land suitability for mangroves that were obtained (with or without the use of the AHP method) with the results of Worthington and Spalding (2019, [102]) ( Figure 10). The results of Worthington and Spalding (2019, [102]) showed that the total area with restoration potential was 303,708 ha. According to our results, the total area with the greatest potential (very suitable) to plant mangroves (with the use of the AHP technique) was 398,000 ha. Without the use of the AHP technique, this area was reduced to 131,756 ha.
Remote Sens. 2020, 12, x FOR PEER REVIEW 17 of 28 There were differences in the potentially suitable land area for mangrove planting (Figure 8). The largest area of potential very suitable land was found with RCP 2.6, which was 330,919 ha, whereas RCP 8.5 had the largest area of potential land area in the suitable category (329,381 ha). RCP 4.5 had the smallest land area in the very suitable category (321,506 ha) when compared with the other two RCPs. Therefore, there was a difference of 9413 ha in the very suitable land area from RCP 2.6 to RCP 4.5. Apart from the differences in the area of the very suitable potential land, the results of the three RCPs were relatively stable (Figure 8). The countries where the land area remained the same were Thailand, the Philippines, Vietnam, Malaysia, Brunei Darussalam, and Singapore. In contrast to the other countries, Indonesia experienced a drastic decrease of 5888 ha in a very suitable land area from RCP 2.6 to RCP 4.5. However, it experienced an increase of 5813 ha in the very suitable land area in RCP 8.5. Therefore, the decrease from RCP 2.6 to RCP 8.5 was 75 ha.

Comparison of our Results with Those of Other Studies
An additional feature of this study was the comparison of the results of the potential land suitability for mangroves that were obtained (with or without the use of the AHP method) with the results of Worthington and Spalding (2019, [102]) ( Figure 10). The results of Worthington and Spalding (2019, [102]) showed that the total area with restoration potential was 303,708 ha. According to our results, the total area with the greatest potential (very suitable) to plant mangroves (with the use of the AHP technique) was 398,000 ha. Without the use of the AHP technique, this area was reduced to 131,756 ha.
In addition, there were differences in the area of land that was determined to be very suitable for mangrove planting in each country ( Figure 11). The sequence of countries obtained in this investigation that had the potential for mangrove planting is also different from the previously reported results of Worthington and Spalding (2019, [102]). In this investigation, the sequence of countries, in order from most suitable to least suitable, was Indonesia, Thailand, the Philippines, Cambodia, Vietnam, Myanmar, Malaysia, Singapore, and Brunei Darussalam. The corresponding sequence of the countries obtained without the use of the AHP technique was Thailand, Indonesia, the Philippines, Cambodia, Vietnam, Malaysia, Myanmar, Brunei Darussalam, and Singapore. According to the results of Worthington and Spalding (2019, [102]), the countries with the highest potentials for mangrove restoration were Indonesia, Myanmar, Thailand, Vietnam, Malaysia, the Philippines, Cambodia, Brunei Darussalam, and Singapore in order of potential. The differences were caused by the variations in the use of environmental parameters during data processing.  In addition, there were differences in the area of land that was determined to be very suitable for mangrove planting in each country ( Figure 11). The sequence of countries obtained in this investigation that had the potential for mangrove planting is also different from the previously reported results of Worthington and Spalding (2019, [102]). In this investigation, the sequence of countries, in order from most suitable to least suitable, was Indonesia, Thailand, the Philippines, Cambodia, Vietnam, Myanmar, Malaysia, Singapore, and Brunei Darussalam. The corresponding sequence of the countries obtained without the use of the AHP technique was Thailand, Indonesia, the Philippines, Cambodia, Vietnam, Malaysia, Myanmar, Brunei Darussalam, and Singapore. According to the results of Worthington and Spalding (2019, [102]), the countries with the highest potentials for mangrove restoration were Indonesia, Myanmar, Thailand, Vietnam, Malaysia, the Philippines, Cambodia, Brunei Darussalam, and Singapore in order of potential. The differences were caused by the variations in the use of environmental parameters during data processing.
Remote Sens. 2020, 12, x FOR PEER REVIEW 18 of 28 Figure 11. Comparison of the results of the land suitability for mangrove planting in Southeast Asia obtained from this study (with and without the use of the analytic hierarchy process (AHP) method) and the potential restoration results obtained from Worthington and Spalding (2019, [102]).
In addition, in this study, human pressure parameters were included. Therefore, an area that was classified as very suitable was selected not only because of the suitability of the environmental parameters but also because of the low human pressure parameter values. A high human pressure parameter value in a country was caused by high population growth and the GDP. In such a scenario, a larger land area will be required for the construction of socio-economic facilities. The results of this study showed that there were some countries with a smaller land area that were very suitable for mangrove planting when compared with other countries, which may have had larger mangrove areas.
In this study, all parameters were derived from remote sensing data and included environmental (hydrodynamic and geomorphology), climate, and human pressure parameters. However, the parameters used by Worthington and Spalding (2019, [102]) were tidal range, sea level rise, projected future sea level rise, sediment change, average size of mangrove loss, and proximity of the lost area to the remaining mangrove forest.

Comparison RCPs Scenario in Terms of the Potential Land Suitability for Mangrove Planting in the Future
In this study, two global climate models were used (namely, the CNRM-CM5.1 and the MIROC5) to predict the potential land suitability for mangrove planting in 2050 and 2070 with the RCP 2.6, 4.5, and 8.5. The different uses of the models created a number of different models of potential land suitability for mangrove planting. As shown in Figure 12, RCP 8.5 had the most suitable potential land area. This was due to the fact that its temperature increases were greater than those for the other RCPs (from 18.5 °C to 30.5 °C). In RCP 2.6 and RCP 4.5, the temperature increase was slightly lesser than that in RCP 8.5, from 18 °C to 29.9 °C, and 18.2 °C to 30.1 °C, respectively. Mangroves live and grow optimally at temperatures of 28 °C to 30 °C. This means that most Southeast Asian countries will have very suitable land area for the mangrove in RCP 8.5 in 2050. In contrast, mangrove land suitability predictions for 2070 showed that the majority of Southeast Asian countries maintained the same area across the three RCPs. This occurred because, in 2070, the lowest temperature was predicted to occur in RCP 4.5 (18.5-30.4 °C). RCP 8.5 was predicted to have the highest temperature in 2070, which ranged from 19.2 °C to 31.3 °C. On the other hand, RCP 2.6 had an intermediate temperature between RCP 4.5 and RCP 8.5 (18.5-30.5 °C). The temperature increase from RCP 2.6 to RCP 8.5 both in 2050 and 2070 means that the majority of Southeast Asian countries would have the same pattern in terms of very suitable land area for mangroves planting (Figure 13). Figure 11. Comparison of the results of the land suitability for mangrove planting in Southeast Asia obtained from this study (with and without the use of the analytic hierarchy process (AHP) method) and the potential restoration results obtained from Worthington and Spalding (2019, [102]).
In addition, in this study, human pressure parameters were included. Therefore, an area that was classified as very suitable was selected not only because of the suitability of the environmental parameters but also because of the low human pressure parameter values. A high human pressure parameter value in a country was caused by high population growth and the GDP. In such a scenario, a larger land area will be required for the construction of socio-economic facilities. The results of this study showed that there were some countries with a smaller land area that were very suitable for mangrove planting when compared with other countries, which may have had larger mangrove areas.
In this study, all parameters were derived from remote sensing data and included environmental (hydrodynamic and geomorphology), climate, and human pressure parameters. However, the parameters used by Worthington and Spalding (2019, [102]) were tidal range, sea level rise, projected future sea level rise, sediment change, average size of mangrove loss, and proximity of the lost area to the remaining mangrove forest.

Comparison RCPs Scenario in Terms of the Potential Land Suitability for Mangrove Planting in the Future
In this study, two global climate models were used (namely, the CNRM-CM5.1 and the MIROC5) to predict the potential land suitability for mangrove planting in 2050 and 2070 with the RCP 2.6, 4.5, and 8.5. The different uses of the models created a number of different models of potential land suitability for mangrove planting. As shown in Figure 12, RCP 8.5 had the most suitable potential land area. This was due to the fact that its temperature increases were greater than those for the other RCPs (from 18.5 • C to 30.5 • C). In RCP 2.6 and RCP 4.5, the temperature increase was slightly lesser than that in RCP 8.5, from 18 • C to 29.9 • C, and 18.2 • C to 30.1 • C, respectively. Mangroves live and grow optimally at temperatures of 28 • C to 30 • C. This means that most Southeast Asian countries will have very suitable land area for the mangrove in RCP 8.5 in 2050. In contrast, mangrove land suitability predictions for 2070 showed that the majority of Southeast Asian countries maintained the same area across the three RCPs. This occurred because, in 2070, the lowest temperature was predicted to occur in RCP 4.5 (18.5-30.4 • C). RCP 8.5 was predicted to have the highest temperature in 2070, which ranged from 19.2 • C to 31.3 • C. On the other hand, RCP 2.6 had an intermediate temperature between RCP 4.5 and RCP 8.5 (18.5-30.5 • C). The temperature increase from RCP 2.6 to RCP 8.5 both in 2050 and 2070 means that the majority of Southeast Asian countries would have the same pattern in terms of very suitable land area for mangroves planting ( Figure 13). Furthermore, the increased air temperature in 2050 and 2070 were not predicted across the temperature limits required for mangrove survival.
year, respectively). The increase in precipitation will also increase the number of areas suitable for planting mangroves in 2070, especially for RCP 8.5. This is consistent with the findings of Eslami-Andargoli et al. (2009, [39]), which increased precipitation, would result in the migration of mangrove forests in land and, consequently, increase the number of areas suitable for planting mangroves. However, the increasing number of inland mangrove areas with RCP 8.5 in 2070 can result in changes to the zonation of mangrove species [17,44,45]. An increase in precipitation could also reduce salinity [44,103,104]. Furthermore, each mangrove species has different hydrodynamic, geomorphological, and climatic criteria. Thus, if there is an increase in precipitation and temperature, some mangrove species would not survive and would eventually become extinct. At the same time, several surviving mangrove species would increase peat production due to decreased salinity and increased freshwater retention [105,106].  planting mangroves in 2070, especially for RCP 8.5. This is consistent with the findings of Eslami-Andargoli et al. (2009, [39]), which increased precipitation, would result in the migration of mangrove forests in land and, consequently, increase the number of areas suitable for planting mangroves. However, the increasing number of inland mangrove areas with RCP 8.5 in 2070 can result in changes to the zonation of mangrove species [17,44,45]. An increase in precipitation could also reduce salinity [44,103,104]. Furthermore, each mangrove species has different hydrodynamic, geomorphological, and climatic criteria. Thus, if there is an increase in precipitation and temperature, some mangrove species would not survive and would eventually become extinct. At the same time, several surviving mangrove species would increase peat production due to decreased salinity and increased freshwater retention [105,106].   According to Duke et al. [45], mangroves will thrive with precipitation levels of 1400 to 3750 mm per year. According to the climate model data of MIROC5 and CNRM-CM5.1 in 2050 and 2070, the amount of precipitation per year in Southeast Asia will range from 0-1257 mm. For RCP 8.5 in 2070, the MIROC5 and CNRM-CM5.1 models had the highest precipitation values (12.574 and 988 mm per year, respectively). The increase in precipitation will also increase the number of areas suitable for planting mangroves in 2070, especially for RCP 8.5. This is consistent with the findings of Eslami-Andargoli et al. (2009, [39]), which increased precipitation, would result in the migration of mangrove forests in land and, consequently, increase the number of areas suitable for planting mangroves. However, the increasing number of inland mangrove areas with RCP 8.5 in 2070 can result in changes to the zonation of mangrove species [17,44,45]. An increase in precipitation could also reduce salinity [44,103,104]. Furthermore, each mangrove species has different hydrodynamic, geomorphological, and climatic criteria. Thus, if there is an increase in precipitation and temperature, some mangrove species would not survive and would eventually become extinct. At the same time, several surviving mangrove species would increase peat production due to decreased salinity and increased freshwater retention [105,106].

Uncertainties in Selecting Land Suitable for Mangrove Replanting
In this study, hydrodynamic, geomorphological, climatic, and socio-economic remote sensing data were used to produce maps of mangrove land suitability. The use of remote sensing data creates uncertainties related to the selection of the most suitable sites for mangrove planting and restoration. This is due to several reasons. First, the data used had different spatial resolutions. This can affect the interpretation of the suitability of the land for mangrove planting. Second, there were many parameters that affected the suitability of the land in terms of mangrove planting. However, most of these parameters have no available data. Some data were available, but their spatial resolutions were so coarse that they could not be used in this study. Owing to the limitations in the availability of these data, several other land suitability parameters were selected for this study. The parameters chosen had a higher level of importance in terms of the suitability of the land for mangrove planting than that of the other parameters. The AHP method was used to determine the importance level of each parameter. The determination of the level of importance for each parameter in the context of the AHP method resulted from expert judgment. However, some experts had conflicting opinions even though the results of the experts' consistency assessments were consistent among the parameters. For example, expert A said that the tidal inundation parameter was the most influential parameter for selecting land suitable for mangrove planting, whereas expert B said that the most important parameter was slope. Both experts assigned values to each parameter consistently using the same rationale. This situation can also affect the land suitability results. Finally, the maps that predict the mangrove land suitability for 2050 and 2070 were created based on the assumption that hydrodynamic, geomorphological, and socio-economic parameters had fixed values based on the current situation even though this may not be the case. Furthermore, the sea-level rise parameters was not considered in this study. However, the sea-level rise is the most pertinent from global warming because it can change the duration of the swamps, their frequency, and salinity [46,107]. Therefore, the predicted results for 2050 and 2070 may be less accurate. Table A3 explains that the CI value of 0.034 meant that the value of the relationship for each different variable had a fairly different similarity index and was also represented as having been sampled from our domain. Then, the sensitivity test was used to test which variables had changed with the first weight and the AHP weight. After the reduction, the differences in the tidal inundation (0.29), slope (0.19), elevation (0.23), air temperature (0.24), and precipitation (0.21) were obtained. In line with this, the slope parameter had the smallest difference in sensitivity with the total difference in sensitivity 5 i=1 (S i ) = 1.16. It can then be interpreted that the use and employment of AHP here was appropriate, as proven by the Wilcoxon Test, which is a non-parametric version of the Student's t-test with value (W = 25) and p-value = 0.007937 less than the significance level α = 5%. The true location shift is not equal to 0. The weight is the most important aspect in determining model performance, after which it is then necessary to pay attention to the parameters that are significant for model construction. If the weight increases, it takes longer time for a model to reach the maximum solution or convergence. Meanwhile, the Wilcoxon value also proved that the weight obtained was correct.
In addition, the results of this study were compared with several similar studies, namely, by Chakraborty et al. (2019, [56]) and Worthington and Spalding (2019, [102]) ( Figure 14). The land suitability of the mangrove in the study by Worthington and Spalding [102] focused on areas lost since 1996. Factors such as tidal ranges, recent sea-level rises, projections of future sea-level rises, and recent sediment changes were added. The time since loss, mean size of patch lost, and proximity of the area lost to the remaining mangrove forests were analyzed using the Delphi method. Meanwhile, the land suitability of mangrove in another study [56] was made based on the climatic, geomorphic, edaphic, and floral conditions and the human interface parameters that were analyzed using the AHP method in the Northern and Central Andaman regions. Although there are some differences in the parameters and the methods used, the results of the mangrove land suitability maps produced by the three studies are very similar. For example, the results of the three studies in the Mayabunder area indicated that this location was relatively unsuitable for mangrove planting, whereas, in the Shyamkund Rangat area, all three studies indicated that most of the area was suitable for mangrove planting. Another example is that all three studies indicated that the Sahari Rangat area was very suitable for planting mangroves (Figure 14). This proves that some of the same parameters used in these three studies have a significant role in selecting suitable locations for mangroves. In addition, there are similarities in several locations suitable for mangrove planting.
Remote Sens. 2020, 12, x FOR PEER REVIEW 21 of 28 three studies are very similar. For example, the results of the three studies in the Mayabunder area indicated that this location was relatively unsuitable for mangrove planting, whereas, in the Shyamkund Rangat area, all three studies indicated that most of the area was suitable for mangrove planting. Another example is that all three studies indicated that the Sahari Rangat area was very suitable for planting mangroves (Figure 14). This proves that some of the same parameters used in these three studies have a significant role in selecting suitable locations for mangroves. In addition, there are similarities in several locations suitable for mangrove planting.

Possible Future Directions
The results of this study can be applied in future research both in the environmental and the socio-economic sectors. One of the applications of this study could be to identify suitable land for planting mangroves in Southeast Asia. The reference map of the mangrove land suitability can help select areas that will be subjected to restoration before field surveys are conducted. However, to restore these areas, consideration should be paid to other mangrove land suitability parameters that were not included in this analysis. This will ultimately help minimize failures when planting mangroves. Part of this application includes an effort to implement sustainable development goals (SDGs) 13 and 14 regarding climate change and protecting the marine ecosystem [108].
A range of other important parameters should be included in future studies that are related to the suitability of land for mangrove planting. For example, determining mangrove species that are suitable to plant in areas subjected to either restoration or rehabilitation is of importance. Sea-level rises are also a restoration consideration in many areas. In addition, parameters for the future population growth models can also be added as barriers for the suitability of the land for future mangrove planting. Another research direction that could be followed is the examination of how much carbon will be absorbed if a number of areas are subjected to restoration. The findings of these studies will help elucidate how much greenhouse gas emissions will be reduced if mangrove restoration is carried out in a number of areas. This can support the implementation of greenhouse gas emission reduction targets in accordance with the Paris agreement [109]. The prediction of future mangrove suitability areas in 2050 and 2070 could also be improved by considering long-term changes in the anthropogenic and naturogenic impact [110], such as in agriculture areas [111,112], urban areas [113], and aquaculture practices [114]. These products could be integrated as one input data to make future analyses more realistic. Improving pixel resolution should also be considered to

Possible Future Directions
The results of this study can be applied in future research both in the environmental and the socio-economic sectors. One of the applications of this study could be to identify suitable land for planting mangroves in Southeast Asia. The reference map of the mangrove land suitability can help select areas that will be subjected to restoration before field surveys are conducted. However, to restore these areas, consideration should be paid to other mangrove land suitability parameters that were not included in this analysis. This will ultimately help minimize failures when planting mangroves. Part of this application includes an effort to implement sustainable development goals (SDGs) 13 and 14 regarding climate change and protecting the marine ecosystem [108].
A range of other important parameters should be included in future studies that are related to the suitability of land for mangrove planting. For example, determining mangrove species that are suitable to plant in areas subjected to either restoration or rehabilitation is of importance. Sea-level rises are also a restoration consideration in many areas. In addition, parameters for the future population growth models can also be added as barriers for the suitability of the land for future mangrove planting. Another research direction that could be followed is the examination of how much carbon will be absorbed if a number of areas are subjected to restoration. The findings of these studies will help elucidate how much greenhouse gas emissions will be reduced if mangrove restoration is carried out in a number of areas. This can support the implementation of greenhouse gas emission reduction targets in accordance with the Paris agreement [109]. The prediction of future mangrove suitability areas in 2050 and 2070 could also be improved by considering long-term changes in the anthropogenic and naturogenic impact [110], such as in agriculture areas [111,112], urban areas [113], and aquaculture practices [114]. These products could be integrated as one input data to make future analyses more realistic. Improving pixel resolution should also be considered to improve future analyses. Several techniques such as a spectral mixture model [115,116] and climate model interpolation [117,118] are key for producing high resolution images of mangrove replanting sites.

Conclusions
A map of the land suitable for mangrove planting in Southeast Asia was produced with a spatial resolution of 250 m, and was divided into four scenarios. The first scenario included the land suitability map for mangrove planting with the use of the AHP method. The second scenario included the AHP method and the influence of human pressure parameters. In the third scenario, all parameters were assumed to have the same level of importance and did not use the AHP model. Finally, in the fourth scenario, it was assumed that all parameters had the same level of importance, and the influence of the human pressure parameters was considered, but the AHP method was not used. The first and the second scenarios showed that approximately 3,960,000 ha and 398,000 ha of land, respectively, appeared to have the potential to be very suitable for mangrove planting. In scenarios three and four, approximately 1,149,781 ha and 131,756 ha of land, respectively, appeared to have the potential to be suitable for mangrove planting. All four scenarios showed that the country with the largest land area suitable for mangrove planting was Indonesia, which accounted for approximately 50%-60% of the total land included in the very suitable category in Southeast Asia. Moreover, this study presented a potential land suitability map for mangrove planting for 2050 and 2070 using two climate models (CNRM-CM5.1 and MIROC5) for each year. Each climate model had three RCP scenarios (2.6, 4.5, and 8.5). The results from both models showed that RCP 8.5 resulted in the largest land area suitable for mangrove planting in 2050 in the majority of the Southeast Asian countries. In 2070, almost all Southeast Asian countries had the same land area suitable for mangrove planting with all three RCPs. Future research could use feature selection in machine learning to overcome the shortcomings of the AHP method.
Author Contributions: L.N.S. and A.D.S. were responsible for the overall design of the study. L.N.S. and R.V. were responsible for the Geographic Information System (GIS) method. W.W. was responsible for processing the hydrography data. R.E.C. validated the design of the study and ran statistical analysis. K.W., S.T. and R.P. supported the model design of the study. L.N.S. wrote the paper. All authors have read and agreed to the published version of the manuscript.  Table A2. The weight of hydrodynamic sub-parameters.

Sub-Parameter Weight
Sea water temperature 0.15 Energy from sea wave and sea tide 0.31 Tidal inundation 0.54 Table A3. Pairwise comparison matrix for determination of the weight of geomorphology sub-parameters.  Table A5. Pairwise comparison matrix for determination of the weight of climate sub-parameters.