Increasing concentrations of atmospheric carbon dioxide are expected to affect carbon assimilation and evapotranspiration (ET), ultimately driving changes in plant growth, hydrology, and the global carbon balance. Direct leaf biochemical effects have been widely investigated, whereas indirect effects, although documented, elude explicit quantification in experiments. Here, we used a mechanistic model to investigate the relative contributions of direct (through carbon assimilation) and indirect (via soil moisture savings due to stomatal closure, and changes in leaf area index) effects of elevated CO 2 across a variety of ecosystems. We specifically determined which ecosystems and climatic conditions maximize the indirect effects of elevated CO 2 . The simulations suggest that the indirect effects of elevated CO 2 on net primary productivity are large and variable, ranging from less than 10% to more than 100% of the size of direct effects. For ET, indirect effects were, on average, 65% of the size of direct effects. Indirect effects tended to be considerably larger in water-limited ecosystems. As a consequence, the total CO 2 effect had a significant, inverse relationship with the wetness index and was directly related to vapor pressure deficit. These results have major implications for our understanding of the CO 2 response of ecosystems and for global projections of CO 2 fertilization, because, although direct effects are typically understood and easily reproducible in models, simulations of indirect effects are far more challenging and difficult to constrain. Our findings also provide an explanation for the discrepancies between experiments in the total CO 2 effect on net primary productivity.he leaf-level response to elevated CO 2 (eCO 2 ) is well known: at current CO 2 levels, photosynthesis of C 3 plants is not saturated, whereas, for C 4 plants, it is close to saturation (e.g., refs. 1-4). If acclimation is limited, leaf-level carbon assimilation of C 3 plants will increase as the CO 2 concentration increases, as shown by, among others, observations in Free-Air CO 2 Enrichment (FACE) experiments (5-7). Concurrently, stomatal conductance decreases consistently with eCO 2 in most species (8-11). Even though the leaf-level responses are well characterized and quantifiable, the ecosystem response to eCO 2 remains considerably more uncertain and difficult to predict (12-17). This discrepancy is not simply a consequence of the uncertainty in scaling up from leaf to canopy and ecosystem but derives from indirect effects and feedbacks that may lead to an amplification or dampening of the direct leaf-level response to eCO 2 .Indirect effects may be related to (i) modifications of plant water status through changes in soil moisture within the root zone, which occur as a consequence of stomatal closure; (ii) changes in Leaf Area Index (LAI), root biomass and depth, and canopy structure; (iii) limitations due to soil nutrient scarcity or plant incapability to take up nutrients at a rate sufficient to support enhanced c...