How drought events during the last century have impacted biomass carbon in Amazonian rainforests.

During the last two decades, inventory data show that droughts have reduced biomass carbon sink of the Amazon forest by causing mortality to exceed growth. However, process-based models have struggled to include drought-induced responses of growth and mortality and have not been evaluated against plot data. A process-based model, ORCHIDEE-CAN-NHA, including forest demography with tree cohorts, plant hydraulic architecture and drought-induced tree mortality, was applied over Amazonia rainforests forced by gridded climate fields and rising CO2 from 1901 to 2019. The model reproduced the decelerating signal of net carbon sink and drought sensitivity of aboveground biomass (AGB) growth and mortality observed at forest plots across selected Amazon intact forests for 2005 and 2010. We predicted a larger mortality rate and a more negative sensitivity of the net carbon sink during the 2015/16 El Niño compared with the former droughts. 2015/16 was indeed the most severe drought since 1901 regarding both AGB loss and area experiencing a severe carbon loss. We found that even if climate change did increase mortality, elevated CO2 contributed to balance the biomass mortality, since CO2 -induced stomatal closure reduces transpiration, thus, offsets increased transpiration from CO2 -induced higher foliage area.

Compensatory water effects link yearly global land CO2 sink changes to temperature.

Large interannual variations in the measured growth rate of atmospheric carbon dioxide (CO2) originate primarily from fluctuations in carbon uptake by land ecosystems. It remains uncertain, however, to what extent temperature and water availability control the carbon balance of land ecosystems across spatial and temporal scales. Here we use empirical models based on eddy covariance data and process-based models to investigate the effect of changes in temperature and water availability on gross primary productivity (GPP), terrestrial ecosystem respiration (TER) and net ecosystem exchange (NEE) at local and global scales. We find that water availability is the dominant driver of the local interannual variability in GPP and TER. To a lesser extent this is true also for NEE at the local scale, but when integrated globally, temporal NEE variability is mostly driven by temperature fluctuations. We suggest that this apparent paradox can be explained by two compensatory water effects. Temporal water-driven GPP and TER variations compensate locally, dampening water-driven NEE variability. Spatial water availability anomalies also compensate, leaving a dominant temperature signal in the year-to-year fluctuations of the land carbon sink. These findings help to reconcile seemingly contradictory reports regarding the importance of temperature and water in controlling the interannual variability of the terrestrial carbon balance. Our study indicates that spatial climate covariation drives the global carbon cycle response.

Rainfall manipulation experiments as simulated by terrestrial biosphere models: Where do we stand?

Changes in rainfall amounts and patterns have been observed and are expected to continue in the near future with potentially significant ecological and societal consequences. Modelling vegetation responses to changes in rainfall is thus crucial to project water and carbon cycles in the future. In this study, we present the results of a new model-data intercomparison project, where we tested the ability of 10 terrestrial biosphere models to reproduce the observed sensitivity of ecosystem productivity to rainfall changes at 10 sites across the globe, in nine of which, rainfall exclusion and/or irrigation experiments had been performed. The key results are as follows: (a) Inter-model variation is generally large and model agreement varies with timescales. In severely water-limited sites, models only agree on the interannual variability of evapotranspiration and to a smaller extent on gross primary productivity. In more mesic sites, model agreement for both water and carbon fluxes is typically higher on fine (daily-monthly) timescales and reduces on longer (seasonal-annual) scales. (b) Models on average overestimate the relationship between ecosystem productivity and mean rainfall amounts across sites (in space) and have a low capacity in reproducing the temporal (interannual) sensitivity of vegetation productivity to annual rainfall at a given site, even though observation uncertainty is comparable to inter-model variability. (c) Most models reproduced the sign of the observed patterns in productivity changes in rainfall manipulation experiments but had a low capacity in reproducing the observed magnitude of productivity changes. Models better reproduced the observed productivity responses due to rainfall exclusion than addition. (d) All models attribute ecosystem productivity changes to the intensity of vegetation stress and peak leaf area, whereas the impact of the change in growing season length is negligible. The relative contribution of the peak leaf area and vegetation stress intensity was highly variable among models.

Effect of climate change, CO2 trends, nitrogen addition, and land-cover and management intensity changes on the carbon balance of European grasslands.

Several lines of evidence point to European managed grassland ecosystems being a sink of carbon. In this study, we apply ORCHIDEE-GM a process-based carbon cycle model that describes specific management practices of pastures and the dynamics of carbon cycling in response to changes in climatic and biogeochemical drivers. The model is used to simulate changes in the carbon balance [i.e., net biome production (NBP)] of European grasslands over 1991-2010 on a 25 km × 25 km grid. The modeled average trend in NBP is 1.8-2.0 g C m(-2)  yr(-2) during the past two decades. Attribution of this trend suggests management intensity as the dominant driver explaining NBP trends in the model (36-43% of the trend due to all drivers). A major change in grassland management intensity has occurred across Europe resulting from reduced livestock numbers. This change has 'inadvertently' enhanced soil C sequestration and reduced N2 O and CH4 emissions by 1.2-1.5 Gt CO2 -equivalent, offsetting more than 7% of greenhouse gas emissions in the whole European agricultural sector during the period 1991-2010. Land-cover change, climate change and rising CO2 also make positive and moderate contributions to the NBP trend (between 24% and 31% of the trend due to all drivers). Changes in nitrogen addition (including fertilization and atmospheric deposition) are found to have only marginal net effect on NBP trends. However, this may not reflect reality because our model has only a very simple parameterization of nitrogen effects on photosynthesis. The sum of NBP trends from each driver is larger than the trend obtained when all drivers are varied together, leaving a residual - nonattributed - term (22-26% of the trend due to all drivers) indicating negative interactions between drivers.

Recent decline in the global land evapotranspiration trend due to limited moisture supply.

More than half of the solar energy absorbed by land surfaces is currently used to evaporate water. Climate change is expected to intensify the hydrological cycle and to alter evapotranspiration, with implications for ecosystem services and feedback to regional and global climate. Evapotranspiration changes may already be under way, but direct observational constraints are lacking at the global scale. Until such evidence is available, changes in the water cycle on land−a key diagnostic criterion of the effects of climate change and variability−remain uncertain. Here we provide a data-driven estimate of global land evapotranspiration from 1982 to 2008, compiled using a global monitoring network, meteorological and remote-sensing observations, and a machine-learning algorithm. In addition, we have assessed evapotranspiration variations over the same time period using an ensemble of process-based land-surface models. Our results suggest that global annual evapotranspiration increased on average by 7.1 ± 1.0 millimetres per year per decade from 1982 to 1997. After that, coincident with the last major El Niño event in 1998, the global evapotranspiration increase seems to have ceased until 2008. This change was driven primarily by moisture limitation in the Southern Hemisphere, particularly Africa and Australia. In these regions, microwave satellite observations indicate that soil moisture decreased from 1998 to 2008. Hence, increasing soil-moisture limitations on evapotranspiration largely explain the recent decline of the global land-evapotranspiration trend. Whether the changing behaviour of evapotranspiration is representative of natural climate variability or reflects a more permanent reorganization of the land water cycle is a key question for earth system science.

The greenhouse gas balance of European grasslands.

The greenhouse gas (GHG) balance of European grasslands (EU-28 plus Norway and Switzerland), including CO2 , CH4 and N2 O, is estimated using the new process-based biogeochemical model ORCHIDEE-GM over the period 1961-2010. The model includes the following: (1) a mechanistic representation of the spatial distribution of management practice; (2) management intensity, going from intensively to extensively managed; (3) gridded simulation of the carbon balance at ecosystem and farm scale; and (4) gridded simulation of N2 O and CH4 emissions by fertilized grassland soils and livestock. The external drivers of the model are changing animal numbers, nitrogen fertilization and deposition, land-use change, and variable CO2 and climate. The carbon balance of European grassland (NBP) is estimated to be a net sink of 15 ± 7 g C m(-2 ) year(-1) during 1961-2010, equivalent to a 50-year continental cumulative soil carbon sequestration of 1.0 ± 0.4 Pg C. At the farm scale, which includes both ecosystem CO2 fluxes and CO2 emissions from the digestion of harvested forage, the net C balance is roughly halved, down to a small sink, or nearly neutral flux of 8 g C m(-2 ) year(-1) . Adding CH4 and N2 O emissions to net ecosystem exchange to define the ecosystem-scale GHG balance, we found that grasslands remain a net GHG sink of 19 ± 10 g C-CO2 equiv. m(-2 ) year(-1) , because the CO2 sink offsets N2 O and grazing animal CH4 emissions. However, when considering the farm scale, the GHG balance (NGB) becomes a net GHG source of -50 g C-CO2 equiv. m(-2 ) year(-1) . ORCHIDEE-GM simulated an increase in European grassland NBP during the last five decades. This enhanced NBP reflects the combination of a positive trend of net primary production due to CO2 , climate and nitrogen fertilization and the diminishing requirement for grass forage due to the Europe-wide reduction in livestock numbers.

Global critical soil moisture thresholds of plant water stress.

During extensive periods without rain, known as dry-downs, decreasing soil moisture (SM) induces plant water stress at the point when it limits evapotranspiration, defining a critical SM threshold (θcrit). Better quantification of θcrit is needed for improving future projections of climate and water resources, food production, and ecosystem vulnerability. Here, we combine systematic satellite observations of the diurnal amplitude of land surface temperature (dLST) and SM during dry-downs, corroborated by in-situ data from flux towers, to generate the observation-based global map of θcrit. We find an average global θcrit of 0.19 m3/m3, varying from 0.12 m3/m3 in arid ecosystems to 0.26 m3/m3 in humid ecosystems. θcrit simulated by Earth System Models is overestimated in dry areas and underestimated in wet areas. The global observed pattern of θcrit reflects plant adaptation to soil available water and atmospheric demand. Using explainable machine learning, we show that aridity index, leaf area and soil texture are the most influential drivers. Moreover, we show that the annual fraction of days with water stress, when SM stays below θcrit, has increased in the past four decades. Our results have important implications for understanding the inception of water stress in models and identifying SM tipping points.

Climate change impacts on tree ranges: model intercomparison facilitates understanding and quantification of uncertainty.

Model-based projections of shifts in tree species range due to climate change are becoming an important decision support tool for forest management. However, poorly evaluated sources of uncertainty require more scrutiny before relying heavily on models for decision-making. We evaluated uncertainty arising from differences in model formulations of tree response to climate change based on a rigorous intercomparison of projections of tree distributions in France. We compared eight models ranging from niche-based to process-based models. On average, models project large range contractions of temperate tree species in lowlands due to climate change. There was substantial disagreement between models for temperate broadleaf deciduous tree species, but differences in the capacity of models to account for rising CO(2) impacts explained much of the disagreement. There was good quantitative agreement among models concerning the range contractions for Scots pine. For the dominant Mediterranean tree species, Holm oak, all models foresee substantial range expansion.

Legislative and functional aspects of different metrics used for ozone risk assessment to forests.

Surface ozone (O3) is a threat to forests by decreasing photosynthesis and, consequently, influencing the strength of land carbon sink. However, due to the lack of continuous surface O3 measurements, observational-based assessments of O3 impacts on forests are largely missing at hemispheric to global scales. Currently, some metrics are used for regulatory purposes by governments or national agencies to protect forests against the negative impacts of ozone: in particular, both Europe and United States (US) makes use of two different exposure-based metrics, i.e. AOT40 and W126, respectively. However, because of some limitations in these metrics, a new standard is under consideration by the European Union (EU) to replace the current exposure metric. We analyse here the different air quality standards set or proposed for use in Europe and in the US to protect forests from O3 and to evaluate their spatial and temporal consistency while assessing their effectiveness in protecting northern-hemisphere forests. Then, we compare their results with the information obtained from a complex land surface model (ORCHIDEE). We find that present O3 uptake decreases gross primary production (GPP) in 37.7% of the NH forested area of northern hemisphere with a mean loss of 2.4% year-1. We show how the proposed US (W126) and the currently used European (AOT40) air quality standards substantially overestimate the extension of potential vulnerable regions, predicting that 46% and 61% of the Northern Hemisphere (NH) forested area are at risk of O3 pollution. Conversely, the new proposed European standard (POD1) identifies lower extension of vulnerability regions (39.6%).

Quantifying the benefits of reducing synthetic nitrogen application policy on ecosystem carbon sequestration and biodiversity.

Synthetic Nitrogen (N) usage in agriculture has greatly increased food supply over the past century. However, the intensive use of N fertilizer is nevertheless the source of numerous environmental issues and remains a major challenge for policymakers to understand, measure, and quantify the interactions and trade-offs between ecosystem carbon and terrestrial biodiversity loss. In this study, we investigate the impacts of a public policy scenario that aims to halve N fertilizer application across European Union (EU) agriculture on both carbon (C) sequestration and biodiversity changes. We quantify the impacts by integrating two economic models with an agricultural land surface model and a terrestrial biodiversity model (that uses data from a range of taxonomic groups, including plants, fungi, vertebrates and invertebrates). Here, we show that the two economic scenarios lead to different outcomes in terms of C sequestration potential and biodiversity. Land abandonment associated with increased fertilizer price scenario facilitates higher C sequestration in soils (+ 1014 MtC) and similar species richness levels (+ 1.9%) at the EU scale. On the other hand, the more extensive crop production scenario is associated with lower C sequestration potential in soils (- 97 MtC) and similar species richness levels (- 0.4%) because of a lower area of grazing land. Our results therefore highlight the complexity of the environmental consequences of a nitrogen reduction policy, which will depend fundamentally on how the economic models used to project consequences.

Disentangling the Impacts of Anthropogenic Aerosols on Terrestrial Carbon Cycle During 1850-2014.

Aerosols have a dimming and cooling effect and change hydrological regimes, thus affecting carbon fluxes, which are sensitive to climate. Aerosols also scatter sunlight, which increases the fraction of diffuse radiation, increasing photosynthesis. There remains no clear conclusion whether the impact of aerosols on land carbon fluxes is larger through diffuse radiation change than through changes in other climate variables. In this study, we quantified the overall physical impacts of anthropogenic aerosols on land C fluxes and explored the contribution from each factor using a set of factorial simulations driven by climate and aerosol data from the IPSL-CM6A-LR experiments during 1850-2014. A newly developed land surface model which distinguishes diffuse and direct radiation in canopy radiation transmission, ORCHIDEE_DF, was used. Specifically, a subgrid scheme was developed to distinguish the cloudy and clear sky conditions. We found that anthropogenic aerosol emissions since 1850 cumulatively enhanced the land C sink by 22.6 PgC. Seventy-eight percent of this C sink enhancement is contributed by aerosol-induced increase in the diffuse radiation fraction, much larger than the effect of the aerosol-induced dimming. The cooling of anthropogenic aerosols has different impacts in different latitudes but overall increases the global land C sink. The dominant role of diffuse radiation changes found in this study implies that future aerosol emissions may have a much stronger impacts on the C cycle through changing radiation quality than through changing climate alone. Earth system models need to consider the diffuse radiation fertilization effect to better evaluate the impacts of climate change mitigation scenarios.

Orbital variations as a major driver of climate and biome distribution during the greenhouse to icehouse transition.

Recent studies suggest increasing sensitivity to orbital variations across the Eocene-Oligocene greenhouse to icehouse climate transition. However, climate simulations and paleoenvironmental studies mostly provide snapshots of the past climate, therefore overlooking the role of this short-term variability in driving major environmental changes and possibly biasing model-data comparisons. We address this problem by performing numerical simulations describing the end-members of eccentricity, obliquity, and precession. The orbitally induced biome variability obtained in our simulations allows to reconcile previous apparent mismatch between models and paleobotanical compilations. We show that precession-driven intermittent monsoon-like climate may have taken place during the Eocene, resulting in biomes shifting from shrubland to tropical forest in the intertropical convergence zone migration region. Our Oligocene simulations suggest that, along with decreased pCO2, orbital variations crucially modulated major faunal dispersal events around the EOT such as the Grande Coupure by creating and fragmenting the biome corridors along several key land bridges.

Climate warming from managed grasslands cancels the cooling effect of carbon sinks in sparsely grazed and natural grasslands.

Grasslands absorb and release carbon dioxide (CO2), emit methane (CH4) from grazing livestock, and emit nitrous oxide (N2O) from soils. Little is known about how the fluxes of these three greenhouse gases, from managed and natural grasslands worldwide, have contributed to past climate change, or the roles of managed pastures versus natural grasslands. Here, global trends and regional patterns of the full greenhouse gas balance of grasslands are estimated for the period 1750 to 2012. A new spatially explicit land surface model is applied, to separate the direct effects of human activities from land management and the indirect effects from climate change, increasing CO2 and regional changes in nitrogen deposition. Direct human management activities are simulated to have caused grasslands to switch from a sink to a source of greenhouse gas, because of increased livestock numbers and accelerated conversion of natural lands to pasture. However, climate change drivers contributed a net carbon sink in soil organic matter, mainly from the increased productivity of grasslands due to increased CO2 and nitrogen deposition. The net radiative forcing of all grasslands is currently close to neutral, but has been increasing since the 1960s. Here, we show that the net global climate warming caused by managed grassland cancels the net climate cooling from carbon sinks in sparsely grazed and natural grasslands. In the face of future climate change and increased demand for livestock products, these findings highlight the need to use sustainable management to preserve and enhance soil carbon storage in grasslands and to reduce greenhouse gas emissions from managed grasslands.

Historical phenology: grape ripening as a past climate indicator.

French records of grape-harvest dates in Burgundy were used to reconstruct spring-summer temperatures from 1370 to 2003 using a process-based phenology model developed for the Pinot Noir grape. Our results reveal that temperatures as high as those reached in the 1990s have occurred several times in Burgundy since 1370. However, the summer of 2003 appears to have been extraordinary, with temperatures that were probably higher than in any other year since 1370.

Vegetation structural change since 1981 significantly enhanced the terrestrial carbon sink.

Satellite observations show that leaf area index (LAI) has increased globally since 1981, but the impact of this vegetation structural change on the global terrestrial carbon cycle has not been systematically evaluated. Through process-based diagnostic ecosystem modeling, we find that the increase in LAI alone was responsible for 12.4% of the accumulated terrestrial carbon sink (95 ± 5 Pg C) from 1981 to 2016, whereas other drivers of CO2 fertilization, nitrogen deposition, and climate change (temperature, radiation, and precipitation) contributed to 47.0%, 1.1%, and -28.6% of the sink, respectively. The legacy effects of past changes in these drivers prior to 1981 are responsible for the remaining 65.5% of the accumulated sink from 1981 to 2016. These results refine the attribution of the land sink to the various drivers and would help constrain prognostic models that often have large uncertainties in simulating changes in vegetation and their impacts on the global carbon cycle.

The large mean body size of mammalian herbivores explains the productivity paradox during the Last Glacial Maximum.

Large herbivores are a major agent in ecosystems, influencing vegetation structure, and carbon and nutrient flows. During the last glacial period, a mammoth steppe ecosystem prevailed in the unglaciated northern lands, supporting a high diversity and density of megafaunal herbivores. The apparent discrepancy between abundant megafauna and the expected low vegetation productivity under a generally harsher climate with a lower CO2 concentration, termed the productivity paradox, requires large-scale quantitative analysis using process-based ecosystem models. However, most of the current global dynamic vegetation models (DGVMs) lack explicit representation of large herbivores. Here we incorporated a grazing module in a DGVM based on physiological and demographic equations for wild large grazers, taking into account feedbacks of large grazers on vegetation. The model was applied globally for present-day and the Last Glacial Maximum (LGM). The present-day results of potential grazer biomass, combined with an empirical land-use map, infer a reduction in wild grazer biomass by 79-93% owing to anthropogenic land replacement of natural grasslands. For the LGM, we find that the larger mean body size of mammalian herbivores than today is the crucial clue to explain the productivity paradox, due to a more efficient exploitation of grass production by grazers with a large body size.

Impact of the 2015/2016 El Niño on the terrestrial carbon cycle constrained by bottom-up and top-down approaches.

Evaluating the response of the land carbon sink to the anomalies in temperature and drought imposed by El Niño events provides insights into the present-day carbon cycle and its climate-driven variability. It is also a necessary step to build confidence in terrestrial ecosystems models' response to the warming and drying stresses expected in the future over many continents, and particularly in the tropics. Here we present an in-depth analysis of the response of the terrestrial carbon cycle to the 2015/2016 El Niño that imposed extreme warming and dry conditions in the tropics and other sensitive regions. First, we provide a synthesis of the spatio-temporal evolution of anomalies in net land-atmosphere CO2 fluxes estimated by two in situ measurements based on atmospheric inversions and 16 land-surface models (LSMs) from TRENDYv6. Simulated changes in ecosystem productivity, decomposition rates and fire emissions are also investigated. Inversions and LSMs generally agree on the decrease and subsequent recovery of the land sink in response to the onset, peak and demise of El Niño conditions and point to the decreased strength of the land carbon sink: by 0.4-0.7 PgC yr-1 (inversions) and by 1.0 PgC yr-1 (LSMs) during 2015/2016. LSM simulations indicate that a decrease in productivity, rather than increase in respiration, dominated the net biome productivity anomalies in response to ENSO throughout the tropics, mainly associated with prolonged drought conditions.This article is part of a discussion meeting issue 'The impact of the 2015/2016 El Niño on the terrestrial tropical carbon cycle: patterns, mechanisms and implications'.

Future productivity and phenology changes in European grasslands for different warming levels: implications for grassland management and carbon balance.

BACKGROUND: Europe has warmed more than the global average (land and ocean) since pre-industrial times, and is also projected to continue to warm faster than the global average in the twenty-first century. According to the climate models ensemble projections for various climate scenarios, annual mean temperature of Europe for 2071-2100 is predicted to be 1-5.5 °C higher than that for 1971-2000. Climate change and elevated CO2 concentration are anticipated to affect grassland management and livestock production in Europe. However, there has been little work done to quantify the European-wide response of grassland to future climate change. Here we applied ORCHIDEE-GM v2.2, a grid-based model for managed grassland, over European grassland to estimate the impacts of future global change. RESULTS: Increases in grassland productivity are simulated in response to future global change, which are mainly attributed to the simulated fertilization effect of rising CO2. The results show significant phenology shifts, in particular an earlier winter-spring onset of grass growth over Europe. A longer growing season is projected over southern and southeastern Europe. In other regions, summer drought causes an earlier end to the growing season, overall reducing growing season length. Future global change allows an increase of management intensity with higher than current potential annual grass forage yield, grazing capacity and livestock density, and a shift in seasonal grazing capacity. We found a continual grassland soil carbon sink in Mediterranean, Alpine, North eastern, South eastern and Eastern regions under specific warming level (SWL) of 1.5 and 2 °C relative to pre-industrial climate. However, this carbon sink is found to saturate, and gradually turn to a carbon source at warming level reaching 3.5 °C. CONCLUSIONS: This study provides a European-wide assessment of the future changes in productivity and phenology of grassland, and their consequences for the management intensity and the carbon balance. The simulated productivity increase in response to future global change enables an intensification of grassland management over Europe. However, the simulated increase in the interannual variability of grassland productivity over some regions may reduce the farmers' ability to take advantage of the increased long-term mean productivity in the face of more frequent, and more severe drops of productivity in the future.

Comparing and evaluating process-based ecosystem model predictions of carbon and water fluxes in major European forest biomes.

Process-based models can be classified into: (a) terrestrial biogeochemical models (TBMs), which simulate fluxes of carbon, water and nitrogen coupled within terrestrial ecosystems, and (b) dynamic global vegetation models (DGVMs), which further couple these processes interactively with changes in slow ecosystem processes depending on resource competition, establishment, growth and mortality of different vegetation types. In this study, four models - RHESSys, GOTILWA+, LPJ-GUESS and ORCHIDEE - representing both modelling approaches were compared and evaluated against benchmarks provided by eddy-covariance measurements of carbon and water fluxes at 15 forest sites within the EUROFLUX project. Overall, model-measurement agreement varied greatly among sites. Both modelling approaches have somewhat different strengths, but there was no model among those tested that universally performed well on the two variables evaluated. Small biases and errors suggest that ORCHIDEE and GOTILWA+ performed better in simulating carbon fluxes while LPJ-GUESS and RHESSys did a better job in simulating water fluxes. In general, the models can be considered as useful tools for studies of climate change impacts on carbon and water cycling in forests. However, the various sources of variation among models simulations and between models simulations and observed data described in this study place some constraints on the results and to some extent reduce their reliability. For example, at most sites in the Mediterranean region all models generally performed poorly most likely because of problems in the representation of water stress effects on both carbon uptake by photosynthesis and carbon release by heterotrophic respiration (Rh ). The use of flux data as a means of assessing key processes in models of this type is an important approach to improving model performance. Our results show that the models have value but that further model development is necessary with regard to the representation of the some of the key ecosystem processes.

Modeling climate change effects on the potential production of French plains forests at the sub-regional level.

We modeled the effects of climate change and two forest management scenarios on wood production and forest carbon balance in French forests using process-based models of forest growth. We combined data from the national forest inventory and soil network survey, which were aggregated over a 50 x 50-km grid, i.e., the spatial resolution of the climate scenario data. We predicted and analyzed the climate impact on potential forest production over the period 1960-2100. All models predicted a slight increase in potential forest yield until 2030-2050, followed by a plateau or a decline around 2070-2100, with overall, a greater increase in yield in northern France than in the south. Gross and net primary productivities were more negatively affected by soil water and atmospheric water vapor saturation deficits in western France because of a more pronounced shift in seasonal rainfall from summer to winter. The rotation-averaged values of carbon flux and production for different forest management options were estimated during four years (1980, 2015, 2045 and 2080). Predictions were made using a two-dimensional matrix covering the range of local soil and climate conditions. The changes in ecosystem fluxes and forest production were explained by the counterbalancing effect of rising CO2 concentration and increasing water deficit. The effect of climate change decreased with rotation length from short rotations with high production rates and low standing biomasses to long rotations with low productivities and greater standing biomasses. Climate effects on productivity, both negative and positive, were greatest on high fertility sites. Forest productivity in northern France was enhanced by climate change, increasingly from west to east, whereas in the southwestern Atlantic region, productivity was reduced by climate change to an increasing degree from west to east.