A spatially varying coefficient model for mapping PM10 air quality at the European scale

Particulate matter (PM) air quality in Europe has improved substantially over the past decades, but it still poses a significant threat to human health. Accurate regional scale maps of PM10 concentrations are needed for monitoring progress in mitigation strategies and monitoring compliance with statutory limit values. Chemistry transport models (CTM) use emission databases and simulate the transport and deposition of pollutants. They deliver such maps but are known to be inaccurate. A promising approach is to use geostatistics to model the relationship between the in situ observations and the CTM. This has been shown to be more accurate than using either observations or CTM's alone. This paper presents a spatially varying coefficients (SVC) geostatistical model as an extension of the standard spatially varying intercept (SVI) geostatistical model. SVC allowed the regression coefficient to vary spatially according to a covariance function, the parameters of which were estimated from the data. It was built as a Bayesian hierarchical model and implemented using Markov chain Monte Carlo. The procedure was applied to Airbase PM10 observations and LOTOS-EUROS simulated PM10 for central, southern and eastern Europe. Model-fit diagnostics showed that SVC delivered a better fit to the data than SVI. Mapping the spatially varying coefficients allowed identification of the locations where the CTM performed well or poorly. This could be used for objective CTM evaluation purposes. The posterior predictive simulations were also used to map median PM10 concentrations as well as the probability of exceeding the 50μgm-3 EU daily PM10 concentration threshold. Although posterior median prediction accuracy was similar for SVI and SVC, SVC better modelled the process and yielded narrower credible intervals. As such, SVC was more appropriate for quantifying uncertainty and for mapping threshold exceedances. The resulting maps may be used to guide air quality assessment and mitigation strategies, including those related to health impacts.

[1]  Bruce Denby,et al.  Comparison of two data assimilation methods for assessing PM10 exceedances on the European scale , 2008 .

[2]  P. Atkinson,et al.  Increased accuracy of geostatistical prediction of nitrogen dioxide in the United Kingdom with secondary data , 2004 .

[3]  Claudio Carnevale,et al.  A comparison of reanalysis techniques: applying optimal interpolation and Ensemble Kalman Filtering to improve air quality monitoring at mesoscale. , 2013, The Science of the total environment.

[4]  R. Reese.  Geostatistics for Environmental Scientists , 2001 .

[5]  Bert Brunekreef,et al.  Development of NO2 and NOx land use regression models for estimating air pollution exposure in 36 study areas in Europe - The ESCAPE project , 2013 .

[6]  C Borrego,et al.  Procedures for estimation of modelling uncertainty in air quality assessment. , 2008, Environment international.

[7]  Michael L. Stein,et al.  Local likelihood estimation for nonstationary random fields , 2011, J. Multivar. Anal..

[8]  Noel A Cressie,et al.  Statistics for Spatial Data, Revised Edition. , 1994 .

[10]  David B. Dunson,et al.  Bayesian Data Analysis , 2010 .

[11]  A. Raftery,et al.  Strictly Proper Scoring Rules, Prediction, and Estimation , 2007 .

[12]  Peter M. Atkinson,et al.  A per-pixel, non-stationary mixed model for empirical line atmospheric correction in remote sensing , 2012 .

[13]  Michael D. Moran,et al.  Comparing emission inventories and model-ready emission datasets between Europe and North America for the AQMEII project , 2012 .

[14]  Andrew O. Finley,et al.  Comparing spatially‐varying coefficients models for analysis of ecological data with non‐stationary and anisotropic residual dependence , 2011 .

[15]  A. Stein,et al.  Statistical mapping of PM10 concentrations over Western Europe using secondary information from dispersion modeling and MODIS satellite observations , 2006 .

[16]  Anna Pederzoli,et al.  Performance criteria to evaluate air quality modeling applications , 2012 .

[17]  C. F. Sirmans,et al.  Spatial Modeling With Spatially Varying Coefficient Processes , 2003 .

[18]  A. Cohen,et al.  Exposure assessment for estimation of the global burden of disease attributable to outdoor air pollution. , 2012, Environmental science & technology.

[19]  Michel Gerboles,et al.  Model quality objectives based on measurement uncertainty. Part II: NO2 and PM10 , 2013 .

[20]  P. Diggle,et al.  Model‐based geostatistics , 2007 .

[21]  Hugo Denier van der Gon,et al.  The origin of ambient particulate matter concentrations in the Netherlands , 2013 .

[22]  B. Brunekreef,et al.  Air pollution and health , 2002, The Lancet.

[23]  Renske Timmermans,et al.  The LOTOS?EUROS model: description, validation and latest developments , 2008 .

[24]  Jan van de Kassteele,et al.  A model for external drift kriging with uncertain covariates applied to air quality measurements and dispersion model output , 2006 .

[25]  Alma Hodzic,et al.  A model inter-comparison study focussing on episodes with elevated PM10 concentrations , 2008 .

[26]  Sw. Banerjee,et al.  Hierarchical Modeling and Analysis for Spatial Data , 2003 .

[27]  Andrew O. Finley,et al.  spBayes for Large Univariate and Multivariate Point-Referenced Spatio-Temporal Data Models , 2013, 1310.8192.

[28]  David J. Lunn,et al.  The BUGS Book: A Practical Introduction to Bayesian Analysis , 2013 .

[29]  Sabine Banzhaf,et al.  Impact of emission changes on secondary inorganic aerosol episodes across Germany , 2013 .

[30]  Alan E. Gelfand,et al.  Model choice: A minimum posterior predictive loss approach , 1998, AISTATS.

[31]  Leon G. Higley,et al.  Forensic Entomology: An Introduction , 2009 .

[32]  Bradley P. Carlin,et al.  Bayesian measures of model complexity and fit , 2002 .

[33]  D. Dockery,et al.  An association between air pollution and mortality in six U.S. cities. , 1993, The New England journal of medicine.

[34]  Timo Mäkelä,et al.  Intercomparison of methods to measure the mass concentration of the atmospheric aerosol during INTERCOMP2000: influence of instrumentation and size cuts , 2004 .

[35]  Spectral tempering to model non-stationary variation of soil properties: Sensitivity to the initial stationary model , 2010 .

[36]  Harald Flentje,et al.  Coupling global chemistry transport models to ECMWF’s integrated forecast system , 2009 .

[37]  Michael D. Moran,et al.  Operational model evaluation for particulate matter in Europe and North America in the context of AQMEII , 2012 .

[38]  Kurt Straif,et al.  The carcinogenicity of outdoor air pollution. , 2013, The Lancet. Oncology.

[39]  Joel Schwartz,et al.  REVIEW OF EPIDEMIOLOGICAL EVIDENCE OF HEALTH EFFECTS OF PARTICULATE AIR POLLUTION , 1995 .