A “toy” model of biogeochemical dynamics on climate gradients

We used a simple “toy” model to aid in the evaluation of the controls of biogeochemical patterns along a climate gradient. The model includes simplified treatments of water balance (precipitation minus Potential Evapotranspiration), leaching, weathering of cation- and P-bearing minerals, N cycling and loss, biomass production, and biological N fixation. We use δ15N as a central integrator of biogeochemical processes, because δ15N integrates multiple pathways of N input, output, and transformation in ecosystems. The model simulated the location and magnitude of a peak in δ15N on a gradient on Kohala Volcano, Hawai‘i which peaked ~  + 14 ‰ in sites receiving ~ 3.5 cm/month average precipitation (− 1300 mm/year water balance); the model also captured a peak in total P in surface soil at intermediate levels of precipitation and water balance, and other biogeochemical features on the gradient. We then applied the model to understanding the patterns of and mechanisms underlying nutrient limitation to net primary production (NPP) and plant biomass on the gradient, testing for the existence and extent of N and P limitation by simulated additions of N and/or P in the model. Both a simulated symbiotic biological N fixer and a simulated non-fixer were limited by P supply across the gradient; the non-fixer was independently limited by N supply in wetter sites. By running the toy model with and without the influence of temperature, we demonstrated that water is the most important factor shaping biogeochemical patterns on this gradient.