A large-scale stochastic spatiotemporal model for Aedes albopictus-borne chikungunya epidemiology

Chikungunya is a viral disease transmitted to humans primarily via the bites of infected Aedes mosquitoes. The virus caused a major epidemic in the Indian Ocean in 2004, affecting millions of inhabitants, while cases have also been observed in Europe since 2007. We developed a stochastic spatiotemporal model of Aedes albopictus-borne chikungunya transmission based on our recently developed environmentally-driven vector population dynamics model. We designed an integrated modelling framework incorporating large-scale gridded climate datasets to investigate disease outbreaks on Reunion Island and in Italy. We performed Bayesian parameter inference on the surveillance data, and investigated the validity and applicability of the underlying biological assumptions. The model successfully represents the outbreak and measures of containment in Italy, suggesting wider applicability in Europe. In its current configuration, the model implies two different viral strains, thus two different outbreaks, for the two-stage Reunion Island epidemic. Characterisation of the posterior distributions indicates a possible relationship between the second larger outbreak on Reunion Island and the Italian outbreak. The model suggests that vector control measures, with different modes of operation, are most effective when applied in combination: adult vector intervention has a high impact but is short-lived, larval intervention has a low impact but is long-lasting, and quarantining infected territories, if applied strictly, is effective in preventing large epidemics. We present a novel approach in analysing chikungunya outbreaks globally using a single environmentally-driven mathematical model. Our study represents a significant step towards developing a globally applicable Ae. albopictus-borne chikungunya transmission model, and introduces a guideline for extending such models to other vector-borne diseases.

Clean air in the Anthropocene

In atmospheric chemistry, interactions between air pollution, the biosphere and human health, often through reaction mixtures from both natural and anthropogenic sources, are of growing interest. Massive pollution emissions in the Anthropocene have transformed atmospheric composition to the extent that biogeochemical cycles, air quality and climate have changed globally and partly profoundly. It is estimated that mortality attributable to outdoor air pollution amounts to 4.33 million individuals per year, associated with 123 million years of life lost. Worldwide, air pollution is the major environmental risk factor to human health, and strict air quality standards have the potential to strongly reduce morbidity and mortality. Preserving clean air should be considered a human right, and is fundamental to many sustainable development goals of the United Nations, such as good health, climate action, sustainable cities, clean energy, and protecting life on land and in the water. It would be appropriate to adopt “clean air” as a sustainable development goal.

Spatiotemporal variability and contribution of different aerosol types to the Aerosol Optical Depth over the Eastern Mediterranean

This study characterizes the spatiotemporal variability and relative contribution of different types of aerosols to the Aerosol Optical Depth (AOD) over the Eastern Mediterranean as derived from MODIS Terra (3/2000-12/2012) and Aqua (7/2002-12/2012) satellite instruments. For this purpose, a 0.1° × 0.1° gridded MODIS dataset was compiled and validated against sunphotometric observations from the AErosol RObotic NETwork (AERONET). The high spatial resolution and long temporal coverage of the dataset allows for the determination of local hot spots like megacities, medium sized cities, industrial zones, and power plant complexes, seasonal variabilities, and decadal averages. The average AOD at 550 nm (AOD₅₅₀) for the entire region is ~ 0.22 ± 0.19 with maximum values in summer and seasonal variabilities that can be attributed to precipitation, photochemical production of secondary organic aerosols, transport of pollution and smoke from biomass burning in Central and Eastern Europe, and transport of dust from the Sahara Desert and the Middle East. The MODIS data were analyzed together with data from other satellite sensors, reanalysis projects and a chemistry-aerosol-transport model using an optimized algorithm tailored for the region and capable of estimating the contribution of different aerosol types to the total AOD₅₅₀. The spatial and temporal variability of anthropogenic, dust and fine mode natural aerosols over land and anthropogenic, dust and marine aerosols over the sea is examined. The relative contribution of the different aerosol types to the total AOD₅₅₀ exhibits a low/high seasonal variability over land/sea areas, respectively. Overall, anthropogenic aerosols, dust and fine mode natural aerosols account for ~ 51 %, ~ 34 % and ~ 15 % of the total AOD₅₅₀ over land, while, anthropogenic aerosols, dust and marine aerosols account ~ 40 %, ~ 34 % and ~ 26 % of the total AOD₅₅₀ over the sea, based on MODIS Terra and Aqua observations.

Spatiotemporal variability and contribution of different aerosol types to the Aerosol Optical Depth over the Eastern Mediterranean

This study characterizes the spatiotemporal variability and relative contribution of different types of aerosols to the Aerosol Optical Depth (AOD) over the Eastern Mediterranean as derived from MODIS Terra (3/2000-12/2012) and Aqua (7/2002-12/2012) satellite instruments. For this purpose, a 0.1° × 0.1° gridded MODIS dataset was compiled and validated against sunphotometric observations from the AErosol RObotic NETwork (AERONET). The high spatial resolution and long temporal coverage of the dataset allows for the determination of local hot spots like megacities, medium sized cities, industrial zones, and power plant complexes, seasonal variabilities, and decadal averages. The average AOD at 550 nm (AOD₅₅₀) for the entire region is ~ 0.22 ± 0.19 with maximum values in summer and seasonal variabilities that can be attributed to precipitation, photochemical production of secondary organic aerosols, transport of pollution and smoke from biomass burning in Central and Eastern Europe, and transport of dust from the Sahara Desert and the Middle East. The MODIS data were analyzed together with data from other satellite sensors, reanalysis projects and a chemistry-aerosol-transport model using an optimized algorithm tailored for the region and capable of estimating the contribution of different aerosol types to the total AOD₅₅₀. The spatial and temporal variability of anthropogenic, dust and fine mode natural aerosols over land and anthropogenic, dust and marine aerosols over the sea is examined. The relative contribution of the different aerosol types to the total AOD₅₅₀ exhibits a low/high seasonal variability over land/sea areas, respectively. Overall, anthropogenic aerosols, dust and fine mode natural aerosols account for ~ 51 %, ~ 34 % and ~ 15 % of the total AOD₅₅₀ over land, while, anthropogenic aerosols, dust and marine aerosols account ~ 40 %, ~ 34 % and ~ 26 % of the total AOD₅₅₀ over the sea, based on MODIS Terra and Aqua observations.

Implementing the US air quality standard for PM2.5 worldwide can prevent millions of premature deaths per year

BACKGROUND

Air pollution by fine aerosol particles is among the leading causes of poor health and premature mortality worldwide. The growing awareness of this issue has led several countries to implement air pollution legislation. However, populations in large parts of the world are still exposed to high levels of ambient particulate pollution. The main aim of this work is to evaluate the potential impact of implementing current air quality standards for fine particulate matter (PM2.5) in the European Union (EU), United States (US) and other countries where PM2.5 levels are high.

METHODS

We use a high-resolution global atmospheric chemistry model combined with epidemiological concentration response functions to investigate premature mortality attributable to PM2.5 in adults ≥30 years and children <5 years. We perform sensitivity studies to estimate the reductions in mortality that could be achieved if the PM2.5 air quality standards of the EU and US and other national standards would be implemented worldwide.

RESULTS

We estimate the global premature mortality by PM2.5 at 3.15 million/year in 2010. China is the leading country with about 1.33 million, followed by India with 575 thousand and Pakistan with 105 thousand per year. For the 28 EU member states we estimate 173 thousand and for the United States 52 thousand premature deaths in 2010. Based on sensitivity analysis, applying worldwide the EU annual mean standard of 25 μg/m(3) for PM2.5 could reduce global premature mortality due to PM2.5 exposure by 17 %; while within the EU the effect is negligible. With the 2012 revised US standard of 12 μg/m(3) premature mortality by PM2.5 could drop by 46 % worldwide; 4 % in the US and 20 % in the EU, 69 % in China, 49 % in India and 36 % in Pakistan. These estimates take into consideration that about 22 % of the global PM2.5 related mortality cannot be avoided due to the contribution of natural PM2.5 sources, mainly airborne desert dust and PM2.5 from wild fires.

CONCLUSIONS

Our results reflect the need to adopt stricter limits for annual mean PM2.5 levels globally, like the US standard of 12 μg/m(3) or an even lower limit to substantially reduce premature mortality in most of the world.

Large-Scale Modelling of the Environmentally-Driven Population Dynamics of Temperate Aedes albopictus (Skuse)

The Asian tiger mosquito, Aedes albopictus, is a highly invasive vector species. It is a proven vector of dengue and chikungunya viruses, with the potential to host a further 24 arboviruses. It has recently expanded its geographical range, threatening many countries in the Middle East, Mediterranean, Europe and North America. Here, we investigate the theoretical limitations of its range expansion by developing an environmentally-driven mathematical model of its population dynamics. We focus on the temperate strain of Ae. albopictus and compile a comprehensive literature-based database of physiological parameters. As a novel approach, we link its population dynamics to globally-available environmental datasets by performing inference on all parameters. We adopt a Bayesian approach using experimental data as prior knowledge and the surveillance dataset of Emilia-Romagna, Italy, as evidence. The model accounts for temperature, precipitation, human population density and photoperiod as the main environmental drivers, and, in addition, incorporates the mechanism of diapause and a simple breeding site model. The model demonstrates high predictive skill over the reference region and beyond, confirming most of the current reports of vector presence in Europe. One of the main hypotheses derived from the model is the survival of Ae. albopictus populations through harsh winter conditions. The model, constrained by the environmental datasets, requires that either diapausing eggs or adult vectors have increased cold resistance. The model also suggests that temperature and photoperiod control diapause initiation and termination differentially. We demonstrate that it is possible to account for unobserved properties and constraints, such as differences between laboratory and field conditions, to derive reliable inferences on the environmental dependence of Ae. albopictus populations.

Present and future projections of habitat suitability of the Asian tiger mosquito, a vector of viral pathogens, from global climate simulation

Climate change can influence the transmission of vector-borne diseases (VBDs) through altering the habitat suitability of insect vectors. Here we present global climate model simulations and evaluate the associated uncertainties in view of the main meteorological factors that may affect the distribution of the Asian tiger mosquito (Aedes albopictus), which can transmit pathogens that cause chikungunya, dengue fever, yellow fever and various encephalitides. Using a general circulation model at 50 km horizontal resolution to simulate mosquito survival variables including temperature, precipitation and relative humidity, we present both global and regional projections of the habitat suitability up to the middle of the twenty-first century. The model resolution of 50 km allows evaluation against previous projections for Europe and provides a basis for comparative analyses with other regions. Model uncertainties and performance are addressed in light of the recent CMIP5 ensemble climate model simulations for the RCP8.5 concentration pathway and using meteorological re-analysis data (ERA-Interim/ECMWF) for the recent past. Uncertainty ranges associated with the thresholds of meteorological variables that may affect the distribution of Ae. albopictus are diagnosed using fuzzy-logic methodology, notably to assess the influence of selected meteorological criteria and combinations of criteria that influence mosquito habitat suitability. From the climate projections for 2050, and adopting a habitat suitability index larger than 70%, we estimate that approximately 2.4 billion individuals in a land area of nearly 20 million km² will potentially be exposed to Ae. albopictus. The synthesis of fuzzy-logic based on mosquito biology and climate change analysis provides new insights into the regional and global spreading of VBDs to support disease control and policy making.

The contribution of outdoor air pollution sources to premature mortality on a global scale

Assessment of the global burden of disease is based on epidemiological cohort studies that connect premature mortality to a wide range of causes, including the long-term health impacts of ozone and fine particulate matter with a diameter smaller than 2.5 micrometres (PM2.5). It has proved difficult to quantify premature mortality related to air pollution, notably in regions where air quality is not monitored, and also because the toxicity of particles from various sources may vary. Here we use a global atmospheric chemistry model to investigate the link between premature mortality and seven emission source categories in urban and rural environments. In accord with the global burden of disease for 2010 (ref. 5), we calculate that outdoor air pollution, mostly by PM2.5, leads to 3.3 (95 per cent confidence interval 1.61-4.81) million premature deaths per year worldwide, predominantly in Asia. We primarily assume that all particles are equally toxic, but also include a sensitivity study that accounts for differential toxicity. We find that emissions from residential energy use such as heating and cooking, prevalent in India and China, have the largest impact on premature mortality globally, being even more dominant if carbonaceous particles are assumed to be most toxic. Whereas in much of the USA and in a few other countries emissions from traffic and power generation are important, in eastern USA, Europe, Russia and East Asia agricultural emissions make the largest relative contribution to PM2.5, with the estimate of overall health impact depending on assumptions regarding particle toxicity. Model projections based on a business-as-usual emission scenario indicate that the contribution of outdoor air pollution to premature mortality could double by 2050.

Abrupt recent trend changes in atmospheric nitrogen dioxide over the Middle East

Nitrogen oxides, released from fossil fuel use and other combustion processes, affect air quality and climate. From the mid-1990s onward, nitrogen dioxide (NO2) has been monitored from space, and since 2004 with relatively high spatial resolution by the Ozone Monitoring Instrument. Strong upward NO2 trends have been observed over South and East Asia and the Middle East, in particular over major cities. We show, however, that a combination of air quality control and political factors, including economical crisis and armed conflict, has drastically altered the emission landscape of nitrogen oxides in the Middle East. Large changes, including trend reversals, have occurred since about 2010 that could not have been predicted and therefore are at odds with emission scenarios used in projections of air pollution and climate change in the early 21st century.

Heat-related cardiovascular mortality risk in Cyprus: a case-crossover study using a distributed lag non-linear model

BACKGROUND

The frequency and intensity of heat waves is projected to increase in many parts of the world, particularly in regions such as the Eastern Mediterranean and Middle East (EMME), where the warming trends are much larger than the global average. The relationship between air temperature and premature mortality is widely recognized, however, it is not well defined in the aforementioned region. The objective of this study is to assess the relationship between cardiovascular mortality risk and air temperature in Cyprus, an island located centrally in the EMME.

METHODS

Daily cardiovascular mortality data and spatially aggregated daily mean, maximum, and minimum temperatures for the period 2004-2010 were analyzed using a case-crossover design combined with a distributed lag non-linear model.

RESULTS

A relationship between high temperatures and cardiovascular mortality was observed for cerebrovascular diseases, ischaemic and other heart diseases; this relationship was exacerbated on days with high temperatures. The highest relative risk was observed on the day of the heat event and remained significantly elevated for another day. The results were consistent regardless whether the minimum, maximum, or mean temperatures were used, although the association seems to be more pronounced with the mean temperatures, which suggests that consecutive high day- and night-time temperatures are the most hazardous.

CONCLUSIONS

The identification of a positive relationship between high temperatures and cardiovascular mortality in Cyprus raises concerns. In view of the projected climate changes and strong increases in extreme heat events in the region, appropriate interventions need to be developed.

Stratospheric sulfur and its implications for radiative forcing simulated by the chemistry climate model EMAC

Multiyear simulations with the atmospheric chemistry general circulation model EMAC with a microphysical modal aerosol module at high vertical resolution demonstrate that the sulfur gases COS and SO₂, the latter from low-latitude and midlatitude volcanic eruptions, predominantly control the formation of stratospheric aerosol. Marine dimethyl sulfide (DMS) and other SO₂ sources, including strong anthropogenic emissions in China, are found to play a minor role except in the lowermost stratosphere. Estimates of volcanic SO₂ emissions are based on satellite observations using Total Ozone Mapping Spectrometer and Ozone Monitoring Instrument for total injected mass and Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) on Envisat or Stratospheric Aerosol and Gases Experiment for the spatial distribution. The 10 year SO₂ and COS data set of MIPAS is also used for model evaluation. The calculated radiative forcing of stratospheric background aerosol including sulfate from COS and small contributions by DMS oxidation, and organic aerosol from biomass burning, is about 0.07W/m². For stratospheric sulfate aerosol from medium and small volcanic eruptions between 2005 and 2011 a global radiative forcing up to 0.2W/m² is calculated, moderating climate warming, while for the major Pinatubo eruption the simulated forcing reaches 5W/m², leading to temporary climate cooling. The Pinatubo simulation demonstrates the importance of radiative feedback on dynamics, e.g., enhanced tropical upwelling, for large volcanic eruptions.

On the temporal and spatial variation of ozone in Cyprus

More than sixteen years (1997-2013) of continuous ozone concentrations at the rural Agia Marina (EMEP, 532 ma.s.l.) station in Cyprus, together with a number of ancillary chemical and meteorological parameters have been analyzed on a multiannual, annual and diurnal basis. The observations reveal a) the presence of a prominent seasonality with maxima observed during summer (54±5 ppbv) and the minima in winter (39±3 ppbv) b) a relatively small diurnal variability with the noon levels (50±9 ppbv) being higher by ~4 pbbv compared to nighttime (46±9 ppbv) and c) a non-significant upward trend over the 16 years of 0.11±0.12 ppbv y(-1). To assess the spatial variability over Cyprus, simultaneous measurements in 2011-2012 have been performed at Inia, Stavrovouni and Cavo Greco, three remote marine monitoring sites located to the west, central and the east of the Island, respectively. Our results show that ambient ozone levels over Cyprus are mostly influenced by regional/transported ozone while the local precursor emissions play a minor role in ozone formation. On an annual basis a net ozone reduction of 1.5 and 1.0 ppbv occurs when the air masses originate from northerly and westerly directions, respectively, while this is 2.4 ppbv during southerly wind. This suggests continuous net ozone loss controlled by surface deposition and photochemical destruction, and highlights the importance of long-range transport in controlling ozone levels in Cyprus.

The role of environmental variables on Aedes albopictus biology and chikungunya epidemiology

Aedes albopictus is a vector of dengue and chikungunya viruses in the field, along with around 24 additional arboviruses under laboratory conditions. As an invasive mosquito species, Ae. albopictus has been expanding in geographical range over the past 20 years, although the poleward extent of mosquito populations is limited by winter temperatures. Nonetheless, population densities depend on environmental conditions and since global climate change projections indicate increasing temperatures and altered patterns of rainfall, geographic distributions of previously tropical mosquito species may change. Although mathematical models can provide explanatory insight into observed patterns of disease prevalence in terms of epidemiological and entomological processes, understanding how environmental variables affect transmission is possible only with reliable model parameterisation, which, in turn, is obtained only through a thorough understanding of the relationship between mosquito biology and environmental variables. Thus, in order to assess the impact of climate change on mosquito population distribution and regions threatened by vector-borne disease, a detailed understanding (through a synthesis of current knowledge) of the relationship between climate, mosquito biology, and disease transmission is required, but this process has not yet been undertaken for Ae. albopictus. In this review, the impact of temperature, rainfall, and relative humidity on Ae. albopictus development and survival are considered. Existing Ae. albopictus populations across Europe are mapped with current climatic conditions, considering whether estimates of climatic cutoffs for Ae. albopictus are accurate, and suggesting that environmental thresholds must be calibrated according to the scale and resolution of climate model outputs and mosquito presence data.

Direct kinetic study of OH and O3 formation in the reaction of CH3C(O)O2 with HO2

The reaction between HO2 and CH3C(O)O2 has three exothermic product channels, forming OH (R3a), peracetic acid (R3b), and acetic acid plus O3 (R3c). The branching ratios of the OH- and ozone-forming reaction channels were determined using a combination of laser-induced fluorescence (LIF, for time-resolved OH concentration measurement) and transient absorption spectroscopy (TAS, for time-resolved O3 concentration measurement) following pulsed laser generation of HO2 and CH3C(O)O2 from suitable precursors. TAS was also used to determine the initial concentration of the reactant peroxy radicals. The data were evaluated by numerical simulation using kinetic models of the measured concentration profiles; a Monte Carlo approach was used to estimate the uncertainties of the rate constants (k3) and branching ratios (α) thus obtained. The reaction channel forming OH (R3a) was found to be the most important with α3a = 0.61 ± 0.09 and α3c = 0.16 ± 0.08. The overall rate coefficient of the title reaction was found to be k3 = (2.1 ± 0.4) × 10(-11) cm(3) molecule(-1) s(-1) for both HO2 and DO2. Use of DO2 resulted in an increase in α3a to 0.80 ± 0.14. Comparison with former studies shows that OH formation via (R3) has been underestimated significantly to date. Possible reasons for these discrepancies and atmospheric implications are discussed.

Direct observation of OH formation from stabilised Criegee intermediates

The first direct observation of OH radicals from syn-CH₃CHOO Criegee intermediates decomposition suggests it affects tropospheric chemistry and ambient OH measurements.