Assessing the space-time trends and detecting the sources of heavy metal accumulation in soils have important consequences in the prevention and treatment of soil heavy metal pollution. In this ...study, we collected soil samples in the eastern part of the Qingshan district, Wuhan city, Hubei Province, China, during the period 2010–2014. The Cd, Cu, Pb and Zn concentrations in soils exhibited a significant accumulation during 2010–2014. The spatiotemporal Kriging technique, based on a quantitative characterization of soil heavy metal concentration variations in terms of non-separable variogram models, was employed to estimate the spatiotemporal soil heavy metal distribution in the study region. Our findings showed that the Cd, Cu, and Zn concentrations have an obvious incremental tendency from the southwestern to the central part of the study region. However, the Pb concentrations exhibited an obvious tendency from the northern part to the central part of the region. Then, spatial overlay analysis was used to obtain absolute and relative concentration increments of adjacent 1- or 5-year periods during 2010–2014. The spatial distribution of soil heavy metal concentration increments showed that the larger increments occurred in the center of the study region. Lastly, the principal component analysis combined with the multiple linear regression method were employed to quantify the source apportionment of the soil heavy metal concentration increments in the region. Our results led to the conclusion that the sources of soil heavy metal concentration increments should be ascribed to industry, agriculture and traffic. In particular, 82.5% of soil heavy metal concentration increment during 2010–2014 was ascribed to industrial/agricultural activities sources.
Using STK and SOA to obtain the spatial distribution of heavy metal concentration increments in soils.
Using PCA-MLR to quantify the source apportionment of soil heavy metal concentration increments.
Display omitted
•Spatiotemporal Kriging was used to obtain the spatiotemporal distribution of soil heavy metals.•Spatial overlay was used to obtain the spatial distribution of soil heavy metals concentration increment.•The PCA-MLR method was used to quantify the source apportionment of soil heavy metal concentration increments.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK, ZRSKP
With rapid economic development, industrialization and urbanization, the ambient air PM2.5 has become a major pollutant linked to respiratory, heart and lung diseases. In China, PM2.5 pollution ...constitutes an extreme environmental and social problem of widespread public concern. In this work we estimate ground-level PM2.5 from satellite-derived aerosol optical depth (AOD), topography data, meteorological data, and pollutant emission using an integrative technique. In particular, Geographically Weighted Regression (GWR) analysis was combined with Bayesian Maximum Entropy (BME) theory to assess the spatiotemporal characteristics of PM2.5 exposure in a large region of China and generate informative PM2.5 space-time predictions (estimates). It was found that, due to its integrative character, the combined BME-GWR method offers certain improvements in the space-time prediction of PM2.5 concentrations over China compared to previous techniques. The combined BME-GWR technique generated realistic maps of space-time PM2.5 distribution, and its performance was superior to that of seven previous studies of satellite-derived PM2.5 concentrations in China in terms of prediction accuracy. The purely spatial GWR model can only be used at a fixed time, whereas the integrative BME-GWR approach accounts for cross space-time dependencies and can predict PM2.5 concentrations in the composite space-time domain. The 10-fold results of BME-GWR modeling (R2 = 0.883, RMSE = 11.39 μg/m3) demonstrated a high level of space-time PM2.5 prediction (estimation) accuracy over China, revealing a definite trend of severe PM2.5 levels from the northern coast toward inland China (Nov 2015–Feb 2016). Future work should focus on the addition of higher resolution AOD data, developing better satellite-based prediction models, and related air pollutants for space-time PM2.5 prediction purposes.
•A combined BME-GWR method was used to study space-time PM2.5 pollution in China.•Cross space-time PM2.5 correlations and environmental factors were jointly processed.•BME-GWR improved PM2.5 prediction accuracy (R2 = 0.883) compared to previous methods.•A severe PM2.5 trend (northern coast to inland China, Nov 2015–Feb 2016) was detected.•BME-GWR is useful for exposure assessment and health risk management purposes.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK, ZRSKP
Improving the understanding and characterization of spatial soil heavy metal distribution is becoming an important component of risk assessment and environmental policy. In this work, 213 soil ...samples collected from Daye (Hubei Province, China) were used as the empirical dataset. First, maps of soil heavy metal distributions, including Cd, Co, Cr, Cu, Mn, Ni, Pb and Zn, were obtained using the ordinary Kriging method. Then, the pollution index (PI) and integrated pollution index (IPI) were calculated based on the ordinary Kriging maps to obtain a comprehensive quantitative pollution characterization of the eight heavy metals in the Daye soil. The results showed that 46.1%, 32.1%, and 0.5% of the soil in the study region are moderately, highly and extremely polluted, respectively. Finally, the one- and two-point stochastic site indicators of IPI were used to assess quantitatively the uncertainties and risks associated with soil heavy metal distributions in the polluted regions. These results showed that the IPI values exceeding a specified threshold increased almost linearly with increasing threshold value, whereas the relative area of excess pollution decreased steadily with increasing threshold. Among the site pairs considered in the study region, about 70% and 26% of them simultaneously experienced moderate and high pollution risk, respectively.
Display omitted
•Site stochastic indicators (SSI) are useful in soil remediation management.•SSI offer a comprehensive way to assess the significance of pollution in soils.•One-point SSI offer global averages of pollution using environmental thresholds.•Two-points SSI offer information regarding the spatial variation of pollution.•SSI quantify and characterize the pollution spatial morphology and evolution.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
The accurate assessment of spatiotemporal rainfall variability is a crucial and challenging task in many hydrological applications, mainly due to the lack of a sufficient number of rain gauges. The ...purpose of the present study is to investigate the spatiotemporal variations of annual and monthly rainfall over Fujian province in China by combining the Bayesian maximum entropy (BME) method and satellite rainfall estimates. Specifically, based on annual and monthly rainfall data at 20 meteorological stations from 2000 to 2012, (1) the BME method with Tropical Rainfall Measuring Mission (TRMM) estimates considered as soft data, (2) ordinary kriging (OK) and (3) cokriging (CK) were employed to model the spatiotemporal variations of rainfall in Fujian province. Subsequently, the performance of these methods was evaluated using cross-validation statistics. The results demonstrated that BME with TRMM as soft data (BME-TRMM) performed better than the other two methods, generating rainfall maps that represented the local rainfall disparities in a more realistic manner. Of the three interpolation (mapping) methods, the mean absolute error (MAE) and root mean square error (RMSE) values of the BME-TRMM method were the smallest. In conclusion, the BME-TRMM method improved spatiotemporal rainfall modeling and mapping by integrating hard data and soft information. Lastly, the study identified new opportunities concerning the application of TRMM rainfall estimates.
Population health attributes (such as disease incidence and prevalence) are often estimated using sentinel hospital records, which are subject to multiple sources of uncertainty. When applied to ...these health attributes, commonly used biased estimation techniques can lead to false conclusions and ineffective disease intervention and control. Although some estimators can account for measurement error (in the form of white noise, usually after de-trending), most mainstream health statistics techniques cannot generate unbiased and minimum error variance estimates when the available data are biased.
A new technique, called the Biased Sample Hospital-based Area Disease Estimation (B-SHADE), is introduced that generates space-time population disease estimates using biased hospital records. The effectiveness of the technique is empirically evaluated in terms of hospital records of disease incidence (for hand-foot-mouth disease and fever syndrome cases) in Shanghai (China) during a two-year period. The B-SHADE technique uses a weighted summation of sentinel hospital records to derive unbiased and minimum error variance estimates of area incidence. The calculation of these weights is the outcome of a process that combines: the available space-time information; a rigorous assessment of both, the horizontal relationships between hospital records and the vertical links between each hospital's records and the overall disease situation in the region. In this way, the representativeness of the sentinel hospital records was improved, the possible biases of these records were corrected, and the generated area incidence estimates were best linear unbiased estimates (BLUE). Using the same hospital records, the performance of the B-SHADE technique was compared against two mainstream estimators.
The B-SHADE technique involves a hospital network-based model that blends the optimal estimation features of the Block Kriging method and the sample bias correction efficiency of the ratio estimator method. In this way, B-SHADE can overcome the limitations of both methods: Block Kriging's inadequacy concerning the correction of sample bias and spatial clustering; and the ratio estimator's limitation as regards error minimization. The generality of the B-SHADE technique is further demonstrated by the fact that it reduces to Block Kriging in the case of unbiased samples; to ratio estimator if there is no correlation between hospitals; and to simple statistic if the hospital records are neither biased nor space-time correlated. In addition to the theoretical advantages of the B-SHADE technique over the two other methods above, two real world case studies (hand-foot-mouth disease and fever syndrome cases) demonstrated its empirical superiority, as well.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
The accurate and informative space-time mapping of air pollutants is a crucial component of many human exposure studies. In the present work, space-time maps of daily distributions of PM
2.5
and NO
2
...concentrations were generated in the severely polluted northern China region using the Bayesian maximum entropy (BME) method. This method can incorporate hard PM
2.5
and NO
2
data (obtained at ground-level monitoring sites), and various kinds of soft (uncertain) data, including satellite data processed in terms of machine learning techniques, meteorological variables, and geographical predictors. The BME maps of space-time PM
2.5
and NO
2
concentrations over northern China generated during the winter season (when severe haze episodes occur frequently) were realistic and informative. As regards their numerical accuracy, for the space-time PM
2.5
estimates, the tenfold cross-validation
R
2
and the RMSE were, respectively, 0.86 and 14.37 μg/m
3
; for the space-time NO
2
estimates, the
R
2
and RMSE values were, respectively, 0.85 and 6.93 μg/m
3
. Lastly, it was shown that the BME method performed better than the mainstream spatiotemporal ordinary kriging technique in terms of the higher
R
2
values of both the predicted PM
2.5
and NO
2
concentration maps.
Typhoid and paratyphoid fever are endemic in Hongta District and their prevalence, at 113 per 100,000 individuals, remains the highest in China. However, the exact sources of the disease and its main ...epidemiological characteristics have not yet been clearly identified.
Numbers of typhoid and paratyphoid cases per day during the period 2006 to 2010 were obtained from the Chinese Center of Disease Control (CDC). A number of suspected disease determinants (or their proxies), were considered for use in spatiotemporal analysis: these included locations of discharge canals and food markets, as well as socio-economic and environmental factors. Results showed that disease prevalence was spatially clustered with clusters decreasing with increasing distance from markets and discharge canals. More than half of the spatial variance could be explained by a combination of economic conditions and availability of health facilities. Temporal prevalence fluctuations were positively associated with the monthly precipitation series. Polluted hospital and residential wastewater was being discharged into rainwater canals. Salmonella bacteria were found in canal water, on farmland and on vegetables sold in markets.
DISEASE TRANSMISSION IN HONGTA DISTRICT IS DRIVEN PRINCIPALLY BY TWO SPATIOTEMPORALLY COUPLED CYCLES: one involving seasonal variations and the other the distribution of polluted farmland (where vegetables are grown and sold in markets). Disease transmission was exacerbated by the fact that rainwater canals were being used for disposal of polluted waste from hospitals and residential areas. Social factors and their interactions also played a significant role in disease transmission.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Improving the spatiotemporal coverage of remote sensing (RS) products, such as sea surface chlorophyll concentration (SSCC), can offer a better understanding of the spatiotemporal SSCC distribution ...for ocean management purposes. In the first part of this work, 834 in-situ SSCC measurements of the SeaBASS-NASA (National Aeronautics and Space Administration) during 2002–2016 served as the empirical dataset. A moving window with ±3 days and ±0.5° centered at each of the in-situ SSCC measurements established a search neighborhood for Moderate Resolution Imaging Spectroradiometer Level 2 (MODIS L2) SSCC and MODIS L2 sea surface temperature (SST) data, and the matched SSCC and SST data were used for building a linear SSCC-SST relationship. The unmatched SST was introduced to the linear model for generating soft SSCC data with uniform distributions. The inherent spatiotemporal dependency of the SSCC distribution was then represented by the Bayesian maximum entropy (BME) method, which incorporated the soft SSCC data as auxiliary variable for SSCC estimation and mapping purposes. The results showed that a 75.3% accuracy improvement of remote SSCC retrieval in terms of R2 can be achieved by BME-based method compared to the original MODIS L2 product. Subsequently, the BME-based method was applied to obtain daily SSCC dataset in Chesapeake Bay (USA) during the period 2010–2019. It was found that the SSCC distribution exhibited a decreasing spatial trend from the upper bay to the outer bay, whereas decreasing and increasing temporal trends were detected during the periods 2011–2014 and 2016–2019, respectively. The generalized Cauchy process was used to quantitatively describe the autocorrelation SSCC function in the Chesapeake Bay. The results showed that the outer bay exhibited the strongest long-range dependence among the four sub-regions, whereas the middle bay exhibited the weakest long-range dependence. Finally, one-point and two-point stochastic site indicators (SSIs) were employed to explore the spatiotemporal SSCC characteristics in Chesapeake Bay. The one-point SSI results showed that nearly 100% of the upper, middle and the lower bay areas experienced a high SSCC level (>5 mg/m3) during the entire study period. The area with SSCC >5 mg/m3 in the outer bay increased a lot during the winter season, but the area with SSCC >10 or 20 mg/m3 decreased significantly in the upper, middle and lower bay. Simultaneously, the SSCC dispersion in these areas was rather small during the winter season. On the other hand, the two-point SSI results showed that although the SSCC levels differ among the four sub-regions, but the SSCC connectivity structures between pairs of points also displayed some similarities in terms of their spatiotemporal dependency. In conclusion, the proposed BME-based method was shown to be a promising remote SSCC mapping technique that exhibited a powerful ability to improve both accuracy and coverage of RS products. The SSIs can be also used to explore the spatiotemporal characteristics of a variety of natural attributes in waters.
Display omitted
•Introducing SST can help improve the accuracy and coverage of remote SSCC by BME.•Modeling at a local scale is necessary to depict the SSCC-SST relationship.•Decreasing trend of SSCC was found from upper to outer bay of Chesapeake Bay (CPB).•SSIs are useful for exploring the spatiotemporal characteristics of SSCC.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
Soil heavy metal pollution can be a serious threat to human health and the environment. The accurate mapping of the spatial distribution of soil heavy metal pollutant concentrations enables the ...detection of high pollution areas and facilitates pollution source apportionment and control. To make full use of auxiliary soil properties information and show that they can improve mapping, a synthesis of the Bayesian Maximum Entropy (BME) theory and the Geographically Weighted Regression (GWR) model is proposed and implemented in the study of the Shanghai City soils (China). The results showed that, compared to traditional techniques, the proposed BME-GWR synthesis has certain important advantages: (a) it integrates heavy metal measurements and auxiliary information on a sound theoretical basis, and (b) it performs better in terms of both prediction accuracy and implementation flexibility (including the assimilation of multiple data sources). Based on the heavy metal concentration maps generated by BME-GWR, we found that the As, Cr and Pb concentration levels are high in the eastern part of Shanghai, whereas high Cd concentration levels were observed in the northwestern part of the city. Organic carbon and pH were significantly correlated with most of the heavy metals in Shanghai soils. We concluded that Cd pollution is mainly the result of agricultural activities, and that the Cr pollution is attributed to natural sources, whereas Pb and As have compound pollution sources. Future studies should investigate the implementation of BME-GWR in the case of space-time heavy metal mapping and its ability to integrate human activity information and soil category variables.
Display omitted
•BME-GWR performs better in terms of prediction accuracy and performance flexibility.•Organic carbon and pH were significantly correlated with most of the heavy metals.•Cr was higher in eastern region, which is attributed to natural sources.•Cd was higher in northwestern region, which is mainly due to agricultural activities.•Pb and As were higher in eastern region having compound pollution sources.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
Remote sensing reflectance (Rrs) values measured by satellite sensors involve large amounts of uncertainty leading to non-negligible noise in remote Chlorophyll-a (Chl-a) concentration estimation. ...This work distinguished between two main stages in the case of estimating distributions of Chl-a within the Gulf of St. Lawrence (Canada). At the model building stage, the retrieval algorithm used both in-situ Chl-a measurements and the corresponding Moderate Resolution Imaging Spectroradiometer (MODIS) L2-level data estimated Rrs at 412, 443, 469, 488, 531, 547, 555, 645, 667, 678 nm at a 1 km spatial resolution during 2004–2013. Through the training and validation of various models and Rrs combinations of the considered eight techniques (including support vector regression, artificial neural networks, gradient boosting machine, random forests, standard CI-OC3M, multiple linear regression, generalized addictive regression, principal component regression), the support vector regression (SVR) technique was shown to have the best performance in Chl-a concentration estimation using Rrs at 412, 443, 488, 531 and 678 nm. The accuracy indicators for both the training (850) and the validation (213) datasets were found to be very good to excellent (e.g., the R2 value varied between 0.7058 and 0.9068). At the space-time estimation stage, this work took a step forward by using the Bayesian maximum entropy (BME) theory to further process the SVR estimated Chl-a concentrations by incorporating the inherent spatiotemporal dependency of physical Chl-a distribution. A 56% improvement was achieved in the reduction of the mean uncertainty of the validation data decreased considerably (from 1.2222 to 0.5322 mg/m3). Then, this novel BME/SVR framework was employed to estimate the daily Chl-a concentrations in the Gulf of St. Lawrence during Jan 1-Dec 31 of 2017 (1 km spatial resolution). The results showed that the daily mean Chl-a concentration varied from 1.6630 to 3.3431 mg/m3, and that the daily mean Chl-a uncertainty reduction of the composite BME/SVR vs. the SVR estimation had a maximum reduction value of 1.0082 and an average reduction value of 0.6173 mg/m3. The monthly spatial Chl-a distribution covariances showed that the highest Chl-a concentration variability occurred during November and that the spatiotemporal Chl-a concentration pattern changed a lot during the period August to November. In conclusion, the proposed BME/SVR was shown to be a promising remote Chl-a retrieval approach that exhibited a significant ability in reducing the non-negligible uncertainty and improving the accuracy of remote sensing Chl-a concentration estimates.
•BME and SVR (ML) were combined to retrieve Chl-a in optically complex waters.•SVR shows great ability to retrieve Chl-a using satellite remote sensing data.•Introducing BME can significantly reduce the uncertainty of Chl-a estimation.•The distribution pattern of Chl-a at GSL varies a lot during the summer season.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP