A Comprehensive Assessment of Climate Change and Coastal Inundation through Satellite-Derived Datasets: A Case Study of Sabang Island, Indonesia

: Climate-change-induced hazards are negatively affecting the small islands across Indonesia. Sabang Island is one of the most vulnerable small islands due to the rising sea levels and increasing coastal inundation which threaten the low-lying coastal areas with and without coastal defences. However, there is still a lack of studies concerning the long-term trends in climatic variables and, consequently, sea level changes in the region. Accordingly, the current study attempts to comprehensively assess sea level changes and coastal inundation through satellite-derived datasets and model-based products around Sabang Island, Indonesia. The ﬁndings of the study show that the temperature (both minimum and maximum) and rainfall of the island are increasing by ~0.01 ◦ C and ~11.5 mm per year, respectively. The trends of temperature and rainfall are closely associated with vegetative growth; an upward trend in the dense forest is noticed through the enhanced vegetation index (EVI). The trend analysis of satellite altimeter datasets shows that the sea level is increasing at a rate of 6.6 mm/year. The DEM-based modelling shows that sea level rise poses the greatest threat to coastal habitations and has signiﬁcantly increased in recent years, accentuated by urbanisation. The GIS-based model results predict that about half of the coastal settlements (2.5 sq km) will be submerged completely within the next 30 years, provided the same sea level rise continues. The risk of coastal inundation is particularly severe in Sabang, the largest town on the island. The results allow regional, sub-regional, and local comparisons that can assess variations in climate change, sea level rise, coastal inundation, and associated vulnerabilities.


Introduction
The earth's temperature has increased by 0.14 • F (0.08 • C) every ten years since 1880.It has almost doubled (increased by 0.32 • F/0.18 • C) since 1981.Globally, 2020 was one of Remote Sens. 2022, 14, 2857 2 of 14 the warmest years on record (+1.12 • C since 1880) [1], whilst 2021 was the sixth warmest year [2].The surface ocean water is also warming.It has increased by ~0.13 • C every ten years over the last hundred years [3].Sea level has risen by 21 to 24 centimetres since 1880, mostly rising in the last few decades [4].In 2020, the global mean sea level rise reached a new high and had increased 9.13 cm from 1993 levels [5].The compounding effect of warming and sea level rise negatively affects coastal systems and human communities.The Asia-Pacific coastal region, in particular, coastal populations in India, Bangladesh, Vietnam, and Indonesia, are highly exposed to extreme coastal hazards that are predicted to worsen in the upcoming years [6].
Indonesia, Far East Asia, is an island nation with an 81,000 km coastline comprising over 17,500 islands.Several extreme weather events, such as intense rainfall, tropical cyclones, and tsunamis, have occurred in the coastal areas of the islands in recent decades, leading to natural disasters [7,8].In April 2021, Tropical Cyclone Seroja caused intense rainfall (up to 100 mm/day over five days), strong winds, and wave setup which caused over 180 deaths in the East Nusa Tenggara province [9].Despite this, the annual population continues to grow at a rate of ~1.1% [10], and towns and cities are widening along the coastal stretches.The average population density in Indonesian coastal areas (80 persons/sq km) is twice the world's average [11].Indonesians prefer to stay near the coast due to the potential livelihood benefits such as fishing.Climate change is, therefore, a serious issue facing Indonesia, where the projected changes in temperatures, rainfall behaviour, and sea level are expected to negatively affect coastal systems and communities [12][13][14][15].In particular, the vulnerability of coastal communities is increasing in many Indonesian islands due to coastal inundation.According to the IPCC [16], the Indonesian sea level will rise ~100 cm by 2100 due to global, national, and regional climatic changes, and this could ultimately lead to exacerbated coastal floods, erosion, and land subsidence.Rapid urbanisation and land use changes, for example, the expansion of monocultures and deforestation [17,18], on Indonesian islands are expected to increase the disaster risks further [19][20][21].
Consequently, they need to employ local adaptation and mitigation measures against sea level rise, with in situ adaptation measures, such as flood proofing property by raising floors [22,23] or income diversification in regions where fishing is impacted by changing storm patterns, preferred to relocation to the mainland [6,11,24].The coastal zones of Indonesia require increased attention from planners and policymakers to address the predicted climate-driven challenges faced by both coastal ecosystems and socio-economic systems.While the Indonesian coasts provide many economic benefits [25], many human habitations are exposed to a variety of natural hazards, such as cyclones, erosion, saltwater intrusion, subsidence, tsunamis, and floods, due to climate change and sea level rise [26].These hazards cannot be prevented, but the disruption and loss of human life can be mediated by implementing appropriate strategies in advance.Such strategies require a systematic assessment of the impacts of climate change and sea level rise on coastal hazards and inundation.
The literature to date on Indonesian coastal risks has focused on the major cities/islands of the country [27][28][29][30].Many small islands in Indonesia have limited resources and are prone to natural disasters and vulnerable to climate change.Sabang Island attracts tourists from all over the world and has experienced land use degradation, rapid development, and coastal inundation [31].Only a few researchers have analysed Sabang Island's vulnerability to coastal hazards via systematic long-term datasets.This literature is not adequate to understand micro-level climate change and coastal inundation.The current study, therefore, evaluates Sabang Island's vulnerability to climate change and coastal inundation using satellite datasets and GIS tools to fill this research gap.

Case Study Area: Sabang Island
Sabang Island/Weh Island, one of the Indonesian islands, is situated in the southwestern part of Indonesia.Geographically, the island is located between 5 • 46 17 and 5 • 54 18 north latitude and 95 • 12 48 and 95 • 22 39 east longitude, with an areal extent of ~131 sq km and a ~95 km long coastline (Figure 1).Sabang is a cluster of five islands (i.e., Pulau Weh, Rubiah Island, Seulako Island, Pulau Klah, and Pulau Rnondo), of which Pulau Weh (i.e., Weh Island) is the largest island, where Sabang, the largest town, with a population of 33,217 in 2015, is located [32].The daily solar radiation of Sabang Island is high throughout the year.The long-term data  show that the climate of Sabang Island is an equatorial type with hot and wet weather characteristics throughout the year.The island experiences an average temperature range of around 28 • C, with an annual rainfall of 2000 mm.The island is made mostly of volcanic rock, and the topographic variation ranges from mean sea level (MSL) to 610 m (highest at a fumarolic volcano).Dense forest is the dominant land cover type followed by the built-up and agricultural land use classes [33].

Case Study Area: Sabang Island
Sabang Island/Weh Island, one of the Indonesian islands, is situated in the southwestern part of Indonesia.Geographically, the island is located between 5°46′17″ and 5°54′18″ north latitude and 95°12′48″ and 95°22′39″ east longitude, with an areal extent of ~131 sq km and a ~95 km long coastline (Figure 1).Sabang is a cluster of five islands (i.e., Pulau Weh, Rubiah Island, Seulako Island, Pulau Klah, and Pulau Rnondo), of which Pulau Weh (i.e., Weh Island) is the largest island, where Sabang, the largest town, with a population of 33,217 in 2015, is located [32].The daily solar radiation of Sabang Island is high throughout the year.The long-term data  show that the climate of Sabang Island is an equatorial type with hot and wet weather characteristics throughout the year.The island experiences an average temperature range of around 28 °C, with an annual rainfall of 2000 mm.The island is made mostly of volcanic rock, and the topographic variation ranges from mean sea level (MSL) to 610 m (highest at a fumarolic volcano).Dense forest is the dominant land cover type followed by the built-up and agricultural land use classes [33].

Data Sources
The long-term (1981-2020) climatic parameters, including minimum temperature, maximum temperature, rainfall, and relative humidity, were extracted from the National Aeronautics and Space Administration (NASA) POWER Release 8 platform [34].The parameters were extracted from various satellite and model-based products which are sufficiently accurate to provide reliable climate datasets for Sabang Island, where surface measurements were nonexistent for a long period [35].All the cloud-free LANDSAT 5 and 8 satellite datasets available  for Sabang Island were used to compute long-term trends in land use/land cover through enhanced vegetation index (EVI) and

Data Sources
The long-term (1981-2020) climatic parameters, including minimum temperature, maximum temperature, rainfall, and relative humidity, were extracted from the National Aeronautics and Space Administration (NASA) POWER Release 8 platform [34].The parameters were extracted from various satellite and model-based products which are sufficiently accurate to provide reliable climate datasets for Sabang Island, where surface measurements were nonexistent for a long period [35].All the cloud-free LANDSAT 5 and 8 satellite datasets available  for Sabang Island were used to compute long-term trends in land use/land cover through enhanced vegetation index (EVI) and normalised difference built-up index (NDBI).In addition, MODIS datasets were also used to extract EVI for analysing the trends in vegetation during the last two decades.The data on sea level anomalies of the island coasts were extracted using the satellite altimeter datasets obtained from Jet Propulsion Laboratory MEaSUREs gridded sea surface height anomaly (SSHA) dataset [36] for the period 1993-2019.Further, the ALOS PALSAR digital elevation model (DEM) was acquired from Alaska Satellite Facility (ASF) to construct a model for coastal inundation scenarios.Table 1 describes the data sources, parameters extracted, and time window of the satellite datasets.

Assessment of Land Use/Land Cover Change
The spatio-temporal assessment of land use/land cover is essential for understanding the interplay between climate and human society.As the case study location is a small island and consists of mainly the forest and built-up land use/land cover classes, we assessed the dynamics of forest and built-up areas through EVI and NDBI, respectively.The indices were calculated from all available cloud-free datasets from LANDSAT 5 (1988LANDSAT 5 ( -2013)), LANDSAT 8 (2014LANDSAT 8 ( -2021)), and MODIS (2000-2021) datasets of the island using Google Earth Engine (GEE).The median index values for the window periods 1988-1992, 1993-1997, 1998-2001, 2004-2009, 2014-2017, and 2018-2021 were extracted and mapped for visualisation.

Enhanced Vegetation Index (EVI)
EVI is primarily used for analysis and quantifying greenness on a relative scale.It is better than normalised difference vegetation index (NDVI) and more sensitive in an area where vegetation is dense [41].EVI values range between −1 and +1, and, for healthy vegetation, the value varies between 0.2 and 0.8.The EVI can be calculated easily with the help of LANDSAT datasets using the following equation: where R rs (NIR), R rs (Red), and Rrs(blue) represent the reflectance of NIR, red and blue bands respectively; C 1 and C 2 are the coefficients used to correct aerosol scattering with the blue band; L is the soil adjustment factor; and G represents the gain factor.In general, C 1 , C 2 , and L values are assigned as 6.0, 7.5, and 1, respectively.

Normalised Difference Built-Up Index (NDBI)
The NDBI is used to discriminate between built-up areas and other land use classes [42].As built-up areas have higher reflectance in the shortwave-infrared (SWIR) region than the near-infrared (NIR) region, the NDBI can be computed by using the following equation [43]: NDBI values vary from −1 to +1, where the positive indicates a built-up region.

Trend Analysis
Some non-parametric methods, such as Sen's slope estimator and Mann-Kendall test, were used in the study to test the statistical significance of trends in climate, land use, and sea level variables.The MAKESENS template (https://en.ilmatieteenlaitos.fi/makesens(accessed on 10 May 2022)) was used to identify the trend Z score and true slope values (annual changes) in time series data.The Mann-Kendall test and Sen's slope estimations were performed using the monthly and annual average datasets.

Mann-Kendall Test
The non-parametric Mann-Kendall test is widely used to detect the significant, monotonic trend pattern (i.e., upward and downward) in atmospheric and hydrological variables [44].The main advantage of the method is that trends are not affected by minimum missing data and outliers.The Mann-Kendall test for time series x 1 , . . .x n can be calculated using: where n represents the number of samples, x is the value of i = 1, 2, . . ., n -1, and j = i + 1.
If S > 0, then the later observations in the time series tend to be larger than those that appeared earlier in the time series, while the reverse is true if S < 0. When the mean value of S is 0, then the variance of S is calculated by the following equation: The value of t varies according to the rank, and f t represents the frequency of rank t.Then, statistical Z values are computed as follows: where the values of Z represent the presence of a significant trend.If the Z value is greater than 0, there is an increasing trend and vice versa.If Z > Zα/2 (α is a significant level), the trend is significant.The standard normal distribution of α (e.g., 0.001, 0.01, 0.05, and 0.1) is used to find the Zα/2 value.

Sen's Slope Estimator
Sen's slope is a simple non-parametric method for estimating the magnitude of time series data [45].
Sen's estimator of the slope is the median of the N values of Q i , and it is ranked from smallest to largest.Then the Sen's slope confidence intervals (upper and lower) are computed as follows: where N is the total number of pairs of time series elements of x i and x j , and se represents the standard error for the Mann-Kendall test.

Modelling of Coastal Inundation
In this study, we considered how future projected sea level rise and the resulting coastal inundation may affect Sabang Island.We used the present mean high tide level as a reference coastline to model the probable extent of coastal inundation due to extreme storm surge events and future sea levels.The study considered a first-order approximation of the inundated coastline due to sea level rise with the help of a radiometrically terrain-corrected (RTC), high-resolution ALOS PALSAR digital elevation model (DEM) dataset.The model showed coastal inundation as a function of present-day topographic and bathymetric data, with a uniformly increased water level rise added to the present-day mean high water tide line.We used satellite altimeter data to understand the sea level rise.Sen's slope estimator of long-term sea level anomalies showed that the sea level of the island is rising at a rate of 6.6 mm per year.The IPCC reports indicate that the global mean sea level will rise 0.61-1.10m by 2100 [46].By considering the local trends and global projections, we applied water level rise/flooding command in Global Mapper software to simulate the scenarios of a 0.5 m (2050) and 1 m (2100) rise from the present sea level/coastline.In addition, 2, 5, and 10 m increases in water level over the mean high water tide line were also attempted to see where the water would reach maximum levels during extreme storm surge events [47].All the modelled inundation results were delineated and used for further analysis and visualisation.

Vulnerability of Built-Up Areas
Low-lying and near-shore coastal settlements are always at high risk of various coastal hazards [48].The risk of coastal inundation can be indirectly estimated with the help of assessing the extent of built-up areas [49].Built-up areas can be used to understand the exposure of populations to sea level rise and associated coastal inundations.They can also be used to estimate economic losses to coastal hazards.In this study, the median NDBI layers of 2018-2021 were used to prepare the built-up layer of Sabang Island.With the help of GIS, the built-up layer was intersected with all predicted coastal inundation layers to estimate the maximum damage exposure [50].By overlaying the modelled coastal inundation layers on a built-up layer, we mapped the vulnerable areas to coastal inundation.Although the actual risk of each inundation scenario is dependent on many factors, this assessment provided a relative assessment of the maximum potential damage exposure of human habitations to sea level rise.

Climate Change in Sabang Island
The long-term data (1981-2020) show that the average annual temperature of the case study island is 28 • C with a diurnal variation of 25-31 • C. The relative humidity is about 80%, and it is comparatively high September-December.The annual average rainfall of the island is about 2000 mm, and it is well distributed throughout the year.Except for February, all the months contribute more than 100 mm of rainfall, and the contribution is at its maximum in May and October-December (Figure 2).
The Mann-Kendall trend analysis and Sen's slope estimate for the long-term climatic data show that the temperature of the island has a strong upward trend with a level of significance α = 0.01.The mean maximum temperature of all the months is increasing, and the annual rate of increase is 0.01 • C per year (Figure 3A).The rate of increase of the maximum temperature is much more pronounced in January (0.02 • C), with a high level of significance (α = 0.001).The minimum temperature also has an increasing trend (Figure 3B), especially in the wetter part of the year (June-December).The Mann-Kendall tests indicate that the variables of mean maximum and minimum temperatures have a positive Z score throughout all the months, with a strong statistical significance (Figure 4).
The rainfall trends in Sabang Island are cyclic but with an upward trend (Figure 3C).Sen's slope estimate shows that rainfall on the island is increasing at a rate of 11.5 mm per year with a statistically significant level (α = 0.01) of Z score (3.16).The month-wise plotting of Z score analysis showed that the rainfall pattern follows a significant upward trend during the comparatively wet months i.e., November-December (Figure 4).This season also has a statistically significant upward trend in minimum and maximum temperature, and increasing trends of air temperature lead to more convection and rainfall during this season.
The long-term data (1981-2020) show that the average annual temperature of the case study island is 28 °C with a diurnal variation of 25-31 °C.The relative humidity is about 80%, and it is comparatively high September-December.The annual average rainfall of the island is about 2000 mm, and it is well distributed throughout the year.Except for February, all the months contribute more than 100 mm of rainfall, and the contribution is at its maximum in May and October-December (Figure 2).The Mann-Kendall trend analysis and Sen's slope estimate for the long-term climatic data show that the temperature of the island has a strong upward trend with a level of significance α = 0.01.The mean maximum temperature of all the months is increasing, and the annual rate of increase is 0.01 °C per year (Figure 3A).The rate of increase of the maximum temperature is much more pronounced in January (0.02 °C), with a high level of significance (α = 0.001).The minimum temperature also has an increasing trend (Figure 3B), especially in the wetter part of the year (June-December).The Mann-Kendall tests indicate that the variables of mean maximum and minimum temperatures have a positive Z score throughout all the months, with a strong statistical significance (Figure 4).
The rainfall trends in Sabang Island are cyclic but with an upward trend (Figure 3C).Sen's slope estimate shows that rainfall on the island is increasing at a rate of 11.5 mm per year with a statistically significant level (α = 0.01) of Z score (3.16).The month-wise plotting of Z score analysis showed that the rainfall pattern follows a significant upward trend during the comparatively wet months i.e., November-December (Figure 4).This season also has a statistically significant upward trend in minimum and maximum temperature, and increasing trends of air temperature lead to more convection and rainfall during this season.
The relative humidity of the island has a slightly decreasing trend but no statistical significance.The level of decrease in relative humidity is high during February, with a statistical level of significance (α = 0.1).The decreasing trend of relative humidity is associated with air temperature, i.e., both minimum and maximum temperatures have an upward trend (Figure 4).The relative humidity of the island has a slightly decreasing trend but no statistical significance.The level of decrease in relative humidity is high during February, with a statistical level of significance (α = 0.1).The decreasing trend of relative humidity is associated with air temperature, i.e., both minimum and maximum temperatures have an upward trend (Figure 4).

Sea Level Rise
The sea level changes of the island coasts were studied indirectly using the satellite altimeter datasets.The averaged sea surface height anomaly data of January over the last three decades (Figure 5) show that the sea level is increasing at a rate of 6.6 mm/year with a statistical level of significance (α = 0.05).The data show that the higher positive anomalies are significant, especially in the last decade, and such a steep rise in sea level will result in more coastal inundation.

Changes in Land Use/Land Cover
The case study island is dominantly covered by dense forest and built-up land use/land cover classes.The spatio-temporal changes in these dominant classes were studied based on EVI and NDBI indices using LANDSAT datasets (Figures 6 and 7).The results show that dense forest cover follows the pattern of rainfall trends, and has increased in the last decade.The areal extent of dense forest cover based on median EVI data in 1988-1992 was about 92 sq km, and it increased to 115 sq km as per the recent median EVI data (2018-2021).The other dominant land use category, that of built-up land use, gradually increased from 1.6 sq km in 1988 to 5.4 sq km in 2021, with an approximately one-fold increase in built-up cover for every decade (Table 2).

Sea Level Rise
The sea level changes of the island coasts were studied indirectly using the satellite altimeter datasets.The averaged sea surface height anomaly data of January over the last three decades (Figure 5) show that the sea level is increasing at a rate of 6.6 mm/year with a statistical level of significance (α = 0.05).The data show that the higher positive anomalies are significant, especially in the last decade, and such a steep rise in sea level will result in more coastal inundation.The analysis shows that the sea level is significantly increasing (α = 0.01) at a rate of 6.6 mm/year.

Changes in Land Use/Land Cover
The case study island is dominantly covered by dense forest and built-up land use/land cover classes.The spatio-temporal changes in these dominant classes were studied based on EVI and NDBI indices using LANDSAT datasets (Figures 6 and 7).The results show that dense forest cover follows the pattern of rainfall trends, and has increased in the last decade.The areal extent of dense forest cover based on median EVI data in 1988-1992 was about 92 sq km, and it increased to 115 sq km as per the recent median EVI data (2018-2021).The other dominant land use category, that of built-up land use, gradually increased from 1.6 sq km in 1988 to 5.4 sq km in 2021, with an approximately one-fold increase in built-up cover for every decade (Table 2).The analysis shows that the sea level is significantly increasing (α = 0.01) at a rate of 6.6 mm/year.

Coastal Inundation and Built-Up Areas
The coastal inundation on the island was assessed based on the present coastline, topography configuration, rate of sea level rise, and exposure to the built-up area.Although the data points were limited (only four points for the entire island), they were valuable for showing a general upward pattern of sea level with a rate of 6.6 mm/year.The water rise model was constructed based on the present coastline and DEM for two scenarios: (1) probable future sea levels, i.e., 0.2 m in 2050 and 1 m in 2100, and (2) extreme storm surge/tsunami events, i.e., a 5 m and 10 m rise of wave heights.Although the island's slope is generally steep, the modelled results signify that about half of the coastal settlements (2.5 sq km) will be submerged due to a rise in sea level of just 0.2 m from the present coastline (Table 3).The vulnerability of coastal inundation is very significant and severe in Sabang town, the largest cluster of settlements on the island (Figure 7), due to the growth of the built-up area over the coastal plain topography (Figure 8).The other important coastal towns that have greater vulnerability to coastal inundation are Sukajaya and Ganpong Iboih.These towns are not only exposed to the risk of sea level rise but are also prone to possible storm surge/tsunami events.ments (2.5 sq km) will be submerged due to a rise in sea level of just 0.2 m from the present coastline (Table 3).The vulnerability of coastal inundation is very significant and severe in Sabang town, the largest cluster of settlements on the island (Figure 7), due to the growth of the built-up area over the coastal plain topography (Figure 8).The other important coastal towns that have greater vulnerability to coastal inundation are Sukajaya and Ganpong Iboih.These towns are not only exposed to the risk of sea level rise but are also prone to possible storm surge/tsunami events.

Discussion
Rising sea levels and their subsequent effect on the coastal zone are a concern at a global level.The recent acceleration in sea level rise due to anthropogenic-induced climate change means that the negative impacts of climate-change-induced hazards are an imminent threat to coastal communities worldwide.Few studies have assessed the impact of sea level rise on low-lying parts of the Indonesian islands using historic observation

Discussion
Rising sea levels and their subsequent effect on the coastal zone are a concern at a global level.The recent acceleration in sea level rise due to anthropogenic-induced climate change means that the negative impacts of climate-change-induced hazards are an imminent threat to coastal communities worldwide.Few studies have assessed the impact of sea level rise on low-lying parts of the Indonesian islands using historic observation records from tide gauge data and sea level anomaly altimetry satellite data [51] and secondary data integrated to form layers in the geospatial environment for inundation modelling [25,28].Further, quantifying the impact of climate change and sea level rise on the coastal environment and ecosystems of Indonesian islands has only been attempted with coarser spatial and temporal resolution datasets.Comprehensive vulnerability assessments are required to understand the impacts of climate-change-induced hazards on small, low-lying islands to inform adaptation strategies, mitigate losses, and improve the resilience of communities to future change [52].
In previous studies, the quantification of potential risk zones was attempted, but they utilised single-term data.Rather, in this study, we utilised the long-term change of multiple environmental factors, such as minimum temperature, maximum temperature, rainfall, and forest cover, and also population dynamics for the delineation of potential risk zones for coastal inundation due to relative sea level rise.A long-term dataset was considered (1988-2021) for climatic and land use variables.It was observed that the temperature (both minimum and maximum) and rainfall of the island showed a strong upward trend (α = 0.01).The Sen's slope estimate indicated that the maximum temperature and rainfall of the island are increasing by about 0.01 • C and 11.5 mm, respectively, per year.It was also observed that the positive trends of temperature and rainfall are closely associated with vegetative growth, and an upward trend in the dense forest was noticed through the enhanced vegetation index (EVI).The sea level dynamics of the island coasts were studied indirectly using the satellite altimeter datasets, and the results show that the sea level is increasing at a rate of 6.6 mm/year with a statistical level of significance (α = 0.05).The positive sea level anomalies are significant in recent years, resulting in more coastal inundation in the low-lying areas.From the NDBI change detection analysis, it was observed that the low-lying coastal areas are being converted into built-up areas, expanding three-fold since 1988, and possess greater risks of coastal inundation.The coastal inundation model utilising the ALOS PALSAR DEM indicated that about half of the coastal settlements (2.5 sq km) would be submerged completely within the next 30 years if the same sea level rise continues.
Through this study, we proposed an effective spatio-temporally transferable methodological framework for the impact assessment of coastal flooding in smaller islands.This study could be improved by including the in situ measurements of tidal variations and a focused group survey to understand the perspective of the local community on climate change and the prevailing issues due to sea level rise.Improved observation records of physical phenomena, including tides, pressure, rainfall, wind, and waves, would also allow for dynamic flood modelling to better understand which areas may be most vulnerable to climate-change-induced hazards [27,53,54].

Conclusions
In this study, we chose Sabang Island, Indonesia, to assess climate change, sea level rise, and coastal inundations.Here, we estimated the trends in climatic variables and sea level rise by utilising long-term (1988-2021), satellite-derived data products and statistical analysis.It was observed that the maximum temperature and rainfall of the island are increasing by about 0.01 • C and 11.5 mm, respectively, per year, with a statistical significance (α = 0.01).The satellite altimeter datasets showed that the sea level is increasing at a rate of 6.6 mm/year.The topographic geospatial modelling of sea level rise scenarios based on IPCC projections and local sea level trends showed that coastal inundation in low-lying areas is a serious issue on the island.About half of the coastal settlements on the islands (2.5 sq km) will be submerged completely within the next 30 years if the same trends in sea level rise continue.The study formed baseline datasets for understanding vulnerability to coastal hazards and gives an insight for stakeholders and policymakers to frame policies for sustainable coastal zone management.Further, the methodological framework we proposed through this study is spatio-temporally transferable and also utilised mainly open-source data.So, the study could be extended to other smaller islands of this region on a larger scale with reasonable accuracy.

Figure 2 .
Figure 2. Climatic diagram of Sabang Island based on long-term climate data (1981-2020).The bar graphs denote the month-wise average amount of rainfall (in mm).The dark-red line represents the average temperature (in °C), and the dotted, blue line shows monthly variations of relative humidity (%).

Figure 2 .
Figure 2. Climatic diagram of Sabang Island based on long-term climate data (1981-2020).The bar graphs denote the month-wise average amount of rainfall (in mm).The dark-red line represents the average temperature (in • C), and the dotted, blue line shows monthly variations of relative humidity (%).Remote Sens. 2022, 14, x FOR PEER REVIEW 8 of 15

Figure 3 .
Figure 3. Climatic trends in Sabang Island (1981-2020) of (A) annual average maximum temperature, (B) annual average minimum temperature, (C) annual average rainfall, and (D) annual average relative humidity.A statistically significant upward trend is noticed in maximum and minimum temperatures.The rainfall pattern follows a cyclic pattern with an overall increasing trend.The relative humidity follows a long-term (~15 years) cyclic trend.

Figure 3 .
Figure 3. Climatic trends in Sabang Island (1981-2020) of (A) annual average maximum temperature, (B) annual average minimum temperature, (C) annual average rainfall, and (D) annual average relative humidity.A statistically significant upward trend is noticed in maximum and minimum temperatures.The rainfall pattern follows a cyclic pattern with an overall increasing trend.The relative humidity follows a long-term (~15 years) cyclic trend.

Figure 3 .
Figure 3. Climatic trends in Sabang Island (1981-2020) of (A) annual average maximum temperature, (B) annual average minimum temperature, (C) annual average rainfall, and (D) annual average relative humidity.A statistically significant upward trend is noticed in maximum and minimum temperatures.The rainfall pattern follows a cyclic pattern with an overall increasing trend.The relative humidity follows a long-term (~15 years) cyclic trend.

Figure 4 .
Figure 4. Monthly and annual distribution of Mann-Kendall test Z score.The long-term data (1981-2020) show that the maximum and minimum temperature (red symbols) of Sabang Island has an upward trend for all the months with a statistically significant Z score.Except for May and September, rainfall (blue symbols) also has an upward trend for all the months.Relative humidity (green symbols) shows an upward trend during the months April-August.The annual average data show

Figure 4 .
Figure 4. Monthly and annual distribution of Mann-Kendall test Z score.The long-term data show that the maximum and minimum temperature (red symbols) of Sabang Island has an upward trend for all the months with a statistically significant Z score.Except for May and September, rainfall (blue symbols) also has an upward trend for all the months.Relative humidity (green symbols) shows an upward trend during the months April-August.The annual average data show that maximum-minimum temperature and rainfall follow an upward trend with a statistically significant Z score, while relative humidity shows a downward trend with an insignificant Z score.

Figure 5 .
Figure 5. Trend of sea surface height anomaly (SSHA) on the coasts of Sabang Island (1993-2019).The analysis shows that the sea level is significantly increasing (α = 0.01) at a rate of 6.6 mm/year.

Figure 5 .
Figure 5. Trend of sea surface height anomaly (SSHA) on the coasts of Sabang Island (1993-2019).The analysis shows that the sea level is significantly increasing (α = 0.01) at a rate of 6.6 mm/year.

15 Figure 6 .
Figure 6.Median values of enhanced vegetation index (EVI) from 1988 to 2021.The dark green colour denotes dense vegetation.

Figure 6 .
Figure 6.Median values of enhanced vegetation index (EVI) from 1988 to 2021.The dark green colour denotes dense vegetation.

Figure 6 .
Figure 6.Median values of enhanced vegetation index (EVI) from 1988 to 2021.The dark green colour denotes dense vegetation.

Figure 7 .
Figure 7. Median values of normalised difference built-up index (NDBI) from 1988 to 2021.The built-up areas increased gradually from 1988 to 2021, and the changes are apparent after 2001.

Figure 7 .
Figure 7. Median values of normalised difference built-up index (NDBI) from 1988 to 2021.The built-up areas increased gradually from 1988 to 2021, and the changes are apparent after 2001.

Figure 8 .
Figure 8. Predicted coastal inundation on Sabang Island.The blue shades represent different water level rise scenarios.The latest built-up area is overlaid to allow an understanding of the built-up elements at risk.

Figure 8 .
Figure 8. Predicted coastal inundation on Sabang Island.The blue shades represent different water level rise scenarios.The latest built-up area is overlaid to allow an understanding of the built-up elements at risk.

Table 1 .
Sources of satellite datasets used in the study.

Table 2 .
Dynamics of built-up areas and dense forests of Sabang Island.

Table 2 .
Dynamics of built-up areas and dense forests of Sabang Island.

Table 3 .
The area of built-up elements at risk under different predicted coastal inundation levels.

Table 3 .
The area of built-up elements at risk under different predicted coastal inundation levels.
Predicted Coastal Inundation (meters)Built-Up Are (sq km)