<p><strong>Abstract.</strong> We examined the net terrestrial carbon flux to the atmosphere (<i>F</i>_TA) simulated by nine models from the TRENDY dynamic global vegetation model project for its seasonal cycle and amplitude trend during 1961–2012. While some models exhibit similar phase and amplitude compared to atmospheric inversions, with spring drawdown and autumn rebound, others tend to rebound early in summer. The model ensemble mean underestimates the magnitude of the seasonal

more »

... cle by 40<span class="thinspace"></span>% compared to atmospheric inversions. Global <i>F</i>_TA amplitude increase (19<span class="thinspace"></span>±<span class="thinspace"></span>8<span class="thinspace"></span>%) and its decadal variability from the model ensemble are generally consistent with constraints from surface atmosphere observations. However, models disagree on attribution of this long-term amplitude increase, with factorial experiments attributing 83<span class="thinspace"></span>±<span class="thinspace"></span>56<span class="thinspace"></span>%, −3<span class="thinspace"></span>±<span class="thinspace"></span>74 and 20<span class="thinspace"></span>±<span class="thinspace"></span>30<span class="thinspace"></span>% to rising CO₂, climate change and land use/cover change, respectively. Seven out of the nine models suggest that CO₂ fertilization is the strongest control – with the notable exception of VEGAS, which attributes approximately equally to the three factors. Generally, all models display an enhanced seasonality over the boreal region in response to high-latitude warming, but a negative climate contribution from part of the Northern Hemisphere temperate region, and the net result is a divergence over climate change effect. Six of the nine models show that land use/cover change amplifies the seasonal cycle of global <i>F</i>_TA: some are due to forest regrowth, while others are caused by crop expansion or agricultural intensification, as revealed by their divergent spatial patterns. We also discovered a moderate cross-model correlation between <i>F</i>_TA amplitude increase and increase in land carbon sink (<i>R</i>² = <span class="thinspace"></span>0.61). Our results suggest that models can show similar results in some benchmarks with different underlying mechanisms; therefore, the spatial traits of CO₂ fertilization, climate change and land use/cover changes are crucial in determining the right mechanisms in seasonal carbon cycle change as well as mean sink change.</p>