Newer
Older
---
title: "Malaria Suitability Indicator"
author: "Earth Sciences department, Barcelona Supercomputing Center (BSC)"
date: "`r Sys.Date()`"
output: rmarkdown::html_vignette
vignette: >
%\VignetteEngine{knitr::knitr}
knitr::opts_chunk$set(warning=FALSE, message=FALSE)
```
Climate Suitability for Malaria Transmission
-----------------------------
## Introduction
Climate sensitive infectious diseases are a subject of major concern in the context
of a changing climate. Malaria is an infectious
disease caused by parasites of the genus _Plasmodium_ and transmitted by _Anopheles_
mosquitoes (Boualam et al., 2021). Eradication efforts managed to achieve complete
elimination of endemic circulation of malaria in Europe, although sporadic transmission
events are frequently reported in travelers coming from endemic areas. Studies
evaluating drivers of malaria elimination found that high socio-economic
conditions and increased life expectancy are key determinants in keeping this disease
out of the continent (Fischer et al., 2020). On the
other hand, the risk of malaria resurgence is present as climatic conditions
become more suitable for both the parasite and mosquitoes transmission (Zhao et al., 2016).
The [Lancet Countdown in Europe collaboration]("https://www.lancetcountdown.org/europe/")
was launched in 2021 with the purpose of using a variety of indicators to monitor
trends in the impact of climate on dynamic of different health processes. As
part of this initiative, the Climate Suitability for Malaria Transmission indicator
was adapted from the [Global Lancet Countdown Collaboration]("https://www.lancetcountdown.org/")
(Romanello et al., 2021)
to track the annual number of months suitable of malaria transmission between 1950
and 2021.This indicator is a threshold based model that
overlaps the climatic and environmental requirements of _Anopheles_ mosquitoes and
_Plasmodium vivax_, which was the endemic parasite in the continent
until the mid 1970s (Boualam et al., 2021).
Although the `CSIndicators` was initially developed for providing a set of tools
useful in the context of agricultural practices, it is also possible to expand
its applications to the computation of other type of indicators. In this vignette,
the Climate Suitability for Malaria indicator will be computed using the `TotalTimeExceedingThreshold()`
function. The output will be an multidimensional array containing the number of
In addition to the `CSIndicators` package, the `CSTools` and `s2dv` packages will
be used for data extraction and visualization.
source("https://earth.bsc.es/gitlab/nperez/Flor/-/raw/master/CST_Load_devel_from_s2dv.R")
```
The number of months suitable for malaria transmission between 2010 and 2020 will
be computed for the Iberian Peninsula.
knitr::include_url("https://es.wikipedia.org/wiki/Pen%C3%ADnsula_ib%C3%A9rica#/media/Archivo:Espa%C3%B1a_y_Portugal.jpg")
The Malaria indicator is computed with climate variables extracted from the
[ERA5-Land reanalysis dataset]("https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-land?tab=overview") (Muñoz-Sabater et al., 2021), which span between 1950
and present day, at a 9km resolution. These variables are 2 meter
temperature (`tas`), 2 meter dew point temperature (`tdps`) and total precipitation
(`prlr`). These data are produced
by the Copernicus service and extracted and reprocessed by the **Earth Sciences department**
from the **Barcelona Supercomputing Center**.
The function `CST_Load` from the `CSTools` package allows the user to access the
processed data and delivers the data in the form of an `s2dv_cube` object. These
multidimensional objects store data and metadata in several elements. For further
information about these objects visit [Data retrieval and storage](https://cran.r-project.org/package=CSTools/vignettes/Data_Considerations.html).
`CST_Load` has its requirements for extracting data. First, the path to the files
is introduced using *whitecards* delimited by dollar signs. Second, the temporal
extent and resolution of the output con be provided with `sdates`, which takes a
vector of all dates of interest. The spatial extent is inputted within `lonmax`,
`lonmin`, `latmax` and `latmin.`
```{r}
# Path to file locations using whitecards
path_ERA5_CDS <- list(path='/esarchive/recon/ecmwf/era5land/$STORE_FREQ$_mean/$VAR_NAME$_f1h/$VAR_NAME$_$YEAR$$MONTH$.nc')
# Temporal extent and resolution:
# We will work with all months from 2010 to 2020
year_in <- 2010 # initial year
year_fi <- 2020 # last year
month_in <- 1 # first month
month_fi <- 12 # last month
sdates <- paste0(year_in:year_fi, '0101')
# Extract data
vars <- c('tas', 'tdps', 'prlr')
out <- NULL
for(var in vars){
out[[var]] <- CST_Load_s2dv(var = var,
exp = NULL, # We only require observed data
obs = list(path_ERA5_CDS),
sdates = sdates,
lonmax = 5, lonmin = 350,
latmax = 45, latmin = 35,
storefreq = 'monthly',
leadtimemin = month_in,
leadtimemax = month_fi,
output = "lonlat")
# library(zeallot)
# c(tas, tdps, prlr) %<-% CST_Load_s2dv(var = c('tas', 'tdps', 'prlr'),
# exp = NULL, # We only require observed data
# obs = list(path_ERA5_CDS),
# sdates = sdates,
# lonmax = 355, lonmin = 350,
# latmax = 45, latmin = 43,
# storefreq = 'monthly',
# leadtimemin = month_in,
# leadtimemax = month_fi,
# output = "lonlat")
```
The output object is a list with the three `s2dv_cube` objects with dimensions
And it is possible to have a look at the data:
```{r}
# Sneak peek of a sample of the temperature data in Jan 2010
## 2. Data transformation
### 2.1 Temperature and dew point temperature
First, it is necessary to transform units from Kelvin (K) to Degrees Celsius (C)
```{r}
for(var in c('tas', 'tdps')){
out[[var]]$data <- out[[var]]$data - 273.15
}
summary(out$tas$data)
```
Then, it is possible to calculate relative humidity (hurs) using temperature ($T$)
and dew point temperature ($T_d$) with the August-Roche-Magnus equation (Alduchov 1996):
RH=100*\frac{exp(\frac{aT_d}{b+T_d})}{exp(\frac{aT}{c+T})}
$$
Where _a_ and _b_ are the coefficients 17.625 and 243.04, respectively.
# Create a new s2dv object with relative humidity
out$hurs <- s2dv_cube(data=100*(exp((17.625*out$tdps$data)/(243.04+out$tdps$data)) /
exp((17.625*out$tas$data) / (243.04+out$tas$data))),
lon=out$tas$lon,
lat=out$tas$lat,
Variable=list(varName="rh", level=NULL),
Dates=out$tas$Dates,
Datasets=out$tas$Datasets,
when=Sys.time(),
source_files="see source files of tas, tdps and pr")
attr(out$hurs$Variable, "units") <- "%"
attr(out$hurs$Variable, "longname") <- "near-surface relative humidity"
summary(out$hurs$data)
Precipitation is extracted from the hourly reanalaysis dataset and delivered in
monthly meters per second. The next step is to transform it into accumulated mm
per month
```{r}
# To compute the accumulated precipitation per month it is necessary to multiply by:
# - total number of seconds in one hour: 3600
# - total number of hours in a day: 24
# - average number of days in a month: 365.25/12=30.44
The suitability for malaria transmission indicator estimates the number of months
suitable for transmission of _P. vivax_ according to the following thresholds:
* Accumulated monthly precipitation above 80 mm (Lyon et al., 2017),
* Monthly temperature between 14.5ºC and 33ºC (Grover-Kopec et al., 2006), and
* Monthly relative humidity above 60% (Grover-Kopec et al., 2006)
# Create a new s2dv object with malaria suitability
out$malariaSuit <- s2dv_cube(data=ifelse(out$tas$data >= 14.5 &
out$tas$data <= 33 &
out$hurs$data >= 60 &
out$prlr$data >= 80, 1, 0),
lon=out$tas$lon,
lat=out$tas$lat,
Variable=list(varName="malar_suit", level=NULL),
Dates=out$tas$Dates,
Datasets=out$tas$Datasets,
when=Sys.time(),
source_files="see source files of tas, tdps and pr")
attr(out$malariaSuit$Variable, "units") <- "none"
attr(out$malariaSuit$Variable, "longname") <- "suitability for malaria transmission"
# Compute number of number of months per year that are suitable
# (i.e. suitability = 1)
malariaInd <- CSIndicators::CST_TotalTimeExceedingThreshold(out$malariaSuit,
threshold=0.5)
Using the function `PlotEquiMap()` from the `s2dv` package, it is possible to
have a visual inspection of the number of suitable months of a particular year.
In this, the first and last years of the time series.
# Source master branch for PlotLayout()
source('https://earth.bsc.es/gitlab/es/s2dv/-/raw/master/R/PlotLayout.R')
# Until the dependency with s2dverification in CSTools is updated, it is necessary
# to make these adjustments
source('https://earth.bsc.es/gitlab/es/s2dv/-/raw/master/R/PlotEquiMap.R')
.FilterUserGraphicArgs <- s2dv:::.FilterUserGraphicArgs
.KnownLonNames <- s2dv:::.KnownLonNames
.KnownLatNames <- s2dv:::.KnownLatNames
ColorBar <- s2dv::ColorBar
clim.palette <- s2dv::clim.palette
.IsColor<- s2dv:::.IsColor
.SelectDevice <- s2dv:::.SelectDevice
var=malariaInd$data[1,1,c(1,dim(malariaInd$data)[[3]]),,],
nrow=1, ncol=2,
lon=malariaInd$lon,
lat=malariaInd$lat,
filled.continents=FALSE,
brks=seq(0, 12, by=1),
toptitle="Number of months per year suitable for malaria transmission in the Iberian Peninsula",
titles=c(as.character(year_in), as.character(year_fi)),
title_scale=0.4,
coast_width=2,
filled.oceans=TRUE,
country.borders=TRUE,
intylat=1,
intxlon=1,
fileout = "MI_MapSuitability_2022-8.png",
width = 7, height = 4)
The annual mean number of months suitable for malaria transmission can be extracted
with the `MeanDims()` function from the `s2dv` package. Then, a plot can be created
with functions from the `dplyr` and `ggplot2` packages. Additionally, the `ggplot()`
function allows the user to compute a linear trend in the data.
```{r fig.width=7, fig.height=4}
library(dplyr)
library(ggplot2)
pt <- data.frame(year=seq(year_in, year_fi, by=1),
ind=c(MeanDims(malariaInd$data, c('lon', 'lat'), na.rm=TRUE))) %>%
ggplot(aes(x=year, y=ind)) +
geom_line() +
scale_x_continuous(breaks=seq(year_in, year_fi, by=1)) +
geom_smooth(method='lm', formula= y~x, alpha=0.3, linetype="dashed", size=0.3, se=FALSE) +
labs(x="Year", y="Mean number of months suitable \n for malaria transmission") +
theme_bw()
first <- mean(malariaInd$data[1,1,1,,], na.rm=TRUE)
last <- mean(malariaInd$data[1,1,dim(malariaInd$data)[[3]],,], na.rm=TRUE)
change <- round((last-first)/first * 100, 0)
The number of months suitable for malaria transmission appear to have increased
in the Iberian Peninsula between 2010 and 2020. When comparing the year 2020 with
the year 2010, there was a `r change`% increase in the number of months suitable for
* Alduchov O. Improved Magnus Form Approximation of Saturation Vapor Pressure. _J Appl Meteorl Climatol. 1996; 35(4),601-609._ [URL]("https://journals.ametsoc.org/view/journals/apme/35/4/1520-0450_1996_035_0601_imfaos_2_0_co_2.xml")
* Benali A, Nunes JP, Freitas FB, et al. Satellite-derived estimation of environmental suitability for malaria vector development in Portugal. _Remote Sens Environ. 2014;145:116-130._ DOI:10.1016/J.RSE.2014.01.014
* Boualam MA, Pradines B, Drancourt M, Barbieri R. Malaria in Europe: a historical perspective. _Front Med. 2021; 8:876._ DOI:10.3389/FMED.2021.691095/BIBTEX
* Fischer L, Gültekin N, Kaelin M, Fehr J, Schlagenhauf P. Rising temperature and its impact on receptivity to malaria transmission in Europe: a systematic review. _Travel Med Infect Dis. 2020;36:101815._ DOI:10.1016/J.TMAID.2020.101815
* Grover-Kopec EK, Blumenthal MB, Ceccato P, Dinku T, Omumbo JA, Connor SJ. Web-based climate information resources for malaria control in Africa. _Malar J. 2006;5(1):1-9._ DOI:10.1186/1475-2875-5-38/FIGURES/5
* Lyon B, Dinku T, Raman A, Thomson MC. Temperature suitability for malaria climbing the Ethiopian Highlands. _Environ Res Lett. 2017;12(6):064015._ DOI:10.1088/1748-9326/AA64E6
* Muñoz-Sabater J, Dutra E, Agustí-Panareda A, et al. ERA5-Land: A state-of-the-art global reanalysis dataset for land applications. _Earth Syst Sci Data. 2021;13(9):4349-4383._ DOI:10.5194/ESSD-13-4349-2021
* Romanello M, McGushin A, Di Napoli C, et al. The 2021 report of the Lancet Countdown on health and climate change: code red for a healthy future. _Lancet. 2021;398(10311):1619-1662._ DOI:https://doi.org/10.1016/S0140-6736(21)01787-6
* Zhao X, Smith DL, Tatem AJ. Exploring the spatiotemporal drivers of malaria elimination in Europe. _Malar J. 2016;15(1):1-13._ DOI:10.1186/S12936-016-1175-Z/TABLES/5