<p>Terrestrial vegetation growth is hypothesised to increase under elevated atmospheric CO<sub>2</sub>, a process known as the CO<sub>2</sub> fertilisation effect. However, the magnitude of this effect and its long-term sustainability remains uncertain. One of the main limitations to the CO2&#160; fertilisation effect is nutrient limitation to plant growth, in particular nitrogen (N) in temperate and boreal ecosystems. Recent studies have suggested that decreases in observed foliar N content (N%) and &#948;<sup>15</sup>N indicate widespread nitrogen limitation with increasing CO<sub>2</sub>&#160; concentrations. However, the factors driving these two variables, and especially the foliar &#948;<sup>15</sup>N values, are complex and can be caused by a number of processes. On one hand, if the observed trends reflect nutrient limitation, this limitation can be caused by either CO<sub>2</sub> or warming driven growth. On the other hand, it is possible that nutrient limitation does not occur to its full extent due to plant plastic responses to alleviate nutrient limitation, causing a decrease in N%, but changes in the anthropogenic N deposition 15N signal cause the observed &#948;<sup>15</sup>N trend. In reality, it is likely that all these factors contribute to the observed trends. To understand ecosystem dynamics it is important to disentangle the processes behind these signals which is very difficult based on observational datasets only.</p><p>We use a novel land surface model to explore the causes behind the observed trends in foliar N% and &#948;<sup>15</sup>N. The QUINCY (QUantifying Interactions between terrestrial Nutrient CYcles and the climate system) model&#160; has the unique capacity to track ecologically relevant isotopic composition, including <sup>15</sup>N in plant and soil pools. The model also includes a realistic representation of plant plastic acclimation processes, specifically a representation of nitrogen allocation to and inside the canopy in response to nitrogen availability, so implicitly to changes in CO<sub>2&#160;</sub> concentrations. We test the different hypotheses above behind the observed changes in N% and &#948;<sup>15</sup>N separately and quantify the contribution of each of the factors towards the observed trend. We then test the different hypotheses against existing observations of N% and &#948;<sup>15</sup>N from the ICP Forests database and other published datasets such as the global dataset of Craine et al. 2018.</p><p>Our study showcases the use of an isotope-enabled land surface model in conjunction with long-term observations to strengthen our understanding of the ecosystem processes behind the observed trends.</p>