Production of high quality interpolation maps of heavy metals is important for risk assessment of environmental pollution. In this paper, the spatial correlation characteristics information obtained from Moran’s I analysis was used to supplement the traditional geostatistics. According to Moran’s I analysis, four characteristics distances were obtained and used as the active lag distance to calculate the semivariance. Validation of the optimality of semivariance demonstrated that using the two distances where the Moran’s I and the standardized Moran’s I, Z(I) reached a maximum as the active lag distance can improve the fitting accuracy of semivariance. Then, spatial interpolation was produced based on the two distances and their nested model. The comparative analysis of estimation accuracy and the measured and predicted pollution status showed that the method combining geostatistics with Moran’s I analysis was better than traditional geostatistics. Thus, Moran’s I analysis is a useful complement for geostatistics to improve the spatial interpolation accuracy of heavy metals.

This study explored the spatial pattern of heavy metals in Beijing agricultural soils using Moran’s I statistic of spatial autocorrelation. The global Moran’s I result showed that the spatial dependence of Cr, Ni, Zn, and Hg changed with different spatial weight matrixes, and they had significant and positive global spatial correlations based on distance weight. The spatial dependence of the four metals was scale-dependent on distance, but these scale effects existed within a threshold distance of 13 km, 32 km, 50 km, and 29 km, respectively for Cr, Ni, Zn, and Hg. The maximal spatial positive correlation range was 57 km, 70 km, 57 km, and 55 km for Cr, Ni, Zn, and Hg, respectively and these were not affected by sampling density. Local spatial autocorrelation analysis detected the locations of spatial clusters and spatial outliers and revealed that the pollution of these four metals occurred in significant High-high spatial clusters, Low-high, or even High-low spatial outliers. Thus, three major areas were identified and should be receiving more attention: the first was the northeast region of Beijing, where Cr, Zn, Ni, and Hg had significant increases. The second was the southeast region of Beijing where wastewater irrigation had strongly changed the content of metals, particularly of Cr and Zn, in soils. The third area was the urban fringe around city, where Hg showed a significant increase.

Multifaceted regime shifts of Earth's surface are ongoing dramatically and – in turn – considerably alter the global carbon budget, energy balance and biogeochemical cycles. Sustainably managing terrestrial ecosystems necessitates a deeper comprehension of the diverse and dynamic nature of multicomponent information within these environments. However, comprehensive records of global-scale fractional vegetation and soil information that encompass these structural and functional complexities remain limited. Here, we provide a globally comprehensive record of monthly vegetation and soil fractions during the period 2001–2022 using a spatiotemporally adaptive spectral mixture analysis framework. This product is designed to continuously represent Earth's terrestrial surface as a percentage of five physically meaningful vegetation and soil endmembers, including photosynthetic vegetation (PV), nonphotosynthetic vegetation (NPV), bare soil (BS), ice or snow (IS) and dark surface (DA), with high accuracy and low uncertainty compared to the previous vegetation index and vegetation continuous-field product as well as traditional fully constrained linear spectral mixture models. We also adopt nonparametric seasonal Mann–Kendall tested fractional dynamics to identify shifts based on interactive changes in these fractions. Our results – superior to previous portrayals of the greening planet – not only report a +9.35 × 10⁵ km² change in photosynthetic vegetation, but also explore decreases in nonphotosynthetic vegetation (−2.19 × 10⁵ km²), bare soil (−5.14 × 10⁵ km²) and dark surfaces (−2.27 × 10⁵ km²). In addition, interactive changes in these fractions yield multifaceted regime shifts with important implications, such as a simultaneous increase in PV and NPV in central and southwestern China during afforestation activities, an increase in PV in cropland of China and India due to intensive agricultural development, a decrease in PV and an increase in BS in tropical zones resulting from deforestation. These advantages emphasize that our dataset provides locally relevant information on multifaceted regime shifts at the required scale, enabling scalable modeling and effective governance of future terrestrial ecosystems. The data about five fractional surface vegetation and soil components are available in the Science Data Bank (https://doi.org/10.57760/sciencedb.13287, Sun and Sun, 2023).

Understanding the effects of intensive agricultural land use activities on water resources is essential for natural resource management and environmental improvement. In this paper, multi-scale nested watersheds were delineated and the relationships between two representative water quality indexes and agricultural land use intensity were assessed and quantified for the year 2000 using multi-scale regression analysis. The results show that the log-transformed nitrate-nitrogen (NO₃-N) index exhibited a relationship with chemical fertilizer input intensity and several natural factors, including soil loss, rainfall and sunlight at the first order watershed scale, while permanganate index (COD_Mn) had a positive relationship with another two input intensities of pesticides and agricultural plastic mulch and organic manure at the fifth order watershed scale. The first order watershed and the fifth order watershed were considered as the watershed adaptive response units for NO₃-N and COD_Mn, respectively. The adjustment of agricultural input and its intensity may be carried out inside the individual watershed adaptive response unit. The multiple linear regression model demonstrated the cause-and-effect relationship between agricultural land use intensity and stream water quality at multiple scales, which is an important factor for the maintenance of stream water quality.