Understanding changes in mangrove ecosystems driven by human activities, climate change, and environmental variations is essential for effective ecological management. This study focuses on the spatiotemporal variability of the Normalized Difference Vegetation Index (NDVI) and examines its responses to parameters such as sea level (SL), Potential Evapotranspiration (PET), rainfall (RF), Standardized Precipitation Index (SPI-1 month), soil moisture (SM), minimum temperature (TN), and maximum temperature (TX) within the study area. Trends, relative influences, spatial autocorrelation, and relationships between NDVI and climatic-environmental variables, as well as partial correlations, were analyzed using the Mann-Kendall monotonic trend test (MKMT), Relative Weight Analysis (RWA), partial correlation coefficients (PCC), and Multiple Linear Regression (MLR) methods.The spatiotemporal patterns of NDVI reveal a reduction in bare soil and an increase in sparse and dense vegetation from 1987 to 2022. Nevertheless, zones of degradation were observed, particularly in southern Godoria in 2022 compared to 1987, as indicated by NDVI. A notable deterioration in NDVI (> 0.2) was recorded from 2000 to 2012, while the overall interannual trend shows a slight decline.Additionally, analyses with Mann-Kendall and Theil-Sen slope reveal that TN, TX, PET, and SPI-1 show increasing trends, though not statistically significant, while SM and LST show decreasing trends. For environmental variables, SL indicates an upward trend. Further, partial correlation analysis identifies SL, TN, SPI-1, TX, and PET as the primary climatic factors controlling vegetation dynamics during the JJAS season, with PCC values of -0.89, 0.87, 0.77, -0.76, -0.75, and 0.86 with NDVI, respectively.These findings highlight the significant influence of select environmental factors on the spatiotemporal dynamics of mangrove vegetation, providing insights critical for conservation and management efforts.
 
2.1.1.  Partial
correlation
The bivariate correlation coefficients may
not effectively represent the complex relationships among variables in
multivariate correlation analysis, given that multiple factors can influence
these relationships. Therefore, partial correlation coefficients were computed
to assess the spatiotemporal strength and direction of the linear relationship
between NDVI and each climate variable, while controlling for the effects of
the other climate variables (i.e, sea level, PET, SM, SPI, LST, TN and TX). The
strongest correlation is close to 1, while the weakest is below 0.5. Thus, the
partial correlation can be calculated as follow (Cheng et al., 2017):
 
                                                                                                 
(1)
Where x, y, and z represent three distinct
variables. Rxy,z signifies the partial correlation between variables
x and y while accounting for the influence of variable z. Similarly, Rxy
denotes the linear correlation coefficient between x and y, with Rxz
and Ryz conveying analogous interpretations.
To explore the spatial autocorrelation of NDVI data, we
used the "global and local autocorrelation analysis based on Moran’s I
statistics." This method enables the evaluation of the average spatial
differences between individual cells and their adjacent neighbors, thereby
characterizing the spatial attributes of a specific property across the entire
study area through global spatial autocorrelation analysis. In Moran’s
statistics, the normalized z-score can range from -1 to +1. A Moran's I value
exceeding 0 indicates a positive correlation, which suggests a clustering
pattern, while a value below 0 points to a negative correlation, reflecting a
dispersed arrangement. The calculation of Moran’s I statistics for examining
spatial autocorrelation is provided by Xu et al.
(2015):
 (2)
Where N represents the number of observations, xi
denotes the observed value for cell i, xj  indicates the observed
value for cell j, 
 
 
  
  
  
  
  
  
  
  
  
  
  
  
 
 
 
 
 is the average of cell
i or cell j, and wij is the weight assigned to the relationship
between cells i and j.
While the global spatial autocorrelation through
Moran's I statistics reveals the overall clustering pattern, it does not allow for
the assessment of spatial association patterns across multiple locations. In
contrast, Local Spatial Autocorrelation focuses on the significance of local
statistics at each individual location and identifies the presence of spatial
clusters, a capability that global spatial autocorrelation lacks. The mathematical
equation of local spatial autocorrelation using Moran's I is described by (Anselin (2010)).

 (3)
 
  | 
     
  
   
   | 
 
 
  | 
   Fig. 1. (a)Estimation of relative importance of climate variables as predictors of
  NDVI. (b) Pearson correlation coefficient between NDVI and climate variables. 
   | 
 
 
 
Table 1. Mann Kendall Trend and
Theil-sen slope statistics of NDVI and climate variables from 1987 to 2022. 
 
Abdi-Basid ADAN, 2024
The Abdi-Basid Courses Institute (TABCI)