Spatiotemporal analysis of the PM10 concentration over the Taranto area

In this paper, an analysis of air quality data is provided for the municipal area of Taranto (southern Italy) characterized by high environmental risks as formally decreed by the Italian government in the 1990s with two administrative measures. This is due to the massive presence of industrial sites with elevated environmental impact activities along the NW boundary of the city conurbation. The aforementioned activities have effects on the environment and on public health, as a number of epidemiological researches concerning this area reconfirm. The present study is focused on particulate matter as measured by PM10 concentrations at 13 monitoring stations, equipped with analogous instruments based on the Beta absorption technology, either reporting hourly, two-hourly, or daily measurements. Daily estimates of the PM10 concentration surfaces are obtained in order to identify areas of higher concentration (hot spots), possibly related to specific anthropic activities. Preliminary analysis involved addressing several data problems: (1) due to the use of two different validation techniques, a calibration procedure was devised to allow for data comparability; (2) imputation techniques were considered to cope with the large number of missing data, due to both different working periods and occasional malfunctions of PM10 sensors; and (3) reliable weather covariates (wind speed and direction, pressure, temperature, etc.) were obtained and considered within the analysis. Spatiotemporal modelling was addressed by a Bayesian kriging-based model proposed by Le and Zidek (2006) characterized by the use of time varying covariates and a semiparametric covariance structure. Advantages and disadvantages of the model are highlighted and assessed in terms of fit and performance. Estimated daily PM10 concentration surfaces are suitable for the interpretation of time trends and for identifying concentration peaks within the urban area.