Remote Sensing to Study Mangrove Fragmentation and Its Impacts on Leaf Area Index and Gross Primary Productivity in the South of Peninsular Malaysia

: Mangrove is classiﬁed as an important ecosystem along the shorelines of tropical and subtropical landmasses, which are being degraded at an alarming rate despite numerous international treaties having been agreed. Iskandar Malaysia (IM) is a fast-growing economic region in southern Peninsular Malaysia, where three Ramsar Sites are located. Since the beginning of the 21st century (2000–2019), a total loss of 2907.29 ha of mangrove area has been estimated based on medium-high resolution remote sensing data. This corresponds to an annual loss rate of 1.12%, which is higher than the world mangrove depletion rate. The causes of mangrove loss were identiﬁed as land conversion to urban, plantations, and aquaculture activities, where large mangrove areas were shattered into many smaller patches. Fragmentation analysis over the mangrove area shows a reduction in the mean patch size (from 105 ha to 27 ha) and an increase in the number of mangrove patches (130 to 402), edge, and shape complexity, where smaller and isolated mangrove patches were found to be related to the rapid development of IM region. The Moderate Resolution Imaging Spectro-radiometer (MODIS) Leaf Area Index (LAI) and Gross Primary Productivity (GPP) products were used to inspect the impact of fragmentation on the mangrove ecosystem process. The mean LAI and GPP of mangrove areas that had not undergone any land cover changes over the years showed an increase from 3.03 to 3.55 (LAI) and 5.81 g C m − 2 to 6.73 g C m − 2 (GPP), highlighting the ability of the mangrove forest to assimilate CO 2 when it is not disturbed. Similarly, GPP also increased over the gained areas (from 1.88 g C m − 2 to 2.78 g C m − 2 ). Meanwhile, areas that lost mangroves, but replaced them with oil palm, had decreased mean LAI from 2.99 to 2.62. In fragmented mangrove patches an increase in GPP was recorded, and this could be due to the smaller patches (<9 ha) and their edge effects where abundance of solar radiation along the edges of the patches may increase productivity. The impact on GPP due to fragmentation is found to rely on the type of land transformation and patch characteristics (size, edge, and shape complexity). The preservation of mangrove forests in a rapidly developing region such as IM is vital to ensure ecosystem, ecology, environment, and biodiversity conservation, in addition to providing economical revenue and supporting human activities. results with fragmentation analysis and mangrove LAI and GPP data to reveal the impact of fragmentation on mangrove productivity in a region that hosts important mangrove (Ramsar) sites in Malaysia. Such results are expected to provide standard guidelines for a stronger policy formulation to protect the remaining mangroves in the region and in other mangrove areas. Further studies on how to prevent mangrove loss, how to reconnect the fragmented mangrove forests, and how to improve the carbon stock estimation to prove the value of mangrove to the ecosystem processes are equally important, now and in the future.


Introduction
Changes in the extent of mangrove cover has resulted in fragmentation that can greatly affect its functions and services to the environment and human kind [1]. Fragmentation of mangrove areas refers to the sub-division of a continuous landscape into smaller units, Remote sensing and Geographic Information System (GIS) are cost-and time-effective technologies that can be used to support large scale mangrove mapping, monitoring its changes [34] and analyzing the landscape's fragmentation [1,2,29]. The launch of the Landsat mission promoted the creation of global and regional mangrove maps at 30 m spatial resolution. These maps can be further improved using freely available satellite images with better spatial resolution, such as that of Sentinel-2. Digital image classification algorithms have been widely developed for obtaining land use and land cover using remotely sensed images. Pixel-based classification involving the unsupervised or supervised method are very common and useful. In the last decade, object-oriented, neural network, decision tree, and machine learning based classification algorithms, with advantages of improved accuracy and reduced computational time, have become more popular. Reviews by Kuenzer et al. [35] and Pham, Xia, Ha, Bui, Le, and Takeuchi [25] concluded that most of the classification algorithms can give overall accuracy above 80%, while, with improvements such as training samples selection, satellite data quality etc., a higher classification accuracy (~90%) can be achieved. Better classification results can improve the mangrove extent mapping and the change detection analysis to accurately detect the gain and loss of mangroves.
This study aims to describe the mangrove loss and fragmentation in the Iskandar Malaysia region over the 2000-2019 period and analyze the loss and fragmentation effects on mangrove LAI and GPP values. More specifically, we intend to (1) document the loss and gain of mangroves over a period of two decades (2000-2019), (2) analyze mangrove fragmentation, and (3) relate and analyze the loss/fragmentation impacts on the LAI and GPP values in order to define the function of mangroves to assimilate CO 2 in relation to mangrove cover changes and fragmentation in the study area. It is hypothesized that a higher loss of mangrove canopy cover and fragmented areas will decrease the LAI and GPP of the ecosystem.

Study Area
Malaysia constitutes only 4% of global mangrove cover, but at least 70 true and associated species from 28 families have been identified in the country, and the majority of the mangrove trees are from the Rhizophoraceae family [36]. More than half of Malaysian mangroves are distributed in the state of Sabah (60%), 22% in Sarawak, and the remaining 18% are found in Peninsular Malaysia [36]. In Peninsular Malaysia, the majority of mangroves are located along the west coast, and the state of Johor alone holds 30% of the total mangroves ( Figure 1). Three of the seven Ramsar Sites in Malaysia, namely Pulau Kukup (647 ha), Tanjung Piai (526 ha), and Sungai Pulai (9126 ha) are located in the Iskandar Malaysia region in Johor (Johor National Parks, https://www.johornationalparks.gov.my/ v3/RAMSAR-site/, assessed on 12 December 2020).
Iskandar Malaysia (IM) is a fast-growing national special economic region with a total area of 221,700 ha, located in the southern part of Johor state, Peninsular Malaysia ( Figure 1). Specific zones where rapid development has taken place in IM are the Port of Tanjung Pelepas, Tanjung Bin Power Plant, the Port of Pasir Gudang, and the rapid reclamation zone of Forest City, as shown in Figure 1. These activities are mainly caused by population growth that requires more land to be developed for housing establishments, transition to agricultural lands (particularly oil palm), and the building of infrastructure such as ports, dams, and waterfront cities may also lead to fragmentation and impact the carbon sink potential and primary productivity of the mangrove forests in IM. Thus, it is timely to investigate the impact of fragmentation on the primary productivity of the mangrove ecosystem. Although sustainable development policies are employed in the region, the decrease of mangrove forests was found to be at an alarming rate (33%) between 2005 and 2014 [32]. Nevertheless, the replantation and conservation efforts since 2014 may have also increased mangrove coverage; thus, it is important to monitor the changes and analyze their impact on mangrove services. the decrease of mangrove forests was found to be at an alarming rate (33%) between 2005 and 2014 [32]. Nevertheless, the replantation and conservation efforts since 2014 may have also increased mangrove coverage; thus, it is important to monitor the changes and analyze their impact on mangrove services.

Data and Methodology
The methodology flow-chart ( Figure 2) for this study and details on the data processing and the methods to achieve the aim of the study are described in the following sections. This study utilized medium-resolution remote sensing satellite imageries to monitor changes in mangrove cover and fragmentation, covering the years 2000 to 2019. Support Vector Machine (SVM) was selected as the classifier for the image classification process. A patch analysis plugin was implemented for fragmentation analysis based on ArcGIS software. The land cover, LAI, and GPP Products of Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 006 were used to define the impact of loss and fragmentation on mangrove forests in IM. The MODIS Land Cover Product is based on a supervised decision-tree classification method using inputs from other MODIS products, including the land/water mask, nadir bidirectional reflectance distribution function (BRDF)-adjusted reflectance, Enhanced Vegetation Index (EVI), snow cover, land surface temperature, terrain elevation information, etc. [37]. The MODIS LAI product incorporated the MODIS Land Cover Product and MODIS Surface Reflectance Product into the radiative transfer model to derive LAI [38]. The MODIS GPP product is derived by relating the absorbed photosynthetically active radiation (APAR) calculated from LAI with the energy conversion efficiency (ε) [26]. More information of the MODIS products and the method to access these products are described below.

Data and Methodology
The methodology flow-chart ( Figure 2) for this study and details on the data processing and the methods to achieve the aim of the study are described in the following sections. This study utilized medium-resolution remote sensing satellite imageries to monitor changes in mangrove cover and fragmentation, covering the years 2000 to 2019. Support Vector Machine (SVM) was selected as the classifier for the image classification process. A patch analysis plugin was implemented for fragmentation analysis based on ArcGIS software. The land cover, LAI, and GPP Products of Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 006 were used to define the impact of loss and fragmentation on mangrove forests in IM. The MODIS Land Cover Product is based on a supervised decision-tree classification method using inputs from other MODIS products, including the land/water mask, nadir bidirectional reflectance distribution function (BRDF)-adjusted reflectance, Enhanced Vegetation Index (EVI), snow cover, land surface temperature, terrain elevation information, etc. [37]. The MODIS LAI product incorporated the MODIS Land Cover Product and MODIS Surface Reflectance Product into the radiative transfer model to derive LAI [38]. The MODIS GPP product is derived by relating the absorbed photosynthetically active radiation (APAR) calculated from LAI with the energy conversion efficiency (ε) [26]. More information of the MODIS products and the method to access these products are described below. . This data is registered on the Worldwide Reference System-2 (WRS-2), and one tile (Path 125, Row 59) is required to cover the whole Iskandar Malaysia zone. A sentinel-2 Level-1C multispectral image dated June 28, 2019, was chosen to present the latest conditions of the study area because no Landsat image with acceptable cloud cover was available. The image was downloaded from the Copernicus Open Access Hub (https://scihub.copernicus.eu/ (accessed on 15 March 2021)). For both Landsat 7 and Sentinel-2, only images with less than 20% cloud cover were used. The Semi-Automatic Classification plug-in (SCP) provided in the QGIS software [39] was used to perform the preprocessing of Landsat 7 and the Sentinel-2 images.
Preprocessing, including radiometric correction and atmospheric correction based on the Dark Object Subtraction 1 (DOS 1) model and pan-sharpening model, was applied to the Landsat 7 image. The radiometric and atmospheric correction transferred the digital numbers (DN) to radiance then to reflectance. The pan-sharpening function, based on intensity computation, fused the 30 m Landsat multispectral bands (6 bands covering the visible (blue, green, red), near infrared, and shortwave infrared) with the panchromatic band to produce a 15 m spatial resolution reflectance image. The Sentinel-2 image underwent similar preprocessing steps (radiometric and atmospheric correction using DOS 1) based on the SCP. The visible bands (band 2-blue, 3-green, 4-red) and near infrared (band 8) of 10 m spatial resolution were stacked in this case. The output was a 4-band multispectral reflectance image at 10 m spatial resolution. A thresholding method based on the reflectance of the blue band (>0.25 for Landsat 7 and >0.21 for Sentinel-2) was applied to detect and mask out any cloud on the image. Note that, for both Landsat and Sentinel images, only a very low cloud cover was detected over the mangrove areas, which is the main focus of this study.
Three latest MODIS products (Collection 006) at 500 m spatial resolution were downloaded, namely the MOD15A2H (Leaf Area Index, LAI) [38], the MOD17A2H (Gross Pri- Preprocessing, including radiometric correction and atmospheric correction based on the Dark Object Subtraction 1 (DOS 1) model and pan-sharpening model, was applied to the Landsat 7 image. The radiometric and atmospheric correction transferred the digital numbers (DN) to radiance then to reflectance. The pan-sharpening function, based on intensity computation, fused the 30 m Landsat multispectral bands (6 bands covering the visible (blue, green, red), near infrared, and shortwave infrared) with the panchromatic band to produce a 15 m spatial resolution reflectance image. The Sentinel-2 image underwent similar preprocessing steps (radiometric and atmospheric correction using DOS 1) based on the SCP. The visible bands (band 2-blue, 3-green, 4-red) and near infrared (band 8) of 10 m spatial resolution were stacked in this case. The output was a 4-band multispectral reflectance image at 10 m spatial resolution. A thresholding method based on the reflectance of the blue band (>0.25 for Landsat 7 and >0.21 for Sentinel-2) was applied to detect and mask out any cloud on the image. Note that, for both Landsat and Sentinel images, only a very low cloud cover was detected over the mangrove areas, which is the main focus of this study.
Three latest MODIS products (Collection 006) at 500 m spatial resolution were downloaded, namely the MOD15A2H (Leaf Area Index, LAI) [38], the MOD17A2H (Gross Primary Production, GPP) [26], and the MCD12Q1 (Land cover type, LC) [37]. Only one tile (H28.V08) is required to cover the whole IM region. Both LAI and GPP are 8-day composite products. of the LAI and GPP products, pixels with a QC bit-converted value below 48 were kept and extracted, referring to the MODIS product user guides [26,38,40]. The products were re-projected into WGS 84 projection, and the QC considered the sensor, detectors, and cloud condition issues of the product. For the LC product, we selected the International Geosphere-Biosphere Program (IGBP) layer for this study. According to the IGBP layer, the permanent wetlands (PW) pixels (LC code 11) were identified. The corresponding LAI and GPP values located in the identified PW pixels were extracted with the strict consideration of the QC. The mean and standard deviation of the LAI and GPP at 8-day intervals were computed for further analysis.

Image Classification and Assessment
Support Vector Machine (SVM) classification has been one of the most applied classifiers in recent years. SVM uses the machine learning technique for multispectral and hyperspectral remotely sensed data classification. SVM relies on n-dimensional spectral space to create a hyperplane and identify and separate the classes based on user optimized parameters, including the selection of kernel type and function and the penalty parameter [41]. A total of 7 land use/land cover classes were detected and discriminated from the Landsat and Sentinel images based on the SVM, as follows: terrestrial forest, mangrove forest, oil palm plantation, rubber plantation, urban, water, and others (orchard, shrub, bush, aquaculture farms, etc.). The details of the parameters needed to run the SVM classifier are listed in Table 1. The classification was run based on a total of 140 training data (20 polygons per class with 50-60 pixels per polygon), where different sets of training data were selected during Landsat 7 and Sentinel 2 classification. Figure A1 in the Appendix A shows the training polygons for the Landsat and Sentinel classification. The accuracy assessment of the classified land covers/uses was based on the sample points selected from the Google Earth images as the "ground truth" samples. Google Earth is notably one of the most complete and useful datasets with acknowledged accuracy, and is also freely available for various types of research, including classification and its validation [42][43][44]. Classification assessment based on the error matrix was performed. The overall accuracy (OA) is the ratio between the correctly classified pixels to the total number of reference pixels. The producer accuracy (PA) is the accuracy of the map from the perspective of the map maker, computed by dividing the correctly classified samples of the class to the sum of the reference samples for the class. The user accuracy (UA) is the accuracy from the point of view of a map user, calculated by dividing the correct classification pixels to the total number of pixels mapped as the class. The kappa coefficient (KC) is a value computed based on the KHAT statistic to measure the agreement and accuracy of the classification. In addition, the Z-statistic has been computed to describe the significance of the classification results. The function of the Z-statistic is to determine if the agreement between the classification map significantly agrees with the reference data [45].

Mangrove Cover Change and Fragmentation
After the accuracy assessment of the 2000 and 2019 classification results, the classified mangrove areas were extracted from the classified map. Because the classification maps are in raster form, we converted the maps into vector polygons using the Raster to Polygon function in the ArcGIS software. Change detection was performed based on the polygons using ArcGIS software to detect the unchanged, gained, and lost mangrove areas in IM The Patch Analyst extension in ArcGIS 10.5 was used as the tool for fragmentation analysis. Patch Analyst is based on the FRAGSTATS [46], which can be used for analysis and to provide attributes of the fragmented patches. Several indices were selected in this study to represent the fragmentation of the mangrove area in IM, including Class Area  Table 2. Because mangrove is the focus of the study, the fragmentation analysis was performed for mangrove areas only. Table 2. List of metrics in the Patch Analyst extension used for fragmentation analysis [46]. LAI and GPP are both useful for studying the coverage, health, and potential of carbon assimilation and storage [47]. Mangrove forests have been proven as the highest carbon stock medium; nevertheless, the role of mangrove forest in carbon assimilation and its processes from local to global scale is less understood. This is due to the difficulties of data collection in the mangrove forests. Plant leaves through the photosynthesis process convert the atmospheric CO 2 into nutrients and biomass to support plant growth. The use of LAI to estimate the ecosystem carbon uptake by linking it with GPP is common [48,49]. Remote sensing products were applied in this study to infer the leaf density and productivity of mangrove areas, aimed at detecting the impact of mangrove changes to the products (LAI and GPP) relevant to carbon stock.

Group Fragmentation Index Definition
To identify the impact of mangrove loss and fragmentation over years, the MODIS LAI and GPP of the permanent wetlands (PW) (based on MODIS IGBP land cover classification) were extracted. LAI can be used to describe the leaf density and therefore the vegetation condition/health/ photosynthesis rate and light use efficiency of the mangrove trees [33,34]. We chose LAI as the vegetation index because it is often widely applied to high leaf and tree density zones, including mangrove forests [50], and LAI has been proven to have a strong relationship with Normalized Difference Vegetation Index (NDVI) at different spatial resolutions [51,52]. Meanwhile, MODIS GPP is an indication of the amount of CO 2 assimilated by the ecosystem; thus, it can be used to identify carbon gain/loss due to any changes occurring in the ecosystem. Here, we used GPP rather than NPP because (1) GPP has more direct relations to the photosynthesis rate, leaves, and biomass growth, (2) NPP requires the input of respiration (data is very hard to acquire and not available in the study site), and (3) MODIS only distribute the NPP product at an annual timespan, while we aim to compare GPP and LAI at finer temporal scales (8-days or gain had occurred in LAI or carbon intake by mangroves (to the selected permanent wetlands pixels).

Classification Accuracy and Mangrove Cover Changes
Using the SVM classifier, seven land use/land cover classes, including forest, oil palm, mangrove, rubber, urban, water, and "others", were classified. The class "others" includes orchard, shrub, bush, farm, and aquaculture. Table 3 reports the accuracy of the classification from the two satellite images, including the producer accuracy (PA), user accuracy (UA), overall accuracy (OA) and the kappa coefficient (KC). The overall accuracy of both images is~90%, with a KC of 0.85 and 0.88 for the 2000 and 2019 classifications, respectively. Both the Landsat 7 and Sentinel-2A classification results showed a high accuracy of around 90%, which may be because of their finer spatial resolution and multi spectral bands, and is also due to the machine learning technique of SVM. On the other hand, the pan-sharpened Landsat image also maintained high classification accuracy when more bands (more information) were considered during the classification process. Zstatistics were computed for both classification results, where both values were higher than 1.96, indicating that the classification is meaningful and significant and the classification results are better than random classification. Because the general land use land cover changes are similar to the previous findings in Kanniah, Sheikhi, Cracknell, Goh, Tan, Ho, and Rasli [32], we omitted the land use maps in this study and focused only on the mangrove cover in the following results and discussions. In this study, we further found an expansion of the urban area, especially at the eastern part of Sungai Pulai, where the new Nusajaya township is located, and in the central part of IM, where the Kulai-Senai township has expanded due to increased population. Rubber conversion to oil palm plantation and urban area is still found (in 2019 image) surrounding the Sultan Iskandar Reservoir in the eastern part of IM, showing similar results to [32] in their 2014 satellite data.
The use of the SVM classifier showed great capability by clearly distinguishing mangroves from forests and water bodies, which often can be easily misclassified. Both the user and producer accuracies are high for mangrove classification in the 2000 and 2019 images (Table 3). From the classification results, the loss of mangrove area is noticeable, particularly where it has been converted to urban area and oil palm plantations. Interestingly, the conversion of mangrove to aquaculture land use (classified as "others") was found at multiple locations, which showed the recent development and expansion of the economically valuable aquaculture activities in the IM region (Table 4 and Figure 4b). The performance of SVM classification in this study supports the accuracy of mangrove loss and gain detection (Table 4 and Figure 4), as well as clearly identifying fragmented mangrove patches (Table 6). Based on the classification results, the total mangrove areas in 2000 and 2019 were 13,612.25 ha and 10,704.96 ha, respectively. Figure 3 shows the mangrove areas, their gain, and their loss, and unchanged mangrove areas between 2000 and 2019. The total unchanged mangrove areas over the 19-year period was 8483.44 ha, and the total loss was 5128.80 ha. Despite the loss, some new (2171.38 ha) mangrove areas were detected. This means that a net loss of 2907.29 ha (21.35%) has occurred in IM over the last 19 years. This is equal to a 1.12% annual mangrove loss in the IM region, which is higher than Malaysia's mean depletion rate of 0.41%, based on the estimation of Hamilton and Casey [53]. The southwestern part of the IM, where the Pulau Kukup, Tanjung Piai, and Sungai Pulai mangrove forests (Ramsar sites) are located (Figure 3c), faced significant loss over the study period. On the other hand, at the eastern side of IM, new patches of mangroves were noticed. Figure 4a shows the pinned locations of gain and loss of mangroves over the IM region. Table 4 lists the locations and types of transition from (loss) and to (gain) mangrove area, with Table 5 showing the conditions of some selected loss and gain areas, including when these transitions occurred. During the identification of the gain and loss locations, we also identified around 40 sites in IM that were used for aquaculture activities (Figure 4b).
were noticed. Figure 4a shows the pinned locations of gain and loss of mangroves ov the IM region. Table 4 lists the locations and types of transition from (loss) and to (ga mangrove area, with Table 5 showing the conditions of some selected loss and gain are including when these transitions occurred. During the identification of the gain and lo locations, we also identified around 40 sites in IM that were used for aquaculture activit (Figure 4b).

Fragmentation Analysis
Fragmentation results in Table 6 show that the IM region underwent significant fragmentation over the two decades of the study period. Since we considered only one class (mangrove) in the fragmentation analysis, the area indices, including class area (CA) and total landscape area (TLA), are the same, which is 13,612 ha and 10,705 ha for years 2000 and 2019, respectively. There was a significant loss of mangrove area (21%) during the study period; thus, the analysis of the remaining fragmentation metrics will also be focused on the influence of the area loss on each metric.                The patch size index provides information about the Number of Patches (NP) and the Mean Patch Size (MPS). The patch size is one of the fragmentation parameters that has been studied most [17]. Increased NP from 130 to 402 and decreased MPS from 104.71 to 26.63 (reduction of −78.08 ha) between 2000 and 2019 indicates that the mangrove areas in the IM region had been broken down into smaller patches. The changes of minimum patch size (3.08 ha to 0.001 ha) and maximum patch size (6883.01 ha to 4067.96 ha) from years 2000 to 2019 showed a similar trend to the MPS values, indicating that the mangrove areas had been fragmented into smaller patches throughout the entire study area. In a previous study, Bryan-Brown, Connolly, Richards, Adame, Friess, and Brown [1] reported a smaller reduction of mean patch size of −7.2 ha for Malaysia between 2000 and 2012, which infers that the fragmentation in the IM region is significantly more serious than the country's average. The higher MPS could also be due to the more rapid development in IM compared to other cities located near to mangrove forests in other parts of Malaysia. Coincidentally, the largest mangrove patch in the year 2019 (4067.96 ha) was found located in the largest mangrove patch of the year 2000 (6883.01 ha). Analysis of this large mangrove patch showed that it had fragmented into 60 smaller patches (0.001 ha to 4067.96 ha, with a mean of 94.29 ha), with a loss of 1131.63 ha over the study period. This specific patch has been selected for detailed fragmentation analysis on its effect of LAI and GPP in the later part of this manuscript.

Fragmentation Analysis
Similarly, the Edge indices (TE and ED) also increased from 2000 to 2019, primarily because of the separation of large patches into small patches. Figure A2 in the m shows the fragmentation of mangrove over a selected area within the study area. The increases in TE (from 1200.5 km to 1365 km) and ED (from 88.2 m/ha to 127.5 m/ha) demonstrate that the fragmentation is a result of man-made disturbances rather than natural effects. This can be explained by natural shaped patches normally having fewer edges than man-made urban areas and agriculture boundaries. The values of patch size and edge indices obtained in this study show typical mangrove forest fragmentation over a rapidly developing region, where the disturbance of a large mangrove area into smaller patches occurred due to the transformation of mangroves into urban, agriculture, and aquaculture lands. Generally, the decrease in any forest habitat area or MPS has been shown to create fragments of different sizes (large and small) and isolation among the fragmented patches. Hence, this can result in species extinction and biodiversity loss due to decreased colonization and population size [17,54].
The Shape metrics (MSI, AWMSI, and MPAR) explain the changes of the mangrove areas shape from regular to more complex patches. An MSI of 1 refers to a circular shape, and higher values indicate more irregular shape patterns. The MSI calculated in this study changed from 2.79 in 2000 to 1.96 in 2019, meaning that the complexity of the shape of the patches was reduced. This phenomenon indicates that the mangrove patches change from natural shaped patches (complex in shape) to man-made shape (regular in shape), due to anthropogenic activities [46,55]. The AWMSI refers to the area weighted MSI, which highly depends on the area of patches and the NP. Slightly lower AWMSI was found in 2019 compared to 2000 due to a decrease in the sum of patch area and a large increase of NP. This indicates that the MSI was divided by a relatively larger value in 2019 than in 2000, meaning that there were more fragmented patches in the region in 2019. MPAR, which considers the mean perimeter-to-area ratio shows an increased value for 2019 compared to 2000. Smaller patches with decreasing areas and relatively longer perimeters are the main reason for the large increase in MPAR. Finally, the MPFD, which has values between 1 (simple shaped) to 2 (complex shaped) have increased between 2019 and 2000, showing that the fragmented patches of mangrove in 2019 were more irregular compared to larger mangrove patches in 2000.  Table 7 summarizes the mean and standard deviation values of LAI and GPP for these cases. Note that there are some missing data in these graphs due to unsuccessful retrieval of the MODIS products on certain dates that could be due to cloud cover, failure of detectors, or failure of the main radiative transfer algorithm used to retrieve the LAI/GPP values. The analysis starts from day of year (DOY) 57 because the MODIS LAI and GPP products are only available after February 2000.      Figure 6a shows the mean LAI extracted from the MODIS 8-day product overlaid with all of the PW pixels found in the classified mangrove area. The yearly mean LAI was 3.18 ± 1.56 (mean ± standard deviation) with minimum and maximum values of 0.10 and 6.90, respectively, for the year 2000; while it was 3.46 ± 2.09 (mean ± standard deviation) with minimum of 0.10 and maximum of 6.60 for the year 2019. Figure 6b shows that the values extracted from the PW pixels on the unchanged mangrove area also showed an increasing trend, increasing from 3.03 to 3.55 for the years 2000 and 2019, respectively. The increases of mean LAI for all PW and unchanged areas indicate that the mangrove leaf area has grown over the years. Previous analysis (Section 3.1 and Figure 3c) show that relatively less mangrove gain was noticed when compared to area of mangrove loss. Based on the values in Figure 6c,d, the mean LAI extracted from the mangrove gain increased from 2.81 to 3.51, while the mangrove loss area recorded mean LAI decrease from 2.99 to 2.62, respectively. The loss of native mangrove trees, which have dense leaves, to oil palm, which have a less dense canopy, may be the main reason for the reduced LAI values. The mean LAI values were found to correspond to the growth or loss of the mangrove in the study area. Note that the standard deviation values of the LAI are relatively high, which might be due to the changes in LAI values but are also dependent on species composition, forest structure, and vegetative phenology, which are all intrinsically related to temperature and rainfall, and this may need to be examined further.

LAI and GPP Analysis
GPP serves as one of the most important parameters in global carbon analysis, ecosystem assessment, and monitoring of environmental changes [56,57]. To date, studies relating to mangrove ecosystem processes (including Net Primary Production) have been few, and the majority of them depend on ground-based flux data (i.e., FLUXNET), which are restricted to certain locations [21,58,59]. Introducing remote sensing products such as MODIS GPP could be a solution that will help in improving and extending the study of mangrove carbon processes at a larger scale. Figure 7 shows the GPP values extracted in this study based on the five cases. For all PW pixels (Figure 7a), a reduction in mean GPP from 6.82 g C m −2 to 5.49 g C m −2 is noticed. This reduction could be due to the lower GPP contributions by the newly gained mangroves in the PW pixels, which may have lowered the productivity of the whole area. Figure 7b shows the GPP values of the area with no change. The mean GPP of these pixels increased from 5.81 g C m −2 to 6.73 g C m −2 . The increase of mean GPP over the no change area shows the great ability of mangrove forests to assimilate CO 2 . Over the mangrove gain areas, the mean GPP shows an increasing trend from 1.88 g C m −2 to 2.78 g C m −2 . Most of the mangrove gained areas were found to be mangroves regrown from abandoned aquaculture ponds. Recent studies [60,61] confirmed that the regrowth of the Rhizophoraceae family on the abandoned ponds is feasible, thus supporting the increase of GPP over the gained areas. For the mangrove loss area, mean GPP increased only slightly from 6.38 g C m −2 to 6.68 g C m −2 . Based on the classification results, most of the mangrove losses are due to transformation of the mangrove areas to oil palm cultivation. The fast-growing oil palm cultivation surrounding the mangrove areas, which are mixed in the larger MODIS pixels (500m), could be the reason for increasing GPP despite the loss of mangrove.
The effect of fragmentation on CO 2 assimilation and biomass accumulation by ecosystems is dependent on the size and location of the fragmented patches [62]. To analyze the impact of mangrove fragmentation on the LAI and GPP, we selected the patch of mangrove area with the largest loss (6883.06 ha in 2000), with a loss 1131.68 ha (16.54%) of mangroves and which had fragmented into 60 smaller patches between 2019 and 2000. The mean LAI and GPP values extracted from the MODIS products (Figures 6e and 7e) show a slight increase from 3.11 to 3.41 and a large change from 4.65 g C m −2 to 6.37 g C m −2 , respectively. Of the fragmented patches, six are of a very large size-total of 5704.50 ha (82.88% of the total area) compared to the MPS (26.63- Table 6), one patch is 17.27 ha, and the remaining 53 patches are smaller than 9 ha. According to Laurance et al. [63], the interior of relatively large fragments tends to change less following fragmentation than it does at the forest margins. This may be due to the edge effects, such as altered microclimate, which do not penetrate into the core habitat of the forest interior. Peh, Lin, Luke, Foster, and Turner [18] concluded that productivity is stable in fragments of more than 9 ha. Therefore, despite serious fragmentation over this selected patch, the increases of LAI and GPP may be supported by the impact of fragmentation on productivity, and the survival rate of trees/plants may be higher in relatively smaller patches because the abundance of solar radiation along the edges of the patches may increase the productivity [63]. In smaller patches (<9 ha), the microclimate (temperature and rainfall) may also be changed, which may affect the productivity.

Mangrove Classification and Loss
SVM classification has been widely explored on different types of remote sensing data and at global, regional, and local scales with high accuracy [64]. Mountrakis et al. [65] reviewed land cover classification using different techniques and concluded that SVM is superior compared to other classification techniques. Classification works related to mangrove area identification and mapping using SVM classifiers showed high accuracy [41,66,67]. SVM proved to be one of the best classifiers where high accuracy (~90%) was achieved based on the Landsat 7 ETM+ pan-sharpened image and the Sentinel-2A image. While several global mangrove maps are available at 30 m [34,53,68,69], the use of higher spatial resolution satellite images (15 m and 10 m) in this study aimed to increase the accuracy of mangrove extent mapping. The classification result showed clear delineation between mangrove, water bodies, and other land use and land cover, providing high quality mangrove cover in 2000 and 2019 for mangrove change identification. In a previous study, Kanniah, Sheikhi, Cracknell, Goh, Tan, Ho, and Rasli [32] performed land use classification of the IM region from 1989 to 2014 using multiple Landsat images. They showed that, during the 21st century, urban area significantly increased, converting from oil palm, rubber, mangrove, and forest.
The spatial pattern of mangrove changes in IM shows a significant loss in the southwestern part, where the three Ramsar sites are located (Figure 3c). The Ramsar sites have been protected since 2003; thus, the loss may be due to a high input of sediments generated via construction and land reclamation activities around the Sungai Pulai estuary [70], where the development of the Tanjung Bin Power Plant, the Forest City, and the expansion area of the Tanjung Pelepas Port (Figure 1) are located. Reports including the Marine Environment Protection Committee (MEPC) [71] stated that the increasing traffic of ships and vessels in this region generates wave currents that can destroy mangroves in coastal areas [72].
Within the mangrove region of IM, 40 sites with aquaculture activities were detected from the satellite images. Aquaculture conversion from mangrove forests is very popular, as can be seen in Indonesia, where there has been a 40% mangrove loss over the past three decades [73]. Figure 4b shows the locations of these aquaculture sites. Most of the aquaculture sites are located on the "others" land use, at the boundary of the mangrove area. Thus, we can conclude that aquaculture activities contributed greatly to the loss of mangrove areas during the study period. At the eastern side of IM (Figure 4b), more aquaculture sites were found compared to the western site (where the Ramsar sites are located). One of the reasons for this could be that the Ramsar sites are protected from these conversion activities while, on the eastern side, where the Johor River flows, the accessibility to water from the Johor River and the richness of nutrients at the site have encouraged the conversion of mangrove areas to aquaculture use.

Fragmentation
Based on all the fragmentation indices, the fragmentation of mangrove in the IM region indicates a serious impact on the ecosystem as the fragmented mangrove forests become ecologically and hydrologically isolated. A reduction of Class Area (CA) and Total Landscape Area (TLA) can decrease the movement zone of animals and discourage flora and fauna succession within or between fragments. An increased Number of Patches (NP) and decreased mean patch size (MPS) mean that increasingly isolated mangrove patches lead to a loss of connectivity between patches, thus prolonging the migration of animals to surrounding habitats (due to the decreased existing habitat size or increased species diversity), which is one of the factors that stimulates species extinction [74]. The fragmentation analysis of the edge and shape indices reveal that human-induced disturbance is significantly influencing the mangrove forests. The decreased MSI value showed a reduction of shape complexity, which could be due to the transformation of mangrove areas (irregular shaped) into regularly shaped/spaced patches of urban, agricultural, and aquaculture uses. The changes of edge and shape complexity could ease human access to the fragmented patches and exploit the remaining mangrove forests. This could accelerate the loss and fragmentation of mangrove areas in the future. Fragmentation analysis can support decision-making on the conservation of mangrove area. The impact of fragmentation to the mangrove ecosystem should be the focus of future study. Identification of the vulnerable mangrove patches and/or areas for immediate actions and further regrowth/preservation plans must be implemented to ensure the longevity of mangrove forests in the IM, as well as in Malaysia as a whole.

Mangrove LAI and GPP
The measurement of leaf area and biomass of mangrove at different scales (local, regional, and global) is very scarce. Traditional destructive methods require a high demand of manpower, cost, and time, and the allometric equations generated for biomass calculation are usually site specific [75]. The majority of the scientific data collection based on flux tower has only built on inland forests, and this kind of data is less available for wetland forests [76]. Therefore, we propose the initiative to relate MODIS LAI and GPP to show their capability of detecting leaf area depletion and GPP rate changes of a selected area with active mangrove loss and fragmentation.
The LAI and GPP values extracted from MODIS products showed reliable values when compared to previous studies. Clough et al. [77] reported the LAI of the main mangrove species (Rhizophora apiculata) in Malaysia to be in the range of 3.10-5.10, which is within the range reported in IM (mean LAI from All PW pixels, 3.18-3.46, respectively in 2000 and 2019). In another study, Shrestha et al. [78] showed MODIS GPP to have a good match with the modelled GPP of a mangrove area in India, which can be useful for carbon assimilation estimation. The mean GPP of mangroves in the IM region is 5.49 g C m −2 (based on 2019 data), and this value is close to the GPP value of the tropical forests in Asia and the Americas, which are in the range of 8.22 and 10.96 g C m −2 . Nevertheless, it should be noted that the sparse resolution of MODIS products (500 m) could have a mixing pixel issue, where the contribution of oil palm plantation can affect the accuracy of the LAI and GPP. Thus, data with higher spatial resolution can improve future analysis.
Mangrove forests serve multiple functions that are significant in terms of ecology and economy, and these functions will be lost or reduced due to mangrove fragmentation. The unique forms and rich biodiversity of mangrove forests can provide tourism and recreational activities, including wildlife-watching and fishing [79]. Mangrove forests can be transformed as habitats and nursery grounds for a variety of flora and fauna [57,80], while mangrove trees are also functional as a source of food, fuel, and building materials [81]. Mangrove forests are also effective at protecting shorelines, thus preventing soil erosion, shoreline depletion, and floods [82]. Recent studies have revealed the potential of mangrove ecosystems in sequestrating CO 2 from the atmosphere, thereby reducing the impact of global warming and climate change [21,83]. The complex structure of mangroves with trunk, branches, leaves, various types of roots, and soil enable the ecosystem to store a high volume of carbon (450 Mg ha −1 ) [83]. Mangrove forests have an equal global averaged biomass (247 t ha −1 ) to humid evergreen tropical forests [57,84]. Significant deforestation and degradation of mangrove forests also leads to the loss of a relatively large carbon pool [9,85].

Conclusions
Mangrove forest loss has become a significant environmental issue over the past few decades, mainly due to the disruption of their ecosystem functions and services. Although various mangrove conservation and preservation actions have been initiated worldwide, mangrove losses and fragmentation are still ongoing in order to fulfil the insatiable needs of human beings. In Iskandar Malaysia, mass development has been ongoing and has destroyed the mangrove cover and fragmented the landscape. Such changes have been shown in this study to affect the areal extent of mangroves, their leaf area, and thereby the potential of the ecosystem to assimilate and accumulate carbon dioxide from the atmosphere. In the last two decades, a net loss of 2907 ha of mangrove has been estimated in this study using freely available medium-high resolution remote sensing data. It is noteworthy that the annual loss rate of 1.12%, is higher than the world's mangrove loss rate, and the study identified urbanization, oil palm plantations, and aquaculture activities as the main reasons for the land conversions in this region. Land cover changes of mangroves involved a reduction in the mean patch size from 105 ha to 27 ha and an increase in the number of mangrove patches from 130 to 402. Other fragmentation indices, such as edge and shape complexity, also showed increased values. The impact of fragmentation on mangroves' primary productivity showed an increase of 37%, which could be related to the smaller patches (< 9 ha) that can allow more sunlight to penetrate the edges of the patches and hence increase the productivity. Thus, it can be concluded that the impact of fragmentation on mangrove productivity is also dependent on the fragmented patch characteristics. Despite most of the mangrove forests in IM being declared as Ramsar Sites, the land use change near these protected forests could aggravate other natural forms of destruction such as soil erosion and tidal changes, which could further affect small patches of mangroves. Nevertheless, the regrowth of mangrove trees has also been seen scattered over the IM region, which could be promoted in the immediate future. Accuracy in the monitoring of the LAI and GPP can assist in precision planting or rehabilitation of species and habitat. It can also be used as a powerful response tool for resource managers. This study introduced a novel approach, combining the mangrove classification results with fragmentation analysis and mangrove LAI and GPP data to reveal the impact of fragmentation on mangrove productivity in a region that hosts important mangrove (Ramsar) sites in Malaysia. Such results are expected to provide standard guidelines for a stronger policy formulation to protect the remaining mangroves in the region and in other mangrove areas. Further studies on how to prevent mangrove loss, how to reconnect the fragmented mangrove forests, and how to improve the carbon stock estimation to prove the value of mangrove to the ecosystem processes are equally important, now and in the future.

Acknowledgments:
We acknowledge the Ministry of Education, Malaysia and Universiti Teknologi Malaysia (through research grants Q.J130000.2452.08G51, Q.J130000.3052.02M01, and Q.J130000.3052.02M11) for providing funding to conduct the study. We thank Nazarin Ezzaty and Fateen Nabilah for their help in downloading and processing the satellite data and the reviewers for their constructive comments and suggestions to improve the quality of this work.

Conflicts of Interest:
The authors declare no conflict of interest.