The estimation of pollution fields, especially in densely populated areas, is an important application in the field of environmental science due to the significant effects of air pollution on public health. In this paper, we investigate the spatial distribution of three air pollutants in Tehran’s atmosphere: carbon monoxide (CO), nitrogen dioxide (NO2), and atmospheric particulate matters less than 10 μm in diameter (PM10μm). To do this, we use four geostatistical interpolation methods: Ordinary Kriging, Universal Kriging, Simple Kriging, and Ordinary Cokriging with Gaussian semivariogram, to estimate the spatial distribution surface for three mentioned air pollutants in Tehran’s atmosphere. The data were collected from 21 air quality monitoring stations located in different districts of Tehran during 2012 and 2013 for 00UTC. Finally, we evaluate the Kriging estimated surfaces using three statistical validation indexes: mean absolute error (MAE), root mean square error (RMSE) that can be divided into systematic and unsystematic errors (RMSES, RMSEU), and D-Willmot. Estimated standard errors surface or uncertainty band of each estimated pollutant surface was also developed. The results indicated that using two auxiliary variables that have significant correlation with CO, the ordinary Cokriginga scheme for CO consistently outperforms all interpolation methods for estimating this pollutant and simple Kriging is the best model for estimation of NO2 and PM10. According to optimal model, the highest concentrations of PM10 are observed in the marginal areas of Tehran while the highest concentrations of NO2 and CO are observed in the central and northern district of Tehran.

Finding an appropriate satellite image as simultaneous as possible with the sampling time campaigns is challenging. Fusion can be considered as a method of integrating images and obtaining more pixels with higher spatial, spectral and temporal resolutions. This paper investigated the impact of Landsat 8-OLI and Sentinel-2A data fusion on prediction of several toxic elements at a mine waste dump. The 30 m spatial resolution Landsat 8-OLI bands were fused with the 10 m Sentinel-2A bands using various fusion techniques namely hue-saturation-value (HSV), Brovey, principal component analysis (PCA), Gram-Schmidt (GS), wavelet, and area-to-point regression kriging (ATPRK). ATPRK was the best method preserving both spectral and spatial features of Landsat 8-OLI and Sentinel-2A after fusion. Furthermore, the partial least squares regression (PLSR) model developed on genetic algorithm (GA)-selected laboratory visible-near infrared-shortwave infrared (VNIR–SWIR) spectra yielded more accurate prediction results compared to the PLSR model calibrated on the entire spectra. It was hence, applied to both individual sensors and their ATPRK-fused image. In case of the individual sensors, except for As, Sentinel-2A provided more robust prediction models than Landsat 8-OLI. However, the best performances were obtained using the fused images, highlighting the potential of data fusion to enhance the toxic elements’ prediction models.

The traditional incremental lifetime cancer risk (ILCR) model of urban soil polycyclic aromatic hydrocarbon (PAH) health risk assessment has a large spatial scale and commonly calculates relevant statistics by regarding the whole area as a geographic unit but fails to consider the high heterogeneity of the PAH distribution and differences in population susceptibility and density in an area. Therefore, the risk assessment spatial performance is insufficient and does not reflect the characteristics of cities, which are centered on human activities and serve the needs of humans, thus making it difficult to effectively support PAH prevention and treatment measures in cities. Here, the random forest model combined with the kriging residual model (RFerr-K) is used to estimate high-precision PAH distributions, separately considering the exposure characteristics of children and adults with different susceptibilities, and kindergarten point-of-interest (POI) and population density index (PDI) data were used to estimate the distributions of the kindergarten children and adults in the study area. Through the refined expression of these three dimensions, a new spatially explicit model of the incremental lifetime cancer-causing population distribution (MapPILCR) was constructed, and the risk threshold range delineation method was proposed to accurately identify regional risk levels. The results showed that the RFerr-K model significantly improves the accuracy of PAH prediction. The susceptibility index (SI) of children is 45% higher than that of adults, and POI and PDI data can be used effectively in population distribution estimation. The MapPILCR model provides a useful method for the spatially explicit assessment of the cancer risk of urban populations to inspire urban pollution grid management.

Different prediction models have important effects on the accuracy of spatial distribution simulations of heavy metals in soil. This study proposes a model (RFOK) combining a random forest (RF) with ordinary kriging (OK), multi-source environmental data such as terrain elements, site environmental elements, and remote sensing data were incorporated to predict the spatial distribution of heavy arsenic (As) in soil of a certain large arsenic slag site. The predictions results of RFOK were compared with those obtained using the RF, OK, inverse distance weighted (IDW), and stepwise regression (STEPREG) models for assessment of prediction accuracy. The results showed that arsenic pollution was widely distributed and the center of the site, including arsenic slag stacking area and production area were seriously polluted. The overall spatial distribution of arsenic pollution simulated by the five models was similar, but the IDW, RF, OK, and STEPREG showed less spatial variation of soil pollution, while RFOK simulation can better express the characteristics of details in change. The cross-validation results showed that RFOK had the lowest root-mean-square error (RMSE), mean absolute error (MAE), and mean relative error (MRE) relative to the other four models, followed by RF, OK, IDW, and STEPREG. The RMSE, MAE and MRE of RFOK decreased by 62.2%, 64.3% and 68.7%, respectively, relative to the RF model with the second highest accuracy. Compared with the traditional spatial distribution prediction model, the RFOK model proposed in this study has excellent spatial distribution prediction ability for soil heavy metal pollution with large spatial variation characteristics, which can fully explain the nonlinear relationship between pollutant content and its environmental impact elements.

Nitrogen is one of the most significant pollutants in the Yangtze River estuary (YRE), China. Reliable estimation of nitrogen concentration in the water is crucial for assessment of the water quality of the estuary. Because ocean fronts exist in the YRE, which divide water masses into different regions, it is necessary to account for the heterogeneity of the water surface when predicting nitrogen concentrations. A new geostatistical method, called spatiotemporal point mean of surface with non-homogeneity (ST-PMSN), is proposed to model the non-stationary spatiotemporal random process of nitrogen concentrations between 2004 and 2013 in the YRE. The method considers the spatiotemporal correlation of surface water nitrogen and uses information from both sides of a boundary for heterogeneous water masses. Comparing with several other interpolating methods, including spatial ordinary kriging (OK), stratified ordinary kriging (SOK), point mean of surface with non-homogeneity (P-MSN), spatiotemporal ordinary kriging (STK), and stratified spatiotemporal ordinary kriging (SSTK), the cross-validation results show that ST-PMSN has the highest accuracy, followed by SSTK, STK, P-MSN, SOK, and OK in descending order. ST-PMSN is therefore demonstrated to be effective in estimating the nitrogen pollutant concentrations in a stratified estuary. According to interpolated nitrogen concentrations in the YRE, water quality has generally deteriorated—with fluctuations—from 2004 to 2013. The average annual reduction in area of water quality of Grades I and II from 2004 to 2013 was 1.10%. At the same time, the average annual increase in area of water quality of Grades III and IV was 0.89% and that of Grade V was 0.21%. The results of this study provide a new and more accurate interpolating method for assessing the pollutant concentration in the marine and offers guidance for more precise classification of water quality in the YRE.

Heavy metal soil contamination is associated with potential toxicity to humans or ecotoxicity. Scholars have increasingly used a combination of geographical information science (GIS) with geostatistical and multivariate statistical analysis techniques to examine the spatial distribution of heavy metals in soils at a regional scale. A review of such studies showed that most soil sampling programs were based on grid patterns and composite sampling methodologies. Many programs intended to characterize various soil types and land use types. The most often used sampling depth intervals were 0–0.10 m, or 0–0.20 m, below surface; and the sampling densities used ranged from 0.0004 to 6.1 samples per km², with a median of 0.4 samples per km². The most widely used spatial interpolators were inverse distance weighted interpolation and ordinary kriging; and the most often used multivariate statistical analysis techniques were principal component analysis and cluster analysis. The review also identified several determining and correlating factors in heavy metal distribution in soils, including soil type, soil pH, soil organic matter, land use type, Fe, Al, and heavy metal concentrations. The major natural and anthropogenic sources of heavy metals were found to derive from lithogenic origin, roadway and transportation, atmospheric deposition, wastewater and runoff from industrial and mining facilities, fertilizer application, livestock manure, and sewage sludge. This review argues that the full potential of integrated GIS and multivariate statistical analysis for assessing heavy metal distribution in soils on a regional scale has not yet been fully realized. It is proposed that future research be conducted to map multivariate results in GIS to pinpoint specific anthropogenic sources, to analyze temporal trends in addition to spatial patterns, to optimize modeling parameters, and to expand the use of different multivariate analysis tools beyond principal component analysis (PCA) and cluster analysis (CA).

Exposure of the general population to airborne metals remains poorly estimated despite the potential health risks. Passive moss biomonitoring can proxy air quality at fine resolution over large areas, mainly in rural areas. We adapted the technique to urban areas to develop fine concentration maps for several metals for Constances cohort's participants. We sampled Grimmia pulvinata in 77 and 51 cemeteries within ∼50 km of Paris and Lyon city centers, respectively. We developed land-use regression models for 14 metals including cadmium, lead, and antimony; potential predictors included the amount of urban, agricultural, forest, and water around cemeteries, population density, altitude, and distance to major roads. We used both kriging with external drift and land use regression followed by residual kriging when necessary to derive concentration maps (500 × 500 m) for each metal and region. Both approaches led to similar results. The most frequent predictors were the amount of urban, agricultural, or forest areas. Depending on the metal, the models explained part of the spatial variability, from 6% for vanadium in Lyon to 84% for antimony in Paris, but mostly between 20% and 60%, with better results for metals emitted by human activities. Moss biomonitoring in cemeteries proves efficient for obtaining airborne metal exposures in urban areas for the most common metals.

Nitrogen dioxide (NO₂) is an important air pollutant that causes direct harms to the environment and human health. Ground NO₂ mapping with high spatiotemporal resolution is critical for fine-scale air pollution and environmental health research. We thus developed a spatiotemporal regression kriging model to map daily high-resolution (3-km) ground NO₂ concentrations in China using the Tropospheric Monitoring Instrument (TROPOMI) satellite retrievals and geographical covariates. This model combined geographically and temporally weighted regression with spatiotemporal kriging and achieved robust prediction performance with sample-based and site-based cross-validation R² values of 0.84 and 0.79. The annual mean and standard deviation of ground NO₂ concentrations from June 1, 2018 to May 31, 2019 were predicted to be 15.05 ± 7.82 μg/m³, with that in 0.6% of China’s area (10% of the population) exceeding the annual air quality standard (40 μg/m³). The ground NO₂ concentrations during the coronavirus disease (COVID-19) period (January and February in 2020) was 14% lower than that during the same period in 2019 and the mean population exposure to ground NO₂ was reduced by 25%. This study was the first to use TROPOMI retrievals to map fine-scale daily ground NO₂ concentrations across all of China. This was also an early application to use the satellite-estimated ground NO₂ data to quantify the impact of the COVID-19 pandemic on the air pollution and population exposures. These newly satellite-derived ground NO₂ data with high spatiotemporal resolution have value in advancing environmental and health research in China.

Previous ecological risk assessments were mainly concentration-oriented rather than source-oriented. Moreover, land use is usually related to source emissions but was rarely used to improve the source apportionment accuracy. In this study, the land-use effects of heavy metals (HMs) in surface (0–20 cm) and subsurface (20–40 cm) soils were first explored using ANOVA in a suburb of Changzhou City, China; next, based on robust absolute principal component scores-robust geographically weighted regression (RAPCS/RGWR), this study proposed RAPCS/RGWR with land-use type (RAPCS/RGWR-LUT) and compared its source apportionment accuracy with those of basic RAPCS/RGWR and commonly-used absolute principal component scores/multiple linear regression (APCS/MLR); then, the source-oriented ecological risks were apportioned based on RAPCS/RGWR-LUT and Hakanson potential ecological risk index method; finally, this study proposed robust residual kriging with land-use type (RRK) for spatially predicting the source-oriented ecological risks, and compared its spatial prediction accuracy with those of robust ordinary kriging (ROK) and traditionally-used ordinary kriging (OK). Results showed that: (i) by incorporating land-use effects, RAPCS/RGWR-LUT obtained higher source apportionment accuracy than RAPCS/RGWR and APCS/MLR; (ii) the two most important external input sources of the ecological risks were 'atmospheric deposition' (PERIₛᵤᵣfₐcₑ = 47.11 and PERIₛᵤbₛᵤᵣfₐcₑ = 35.27) and 'agronomic measure' (PERIₛᵤᵣfₐcₑ = 28.93 and PERIₛᵤbₛᵤᵣfₐcₑ = 20.37); (iii) the biggest ecological risk factor was soil Cd (ERₛᵤᵣfₐcₑ = 57.14 and ERₛᵤbₛᵤᵣfₐcₑ = 47.62), which was mainly contributed by 'atmospheric deposition' (ERₛᵤᵣfₐcₑ=33.14 and ERₛᵤbₛᵤᵣfₐcₑ=25.71); (iv) RRK obtained higher spatial prediction accuracy than ROK and OK; (v) the high-risk areas derived from 'atmospheric deposition' were mainly located in the southwest of the study area, and the high-risk areas derived from 'agronomic measure' were scattered in the agricultural land in the north and south of the study area. The above information provided effective spatial decision support for reducing the source-oriented input of the ecological risks of soil HMs in a large-scale area.

The interrelationships among long-term ambient air pollution exposure, emotional distress and cognitive decline in older adulthood remain unclear. Long-term exposure may impact cognitive performance and subsequently impact emotional health. Conversely, exposure may initially be associated with emotional distress followed by declines in cognitive performance. Here we tested the inter-relationship between global cognitive ability, emotional distress, and exposure to PM₂.₅ (particulate matter with aerodynamic diameter <2.5 μm) and NO₂ (nitrogen dioxide) in 6118 older women (aged 70.6 ± 3.8 years) from the Women’s Health Initiative Memory Study. Annual exposure to PM₂.₅ (interquartile range [IQR] = 3.37 μg/m³) and NO₂ (IQR = 9.00 ppb) was estimated at the participant’s residence using regionalized national universal kriging models and averaged over the 3-year period before the baseline assessment. Using structural equation mediation models, a latent factor capturing emotional distress was constructed using item-level data from the 6-item Center for Epidemiological Studies Depression Scale and the Short Form Health Survey Emotional Well-Being scale at baseline and one-year follow-up. Trajectories of global cognitive performance, assessed by the Modified-Mini Mental State Examination (3MS) annually up to 12 years, were estimated. All effects reported were adjusted for important confounders. Increases in PM₂.₅ (β = -0.144 per IQR; 95% CI = −0.261; −0.028) and NO₂ (β = −0.157 per IQR; 95% CI = −0.291; −0.022) were associated with lower initial 3MS performance. Lower 3MS performance was associated with increased emotional distress (β = −0.008; 95% CI = −0.015; −0.002) over the subsequent year. Significant indirect effect of both exposures on increases in emotional distress mediated by exposure effects on worse global cognitive performance were present. No statistically significant indirect associations were found between exposures and 3MS trajectories putatively mediated by baseline emotional distress. Our study findings support cognitive aging processes as a mediator of the association between PM₂.₅ and NO₂ exposure and emotional distress in later-life.