Last week I looked at areas of mangrove forests around the globe, and where these forests were dying off or being deforested fastest. This week I want to continue with mangroves, and ask: How does a country’s loss of mangrove forests relate to other environmental performance factors?
The Environmental Performance Index (EPI) quantifies how a country is doing with regards to environmental policy and ecosystem health, is a score from 0-100. It is put out by Yale University, who make the data used in their calculation available here. This dataset contains an amazing amount of environmental statistics for each country, all in one place. It does not, however, incorporate data on mangrove loss.
If mangrove loss is caused by poor environmental stewardship, then we would expect that larger mangrove loss would correlate with lower EPI. The following plot looks at this relationship. We’re only plotting the countries that have mangrove forests to lose.
It doesn’t really look like there’s any correlation. Let’s try putting it on a log scale (and to do this we’ll switch to Percent.Loss = -Percent.Change, and replace 0 change with 0.01).
It really does not seem like the Environmental Performance Index correlates with 2000-2012 mangrove forest loss levels. There are countries like Bangladesh, with a low EPI of 42% but barely any mangrove loss. There’s Myanmar, with a low EPI of 49% and a very high 8.4% loss of mangroves over those 12 years. On the other hand, the US has a relatively high EPI of 85%, but lost a significant 2.7% of its mangroves. And there are countries with no recorded loss of mangrove area with EPIs ranging from 87% (Singapore) to 28% (Somalia).
One explanation might be that the causes of mangrove deforestation are more atmospheric, oceanic and environmental, rather than human-related.
Let’s try one last thing: plotting absolute mangrove area loss rather than percent loss.
This seems like gibberish: increasing your mangrove deforestation in 2000-2012 will increase your EPI score in 2016! The p-value for the linear coefficient is 0.0183.
x <- log(-df$Amount.Lost+0.01)
y <- df$EPI.Score.2016
fit <- lm(y~x)
summary(fit)$coefficients
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 63.503241 1.4452413 43.939541 1.063372e-60
## x 1.125302 0.4678195 2.405418 1.829766e-02
Next week I’ll do a more interesting analysis, where I unpack the single EPI score into the 33 environmental factors of “environmental health” and “ecosystem vitality” that go into its calculation, and see which relate to a country’s mangrove loss.