1. Introduction
Mangrove forests play a crucial role in providing valuable ecosystem services, contributing to a staggering annual value of at least US $1.6 billion (Polidoro, Carpenter et al. 2010). These services encompass a wide range of benefits that support coastal livelihoods on a global scale (Dahdouh-Guebas, Jayatissa et al. 2005, Duke, Meynecke et al. 2007, Ellison 2008, Abrantes, Johnston et al. 2015). However, Mangrove forests are disappearing worldwide by 1 to 2% per year (Duke, Meynecke et al. 2007). Clearing for coastal development, expansion of aquaculture, logging for timber, and fuel production are among the primary drivers behind this concerning trend (Daru, Yessoufou et al. 2013, Kirui, Kairo et al. 2013, Yessoufou and Stoffberg 2016). It has been shown that over 40% of the assessed vertebrate species endemic to Mangrove Forests are currently facing global threats to their survival (Luther and Greenberg 2009). As a result, urgent attention is needed to monitor the state of Mangrove forests and get a better understanding of Mangrove Forest dynamics in response to disturbances.
The response of an ecosystem, such as the Mangrove ecosystem, to disturbance is not always a gradual process; it can lead to sudden and irreversible changes (Scheffer 1990, Scheffer 2001, Scheffer, Carpenter et al. 2001). In fact, even gradual changes in the environment may not result in a corresponding gradual response from the ecosystem. Instead, they can trigger sudden, unpredictable, and irreversible shifts known as regime shifts (Capon, Lynch et al. 2015). Apart from the growing evidence of critical changes occurrences in different ecosystems (Barbier, Koch et al. 2008, Guttal and Jayaprakash 2008, Lenton 2011, Verbesselt, Umlauf et al. 2016, Alibakhshi, Groen et al. 2017), the need to enhance the understanding of generating early warning signals of critical transition in ecosystems is highlighted.
Abrupt changes in the state of an ecosystem can occur when ecosystems are unable to cope with the effects of disturbances, leading to slower resilience and a reduced ability to recover (Carpenter and Brock 2006, Carpenter, Brock et al. 2008, Scheffer, Bascompte et al. 2009, Carpenter and Brock 2011, Carpenter, Cole et al. 2011). Disturbances can push the state of an ecosystem to a state that is near a critical threshold. Once a critical threshold is reached, even a small disturbance can trigger a significant transition to a new state, where it is challenging and sometimes impossible to return to the previous state (Scheffer, Carpenter et al. 2001, Carpenter and Brock 2006, Carpenter, Brock et al. 2008, Scheffer, Bascompte et al. 2009, Carpenter and Brock 2011, Carpenter, Cole et al. 2011). This knowledge is crucial for assessing the current state of an ecosystem and determining whether it is approaching a critical transition. Understanding the current state of an ecosystem is particularly crucial in Mangrove Forests, as it can inform conservation efforts and has abundant resources and socio-economic impacts for locals (Martínez, Intralawan et al. 2007). Detecting the state of an ecosystem can prevent irreversible changes and facilitate early interventions (Lenton, Held et al. 2008, Hirota, Holmgren et al. 2011, Alibakhshi 2023).
Various methods have been developed to quantify the state of an ecosystem (Carpenter and Brock 2006, Dakos, Carpenter et al. 2012, Kéfi, Guttal et al. 2014). For example, an increasing trend in autocorrelation and standard deviation of the state variables over time (Dakos, Van Nes et al. 2012) can serve as reliable early warning signals, indicating an approaching critical transition in the ecosystem (Dakos, Carpenter et al. 2012, Dakos, Van Nes et al. 2012, Lenton, Livina et al. 2012, Dakos, Carpenter et al. 2015). Furthermore, disruptions or disturbances within the ecosystem can lead to alterations in the asymmetrical distribution of the state variable time series, resulting in an increased skewness, which can be served as early warning signals of critical transition (Guttal and Jayaprakash 2008, Dakos, Carpenter et al. 2012).
Despite the availability of various methods for assessing early warning signals (Eslami-Andergoli, Dale et al. 2015), the anticipation of critical transitions in ecosystems poses significant challenges. Hence, the application of the methods to assess the state of ecosystems and identify early warning signals of critical transition in the state of ecosystems is still limited. The complexity and diversity of ecological systems, coupled with the need for high-frequency observations and comprehensive time series data on relevant environmental variables, present significant obstacles (Dakos, Carpenter et al. 2015). To address this problem, this study is using the high spatio-temporal resolution of satellite images. Satellite images offer high-resolution, comprehensive coverage, and multispectral capabilities, making them a powerful tool for the aim of this study.
To effectively detect the approach of an ecosystem towards a critical threshold, it is crucial to select the right state variables, which is variable that can accurately present the state of an ecosystem and measure the proximity of the ecosystem to critical conditions (Alibakhshi, Groen et al. 2017). For example, a research study has demonstrated that the utilization of remotely sensed indicators capable of simultaneously capturing variations in both water and vegetation can provide a more comprehensive understanding of the state of an ecosystem compared to relying solely on vegetation-based or water-based indices (Alibakhshi, Groen et al. 2017).
This study aims to investigate the potential of time series analysis of remote sensing data in detecting regime shifts in Mangrove Forest ecosystems and establishing an early warning system. As two main components of Mangrove Forests are vegetation and water, a vegetation-based indicator, namely the Normalized Difference Vegetation Index (NDVI), a water-based indicator, namely the Modified Normalized Difference Water Index (MNDW), and a vegetation-water-based indicator, namely Modified Vegetation Water Ratio (MVWR) has been employed as remotely sensed state variables. Additionally, changes in land-use and land cover-change are quantified using Landsat data from 1996, 2002, 2008, and 2014 to assess the reasons behind the reduced resilience and possible early warning signals of upcoming critical transition in study sites. Mangrove Forest ecosystems in Qeshm Island and Gabrik are chosen as case studies, where Qeshm Island is selected as a representative of an unhealthy ecosystem, while Gabrik serves as the reference site for comparison. The application of critical slowing down in monitoring Mangrove Forests has been rarely tested (Wang, Zhang et al. 2023). This is the first study that explores the potential of remote sensing to explore early warning signals of impending critical transition in Mangrove forests in Iran.
2. Material and Methods
2.1. Study Area
This study focuses on Qeshm Island (Mazraeh and Pazhouhanfar 2018, Kourosh Niya, Huang et al. 2019), located in the southern region of Iran, between the Persian Gulf and the Oman Sea. Qeshm Island, with a total area of 1667 km2, has undergone substantial and, at times, drastic changes in its land-use and land cover-change time (Kourosh Niya, Huang et al. 2019). Gabrik, with a total area of 2496 km2, has been selected as a reference case study, representing a healthy ecosystem. Gabrik shares identical climatic conditions with Qeshm Island and is also characterized by the presence of Mangrove Forest ecosystems (Zahed, Rouhani et al. 2010, Naderloo, Türkay et al. 2013).
The geographical coordinates of Qeshm Island range from 55° 15’ 38” to 56° 16’ 52” E and 26° 32’ 20” to 27° 00’ 00” N, and Gabrik is located southeast of Qeshm Island with geographical coordinates ranging from 55° 15’ 38” to 56° 16’ 52” E and 26° 32’ 20” to 27° 00’ 00” N (
Figure 1). Based on ERA5 data (Muñoz Sabater 2019), the mean temperature in Qeshm Island is 301 Kelvin (28 degrees Celsius), and the mean precipitation is 36 mm from February 18, 2000 to July 31, 2021. In Gabrik, the mean temperature is 301 Kelvin (28 degrees Celsius), with an average precipitation of 28 mm from February 18, 2000 to July 31, 2021. These show the similarity of climatic conditions and local weather patterns of Qeshm Island and Gabrik study sites (
Figure 1). According to the population dataset provided by WorldPop (Sorichetta, Hornby et al. 2015), the estimated population residing in the region of Gabrik is reported to be remarkably low, with a mere 74 individuals. On the other hand, a significantly higher population count of 957 individuals are living in Qeshm Island.
Figure 1.
Monthly temperature in Kelvin (°K) and precipitation millimetres (mm) for Qeshm Island and Gabrik from February 18, 2000, to July 31, 2021.
Figure 1.
Monthly temperature in Kelvin (°K) and precipitation millimetres (mm) for Qeshm Island and Gabrik from February 18, 2000, to July 31, 2021.
Figure 2.
Location of the study sites Qeshm Island and Gabrik (a:c).
Figure 2.
Location of the study sites Qeshm Island and Gabrik (a:c).
2.2. Materials and data collection
2.2.1. State variables
The Normalized Difference Vegetation Index has been obtained from the latest version of the Moderate Resolution Imaging Spectroradiometer (MODIS) product, specifically MOD09A1 (version 006). The 006 version of MOD09A1 has undergone algorithm enhancements aimed at improving its accuracy (Didan, 2015). The data utilized in this study was acquired from the Google Earth Engine platform (Gorelick et al., 2017), which offers a comprehensive repository of geospatial data. The dataset utilized herein possesses a spatial resolution of 500 meters with a temporal resolution of 16 days, enabling comprehensive analysis from February 18, 2000, to July 31, 2021. The MOD09A1 product is derived from atmospherically corrected bi-directional surface reflectance. The Red and Near-Infrared (NIR) bands, specifically within the wavelengths of 0.620 μm to 0.670 μm and 0.841 μm to 0.876 μm, respectively has been used to calculate NDVI (Eq. 1). The NDVI values range from -1 to +1, where a value of -1 indicates bare land, while a value of +1 indicates dense forest.
Furthermore, the Modified Normalized Difference Water Index (MNDWI) is a useful indicator for identifying and quantifying water bodies within an ecosystem (Xu 2006). MNDWI is calculated using MOD09A1 (version 006), which is described earlier. Among the numerous remotely sensed water indices available (Mozumder, Tripathi et al. 2014, Rokni, Ahmad et al. 2014, Li, Chen et al. 2015), MNDWI has been recognized as the most accurate indicator for extracting water area variations (Ji, Zhang et al. 2009, Chen, Huang et al. 2013). MNDWI, similar to NDVI, is a dimensionless index that ranges from -1 to 1. When MNDWI values are below 0, it indicates low water content, which includes soil and vegetation. On the other hand, MNDWI values above 0 indicate high water content and varying water levels. The Green and shortwave infrared (SWIR) bands, specifically within the wavelengths of 1.23 μm to 1.25 μm and 0.55 μm to 0.57 μm, respectively has been used to calculate MNDWI) (Eq. 2).
Additionally, the Modified Vegetation Water Ratio (MVWR) was used to measure early warning signs in aquatic ecosystems (Alibakhshi, Groen et al. 2017) (Eq. 3). The MVWR is sensitive to changes in vegetation water content, which are the main component of Mangrove Forests. It effectively captures variations in water availability and reflects hydrological dynamics, such as seasonal fluctuations and long-term shifts. Moreover, the MVWR proves valuable in assessing vegetation health and stress levels (Tehrani and Janalipour 2021).
2.2.2. Map of mangrove Forest ecosystems and land-use and land cover-change
The land-use and land cover-changes map was calculated with a spatial resolution of 30 meters (Kourosh Niya, Huang et al. 2019), by using Landsat satellite imagery for the years 1996, 2002, 2008, and 2014. The maps provide a classification map of study sites, including six distinct land use classes: agriculture, bare-land, built-up, dense-vegetation, mangrove (Mangrove Forest ecosystem), and waterbody.
In addition, the Global Mangrove Forest Distribution dataset was employed for the year 2000 which provides a map of the world’s Mangrove Forest ecosystems as of the year 2000 with a spatial resolution of approximately 30 meters (Giri, Ochieng et al. 2011, Giri, Ochieng et al. 2013). The data compilation involved analyzing over 1,000 Landsat Thematic Mapper (TM) scenes using a hybrid approach of classification techniques.
2.3. Methods for exploring early warning signals of critical transition
In this study, we only used a large variety of freely available remotely sensed data. All statistical analyses and visualizations were performed in R statistical software (Pinheiro, Bates et al. 2000) QGIS (Qgis 2016), and Google Earth Engine (Gorelick, Hancher et al. 2017).
First Mangrove Forest ecosystems in each study site were delineated using the Global Mangrove Forest Distribution maps (
Section 2.2.2). Second, 100 points in each study site were randomly selected. From each point, the time series of the state variables (
Section 2.2.1) was extracted using a mean function for the period from February 18, 2000, to July 31, 2021. Finally, autocorrelation, skewness, and standard deviation were applied to detect early warning signals of critical transition (Dakos, Carpenter et al. 2012). More specifically, autocorrelation refers to the degree of correlation between the values of the same variables in different observations in the data. The concept of autocorrelation is usually used in time series data to calculate the correlation of value to itself in different observations over a consecutive time. For autocorrelation calculation, the first step is to define a lag operator, which is represented by (
t) and is a time display. The autocorrelation function (ACF) can be calculated using equation (4):
Standard deviation (SD) measures the degree of variability or distribution for a set of data relative to the mean of the same data. SD is obtained from the variance as shown in equation (5):
Skewness is the third statistical index used in this study, which is calculated using equation (6):
Prior to applying metric-based models, the data needs to undergo detrending and smoothing procedures to mitigate the impact of nonstationary conditions on the leading indicators (Dakos, Carpenter et al. 2012). Several detrending approaches are commonly employed, such as Gaussian, Linear, Loess filters, and first-differencing (Lenton, Livina et al. 2012). These methods detrend the data within a rolling window. In this study, a sensitivity analysis was conducted using Kendall’s τ, a nonparametric statistic measuring the association between indicators and time, to identify the optimal size of the rolling window and bandwidth for the Gaussian filter (Bevan and Kendall 1971). Kendall’s τ ranges from -1 to +1, where higher values indicate stronger trends, aiming to identify the detrending settings that best capture trends in the leading indicators. To achieve this, the leading indicators for various rolling window sizes (ranging from 25% to 75% of the time series length) and bandwidths (ranging from 25% to 75% of the time series length for the Gaussian filter) with increments of 10%, using Gaussian, Loess, Linear filters, and first-differencing approaches was calculated.
To ensure that the observed trends in the leading indicators were not due to random chance, 1000 surrogate datasets were generated. These datasets were created by fitting the best linear autoregressive moving average model (ARMA) based on AIC to the residuals obtained after detrending the data. Each surrogate dataset had the same length as the residual time series. Following previous research, the trend estimations from the original data with those from the surrogate data, had similar correlation structures and probability distributions. Kendall’s τ was employed to estimate trends in autocorrelation, skewness, and standard deviation was compared. The probability of finding a trend by chance was measured by comparing Kendall’s τ of the original data with the number of cases in which the statistic was equal to or smaller than the estimates of the simulated records, denoted as P(τ∗≤τ) (Dakos, Carpenter et al. 2012). Additionally, when a clear breakpoint was present in the time series, Kendall’s τ for the period before the breakpoint was reported to ensure robust trend estimations.
3. Results
The Results include an analysis of land-use and land cover-change in Qeshm Island and Gabrik, highlighting changes in agriculture, bare land, built-up areas, mangroves, vegetation, and water bodies (
Section 3.1). Time series of state variables, including NDVI, MNDWI, and MVWR are represented, which show temporal dynamics and seasonal variations of Mangrove Forests in study sites (
Section 3.2). Finally, the early warning signals observed time series of NDVI, MNDWI, and MVWR have been presented (
Section 3.3).
3.1. land-use and land cover-change
Qeshm Island experienced various changes in land-use and land cover-change during the study period (
Table 1). The area dedicated to agriculture witnessed a marginal increase of 338 units (0.65% change). However, there was a significant decrease in bare land, with a change of -57,360 units (-3.65% change). On the other hand, built-up areas expanded dramatically, showing an increase of 35,759 units (64.40% change). Mangroves and vegetation also exhibited positive growth, with changes of 6,049 units (9.07% change) and 6,880 units (70.30% change), respectively. Water bodies slightly expanded, with a change of 8,336 units (0.82% change).
In Gabrik agriculture experienced substantial growth (
Table 1), with a change of 4,708 units (15.22% change). Bare land decreased, albeit to a lesser extent, with a change of -8,629 units (-0.52% change). The built-up areas expanded modestly, with a change of 215 units (4.04% change). Mangroves and vegetation also showed positive growth, with changes of 790 units (9.47% change) and 2,628 units (3.22% change), respectively. Water bodies remained relatively stable, with a minimal change of 288 units (0.04% change).
Figure 3.
Land-use maps of Qeshm Island (A) and Gabrik (B) were extracted from Landsat images from 1996 to 2014.
Figure 3.
Land-use maps of Qeshm Island (A) and Gabrik (B) were extracted from Landsat images from 1996 to 2014.
3.2. Time series of state variables
The extracted data for the NDVI, MNDVI, and MVWR indices represent temporal dynamics and seasonal variations in study sites (Figure 5). The data presented in the table consists of four variables: time, NDVI, MNDWI, and MVWR. The time column represents the date of the measurements. The NDVI values range from 0.01 to 0.36, indicating the density of green vegetation cover in the study area. The MNDWI values range from -0.01 to 0.36, representing the presence of water bodies. The MVWR values range from -0.242 to 0.263, indicating the vegetation’s water content. In Gabrik, the NDVI values range from -0.078 to 0.35, indicating the density of green vegetation cover in the study area. The MNDWI values range from -0.10 to 0.27, representing the presence of water bodies. The MVWR values range from -0.25 to 0.27, indicating the vegetation’s water content.
Qeshm Island has a slightly higher average NDVI value of 0.20, indicating a relatively denser green vegetation cover compared to Gabrik, which has an average NDVI of 0.17. This suggests that Qeshm Island may have a higher overall vegetation density. In terms of water presence, Gabrik exhibits a lower average MNDWI value of 0.04, suggesting a relatively lower presence of water bodies compared to Qeshm Island, which has an average MNDWI of 0.10. This indicates that Qeshm Island may have more abundant water bodies within its study area. Regarding vegetation water content, Gabrik and Qeshm Island have similar average MVWR values of 0.10 and 0.09, respectively.
3.3. Early warning signals of a critical transition
Comparing the values obtained for autocorrelation, Qeshm exhibited significantly higher values (0.411, 0.817, and 0.181) compared to Gabrik (0.037, -0.046, and 0.004), indicating a stronger temporal dependency within the ecosystem. Conversely, Gabrik displayed relatively lower autocorrelation values. Furthermore, the analysis of standard deviation demonstrated notably higher values for Qeshm (0.535, 0.827, and -0.207) compared to Gabrik (-0.689, 0.652, and -0.483). Regarding skewness, Qeshm exhibited positive values (0.631, 0.401), while Gabrik displayed negative values (-0.199, -0.191) for NDVI and MNDWI indices. These contrasting skewness values indicate asymmetry in the distribution of the variables.
4. Discussion
This study illuminates the provision of early warning signals for abrupt changes in the Mangrove Forest ecosystem, offering valuable information about the impact of land-use and land cover-change on ecosystem dynamics. Mangrove Forest ecosystems are crucial for biodiversity conservation, harboring a diverse range of plant and animal species and providing vital habitats and breeding grounds for numerous marine and avian organisms, contributing to the overall ecological resilience of coastal ecosystems (Nagelkerken, Van der Velde et al. 2000, Lugendo, Nagelkerken et al. 2007, Kathiresan 2012). The results showed although the study sites, Qeshm Island and Gabrik share similar climate, geography, and annual rainfall (
Figure 1,
Figure 2), the analysis of remote sensing indices, including NDVI, MNDWI, and MVWR, revealed notable differences between the two study areas (
Figure 4). Comparing Qeshm Island and Gabrik, it is evident that the expansion of built-up areas was more prominent in Qeshm Island than in Gabrik (
Table 1) which can explain the differences in the remotely sensed indices that are representing the variations in vegetation and water (
Figure 4) and observed early warning signals (
Figure 5) in Mangrove Forest ecosystems.
The results revealed that in Qeshm Island, autocorrelation has increased in all three indices compared to Gabrik (
Figure 5), indicating early warning signals of upcoming critical transition and also loss of resilience in coping with the stress of disturbances. In addition, Qeshm Island exhibited higher standard deviation and skewness values, than Gabrik, indicating greater variability and asymmetry in the distribution of index values. These findings suggest that Qeshm Island experiences more pronounced fluctuations which presents a high potential for ecosystem state change occurrences compared to Gabrik (
Figure 5). The observed early warnings could be attributed to various factors such as land-use and land cover-change in the area (
Table 1). The expansion of built-up areas in Qeshm Island, as evidenced by the significant increase in urban development, raises concerns about its potential impacts on the island’s ecological balance (
Figure 3 and
Table 2). Previous studies also have shown that Qeshm Island has undergone a lot of changes, from land-use and land cover-change (Toosi, Soffianian et al. 2019, Tajbakhsh, Karimi et al. 2020), to petroleum pollution (Ebrahimi-Sirizi and Riyahi-Bakhtiyari 2013). In contrast, Gabrik, being the reference site, showed relatively lower values for autocorrelation, standard deviation, and skewness. This indicates a more stable and consistent pattern of index values, implying a healthier and less perturbed ecosystem compared to Qeshm Island (Figs 4 and 5). The lower variability and symmetry in the distribution of data imply a more balanced and stable state in terms of land cover and ecological conditions (Guttal and Jayaprakash 2008).
The results showed that NDVI outperformed MVWR and MNDWI as a robust indicator of ecosystem dynamics for several reasons (
Figure 5). The NDVI, which provides information about vegetation cover (Alatorre, Sánchez-Carrillo et al. 2016, Li, Jia et al. 2019, Alibakhshi 2020, Cabello, Germentil et al. 2021), proved to be highly sensitive in capturing changes in Mangrove Forest ecosystems. The NDVI has been widely used and validated in numerous studies for assessing vegetation dynamics (Verbesselt, Hyndman et al. 2010, Verbesselt, Hyndman et al. 2010, Ruan, Yan et al. 2022, Tran, Reef et al. 2022). It is based on the principle that healthy vegetation, such as Mangrove Forest, reflects more near-infrared (NIR) radiation and absorbs more red light (Eq. 1), resulting in higher NDVI values (Shimu, Aktar et al. 2019). The time series of NDVI has explained 62 % of the global mangrove loss due to land-use and land cover-change (Goldberg, Lagomasino et al. 2020). In addition, in this study, NDVI demonstrated consistent increases in standard deviation, autocorrelation, and skewness in the unhealthy study site of Qeshm, indicating significant variations in vegetation density and dynamics (
Table 2). This behaviour aligns with the expected response of an ecosystem undergoing land cover changes, providing valuable information about the temporal patterns and health status of the vegetation. NDVI provided more robust signals compared with MVWR and MNDWI that exhibited limitations and potential false alarms in assessing ecosystem dynamics. The negative skewness and negative standard deviation observed in the MVWR index in Qeshm suggest unreliable signals and interpretations regarding vegetation water content. Similarly, the high autocorrelation and skewness values observed in the MNDWI index raise concerns about its reliability in detecting water presence. This discrepancy may be attributed to the complex interaction of various factors in Mangrove Forests such as vegetation health, type, and forest structure, which have less effects on MNDWI values and lead to inaccurate generation of warning signals. These findings indicate potential inconsistencies and limitations in using MNDWI and MDNDWI as standalone indicators for assessing ecosystem health. In contrast, NDVI’s ability to capture changes in vegetation density and its consistent patterns between the study areas make it a superior indicator. In sum, the varying performance of state variables in this study can be attributed to the specific characteristics and dynamics of the Mangrove Forest ecosystem.
Despite its valuable findings, this study has certain limitations that should be acknowledged. Firstly, the analysis solely relies on remote sensing data, which may have limitations in accurately capturing certain ecological processes and dynamics at a finer spatial scale. Additionally, the study focused on a specific region (Qeshm Island and Gabrik), limiting the generalizability of the results to other Mangrove Forest ecosystems. However, it should be noted that, due to the current political situation in some countries such as Iran, accessing field data is difficult, and thus the common understanding of ecosystem dynamics is strongly based on satellite data.