Soil pH regulates the capacity of soils to store and supply nutrients, and thus contributes substantially to controlling productivity in terrestrial ecosystems 1 . However, soil pH is not an independent regulator of soil fertility-rather, it is ultimately controlled by environmental forcing. In particular, small changes in water balance cause a steep transition from alkaline to acid soils across natural climate gradients 2,3 . Although the processes governing this threshold in soil pH are well understood, the threshold has not been quantified at the global scale, where the influence of climate may be confounded by the effects of topography and mineralogy. Here we evaluate the global relationship between water balance and soil pH by extracting a spatially random sample (n = 20,000) from an extensive compilation of 60,291 soil pH measurements. We show that there is an abrupt transition from alkaline to acid soil pH that occurs at the point where mean annual precipitation begins to exceed mean annual potential evapotranspiration. We evaluate deviations from this global pattern, showing that they may result from seasonality, climate history, erosion and mineralogy. These results demonstrate that climate creates a nonlinear pattern in soil solution chemistry at the global scale; they also reveal conditions under which soils maintain pH out of equilibrium with modern climate.Climate controls many aspects of soil chemistry, affecting soil pH (ref. 4). Alkaline soils are known to be common in arid climates, while acid soils are known to be common in humid climates 1 . Surprisingly, however, the global-scale mechanisms governing this pattern remain broadly defined, and untested by direct observation. What are the dominant chemical equilibria that constrain soil pH? What aspect of climate defines the transition between alkaline and acid soils, and is the transition linear? The answers to these questions are fundamental to understanding soil development and surface geochemistry at the global scale. Furthermore, achieving this understanding may prove essential for representing soils in models of the terrestrial biosphere, given that soil pH controls many aspects of soil fertility 5,6 . Here we illustrate that simple geochemical and hydrological concepts can be used to build a mechanistic understanding of soil pH at the global scale.Interpretations of acid-titration experiments indicate that the soil pH is typically most strongly buffered by equilibrium with two secondary minerals: calcite (CaCO 3 ), or gibbsite (Al(OH) 3 ) 7,8 Fig. 1).Local studies of climate gradients have shown that the relative importance of these two buffers is determined by leaching, which removes Ca 2+ from the soil 2,3,8,9 . In climates where evaporative demand exceeds precipitation, leaching rates are low, and dissolved Ca 2+ accumulates as CaCO 3 -buffering soil pH near 8.2 (ref. 4). Conversely, in climates where precipitation exceeds evaporative demand, water leaches through the soil, removing Ca 2+ and allowing accumulation of relatively immobile...

Abstract. Stomatal conductance (g s ) affects the fluxes of carbon, energy and water between the vegetated land surface and the atmosphere. We test an implementation of an optimal stomatal conductance model within the Community Atmosphere Biosphere Land Exchange (CABLE) land surface model (LSM). In common with many LSMs, CABLE does not differentiate between g s model parameters in relation to plant functional type (PFT), but instead only in relation to photosynthetic pathway. We constrained the key model parameter "g 1 ", which represents plant water use strategy, by PFT, based on a global synthesis of stomatal behaviour. As proof of concept, we also demonstrate that the g 1 parameter can be estimated using two long-term average bioclimatic variables: (i) temperature and (ii) an indirect estimate of annual plant water availability. The new stomatal model, in conjunction with PFT parameterisations, resulted in a large reduction in annual fluxes of transpiration (∼ 30 % compared to the standard CABLE simulations) across evergreen needleleaf, tundra and C4 grass regions. Differences in other regions of the globe were typically small. Model performance against upscaled data products was not degraded, but did not noticeably reduce existing model-data biases. We identified assumptions relating to the coupling of the vegetation to the atmosphere and the parameterisation of the minimum stomatal conductance as areas requiring further investigation in both CABLE and potentially other LSMs.We conclude that optimisation theory can yield a simple and tractable approach to predicting stomatal conductance in LSMs.

Abstract. We implement a new stomatal conductance scheme, based on the optimality approach, within the Community Atmosphere Biosphere Land Exchange (CABLEv2.0.1) land surface model. Coupled land–atmosphere simulations are then performed using CABLEv2.0.1 within the Australian Community Climate and Earth Systems Simulator (ACCESSv1.3b) with prescribed sea surface temperatures. As in most land surface models, the default stomatal conductance scheme only accounts for differences in model parameters in relation to the photosynthetic pathway but not in relation to plant functional types. The new scheme allows model parameters to vary by plant functional type, based on a global synthesis of observations of stomatal conductance under different climate regimes over a wide range of species. We show that the new scheme reduces the latent heat flux from the land surface over the boreal forests during the Northern Hemisphere summer by 0.5–1.0 mm day−1. This leads to warmer daily maximum and minimum temperatures by up to 1.0 °C and warmer extreme maximum temperatures by up to 1.5 °C. These changes generally improve the climate model's climatology of warm extremes and improve existing biases by 10–20 %. The bias in minimum temperatures is however degraded but, overall, this is outweighed by the improvement in maximum temperatures as there is a net improvement in the diurnal temperature range in this region. In other regions such as parts of South and North America where ACCESSv1.3b has known large positive biases in both maximum and minimum temperatures (~ 5 to 10 °C), the new scheme degrades this bias by up to 1 °C. We conclude that, although several large biases remain in ACCESSv1.3b for temperature extremes, the improvements in the global climate model over large parts of the boreal forests during the Northern Hemisphere summer which result from the new stomatal scheme, constrained by a global synthesis of experimental data, provide a valuable advance in the long-term development of the ACCESS modelling system.

Abstract. We implement a new stomatal conductance model, based on the optimality approach, within the Community Atmosphere Biosphere Land Exchange (CABLE) land surface model. Coupled land-atmosphere simulations are then performed using CABLE within the Australian Community Climate and Earth Systems Simulator (ACCESS) with prescribed sea surface temperatures. As in most land surface models, the default stomatal conductance scheme only accounts for differences in model parameters in relation to the photosynthetic pathway, but not in relation to plant functional types. The new scheme allows model parameters to vary by plant functional type, based on a global synthesis of observations of stomatal conductance under different climate regimes over a wide range of species. We show that the new scheme reduces the latent heat flux from the land surface over the boreal forests during the Northern Hemisphere summer by 0.5 to 1.0 mm day-1. This leads to warmer daily maximum and minimum temperatures by up to 1.0 °C and warmer extreme maximum temperatures by up to 1.5 °C. These changes generally improve the climate model's climatology and improve existing biases by 10–20 %. The change in the surface energy balance also affects net primary productivity and the terrestrial carbon balance. We conclude that the improvements in the global climate model which result from the new stomatal scheme, constrained by a global synthesis of experimental data, provide a valuable advance in the long-term development of the ACCESS modelling system.

No abstract