Mbouzi

Author

Alix Biton

Published

March 18, 2026

Chissouia Mbouzi

Présentation

L’îlot Chissouia Mbouzi est une réserve naturelle de 82 hectares de terres avec un sommet a 153m d’altitude donc un relief visible. Il y a 11 hectares de forêt sèche primaire.

Comme on peut le voir visuellement sur ces premières cartes, les dégâts observés apparus entre décembre 2024 et janvier 2025 sont très visibles sur les zones exposées au vent. Ici l’image “Exposition topographique au vent de Mbouzi” recense l’exposition topographique de l’ilôt lorsque la vitesse était maximale en chaque point. On peut dire qu’en chaque point, on a le moment le plus vulnérable.

Premiers tests

terra 1.9.1

Avant de commencer les statistiques, il faut faire des statistiques basiques afin de vérifier quantitativement si les différences sont significatives.

Pour rappel: les valeurs négatives correspondent à une perte de végétation et les valeurs positives correspondent une végétalisation de la zone (possibilité de trouées attention)

# few key stats
statSummary <- data.frame(
  meanDec = mean(v_dec),
  meanJan = mean(v_jan),
  meanChange = mean(v_diff),
  sdChange = sd(v_diff),
  lossPercentage = (sum(v_diff < 0) / length(v_diff)) * 100
)

print(statSummary)

    meanDec   meanJan  meanChange   sdChange lossPercentage
1 0.6253096 0.5700586 -0.05525097 0.09354879        70.6541

t.test(v_dec, v_jan, paired = TRUE)


    Paired t-test

data:  v_dec and v_jan
t = 56.093, df = 9019, p-value < 2.2e-16
alternative hypothesis: true mean difference is not equal to 0
95 percent confidence interval:
 0.05332015 0.05718178
sample estimates:
mean difference 
     0.05525097

hist(v_diff, breaks=50, col="lightblue2", main="Distribution de la différence NDVI",
     xlab="Δ NDVI (Jan - Dec)", ylab="Fréquence")
abline(v=0, col="red", lwd=2, lty=2) # 0 abline

pairs(dtaMbouzi[, c("perte_ndvi", "pente", "topoexp")], 
      main = "Corrélations Topographie / NDVI",
      pch = ".")

cor(dtaMbouzi, use = "complete.obs")

           perte_ndvi      pente    topoexp
perte_ndvi  1.0000000  0.2618887 -0.4985377
pente       0.2618887  1.0000000 -0.2045091
topoexp    -0.4985377 -0.2045091  1.0000000

plot(dtaMbouzi$topoexp,dtaMbouzi$perte_ndvi,
     col = "cornflowerblue")
abline(lm(perte_ndvi ~ topoexp, data = dtaMbouzi), col = "orange", lwd = 2, lty = 2)

D’après ces tests, on peut dire qu’il y a une différence significative entre le changement de NDVI et l’exposition topographique de la zone de Mbouzi.