Future phosphorus (P) shortages could seriously affect terrestrial productivity and food security. We investigated the changes in topsoil available P (AP) and total P (TP) in China\'s forests, grasslands, paddy fields, and upland croplands during the 1980s-2010s based on substantial repeated soil P measurements (63,220 samples in the 1980s, 2000s, and 2010s) and machine learning techniques. Between the 1980s and 2010s, total soil AP stock increased with a small but significant rate of 0.13 kg P ha-1  year-1 , but total soil TP stock declined substantially (4.5 kg P ha-1  year-1 ) in the four ecosystems. We quantified the P budgets of soil-plant systems by harmonizing P fluxes from various sources for this period. Matching trends of soil contents over the decades with P budgets and fluxes, we found that the P-surplus in cultivated soils (especially in upland croplands) might be overestimated due to the great soil TP pool compared to fertilization and the substantial soil P losses through plant uptake and water erosion that offset the P additions. Our findings of P-deficit in China raise the alarm on the sustainability of future biomass production (especially in forests), highlight the urgency of P recycling in croplands, and emphasize the critical role of country-level basic data in guiding sound policies to tackle the global P crises.

Soil organic carbon (SOC), total nitrogen (TN) and total phosphorus (TP) are important indicators reflecting soil quality, and they can be used to effectively evaluate the effect of soil remediation. Many studies have evaluated the content of SOC, TN and TP in different ecosystems. However, after constructing protected forests for ecological restoration in the ecologically fragile coastal zone, the spatial distribution and influencing mechanism of SOC, TN and TP content is still uncertain. In this study, the spatial heterogeneity and influencing factors of SOC, TN and TP in surface (0-20 cm) soil were analyzed by traditional analysis and geostatistics. A total of 39 soil samples were collected under the coastal zone protected forest types including Quercus acutissima Carruth (QAC), Pinus thunbergii Parl (PTP), mixed PTP and QAC (QP) and Castanea mollissima BL (CMB) in the coastal zone protected forests in northern China. The results show that SOC, TN and TP content were defined as moderate variation, and they also show significant changes under different protected forest types (P < 0.05). The semivariance results indicate that SOC, TN and TP all exhibited strong spatial dependence class, with Range of 224 m, 229 m and 282 m respectively, which were more than the sampling scale of 200 m. The spatial prediction results showed that SOC, TN and TP content all appear in large areas of extremely low value in CMB, and its cross validation results showed that using vegetation and terrain factors as covariates in the spatial prediction of SOC, TN and TP can improve the prediction accuracy. The results of correlation analysis showed that the influencing factor for SOC and TN, and TP were NDVI and topographical changes, respectively. In general, vegetation and terrain factors as auxiliary factors can improved the accuracy of soil C-N-P spatial distribution prediction after afforestation in coastal zone.

Ecological restoration is widespread in the karst region, southwest China, but the impacts of different restoration strategies on soil fertility indices have rarely been compared. Here soil nutrients and microbial communities were measured 16 years after agricultural abandonment in a karst area, southwest China. Three restoration strategies were included, i.e., i) restoration with an economic tree species Toona sinensis (TS), ii) restoration with Guimu-1 hybrid elephant grass (GG), iii) restoration with a combination of Zenia insignis and Guimu-1 hybrid elephant grass (ZG). Cropland under maize-soybean rotation (CR) was used as reference. Soil organic carbon level was more than doubled in TS, and that in GG and ZG was elevated by about 50% relative to CR. Soil total nitrogen concentration in GG was not significantly different from CR, but that in TS and ZG was increased by 93% and 55% relative to CR. Similar to nitrogen, soil total phosphorus concentration in GG was not changed relative to CR, but that in TS and ZG were significantly increased. Microbial biomass carbon and nitrogen concentrations were significantly increased in TS and GG by 124% and 82%, respectively, compared to CR, but those in ZG and CR were similar. The abundance of total PLFAs (phospholipid fatty acids) was significantly increased by 55-69% following agricultural abandonment, and there was no significant difference among the three restoration strategies. The patterns of the other microbial groups and the ratio of fungal to bacterial (F:B) PLFAs were largely similar to that of total PLFAs. Soil organic carbon was identified as the primary factor affecting the abundance of soil microbial communities. Our findings suggest that the three restoration strategies, particularly TS are efficient in improving soil fertility.

Soil total nitrogen (STN) and total phosphorus (STP) are important indicators of soil nutrients and the important indexes of soil fertility and soil quality evaluation. Using geographic information system (GIS) and geostatistics, the spatial heterogeneity distribution of STN and STP in the Yaoxiang watershed in a hilly area of northern China was studied. The results showed that: (1) The STN and STP contents showed a declining trend with the increase in soil depth; the variation coefficients (Cv) of STN and STP in the 0- to 10-cm soil layer (42.25% and 14.77%, respectively) were higher than in the 10- to 30-cm soil layer (28.77% and 11.60%, respectively). Moreover, the Cv of STN was higher than that of STP. (2) The maximum C0/(C0 + C1) of STN and STP in the soil layers was less than 25%, this indicated that a strong spatial distribution autocorrelation existed for STN and STP; and the STP showed higher intensity and more stable variation than the STN. (3) From the correlation analysis, we concluded that the topographic indexes such as elevation and slope direction all influenced the spatial distribution of STN and STP (correlation coefficients were 0.688 and 0.518, respectively). (4) The overall distribution of STN and STP in the Yaoxiang watershed decreased from the northwest to the southeast. This variation trend was similar to the watershed DEM trend and was significantly influenced by vegetation and topographic factors. These results revealed the spatial heterogeneity distribution of STN and STP, and addressed the influences of forest vegetation coverage, elevation, and other topographic factors on the spatial distribution of STN and STP at the watershed scale.