- Research
- Open Access

# Modeling the role of environmental variables on the population dynamics of the malaria vector *Anopheles gambiae* sensu stricto

- Paul E Parham
^{1}Email author, - Diane Pople
^{2}, - Céline Christiansen-Jucht
^{2}, - Steve Lindsay
^{3}, - Wes Hinsley
^{2}and - Edwin Michael
^{4}

**11**:271

https://doi.org/10.1186/1475-2875-11-271

© Parham et al.; licensee BioMed Central Ltd. 2012

**Received:**4 April 2012**Accepted:**31 July 2012**Published:**9 August 2012

## Abstract

### Background

The impact of weather and climate on malaria transmission has attracted considerable attention in recent years, yet uncertainties around future disease trends under climate change remain. Mathematical models provide powerful tools for addressing such questions and understanding the implications for interventions and eradication strategies, but these require realistic modeling of the vector population dynamics and its response to environmental variables.

### Methods

Published and unpublished field and experimental data are used to develop new formulations for modeling the relationships between key aspects of vector ecology and environmental variables. These relationships are integrated within a validated deterministic model of *Anopheles gambiae* s.s. population dynamics to provide a valuable tool for understanding vector response to biotic and abiotic variables.

### Results

A novel, parsimonious framework for assessing the effects of rainfall, cloudiness, wind speed, desiccation, temperature, relative humidity and density-dependence on vector abundance is developed, allowing ease of construction, analysis, and integration into malaria transmission models. Model validation shows good agreement with longitudinal vector abundance data from Tanzania, suggesting that recent malaria reductions in certain areas of Africa could be due to changing environmental conditions affecting vector populations.

### Conclusions

Mathematical models provide a powerful, explanatory means of understanding the role of environmental variables on mosquito populations and hence for predicting future malaria transmission under global change. The framework developed provides a valuable advance in this respect, but also highlights key research gaps that need to be resolved if we are to better understand future malaria risk in vulnerable communities.

## Keywords

- Malaria
*Anopheles gambiae*s.s.- Temperature
- Rainfall
- Density-dependence
- Mathematical modeling
- Climate change

## Background

Among the potential effects of climate change on human health, the impact on infectious diseases has attracted increasing attention in recent years[1]. Vector-borne diseases (VBDs) are likely to be particularly vulnerable given the poikilothermic nature of vector survival and development, as well as the effects of temperature on pathogen development. Although the link between climatic variables and transmission has attracted interest for VBDs such as dengue and schistosomiasis, the combined global mortality of these diseases is less than 7% of that due to malaria[2], and this, combined with the significant effects of climatic variables on multiple stages of the transmission cycle, has led to malaria remaining an important focus of ongoing debate regarding climate change and VBDs[3, 4].

In the context of better understanding the role of weather and climate on transmission, two modeling approaches are possible. Statistical models use empirical relationships between climatic variables and past (or current) disease incidence (or prevalence) to predict future disease trends[5, 6]. Mechanistic models, on the other hand, adopt a process-based approach, incorporating known biological, epidemiological and entomological relationships affecting vector and pathogen vital rates and formulating mathematically how these combine[7–9]. Both types of model have important roles to play in improving our understanding of climate-driven transmission changes, but the focus here is on exploiting the explanatory power of the latter.

A vital component in developing reliable VBD transmission models is establishing a realistic model of the vector population dynamics, yet only a few studies have explicitly modeled and parameterized the impact of climatic drivers on vector vital rates[8, 10–12]. While these studies have greatly improved our understanding of the relative importance of temperature, rainfall and relative humidity (RH) on vector populations, they also highlight the need to develop a comprehensive mathematical framework for analysing how a range of environmental factors, arising at different spatial scales, combine at the level of breeding sites to affect stage-specific vector abundance in malaria-affected regions.

This work aims to provide such a framework by formulating and parameterizing environment-vector relationships through surveying and modeling relevant experimental and field data, and incorporating these relationships within a low-dimensional, deterministic mathematical framework. Model simplicity permits ease of integration into malaria transmission models and the model is calibrated and validated against longitudinal *Anopheles gambiae* abundance data from Tanzania[13]. The model also highlights where further experimental and modeling work is required to improve parameterization, in addition to developing a framework readily generalized to different *Anopheles* species and other disease vectors.

## Methods

*An. gambiae*s.s. development and mortality depends on the life cycle stage and that field data available to parameterize mathematical models is often collected daily, a stage-structured, discrete-time model (with a daily time-step) is motivated. An alternative framework is based on physiological, rather than chronological, age and this has been adopted elsewhere[7, 8, 10]. In physiological age-structured models, progression through the life cycle is dependent on temperature conditions within a time-step and the minimum temperature for physiological development. However, while processes such as age-dependent mortality, heterogeneities in larval instars, and oviposition differences between gonotrophic cycles are more naturally incorporated within such approaches, there are several drawbacks of relevance to this article. For a general physiological age-structured model of the form

where$\mathbf{n}={\left({n}_{1}{n}_{2}{n}_{3}{n}_{4}\right)}^{\text{T}}$ and **M** is the projection matrix, the high-dimensional nature of **M** increases by an order of magnitude as temperature measurements become more precise. The dependence of development on other factors (such as RH for adults) also increases the complexity of **M**, as well as making an implicit assumption about the linearity of development with temperature that is often violated. Thus, a low-dimensional approach is instead adopted here, providing a simple, structurally-parsimonious, deterministic model that more transparently illustrates the basic structure that may be built upon in future model development, is considerably easier to construct, analyse and interpret, and may be readily appended to malaria transmission models.

*An. gambiae*s.s. pass through three distinct aquatic stages (eggs, larvae (instars L1 to L4) and pupae) prior to adult development. Let

*n*

_{ i }(

*t*) represent the number of vectors in state

*i*(where

*i*= 1, 2, 3 and 4 refers to eggs, larvae, pupae, and adults respectively). The exposed nature of breeding sites results in considerable vulnerability to environmental influences and the impacts of rainfall, temperature, and biotic effects on immature survival and development are considered here. For immature stages, the daily survival probability

*p*

_{ i }of stage

*i*is assumed to be determined by (independent) factors attributable to the mean daily water temperature

*T*

_{ W }(°C), cumulative daily rainfall

*R*

_{ t }(mm), prolonged periods of desiccation

*D*(days), and density-dependence

*DD*, so that

*i*= 1, 2, 3), while, for adults,${p}_{4}={p}_{4}\left({T}_{A},RH\right)$where

*T*

_{ A }is the mean daily air temperature (°C) and

*RH*the relative humidity (%). If

*n*

_{ i }

*(t)*represents the number of (female)

*An. gambiae*s.s. in stage

*i*at the breeding site at time

*t*, then

where *F*_{
4
} is the average number of eggs laid per day per female adult, *P*_{
i
} is the proportion of vectors surviving and remaining in stage *i* in *t* to *t +* 1, and *G*_{
i
} the proportion surviving and progressing from stage *i* in *t* to *t* + 1. To calculate *P*_{
i
} and *G*_{
i
}, the expressions from[14] are used, namely${P}_{i}=\left(\frac{1-{p}_{i}^{{d}_{i}-1}}{1-{p}_{i}^{{d}_{i}}}\right){p}_{i}$ and${G}_{i}=\frac{{p}_{i}^{{d}_{i}}\left(1-{p}_{i}\right)}{1-{p}_{i}^{{d}_{i}}}$(for all values of *i*), where 0 ≤ *p*_{
i
} ≤ 1 is given by (2) and *d*_{
i
} > 1 the average duration spent in stage *i*. To parameterize the model, the literature is reviewed to source relevant data, as well as using previously unpublished data, to develop, where appropriate, functional forms for *F*_{
4
}, *d*_{
i
} and the components of *p*_{
i
} in (2). The resultant population model (3) is then calibrated and validated against vector abundance data from[15].

## Results and discussion

### Modelling breeding site hydrodynamics

*V*

_{ t }be the volume of water (ml) within the site at time

*t*given a fixed site opening of surface area

*A*

_{ T }(mm

^{2}),

*A’*the exposed surface area of water within the site after rainfall (mm

^{2}) (where

*A’*≤

*A*

_{ T }) (which is then used to calculate the evaporation

*E*

_{ t }from the site at the end of day

*t*), and

*h’*the water depth after all daily rainfall (mm) (see Figure1).

*E*

_{ t }is the evaporation from the site on day

*t*(mm). Since the total volume of water (existing volume plus new rainfall) on day

*t*is$V\text{'}={V}_{t}+\frac{{A}_{T}{R}_{t}}{1000}$ and$A\text{'}=1000\left(\frac{3V\text{'}}{h\text{'}}\right)\text{,}$ substituting into (4) gives

*h’*, consideration of the geometry of the cone before and after rainfall on day

*t*gives, using similar triangles,

*V’*/

*V*

_{ 0 }= (

*h’*/

*h*

_{ 0 })

^{3}(where

*V*

_{ 0 }and

*h*

_{ 0 }are the initial volume and depth of water respectively). Rearranging for

*h’*, using the expression above for

*V’*, and substituting into (5) gives

*E*

_{ t }, the standard FAO Penman-Monteith method is used to first calculate the daily reference crop evapotranspiration

*ET*

_{ 0 }(mm/day)[17] as

Here, Δ is the slope of the vapour pressure curve (kPa°C^{-1}) (which depends on *T*_{
A
}), *R*_{
n
} the daily net radiation transferred to the breeding site (MJm^{-2} day^{-1}) (which, for a given location and day number, depends on the daily cloud fraction *CF* (through its relationship with the number of sunshine hours per day), dew-point temperature *T*_{
DP
} (°C), minimum daily temperature *T*_{
min
} (°C) and maximum daily temperature *T*_{
max
} (°C)), *G* the soil heat flux (MJm^{-2} day^{-1}), *γ* the psychrometric constant (kPa°C^{-1}) (constant for a given site), *U*_{
2
} the wind speed at 2 m (ms^{-1}), *e*_{
s
} the saturation vapour pressure (kPa) (dependent on *T*_{
min
} and *T*_{
max
}), and *e*_{
a
} the actual vapour pressure (kPa) (dependent on *T*_{
DP
}). The climatic variables *R*_{
t
}, *T*_{
A
}, *T*_{
DP
} and *CF* are readily available from the ECMWF ERA-40 re-analysis dataset[18], while *U*_{
2
} may be approximated from *U*_{
10
} (the wind speed at 10 m, available from ERA-40) using the conversion *U*_{
2
} = 0.748*U*_{
10
}[17]. The outgoing heat conduction between the water body and surrounding soil *G* is typically negligible compared to *R*_{
n
}[17] and, as in[17], is neglected here.

*ET*

_{ 0 }, however, due to differences in the reflectivity, heat capacity and typical microclimatic conditions of water bodies compared to crops. Pan evaporation

*E*

_{ pan }, the evaporation rate from pans filled with water and sunken into the ground, is more akin to breeding site conditions and hence

*E*

_{ t }can be estimated as

*K*

_{ p }is an empirically-derived pan coefficient (dimensionless) that depends on the type of pan, breeding site surroundings,

*RH*(obtained from$RH=100exp\left(17.27{T}_{\mathit{DP}}/\left(237.3+{T}_{\mathit{DP}}\right)-17.27{T}_{A}/\left(237.3+{T}_{A}\right)\right)$ and

*U*

_{ 2 }. Although immature

*An*.

*gambiae*s.s. typically prefer clear water, examples of breeding within turbid waters also exist[19], but the turbidity of water does not typically affect

*ET*

_{ 0 }(and hence

*E*

_{ t }) by more than 5%[17], so this is ignored here. Daily values of

*K*

_{ p }are estimated using the empirical tables for Colorado sunken pans (with 1 m radius dry fetch) in[17] based on daily values of

*RH*and

*U*

_{ 2 }. A summary of model parameters is given in Table1.

**Key model variables, parameters, and climatic variables**

Quantity | Definition |
---|---|

| The number of |

| The daily survival probability of stage |

| The average duration spent in stage |

| The volume of the breeding site on day |

| Evaporation from the breeding site on day |

| The number of consecutive days without water in the breeding site (days) |

| Daily mean air temperature (°C) |

| Daily mean water temperature in the breeding site (°C) |

| Total rainfall on day |

| Dew-point temperature (°C) |

| Relative humidity (%) (can be calculated from knowledge of |

| Cloud fraction |

| Minimum daily temperature (°C) |

| Maximum daily temperature (°C) |

| Wind speed at 2 m (ms |

### Environmental influences on immature development

#### Rainfall

*An. gambiae*s.s. prefer to breed in small, shallow, temporary rain pools or stagnant bodies of water fully exposed to the sun (such as hoof marks, tyre tracks or other pools created during land use changes)[21], while other species within the

*An. gambiae*complex differ in their preference for freshwater, brackish and saline water[19]. To capture the dependence of oviposition behaviour on environmental conditions, let

*N*

_{ EP }and

*N*

_{ EO }be the number of eggs per female per oviposition produced and laid (respectively), so that

*t*.

*An. gambiae*s.s. oviposition may be influenced by two signals – a chemical cue directing the suitability of habitat water for oviposition and the existing density of juveniles present[22–25]. Dependence on the latter is quantified using the oviposition index

*OI*introduced in[26]. Using[24] and refitting to find

*OI*as a function of the number of immature per ml

*ρ*

_{ t }(using data on L1 and L2 instars) gives

*n*

_{ i }

*(t)*represents the number of vectors in stage

*i*, the relevant density is${\rho}_{t}={n}_{2}\left(t\right)/{V}_{t}$. Hence, since$OI=\left({N}_{t}-{N}_{s}\right)/\left({N}_{T}+{N}_{s}\right)$, where

*N*

_{ T }and

*N*

_{ S }are the number of eggs laid in the test substrate (pool water with larvae) and control substrate (pool water without larvae) respectively,

*p*

_{ i }(

*R*

_{ t }) represent the daily survival probability of immatures in stage

*i*given rainfall

*R*

_{ t }. It is assumed that${p}_{i}\left({R}_{t}\right)=exp\left(-{\sigma}_{i}{R}_{t}\right)$ (

*i*= 1, 2, 3) where

*σ*

_{ i }quantifies the decrease in survival of stage

*i*. Given the focus on L1 and L4 larvae in[27] and the absence of data elsewhere on egg and pupal mortality due to rainfall, eggs and pupae are assumed to respond similarly to L1 and L4 larvae respectively (although pupal response may differ from L4 larvae in reality due to their ventral air space that aids buoyancy, yet significantly increases mortality if this hydrostatic balance is disrupted[28]). Assuming average L1 and L4 losses of 17.5% and 4.8% per night respectively over the study period with 207 mm rainfall across 26 rainfall nights (K. P. Paaijmans, pers. comm.) gives${\sigma}_{1}={\mathrm{0.0242mm}}^{-1}$ and${\sigma}_{3}={\mathrm{0.00618mm}}^{-1}$. Given that the model here does not distinguish between larval instars, the average duration spent in each instar (as a function of

*T*

_{ W }) is accounted for by interpolating between L1 and L4 mortalities in[27] to determine L2 and L3 survival, whereupon averaging over all temperatures gives${\sigma}_{2}={\mathrm{0.0127mm}}^{-1}$ (Figure2a).

The prolonged absence of water also affects immature longevity; anopheline egg survival in desiccating conditions is two to three weeks[29], while *An. gambiae* s.l. eggs are viable for up to 12 days without water[30]. To model the decrease in egg viability in dry habitats, the findings of[31] are used, which demonstrate that the duration of exposure to desiccating conditions is a better measure of egg viability than soil moisture content. If *p*_{
i
}(*D*) is the daily survival probability of stage *i* given *D* days without water, the functional form *p*_{
i
}(*D*) = 2exp(−*ω*_{
i
}*D*)/(1 + exp(−*ω*_{
i
}*D*)) (*i* = 1, 2, 3) is fitted, where *ω*_{
i
} quantifies the sensitivity of stage *i* to desiccation and the functional form ensures that survival is near unity when *D* is small and approaches zero as desiccation increases. Least-squares estimation using field populations under medium-moisture conditions gives *ω*_{
i
} = 0.405days^{−1} (*R*^{2} > 0.99). Survival of larvae and pupae may be similarly parameterized using[29], which demonstrates that L4 larvae survive significantly better than L1, L2 and L3 instars in such conditions – weighting by the average duration in each instar stage gives *ω*_{2} = 0.855days^{−1} (*R*^{2} = 0.97). In the absence of data on pupal survival, pupae are assumed to demonstrate a similar response to L4 larvae, whereupon using[29] gives *ω*_{2} = 0.602days^{−1} (*R*^{2} = 0.94) (Figure2b).

#### Temperature

*p*

_{ i }(

*T*

_{ W }) and stage duration

*d*

_{ i }(

*T*

_{ W }) for each

*i*. For all three stages, age-independent mortality is assumed and hence${p}_{i}\left({T}_{W}\right)=exp\left(-1/{d}_{i}\left({T}_{W}\right)\right)$ (Figure3a and3b).

*An. gambiae*s.s. eggs survive more than five hours at or above 41°C, with survival decreasing exponentially beyond 40°C. For egg development time

*d*

_{ 1 }(

*T*

_{ W }), the functional form of[33], with the corrected coefficients of Bayoh and Lindsay (unpublished data) (Table2), is adopted.

**Average duration**
d
_{
i
}
**(**
T
_{
W
}
**) of immature stage**
i
**at water temperature**
T
_{
W
}
**(from Bayoh and Lindsay (unpublished data))**

Parameter | Functional form |
---|---|

| $1.011+20.212{\left(1+{\left(\frac{{T}_{W}}{12.096}\right)}^{4.839}\right)}^{-1}$ |

| $8.130+13.794{\left(1+{\left(\frac{{T}_{W}}{12.096}\right)}^{4.839}\right)}^{-1}$− |

| $8.560+20.654{\left(1+{\left(\frac{{T}_{W}}{19.759}\right)}^{6.827}\right)}^{-1}-{d}_{2}\left({T}_{W}\right)$ |

Of the juvenile stages, larval survival demonstrates the strongest dependence on temperature and the effect of competition between *An. gambiae* s.s. and *Anopheles arabiensis* on temperature-dependent survival has been examined[34]. The relationship between survival, development and water temperature, and age-dependent mortality, for *An. gambiae* s.s. is considered in[35]. Larval duration is parameterized as a function of *T*_{
W
} in[5], but this is *An. gambiae* s.l., rather than *An. gambiae* s.s. Moreover, this parameterization is based only on temperatures between 23.0 and 32.8°C and extrapolating to temperature extremes gives inconsistent results with experimental findings in[33] (such as development times around 30 days at 18°C in the former compared to 15 days in the latter). While[10] provides a literature survey of larval development times as a function of *T*_{
W
}, eight of the twelve data points for *An. gambiae* s.s. are calculated from[33] on the assumption of eggs and pupae developing within one day, which is inconsistent with experimental data in the latter. The revised coefficients from Bayoh and Lindsay (unpublished data) are therefore used to determine *d*_{
2
}(*T*_{
W
}).

Aside from the work of[36] on the effects of temperatures from 21.2 to 29.5°C on *An. gambiae* s.l. pupal mortality and[33], there is little experimental data to parameterize pupal development and survival. The latter, with the corrected values in Bayoh and Lindsay (unpublished data), are therefore used to parameterize *d*_{
3
}(*T*_{
W
}).

Finally, it is important to note the importance of using water temperature to calculate juvenile survival and development, rather than air temperature. The difference between mean daily water and air temperatures is typically around 3-6°C depending on factors such as breeding site dimensions, microclimate and weather conditions[32, 37]. To account for this, it is assumed that${T}_{W}={T}_{A}+\Delta T\text{,}$where$\Delta T>0$ is assumed to capture all thermodynamic processes taking place at breeding sites leading to a difference between mean water and air temperatures. Lower and upper temperature thresholds for juveniles are taken from[33].

#### Predation and density-dependence

*Anopheles*development and ecology[38]. Here, only within-stage density-dependent mortality is assumed and the potential effects of juvenile density on adult longevity or fertility are not considered. Figure4 demonstrates the dependence of larval survival on existing larval density (H Tsila, unpublished data), while field populations of

*Anopheles*larvae typically demonstrate low densities (for species that do not breed in tree holes or containers) – some field estimates suggest densities of less than 0.3/ml in rice fields, pools and small ponds[36, 39], while others suggest densities around 0.02-0.06 larvae/ml and 1.5 larvae/ml[40, 41] (respectively). Comparing these estimates with Figure4 suggests that larval densities in field populations occur in regimes where intra-species competition for resources is minimal, suggesting that density-dependent mortality is most likely due to predation, although cannibalism may also occur[42].

*An. gambiae*s.l., with larvae typically distributed negative binomially[43]. To model the effects of density-dependence, these observations are incorporated within application of the framework of[44] for developing first principles population models given knowledge of intra-species competition and spatial distribution. For larval populations following a negative binomial distribution and demonstrating predominantly contest competition (given the dominance of predation, also consistent with findings elsewhere such as[40]),[44] demonstrates that if

*X*

_{ t }is the population size at time

*t*

*m*is the number of resource sites across which the population is distributed,

*λ*the aggregation parameter of the negative binomial distribution and

*b*a positive constant. The value of

*λ*is calculated by averaging the aggregation parameters from the five experiments in[43] for which the negative binomial provides the best fit to obtain$\lambda =1.5$. To determine

*b*, consider, without loss of generality, an arbitrary one litre volume of water within a breeding site and divide this into 1ml blocks (so$m=1000$). The observed difference in juvenile mortality between field data and the contribution from temperature and rainfall is attributed to density-dependence. Given that the datasets used consider survival from L1 instars and that the duration of field studies is often longer than the development time from L1 to pupae, it is assumed that predation reduces larval and pupal survival and acts identically on both stages. Since the population affected is${n}_{2}\left(t\right)+{n}_{3}\left(t\right)$, the number of larvae and pupae per litre is${N}_{23}\left(t\right)=1000\left({n}_{2}\left(t\right)+{n}_{3}\left(t\right)\right)/{V}_{t}$ and${p}_{i}\left(DD\right)={X}_{t+1}/{X}_{t\text{,}}$ the daily larval and pupal survival probability due to density-dependence,

*p*

_{ 2 }(

*DD*) and

*p*

_{ 3 }(

*DD*), is

for *i* = 2 and 3. If *n*_{
2
}(*t*) + *n*_{
3
}(*t*) = 0 or *V*_{
t
} = 0, *p*_{
i
}(*DD*) is assumed to be unity. Predation on eggs is assumed to be negligible by comparison and anopheline rice-field survival data from[39],[41] and[45] is used to provide seven independent datasets to fit *b* at Δ*T = 3°C,* 4°C, 5°C and 6°C. For each dataset, air temperature and rainfall data from the nearest meteorological station (using[46] and where missing values are interpolated) are used to calculate the daily survival and development of larvae and pupae due to climatic influences (assuming fixed vector density and assuming no desiccation effects for rice fields) and estimate the additional mortality required to agree with the study data (attributed to *p*_{
i
}(*DD*)). Two approaches are adopted, namely to (a) calculate the number of juveniles remaining after a fixed number of days (determined by the study design), and (b) track the number of cohort larvae and pupae until less than 0.05% of the original population remain. For method (a), where experimental dates are not specified, *b* is calculated for a range of plausible start dates and the average computed. No significant difference in calculating *b* using these methods is found and$b=0.89$ for$\Delta T=3\mathrm{\xb0C}$ and 0.88 for$\Delta T=4\mathrm{\xb0C}$, 5°C and 6°C.

#### Environmental influences on adult development

*Anopheles*is sensitive to temperature and RH, although few experimental studies have examined this in detail and[11] have recently undertaken a review of parameterization work to date. Although the fitting of[47] has been used in work examining the effects of climatic variables on malaria transmission (such as[8]), this parameterization is inconsistent with[48] demonstrating that

*An. gambiae*s.s. cannot survive longer than one day at 40°C. Thus, the majority of modeling studies to date investigating malaria transmission under changing environmental conditions[5, 7, 9, 49] use

despite its basis on fitting a three-parameter function to three data points in the range 9-40°C[9] (with the 40°C point inconsistent with[48]). This relationship assumes no adverse effects of RH on mortality, which is unlikely given that RH < 50% leads to significantly reduced survival[50]. Field observations of *An. gambiae* adults are only approximately consistent with (18), but reflect the relatively high survival at 22-30°C[36].

*p*

_{ 4 }(

*T*

_{ A }

*,RH*) fitted given mean female survival times at 5-40°C inclusive (in 5°C intervals) and 40-100% RH (at 20% intervals) to obtain

*G*

_{ c }is also temperature-dependent and[9] parameterizes this as

*G*

_{ c }(

*T*

_{ A }) =

*D*

_{ M }/(

*T*

_{ A }

*T*

_{ M }) where

*D*

_{ M }= 36.5°C days and

*T*

_{ M }= 9.9°C. An alternative functional form is given by[11] as

*D*

_{ E }= 37.1°C days and

*T*

_{ E }= 7.7°C. Comparing these formulations, the latter gives longer gonotrophic cycles at temperatures above 17.6°C, the regime generally of interest. On the basis of better agreement with[51], (20) is adopted (Figure5b). To calculate

*F*

_{ 4 },

*d*

_{4}/

*G*

_{ c }(

*T*

_{ A }) gives the average number of ovipositions per adult across her lifetime. If the average number of eggs per oviposition is

*N*

_{ EO }=

*f*

_{ t }

*N*

_{ EP }, the average number of eggs laid over her lifetime is

*d*

_{ 4 }

*N*

_{ EO }/

*G*

_{ c }(

*T*

_{ A }), so that

*F*

_{4}, the average number laid per day, is

Given the absence of age-structure in this model, each gonotrophic cycle is assumed to be of equal duration for all adults and produce the same number of eggs, although studies have shown variation in both[52].

No direct influences of rainfall on adult survival are assumed (with indirect effects through changes in RH captured by (19)) and adult survival is assumed to be density-independent following[53] and the weak, but statistically significant, relationship between adult density and survivorship in[15]. There is some evidence of predation on adult *An. gambiae* s.l. at oviposition sites, with the severity potentially depending on the type of site[39], but there are few quantitative studies in this respect.

#### Model calibration and validation

To assess performance, the model is calibrated and validated against longitudinal *An. gambiae* s.l. abundance data from[13] collected in an environment free of vector controls. Data on *T*_{
A
}, *T*_{
DP
} (for calculation of *RH*), (low) cloud fraction *CF*, and the horizontal and vertical components of 10 m wind speed (to calculate *U*_{
2
}) are taken from the ERA-40 re-analysis dataset[18] for the rural community in Masaika, Tanzania (5 16' 0'' S, 38 49' 60'' E) (with the nearest ERA-40 point at 5^{o}, 0’ 0” S, 37^{o} 30’ 0” E). Rainfall data from the Maji Depot Tanga Rainfall station (at 5 4' 58'' S, 39 5' 21'' E), approximately 35 km from Masaika, is used when available (see[13]), with missing data taken from[18]. Since the daily values of *T*_{
min
} and *T*_{
max
} are not available from[18], we derive empirical relationships between *T*_{
A
} and these variables using data from the nearest meteorological station (Tanga at 5^{o} 4’ 48” S, 39^{o} 4’ 12” E), approximately 34 km from the study site, and apply these relationships (*T*_{
min
} = 0.724*T*_{
A
} + 14.4, with *R*^{
2
} = 0.53, and *T*_{
max
} = 0.728*T*_{
A
} + 28.3, with *R*^{
2
} = 0.61) to ERA-40 data on *T*_{
A
} to estimate the associated values of *T*_{
min
} and *T*_{
max
}.

*An*.

*gambiae*s.l. caught in CDC light traps; further details on mosquito collection and experimental procedures are given in[13]. The model is calibrated using the first twenty months of complete monthly data (August 1998 to March 2000 inclusive) and validated over the subsequent twenty months (April 2000 to November 2001 inclusive). The variation in

*T*

_{ A }and

*RH*over the calibration and validation periods is shown in Figure6. Data on the average number of

*An*.

*gambiae*s.l. caught per trap is available for each weekday in the period (aside from short breaks for public holidays), but not at weekends. Given the daily time-step nature of the model, weekend abundance is estimated using linear interpolation and these values appended to the weekday values. This data is then aggregated by month and the model fitted at this scale. A minimum 365 day burn-in period is applied to remove early model transients.

For model calibration, the average number of adult *An. gambiae* s.l. per light trap is fitted to model output after the burn-in period. To account for the difference in scale between data and the model, the scaled fecundity$\stackrel{\xc2\xaf}{{F}_{4}}={\alpha}_{1}{F}_{4}$ and adult *An*. *gambiae* s.s. abundance$\stackrel{\xc2\xaf}{{n}_{4}}={\alpha}_{2}{n}_{4}$are defined and just three parameters fitted over the calibration period – the scale parameters *α*_{1} and *α*_{2}, and Δ*T*. All other parameters are derived from parameterizations in this paper and local breeding site properties (altitude and latitude). It is assumed that${N}_{\mathit{EP}}=120$ (based on model calibration in[7]) and breeding site dimensions consistent with the characteristics of typical *An*. *gambiae* s.s. habitats (in the presence of multiple *An*. *gambiae* s.l. species given the collection of multiple *Anopheles* species in data collection in[13]) reported in[21], namely *A*_{
T
} = 1.79 x 10^{6}mm^{2} and *h*_{
0
} = 97mm. An initial water volume of 1 litre is assumed (*V*_{
0
} = 1000ml). Model fit to data is found to be independent of the initial conditions, so 100 mosquitoes are arbitrarily initially assumed to be in each lifecycle stage.

*R*

^{ 2 }= 0.84). Running the model for a further 20 months with these parameters and assessing the goodness-of-fit gives

*R*

^{ 2 }= 0.50 across the validation period (Figure7). The model is encouragingly able to capture the overall decline in

*An. gambiae*s.l. abundance in Masaika reported in[13] across the calibration and validation periods, as well as the general seasonal trend (although the timing of the two abundance peaks in the validation period are underestimated by one month in both cases, as well as the magnitude of the peaks). The water volume within the breeding site over time (with these best-fit parameters) is shown in Figure8, while the immature population dynamics, and estimated daily water temperature, are plotted in Figure9. Alternatively fitting the model across the entire 40 months of data (Figure10) gives

*R*

^{ 2 }= 0.70 (with${\alpha}_{1}=280.486$,${\alpha}_{2}=0.026$ and$\Delta T=6.1\mathrm{\xb0C}$) and, in this case, the timing of two of the three seasonal peaks are correctly predicted, as well as the approximate severity of these peaks. The fitted values of$\Delta T=6.9\mathrm{\xb0C}$ and 6.1°C are slightly greater than typical Δ

*T*values observed in the field (values in[37], for example, lie in the range 4.0-6.1°C on clear days for three different sized pools), and this reflects the simplified nature of the${T}_{W}={T}_{A}+\Delta T$ formulation and fitting a single mean value of Δ

*T*across annual timescales. Future refinements will improve this component of the model by calculating Δ

*T*from thermodynamical principles based on daily weather conditions and this is expected to further improve model fit. Nonetheless, the model offers the potential for mechanistic insight into vector response to temperature, rainfall, RH, wind speed and cloudiness, and hence how future changes in these variables may affect mosquito dynamics. The results suggest that the observed decline in vector numbers (and malaria) reported in[13] could, in turn, be due to long-term changes in environmental conditions. Further model analysis (such as application of the methods of matrix population modeling[54]) will provide valuable insight into the dominant environmental variables influencing the observed changes in vector numbers, as well as furthering our understanding of the dominant drivers on short and long-term timescales. While such analysis is beyond the scope of this paper and will follow in a forthcoming article, these results highlight the explanatory power of validated mathematical models and their role in evaluating the effects of temporal changes in weather and climate on vector dynamics and, ultimately, disease transmission.

## Conclusions

Along with *An. arabiensis* and *Anopheles funestus*, *An. gambiae* s.s. is one of the principal malaria vectors in Africa[19] and understanding its ecology and dynamics is vital in better understanding the associated impact on malaria transmission and the prospects for eradication[55], as well as the effectiveness of vector controls in different communities and settings. Vector population dynamics are driven by a range of biotic and abiotic factors and clarifying the role of both is key, particularly in the context of how climate change may influence the future spread and distribution of VBDs. Here, a useful framework for understanding how changes in rainfall, temperature, RH, wind speed and cloudiness (both mean values and temporal variability), and density-dependence, at breeding sites may influence vector abundance is presented. By calibrating and validating the model against longitudinal abundance data, this framework is shown to be capable of reproducing the observations in[13] on long-term timescales, suggesting a mechanistic underpinning of mosquito dynamics in terms of environmental variables, an important result given the ongoing debate regarding the link between malaria transmission and climatic changes in Africa[3, 4]. This work also highlights the power of mathematical models in addressing key questions surrounding the role of environmental variables, compared to the multitude of other ecological, epidemiological, socioeconomic and demographic factors, on disease transmission[1]. An important advance of this work is the construction of a modeling framework enabling the linkage of climatic events at large spatial scales to processes at the localized scale of vector breeding sites, enabling assessments of how climatic phenomena at different scales may affect disease transmission in host communities.

Model reliability may be enhanced with improved parameterization and future experimental and modeling research will lead to further understanding of species-specific *Anopheles* population dynamics and their response to environmental variables. These include (i) improving our understanding of *Anopheles* oviposition behaviour, (ii) better quantifying the role of rainfall and temperature on egg, larval and pupal survival, as well as the role of heterogeneities, such as body size, that might influence response, (iii) improved modeling of the relationship between air and water temperatures at breeding sites, (iv) improving our understanding of density-dependent effects on juvenile and adult development and survival (including intra-specific competition, inter-specific interactions between species, cannibalistic tendencies, and predation, as well as their dependence on climatic variables), (v) assessing evidence for age-dependent mortality in juveniles and adults, and (vi) better understanding variability in gonotrophic cycles.

New longitudinal vector studies that simultaneously measure changes in environmental variables are also required to improve the validity and reliability of vector models, which will not only further our understanding of dominant factors driving mosquito dynamics, but will also improve our understanding of the implications for VBD transmission. Nonetheless, the approach here not only provides a useful framework for *An. gambiae* s.s. modeling, but its structure may be readily applied to other *Anopheles* species with suitable parameterization, as well as other vectors (such as *Aedes* or *Culex*). This will ultimately enable a better understanding of the response of a variety of VBDs to environmental change, an important question given the likely influences of weather and climate on many regions of VBD risk over the coming decades.

## Declarations

### Acknowledgements

The authors would like to thank Henri Tsila for providing unpublished data that formed the basis of Figure4. PEP and EM would like to thank the Grantham Institute of Climate Change at Imperial College London for funding this research. EM also acknowledges the Eck Institute for Global Health at the University of Notre Dame for part funding this work. None of the funding bodies mentioned contributed to the design, collection, analysis or interpretation of data, nor the writing of the manuscript or decision to submit. The authors would like to thank the two anonymous reviewers whose comments greatly improved this manuscript.

## Authors’ Affiliations

## References

- McMichael A, Woodruff R, Hales S: Climate change and human health: present and future risks. Lancet. 2006, 367: 859-10.1016/S0140-6736(06)68079-3.View ArticlePubMedGoogle Scholar
- WHO: The global burden of disease: 2004 update. 2008, World Health Organization, GenevaGoogle Scholar
- Hay SI, Cox J, Rogers DJ, Randolph SE, Stern DI, Shanks GD, Myers MF, Snow RW: Climate change and the resurgence of malaria in the East African highlands. Nature. 2002, 415: 905-909. 10.1038/415905a.PubMed CentralView ArticlePubMedGoogle Scholar
- Pascual M, Ahumada JA, Chaves LF, Rodo X, Bouma M: Malaria resurgence in the East African highlands: temperature trends revisited. Proc Natl Acad Sci USA. 2006, 103: 5829-5834. 10.1073/pnas.0508929103.PubMed CentralView ArticlePubMedGoogle Scholar
- Craig MH, Snow RW, le Sueur D: A climate-based distribution model of malaria transmission in sub-Saharan Africa. Parasitol Today. 1999, 15: 105-111. 10.1016/S0169-4758(99)01396-4.View ArticlePubMedGoogle Scholar
- Rogers DJ, Randolph SE: The global spread of malaria in a future, warmer world. Science. 2000, 289: 1763-1766.View ArticlePubMedGoogle Scholar
- Ermert V, Fink A, Jones A, Morse A: Development of a new version of the Liverpool Malaria Model. II. Calibration and validation for West Africa. Malar J. 2011, 10: 62-10.1186/1475-2875-10-62.PubMed CentralView ArticlePubMedGoogle Scholar
- Hoshen MB, Morse AP: A weather-driven model of malaria transmission. Malar J. 2004, 3: 32-10.1186/1475-2875-3-32.PubMed CentralView ArticlePubMedGoogle Scholar
- Martens P: Health and climate change: modeling the impacts of global warming and ozone depletion. 1998, Earthscan, Oxford, UK, 176-Google Scholar
- Depinay J-M, Mbogo C, Killeen G, Knols B, Beier J, Carlson J, Dushoff J, Billingsley P, Mwambi H, Githure J, Toure A, McKenzie F: A simulation model of African Anopheles ecology and population dynamics for the analysis of malaria transmission. Malar J. 2004, 3: 29-10.1186/1475-2875-3-29.PubMed CentralView ArticlePubMedGoogle Scholar
- Ermert V, Fink A, Jones A, Morse A: Development of a new version of the Liverpool Malaria Model. I. Refining the parameter settings and mathematical formulation of basic processes based on a literature review. Malar J. 2011, 10: 35-10.1186/1475-2875-10-35.PubMed CentralView ArticlePubMedGoogle Scholar
- White MT, Griffin JT, Churcher TS, Ferguson NM, Basanez MG, Ghani AC: Modelling the impact of vector control interventions on
*Anopheles gambiae*population dynamics. Parasit Vectors. 2011, 4: 153-10.1186/1756-3305-4-153.PubMed CentralView ArticlePubMedGoogle Scholar - Meyrowitsch DW, Pedersen EM, Alifrangis M, Scheike TH, Malecela MN, Magesa SM, Derua YA, Rwegoshora RT, Michael E, Simonsen PE: Is the current decline in malaria burden in sub-Saharan Africa due to a decrease in vector population?. Malar J. 2011, 10: 188-10.1186/1475-2875-10-188.PubMed CentralView ArticlePubMedGoogle Scholar
- Crouse DT, Crowder LB, Caswell H: A stage-based population model for loggerhead sea turtles and implications for conservation. Ecology. 1987, 68: 1412-1423. 10.2307/1939225.View ArticleGoogle Scholar
- Diuk-Wasser MA, Toure MB, Dolo G, Bagayoko M, Sogoba N, Traore SF, Manoukis N, Taylor CE: Vector abundance and malaria transmission in rice-growing villages in Mali. Am J Trop Med Hyg. 2005, 72: 725-731.PubMed CentralPubMedGoogle Scholar
- Koenraadt CJM, Githeko AK, Takken W: The effects of rainfall and evapotranspiration on the temporal dynamics of Anopheles gambiae s.s. and Anopheles arabiensis in a Kenyan village. Acta Trop. 2004, 90: 141-153. 10.1016/j.actatropica.2003.11.007.View ArticlePubMedGoogle Scholar
- Allen RG, Pereira LS, Raes D, Smith M: Crop evapotranspiration – Guidelines for computing crop water requirements. FAO Irrigation and drainage paper 56. Rome. 1998, FAO – Food and Agriculture Organization of the United Nations, Google Scholar
- ECMWF: data-portal.ecmwf.int.Google Scholar
- Sinka ME, Bangs MJ, Manguin S, Coetzee M, Mbogo CM, Hemingway J, Patil AP, Temperley WH, Gething PW, Kabaria CW, Okara RM, Van Boeckel T, Godfray HC, Harbach RE, Hay SI: The dominant Anopheles vectors of human malaria in Africa, Europe and the Middle East: occurrence data, distribution maps and bionomic precis. Parasit Vectors. 2010, 3: 117-10.1186/1756-3305-3-117.PubMed CentralView ArticlePubMedGoogle Scholar
- Bi P, Tong S, Donald K, Parton K, Ni J: Climatic variables and transmission of malaria: a 12-year data analysis in Shuchen County. China. Public Health Reports. 2003, 118: 65-10.1016/S0033-3549(04)50218-2.View ArticlePubMedGoogle Scholar
- Gimnig JE, Ombok M, Kamau L, Hawley WA: Characteristics of larval anopheline (Diptera: Culicidae) habitats in Western Kenya. J Med Entomol. 2001, 38: 282-288. 10.1603/0022-2585-38.2.282.View ArticlePubMedGoogle Scholar
- McCrae AW: Oviposition by African malaria vector mosquitoes. II. Effects of site tone, water type and conspecific immatures on target selection by freshwater
*Anopheles gambiae*Giles, sensu lato. Ann Trop Med Parasitol. 1984, 78: 307-318.PubMedGoogle Scholar - Munga S, Minakawa N, Zhou G, Barrack OO, Githeko AK, Yan G: Effects of larval competitors and predators on oviposition site selection of
*Anopheles gambiae*sensu stricto. J Med Entomol. 2006, 43: 221-224. 10.1603/0022-2585(2006)043[0221:EOLCAP]2.0.CO;2.View ArticlePubMedGoogle Scholar - Sumba LA, Ogbunugafor CB, Deng AL, Hassanali A: Regulation of oviposition in Anopheles gambiae s.s.: role of inter- and intra-specific signals. J Chem Ecol. 2008, 34: 1430-1436. 10.1007/s10886-008-9549-5.View ArticlePubMedGoogle Scholar
- Sumba LA, Okoth K, Deng AL, Githure J, Knols BG, Beier JC, Hassanali A: Daily oviposition patterns of the African malaria mosquito
*Anopheles gambiae*Giles (Diptera: Culicidae) on different types of aqueous substrates. J Circadian Rhythms. 2004, 2: 6-10.1186/1740-3391-2-6.PubMed CentralView ArticlePubMedGoogle Scholar - Kramer WL, Mulla MS: Oviposition attractants and repellents of mosquitoes: oviposition responses of Culex mosquitoes to organic infusions. Environ Entomol. 1979, 8: 1111-1117.View ArticleGoogle Scholar
- Paaijmans KP, Wandago MO, Githeko AK, Takken W: Unexpected high losses of
*Anopheles gambiae*larvae due to rainfall. PLoS One. 2007, 2: e1146-10.1371/journal.pone.0001146.PubMed CentralView ArticlePubMedGoogle Scholar - Romoser W, Lerdthusnee K, Berry R, Kittayapong P: Effect of mechanical shock on hydrostatic balance and survival of mosquito pupae. J Am Mosq Control Assoc. 1994, 10: 10-PubMedGoogle Scholar
- Koenraadt CJ, Paaijmans KP, Githeko AK, Knols BG, Takken W: Egg hatching, larval movement and larval survival of the malaria vector Anopheles gambiae in desiccating habitats. Malar J. 2003, 2: 20-10.1186/1475-2875-2-20.PubMed CentralView ArticlePubMedGoogle Scholar
- Beier JC, Copeland R, Oyaro C, Masinya A, Odago WO, Oduor S, Koech DK, Roberts CR:
*Anopheles gambiae*complex egg-stage survival in dry soil from larval development sites in western Kenya. J Am Mosq Control Assoc. 1990, 6: 105-109.PubMedGoogle Scholar - Shililu JI, Grueber WB, Mbogo CM, Githure JI, Riddiford LM, Beier JC: Development and survival of
*Anopheles gambiae*eggs in drying soil: influence of the rate of drying, egg age, and soil type. J Am Mosq Control Assoc. 2004, 20: 243-247.PubMedGoogle Scholar - Huang J, Walker ED, Vulule J, Miller JR: Daily temperature profiles in and around Western Kenyan larval habitats of
*Anopheles gambiae*as related to egg mortality. Malar J. 2006, 5: 87-10.1186/1475-2875-5-87.PubMed CentralView ArticlePubMedGoogle Scholar - Bayoh MN, Lindsay SW: Effect of temperature on the development of the aquatic stages of
*Anopheles gambiae*sensu stricto (Diptera: Culicidae). Bull Entomol Res. 2003, 93: 375-381.PubMedGoogle Scholar - Kirby MJ, Lindsay SW: Effect of temperature and inter-specific competition on the development and survival of
*Anopheles gambiae*sensu stricto and An. arabiensis larvae. Acta Trop. 2009, 109: 118-123. 10.1016/j.actatropica.2008.09.025.View ArticlePubMedGoogle Scholar - Bayoh MN, Lindsay SW: Temperature-related duration of aquatic stages of the Afrotropical malaria vector mosquito
*Anopheles gambiae*in the laboratory. Med Vet Entomol. 2004, 18: 174-179. 10.1111/j.0269-283X.2004.00495.x.View ArticlePubMedGoogle Scholar - Armstrong JA, Bransby-Williams WR: The maintenance of a colony of
*Anopheles gambiae*, with observations on the effects of changes in temperature. Bull World Health Organ. 1961, 24: 427-435.PubMed CentralPubMedGoogle Scholar - Paaijmans K, Jacobs A, Takken W, Heusinkveld B, Githeko A, Dicke M, Holtslag A: Observations and model estimates of diurnal water temperature dynamics in mosquito breeding sites in western Kenya. Hydrological Processes. 2008, 22: 4789-4801. 10.1002/hyp.7099.View ArticleGoogle Scholar
- Tsila HaM J, Foko Dadji GA: Adaptative responses of
*Anopheles gambiae*in crowding larvae conditions in laboratory. Asian J Biol Sci. 2011, 4: 259-265.View ArticleGoogle Scholar - Service MW: Mortalities of the immature stages of species B of the
*Anopheles gambiae*complex in Kenya: comparison between rice fields and temporary pools, identification of predators, and effects of insecticidal spraying. J Med Entomol. 1977, 13: 535-545.View ArticlePubMedGoogle Scholar - Gimnig JE, Ombok M, Otieno S, Kaufman MG, Vulule JM, Walker ED: Density-dependent development of
*Anopheles gambiae*(Diptera: Culicidae) larvae in artificial habitats. J Med Entomol. 2002, 39: 162-172. 10.1603/0022-2585-39.1.162.View ArticlePubMedGoogle Scholar - Mogi M, Miyagi I, Cabrera BD: Development and survival of immature mosquitoes (Diptera: Culicidae) in Philippine rice fields. J Med Entomol. 1984, 21: 283-291.View ArticlePubMedGoogle Scholar
- Koenraadt CJ, Takken W: Cannibalism and predation among larvae of the
*Anopheles gambiae*complex. Med Vet Entomol. 2003, 17: 61-66. 10.1046/j.1365-2915.2003.00409.x.View ArticlePubMedGoogle Scholar - Service MW: Studies on sampling larval populations of the
*Anopheles gambiae*complex. Bull World Health Organ. 1971, 45: 169-180.PubMed CentralPubMedGoogle Scholar - Brannstrom A, Sumpter DJ: The role of competition and clustering in population dynamics. Proc Biol Sci. 2005, 272: 2065-2072. 10.1098/rspb.2005.3185.PubMed CentralView ArticlePubMedGoogle Scholar
- Mogi M, Okazawa T, Miyagi I, Sucharit S, Tumrasvin W, Deesin T, Khamboonruang C: Development and survival of anopheline immatures (Diptera: Culicidae) in rice fields in northern Thailand. J Med Entomol. 1986, 23: 244-250.View ArticlePubMedGoogle Scholar
- NCDC:www.ncdc.noaa.gov,
- Lindsay SW, Birley MH: Climate change and malaria transmission. Ann Trop Med Parasitol. 1996, 90: 573-588.PubMedGoogle Scholar
- Kirby MJ, Lindsay SW: Responses of adult mosquitoes of two sibling species, Anopheles arabiensis and A-gambiae s.s. (Diptera : Culicidae), to high temperatures. Bull Entomol Res. 2004, 94: 441-448.PubMedGoogle Scholar
- Parham PE, Michael E: Modeling the effects of weather and climate change on malaria transmission. Environ Health Perspect. 2010, 118: 620-626.PubMed CentralView ArticlePubMedGoogle Scholar
- Essential malariology. Edited by: Warrell DA, Gilles HM. 2002, xii + 348-4Google Scholar
- Gillies MT, Wilkes TJ: A study of the age-composition of populations of
*Anopheles gambiae*Giles and A.*funestus*Giles in North-Eastern Tanzania. Bull Entomol Res. 1965, 56: 237-262. 10.1017/S0007485300056339.View ArticlePubMedGoogle Scholar - Clements A, Paterson G: The analysis of mortality and survival rates in wild populations of mosquitoes. J Appl Ecology. 1981, 18: 373-399. 10.2307/2402401.View ArticleGoogle Scholar
- Charlwood J, Smith T, Kihonda J, Heiz B, Billingsley P, Takken W: Density independent feeding success of malaria vectors (Diptera: Culicidae) in Tanzania. Bull Entomol Res. 1995, 85: 29-35. 10.1017/S0007485300051981.View ArticleGoogle Scholar
- Caswell H: Matrix population models: Construction, analysis, and interpretation. 2001, Sinauer Associates, Inc, Sunderland, 2Google Scholar
- Ferguson HM, Dornhaus A, Beeche A, Borgemeister C, Gottlieb M, Mulla MS, Gimnig JE, Fish D, Killeen GF: Ecology: a prerequisite for malaria elimination and eradication. PLoS Med. 2010, 7: e1000303-10.1371/journal.pmed.1000303.PubMed CentralView ArticlePubMedGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Please note that comments may be removed without notice if they are flagged by another user or do not comply with our community guidelines.