Monitoring Mangrove Forest Degradation and Regeneration: Landsat Time Series Analysis of Moisture and Vegetation Indices at Rabigh Lagoon, Red Sea

: Rabigh Lagoon, located on the eastern coast of the Red Sea, is an ecologically rich zone in Saudi Arabia, providing habitat to Avicennia marina mangrove trees. The environmental quality of the lagoon has been decaying since the 1990s mainly from sedimentation, road construction, and camel grazing. However, because of remedial measures, the mangrove communities have shown some degree of restoration. This study aims to monitor mangrove health of Rabigh Lagoon during the time it was under stress from road construction and after the road was demolished. For this purpose, time series of EVI (Enhanced Vegetation Index), MSAVI (Modiﬁed, Soil-Adjusted Vegetation Index), NDVI (Normalized Difference Vegetation Index), and NDMI (Normalized Difference Moisture Index) have been used as a proxy to plant biomass and indicator of forest disturbance and recovery. Long-term trend patterns, through linear, least square regression, were estimated using 30 m annual Landsat surface-reﬂectance-derived indices from 1986 to 2019. The outcome of this study showed (1) a positive trend over most of the study region during the evaluation period; (2) most trend slopes were gradual and weakly positive, implying subtle changes as opposed to abrupt changes; (3) all four indices divided the times series into three phases: degraded mangroves, slow recovery, and regenerated mangroves; (4) MSAVI performed best in capturing various trend patterns related to the greenness of vegetation; and (5) NDMI better identiﬁed forest disturbance and recovery in terms of water stress. Validating observed patterns using only the regression slope proved to be a challenge. Therefore, water quality parameters such as salinity, pH/dissolved oxygen should also be investigated to explain the calculated trends.


Introduction
Mangrove swamps are a common feature in coastal areas that are exposed to the daily fluctuation of tides. They are predominantly found in muddy substrate with low wave energy and hypersaline environment [1][2][3][4][5][6][7]. Mangroves serve as an important habitat for fish communities and other benthic flora and fauna [8]. Moreover, they play a key role in shoreline protection and waste assimilation through the purification of marine water and surrounding air. Mangroves are also important for carbon sequestration since they are a major sink for carbon by storing it in the sediments [9,10]. In addition, mangrove trees have shown great tolerance to low levels of dissolved oxygen in coastal environments [11][12][13][14].
The spatial distribution of mangrove trees along the Red Sea coast of Saudi Arabia was found in 104 locations with an estimated total area of 3452 hectares (ha) [15]. Two major species common to this area are Avicennia marina and Rhizophora mucronata. Lagoons share 15% of the globes' coastal areas. Mangrove swamps of Red Sea are essential for various ecological functions such as nursery for commercial fish species [16], protection of coral reefs [17], and provision of nesting sites for bird species like Goliath heron and Reef heron [18]. There has been an observed degradation of mangrove trees in the Red Sea region. Some of the major drivers of degradation are the clearing of large areas of mangrove for hardwood, shrimp farming, over grazing by camels, construction activities that alter the water level, and tidal flow. For instance, building dams has significantly lowered fresh water flow into the swamps, thereby, leading to a substantial rise in the level of salinity [19]. In addition to that, mangrove species of this region grow in a particularly harsh environment including: salinity levels of more than 40 ppt, sea surface temperature of more than 31 • C in summers, and no permanent source of freshwater recharge [20].
In the context of the Rabigh lagoon, in 1987, a new road was constructed in the northwest of the study area crossing the mouth of the bay. This road obstructed the recurrence of water flow within the swamp. As a result, the amount of freshwater flow decreased and salinity of the bay area gradually elevated, thereby reducing the productivity of the mangrove trees. By the end of 2012, the road was flooded by sea water, completely opening the bay, and consequently improving the health of mangroves due to the enhancement of water movement.
In order to conserve ecologically sensitive regions, it is inevitable to first monitor their status. Remotely sensed data, especially multi-spectral Landsat data, have been widely used for such purposes [21][22][23][24][25]. Essentially the spectral information from satellite imagery is used to correlate it with the biophysical properties of vegetation cover. Several studies have indicated that the NDVI (Normalized Difference Vegetation Index) of mangrove represents high correlation with biomass and leaf area index [26][27][28][29][30]. Other indices such as Enhanced Vegetation Index (EVI) and Soil Adjusted Vegetation Index (SAVI) are also used to monitor plant health as they are most robust than NDVI. For instance, EVI is sensitive to background canopy variations and does not saturate over high biomass. SAVI accounts for background soil component that can interfere with the vegetation signal [31]. Normalized Difference Moisture Index (NDMI) has been used to detect forest disturbance and recovery since it detects variation in the moisture content of the vegetation. It can differentiate between the moisture content of soil and vegetation and is less sensitive to atmospheric scattering; SWIR (shortwave infrared) used for NDMI can penetrate thin clouds as well [32].
Alamahasheer et al. [33] studied changes in Red Sea mangroves using multi-temporal Landsat imagery, Elsebaie et al. [34] used integrated remote sensing and GIS (geographic information system) approach to identify suitable plantation sites for mangroves on the southern coast of the Red Sea, Kumar et al. [35] surveyed distribution of mangroves along the coast of theRed Sea using NDVI, and Monsef and Smith [36] estimated mangrove cover in the Red Sea using Landsat spectral indices. This study aims to monitor the status of mangrove trees in Rabigh Lagoon by comparing spatial and temporal patterns during the time of obstructed water circulation and after the water circulation was resumed. To our knowledge, this is the first study to use hyper temporal data and trend estimation, over this region, in the form of annual images for the past 34 years.

Study Area
Rabigh Lagoon is located on the northeast coast of Red Sea. It, also called as Sharm El Kharar, extends between longitudes 38 • E and 39 • E and latitudes 22 • N and 23 • N. Its origin dates to late Pleistocene because of erosion, followed by flooding due to postglacial sea level rise in early Holocene transgression [37][38][39]. The study region is only 16 km from Rabigh city and occupies an area of about 75 km 2 ( Figure 1) with many patches of mangrove communities spreading over an area of 136.7 ha [35]. Surface elevation of the area differs widely from −10 to 58.5 m above sea level. Sediment texture is mostly mud, gravelly sand, and sandy mud. The prevailing climate is dry and tropical. In winters, there is also influx of freshwater from wadis: Rabigh, Rehab, Murayykh, and Al-Khariq. Tidal range is 20 to 30 cm. Water exchange between the lagoon and the sea takes place mainly via the following process: Water, having temperature of 25 to 30 • C and salinity of 39‰, enters the lagoon as surface inflow. After circulating the lagoon, it exits in the form of subsurface flow with relatively similar temperature and a salinity of 39.8‰ and 40.5‰ during winter and summer, respectively [40]. On average, it takes 15 days for the water to circulate the lagoon, with the speed of flow of 50 cm per second at the entrance 5-20 cm per second inside the lagoon. Avicennia marina is the only mangrove species that is found in the lagoon [41]. They mostly have a stunted growth reaching to a height of 2 m. The Rabigh lagoon has suffered from climatic changes and global warming in the form of increased water temperature and sea level rise [42].
Rabigh city and occupies an area of about 75 km 2 ( Figure 1) with many patches of mangrove communities spreading over an area of 136.7 ha [35]. Surface elevation of the area differs widely from −10 to 58.5 m above sea level. Sediment texture is mostly mud, gravelly sand, and sandy mud. The prevailing climate is dry and tropical. In winters, there is also influx of freshwater from wadis: Rabigh, Rehab, Murayykh, and Al-Khariq. Tidal range is 20 to 30 cm. Water exchange between the lagoon and the sea takes place mainly via the following process: Water, having temperature of 25 to 30 °C and salinity of 39‰, enters the lagoon as surface inflow. After circulating the lagoon, it exits in the form of subsurface flow with relatively similar temperature and a salinity of 39.8‰ and 40.5‰ during winter and summer, respectively [40]. On average, it takes 15 days for the water to circulate the lagoon, with the speed of flow of 50 cm per second at the entrance 5-20 cm per second inside the lagoon. Avicennia marina is the only mangrove species that is found in the lagoon [41]. They mostly have a stunted growth reaching to a height of 2 m. The Rabigh lagoon has suffered from climatic changes and global warming in the form of increased water temperature and sea level rise [42].

Trend Analysis Using Time Series
Time series are data points distributed across equal intervals. Time series of vegetation indices have been used for ecological studies such as monitoring forest and land degradation. This study used 30 m annual images, from 1986 to 2019, generated from sensors onboard Landsat 5 (1986 to 2002), Landsat 7 (2003 to 2012), and Landsat 8 (2013 to 2019) to derive the long-term trend component over mangrove forests in the study region. The different indices used were as follows (Table 1): NDVI, EVI, MSAVI, and NDMI. Images

Trend Analysis Using Time Series
Time series are data points distributed across equal intervals. Time series of vegetation indices have been used for ecological studies such as monitoring forest and land degradation. This study used 30 m annual images, from 1986 to 2019, generated from sensors onboard Landsat 5 (1986 to 2002), Landsat 7 (2003 to 2012), and Landsat 8 (2013 to 2019) to derive the long-term trend component over mangrove forests in the study region. The different indices used were as follows (Table 1): NDVI, EVI, MSAVI, and NDMI. Images for all years belonged to the summer season, which is June to September, depending on data availability and minimum cloud cover. Data were obtained from the online archive of USGS EROS (United States Geological Survey Earth Resources Observation and Science) Center Science Processing Architecture. The final product (normalized indices) downloaded, with path row 170/044 on WRS-2 (World Reference system), was derived using surface reflectance (tier 1) products from Landsat 5 TM (Thematic Mapper), Landsat 7 ETM (Enhanced Thematic Mapper), and Landsat 8 OLI (Operational Land Imager). Surface reflectance products are radiometrically and atmospherically corrected and come along with cloud, shadow, water, and snow masks developed using the CFMASK (C Function of Mask) algorithm. More information about Landsat surface reflectance product characteristics can be found in Schmidt et al. [43] and Vermote et al. [44] findings. All images were rescaled to fall in the range between −1 to 1 by multiplying the values with 0.0001. Images from Landsat 7 between the years 2003 and 2012 with SLC (Scan Line Corrector) failure were gap filled using linear interpolation. All non-vegetated areas and water bodies were masked out from further analysis using a threshold median EVI of 0.08. This threshold was determined by comparing EVI values of vegetation with non-vegetated surface around Rabigh Lagoon. Figure 2 shows spatial distribution of prominent vegetation communities (dark green tones). data availability and minimum cloud cover. Data were obtained from the online archive of USGS EROS (United States Geological Survey Earth Resources Observation and Science) Center Science Processing Architecture. The final product (normalized indices) downloaded, with path row 170/044 on WRS-2 (World Reference system), was derived using surface reflectance (tier 1) products from Landsat 5 TM (Thematic Mapper), Landsat 7 ETM (Enhanced Thematic Mapper), and Landsat 8 OLI (Operational Land Imager). Surface reflectance products are radiometrically and atmospherically corrected and come along with cloud, shadow, water, and snow masks developed using the CFMASK (C Function of Mask) algorithm. More information about Landsat surface reflectance product characteristics can be found in Schmidt et al. [43] and Vermote et al. [44] findings. All images were rescaled to fall in the range between −1 to 1 by multiplying the values with 0.0001. Images from Landsat 7 between the years 2003 and 2012 with SLC (Scan Line Corrector) failure were gap filled using linear interpolation. All non-vegetated areas and water bodies were masked out from further analysis using a threshold median EVI of 0.08. This threshold was determined by comparing EVI values of vegetation with non-vegetated surface around Rabigh Lagoon. Figure 2 shows spatial distribution of prominent vegetation communities (dark green tones).  Workflow in Figure 3 shows how the images were processed and analyzed. Trend estimation was carried out on four (EVI, MSAVI, NDVI, and NDMI) annual raster time series covering the period of the last 34 years. Temporal changes in indices were interpreted as follows: increasing EVI, MSAVI, and NDVI meant the canopy is more productive (green) and in a healthy state. This is due to higher reflection of NIR (near infrared) by the spongy mesophyll of the all the leaves in the canopy and absorption of red wavelength by chloro-phyll. In the case of NDMI, the SWIR1 reflection tells about the vegetation water content and spongy mesophyll structure. Higher water content would mean more absorption of SWIR1, thus leading to higher NDMI [48]. Slope of the yearly trend was calculated by linear least square regression where the index values act as response variable against the predictor variable, which in this case is time. Statistical significance of the regression coefficient was determined through Mann-Kendall trend test. The analysis was done using the "greenbrown" library [49] in RStudio (version 1.2.5033, Affero General Public License v3, Boston, MA, USA). The trend calculation algorithm in the "greenbrown" library offers the possibility to compute breaks in the time series by controlling the parameter h. Statistically, a breakpoint in the time series occurs when there is a structural change in the regression parameters before and after the break [50]. Ordinary Least Square Moving Sum (OLS-MOSUM) test is applied to detect whether any breaks exist in the time series. The detected break in a signal might depict an abrupt change such as clearcutting, logging, forest fire, landslide, and land management practices that might cause the index values to change with a large magnitude. However, breaks can also be induced by data artifacts [51]. For this study, the parameter h, which is the minimum time between detected breaks, was set to 5 years (0.15) for a 34-year time span. This means only the most significant abrupt change, if it exists, in the five-year period will be recorded. The regression analysis resulted in four trend maps showing areas with monotonic upward (greening) and downward (browning) patterns. Figure 3 shows how the images were processed and analyzed. Trend estimation was carried out on four (EVI, MSAVI, NDVI, and NDMI) annual raster time series covering the period of the last 34 years. Temporal changes in indices were interpreted as follows: increasing EVI, MSAVI, and NDVI meant the canopy is more productive (green) and in a healthy state. This is due to higher reflection of NIR (near infrared) by the spongy mesophyll of the all the leaves in the canopy and absorption of red wavelength by chlorophyll. In the case of NDMI, the SWIR1 reflection tells about the vegetation water content and spongy mesophyll structure. Higher water content would mean more absorption of SWIR1, thus leading to higher NDMI [48]. Slope of the yearly trend was calculated by linear least square regression where the index values act as response variable against the predictor variable, which in this case is time. Statistical significance of the regression coefficient was determined through Mann-Kendall trend test. The analysis was done using the "greenbrown" library [49] in RStudio (version 1.2.5033, Affero General Public License v3, Boston, MA, USA). The trend calculation algorithm in the "greenbrown" library offers the possibility to compute breaks in the time series by controlling the parameter h. Statistically, a breakpoint in the time series occurs when there is a structural change in the regression parameters before and after the break [50]. Ordinary Least Square Moving Sum (OLS-MOSUM) test is applied to detect whether any breaks exist in the time series. The detected break in a signal might depict an abrupt change such as clearcutting, logging, forest fire, landslide, and land management practices that might cause the index values to change with a large magnitude. However, breaks can also be induced by data artifacts [51]. For this study, the parameter h, which is the minimum time between detected breaks, was set to 5 years (0.15) for a 34-year time span. This means only the most significant abrupt change, if it exists, in the five-year period will be recorded. The regression analysis resulted in four trend maps showing areas with monotonic upward (greening) and downward (browning) patterns.  All the spatial mapping was carried in ArcGIS Pro (version 2.6.3, Esri, Redlands, CA, USA). Using the trend maps, pixel-wise time series were extracted for each index and were analyzed as follows: to depict the health of mangroves, each time series was divided into two periods; during road obstruction (disturbed period) and after demolition of the road (recovery period). The disturbed period was represented by the years 1987 to 2012, whereas recovery period was represented by years 2013 to 2019. The median was calculated for each period per site and subsequently compared. If the difference between the median index values between the two periods was positive that was interpreted as recovery. Additionally, for easy comparison, the differences between disturbed and recovery periods for all twelve sites were further averaged for each index. Similar calculations were made for the median of raster images. Kruskal-Wallis test was also performed in RStudio to identify any significant differences between the medians of both periods using median values for each site per index. Kruskal-Wallis is a non-parametric approach to detect any significant differences in a continuous dependent variable (median index values) against a categorical independent variable (period of disturbance and post disturbance recovery).

Results
Trend estimations from all four indices were analyzed for their spatial explicitness and temporal development. Figure 4 shows the spatial trend patterns computed by the regression algorithm. Only significant trends with p-values less than and equal to 0.05 (95% significance level) are shown. All trend slopes lie within the range of −0.01 to 0.01. Areas on trend maps ( Figure 4) with bright orange/yellow to green colors indicate the range of weak positive (<0.01) to strong positive (0.01) trend, while areas with dark orange to red indicate the range for weak negative (<−0.01) to strong negative trends (−0.01). According to the histograms in Figure 5a-d, most pixels had a weak positive trend for all indices falling in the range of 0.001 to 0.005. MSAVI showed higher variation in the calculation of trend as can be seen by a greater number of bars. In addition to that, the number of pixels with significant positive trends (more than 0.001) was the highest for MSAVI (Figure 5b). This presents a good case for using MSAVI over regions with sparse vegetation since it can prevent the influence of background soil. The number of pixels with negative trends (slope value more than −0.001) was highest for NDMI, indicating water-stressed mangroves ( Figure 5d). The increasing trend shown in Figure 4 was mostly observed over natural areas with mangrove forests depicting subtle changes. The vegetation around Rabigh lagoon is dominated by mangrove trees [41]. Just near the coast, on the south west of the lagoon, there is a cluster of pixels with a strong positive trend over a region that has undergone modification and land use in terms of managed tree plantation on an area that was previously bare soil. Stronger greening patterns were observed for mangroves inside the lagoon compared to those for the mangroves just at the north western border of the lagoon. Of all four indices, the NDVI was the least reliable. Although overall patterns were similar, the NDVI does not perform well in sparsely vegetated areas, especially where highly heterogeneous pixels are influenced by surrounding bright soil.
Different plots ( Figure 6) of calculated trend slopes for all four indices show correlation between EVI, MSAVI, NDVI, and NDMI and their statistical distribution. Overall correlation values between EVI, MSAVI, and NDVI are quite high, up to 0.98. On the other hand, lower values of correlation between NDMI with other indices were observed, with the maximum reaching up to 0.58. Nonetheless, all indices showed a dominant greening pattern. The box plots and QQ (quantile-quantile) plots show that all indices follow similar statistical distribution. The prominent presence of extreme values (especially positive) is indicated by data points falling on the line in the middle and deviating away from the line (above blue line for positive values and below blue line for negative values) forming heavy tails in the QQ plots. This is also reflected by the data points falling outside the upper and lower whisker in the box plot.       Figure 7 shows the extracted time series over various sites with mangrove trees. A total of 12 time series were extracted. All the graphs show a very subtle development of index values. NDMI shows lower values than the other indices, which reflects water stress for plants largely triggered by obstructed water flow. For all indices, the end of time series can be seen showing higher values, and they stay above the values of those in the beginning of the time series (discussed further in Figure 8). That is why there is an observed monotonic upward trend. Low index values in 2005 were accounted for by the poor quality according to the "pixel quality assurance" band. Therefore, it was considered as noise. Outliers were considered as noise in this study. As long as noise did not result in break detection, cause regression parameters to be unstable, or lead to trend shifts, it was ignored from further investigation. In addition, any rise and fall in the signal between consecutive years was not interpreted as a trend because it was assumed that the mangrove ecosystem takes time to respond to external factors. Moreover, the lag period between the cause and effect has to be kept in mind when detecting vegetation response with satellite instruments. Local scale differences between the time series of the VIs (vegetation index) and moisture index can be attributed to the spectral differences where the VI used the IR (infrared) and visible (red and blue) part of the spectrum for highlighting the greenness of the vegetation. On the other hand, NDMI use only the infrared waves (IR and SWIR) characterizing the wetness of the vegetation.   Figure 8). That is why there is an observed monotonic upward trend. Low index values in 2005 were accounted for by the poor quality according to the "pixel quality assurance" band. Therefore, it was considered as noise. Outliers were considered as noise in this study. As long as noise did not result in break detection, cause regression parameters to be unstable, or lead to trend shifts, it was ignored from further investigation. In addition, any rise and fall in the signal between consecutive years was not interpreted as a trend because it was assumed that the mangrove ecosystem takes time to respond to external factors. Moreover, the lag period between the cause and effect has to be kept in mind when detecting vegetation response with satellite instruments. Local scale differences between the time series of the VIs (vegetation index) and moisture index can be attributed to the spectral differences where the VI used the IR (infrared) and visible (red and blue) part of the spectrum for highlighting the greenness of the vegetation. On the other hand, NDMI use only the infrared waves (IR and SWIR) characterizing the wetness of the vegetation.  Figure 9 shows the spatial pattern of median differences for the whole study site, which also came out to be positive. The average differences (Table 2) between the two periods were similar for EVI and MSAVI, whereas for NDVI, it was the highest. The NDVI usually shows higher values because it tends to saturate over high biomass; therefore, the larger difference was not interpreted as a higher recovery. In the case of NDMI, the difference was only slightly larger than the difference of EVI and MSAVI. This could be a random occurrence. However, it is worthy to mention that NDMI is better at identifying forest disturbance and recovery because SWIR is more influenced by canopy moisture content [52]. For some sites (38.8361426° E 22.9962740° N, 38.9113980° E 22.8789328° N, 38.9140997° E 22.8773751° N), NDMI shows negative values even in the recovery period, which could mean higher water stress. In addition, negative NDMI values also mean low canopy cover, which was illustrated by visualizing high-resolution Google Earth imagery.   In Figure 8, median values for each selected site can be seen changing during the two classified periods of time series. According to the output of Kruskal-Wallis test with p-value of 0.0006918 (against significance level of 0.05), there are significant differences in the median values of disturbance period from the recovery period. For all sites during the disturbance period (1987-2012), median values are lower, showing mangrove communities in an unhealthy state. Post disturbance (2013-2019), after demolishing the road, a gradual recovery of the mangroves was observed, which was reflected in the increased median index values compared to the values in disturbed state. Figure 9 shows the spatial pattern of median differences for the whole study site, which also came out to be positive. The average differences (Table 2) between the two periods were similar for EVI and MSAVI, whereas for NDVI, it was the highest. The NDVI usually shows higher values because it tends to saturate over high biomass; therefore, the larger difference was not interpreted as a higher recovery. In the case of NDMI, the difference was only slightly larger than the difference of EVI and MSAVI. This could be a random occurrence. However, it is worthy to mention that NDMI is better at identifying forest disturbance and recovery because SWIR is more influenced by canopy moisture content [52]. For some sites ( N), NDMI shows negative values even in the recovery period, which could mean higher water stress. In addition, negative NDMI values also mean low canopy cover, which was illustrated by visualizing high-resolution Google Earth imagery.     This study relates degradation and recovery patterns in response to obstructed water recharge, resulting in disturbed levels of salinity, dissolved oxygen, and pH. Tidal and water chemistry are two main parameters considered to effect mangrove productivity [53]. Therefore, it is helpful to consider these factors when interpreting trend patterns. Under higher-salinity and low-nutrient conditions, mangroves increase the process of transpiration to excrete the salt, which could lead to less growth or productivity. Figure 10 shows the aforementioned parameters of the lagoon for the year 2014. The data were taken from the study of Youssef and Sorogy [54] and mapped to visualize the spatial pattern. Since no such data are available during the disturbance period, a direct comparison cannot be made. Nevertheless, values in Figure 10 do suggest that mangroves are returning to a healthy state because the parameters in the surrounding lagoon water are within the normal range (mangroves in Rabigh are acclimatized to 40 ppt salinity) [20]. recharge, resulting in disturbed levels of salinity, dissolved oxygen, and pH. Tidal and water chemistry are two main parameters considered to effect mangrove productivity [53]. Therefore, it is helpful to consider these factors when interpreting trend patterns. Under higher-salinity and low-nutrient conditions, mangroves increase the process of transpiration to excrete the salt, which could lead to less growth or productivity. Figure  10 shows the aforementioned parameters of the lagoon for the year 2014. The data were taken from the study of Youssef and Sorogy [54] and mapped to visualize the spatial pattern. Since no such data are available during the disturbance period, a direct comparison cannot be made. Nevertheless, values in Figure 10 do suggest that mangroves are returning to a healthy state because the parameters in the surrounding lagoon water are within the normal range (mangroves in Rabigh are acclimatized to 40 ppt salinity) [20]. Rainfall and temperature data, provided by CRUTS (Climate Research Unit gridded Time Series), were also used to assess their impact on degradation and regeneration patterns. The data were produced through the interpolation of observations from weather stations. Monthly data was downloaded, and the trend was estimated through linear least square regression on annually aggregated series against a 95% significance level. According to Figure 11, an increasing trend in temperature and non-significant trend in precipitation was observed. Whereas it could be contemplated that a combination of both climatic parameters could have had an influence on mangroves, this cause-and-effect relation must be proceeded on with extreme caution. Firstly, the data had a very coarse spatial resolution. For pixel-to-pixel correlation, climatic and index data must be of the same spatial resolution. Rainfall and temperature were recorded from the monitoring station in Rainfall and temperature data, provided by CRUTS (Climate Research Unit gridded Time Series), were also used to assess their impact on degradation and regeneration patterns. The data were produced through the interpolation of observations from weather stations. Monthly data was downloaded, and the trend was estimated through linear least square regression on annually aggregated series against a 95% significance level. According to Figure 11, an increasing trend in temperature and non-significant trend in precipitation was observed. Whereas it could be contemplated that a combination of both climatic parameters could have had an influence on mangroves, this cause-and-effect relation must be proceeded on with extreme caution. Firstly, the data had a very coarse spatial resolution. For pixel-to-pixel correlation, climatic and index data must be of the same spatial resolution. Rainfall and temperature were recorded from the monitoring station in Jeddah that is at a distance of 141 km from Rabigh. While such datasets might be useful to monitor the impacts of major climatic events and general trend patterns in the context of globally rising temperatures, it might not always be helpful in reflecting a spatially and temporally explicit situation. This is especially the case where a combination (multivariate) of factors have a nonlinear effect and linear statistical correlations might not be enough.
Jeddah that is at a distance of 141 km from Rabigh. While such datasets might be useful to monitor the impacts of major climatic events and general trend patterns in the context of globally rising temperatures, it might not always be helpful in reflecting a spatially and temporally explicit situation. This is especially the case where a combination (multivariate) of factors have a nonlinear effect and linear statistical correlations might not be enough. shown. The symbol * shows the trend is significant, with three * meaning highly signif-icant having p value of less than 0.000001.

Discussion
The deterioration of mangrove species has been observed at a regional and global scale. According to Al Shiekh [55], there has been a loss of 46% of vegetation cover be-tween 1987 and 2002 in the Red Sea coast at the Jazan region, Saudi Arabia. The reasons discussed for this loss are demographic growth and aquaculture. On the contrary, a recent study by Almahasheer et al. [33] reported an expansion of 12% in mangrove area along the Red Sea. This expansion was observed over 41 years from 1972 to 2013 based on multi-temporal Landsat data analysis. Although our study investigated only a small part of the Red Sea, it does correlate with the findings of Almahasheer et al. [33] such that our study showed a monotonic upward trend for mangroves in the form of regeneration after the road in the north west of the lagoon was dismantled. Our study used the approach of "intercomparison" for indirect validation of the observed trends. Intercomparison in-volves [56] using multiple data products to draw out simplified estimations of major spa-tial and temporal patterns. The multiple data products in this study were the four indices using different band combination to assess the temporal change in the spectral response.
Mangroves show a lot of tolerance to various stress factors they are exposed to [57]. The generated time series (Figure 7) of annual records of EVI, MSAVI, NDVI, and NDMI over twelve sites of mangrove habitats in the Rabigh lagoon from 1986 to 2019 showed an interesting pattern of ecological dynamics of mangrove population. There are very subtle changes taking place in the region, where at some sites, mangroves tend to be stable through a loss and gain strategy, and at others, they show a declining tendency. According to Vovides et al. [58], mangrove regeneration post disturbance in arid and semi-arid re-gions could prove to be a challenge. This is especially the case where plantation activities alone, without eliminating the stress factor, will not bring desired results. Twilley et al. [59] investigated similar hydrological stress on mangroves in CGSM lagoon (Cienaga shown. The symbol * shows the trend is significant, with three * meaning highly significant having p value of less than 0.000001.

Discussion
The deterioration of mangrove species has been observed at a regional and global scale. According to Al Shiekh [55], there has been a loss of 46% of vegetation cover between 1987 and 2002 in the Red Sea coast at the Jazan region, Saudi Arabia. The reasons discussed for this loss are demographic growth and aquaculture. On the contrary, a recent study by Almahasheer et al. [33] reported an expansion of 12% in mangrove area along the Red Sea. This expansion was observed over 41 years from 1972 to 2013 based on multi-temporal Landsat data analysis. Although our study investigated only a small part of the Red Sea, it does correlate with the findings of Almahasheer et al. [33] such that our study showed a monotonic upward trend for mangroves in the form of regeneration after the road in the north west of the lagoon was dismantled. Our study used the approach of "intercomparison" for indirect validation of the observed trends. Intercomparison involves [56] using multiple data products to draw out simplified estimations of major spatial and temporal patterns. The multiple data products in this study were the four indices using different band combination to assess the temporal change in the spectral response.
Mangroves show a lot of tolerance to various stress factors they are exposed to [57]. The generated time series (Figure 7) of annual records of EVI, MSAVI, NDVI, and NDMI over twelve sites of mangrove habitats in the Rabigh lagoon from 1986 to 2019 showed an interesting pattern of ecological dynamics of mangrove population. There are very subtle changes taking place in the region, where at some sites, mangroves tend to be stable through a loss and gain strategy, and at others, they show a declining tendency. According to Vovides et al. [58], mangrove regeneration post disturbance in arid and semi-arid regions could prove to be a challenge. This is especially the case where plantation activities alone, without eliminating the stress factor, will not bring desired results. Twilley et al. [59] investigated similar hydrological stress on mangroves in CGSM lagoon (Cienaga Grande de Santa Marta) Colombia, brought on by coastal highway construction. The authors highlighted that the recovery rate of the mangrove wetlands depends on the intensity and persistence of the disturbance, type and age of species, and natural tree mortality, where adult trees after dying leave gaps in the canopy giving space for new seedlings to grow. Based on the study of Twilley et al. [59], natural hydrological restoration and mangrove plantations were the main assessed factors important for mangrove recovery.
In the region of Rabigh, the prominent limiting factors for mangrove growth are optimum hydrological conditions in the form of mainly water circulation from the open sea and limited rainfall. A sand pathway, hindering water circulation, has been present between 1987 and 2012. Since mangrove forests are wetlands and have a deep root system, persistent water obstruction may take a few years to have a prolonged and visible (such as stunted growth and yellowing of leaves) and detectable impact on mangrove's health [60]. This is so because mangrove roots are efficient in absorbing water from sediments. Moreover, since Rabigh Lagoon is surrounded by inland desert, mangrove species there are adapted to dry and warm conditions and will show resilience to drought-like conditions. It should be kept in mind though that at the physiological level, degradation might be immediately detected. The lag discussed in this paper is solely in terms of the satellite detection of vegetation status when viewing the canopy greenness and wetness. The sand passageway was partly flooded by water in 1990, after which the lagoon was connected to the sea through manually laid downpipes (Figure 12). With the passage of time, these pipes were blocked by sedimentation. As a result, there was a slow and reduced water re-charge into the lagoon. During this period, the mangroves were gradually recovering while the degradation processes were still active. Grande de Santa Marta) Colombia, brought on by coastal highway construction. The au thors highlighted that the recovery rate of the mangrove wetlands depends on the inten sity and persistence of the disturbance, type and age of species, and natural tree mortality where adult trees after dying leave gaps in the canopy giving space for new seedlings t grow. Based on the study of Twilley et al. [59], natural hydrological restoration and man grove plantations were the main assessed factors important for mangrove recovery.
In the region of Rabigh, the prominent limiting factors for mangrove growth are op timum hydrological conditions in the form of mainly water circulation from the open se and limited rainfall. A sand pathway, hindering water circulation, has been present be tween 1987 and 2012. Since mangrove forests are wetlands and have a deep root system persistent water obstruction may take a few years to have a prolonged and visible (suc as stunted growth and yellowing of leaves) and detectable impact on mangrove's healt [60]. This is so because mangrove roots are efficient in absorbing water from sediment Moreover, since Rabigh Lagoon is surrounded by inland desert, mangrove species ther are adapted to dry and warm conditions and will show resilience to drought-like cond tions. It should be kept in mind though that at the physiological level, degradation migh be immediately detected. The lag discussed in this paper is solely in terms of the satellit detection of vegetation status when viewing the canopy greenness and wetness. The san passageway was partly flooded by water in 1990, after which the lagoon was connecte to the sea through manually laid downpipes ( Figure 12). With the passage of time, thes pipes were blocked by sedimentation. As a result, there was a slow and reduced water re charge into the lagoon. During this period, the mangroves were gradually recoverin while the degradation processes were still active. In 2012, the sand passageway was completely demolished, and the lagoon was re connected with the open sea. As a result, there was unobstructed natural water circulatio and also an increased nutrient supply. This is reflected by increasing index values for a indices in all twelve sites, especially from 2013 onward where the values are the highes during the whole evaluation period. Three sites in particular were visited as part of a In 2012, the sand passageway was completely demolished, and the lagoon was reconnected with the open sea. As a result, there was unobstructed natural water circulation and also an increased nutrient supply. This is reflected by increasing index values for all indices in all twelve sites, especially from 2013 onward where the values are the highest during the whole evaluation period. Three sites in particular were visited as part of an ongoing project in Rabigh Lagoon since the 1990s and therefore are presented as case studies in the following sections. Figures 13-15 show results from pixel-based trend analysis for EVI, MSAVI, NDVI, and NDMI for the three locations. They are further discussed below. All three sites present a case study of disturbance (due to road construction) and recovery after the road was demolished and the effect of various stress factors was diminished.
under high environmental stress. Slope values for EVI, MSAVI, and NDVI range from 0.0015 to 0.003. The MK (Mann-Kendall) tau correlation coefficient came out to be the highest for EVI, MSAVI, and NDVI, which meant that around 35% variability in VI is explained by time. However, in the case of NDMI, low tau and slope values indicate a weak relation between trend and time and weak trend, respectively. This could indicate that the moisture content of mangroves at this particular location might not be a driving factor for degradation and recovery patterns as much as chlorophyll content is. No breaks were detected, which means there has been a gradual upward trend.

South of Lagoon
This site is present in the second island toward the southern part of the lagoon where freshwater flow is mainly supplied by rain and flood. According to Figure 14, the slope values lie between 0.0015 and 0.003. Only gradual change was observed. MK tau values showed that 20% to 30% variation in index values was explained by time. This site was under sedimentation process, receiving more mud and clay from runoff. The increase of sedimentary conditions may have affected water circulation and increased the deposition of organic matter causing lower index values during the disturbance period.

Extreme Northwest of Lagoon
Located in the Camel Island facing the third entrance of water, it is relatively a small channel. It can be seen from Figure 15 that no abrupt changes were detected by the regression algorithm. At low tide, camels used to cross the shallow water toward the island for grazing. This site has suffered the most degradation, as the index values are lower compared to the other two sites. The calculated trend for EVI, MSAVI, and NDVI was in the range of 0.001. On the contrary, slope value (0.002) and tau value (25%) were the highest for NDMI for this site, suggesting that moisture was the dominant limiting factor. Figure  16 shows photos taken from the field visit conducted at this particular location during various stages of an ongoing project. The temporal variability of mangrove vegetation can be visualized; there is evident degradation in the 1990s followed by a return to a healthy state after natural water flow was resumed.

Extreme Northwest of Lagoon
Located in the Camel Island facing the third entrance of water, it is relatively a small channel. It can be seen from Figure 15 that no abrupt changes were detected by the regression algorithm. At low tide, camels used to cross the shallow water toward the island for grazing. This site has suffered the most degradation, as the index values are lower compared to the other two sites. The calculated trend for EVI, MSAVI, and NDVI was in the range of 0.001. On the contrary, slope value (0.002) and tau value (25%) were the highest for NDMI for this site, suggesting that moisture was the dominant limiting factor. Figure  16 shows photos taken from the field visit conducted at this particular location during various stages of an ongoing project. The temporal variability of mangrove vegetation can be visualized; there is evident degradation in the 1990s followed by a return to a healthy state after natural water flow was resumed.

North of Lagoon
This site, as shown in Figure 13, is in the northern part of the lagoon facing a skewed, deep entrance to the open sea. Due to this, only the tiny deep channel allows water circulation between the lagoon and the open sea. The habitat at this location is mixed between coral reefs and mangroves. Mangroves are scattered in patches. This area has suffered from extreme wind and tides that led to massive erosion; therefore, mangrove trees were under high environmental stress. Slope values for EVI, MSAVI, and NDVI range from 0.0015 to 0.003. The MK (Mann-Kendall) tau correlation coefficient came out to be the highest for EVI, MSAVI, and NDVI, which meant that around 35% variability in VI is explained by time. However, in the case of NDMI, low tau and slope values indicate a weak relation between trend and time and weak trend, respectively. This could indicate that the moisture content of mangroves at this particular location might not be a driving factor for degradation and recovery patterns as much as chlorophyll content is. No breaks were detected, which means there has been a gradual upward trend.

South of Lagoon
This site is present in the second island toward the southern part of the lagoon where freshwater flow is mainly supplied by rain and flood. According to Figure 14, the slope values lie between 0.0015 and 0.003. Only gradual change was observed. MK tau values showed that 20% to 30% variation in index values was explained by time. This site was under sedimentation process, receiving more mud and clay from runoff. The increase of sedimentary conditions may have affected water circulation and increased the deposition of organic matter causing lower index values during the disturbance period.

Extreme Northwest of Lagoon
Located in the Camel Island facing the third entrance of water, it is relatively a small channel. It can be seen from Figure 15 that no abrupt changes were detected by the regression algorithm. At low tide, camels used to cross the shallow water toward the island for grazing. This site has suffered the most degradation, as the index values are lower compared to the other two sites. The calculated trend for EVI, MSAVI, and NDVI was in the range of 0.001. On the contrary, slope value (0.002) and tau value (25%) were the highest for NDMI for this site, suggesting that moisture was the dominant limiting factor. Figure 16 shows photos taken from the field visit conducted at this particular location during various stages of an ongoing project. The temporal variability of mangrove vegetation can be visualized; there is evident degradation in the 1990s followed by a return to a healthy state after natural water flow was resumed. For all three sites, no breaks were found as result of abrupt changes in the index values. This indicates there is no clearcutting and large-scale urbanization that would cause a sudden drop in the index value with at least a change in magnitude exceeding 0.2 units. Mostly the changes are subtle and persistent. That is why despite the curve having a period of low values, the slope line shows up as a long-term pattern of increasing trend. This For all three sites, no breaks were found as result of abrupt changes in the index values. This indicates there is no clearcutting and large-scale urbanization that would cause a sudden drop in the index value with at least a change in magnitude exceeding 0.2 units. Mostly the changes are subtle and persistent. That is why despite the curve having a period of low values, the slope line shows up as a long-term pattern of increasing trend. This study only investigated EVI, MSAVI, NDVI, and NDMI as an indicator for mangrove disturbance and recovery. Other parameters such as water salinity gradients, dissolved oxygen, and other water quality parameters could not be analyzed with the same temporal frequency as the indices. A change in population or community structure occurring at the boundaries of habitats is also an essential aspect that is difficult to monitor with medium-resolution imagery. Additionally noteworthy here is the influence of marine erosion on the stability of mangroves, especially for smaller patches of mangroves, as is the case for Rabigh. It is suggested to consider these additional drivers in future studies regarding mangroves of the Red Sea and particularly in the area of Rabigh lagoon. In addition, to better understand vegetation recovery, experimental setups or projects involving mangrove restoration could be devised to study and prioritize drivers essential for regeneration processes.

Conclusions
This study attempted to monitor mangrove disturbance and the post-disturbance recovery pattern in Rabigh lagoon along the coast of the Red Sea through trend analysis of various indices. Moreover, observed patterns were analyzed in the context of road construction that hindered water circulation through the lagoon. Our aim was to capture long-lived processes as opposed to temporary and short-lived changes. The outcome of this study showed that time series analysis of selected vegetation and moisture indices (EVI, MSAVI, NDVI and NDMI) was able to capture a pattern of disturbance and recovery. Even if the changes were gradual, they could still be identified. It is suggested to use MSAVI and NDMI to assess ecological disturbance and recovery, particularly for sparsely vegetated and arid areas. However, validating observed patterns using only the regression slope proved to be a challenge. This is so because the response of a natural ecosystem to gradual changes is complex. Integrating rain and temperature data in our interpretation scheme was only partially helpful, as it aided in explaining general patterns, but climatic anomalies were not able to account for extreme values in the index time series. Future studies must integrate the estimation of surface water parameters at constant time intervals to relate various stress factors with mangrove response. Furthermore, extensive and regular field visits to visualize the condition of mangrove stands could also help validate the observed patterns.