Vegetation and peatland fires cause poor air quality and thousands of premature deaths across densely populated regions in Equatorial Asia. Strong El-Niño and positive Indian Ocean Dipole conditions are associated with an increase in the frequency and intensity of wildfires in Indonesia and Borneo, enhancing population exposure to hazardous concentrations of smoke and air pollutants. Here we investigate the impact on air quality and population exposure of wildfires in Equatorial Asia during Fall 2015, which were the largest over the past two decades. We performed high-resolution simulations using the Weather Research and Forecasting model with Chemistry based on a new fire emission product. The model captures the spatio-temporal variability of extreme pollution episodes relative to space- and ground-based observations and allows for identification of pollution sources and transport over Equatorial Asia. We calculate that high particulate matter concentrations from fires during Fall 2015 were responsible for persistent exposure of 69 million people to unhealthy air quality conditions. Short-term exposure to this pollution may have caused 11,880 (6,153–17,270) excess mortalities. Results from this research provide decision-relevant information to policy makers regarding the impact of land use changes and human driven deforestation on fire frequency and population exposure to degraded air quality.

Background Exposure to poor air quality leads to increased premature mortality from cardiovascular and respiratory diseases. Among the far-reaching implications of the ongoing COVID-19 pandemic, a substantial improvement in air quality was observed worldwide after the lockdowns imposed by many countries. We aimed to assess the implications of different lockdown measures on air pollution levels in Europe and China, as well as the short-term and long-term health impact. Methods For this modelling study, observations of fine particulate matter (PM 2·5 ) concentrations from more than 2500 stations in Europe and China during 2016–20 were integrated with chemical transport model simulations to reconstruct PM 2·5 fields at high spatiotemporal resolution. The health benefits, expressed as short-term and long-term avoided mortality from PM 2·5 exposure associated with the interventions imposed to control the COVID-19 pandemic, were quantified on the basis of the latest epidemiological studies. To explore the long-term variability in air quality and associated premature mortality, we built different scenarios of economic recovery (immediate or gradual resumption of activities, a second outbreak in autumn, and permanent lockdown for the whole of 2020). Findings The lockdown interventions led to a reduction in population-weighted PM 2·5 of 14·5 μg m −3 across China (−29·7%) and 2·2 μg m −3 across Europe (−17·1%), with unprecedented reductions of 40 μg m −3 in bimonthly mean PM 2·5 in the areas most affected by COVID-19 in China. In the short term, an estimated 24 200 (95% CI 22 380–26 010) premature deaths were averted throughout China between Feb 1 and March 31, and an estimated 2190 (1960–2420) deaths were averted in Europe between Feb 21 and May 17. We also estimated a positive number of long-term avoided premature fatalities due to reduced PM 2·5 concentrations, ranging from 76 400 (95% CI 62 600–86 900) to 287 000 (233 700–328 300) for China, and from 13 600 (11 900–15 300) to 29 500 (25 800–33 300) for Europe, depending on the future scenarios of economic recovery adopted. Interpretation These results indicate that lockdown interventions led to substantial reductions in PM 2·5 concentrations in China and Europe. We estimated that tens of thousands of premature deaths from air pollution were avoided, although with significant differences observed in Europe and China. Our findings suggest that considerable improvements in air quality are achievable in both China and Europe when stringent emission control policies are adopted. Funding None.

Global climate models aim to reproduce physical processes on a global scale and predict quantities such as temperature given some forcing inputs. We consider climate ensembles made of collections of such runs with different initial conditions and forcing scenarios. The purpose of this work is to show how the simulated temperatures in the ensemble can be reproduced (emulated) with a global space/time statistical model that addresses the issue of capturing nonstationarities in latitude more effectively than current alternatives in the literature. The model we propose leads to a computationally efficient estimation procedure and, by exploiting the gridded geometry of the data, we can fit massive data sets with millions of simulated data within a few hours. Given a training set of runs, the model efficiently emulates temperature for very different scenarios and therefore is an appealing tool for impact assessment.

In multivariate or spatial extremes, inference for max-stable processes observed at a large collection of locations is a very challenging problem in computational statistics, and current approaches typically rely on less expensive composite likelihoods constructed from small subsets of data. In this work, we explore the limits of modern state-of-the-art computational facilities to perform full likelihood inference and to efficiently evaluate high-order composite likelihoods. With extensive simulations, we assess the loss of information of composite likelihood estimators with respect to a full likelihood approach for some widely-used multivariate or spatial extreme models, we discuss how to choose composite likelihood truncation to improve the efficiency, and we also provide recommendations for practitioners. This article has supplementary material online.

The authors describe a new approach for emulating the output of a fully coupled climate model under arbitrary forcing scenarios that is based on a small set of precomputed runs from the model. Temperature and precipitation are expressed as simple functions of the past trajectory of atmospheric CO 2 concentrations, and a statistical model is fit using a limited set of training runs. The approach is demonstrated to be a useful and computationally efficient alternative to pattern scaling and captures the nonlinear evolution of spatial patterns of climate anomalies inherent in transient climates. The approach does as well as pattern scaling in all circumstances and substantially better in many; it is not computationally demanding; and, once the statistical model is fit, it produces emulated climate output effectively instantaneously. It may therefore find wide application in climate impacts assessments and other policy analyses requiring rapid climate projections.

We introduce a nonstationary spatio-temporal model for gridded data on the sphere. The model specifies a computationally convenient covariance structure that depends on heterogeneous geography. Widely used statistical models on a spherical domain are nonstationary for different latitudes, but stationary at the same latitude (axial symmetry). This assumption has been acknowledged to be too restrictive for quantities such as surface temperature, whose statistical behavior is influenced by large scale geographical descriptors such as land and ocean. We propose an evolutionary spectrum approach that is able to account for different regimes across the Earth's geography, and results in a more general and flexible class of models that vastly outperforms axially symmetric models and captures longitudinal patterns that would otherwise be assumed constant. The model can be estimated with a multi-step conditional likelihood approximation that preserves the nonstationary features while allowing for easily distributed computations: we show how the model can be fit to more than 20 million data points in less than one day on a state-of-the-art workstation. The resulting estimates from the statistical model can be regarded as a synthetic description (i.e. a compression) of the space-time characteristics of an entire initial condition ensemble.

One of the main challenges when working with modern climate model ensembles is the increasingly larger size of the data produced, and the consequent difficulty in storing large amounts of spatio-temporally resolved information. Many compression algorithms can be used to mitigate this problem, but since they are designed to compress generic scientific data sets, they do not account for the nature of climate model output and they compress only individual simulations. In this work, we propose a different, statistics-based approach that explicitly accounts for the space-time dependence of the data for annual global three-dimensional temperature fields in an initial condition ensemble. The set of estimated parameters is small (compared to the data size) and can be regarded as a summary of the essential structure of the ensemble output; therefore, it can be used to instantaneously reproduce the temperature fields in an ensemble with a substantial saving in storage and time. The statistical model exploits the gridded geometry of the data and parallelization across processors. It is therefore computationally convenient and allows to fit a non-trivial model to a data set of one billion data points with a covariance matrix comprising of 10 18 entries.

No abstract