The spatial characteristics, pollution sources, and risks of soil heavy metals were analyzed on Hainan Island. The results showed that the heavily polluted points accounted for 0.56%, and the number of mildly and above polluted points accounted for 15.27%, respectively, which were mainly distributed in the northern part of the study area. The principal component analysis–absolute principal component score–multiple linear regression (APCS-MLR) and the positive matrix factorization (PMF) revealed four sources of heavy metals: agricultural pollution sources for cadmium, (Cd), industrial and mining pollution sources for arsenic, (As), transportation pollution sources for zinc and lead (Zn and Pb), and natural pollution sources for chromium, nickel, and copper (Cr, Ni, and Cu). The human health risk assessment indicated that the average non-carcinogenic risk (HI) for both adults and children was within the safe threshold (<1), whereas Cr and Ni posed a carcinogenic risk (CR) to human health. In addition, the total non-carcinogenic risk (THI) indicated that heavy metals posed a potential non-carcinogenic risk to children, while the total carcinogenic risk (TCR) remained relatively high, mainly in the northern part of the study area. The results of the Monte Carlo simulation showed that the non-carcinogenic risk (HI) for all heavy metals was 1, indicating a potential health risk above the safe threshold. Meanwhile, nearly 100% and 99.94% of the TCR values exceeded 1 × 10⁻⁴ for children and adults, indicating that Cr and Ni are priority heavy metals for control. The research results provide the necessary scientific basis for the prevention and control of heavy metals in agricultural soils.

Terrestrial ecosystem carbon storage (TECS) could significantly affect the concentration of atmospheric CO₂, which is critical for climate change prediction. Along these lines, the Integrated Valuation of Ecosystem Services and Trade-offs model was employed to determine the TECS of Hainan Island (HN) from 2015 to 2050 accurately. Besides, the Future Land-use Simulation model combined with natural and anthropogenic factors was used to forecast the land-use types from 2025 to 2050 in HN by considering different Shared-socioeconomic pathway-Rrepresentative concentration pathway (SSP-RCP) scenarios. Finally, the geographical detector explored the influence mechanism concerning the TECS. Under the SSP1-RCP1.9 scenario, the TECS of HN will be gradually increased to 388.10 million tons in 2050, mainly due to the increase in forest areas and the fact that the majority of grassland in the western part of HN is being converted into forest. Under different SSP-RCP scenarios except for SSP1-RCP1.9, HN's TECS is expected to gradually decrease from 2015 to 2050, mainly due to the loss of grassland and forest in coastal low-altitude areas. From the single/pair factor perspective influenced mechanism concerning the TECS, the elevation (DEM) and DEM∩Slope were found to be the dominant single/pair factor under the SSP1-RCP1.9, SSP1-RCP2.6 and SSP2-RCP4.5 scenarios. The least distance to residential area (LDP) and LDP∩LDR (i.e. LDP and least distance to roads or railways) were found to be the dominant factors under the SSP3-RCP7.0, SSP4-RCP3.4, SSP4-RCP6.0, SSP5-RCP3.4 and SSP5-RCP8.5 scenarios. Besides, the pair factors provided a higher determinant power for TECS than a single factor. Given the results of the TECS and the influence mechanism concerning the TECS under different SSP-RCP scenarios, we suggest reasonably planning the transportation network and limiting the disorderly expansion of construction land.

Conventional fertilization practices can lead to many ecological problems, such as nutrient imbalance, soil acidity, and reduced soil fertility, in natural rubber plantations. To address these challenges, a field investigation was strategically carried out to substitute inorganic fertilizer with organic fertilizer, consisting of six treatments: no fertilization (CK), inorganic fertilizer (NPK), 25% replacement of inorganic through organic (25% manure (M)), 50% replacement of inorganic through organic (50% manure (M)), 75% replacement of inorganic through organic (75% manure (M)), and 100% organic fertilizer (100% manure). The soil physicochemical properties (soil organic carbon (SOC), total nitrogen (TN), mineral nitrogen (N), ammonium nitrogen (NH₄⁺-N), and nitrate nitrogen (NO₃⁻-N)), C:N, pH, and the carbon- and nitrogen-converting enzymes β-1,4-glucosidase (BG), N-acetylglucosaminidase (NAG) and L-leucine aminopeptidase (LAP) were all determined. The partial substitution of inorganic fertilizer with organic fertilizer (i.e., 75% M at surface soil layer) showed higher SOC (14.52 g·kg⁻¹), TN (1.06 g·kg⁻¹), N (20.07 mg·kg⁻¹), C:N (14.63), NH₄⁺-N (10.63 mg·kg⁻¹), and NO₃⁻-N (11.06 mg·kg⁻¹) than NPK and CK. This increase in physicochemical properties after partial replacement of inorganic with organic fertilizer resulted from higher carbon and nitrogen enzyme activities (BG (143.17·nmol·g⁻¹·h⁻¹), NAG (153.96 nmol·g⁻¹·h⁻¹), and LAP (153.48 nmol·g⁻¹·h⁻¹)) compared to NPK and CK. Further, the Pearson correlation and redundancy analysis (RDA) analyses confirmed a significant positive correlation between SOC, N, and soil enzymes. This study presents a new strategy for assessing the impact of partially replacing inorganic fertilizer with organic fertilizer in rubber plantations in tropical regions, mainly by modifying the soil nutrient composition.

Introduction: Soil organic carbon (SOC) accumulation changed with forest succession and hence impacted the SOC storage. However, the variation and underlying mechanisms about SOC during tropical forest succession are not fully understood.Methods: Soil samples at four depths (0–10 cm, 10–20 cm, 20–40 cm and 40–60 cm), litter, and roots of 0–10 cm and 10–20 cm were collected from three forest succession stages (plantation forest, secondary forest, and old– growth forest) in the Jianfengling (JFL) National Nature Reserve in Hainan Island, China. The SOC, soil enzyme activities, physiochemical properties, the biomass of litter and roots were analyzed.Results: Results showed that forest succession significantly increased SOC at 0–10 cm and 10–20 cm depth (from 23.00 g/kg to 33.70 g/kg and from 14.46 g/kg to 22.55 g/kg, respectively) but not at a deeper depth (20–60 cm). SOC content of the three forest succession stages decreased with increasing soil depth and bulk density (BD). With forest succession from plantation to secondary and old–growth forest, the soil pH at 0–10 cm and 10–20 cm depth decreased from 5.08 to 4.10 and from 5.52 to 4.64, respectively. Structural equation model (SEM) results showed that the SOC at depths of 0–20 cm increased with total root biomass but decreased with increasing soil pH value. The direct positive effect of soil TP on SOC was greater than the indirect negative effect of decomposition of SOC by soil acid phosphatase (AP).Discussion: To sum up, the study highlighted there was soil P– limited in tropical forests of JFL, and the increase in TP and total root biomass inputs were main factors favoring SOC sequestration during the tropical forest succession. In addition, soil acidification is of great importance for SOC accumulation in tropical forests for forest succession in the future. Therefore, forest succession improved SOC accumulation, TP and roots contributed to soil C sequestration.