Evaluating two land surface models for Brazil using a full carbon cycle benchmark with uncertainties

Abstract

AbstractForecasts of tropical ecosystem C cycling diverge among models due to differences in simulation of internal processes such as turnover, or transit times, of carbon pools. Estimates of these processes for the recent past are needed to test model representations, and so build confidence in model forecasts within and across biomes. Here, we evaluate carbon cycle process representation in two land surface models [Joint UK Land Environment Simulator (JULES) and Integrated Model of Land Surface Processes (INLAND)] for the period 2001–10 across Brazilian biomes. Model outputs are evaluated using the ILAMB system. Probabilistic benchmarking data were created using the carbon data model framework that assimilates observational times series of leaf area index and maps of woody biomass and soil C. New custom uncertainty metrics assess if models are within benchmark uncertainties. Simulations are better in homogeneous areas of vegetation type, and are less robust at ecotones between biomes, likely due to disturbance effects and parameter errors. Gross biosphere‐atmosphere fluxes are robustly modelled across Brazil. However, benchmark uncertainty is too high on net ecosystem exchange to provide an accurate evaluation of the models. The LSMs have significant differences in internal carbon allocation and the dynamics of the different C pools. JULES models dead C stocks more accurately while living C stocks are best resolved for INLAND. JULES' over‐estimate of the C wood pool results from over‐estimation of both inputs to wood and the transit time of wood. INLAND's under‐estimate of dead C stocks arises from an under‐estimate of the transit time of dead organic matter. The models are better at simulating annual averages than seasonal variation of fluxes. Analyses of monthly net C exchanges show that INLAND correctly simulates seasonality, but over‐estimates amplitudes, whereas JULES correctly simulates the annual amplitudes, but is out of phase with the benchmark.