Soil contamination by potentially toxic elements (PTEs) is one of the greatest threats to environmental degradation. Knowing where PTEs accumulated in soil can mitigate their adverse effects on plants, animals, and human health. We evaluated the potential of using long-term remote sensing images that reveal the bare soils, to detect and map PTEs in agricultural fields. In this study, 360 soil samples were collected at the superficial layer (0–20 cm) in a 2574 km² agricultural area located in São Paulo State, Brazil. We tested the Soil Synthetic Image (SYSI) using Landsat TM/ETM/ETM+, Landsat OLI, and Sentinel 2 images. The three products have different spectral, temporal, and spatial resolutions. The time series multispectral images were used to reveal areas with bare soil and their spectra were used as predictors of soil chromium, iron, nickel, and zinc contents. We observed a strong linear relationship (−0.26 > r > −0.62) between the selected PTEs and the near infrared (NIR) and shortwave infrared (SWIR) bands of Sentinel (ensemble of 4 years of data), Landsat TM (35 years data), and Landsat OLI (4 years data). The clearest discrimination of soil PTEs was obtained from SYSI using a long term Landsat 5 collection over 35 years. Satellite data could efficiently detect the contents of PTEs in soils due to their relation with soil attributes and parent materials. Therefore, distinct satellite sensors could map the PTEs on tropics and assist in understanding their spatial dynamics and environmental effects.

Pollution of riverine ecosystems through the multidimensional impact of human footprints around the world poses a serious challenge. Research studies that communicate potential repercussions of landscape structure metrics on snowmelt riverine water quality particularly, in climatically fragile Himalayan watersheds are very scarce. Though, worldwide, grasping the influence of land-use practices on water quality (WQ) has received renewed attention yet, the relevance of spatial scale linked to landscape pattern is still elusive due to its heterogenic nature across diverse geomorphic regions. In this work, therefore, we tried to capture the insights on landscape-aquascape interface by juxtapositioning the impacts of landscape structure pattern on snowmelt stream WQ of the whole Jhelum River Basin (JRB) under three varying spatial scales viz., watershed scale, riparian corridor (1000 m wide) and reach buffer (500 m wide). The percentage of landscape pattern composition and configuration metrics in the JRB were computed in GIS utilizing Landsat-8 OLI/TIRS satellite image having 30 m resolution. To better explicate the influence of land-use metrics on riverine WQ with space and time, we used Redundancy analysis (RDA) and multilinear regression (MLR) modeling. MLR selected land-use structure metrics revealed the varied response of WQ parameters to multi-scale factors except for total faecal coliform bacteria (TC) which showed perpetual presence. The reach-scale explained slightly better (76%) variations in WQ than riparian (75%) and watershed (70%) scales. Likewise, across seasonal scale, autumn (75%), winter (83%), and summer (77%) captured the most WQ variation at catchment, riparian, and reach scales respectively. We observed impairing WQ linkages with agriculture, built-up and barren rocky areas across watersheds, besides, pastures in riparian buffer areas, and fragmentation of landscape patches at the reach scale. Due to little appearance of spatial scale differences, a multi scale perspective landscape planning is emphasized to ensure future sustainability of Kashmir Himalayan water resources.

The Mugan Plain is the most productive area in the Republic of Azerbaijan, but a previous study confirmed trace metal and metalloid (TM&M) contamination with Cr, Ni and Pb, and the potential ecological risk of As was estimated. However, no industrial activity was previously reported in this area; thus, a source apportionment model using positive matrix factorization (PMF) was employed to identify pollution sources, and a human health risk assessment was conducted to evaluate noncarcinogenic and carcinogenic risks. Surface soil samples were collected from 349 sites, and six major elements (Si, Ca, Cl, P, S and Sr) and 8 TM&Ms (As, Cd, Cr, Co, Cu, Ni, Pb and Zn) were analyzed by X-ray fluorescence and employed for further apportionment and risk assessment. As a result, the PMF model showed 7 factors, assigned to natural activity (12.9%), dry riverbed (13.6%), surface accumulation (3.1%), desalinization activity (3.2%), residential activity (12.3%), fossil fuel combustion (35.5%) and agricultural activity (19.3%). The PMF model characterized certain areas with desalinization activity in the previous Soviet period and with surface accumulation of salt, and these findings were confirmed by additional field surveys and historical Landsat satellite images. The risk assessment results showed that there was no risk for the adults, while for children, there was a noncarcinogenic risk, but no carcinogenic risk. Dermal contact was estimated to be the primary pathway, and Ni and As were identified as the most problematic TM&Ms for noncarcinogenic and carcinogenic risks, respectively. According to the results, fossil fuel combustion associated with heating and vehicle transportation was estimated to be the main source of pollution, contributing 42.6% of the noncarcinogenic and 48.0% of the carcinogenic risks. These results can provide scientific guidance to understand and prevent the risk of TM&Ms on the Mugan Plain.

Ecosystem services are present everywhere, green vegetation coverage (or green areas) is one of the primary sources of ecosystem services considering urban areas sustainability and peoples urban life quality. Urban vegetation cover loss decreases the capacity of nature to provision ecosystem services; the loss of urban vegetation is also observed within the Amazon. This study aims at identifying urban vegetation loss and relate it to the provision of ecosystem services of reduction of air quality, reduction of air pollution, and climate regulation. Urban vegetation coverage loss was calculated using NDVI on LANDSAT 5 imagery over a 23-year period from 1986 to 2009. NDVI thresholds were arbitrarily selected, and complemented by in locus observation, to establish guidelines for quantitative (area) and qualitative (density) evolution of green cover, divided in six different categories, named as water, bare soil, poor vegetation, moderate vegetation, dense vegetation and very dense vegetation. Data on air pollution, noise pollution and temperature were outsourced from previous works. Measurement show a significant loss of very dense, dense and moderate vegetation coverage and an increase in poor vegetation and bare soil areas, in accordance to increase in air and noise pollution, and local temperature, and provides positive refashions between the loss of urban green coverage and decrease in ecosystem services.

Natural environments have been linked to decreased risk of cardiovascular disease (CVD) and respiratory disease (RSD) mortality. However, few cohort studies have looked at associations of natural environments with CVD or RSD hospitalization. The aim of this study was to evaluate these associations in a cohort of U.S. Medicare beneficiaries (∼63 million individuals). Our open cohort included all fee-for-service Medicare beneficiaries (2000–2016), aged ≥65, living in the contiguous U.S. We assessed zip code-level park cover based on the United States Geological Survey Protected Areas Database, average greenness (Normalized Difference Vegetation Index, NDVI), and percent blue space cover based on Landsat satellite images. Cox-equivalent Poisson models were used to estimate associations of the exposures with first CVD and RSD hospitalization in the full cohort and among those living in urban zip codes (≥1000 persons/mile²). NDVI was weakly negatively correlated with percent park cover (Spearman ρ = −0.23) and not correlated with percent blue space (Spearman ρ = 0.00). After adjustment for potential confounders, percent park cover was not associated with CVD or RSD hospitalization in the full or urban population. An IQR (0.27) increase in NDVI was negatively associated with CVD (HR: 0.97, 95%CI: 0.96, 0.97), but not with RSD hospitalization (HR: 0.99, 95%CI: 0.98, 1.00). In urban zip codes, an IQR increase in NDVI was positively associated with RSD hospitalization (HR: 1.02, 95%CI: 1.00, 1.03). In stratified analyses, percent park cover was negatively associated with CVD and RSD hospitalization for Medicaid eligible individuals and individuals living in low socioeconomic status neighborhoods in the urban population. We observed no associations of percent blue space cover with CVD or RSD hospitalization. This study suggests that natural environments may benefit cardiorespiratory health; however, benefits may be limited to certain contexts and certain health outcomes.

Dissolved oxygen (DO) is an effective indicator for water pollution. However, since DO is a non-optically active parameter and has little impact on the spectrum captured by satellite sensors, research on estimating DO by remote sensing at multiple spatiotemporal scales is limited. In this study, the support vector regression (SVR) models were developed and validated using the remote sensing reflectance derived from both Landsat and Moderate Resolution Imaging Spectroradiometer (MODIS) data and synchronous DO measurements (N = 188) and water temperature of Lake Huron and three other inland waterbodies (N = 282) covering latitude between 22–45 °N. Using the developed models, spatial distributions of the annual and monthly DO variability since 1984 and the annual monthly DO variability since 2000 in Lake Huron were reconstructed for the first time. The impacts of five climate factors on long-term DO trends were analyzed. Results showed that the developed SVR-based models had good robustness and generalization (average R² = 0.91, root mean square percentage error = 2.65%, mean absolute percentage error = 4.21%), and performed better than random forest and multiple linear regression. The monthly DO estimates by Landsat and MODIS data were highly consistent (average R² = 0.88). From 1984 to 2019, the oxygen loss in Lake Huron was 6.56%. Air temperature, incident shortwave radiation flux density, and precipitation were the main climate factors affecting annual DO of Lake Huron. This study demonstrated that using SVR-based models, Landsat and MODIS data could be used for long-term DO retrieval at multiple spatial and temporal scales. As data-driven models, combining spectrum and water temperature as well as extending the training set to cover more DO conditions could effectively improve model robustness and generalization.

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.

Intra-urban pollution monitoring requires fine particulate (PM₂.₅) concentration mapping at ultrahigh-resolution (dozens to hundreds of meters). However, current PM₂.₅ concentration estimation, which is mainly based on aerosol optical depth (AOD) and meteorological data, usually had a low spatial resolution (kilometers) and severe spatial missing problem, cannot be applied to intra-urban pollution monitoring. To solve these problems, top-of-atmosphere reflectance (TOAR), which contains both the information about land and atmosphere and has high resolution and large spatial coverage, may be efficiently used for PM₂.₅ estimation. This study aims to systematically evaluate the feasibility of retrieving ultrahigh-resolution PM₂.₅ concentration at a large scale (national level) from TOAR. Firstly, we make a detailed discussion about several important but unsolved theoretic problems on TOAR-based PM₂.₅ retrieval, including the band selection, scale effect, cloud impact, and mapping quality evaluation. Secondly, four types and eight retrieval models are compared in terms of quantitative accuracy, mapping quality, model generalization, and model efficiency, with the pros and cons of each type summarized. Deep neural network (DNN) model shows the highest retrieval accuracy, and linear models were the best in efficiency and generalization. As a compromise, ensemble learning shows the best overall performance. Thirdly, using the highly accurate DNN model (cross-validated R² equals 0.93) and through combining Landsat 8 and Sentinel 2 images, a 90 m and ∼4-day resolution PM₂.₅ product was generated. The retrieved maps were used for analyzing the fine-scale interannual pollution change inner the city and the pollution variations during novel coronavirus disease 2019 (COVID-19). Results of this study proves that ultrahigh resolution can bring new findings of intra-urban pollution change, which cannot be observed at previous coarse resolution. Lastly, some suggestions for future ultrahigh-resolution PM₂.₅ mapping research were given.

Due to rapid urbanization in China, lead (Pb) continues to accumulate in urban topsoil, resulting in soil degradation and increased public exposure. Mapping Pb concentrations in urban topsoil is therefore vital for the evaluation and control of this exposure risk. This study developed spatial models to map Pb concentrations in urban topsoil using proximal and remote sensing data. Proximal sensing reflectance spectra (350–2500 nm) of soils were pre-processed and used to calculate the principal components as landscape factors to represent the soil properties. Other landscape factors, including vegetation and land-use factors, were extracted from time-sequential Landsat images. Two hybrid statistical approaches, regression kriging (RK) and geographically weighted regression (GWR), were adopted to establish prediction models using the landscape factors. The results indicated that the use of landscape factors derived from combined remote and proximal sensing data improved the prediction of Pb concentrations compared with useing these data individually. GWR obtained better results than RK for predicting soil Pb concentration. Thus, joint proximal and remote sensing provides timely, easily accessible, and suitable data for extracting landscape factors.

This study applies remote sensing technology to assess and examine the spatial and temporal Brightness Temperature (BT) profile in the city of Tel-Aviv, Israel over the last 30 years using Landsat imagery. The location of warmest and coldest zones are constant over the studied period. Distinct diurnal and temporal BT behavior divide the city into four different segments. As an example of future application, we applied mixed regression models with daily random slopes to correlate Landsat BT data with monitored air temperature (Tair) measurements using 14 images for 1989–2014. Our preliminary results show a good model performance with R² = 0.81. Furthermore, based on the model's results, we analyzed the spatial profile of Tair within the study domain for representative days.