The black soil region in Northeast China is one of the major grain-producing areas of China. Soil erosion in the black soil region caused by natural factors and anthropic activities has attracted much attention, especially in a rolling hilly region. Compared with natural factors, the land use factor of cropland encompasses the most easily optimized measures. Jiutai County of Changchun City, located in the hilly areas of Northeast China, was taken as an example to calculate the soil erosion modulus using the Revised Universal Soil Loss Equation model. The overall soil erosion status of cultivated land in the study area was mainly slight and light, the proportion of cultivated land affected by extremely intensive and severe erosion was relatively small, and the average soil erosion modulus was 7.09 t·hm⁻²·a⁻¹ in 2019. In view of spatial distribution characteristics of soil erosion revealed by the spatial aggregation and hot spot analysis, the most serious soil erosion intensity was concentrated in the southeast and northeast sloping farmland over 8°. With the increase in elevation and topographic slope, the proportion of slight and light soil erosion gradually decreased, which was closely related to the increase in soil erodibility caused by the space–time migration of soil organic carbon caused by the interaction of hydraulic and tillage erosion in complex topographic areas. The Geographically Weighted Regression model was introduced to explore the driving factors and superposition mechanism of farmland soil erosion in the hilly region of Northeast China. Based on the relationship between soil erosion and landscape fragmentation, landscape fragmentation was an important driving force promoting soil erosion, sediment yield, and sediment transport. This paper is committed to providing a basis for accurately deploying regional soil and water conservation measures and formulating macro land management policies.

Accurate mapping of farmland soil organic carbon (SOC) provides valuable information for evaluating soil quality and guiding agricultural management. The integration of natural factors, agricultural activities, and landscape patterns may well fit the high spatial variation of SOC in low-relief farmlands. However, commonly used prediction methods are global models, ignoring the stratified heterogeneous relationship between SOC and environmental variables and failing to reveal the determinants of SOC in different subregions. Using 242 topsoil samples collected from Jianghan Plain, China, this study explored the stratified heterogeneous relationship between SOC and natural factors, agricultural activities, and landscape metrics, determined the dominant factors of SOC in each stratum, and predicted the spatial distribution of SOC using the Cubist model. Ordinary kriging, stepwise linear regression (SLR), and random forest (RF) were used as references. SLR and RF results showed that land use types, multiple cropping index, straw return, and percentage of water bodies are global dominant factors of SOC. Cubist results exhibited that the dominant factors of SOC vary in different cropping systems. Compared with the SOC of paddy fields, the SOC of irrigated land was more affected by irrigation-related factors. The effect of straw return on SOC was diverse under different cropping intensities. The Cubist model outperformed the other models in explaining SOC variation and SOC mapping (fitting R² = 0.370 and predicted R² = 0.474). These results highlight the importance of exploring the stratified heterogeneous relationship between SOC and covariates, and this knowledge provides a scientific basis for farmland zoning management. The Cubist model, integrating natural factors, agricultural activities, and landscape metrics, is effective in explaining SOC variation and mapping SOC in low-relief farmlands.

Accurately mapping soil organic carbon (SOC) is conducive to evaluating carbon storage and soil quality. However, the high spatial heterogeneity of SOC caused by river-related factors and agricultural management brings challenges to digital soil mapping in floodplain farmland. Moreover, current studies focus on the non-linear relationship between SOC and covariates, but ignore the effective range of environmental variables on SOC, which prevents the revelation of the SOC differentiation mechanism. Using the 375 samples collected from the Jiangchang Town near Han River, we aim to determine the main controlling factors of SOC, reveal the effective range of environmental variables, and obtain the spatial map of SOC by using the gradient boosting decision tree (GBDT) model and partial dependence plots. Linear regression was used as a reference. Results showed that GBDT outperformed linear regression. GBDT results show that the distance from the river was the most important SOC factor, confirming the importance of the Han River to the SOC pattern. The partial dependence plots indicate that all environmental variables have their effective ranges, and when their values are extremely high or low, they do not respond to changes in SOC. Specifically, the influential ranges of rivers, irrigation canals, and rural settlements on SOC were within 4000, 200, and 50 m, respectively. The peak SOC was obtained with high clay (≥31%), total nitrogen (≥1.18 g/kg), and total potassium contents (≥11.1 g/kg), but it remained steady when these covariates further increased. These results highlight the importance of revealing the effective range of environmental variables, which provides data support for understanding the spatial pattern of SOC in floodplain farmland, achieving carbon sequestration in farmland and precision agriculture. The GBDT with the partial dependence plot was effective in SOC fitting and mapping.

Exploring the local influencing factors and sources of soil arsenic (As) is crucial for reducing As pollution, protecting soil ecology, and ensuring human health. Based on geographically weighted regression (GWR), multiscale GWR (MGWR) considers the different influence ranges of explanatory variables and thus adopts an adaptative bandwidth. It is an effective model in many fields but has not been used in exploring local influencing factors and sources of As. Therefore, using 200 samples collected from the northeastern black soil zone of China, this study examined the effectiveness of MGWR, revealed the spatial non-stationary relationship between As and environmental variables, and determined the local impact factors and pollution sources of As. The results showed that 49% of the samples had arsenic content exceeding the background value, and these samples were mainly distributed in the central and southern parts of the region. MGWR outperformed GWR with the adaptative bandwidth, with a lower Moran’s I of residuals and a higher R² (0.559). The MGWR model revealed the spatially heterogeneous relationship between As and explanatory variables. Specifically, the road density and total nitrogen, clay, and silt contents were the primary or secondary influencing factors at most points. The distance from an industrial enterprise was the secondary influencing factor at only a few points. The main pollution sources of As were thus inferred as traffic and fertilizer, and industrial emissions were also included in the southern region. These findings highlight the importance of considering adaptative bandwidths for independent variables and demonstrate the effectiveness of MGWR in exploring local sources of soil pollutants.