2008, Vetter, M., Churkina, G., Jung, M., Reichstein, M., Zaehle, S., Bondeau, A., Chen, Y., Ciais, P., Feser, F., Freibauer, A., Geyer, R., Jones, C., Papale, D., Tenhunen, J., Tomelleri, E., Trusilova, K., Viovy, N., Heimann, M.

Globally, the year 2003 is associated with one of the largest atmospheric CO2 rises on record. In the same year, Europe experienced an anomalously strong flux of CO2 from the land to the atmosphere associated with an exceptionally dry and hot summer in Western and Central Europe. In this study we analyze the magnitude of this carbon flux anomaly and key driving ecosystem processes using simulations of seven terrestrial ecosystem models of different complexity and types (process-oriented and diagnostic). We address the following questions: (1) how large were deviations in the net European carbon flux in 2003 relative to a short-term baseline (1998–2002) and to longer-term variations in annual fluxes (1980 to 2005), (2) which European regions exhibited the largest changes in carbon fluxes during the growing season 2003, and (3) which ecosystem processes controlled the carbon balance anomaly . In most models the prominence of 2003 anomaly in carbon fluxes declined with lengthening of the reference period from one year to 16 years. The 2003 anomaly for annual net carbon fluxes ranged between 0.35 and –0.63 Pg C for a reference period of one year and between 0.17 and –0.37 Pg C for a reference period of 16 years for the whole Europe. In Western and Central Europe, the anomaly in simulated net ecosystem productivity (NEP) over the growing season in 2003 was outside the 1σ variance bound of the carbon flux anomalies for 1980–2005 in all models. The estimated anomaly in net carbon flux ranged between –42 and –158 Tg C for Western Europe and between 24 and –129 Tg C for Central Europe depending on the model used. All models responded to a dipole pattern of the climate anomaly in 2003. In Western and Central Europe NEP was reduced due to heat and drought. In contrast, lower than normal temperatures and higher air humidity decreased NEP over Northeastern Europe. While models agree on the sign of changes in simulated NEP and gross primary productivity in 2003 over Western and Central Europe, models diverge in the estimates of anomalies in ecosystem respiration. Except for two process models which simulate respiration increase, most models simulated a decrease in ecosystem respiration in 2003. The diagnostic models showed a weaker decrease in ecosystem respiration than the process-oriented models. Based on the multi-model simulations we estimated the total carbon flux anomaly over the 2003 growing season in Europe to range between –0.02 and –0.27 Pg C relative to the net carbon flux in 1998–2002.

2015, Frieler, K., Levermann, A., Elliott, J., Heinke, J., Arneth, A., Bierkens, M.F.P., Ciais, P., Clark, D.B., Deryng, D., Döll, P., Falloon, P., Fekete, B., Folberth, C., Friend, A.D., Gellhorn, C., Gosling, S.N., Haddeland, I., Khabarov, N., Lomas, M., Masaki, Y., Nishina, K., Neumann, K., Oki, T., Pavlick, R., Ruane, A.C., Schmid, E., Schmitz, C., Stacke, T., Stehfest, E., Tang, Q., Wisser, D., Huber, V., Piontek, F., Warszawski, L., Schewe, J., Lotze-Campen, H., Schellnhuber, H.J.

Climate change and its impacts already pose considerable challenges for societies that will further increase with global warming (IPCC, 2014a, b). Uncertainties of the climatic response to greenhouse gas emissions include the potential passing of large-scale tipping points (e.g. Lenton et al., 2008; Levermann et al., 2012; Schellnhuber, 2010) and changes in extreme meteorological events (Field et al., 2012) with complex impacts on societies (Hallegatte et al., 2013). Thus climate change mitigation is considered a necessary societal response for avoiding uncontrollable impacts (Conference of the Parties, 2010). On the other hand, large-scale climate change mitigation itself implies fundamental changes in, for example, the global energy system. The associated challenges come on top of others that derive from equally important ethical imperatives like the fulfilment of increasing food demand that may draw on the same resources. For example, ensuring food security for a growing population may require an expansion of cropland, thereby reducing natural carbon sinks or the area available for bio-energy production. So far, available studies addressing this problem have relied on individual impact models, ignoring uncertainty in crop model and biome model projections. Here, we propose a probabilistic decision framework that allows for an evaluation of agricultural management and mitigation options in a multi-impact-model setting. Based on simulations generated within the Inter-Sectoral Impact Model Intercomparison Project (ISI-MIP), we outline how cross-sectorally consistent multi-model impact simulations could be used to generate the information required for robust decision making. Using an illustrative future land use pattern, we discuss the trade-off between potential gains in crop production and associated losses in natural carbon sinks in the new multiple crop- and biome-model setting. In addition, crop and water model simulations are combined to explore irrigation increases as one possible measure of agricultural intensification that could limit the expansion of cropland required in response to climate change and growing food demand. This example shows that current impact model uncertainties pose an important challenge to long-term mitigation planning and must not be ignored in long-term strategic decision making.

2014, Van Oijen, M., Balkovi, J., Beer, C., Cameron, D.R., Ciais, P., Cramer, W., Kato, T., Kuhnert, M., Martin, R., Myneni, R., Rammig, A., Rolinski, S., Soussana, J.-F., Thonicke, K., Van der Velde, M., Xu, L.

We analyse how climate change may alter risks posed by droughts to carbon fluxes in European ecosystems. The approach follows a recently proposed framework for risk analysis based on probability theory. In this approach, risk is quantified as the product of hazard probability and ecosystem vulnerability. The probability of a drought hazard is calculated here from the Standardized Precipitation–Evapotranspiration Index (SPEI). Vulnerability is calculated from the response to drought simulated by process-based vegetation models. We use six different models: three for generic vegetation (JSBACH, LPJmL, ORCHIDEE) and three for specific ecosystems (Scots pine forests: BASFOR; winter wheat fields: EPIC; grasslands: PASIM). The periods 1971–2000 and 2071–2100 are compared. Climate data are based on gridded observations and on output from the regional climate model REMO using the SRES A1B scenario. The risk analysis is carried out for ~ 18 000 grid cells of 0.25 × 0.25° across Europe. For each grid cell, drought vulnerability and risk are quantified for five seasonal variables: net primary and ecosystem productivity (NPP, NEP), heterotrophic respiration (Rh), soil water content and evapotranspiration. In this analysis, climate change leads to increased drought risks for net primary productivity in the Mediterranean area: five of the models estimate that risk will exceed 15%. The risks increase mainly because of greater drought probability; ecosystem vulnerability will increase to a lesser extent. Because NPP will be affected more than Rh, future carbon sequestration (NEP) will also be at risk predominantly in southern Europe, with risks exceeding 0.25 g C m−2 d−1 according to most models, amounting to reductions in carbon sequestration of 20 to 80%.

2014, Nishina, K., Ito, A., Beerling, D.J., Cadule, P., Ciais, P., Clark, D.B., Friend, A.D., Kahana, R., Kato, E., Keribin, R., Lucht, W., Lomas, M., Rademacher, T.T., Pavlick, R., Schaphoff, S., Vuichard, N., Warszawaski, L., Yokohata, T.

Soil organic carbon (SOC) is the largest carbon pool in terrestrial ecosystems and may play a key role in biospheric feedbacks with elevated atmospheric carbon dioxide (CO2) in a warmer future world. We examined the simulation results of seven terrestrial biome models when forced with climate projections from four representative-concentration-pathways (RCPs)-based atmospheric concentration scenarios. The goal was to specify calculated uncertainty in global SOC stock projections from global and regional perspectives and give insight to the improvement of SOC-relevant processes in biome models. SOC stocks among the biome models varied from 1090 to 2650 Pg C even in historical periods (ca. 2000). In a higher forcing scenario (i.e., RCP8.5), inconsistent estimates of impact on the total SOC (2099–2000) were obtained from different biome model simulations, ranging from a net sink of 347 Pg C to a net source of 122 Pg C. In all models, the increasing atmospheric CO2 concentration in the RCP8.5 scenario considerably contributed to carbon accumulation in SOC. However, magnitudes varied from 93 to 264 Pg C by the end of the 21st century across biome models. Using the time-series data of total global SOC simulated by each biome model, we analyzed the sensitivity of the global SOC stock to global mean temperature and global precipitation anomalies (ΔT and ΔP respectively) in each biome model using a state-space model. This analysis suggests that ΔT explained global SOC stock changes in most models with a resolution of 1–2 °C, and the magnitude of global SOC decomposition from a 2 °C rise ranged from almost 0 to 3.53 Pg C yr−1 among the biome models. However, ΔP had a negligible impact on change in the global SOC changes. Spatial heterogeneity was evident and inconsistent among the biome models, especially in boreal to arctic regions. Our study reveals considerable climate uncertainty in SOC decomposition responses to climate and CO2 change among biome models. Further research is required to improve our ability to estimate biospheric feedbacks through both SOC-relevant and vegetation-relevant processes.

2016, Jones, C.D., Ciais, P., Davis, S.J., Friedlingstein, P., Gasser, T., Peters, G.P., Rogelj, J., van Vuuren, D.P., Canadell, J.G., Cowie, A., Jackson, R.B., Jonas, M., Kriegler, E., Littleton, E., Lowe, J.A., Milne, J., Shrestha, G., Smith, P., Torvanger, A., Wiltshire, A.

Natural carbon sinks currently absorb approximately half of the anthropogenic CO2 emitted by fossil fuel burning, cement production and land-use change. However, this airborne fraction may change in the future depending on the emissions scenario. An important issue in developing carbon budgets to achieve climate stabilisation targets is the behaviour of natural carbon sinks, particularly under low emissions mitigation scenarios as required to meet the goals of the Paris Agreement. A key requirement for low carbon pathways is to quantify the effectiveness of negative emissions technologies which will be strongly affected by carbon cycle feedbacks. Here we find that Earth system models suggest significant weakening, even potential reversal, of the ocean and land sinks under future low emission scenarios. For the RCP2.6 concentration pathway, models project land and ocean sinks to weaken to 0.8 ± 0.9 and 1.1 ± 0.3 GtC yr−1 respectively for the second half of the 21st century and to −0.4 ± 0.4 and 0.1 ± 0.2 GtC yr−1 respectively for the second half of the 23rd century. Weakening of natural carbon sinks will hinder the effectiveness of negative emissions technologies and therefore increase their required deployment to achieve a given climate stabilisation target. We introduce a new metric, the perturbation airborne fraction, to measure and assess the effectiveness of negative emissions.

2012, Souty, F., Brunelle, T., Dumas, P., Dorin, B., Ciais, P., Crassous, R., Müller, C., Bondeau, A.

Interactions between food demand, biomass energy and forest preservation are driving both food prices and land-use changes, regionally and globally. This study presents a new model called Nexus Land-Use version 1.0 which describes these interactions through a generic representation of agricultural intensification mechanisms within agricultural lands. The Nexus Land-Use model equations combine biophysics and economics into a single coherent framework to calculate crop yields, food prices, and resulting pasture and cropland areas within 12 regions inter-connected with each other by international trade. The representation of cropland and livestock production systems in each region relies on three components: (i) a biomass production function derived from the crop yield response function to inputs such as industrial fertilisers; (ii) a detailed representation of the livestock production system subdivided into an intensive and an extensive component, and (iii) a spatially explicit distribution of potential (maximal) crop yields prescribed from the Lund-Postdam-Jena global vegetation model for managed Land (LPJmL). The economic principles governing decisions about land-use and intensification are adapted from the Ricardian rent theory, assuming cost minimisation for farmers. In contrast to the other land-use models linking economy and biophysics, crops are aggregated as a representative product in calories and intensification for the representative crop is a non-linear function of chemical inputs. The model equations and parameter values are first described in details. Then, idealised scenarios exploring the impact of forest preservation policies or rising energy price on agricultural intensification are described, and their impacts on pasture and cropland areas are investigated.

2015, Nishina, K., Ito, A., Falloon, P., Friend, A.D., Beerling, D.J., Ciais, P., Clark, D.B., Kahana, R., Kato, E., Lucht, W., Lomas, M., Pavlick, R., Schaphoff, S., Warszawaski, L., Yokohata, T.

We examined the changes to global net primary production (NPP), vegetation biomass carbon (VegC), and soil organic carbon (SOC) estimated by six global vegetation models (GVMs) obtained from the Inter-Sectoral Impact Model Intercomparison Project. Simulation results were obtained using five global climate models (GCMs) forced with four representative concentration pathway (RCP) scenarios. To clarify which component (i.e., emission scenarios, climate projections, or global vegetation models) contributes the most to uncertainties in projected global terrestrial C cycling by 2100, analysis of variance (ANOVA) and wavelet clustering were applied to 70 projected simulation sets. At the end of the simulation period, changes from the year 2000 in all three variables varied considerably from net negative to positive values. ANOVA revealed that the main sources of uncertainty are different among variables and depend on the projection period. We determined that in the global VegC and SOC projections, GVMs are the main influence on uncertainties (60 % and 90 %, respectively) rather than climate-driving scenarios (RCPs and GCMs). Moreover, the divergence of changes in vegetation carbon residence times is dominated by GVM uncertainty, particularly in the latter half of the 21st century. In addition, we found that the contribution of each uncertainty source is spatiotemporally heterogeneous and it differs among the GVM variables. The dominant uncertainty source for changes in NPP and VegC varies along the climatic gradient. The contribution of GVM to the uncertainty decreases as the climate division becomes cooler (from ca. 80 % in the equatorial division to 40 % in the snow division). Our results suggest that to assess climate change impacts on global ecosystem C cycling among each RCP scenario, the long-term C dynamics within the ecosystems (i.e., vegetation turnover and soil decomposition) are more critical factors than photosynthetic processes. The different trends in the contribution of uncertainty sources in each variable among climate divisions indicate that improvement of GVMs based on climate division or biome type will be effective. On the other hand, in dry regions, GCMs are the dominant uncertainty source in climate impact assessments of vegetation and soil C dynamics.

2015, Yue, C., Ciais, P., Cadule, P., Thonicke, K., van Leeuwen, T.T.

Carbon dioxide emissions from wild and anthropogenic fires return the carbon absorbed by plants to the atmosphere, and decrease the sequestration of carbon by land ecosystems. Future climate warming will likely increase the frequency of fire-triggering drought, so that the future terrestrial carbon uptake will depend on how fires respond to altered climate variation. In this study, we modelled the role of fires in the global terrestrial carbon balance for 1901–2012, using the ORCHIDEE global vegetation model equipped with the SPITFIRE model. We conducted two simulations with and without the fire module being activated, using a static land cover. The simulated global fire carbon emissions for 1997–2009 are 2.1 Pg C yr−1, which is close to the 2.0 Pg C yr−1 as estimated by GFED3.1. The simulated land carbon uptake after accounting for emissions for 2003–2012 is 3.1 Pg C yr−1, which is within the uncertainty of the residual carbon sink estimation (2.8 ± 0.8 Pg C yr−1). Fires are found to reduce the terrestrial carbon uptake by 0.32 Pg C yr−1 over 1901–2012, or 20% of the total carbon sink in a world without fire. The fire-induced land sink reduction (SRfire) is significantly correlated with climate variability, with larger sink reduction occurring in warm and dry years, in particular during El Niño events. Our results suggest a "fire respiration partial compensation". During the 10 lowest SRfire years (SRfire = 0.17 Pg C yr−1), fires mainly compensate for the heterotrophic respiration that would occur in a world without fire. By contrast, during the 10 highest SRfire fire years (SRfire = 0.49 Pg C yr−1), fire emissions far exceed their respiration partial compensation and create a larger reduction in terrestrial carbon uptake. Our findings have important implications for the future role of fires in the terrestrial carbon balance, because the capacity of terrestrial ecosystems to sequester carbon will be diminished by future climate change characterized by increased frequency of droughts and extreme El Niño events.

2014, Yue, C., Ciais, P., Cadule, P., Thonicke, K., Archibald, S., Poulter, B., Hao, W.M., Hantson, S., Mouillot, F., Friedlingstein, P., Maignan, F., Viovy, N.

Fire is an important global ecological process that influences the distribution of biomes, with consequences for carbon, water, and energy budgets. Therefore it is impossible to appropriately model the history and future of the terrestrial ecosystems and the climate system without including fire. This study incorporates the process-based prognostic fire module SPITFIRE into the global vegetation model ORCHIDEE, which was then used to simulate burned area over the 20th century. Special attention was paid to the evaluation of other fire regime indicators such as seasonality, fire size and fire length, next to burned area. For 2001–2006, the simulated global spatial extent of fire agrees well with that given by satellite-derived burned area data sets (L3JRC, GLOBCARBON, GFED3.1), and 76–92% of the global burned area is simulated as collocated between the model and observation, depending on which data set is used for comparison. The simulated global mean annual burned area is 346 Mha yr−1, which falls within the range of 287–384 Mha yr−1 as given by the three observation data sets; and is close to the 344 Mha yr−1 by the GFED3.1 data when crop fires are excluded. The simulated long-term trend and variation of burned area agree best with the observation data in regions where fire is mainly driven by climate variation, such as boreal Russia (1930–2009), along with Canada and US Alaska (1950–2009). At the global scale, the simulated decadal fire variation over the 20th century is only in moderate agreement with the historical reconstruction, possibly because of the uncertainties of past estimates, and because land-use change fires and fire suppression are not explicitly included in the model. Over the globe, the size of large fires (the 95th quantile fire size) is underestimated by the model for the regions of high fire frequency, compared with fire patch data as reconstructed from MODIS 500 m burned area data. Two case studies of fire size distribution in Canada and US Alaska, and southern Africa indicate that both number and size of large fires are underestimated, which could be related with short fire patch length and low daily fire size. Future efforts should be directed towards building consistent spatial observation data sets for key parameters of the model in order to constrain the model error at each key step of the fire modelling.