Next Article in Journal
A Novel Scheme for Merging Active and Passive Satellite Soil Moisture Retrievals Based on Maximizing the Signal to Noise Ratio
Previous Article in Journal
Assessment of Urban Land Surface Temperature and Vertical City Associated with Dengue Incidences
Previous Article in Special Issue
Estimation of Hourly near Surface Air Temperature Across Israel Using an Ensemble Model
 
 
Correction published on 9 September 2021, see Remote Sens. 2021, 13(18), 3588.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Satellite-Based Spatio-Temporal Machine Learning Model to Reconstruct Daily PM2.5 Concentrations across Great Britain

1
Department of Public Health, Environments and Society, London School of Hygiene & Tropical Medicine, London WC1H 9SH, UK
2
The Centre on Climate Change and Planetary Health, London School of Hygiene & Tropical Medicine, London WC1H 9SH, UK
3
European Centre for Medium-Range Weather Forecast (ECMWF), Shinfield Rd, Reading RG2 9AX, UK
4
Institute of Social and Preventive Medicine, University of Bern, 3012 Bern, Switzerland
5
Oeschger Center for Climate Change Research, University of Bern, 3012 Bern, Switzerland
6
Department of Epidemiology, Lazio Regional Health Service, 00147 Rome, Italy
7
Swiss Tropical and Public Health Institute, Socinstrasse 57, 4051 Basel, Switzerland
8
University of Basel, Petersplatz 1, 4051 Basel, Switzerland
9
Department of Geography and Environmental Development, Ben-Gurion University of the Negev, Beer Sheva P.O. Box 653, Israel
10
UK Centre for Ecology & Hydrology, Bush Estate, Penicuik, Edinburgh, Midlothian EH26 0QB, UK
11
University of Exeter Medical School, Knowledge Spa, Truro TR1 3HD, UK
12
Centre for Statistical Methodology, London School of Hygiene & Tropical Medicine, London WC1E 7HT, UK
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(22), 3803; https://0-doi-org.brum.beds.ac.uk/10.3390/rs12223803
Submission received: 18 September 2020 / Revised: 30 October 2020 / Accepted: 2 November 2020 / Published: 20 November 2020

Abstract

:
Epidemiological studies on the health effects of air pollution usually rely on measurements from fixed ground monitors, which provide limited spatio-temporal coverage. Data from satellites, reanalysis, and chemical transport models offer additional information used to reconstruct pollution concentrations at high spatio-temporal resolutions. This study aims to develop a multi-stage satellite-based machine learning model to estimate daily fine particulate matter (PM2.5) levels across Great Britain between 2008–2018. This high-resolution model consists of random forest (RF) algorithms applied in four stages. Stage-1 augments monitor-PM2.5 series using co-located PM10 measures. Stage-2 imputes missing satellite aerosol optical depth observations using atmospheric reanalysis models. Stage-3 integrates the output from previous stages with spatial and spatio-temporal variables to build a prediction model for PM2.5. Stage-4 applies Stage-3 models to estimate daily PM2.5 concentrations over a 1 km grid. The RF architecture performed well in all stages, with results from Stage-3 showing an average cross-validated R2 of 0.767 and minimal bias. The model performed better over the temporal scale when compared to the spatial component, but both presented good accuracy with an R2 of 0.795 and 0.658, respectively. These findings indicate that direct satellite observations must be integrated with other satellite-based products and geospatial variables to derive reliable estimates of air pollution exposure. The high spatio-temporal resolution and the relatively high precision allow these estimates (approximately 950 million points) to be used in epidemiological analyses to assess health risks associated with both short- and long-term exposure to PM2.5.

Graphical Abstract

1. Introduction

The World Health Organization estimates in seven million global deaths associated with air pollution (both outdoor and household) every year, emphasising that exposure to particulate matter (PM) is among the greatest causes of concern [1]. Fine particles with an aerodynamic diameter smaller than 2.5 µm (PM2.5) can penetrate the human circulatory system through the lungs and provoke multiple adverse health outcomes, including mortality [2], hospital admissions [3], lung dysfunction [4], cardiovascular diseases [5], and allergic reactions [6]. Usually, epidemiological studies collect air quality (AQ) data from ground monitors to quantify both short-term and long-term PM2.5 exposure associated with acute and chronic health effects, respectively. The limitation in this health assessment approach is the lack of continuous temporal records of PM2.5 and the limited spatial distribution of the monitors. Great Britain is an example of countries with very limited spatio-temporal coverage of PM2.5, whereby the monitoring network is densely located only in major cities and widespread measurements of PM2.5 only started from 2010.
Remote sensing observations of aerosol optical depth (AOD) obtained from satellites, which measures how much direct sunlight has been scattered and absorbed by aerosol particles suspended in the atmosphere, has recently been proposed as an alternative to measuring PM variability for epidemiological purposes [7]. However, while offering the advantage of global coverage and relatively high spatio-temporal resolution, the use of AOD for PM2.5 exposure assessments presents limitations, for instance the fact that it represents the total atmospheric column concentration of the aerosol rather than surface values [8]. Unsurprisingly, early studies based only on satellite-AOD achieved very low performances in predicting PM2.5 [9]. Recent studies have proposed more sophisticated approaches, combining AOD measures with information from other satellite products, reanalysis data, chemical transport models, and geospatial features to improve the prediction of PM2.5. Such studies used various analytical methods, including multiple linear regression [10], land-use regression [11,12], and mixed effect models [13,14,15,16,17]. The last development in this research area is represented by the application of machine learning (ML) algorithms, including various architectures such as random forests [18,19,20,21,22], neural networks [19,21], and gradient boosting [19,23,24]. These have demonstrated higher performances, linked with an ability to model any kind of predictor(s)-response association and to deal better with the potentially complex relationship between PM2.5, spatial, and spatio-temporal predictors [19,25].
The aim of this study is to develop and apply a multi-stage satellite-based ML model to estimate daily concentrations of PM2.5 over a 1 km grid across Great Britain in the period 2008–2018. The analysis is based on a dataset with synchronised information from various data sources, such as several remote sensing satellite products, multiple climate and atmospheric reanalysis databases, chemical transport models, and spatial and spatio-temporal variables. The model is assessed through measures of predictive performance, error, and bias, obtained through cross-validation.

2. Materials and Methods

2.1. Study Area and Period

Great Britain is an island with an extension of 229,462 km2 surrounded by the Atlantic Ocean, Irish Sea, North Sea, and the English Channel. It comprises the countries of England, Scotland, and Wales with a total population in 2018 of almost 65 million [26]. According to the Köppen climate classification, the United Kingdom (Great Britain and Northern Ireland) is defined as having a warm temperate climate, fully humid with a mostly warm summer (cold summer for some parts of Scotland and England) [27]. The study area included 234.429 1 km grid cells (containing a unique identification code, cell-ID) from the original 1 km Great Britain National Grid Squares [28] for a period between 1 January 2008 and 31 December 2018.

2.2. PM2.5 and PM10 Observed Data

Daily PM10 and PM2.5 (µg/m3) measurements in the study period were obtained from five monitoring network sources through the R package openair [29]: Automatic Urban and Rural Network (AURN), Air Quality England (AQE), King College London (KCL), Scotland Air Quality Network (SAQN), and Wales Air Quality Network (WAQN). When monitors from different sources showed the same temporal distributions (i.e., a correlation equal to 1) and were located at approximately the same coordinates, only the AURN monitor was kept. Monitors with less than 18 hours of PM2.5 records per day as well as less than 30 days by year were removed. The final set includes 581 and 183 monitors measuring PM10 and PM2.5 along the study period, respectively. Figure 1 shows the locations of these monitors across Great Britain, illustrating how the network coverage is densely located in major cities, leaving several small cities and rural areas with only a few or no AQ records. Each monitor was indexed using the cell-ID of the 1 km grid cell that contained it.

2.3. Spatially-Lagged and Nearest Monitor PM2.5 Variables

Four spatio-temporal variables were generated from the monitor series of PM2.5 to represent spatially-lagged annual average concentrations. Monitor types were grouped into two classes (background (urban, suburban, and rural) and hotspots (traffic and industrial)) to compute the annual average values of nearby monitors by class using an inverse-distance weighted leave-one-out cross-validated (IDW-LOOCV) approach. Two different weights were applied, namely the inverse distance and the inverse squared distance (in km). The former assigns relatively more weight to distant monitors and therefore can represent a regional background, while the latter captures local differences. Therefore, these four IDW-LOOCV variables were named as follows: (i) Spatially-lagged hotspot-PM2.5 regional, (ii) Spatially-lagged background-PM2.5 regional, (iii) Spatially-lagged hotspot-PM2.5 local, and (iv) Spatially-lagged background-PM2.5 local. To improve the model performance in the spatial domain, two additional spatial variables were generated based on the closest Euclidean distance for each monitor class, named as: (i) Nearest hotspot monitor distance and (ii) Nearest background monitor distance. These six variables were used as additional predictors to capture the heterogeneity across monitors and exploit their spatial autocorrelation, and thus help the model to better categorise the differences in the spatial patterns of measured-PM2.5 series.

2.4. AOD Data: Satellite and Atmospheric Reanalysis Models

Daily satellite-AOD was obtained from the Collection 6 Level-2 gridded product (MCD19A2). These data are generated at a 1 km grid through the Multi-angle Implementation of Atmospheric Correction (MAIAC) algorithm using data from a Moderate Resolution Imaging Spectroradiometer (MODIS) sensor on board both Terra and Aqua Earth Observation satellites [30]. Four layers from the MCD19A2 product were extracted: (i) AOD Blue band (0.47 µm), (ii) AOD Green band (0.55 µm), (iii) AOD uncertainty (i.e., the level of uncertainty based on blue-band surface brightness (reflectance)), (iv) AOD_QA (quality assurance flags to retrieve only the best quality AOD). These layers are generated for each passing time of Terra and Aqua satellites over the area of study and combined by day. The layers AOD-0.47 µm and AOD-0.55 µm were used as the outcome variables after their values were filtered using AOD uncertainty and AOD_QA layers to guarantee high product quality.
This calibration process, together with pixels covered by clouds, removed a large sample of AOD grid cells over Great Britain. To fill the satellite-AOD gaps, the modelled-AOD total column was used from Copernicus Atmosphere Monitoring Service (CAMS) reanalysis provided by the European Centre for Medium-Range Weather Forecasts (ECMWF) [31]. CAMS reanalysis provides every 3-hourly modelled-AOD at five different wavelengths (0.47 µm, 0.55 µm, 0.67 µm, 0.87 µm, and 1.24 µm) with a spatial resolution of approximately 80 km, but the data were downloaded at 10 km, based on an interpolation performed through the ECMWF’s API request. The satellite-AOD and CAMS modelled-AOD were indexed to the closest 1 km grid cell from their pixel centroid.

2.5. Other Spatio-Temporal Predictors

2.5.1. Modelled PM2.5 from Chemical Transport Models

Atmospheric chemistry transport models (ACTMs) incorporate anthropogenic and natural sources of emission, land use, and meteorological conditions to simulate the atmospheric compositions and deposition of various air pollutants (trace gases and particles). Based on the European Modelling and Evaluation Programme (EMEP) ACTM, EMEP4UK, has been developed to represent the UK hourly atmospheric composition at a spatial resolution of approximately 5 km [32]. The description of the EMEP4UK model framework and setup can be found elsewhere [33,34]. Daily EMEP4UK simulations of PM2.5 (µg/m3) concentrations at the surface-level were included to represent ground-level contributions, in contrast to AOD products that refer to the total column of aerosol concentration. Each EMEP4UK 5 km pixel was linked to the closest 1 km grid cell centroid.

2.5.2. Meteorological Variables from Climate Reanalysis Models

Meteorological variables were retrieved from the ECMWF’s climate reanalysis models with the highest spatial resolution available during 2008–2018 and at two sub-day times (0:00 and 12:00). Sea-level pressure and the boundary layer height (BLH) were downloaded from the ERA 5 global reanalysis with a spatial resolution of approximately 30 km [35]. Air temperature at 2 m height and total precipitation were obtained from the ERA 5 Land global reanalysis with a spatial resolution of approximately 9 km [36]. Relative humidity, wind direction, and wind speed were downloaded from the UERRA regional reanalysis at 5.5 km for the MESCAN-SURFEX system [37]. All meteorological variables were indexed to the closest 1 km grid cell to their centroid.

2.5.3. Normalized Difference Vegetation Index

Monthly Normalized difference vegetation index (NDVI) is used to quantify vegetation presence and it ranges from −1 to 1. The NDVI 1 km grid was obtained from MOD13A3 Version 6 Level 3, Terra-MODIS product [38]. Each NDVI pixel was indexed to the closest 1 km grid cell to its pixel centroid and the NDVI values repeated for the days inside the corresponded month.

2.6. Spatial Predictors

2.6.1. Land Variables and Night-Time Light Data from Earth Observation Satellites

Three land predictors were collected from the Copernicus Land Monitoring Service (CLMS) [39] database: elevation, land cover, and impervious surfaces. Elevation data were obtained from the 2011 European Digital Elevation Model (EU-DEM) version 1.1 with a spatial resolution of 25 m. The elevation values were obtained from the mean of all 25 m-pixel values located inside each 1 km grid cell. Land cover data were obtained from the 2012 CORINE Land Cover (CLC) inventory. It was derived from high-resolution ortho-rectified satellites images that mapped all land elements at a spatial resolution ranging from 5 m to 60 m and aggregated into 100 m. Nine predictors were defined by grouping the original 44 CLC classes and each predictor represents its group proportion inside each 1 km Great Britain Grid cell. The imperviousness degree is a binary raster product at a spatial resolution of 100 m, where the value 0 represents natural land cover or water surface and value 1 represents entirely artificial surfaces (i.e., built-up areas). The amount of impervious surfaces was estimated by the proportion of artificial surfaces inside each 1 km Great Britain Grid cell.
Night-time lights data were provided by the visible Infrared Imaging Radiometer Suite (VIIRS) sensor aboard the Suomi-National Polar-orbiting Partnership (Suomi-NPP) satellite. The VIIRS Day/Night band collects cloud-free average radiance values at annual and monthly composites in a spatial resolution of 750 m [40]. The 2015 night-time lights annual mean was computed based on the weighted average of VIIRS pixels inside each 1 km Great Britain Grid cell.

2.6.2. Population Density

The resident population counts data from 2011 that were collected from the Office for National Statistics (England and Wales) and the National Records of Scotland. The smallest geographic unit of the UK census is output areas (OA), with a total of 227.769 OAs polygons for Great Britain. The proportion of each OA’s area inside each 1 km Great Britain Grid cell was extracted to estimate the weighted average population by 1 km grid cell.

2.6.3. Road Density and Distance

Road density and length predictors were derived from the Ordnance Survey (OS) [41] Open Roads product, which offers a geospatial representation of Great Britain’s Road network. Three density predictors were defined by grouping the original eight OS road types, where each predictor was computed as the sum of all roads length inside each 1 km Great Britain Grid cell by group (highway, secondary, and local). Three distance predictors were defined by computing the inverse distance of each 1 km grid centroid from the closest road group (highway, secondary, and local).

2.6.4. Inverse Distance from Airports and Seashore

Information on location and size of airports was derived from the Civil Aviation Authority (CAA) [42], which collects monthly statistics about air traffic movements for more than 60 UK Airports, including aviation activities for terminal passengers, commercial flights, and cargo tonnage. A total of 19 airports across Great Britain were selected based on a minimum of 1% for the annual percentage of passengers at the airport from 2015–2018. For each cell of the 1 km Great Britain grid, the inverse distance from the closest airport was calculated.
The inverse distance from the seashore was computed for each 1 km grid cell using the geographical information about the boundaries of England, Wales, and Scotland provided by the UK Data Service [43].

2.7. Statistical Methods

A four-stage model was developed to obtain daily PM2.5 concentrations for all 234.429 grid cells covering Great Britain. Each stage is described in detail below. Briefly, Stage-1 applies a random forest (RF) algorithm to predict PM2.5 concentrations in monitors with only PM10 records. Stage-2 uses RF models to impute missing satellite-AOD from Terra and Aqua satellites, using modelled-AOD from CAMS. Stage-3 combines the output from Stage-1 and Stage-2 with a list of spatial and spatio-temporal synchronised predictors to estimate PM2.5 concentrations at the locations of the monitors. Stage-4 uses the Stage-3 model to predict daily PM2.5 across the whole of Great Britain.

2.7.1. Random Forest Algorithm

RF is a supervised tree-based design ML algorithm that trains an ensemble of independent decision trees (or forests) in parallel. The final model accuracy is estimated by the performance average of all decision trees. There are two main advantages of the RF architecture (known as the bagging ensemble method): first, it controls the bias-variance trade-off by feeding each tree model with two-thirds of the training set while one-third is left out for validation (i.e., out-of-bag (OOB)) [44]. When all decision trees receive the same amount and list of predictors, they become highly correlated, not solving the variance problem. Therefore, the second positive aspect of RF is that the number of predictors on each tree is less than the full list available and they are randomly selected. This approach changes the predictor placed on the top of the tree, generating different splits and internal nodes. The algorithm is then able to estimate an importance ranking by quantifying the amount of error decreased due to a split of a specific predictor [45].
The performance of the models in each stage was assessed using statistics based on OOB samples and then from a 10-fold cross-validation (CV) procedure based on monitors. In the latter, ten random groups of monitors were defined, and the complete outcome series in each group were predicted using a model fitted in the other nine. This procedure offers a measure of the true predictive ability of the RF model in locations where no monitor is available. Measures of performance were generated by regressing the OOB or cross-validated predicted values on the observed series, and computing the R2, the root mean square error (RMSE), and the intercept and slope of the prediction. These statistics were computed overall using the whole series and then separated into spatial and temporal contributions. The former was computed using the averages of predicted and observed values across the series, and it offers a measure of performance in capturing long-term average PM2.5 concentrations. The latter was computed as daily deviations from the averages, and it quantifies the temporal variability explained by the model.

2.7.2. Stage-1: Increasing PM2.5 Measurements Using Co-located PM10 Monitors

The number of PM2.5 monitors across Great Britain at the beginning of the study period was relatively low, and even if the quantity has increased substantially after 2010, most of the monitors were mostly installed in major cities. Stage-1 aims to increase the number of spatio-temporal ground-level PM2.5 references using observations from co-located monitors measuring PM10, which are available in a higher number of locations and are better distributed across Great Britain. Specifically, in this stage, we fitted an RF model in locations with both PM2.5 and PM10 measurements, using only data from the specific year-model (2008–2018). The RF model for each year y is defined as:
P M 2.5 i , t y =   f ( P M 10 i , t ,   m o n . t y p e i ,   m o n t h t ,   d o w t ,   l a t i   ,   l o n i )
where: P M 2.5 i , t y and P M 10 i , t are the target variable and main predictor, respectively, measured in year y at monitor i on day t; m o n . t y p e i is a categorical variable classifying monitor i (traffic, industrial, urban, suburban, and rural); m o n t h t and d o w t are categorical variables representing the months and day of the week of day t; and l a t i and l o n i define the coordinates of monitor i. The RF model was defined based on the best parameter setting. The optimised parameters were 500 decision trees as the RF ensemble (Ntree = 500) and 4 variables randomly selected to be used on each tree (mtry = 4). This model was eventually used to predict PM2.5 measurements in locations/days in which only PM10 was measured.

2.7.3. Stage-2: Imputing Missing Satellite-AOD from CAMS Modelled-AOD

The percentage of missing satellite-AOD measurements in Great Britain ranged between 87% and 94% during 2008–2018 with the greatest portion during autumn and winter, near to the coast, and in the North of Scotland. Stage-2 imputes satellite-AOD missing for every day and 1 km grid based on an optimised RF model (Ntree = 50 and mtry = 20) and satellite-AOD wavelength (0.47 µm and 0.55 µm), separately in each year within the study period. These RF models were built for each year y as follows:
Satellite - A O D i , t ( z , y ) =   f ( C A M S . A O D 1 , 1 i , t ,     ,   C A M S . A O D 5 , 1 i , t ,   C A M S . A O D 1 , 2 i , t ,   ,   C A M S . A O D 5 , 7 i , t ,   y d a y t ,   l a t i ,   l o n i )
where: Satellite- A O D i , t z ,   y is the target variable representing satellite-AOD (wavelength z, year y) estimates at grid cell i on day t; CAMS.AOD is the main predictor, representing CAMS modelled-AOD estimates at grid cell i, on day t, at five wavelengths (0.47 µm, 0.55 µm, 0.67 µm, 0.865 µm, and 1.24 µm), and at seven sub-day times (3 h, 6 h, 9 h, 12 h, 15 h, 18 h, and 21 h); y d a y t defines the sequence of days in a year from 1 to 365 (366, for leap years); l a t i and l o n i represent the coordinates of grid cell centroid i. This model was eventually used to predict missing satellite-AOD measurements.

2.7.4. Stage-3: Estimating PM2.5 Concentrations Using Spatial and Spatio-Temporal Variables

Stage-3 aims to build a predictive model for daily PM2.5 concentrations, using as target variable the combined set of PM2.5 directly measured from monitors or predicted from Stage-1, AOD measured from satellite instruments or predicted from Stage-2, together with all the spatially and spatio-temporally synchronized predictors described in the previous section. The RF models were fit separately by year y. The optimization process for this stage selected Ntree = 500 and mtry = 20 as the best set of parameter values. In this stage, PM2.5 concentrations were log-transformed to ensure the prediction of non-negative values. The Stage-3 model for each year y is defined as:
log ( P M 2.5 i , t y ) = f ( S P T 1 i , t , S P T 15 i , t ,   S P 1 i ,   ,   S P 27 i )
where: log ( P M 2.5 i , t y ) is the target variable representing the log-PM2.5 concentrations in year y at the monitor located in grid cell i on day t, while S P T i , t and S P i represent the spatio-temporal and spatial predictors, respectively.
Specifically, SPTi,t included: surface-level log-PM2.5 values from the EMEP4UK model; Stage-2 AOD (at 0.47 µm, 0.55 µm); meteorological variables (air temperature, sea pressure, relative humidity, total precipitation, wind speed, and direction); BLH (at 12:00 and 24:00 hours); monthly-averaged NDVI; and day of the year, month, and week. SPi included: spatially–lagged regional and local log-PM2.5 average concentrations from background and hotspot groups; nearest background and hotspot monitor distance; land variables (elevation, land cover in 9 groups, and % of impervious surface); night-time light; road density and inverse distance (each for highway, secondary, and local); population density; and inverse distance from closest airport and seashore.

2.7.5. Stage-4: Reconstructing PM2.5 Time-Series at 1 km Grid

Using the RF models developed by year in Stage-3, daily PM2.5 concentrations for each 1 km grid cell were reconstructed across Great Britain for the whole study period (2008–2018).

3. Results

3.1. Stage-1 Results

Not all monitors in each network provided the full series of daily PM concentrations betwen 2008–2018. However, PM10 was more often available, especially in the years 2008–2009, when there were less than 80 PM2.5 monitors across Great Britain. The imputation process in Stage-1 enabled the expansion from 46 to 269 in 2008 and from 65 to 278 in 2009. Table 1 shows the Stage-1 model performance (i.e., a separate RF model for each year) reported by two CV methods: (i) OOB, the original RF CV algorithm, which presented better accuracy (R2 average of 0.932), and (ii) 10-Fold CV, a targeted sampling, which leaves out monitors with the full set of observations (R2 average of 0.855). As expected, the 10-Fold CV approach resulted in a slightly lower predictive accuracy although it still showed a good performance. There is some indication of bias in the 10-Fold CV, with a slope lower than the expected value of 1 and a slightly negative intercept. The results for 10-Fold CV spatial and temporal domains are in Table S1 in Supplementary Materials.

3.2. Stage-2 Results

The Stage-2 procedure is illustrated in Figure 2, showing missing satellite-AOD, imputed through a combination of multiple modelled-AOD wavelengths and sub-day times. Table 2 shows the performance of Stage-2 models validated using the OOB CV method. The results presented consistently high R2 (ranging from 0.963 to 0.988), low RMSE (varying between 0.007 to 0.010), and almost no bias (intercept zero and slope close to 1 in almost all years and wavelengths).

3.3. Stage-3 Results

Stage-3, the main step of the satellite-based machine learning framework, combines the output of Stage-1 and Stage-2 with a list of spatial and spatio-temporal predictors to estimate PM2.5 at the locations of the monitors. The relative importance of the predictors for the Stage-3 RF models is ranked in Table 3. The list with the top-15 predictors demonstrates the larger contribution of the more informative spatio-temporal variables (EMEP4UK PM2.5, meteorological parameters, BLH, and sea-level pressure). All the proposed spatially-lagged PM2.5 variables were classified as highly important and their ranking positions varied slightly across the displayed years. This suggests the presence of spatial correlations in PM2.5 values that are not entirely captured by the other variables.
Table 4 shows the results of the 10-Fold CV Stage-3 RF models by year. The results indicate a good predictive performance of the model throughout the study period. Overall cross-validated R2 ranged from 0.704 (2008) to 0.821 (2011), with an average of 0.767. The average prediction error is 4.042 µg/m3, with a negligible bias in intercept and slope. The inspection of the spatial and temporal contributions shows that the model performs well generally across the two components, displaying a spatial performance drop only for 2008 (0.486) and 2015 (0.579). The cross-validated spatial R2 ranges from 0.486 (2008) to 0.746 (2017), while the cross-validated temporal R2 from 0.738 (2010) to 0.843 (2011). The two components have an average of 0.658 and 0.795, respectively. The high spatial R2 performance across the years demonstrates that the Stage-3 RF models were able to predict the spatial variation of long-term PM2.5 across Great Britain with good accuracy. In Supplementary Materials, Table S4 shows the Stage-3 results by season, demonstrating that the seasonal patterns were well described by RF models, although with a lower accuracy for the temporal domain in summer, characterized by a lower 10-Fold R2. This drop in performance during summer was also seen by other studies [17,20].

3.4. Stage-4 Results

Stage-4 provides the prediction of daily PM2.5 concentrations for each of the 234.429 1 km grids covering Great Britain. Results indicate an annual PM2.5 average of 9.41 µg/m3 for 2008, 10.17 µg/m3 for 2013, 8.05 µg/m3 for 2018, and 8.84 µg/m3 for 2008–2018 but with a strong spatial and temporal variation. The spatial distribution of annual average PM2.5 concentrations for 2008, 2013, and 2018 are shown in Figure 3, revealing a decrease of pollution levels in recent years across the whole territory, although slightly stronger in England. Table S5 in Supplementary Materials provides the same figures for all the years, confirming the decreasing trend. The spatial comparison suggests that PM2.5 concentrations are lower in Scotland and Wales compared to the more populated southern regions of England, with hotspots located in urban areas such as Liverpool, Manchester, Birmingham, and Greater London. At the bottom of Figure 3, the maps display the corresponding annual average of PM2.5 levels in London, demonstrating the precision of the multi-stage ML model in reconstructing PM2.5 concentrations in 1 km grid cells within urban areas. The maps show local hotspots of high pollution, with a spatial distribution that, however, changes along the study period.
The greatest contribution of this study is in the ability of the satellite-based ML models to reconstruct daily levels of PM2.5 over a 1 km grid across a wide geographical domain. Figure 4 displays the time series plots of observed and predicted PM2.5 concentrations at three monitoring sites and their locations within the geographical domain of Great Britain. While not necessarily representative, the results show that the error varies depending on location, type of monitor, and period. However, generally the plots indicate a very good performance of the ML algorithms in recovering the observed temporal variation in PM2.5 levels, capturing peaks and periods of stable low concentrations.
Figure 5 complements the analysis of daily variations by displaying the spatial distribution of PM2.5 estimations across Great Britain (Top) and London (Bottom) for specific days within the study period. It is interesting to note the wide variation in PM2.5 concentrations between days in the same area and between areas on different days. For instance, the maps in the left panels represent a day with almost no variation and generally low concentrations; the maps in the mid panels display a strong north–south split likely linked to weather conditions; the right panels show a more complex pattern with a wider range of PM2.5 values, a large area with very high pollution concentrations located in east London, and hotspots in highly populated urban areas across England.

4. Discussion

This study presents the first application of satellite-based spatio-temporal ML methods to reconstruct levels of pollution across Great Britain, providing estimates of daily PM2.5 concentrations over a 1 km grid during 2008–2018. The multi-stage ML framework provided significant advantages, allowing the combination of information from multiple data sources, such as air quality monitoring networks, remote sensing satellite products, chemical dispersion models, reanalysis databases, and administrative census data, among others.
The beginning of the study period (2008 and 2009) had a lower quantity of monitors measuring PM2.5 across Great Britain, but this number increased considerably from 2010 after a new measuring network had been established in 2009 [46]. Therefore, Stage-1 was an extremely important step to extend the number of PM2.5 measurements. The implementation of Stage-2 was also relevant to fill the gaps in MAIAC AOD retrievals from satellites, thus maximizing the available information. Stage-3 produced an ML prediction model based on a long list of informative predictors, accounting for potentially complex inter-relationships and functional forms. The Stage-3 ML algorithms offered excellent performance, showing an average cross-validated R2 of 0.767 across the period, with an increased predictive ability in the last years. Stage-4 provided a single PM2.5 estimation for each of the 234.429 1 km grid cells in each of the 4018 days, totalling around 950 million data points. The methodology provides complete spatial coverage, high resolution, and a relative small error of the predictions, and the ability to capture variations in PM2.5 concentrations across both spatial and temporal domains. The model offers a prediction accuracy that makes the output suitable for application in epidemiological studies on the short- and long-term health effects of air pollution. As demonstrated in Figure 4, the methodological approach presented in this study was able to capture the daily PM2.5 variability across different years, locations, and monitor types.
The output from the empirical ML model developed in this study complements existing databases of modelled PM2.5 in the United Kingdom, with some advantages for applications in epidemiological studies. Country-wide maps generated by emission-dispersion models are usually available at a coarser spatial or temporal resolution [32,47,48], and they generally show lower small-scale accuracy when tested against observed monitoring data [49,50]. The spatio-temporal ML models presented here demonstrated comparable predictive performance to similar methods applied in other countries, based either on single-learner ML models [13,20], ensemble ML models [18,19,21], or generalised additive models (GAM) [15]. Ref [21] developed an ensemble ML model (composed by RF, Deep Neural Networks, GAM, Gradient Boosting, K-nearest Neighbour) to estimate PM2.5 for Greater London, reaching a mean 2005–2013 CV spatial-R­2 of 0.396. Using the 2008–2013 period, this study reached a CV spatial-R2 of 0.637 for the whole of Great Britain. Modelling a larger area might have provided more information and a higher spatial variability, improving substantially the CV spatial-R­2. The application of a single learner to model air pollution in both spatial and temporal domains for Great Britain achieved a satisfactory performance. Nonetheless, ensemble model formats can be used in future applications using the same area of study and variables to assess how much performance gain is reached compared to the RF-only learner.
In 2008, a new Air Quality Directive came into force, bringing considerable changes to the following UK annual air quality assessments in 2010, setting an annual mean target of 25 µg/m3 [51]. Several policy controls and emission reductions had been put in place aiming to reduce from 2010 the road traffic PM emission by 83%, off-road mobile machinery by 54%, and energy production by 32% until 2020 [52]. Independently of the temporal aggregation (daily or annual), all maps shown in this study detected this considerable drop in the PM2.5 concentrations from 2010 across Great Britain.
Some limitations of this study must be acknowledged. First, the multi-stage model relies on the extension of the observed series of PM2.5 by predicting values based on co-located PM10 in Stage-1. This step was necessary for the application of the method in the early years characterized by sparse PM2.5 monitoring, which is likely to have contributed to the lower predictive performance in this period. Morover, the cross-validation procedures revealed the presence of some bias in the Stage-1 predictions, particularly in the spatial domain, which was probably linked to limitations in modelling the relative distribution of the two PM components. Second, while the model displayed a good performance throughout the period, the accuracy is worse in the temporal domain in the summer (Table S4), suggesting limitations in capturing the higher temporal variation in this season. Third, the generalization of the prediction model was dependent on the selected locations of the monitors, which may be not representative of the study domain. This can result in an underestimation of the error and potential biases in the predictions in more remote and less represented areas if structural spatial differences are not entirely captured by the model covariates.
The findings of this study also highlight additional limitations not yet discussed in previous studies related to the use of remote sensing satellite products for reconstructing air pollution exposure for epidemiological purposes. The RF importance ranking order showed that Stage-3 was mainly informed by EMEP4UK PM2.5, regional and local interpolated PM2.5 estimations, and meteorological variables, while the contribution from Stage-2 outputs (i.e., predicted-AOD 0.47 µm and 0.55 µm) was very limited. Therefore, direct satellite observations (i.e., AOD, NDVI, elevation, land cover, impervious surfaces, and night-time light) offered a relatively minor contribution to represent ground PM2.5 exposures (Stage-3), while satellite-based products (i.e., indirect satellite observations) such as EMEP4UK and climate reanalysis products played a much more relevant role to derive reliable PM2.5 estimates across Great Britain.
Future directions in remote sensing are pointing to new satellite instruments developed for air pollution monitoring which are likely to provide better resolution and reliability, thereby improving the predictive performance. As examples, the recently launched Copernicus Sentinel-5 Precursor [53] and the two future earth observation missions from European Space Agency (e.g., Copernicus Sentinel-4 [54] and Sentinel-5 [55]) were developed to provide information on atmospheric variables, such as air quality parameters (nitrogen dioxide, ozone, and aerosols). Finally, future research developments related to advanced statistical methodology, as demonstrated by the ML framework proposed here, include the application of geostatistical techniques, the use of alternative single-learner or ensemble ML algorithms, and statistical downscaling methods to increase further the resolution of the predictions.

5. Conclusions

This study developed and applied a multi-stage ML model, combining data for multiple sources, including remote sensing satellite products, climate and atmospheric reanalysis models, chemical transport models, and geospatial features, to generate a complete map of daily PM2.5 concentrations in a 1 km grid across Great Britain between 2008 and 2018. The model showed good performance overall and in both spatial and temporal domains, with an accuracy that is compatible with the use of such reconstructed values as a proxy for PM2.5 exposures in epidemiological studies. In particular, the availability of high-resolution measures that can be linked as such or aggregated at different spatial and/or temporal scales makes the output suitable for investigations on both transient and chronic health risks associated to short and long-term exposures to PM2.5, respectively.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-4292/12/22/3803/s1, Figure S1. Number of 1 km Great Britain Grid cells with daily satellite-AOD available for (a) 2008, (b) 2013, and (c) 2018. The horizontal red lines represent the thresholds: 10.000, 50.000, and 100.000 1 km grid cells. These data refer to satellite-AOD available after the filtering process to remove bad retrievals. Table S1 Stage-1 results for 10-Fold CV spatial and temporal domains. Table S2. Pearson correlation between Satellite-AOD 0.47 µm and five CAMS modelled-AOD wavelengths at the time 12:00. 0.47 µm, 0.55 µm, 0.67 µm, 0.865 µm, and 1.24 µm. Table S3. Relative importance (%) of the predictors in Stage-2. Table S4. Predicted-PM2.5 concentrations obtained from Stage-3 RF models by season were regressed against Stage-1 measured/predicted-PM2.5 concentrations in a linear regression model. The CV–R2 described in three different patterns (overall, spatial, and temporal), RMSE (a measure of the model error, µg/m3), intercept (µg/m3), and slope (µg/m3). In Great Britain, the season is composed of the following months: (i) Winter: December, January, and February, (ii) Spring: March, April, and May, (iii) Summer: June, July, and August, and (iv) Autumn: September, October, and November. Table S5. Distribution of annual averages PM2.5 concentrations for all years.

Author Contributions

The authors’ contribution are as follows: Conceptualization, R.S.d.S. and A.G.; methodology, R.S.d.S., A.G., I.K., M.S. and, K.H.; validation, R.S.d.S. and A.G.; formal analysis, R.S.d.S., A.G., F.S. and P.M.; data curation, R.S.d.S.; writing—original draft preparation, R.S.d.S.; writing—review and editing, all authors; visualization, R.S.d.S.; project administration, A.G.; funding acquisition, A.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Medical Research Council-UK grant ID: MR/M022625/1, the Natural Environment Research Council UK grant ID: NE/R009384/1, and the European Union’s Horizon 2020 Project Exhaustion grant ID: 820655. EMEP4UK Model results and contributions by S.R. and M.V. were funded by grant number NE/R016429/1 as part of the UK-SCAPE programme delivering National Capability.

Acknowledgments

The authors are grateful for the technical support received from the European Centre for Medium-Range Weather Forecasts (ECMWF). The authors also would like to thank the European Space Research Institute (ESRIN) from the European Space Agency (ESA) for their feedback on this study. This research had free and open access to all data sources: (i) NASA EOSDIS Land Processes Distributed Active Archive Center (LP DAAC), (ii) NOAA National Centers for Environmental Information (NCEI), (iii) Copernicus Land Monitoring Service (CLMS), (iv) Copernicus Atmosphere Data Store, (v) Copernicus Climate Data Store, (vi) EMEP4UK (except 2016-2018 which will be available later this year), (vii) UK Department for Environment, Food and Rural Affairs.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Word Health Organization (WHO). Available online: https://www.who.int/health-topics/air-pollution#tab=tab_1 (accessed on 20 March 2020).
  2. Liu, C.; Chen, R.; Sera, F.; Vicedo-Cabrera, A.M.; Guo, Y.; Tong, S.; Coelho, M.S.Z.S.; Saldiva, P.H.N.; Lavigne, E.; Matus, P. Ambient particulate air pollution and daily mortality in 652 cities. N. Engl. J. Med. 2019, 381, 705–715. [Google Scholar] [CrossRef] [PubMed]
  3. Basagaña, X.; Jacquemin, B.; Karanasiou, A.; Ostro, B.; Querol, X.; Agis, D.; Alessandrini, E.; Alguacil, J.; Artiñano, B.; Catrambone, M. Short-term effects of particulate matter constituents on daily hospitalizations and mortality in five South-European cities: Results from the MED-PARTICLES project. Environ. Int. 2015, 75, 151–158. [Google Scholar] [CrossRef] [PubMed]
  4. Raaschou-Nielsen, O.; Beelen, R.; Wang, M.; Hoek, G.; Andersen, Z.J.; Hoffmann, B.; Stafoggia, M.; Samoli, E.; Weinmayr, G.; Dimakopoulou, K. Particulate matter air pollution components and risk for lung cancer. Environ. Int. 2016, 87, 66–73. [Google Scholar] [CrossRef] [PubMed]
  5. Lavigne, E.; Lima, I.; Hatzopoulou, M.; Van Ryswyk, K.; Decou, M.L.; Luo, W.; van Donkelaar, A.; Martin, R.V.; Chen, H.; Stieb, D.M. Spatial variations in ambient ultrafine particle concentrations and risk of congenital heart defects. Environ. Int. 2019, 130, 1–7. [Google Scholar] [CrossRef] [PubMed]
  6. Lavigne, E.; Donelle, J.; Hatzopoulou, M.; Van Ryswyk, K.; van Donkelaar, A.; Martin, R.V.; Chen, H.; Stieb, D.M.; Gasparrini, A.; Crighton, E. Spatiotemporal variations in ambient ultrafine particles and the incidence of childhood asthma. Am. J. Respir. Crit. Care Med. 2019, 199, 1487–1495. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. NASA Earth Observations. Available online: https://neo.sci.gsfc.nasa.gov/view.php?datasetId=MODAL2_M_AER_OD (accessed on 20 March 2020).
  8. Van Donkelaar, A.; Martin, R.V.; Brauer, M.; Kahn, R.; Levy, R.; Verduzco, C.; Villeneuve, P.J. Global estimates of ambient fine particulate matter concentrations from satellite-based aerosol optical depth: Development and application. Environ. Health Perspect. 2010, 118, 847–855. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Koelemeijer, R.; Homan, C.; Matthijsen, J. Comparison of spatial and temporal variations of aerosol optical thickness and particulate matter over Europe. Atmos. Environ. 2006, 40, 5304–5315. [Google Scholar] [CrossRef]
  10. Gupta, P.; Christopher, S.A. Particulate matter air quality assessment using integrated surface, satellite, and meteorological products: Multiple regression approach. J. Geophys. Res. Atmos. 2009, 114, 1–14. [Google Scholar] [CrossRef] [Green Version]
  11. Beckerman, B.S.; Jerrett, M.; Martin, R.V.; van Donkelaar, A.; Ross, Z.; Burnett, R.T. Application of the deletion/substitution/addition algorithm to selecting land use regression models for interpolating air pollution measurements in California. Atmos. Environ. 2013, 77, 172–177. [Google Scholar] [CrossRef]
  12. Vienneau, D.; de Hoogh, K.; Beelen, R.; Fischer, P.; Hoek, G.; Briggs, D. Comparison of land-use regression models between Great Britain and the Netherlands. Atmos. Environ. 2010, 44, 688–696. [Google Scholar] [CrossRef]
  13. De Hoogh, K.; Héritier, H.; Stafoggia, M.; Künzli, N.; Kloog, I. Modelling daily PM2.5 concentrations at high spatio-temporal resolution across Switzerland. Environ. Pollut. 2018, 233, 1147–1154. [Google Scholar] [CrossRef] [PubMed]
  14. Kloog, I.; Koutrakis, P.; Coull, B.A.; Lee, H.J.; Schwartz, J. Assessing temporally and spatially resolved PM2.5 exposures for epidemiological studies using satellite aerosol optical depth measurements. Atmos. Environ. 2011, 45, 6267–6275. [Google Scholar] [CrossRef]
  15. Kloog, I.; Sorek-Hamer, M.; Lyapustin, A.; Coull, B.; Wang, Y.; Just, A.C.; Schwartz, J.; Broday, D.M. Estimating daily PM2.5 and PM10 across the complex geo-climate region of Israel using MAIAC satellite-based AOD data. Atmos. Environ. 2015, 122, 409–416. [Google Scholar] [CrossRef] [Green Version]
  16. Lee, H.; Liu, Y.; Coull, B.; Schwartz, J.; Koutrakis, P. A novel calibration approach of MODIS AOD data to predict PM2. 5 concentrations. Atmos. Chem. Phys. 2011, 11, 7991–8002. [Google Scholar] [CrossRef] [Green Version]
  17. Stafoggia, M.; Schwartz, J.; Badaloni, C.; Bellander, T.; Alessandrini, E.; Cattani, G.; de’ Donato, F.; Gaeta, A.; Leone, G.; Lyapustin, A. Estimation of daily PM10 concentrations in Italy (2006–2012) using finely resolved satellite data, land use variables and meteorology. Environ. Int. 2016, 99, 234–244. [Google Scholar] [CrossRef] [PubMed]
  18. Chen, G.; Li, S.; Knibbs, L.D.; Hamm, N.; Cao, W.; Li, T.; Guo, J.; Ren, H.; Abramson, M.J.; Guo, Y. A machine learning method to estimate PM2.5 concentrations across China with remote sensing meteorological and land use information. Sci. Total Environ. 2018, 636, 52–60. [Google Scholar] [CrossRef]
  19. Di, Q.; Amini, H.; Shi, L.; Kloog, I.; Silvern, R.; Kelly, J.; Sabath, M.B.; Choirat, C.; Koutrakis, P.; Lyapustin, A. An ensemble-based model of PM2.5 concentration across the contiguous United States with high spatiotemporal resolution. Environ. Int. 2019, 130, 1–13. [Google Scholar] [CrossRef]
  20. Stafoggia, M.; Bellander, T.; Bucci, S.; Davoli, M.; de Hoogh, K.; de’ Donato, F.; Gariazzo, C.; Lyapustinf, A.; Michelozzi, P.; Renzi, M. Estimation of daily PM10 and PM2.5 concentrations in Italy, 2013–2015, using a spatiotemporal land-use random-forest model. Environ. Int. 2019, 124, 170–179. [Google Scholar] [CrossRef] [PubMed]
  21. Yazdi, M.D.; Kuang, Z.; Dimakopoulou, K.; Barratt, B.; Suel, E.; Amini, H.; Lyapustin, A.; Katsouyanni, K.; Schwartz, J. Predicting Fine Particulate Matter (PM2.5) in the Greater London Area: An Ensemble Approach using Machine Learning Methods. Remote Sens. 2020, 12, 914. [Google Scholar] [CrossRef] [Green Version]
  22. Wei, J.; Huang, W.; Li, Z.; Xue, W.; Peng, Y.; Sune, L.; Cribb, M. Estimating 1-km-resolution PM2.5 concentrations across China using the space-time random forest approach. Remote Sens. Environ. 2019, 231, 1–14. [Google Scholar] [CrossRef]
  23. Chen, Z.Y.; Zhang, T.H.; Zhang, R.; Zhu, Z.M.; Yang, J.; Chen, P.Y.; Ou, C.Q.; Guo, Y. Extreme gradient boosting model to estimate PM2.5 concentrations with missing-filled satellite data in China. Atmos. Environ. 2019, 202, 180–189. [Google Scholar] [CrossRef]
  24. Zhan, Y.; Luo, Y.; Deng, X.; Chen, H.; Grieneisen, M.L.; Shen, X.; Zhu, L.; Zhang, M. Spatiotemporal prediction of continuous daily PM2.5 concentrations across China using a spatially explicit machine learning algorithm. Atmos. Environ. 2017, 155, 129–139. [Google Scholar] [CrossRef]
  25. Polley, E.C.; Rose, S.; van der Laan, M.J. Super Learning. In Targeted Learning: Causal Inference for Observational and Experimental Data; Van der Laan, M.J., Rose, S., Eds.; Springer: New York, NY, USA, 2011; pp. 43–65. [Google Scholar]
  26. Office for National Statistics (ONS). Available online: https://www.ons.gov.uk/peoplepopulationandcommunity/populationandmigration/populationestimates (accessed on 1 April 2020).
  27. Kottek, M.; Grieser, J.; Beck, C.; Rudolf, B.; Rubel, F. World Map of the Köppen-Geiger climate classification updated. Meteorol. Z. 2006, 15, 259–263. [Google Scholar] [CrossRef]
  28. Digimap. Available online: https://0-digimap-edina-ac-uk.brum.beds.ac.uk/webhelp/os/data_information/os_data_issues/grid_references.htm (accessed on 1 April 2020).
  29. Openair R Package. Available online: https://cran.r-project.org/web/packages/openair/openair.pdf (accessed on 25 May 2020).
  30. Lyapustin, A.; Wang, Y. MCD19A2 MODIS/Terra+Aqua Land Aerosol Optical Depth Daily L2G Global 1km SIN Grid V006. 2018, distributed by NASA EOSDIS Land Processes DAAC. Available online: https://0-doi-org.brum.beds.ac.uk/10.5067/MODIS/MCD19A2.006 (accessed on 28 May 2020).
  31. Bozzo, A.; Remy, S.; Benedetti, A.; Flemming, J.; Bechtold, P.; Rodwell, M.J.; Morcrette, J.J. Implementation of a CAMS-Based Aerosol Climatology in the IFSA; European Centre for Medium-Range Weather Forecasts: Reading, UK, 2017; Volume 801, pp. 1–33. Available online: https://www.ecmwf.int/sites/default/files/elibrary/2017/17219-implementation-cams-based-aerosol-climatology-ifs.pdf (accessed on 3 November 2020).
  32. European Modelling and Evaluation Programme for the UK (EMEP4UK). Available online: http://www.emep4uk.ceh.ac.uk/ (accessed on 13 July 2020).
  33. Vieno, M.; Heal, M.R.; Twigg, M.M.; MacKenzie, I.A.; Braban, C.F.; Lingard, J.J.N.; Ritchie, S.; Beck, R.C.; Móring, A.; Ots, R.; et al. The UK particulate matter air pollution episode of March–April 2014: More than Saharan dust. Environ. Res. Lett. 2016, 11, 12. [Google Scholar] [CrossRef]
  34. Vieno, M.; Dore, A.J.; Stevenson, D.S.; Doherty, R.; Heal, M.R.; Reis, S.; Hallsworth, S.; Tarrason, L.; Wind, P.; Fowler, D.; et al. Modelling surface ozone during the 2003 heat-wave in the UK. Atmos. Chem. Phys. 2010, 10, 7963–7978. [Google Scholar] [CrossRef] [Green Version]
  35. ERA 5 Global Climate Reanalysis. Available online: https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels?tab=overview (accessed on 28 May 2020).
  36. ERA 5 Land Global Climate Reanalysis. Available online: https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-land?tab=overview (accessed on 28 May 2020).
  37. UERRA Regional Reanalysis. Available online: https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-uerra-europe-soil-levels?tab=overview (accessed on 28 May 2020).
  38. Didan, K. MOD13A3 MODIS/Terra Vegetation Indices Monthly L3 Global 1 km SIN Grid V006 [Data set]. NASA EOSDIS LP DAAC, 2015. Available online: https://lpdaac.usgs.gov/products/mod13a3v006/ (accessed on 3 November 2020).
  39. Copernicus Land Monitoring Service (CLMS). Available online: https://land.copernicus.eu/pan-european (accessed on 29 May 2020).
  40. Earth Observation Group (EOG). Available online: https://ngdc.noaa.gov/eog/viirs/download_dnb_composites.html (accessed on 1 July 2020).
  41. Ordnance Survey Open Roads. Available online: https://www.ordnancesurvey.co.uk/documents/os-open-roads-user-guide.pdf. (accessed on 29 May 2020).
  42. Civil Aviation Authority (CAA). Available online: caa.co.uk/home (accessed on 29 May 2020).
  43. UK Data Service. Available online: https://www.ukdataservice.ac.uk/ (accessed on 29 May 2020).
  44. Schneider dos Santos, R. Estimating spatio-temporal air temperature in London (UK) using machine learning and earth observation satellite data. Int. J. Appl. Earth Obs. Geoinf. 2020, 88, 1–10. [Google Scholar] [CrossRef]
  45. James, G.; Witten, D.; Hastie, T.; Tibshirani, R. An introduction to statistical learning; Springer: Berlin, Germany, 2013; 430p. [Google Scholar]
  46. Department for Environment, Food & Rural Affairs (DEFRA). Fine Particulate Matter (PM2.5) in the UK 2012. Available online: https://www.gov.uk/government/publications/fine-particulate-matter-pm2-5-in-the-uk (accessed on 25 May 2020).
  47. DEFRA. Modelled Background Pollution Data. Available online: https://uk-air.defra.gov.uk/data/pcm-data (accessed on 13 July 2020).
  48. Savage, N.H.; Agnew, P.; Davis, L.S.; Ordonez, C. Air quality modelling using the Met Office Unified Model (AQUM OS24-26): Model description and initial evaluation. Geosci. Model Dev. 2013, 6, 353–372. [Google Scholar] [CrossRef] [Green Version]
  49. Hood, C.; MacKenzie, I.; Stocker, J.; Johnson, K.; Carruthers, D.; Vieno, M.; Doherty, R. Air quality simulations for London using a coupled regional-to-local modelling system. Atmos. Chem. Phys. 2018, 18, 11221–11245. [Google Scholar] [CrossRef] [Green Version]
  50. Lin, C.; Heal, M.R.; Vieno, M.; MacKenzie, I.A.; Armstrong, B.G.; Butland, B.K.; Milojevic, A.; Chalabi, Z.; Atkinson, R.W.; Stevenson, D.S.; et al. Spatiotemporal evaluation of EMEP4UK-WRF v4.3 atmospheric chemistry transport simulations of health-related metrics for NO2, O3, PM10, and PM2. 5 for 2001–2010. Geosci. Model Dev. 2017, 10, 1767–1787. [Google Scholar] [CrossRef] [Green Version]
  51. Brookes, D.M.; Stedman, J.R.; Grice, S.E.; Kent, A.J.; Walker, H.L.; Cooke, S.L.; Vincent, K.J.; Lingard, J.J.N.; Bush, T.J.; Abbott, J. UK Air Quality Modelling under the Air Quality Directive (2008/50/EC) for 2010 Covering the Following Air Quality Pollutants: SO2, NOx, NO2, PM10, PM2.5, Lead, Benzene, CO, and Ozone. Report for the Department for Environment, Food and Rural Affairs (Defra), Welsh Government, Scottish Government and the Department of the Environment in Northern Ireland. AEA report. AEAT/ENV/R/3215 Issue 1. 2011. Available online: http://uk-air.defra.gov.uk/reports/cat09/1204301513_AQD2010mapsrep_master_v0.pdf (accessed on 6 July 2020).
  52. Air Quality Expert Group (AQEG). Mitigation of United Kingdom PM2.5 Concentrations 2013. Available online: https://uk-air.defra.gov.uk/assets/documents/reports/cat11/1508060903_DEF-PB14161_Mitigation_of_UK_PM25.pdf (accessed on 6 July 2020).
  53. European Space Agency. Copernicus Sentinel-5 Precursor Mission. Available online: https://sentinel.esa.int/web/sentinel/missions/sentinel-5p. (accessed on 15 October 2020).
  54. European Space Agency. Copernicus Sentinel-4 Mission. Available online: https://sentinel.esa.int/web/sentinel/missions/sentinel-4. (accessed on 16 October 2020).
  55. European Space Agency. Copernicus Sentinel-5 Mission. Available online: https://sentinel.esa.int/web/sentinel/missions/sentinel-5. (accessed on 16 October 2020).
Figure 1. Spatial distribution of 581 PM10 (black star) and 183 PM2.5 (turquoise dots) monitors across Great Britain during the study period.
Figure 1. Spatial distribution of 581 PM10 (black star) and 183 PM2.5 (turquoise dots) monitors across Great Britain during the study period.
Remotesensing 12 03803 g001
Figure 2. Satellite-AOD 0.47 µm values are represented in Figure 2a (mean of all Terra- and Aqua-MODIS passing times). The missing values in grey were imputed through a combination of multiple modelled-AOD wavelengths and sub-day times, represented in Figure 2b by modelled-AOD 0.47 µm at 12:00. The Stage-2 output is shown in Figure 2c, illustrating the full coverage of AOD 0.47 µm based on the combination of measurements and estimations across Great Britain. The maps correspond to values measured or reconstructed on 6 July 2018.
Figure 2. Satellite-AOD 0.47 µm values are represented in Figure 2a (mean of all Terra- and Aqua-MODIS passing times). The missing values in grey were imputed through a combination of multiple modelled-AOD wavelengths and sub-day times, represented in Figure 2b by modelled-AOD 0.47 µm at 12:00. The Stage-2 output is shown in Figure 2c, illustrating the full coverage of AOD 0.47 µm based on the combination of measurements and estimations across Great Britain. The maps correspond to values measured or reconstructed on 6 July 2018.
Remotesensing 12 03803 g002
Figure 3. Stage-4 predicted PM2.5 concentrations across Great Britain (Top) and London (Bottom) for 2008, 2013, and 2018 aggregated by annual means. All plots were built under the same colour scale.
Figure 3. Stage-4 predicted PM2.5 concentrations across Great Britain (Top) and London (Bottom) for 2008, 2013, and 2018 aggregated by annual means. All plots were built under the same colour scale.
Remotesensing 12 03803 g003
Figure 4. Stage-4 daily series of observed and predicted PM2.5 values across Great Britain for 2008, 2013, and 2018: (1) Traffic monitor—located in Dundee (Scotland); (2) Urban monitor—located in London (England); and (3) Rural monitor—located in Sutton Bridge (England).
Figure 4. Stage-4 daily series of observed and predicted PM2.5 values across Great Britain for 2008, 2013, and 2018: (1) Traffic monitor—located in Dundee (Scotland); (2) Urban monitor—located in London (England); and (3) Rural monitor—located in Sutton Bridge (England).
Remotesensing 12 03803 g004
Figure 5. Stage-4 day-specific PM2.5 estimations across Great Britain (Top) and London (Bottom).
Figure 5. Stage-4 day-specific PM2.5 estimations across Great Britain (Top) and London (Bottom).
Remotesensing 12 03803 g005
Table 1. Predicted-PM2.5 concentrations obtained from Stage-1 RF models were regressed against measured-PM2.5 concentrations in a linear regression model. The performance was evaluated using two CV methods (OOB and 10-Fold) together with RMSE (a measure of the model error, µg/m3), intercept (µg/m3), and slope (µg/m3).
Table 1. Predicted-PM2.5 concentrations obtained from Stage-1 RF models were regressed against measured-PM2.5 concentrations in a linear regression model. The performance was evaluated using two CV methods (OOB and 10-Fold) together with RMSE (a measure of the model error, µg/m3), intercept (µg/m3), and slope (µg/m3).
Stage-1
OOB-CV10-Fold CV
R2RMSEInter.SlopeR2RMSEInter.Slope
20080.9181.196−0.4171.0330.7074.9540.8030.886
20090.9211.110−0.3901.0300.7913.9960.4100.937
20100.9191.122−0.4011.0290.8433.4960.0430.983
20110.9491.124−0.2661.0190.9023.439−0.0870.997
20120.9421.058−0.2741.0210.8893.218−0.0350.986
20130.9291.087−0.3681.0280.8473.5840.2180.972
20140.9440.963−0.2671.0220.8913.003−0.0070.995
20150.9330.865−0.2651.0260.8712.662−0.0030.983
20160.9350.896−0.2511.0250.8852.654−0.0500.996
20170.9390.828−0.1961.0220.8952.4300.0100.985
20180.9280.791−0.2481.0280.8862.235−0.0190.993
Mean0.9321.003−0.3041.0260.8553.2430.1170.974
Table 2. Predicted-AOD 0.47 µm and 0.55 µm obtained from Stage-2 RF models were regressed against measured Satellite-AOD 0.47 µm and 0.55 µm in a linear regression model. The performance was evaluated using OOB CV together with RMSE (µg/m3), intercept (µg/m3), and slope (µg/m3).
Table 2. Predicted-AOD 0.47 µm and 0.55 µm obtained from Stage-2 RF models were regressed against measured Satellite-AOD 0.47 µm and 0.55 µm in a linear regression model. The performance was evaluated using OOB CV together with RMSE (µg/m3), intercept (µg/m3), and slope (µg/m3).
Stage-2 OOB CV
Predicted-AOD 0.47 µmPredicted-AOD 0.55 µm
R2RMSEInter.SlopeR2RMSEInter.Slope
20080.9770.010−0.0011.0090.9770.007−0.0011.009
20090.9760.010−0.0011.0100.9760.007−0.0011.010
20100.9680.009−0.0011.0130.9680.007−0.0011.013
20110.9880.010−0.0011.0050.9880.0070.0001.005
20120.9800.010−0.0011.0080.9810.007−0.0011.008
20130.9840.010−0.0011.0070.9840.007−0.0011.006
20140.9700.009−0.0011.0120.9700.007−0.0011.012
20150.9720.009−0.0011.0110.9730.007−0.0011.011
20160.9750.009−0.0011.0100.9750.007−0.0011.010
20170.9630.009−0.0011.0150.9630.007−0.0011.014
20180.9690.010−0.0011.0130.9690.007−0.0011.013
Mean0.9780.010−0.0011.0090.9780.007−0.0011.009
Table 3. Relative importance (%) of the predictors in Stage-3 for the first, middle, and last years.
Table 3. Relative importance (%) of the predictors in Stage-3 for the first, middle, and last years.
Stage-3 Predictors200820132018
EMEP4UK PM2.532.4132.8336.74
Spatially-lagged hotspot-PM2.5 regional2.552.496.77
Wind direction6.357.335.34
Spatially-lagged background-PM2.5 regional1.146.064.93
Day of the year3.224.333.73
Spatially-lagged hotspot-PM2.5 local3.971.653.66
Precipitation6.632.423.25
BLH 0h2.282.912.79
Spatially-lagged background-PM2.5 local0.943.072.75
Month1.762.722.68
2m Air temperature2.602.932.65
Wind speed3.083.752.53
Sea-level pressure3.092.602.49
Relative humidity1.771.561.81
Nearest non-traffic monitor distance3.102.301.78
Note: The top 15 predictors in the RF importance ranking order is determined by the year 2018. The importance is measured by the amount of error reduced due to the splits of a given predictor over all trees used in the RF ensemble.
Table 4. Predicted-PM2.5 concentrations obtained from Stage-3 RF models were regressed against Stage-1 measured/predicted-PM2.5 concentrations in a linear regression model. The CV–R2 (how well the model described the PM2.5 variability in new locations) described in three different patterns (overall, spatial, and temporal), RMSE (µg/m3), intercept (µg/m3), and slope (µg/m3).
Table 4. Predicted-PM2.5 concentrations obtained from Stage-3 RF models were regressed against Stage-1 measured/predicted-PM2.5 concentrations in a linear regression model. The CV–R2 (how well the model described the PM2.5 variability in new locations) described in three different patterns (overall, spatial, and temporal), RMSE (µg/m3), intercept (µg/m3), and slope (µg/m3).
Stage-3
OverallSpatialTemporal
R2RMSEInter.SlopeR2RMSEInter.SlopeR2RMSEInter.Slope
20080.7044.547−1.2511.0640.4862.698−0.7491.0260.7603.6770.0001.074
20090.7424.247−1.1041.0420.6802.255−0.2030.9820.7623.5930.0001.055
20100.7094.330−1.4241.0750.6272.3420.1370.9720.7383.6280.0001.102
20110.8214.421−0.8981.0290.7332.280−0.5091.0030.8433.7560.0001.035
20120.7864.354−0.7491.0270.6612.5270.0730.9660.8233.5520.0001.043
20130.7644.305−1.0931.0470.6372.616−0.5651.0130.7913.6040.0001.061
20140.7844.140−1.0441.0510.6322.292−0.1450.9830.8153.4780.0001.062
20150.7363.792−1.1941.0720.5792.139−0.0260.9690.7763.1270.0001.095
20160.7813.702−0.9801.0500.7251.964−0.5321.0100.7963.1490.0001.061
20170.8163.343−0.9331.0410.7461.720−0.4060.9940.8342.8600.0001.055
20180.7903.275−1.0301.0460.7261.776−0.7451.0150.8072.7750.0001.056
Mean 0.7674.042−1.0641.0490.6582.237−0.3340.9940.7953.3820.0001.063
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Schneider, R.; Vicedo-Cabrera, A.M.; Sera, F.; Masselot, P.; Stafoggia, M.; de Hoogh, K.; Kloog, I.; Reis, S.; Vieno, M.; Gasparrini, A. A Satellite-Based Spatio-Temporal Machine Learning Model to Reconstruct Daily PM2.5 Concentrations across Great Britain. Remote Sens. 2020, 12, 3803. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12223803

AMA Style

Schneider R, Vicedo-Cabrera AM, Sera F, Masselot P, Stafoggia M, de Hoogh K, Kloog I, Reis S, Vieno M, Gasparrini A. A Satellite-Based Spatio-Temporal Machine Learning Model to Reconstruct Daily PM2.5 Concentrations across Great Britain. Remote Sensing. 2020; 12(22):3803. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12223803

Chicago/Turabian Style

Schneider, Rochelle, Ana M. Vicedo-Cabrera, Francesco Sera, Pierre Masselot, Massimo Stafoggia, Kees de Hoogh, Itai Kloog, Stefan Reis, Massimo Vieno, and Antonio Gasparrini. 2020. "A Satellite-Based Spatio-Temporal Machine Learning Model to Reconstruct Daily PM2.5 Concentrations across Great Britain" Remote Sensing 12, no. 22: 3803. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12223803

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop