Current understanding of phosphorus (P) cycling in soils can be enhanced by integrating previously discrete findings concerning P speciation, exchange kinetics, and the underlying biological and geochemical processes. Here, we combine sequential extraction with P K-edge X-ray absorption spectroscopy and isotopic methods (³³P and ¹⁸O in phosphate) to characterize P cycling on a climatic gradient in Hawaii. We link P pools to P species and estimate the turnover times for commonly considered P pools. Dissolved P turned over in seconds, resin-extractable P in minutes, NaOH-extractable inorganic P in weeks to months, and HCl-extractable P in years to millennia. Furthermore, we show that in arid-zone soils, some primary mineral P remains even after 150 ky of soil development, whereas in humid-zone soils of the same age, all P in all pools has been biologically cycled. The integrative information we provide makes possible a more dynamic, process-oriented conceptual model of P cycling in soils.

Soil represents the largest phosphorus (P) stock in terrestrial ecosystems. Determining the amount of soil P is a critical first step in identifying sites where ecosystem functioning is potentially limited by soil P availability. However, global patterns and predictors of soil total P concentration remain poorly understood. To address this knowledge gap, we constructed a database of total P concentration of 5275 globally distributed (semi-)natural soils from 761 published studies. We quantified the relative importance of 13 soil-forming variables in predicting soil total P concentration and then made further predictions at the global scale using a random forest approach. Soil total P concentration varied significantly among parent material types, soil orders, biomes, and continents and ranged widely from 1.4 to 9630.0 (median 430.0 and mean 570.0) mg kg⁻¹ across the globe. About two-thirds (65 %) of the global variation was accounted for by the 13 variables that we selected, among which soil organic carbon concentration, parent material, mean annual temperature, and soil sand content were the most important ones. While predicted soil total P concentrations increased significantly with latitude, they varied largely among regions with similar latitudes due to regional differences in parent material, topography, and/or climate conditions. Soil P stocks (excluding Antarctica) were estimated to be 26.8 ± 3.1 (mean ± standard deviation) Pg and 62.2 ± 8.9 Pg (1 Pg = 1 × 10¹⁵ g) in the topsoil (0–30 cm) and subsoil (30–100 cm), respectively. Our global map of soil total P concentration as well as the underlying drivers of soil total P concentration can be used to constraint Earth system models that represent the P cycle and to inform quantification of global soil P availability. Raw datasets and global maps generated in this study are available at https://doi.org/10.6084/m9.figshare.14583375 (He et al., 2021).

Context: Agroecology combines agronomic and ecological concepts. It relies on the enhancement of biodiversity and related ecosystem services to support agricultural production. It is dependent on biological interactions for the design and management of agricultural systems in agricultural landscapes.Objectives: We review the role of landscape ecology to understand and promote biodiversity, pest regulation and crop pollination for the designing of “agroecology landscapes”. We illustrate the use of landscape ecological methods for supporting agroforestry systems as an example of agroecological development, and we propose pathways to implement agroecology at landscape scale.Methods: The state of the art of how landscape ecology contributes to agroecology development is summarized based on a literature review.Results: Agroecology requires thinking beyond the field scale to consider the positioning, quality and connectivity of fields and semi-natural habitats at larger spatial scales. The spatial and temporal organisation of semi-natural elements and the crop mosaic interact. Understanding this interaction is the pre-requisite for promoting patterns and mechanisms that foster biodiversity and ecosystem service provision. Promoting agroecological practices beyond individual farm borders can be rooted in a bottom-up approach from agroecological lighthouse farms to farm networks to amplify agroecology adoption at the landscape scale.Conclusions: Achieving agricultural landscapes composed of fields and farms following agroecological management requires understanding of biodiversity patterns, biological interactions and mechanisms that determine and boost ecosystem functioning to improve services at landscape scale, involving farmers in a bottom-up and context-specific approach.

Quantification of turnover of inorganic soil phosphorus (P) pools is essential to improve our understanding of P cycling in soil–plant systems and improve representations of the P cycle in land surface models. Turnover can be quantified using mean residence time (MRT); however, to date there is little information on MRT of P in soil P pools. We introduce an approach to quantify MRT of P in sequentially extracted inorganic soil P pools using data from isotope exchange kinetic experiments. Our analyses of 53 soil samples from the literature showed that MRT of labile P (resin- and bicarbonate-extractable P) was on the order of minutes to hours for most soils, MRT in NaOH-extractable P (NaOH-P) was in the range of days to months, and MRT in HCl-extractable P (HCl-P) was on the order of years to millennia. Multiple-regression models were able to capture 54 %–63 % of the variability in MRT among samples and showed that land use was the most important predictor of MRT of P in labile and NaOH pools. MRT of P in HCl-P was strongly dependent on pH, as high-pH soils tended to have longer MRTs. This was interpreted to be related to the composition of HCl-P. Under high pH, HCl-P contains mostly apatite, with a low solubility, whereas under low-pH conditions, HCl-P may contain more exchangeable P forms. These results suggest that current land surface models underestimate the dynamics of inorganic soil P pools and could be improved by reducing model MRTs of the labile and NaOH-P pools, considering soil-type-dependent MRTs rather than universal exchange rates and allowing for two-way exchange between HCl-P and the soil solution.