Historical and future perspectives of global soil carbon response to climate and land-use changes

In this paper, we attempt to analyse the respective influences of land-use and climate changes on the global and regional balances of soil organic carbon (SOC) stocks. Two time periods are analysed: the historical period 1901–2000 and the period 2000–2100. The historical period is analysed using a synthesis of published data as well as new global and regional model simulations, and the future is analysed using models only. Historical land cover changes have resulted globally in SOC release into the atmosphere. This human induced SOC decrease was nearly balanced by the net SOC increase due to higher CO2 and rainfall. Mechanization of agriculture after the 1950s has accelerated SOC losses in croplands, whereas development of carbon-sequestering practices over the past decades may have limited SOC loss from arable soils. In some regions (Europe, China and USA), croplands are currently estimated to be either a small C sink or a small source, but not a large source of CO2 to the atmosphere. In the future, according to terrestrial biosphere and climate models projections, both climate and land cover changes might cause a net SOC loss, particularly in tropical regions. The timing, magnitude, and regional distribution of future SOC changes are all highly uncertain. Reducing this uncertainty requires improving future anthropogenic CO2 emissions and land-use scenarios and better understanding of biogeochemical processes that control SOC turnover, for both managed and un-managed ecosystems.


Introduction
During the past decades, research efforts have been made to improve our knowledge about the effects of climate, atmospheric CO 2 and land-use change on terrestrial C cycle. However, the mechanisms and factors that govern global uptake and release of C from the terrestrial reservoir, and their regional importance, are still poorly quantified. This could only be achieved through C cycle modelling studies as long-term in situ measurements are very sparse and lack global coverage, and remote sensing techniques have limited capability for estimating below canopy processes such as soil respiration and are subject to significant absolute errors in carbon fluxes (Valentini et al., 2000;Baldocchi, 2008). There exist several C cycle modelling studies for understanding and predicting the global terrestrial carbon cycle (McGuire et al., 2001;Cao et al., 2002;Peylin et al., 2005;Zeng et al., 2005;Friedlingstein et al., 2006;Sitch et al., 2008;Piao et al., 2009a,b). However, in these studies, most models either did not incorporate land use changes, or did not separate the contributions of land use change to the global and regional C balances from the contributions of change in CO 2 , temperature and precipitation.
The goal of this study is to analyse implications of changes in climate and land use for the global and regional balance of soil organic carbon (SOC). Two contiguous time period are analysed, the historical period 1901-2000 and the future 2000-2100. Models are used for both periods, and observations for the historical period. We seek answers to the following questions: • What was the impact of rising CO 2 and changing climate on the global SOC budget during the 20th century?
• What was the impact of land cover change and agricultural management change on SOC during the 20th century?
• Which regions have lost or gained SOC in response to land cover and climate change during the 20th century?
• How will the SOC balance change in the 21st century in response to land cover, agricultural management and to climate change?
• What are the uncertainties?
A presentation of the SOC global distribution for large biomes is given in Section 2. The global SOC vulnerability to climate change is presented in Section 3. Historical changes in SOC induced by climate, CO 2 and land cover change are analysed for different regions using a model in Section 4, and future projections in Section 5. We discuss uncertainties in Section 6.

Global distribution of SOC
The global mass of SOC is at least ∼2300 PgC (Jobbágy and Jackson, 2000). It is ∼75% of the total land biosphere C reservoir given that litter and vegetation C stocks amount to 850 PgC (Houghton, 2007). SOC comprises a wide range of compounds, with different age, physical and chemical properties. SOC is formed by input from litter and marginally from charcoal formation after fire (0.005-0.27 PgC yr −1 according to Kuhlbusch and Crutzen, 1995). SOC is removed from the soil by microbial decomposition (heterotrophic respiration), by river export of dissolved organic carbon, by ecosystem disturbance such as fire and by wind and water erosion. Erosion acts generally on longer timescales than the former processes. Table 1 summarizes the SOC distribution among the main biomes. The distribution of C in the soil compared to total ecosystem C stock varies greatly among biomes (from two-third to almost 100%). SOC mass is divided between 'modern' SOC pools summing up to 2300 PgC (Jobbágy and Jackson, 2000), and old (typically >5000 yr) frozen SOC and peat SOC pools. Old SOC pools are mainly distributed in northern permafrost regions. The mass of frozen and peat SOC is huge, as recently appraised (Schuur et al., 2008;Tarnocai et al., 2009). In these regions, total SOC amounts to ∼1600 PgC with 1024 PgC in the top 3 m (Tarnocai et al., 2009). Old SOC pools were formed during the Pleistocene in the case of Yedoma deposits in Eastern Siberia (Zimov et al., 2006), or during the early Holocene for peat deposits in the northern high latitudes (Beilman et al., 2009). This large frozen pool is equivalent to 70% of the 'conventional' C stocks. Finally, extensive peat deposits exist both in permafrost and non-permafrost regions, where decomposition is limited by temperature and/or oxygen availability in saturated organic soils. The CO 2 and CH 4 balance of peatlands is thus acutely sensitive to changes in both hydrology and temperature. Hence, there is a concern that these old SOC pools are vulnerable to future warming at high latitude, as will be seen in the next section.
In the other terrestrial biomes presented in Table 1, SOC stock is determined by the interaction between factors limiting decomposition and C inputs from plants. SOC stocks distribution between these biomes is therefore not obvious. Generally, the mean residence time of SOC, defined as the ratio of the SOC stock to heterotrophic respiration loss (Malhi et al., 2002), is longer at high latitudes for boreal forest and tundra (Bird et al., 1996;  Data were gathered from Jobbagy and Jackson (2000), IPCC (2000), Davidson and Janssens (2006) and Tarnocai et al. (2009). a Jobbagy and Jackson (2000) except for peatlands (IPCC, 2000) and permafrost (Tarnocai et al., 2009). WBGU (German Advisory Council on Global Change, 1988in IPCC, 2000 and MRS (Mooney, Roy andSaugier, 2001 in IPCC, 2000) estimates are given in parentheses. b Jobbagy and Jackson (2000) except for peatlands (IPCC, 2000) and permafrost (Tarnocai et al., 2009). WBGU (German Advisory Council on Global Change, 1988in IPCC, 2000 and IGBP (International Atmosphere-Biosphere Program in IPCC, 2000) estimates are given. c Calculated from WBGU and MRS estimates (IPCC, 2000). d Calculated from Jobbagy and Jackson (2000) except for peatlands and permafrost (Tarnocai et al. 2009). WBGU estimates (IPCC, 2000) are given in parentheses. e Permafrost partly includes peatlands, boreal forests and boreal grasslands. Trumbore, 2000). Primary production is the dominant process determining SOC stocks in tropical forests (Zhao et al., 2004). According to Table 1, the fraction of SOC to total ecosystem C represents roughly 60-80% in forests and more than 80% in grasslands. In tundra soils and peatland and permafrost C deposits, SOC represents nearly 100% of total ecosystem C stocks. It should be noted that SOC stocks estimates presented in Table 1 are generally higher than those given in the IPCC LULUCF report for the top 1 m of soil (IPCC, 2000) because they were estimated down to 3 m (Jobbagy and Jackson, 2000). This difference is particularly apparent for tropical forests (about +480 PgC).

Global vulnerability of SOC to climate
The land biospheric reservoir (soil and living biomass) constitutes a net sink of atmospheric CO 2 at global scale. The size of this sink has increased alongside with increasing emissions, from 1.8 PgC in the 1980s to 2.6 PgC in the 1990s and 3 PgC over (Le Quéré et al., 2009. It is currently considered that land biosphere carbon sink absorbs about 30% of the anthropogenic emissions from fossil fuel combustion and deforestation every year (Canadell et al., 2007;IPCC, 2007). This sink is explained by the increased net primary productivity (NPP) driven by a variety of environmental and management factors (Nemani et al., 2003;Piao et al., 2005;Magnani et al., 2007;Ciais et al., 2008;Lewis et al., 2009). Among the environmental factors, rising CO 2 , warming in 'temperature-limited' ecosystems, increased rainfall in dry ecosystems, extra-nitrogen deposition in 'nitrogen-limited' ecosystems and the combination of these factors are the main causes of increased NPP (Piao et al., 2006). This increase in NPP induces higher soil C stock by increasing C input to soil. However, although soil contains three times more C than the vegetation, it is a two times smaller sink than living biomass (Reichstein, 2008). There are two reasons for that. First, in managed ecosystems, biomass is partly harvested and consequently the soil does not 'benefit' from all the increased primary productivity. Second, in some ecosystems, increased primary productivity adds C into biomass pools where it could be stored for decades, and the soil has therefore not 'seen' the increase in primary productivity yet. This is especially true for the woody biomass pools of forests. Yet, global increase of soil C stock is mostly due to increased foliage and fine root productivity, which are neither harvested nor long lasting biomass pools (Nemani et al., 2003;Piao et al., 2006). However, the ability of soil to behave as a C sink in the future is highly uncertain because SOC decomposition is sensitive to climate and increase with temperature. Generally, soil C respiration rate is multiplied by roughly a factor of two for a 10 • C warming (Raich and Schlesinger, 1992;Davidson and Janssens, 2006;Reichstein et al., 2005) if microbial activity is not limited by substrate availability or soil moisture. The potential loss of soil C under the future warming forecasted by climate models is six times larger than the current soil C sink (Reichstein, 2008). Large uncertainties are associated with this estimate, both in the amount and spatial distribution of warming which depends on economic development scenarios and modelled climate processes, and in the response of soil C to climate change.
In the future, the net balance between extra inputs and extra losses is difficult to assess. Indeed, neither SOC input nor output fluxes are linear function of climate, and these functions are uncertain. It is generally assumed that decomposition increases exponentially with temperature (Lloyd and Taylor, 1994) and productivity is a saturating function of climate. Conceptually, this means that, above a certain degree of warming, the increase of respiration will not be offset by an equivalent productivity increase, and soils will lose C. Moreover, rising atmospheric CO 2 concentration and warming may no more induce productivity increase in ecosystems where nutrient availability is limited and/or where drought accompanies the warming (Angert et al., 2005). This latter response is observed in simulations of the climate-carbon coupled system. For instance, in the widely cited study of Cox et al. (2000), Amazon forest productivity drops, causing a large and uncontrolled loss of soil C in the future century, acting as a positive feedback to global warming. Contrary to the results of this study, most global models of the coupled climate-SOC system predict a net gain of SOC during the 21st century even in the tropics, except for two models including the one used by Cox et al. (see Friedlingstein et al., 2006).
Moreover, in face of transient climate change, an ecosystem may lose soil C rapidly in response to warming while its productivity increases slowly over time. In that case, although the asymptotic value of soil C is theoretically proportional to the new productivity and may thus eventually increase above initial levels, a transient phase of soil C loss lasting several decades, followed by recovery and soil C build-up can be produced.
For soil layers where the C input are virtually zero, such as subsurface peat and permafrost C, the process of decomposition dominates the vulnerability of soil C (Zimov et al., 2006;Khvorostyanov et al., 2008a,b). Permafrost-affected soils span a wide range of soil environments, and thus the huge 'permafrost soil carbon pool' includes carbon stocks which have built up over a variety of timescales and with differing levels of vulnerability to global warming ( Fig. 1). At their surface, permafrostaffected soils have a seasonally thawing active layer, in which SOC balance is, as in temperate soils, sensitive to both litter inputs and respiration fluxes over short timescales. However, the defining feature of permafrost soils is that below the surface there exist soil layers that do not thaw, and in these layers carbon may accumulate for long periods of time, leading to the massive quantities of carbon stored there. Shallow permafrost layers often contain carbon that is transported from the active Fig. 1. Schematic of permafrost carbon cycle. Shallow permafrost carbon is likely to be the most vulnerable to deepening of active layers with warming, whereas deeper permafrost carbon may be vulnerable as well due to processes such as thermokarst or microbial heating. Roots are not figured but may intervene in both active layer and shallow permafrost. layer through cryoturbation, and including this process in terrestrial models leads to a large increase in modelled C storage at high latitudes (Koven et al., 2009). The uncertainties of both the magnitude and vulnerability of these carbon pools increase with depth. Shallow permafrost carbon is clearly highly vulnerable, and field evidence suggests that carbon loss due to thawing of shallow permafrost is already taking place in Alaska (e.g. Schuur et al., 2008). The uncertainty on the quantity of carbon in deeper permafrost stocks, such as Yedoma, is much higher, as is the uncertainty on its vulnerability. Indeed, the thawing of deep permafrost within the next 100-200 years would require additional processes such as expansion of thermokarst lakes formed by meltwater (Walter et al., 2007) or release of heat by soil microorganisms (Khvorostyanov et al., 2008a,b), which are only beginning to be added to models. It is nonetheless interesting to note that most global carbon cycle models  show a net gain in high-latitude terrestrial carbon due to warming, while initial results using a permafrost-enabled version of ORCHIDEE (Koven C., Ringeval B., Ciais P., Friedlingstein P., personal communication) show carbon losses due to warming.

Global vulnerability of SOC to land-use changes
Land cover and land management also strongly impact soil C balance acting on both decomposition and input. Land use changes that began with human agriculture several millennia ago, have escalated during the last centuries. These have consisted mostly of expanding cropland and pastureland and, more recently of urbanization. Today land use change is concentrated in the tropics, where forest clearing is taking place (Houghton, 2003). Tropical forest clearing induced a direct loss of C in biomass of 1.4 PgC yr −1 over 2000-2008, but clearing also causes an indirect loss of soil C (Houghton, 2003) that may last for decades after land use change (Pongratz et al., 2009). The fate of soil C in newly deforested areas converted to cropland or pastureland will depend on agricultural practice. Intensive cropping will quickly empty the former soil C pools within two to three decades (Arrouays et al., 1995;Jolivet et al., 1997;Reeves et al., 1997), whereas sustainable pasture management may stabilize them to a value close to the former tropical forest (Trumbore et al., 1995). In the following, we thus examine the effects of land use change on soil C pools during the last century, and into the future.
To appreciate the timescales involved in the evolution of soil C balance after a land use change, Fig. 2 shows the measured evolution of soil C in so-called 'bare fallow' field experiments. These soils, formerly under cropland or grassland were kept free of vegetation by human intervention during several decades, while SOC decomposition continues. Figure 2 shows that even after 40 years of bare soil conditions, the soil is continuing to lose C, apparently without reaching a new equilibrium. The shape of SOC decay in bare fallow plots cannot be reproduced with a simple 1-pool exponential decay model. Moreover, the SOC decay differs according to the former land use (Fig. 2). For example, the Rothamsted bare fallow implemented on grassland shows an initial large C loss that is not observed in the Askov bare fallow implemented on arable land. This suggests that soil organic matter is composed of two or more pools, with some pools having a decay time of 50-100 years, or more. Moreover, the size of the more stable pool may depend on previous land-use (Balesdent et al., 1988).
Consequently, any change in land use or ecosystem disturbance is expected to have a long lasting effect on the soil C balance. Former land use will continue to exert an influence on today's soil C budget and pre-condition its future evolution as well. Because of the cascade of long timescales involved, one needs to take a long historical perspective to simulate land use and climate change impacts on soil C balance.

Historical land cover, agricultural management and climate change effects on SOC
To gain insights on past land-use change effects on SOC regionally, the results of a global biosphere model forced by land-use and climate changes since 1900 are analysed. This model is a process-based global ecosystem model called ORCHIDEE (ORganizing Carbon and Hydrology In Dynamic EcosystEms model; Krinner et al., 2005). ORCHIDEE was used to quantify the effects of changes in climate, atmospheric CO 2 concentration and changes in land cover area on terrestrial C cycle for the period from 1901 to 2002. SOC stocks were considered to be at steady-state equilibrium in 1901. This is a very crude assumption because human-induced land-use change has occurred before the 20th century particularly in Europe, in the Former Soviet Union and in East Asia (Pongratz et al., 2008). This human-induced land-use change is not negligible because it was estimated that about 5 millions of km 2 were transformed to cropland or pasture between 800 and 1700 AD, leading SOC stocks to a significant disequilibrium in 1900 (Pongratz et al., 2009). For example, they estimated that, whereas vegetation lost about 100 Gt C from 800 to 1900, global SOC stocks increased by about 30 PgC due to additional plant material added to the soil pools from the converted natural vegetation. The steady-state assumption underlying the ORCHIDEE simulation results analysed below implies that the effect of land-cover change on SOC stocks is likely to be overestimated. More details about the model set up and architecture are given in Appendix.

Historical global land cover change and SOC balance
Both land cover and CO 2 + climate changes have a strong impact in SOC balance. Running simulations taking into ac-count CO 2 + climate changes only or CO 2 + climate + land cover changes allow determining the relative importance of land cover change. The results of these two simulations at global scale are presented on Fig. 3a. If only CO 2 and climate drivers are accounted for, NPP and thereby biomass increases ( Fig. 3a). Soil C shows a slight decreasing trend until 1960, due to warming induced acceleration of SOC decomposition in Europe and in temperate North America. Elsewhere, the CO 2 fertilization effect on NPP incorporated in the model increases soil C input in excess to warming increased decomposition. Therefore, the soil C balance is slightly positive over the 20th century. By contrast, when annual land cover changes are incorporated in the study, the soil C stock decreases globally between 1900 and 1960. This land use induced decrease is later reversed, ending up only in a small variation of soil C, SOC = −7.3 PgC by year 2000. The net change represents an insignificant 0.3% decrease of the initial stocks in 1901, assumed in equilibrium. The essential point here is that land-use change has an effect on soil C change between 1901 and 2000 that is equivalent in magnitude to the effect of climate and CO 2 change, but of opposite sign. These results are in agreement with the results of a modelling study by McGuire et al. (2001) who did Fig. 3. Evolution (a) of global soil and vegetation C stocks simulated by ORCHIDEE-LUC in response to combined effects of atmospheric CO 2 concentration and climate, and to combined effects of atmospheric CO 2 concentration, climate and land cover change and (b) of global soil C stocks simulated by the eleven models of the C4MIP intercomparison (Friedlingstein et al., ) from 1901(Friedlingstein et al., to 2002 not distinguish soil carbon pools from the total land biospheric reservoir. Hence, McGuire et al. (2001) estimate with four terrestrial biosphere models that terrestrial ecosystems transitioned during the 20th century from releasing carbon to the atmosphere (∼8.8 PgC from 1920 to 1957) because of land-use changes to storing carbon (∼14.3 PgC from 1958 to 1992) mainly because of the increase in atmospheric CO 2 concentration.
Because the change in NPP consecutive to forest clearing is not as large as the direct loss of biomass ( VEG), the impact of historical land cover change is more drastic on the biomass pools change than on the SOC pools change. Therefore, in our simulation, compared to the effect of CO 2 and climate alone, including land cover change modifies VEG globally from + 140 to −51 PgC and SOC only from +27 to −7.3 PgC. This land cover effect on global SOC stocks (−34.3 PgC) is almost half the range of global SOC changes simulated by the eleven C4MIP models for the 20th century and represents ∼1/3 (MPI) to ∼100 times (UMD) the estimated effect of climate and atmospheric CO 2 changes (Figs 3a and b). Moreover, ORCHIDEE-LUC shows the only negative SOC balance. This result highlights (1) the importance of land-use change comparatively to climate change and (2) the need for consideration of past land-use change in models that were parameterized to be consistent with the current contemporary global land carbon budget. On the other hand, the model does not reproduce secondary forest ecosystem growth unless a given grid point is subject to deforestation and later to reforestation. Thus, the negative impact of deforestation on the SOC balance is probably overestimated (Shevliakova et al., 2009).
When looking at the combined effects of land use, climate and CO 2 changes on the SOC balance at regional scale, we can see in Fig. 4a different behaviour between tropical regions on the Fig. 4. Evolution of regional soil and vegetation C stocks simulated by ORCHIDEE-LUC in response to combined effects of atmospheric CO 2 concentration and climate, and to combined effects of atmospheric CO 2 concentration, climate and land cover change from 1901 to 2002. one hand, and temperate and boreal regions on the other hand. In temperate and boreal regions, the effect of land cover change on the SOC stocks is larger than in the tropics, even though the loss of biomass is greater in the tropics.
More specifically, in the Former Soviet Union (FSU) and North America, forest clearing took place during the first half of the 20th century. This forest clearing has impacted the SOC balance with inertia, owing to the long residence times of SOC (Figs 4c and d). In the FSU, the modelled difference between SOC stocks with climate and CO 2 changes only (S1) and SOC stocks with land-use changes (S2) is still increasing with time. In North America, this difference is rather constant after 1970, indicating that the indirect effect of land cover change on SOC stocks is stabilised. In Europe and China (Figs 4e and f), no major change in forest area took place during the 20th century (most of deforestation being earlier to that date). The SOC difference between S2 and S1 is thus negligible, indicating that land cover change did not impact SOC stocks. However, we will see below that SOC stocks in China and Europe has likely declined during the 20th Century due to intensification of agriculture. In China, large plantations created in the 1980s may contribute to an increase in SOC (Piao et al., 2009a,b), which is not well captured by the coarse scale global land cover maps used to force our model.
In our model simulation, climate and CO 2 have increased SOC stocks everywhere but in Africa (because of drought that decreased NPP, see Fig. 4h) and in Canada (because increased NPP did not compensate warming increased heterotrophic respiration, see Fig. 4g). In contrast, land cover change has diminished SOC stocks everywhere. On average, the effect of land cover change on SOC balance between 1901 and 2000 is opposite to the one of climate and CO 2 . Regions where land cover induced SOC losses dominate over climate and CO 2 -induced SOC accumulation are the USA, the FSU and Europe. These regions experienced a net loss of SOC since 1901. Oppositely, a net SOC gain since 1901 was modelled for tropical Asia, China and South America. Despite deforestation in tropical Asia and South America, the increase in NPP induced by higher CO 2 and rainfall in undisturbed forests (Lewis et al., 2009) has driven the pan-tropical SOC balance positive up to the present.

Agricultural intensity and SOC balance
On top of climate, CO 2 and land cover, land management intensity also impacts SOC stocks, especially in cultivated soils. Several studies reported that changes in the fate of cropping residues (straw and stubble), cultivation and tillage practices have a strong impact on SOC stocks (Paustian et al., 1998(Paustian et al., , 2000West and Post, 2002). In addition, arable soil erosion (which depends on SOC content, cultivation and tillage practices) causes a significant loss of SOC in cultivated soils. However, at globalscale, accumulation of eroded soil particles in floodplains may compensate for this loss at the field scale (van Oost et al., 2007).
Changes in agriculture that took place in North America, China and Europe after 1950 have likely accelerated the SOC losses (Qiu et al., 2009). This has been illustrated by numerous site measurements (Burke et al., 1989;Paul et al., 1997), but few measurements at regional or continental scale (Mäkelä-Kurtto and Sippola, 2002;Sleutel et al., 2003;Bellamy et al., 2005).
In North America, this trend may have been reversed recently, as it has been considered that changes in agricultural practices have increased SOC sequestration during the late 30 years. Indeed, Ogle et al. (2003) estimated that adoption of reduced use of bare fallow, conservation tillage practices on cropland, and conversion of annual cropland to grass and trees in the Conservation Reserve Program (CRP) have resulted in a net gain of 10.8 Tg C yr −1 during the period from 1982 to 1997. Most of this gain was due to setting-aside lands in the CRP and only 1.3 Tg C yr −1 was attributed to change in US agricultural land-use and management. Using the ISAM-2 model, Jain and Yang (2005) estimated a 868 Tg SOC increase (i.e. 43.4 Tg C yr −1 from 1980 to 2000) in North American soils due to no-till practices. Productivity increases induced by irrigation may also have played a significant role in SOC stock increase in the USA. The SOC stock increase induced by irrigation was estimated at 21.3 Tg C in the Great Plains since 1950s using the CENTURY model (Parton et al., 2005).
In China, a continuous decline in SOC levels for a wide scope of Chinese agricultural soils has been observed since the 1950s (Qiu et al., 2009). This was mainly explained by the reduction of crop residue incorporation and manure amendments, which have been the major source of SOC for most Chinese farmlands (Qiu et al., 2009). As in North America, this trend seems to have been reversed recently. New policies have been launched during the past decades to encourage farmers to return more organic matter to the soil and some agricultural fields have shown SOC sequestration (Huang and Sun, 2006;Yan et al., 2007). In addition, intensive agriculture practices, such as application of chemical fertilizer, increased irrigation in arid areas, expansion of straw incorporation, and shallow plowing have also led to an increase in the C sequestration of agricultural ecosystems (Huang and Sun, 2006). Preliminary estimates suggested that along with the increase of agricultural productivity and less removal of biogenic material for societal energy purposes, about 0.025-0.037 PgC yr −1 is potentially fixed in soil (Lal, 2004), which is slightly higher than current estimation of Huang and Sun (2006) and Xie et al. (2007). Huang and Sun (2006) analysed changes of organic carbon stocks in China's agricultural soils during the past two decades, and conclude that the mean carbon sink of agricultural soils was 0.015-0.020 PgC yr −1 . Overall, Piao et al. (2009a,b) estimated in a meta-analysis that topsoil in croplands is a net sink of 0.026 ± 0.011 PgC yr −1 since the early 1980s.
In Europe, SOC change was estimated using three models including specific crop management and practice ORCHIDEE-STICS, RothC, and LPJml. In both ORCHIDEE-STICS and RothC, the simulated soil C dynamics included the effect of past changes in agricultural technology on input and tillage, as well as the effects of climate and atmospheric CO 2 changes. The ORCHIDEE-STICS model was initialized with (reconstructed) ancestral farming practice and crop varieties in 1900 (Gervois et al., 2008). This model provides a net change in SOC, SOC = 0.01 ± 0.06 tC ha −1 yr −1 between 1901 and 2000. This insignificant gain is mostly explained by agricultural intensification and agricultural technology changes, with a small negative contribution of recent climate change over Southern Europe (Gervois et al., 2008). Indeed, in the Iberian Peninsula and a few southern Mediterranean regions, the increase in heterotrophic respiration and a precipitation-induced reduction in soil C inputs slightly overcome positive SOC stock increment induced by agricultural technology changes. At the European scale, a strong loss of SOC is modelled when agriculture became mechanized in the 1950s. But this loss trend was reversed in the 1970s and SOC increased thereafter. Arable lands are today C neutral (+0.16 tC ha −1 yr −1 with a range of 0-0.3 tC ha −1 yr −1 ). This modelled small sink is not supported by regional inventories, available in UK, France, Belgium and Finland, and in some regions of Germany, which rather suggest net SOC losses (Mäkelä-Kurtto and Sippola, 2002;Sleutel et al., 2003;Bellamy et al., 2005;Smith et al., 2005a,b). However, if the net C sink provided by erosion (Lal, 2003;van Oost et al., 2007) is taken into account, SOC losses measured in the inventory become comparable to the model results. This confirms small net changes in the SOC balance of arable lands in Europe.
Over 1990-1999, RothC and LPJml provide Net Biome Productivities (NBP) of −7.6 and 1.3 g C m −2 yr −1 , respectively. It is noteworthy that these results are extremely dependent on assumptions about the management options employed and how these have changed over recent decades. In summary, the three process models integrated over Europe and the inventories available at regional-scale agree on the relatively small absolute magnitude of the net SOC balance of arable lands in Europe. However, each model's sensitivity tests point to the strong sensitivity of NBP to the assumed choice and past history of management options.

Forest management and SOC balance
The impact of forest management change on soil carbon has not been assessed for the past century. ORCHIDEE-FM, a novel version of ORCHIDEE includes a forest management module. Recent results emphasized the importance of forest management on SOC stocks (Bellassen et al., 2010a,b). The evolution of carbon stocks in a managed forest, as simulated by ORCHIDEE-FM, is reported on Fig. 5a. Harvesting wood increases soil carbon in the short term as 50-100 tC ha −1 of slash wood is added to the litter pool as coarse woody debris. This input leads to a peak in soil carbon respiration of about +50% after harvest, which is consistent with measurements (Savage and Davidson, 2001;Pypker and Fredeen, 2002;Pregitzer and Euskirchen, 2004). During decomposition, a part of the coarse woody debris is buried in deep soil horizons, driving a smaller and slower increase in soil carbon per se (Nabuurs et al., 2008). However, in the long term the impact of forest management on SOC stock is negative: as the stems of dead trees are no longer left to feed soil carbon, the long-term litter input is lower under management than in an unmanaged forest. On average, ORCHIDEE-FM simulates a soil carbon stock reduction of about 15 tC ha −1 under management (Fig. 5a) compared to no management (Fig. 5b).

Future changes in climate and land use, and the SOC balance
At this stage, few integrated global model projections combine both land use change and climate change impacts on the SOC balance during the 21st century. We discuss below results from coupled climate-carbon C4MIP models from the study of Friedlingstein et al. (2006). This study only allowed the effects of change in climate and CO 2 to be examined because it did not incorporate land cover change. Hence, results from a coarse resolution terrestrial C model OSCAR (Gitz and Ciais, 2004) that incorporates both land use change and warming (∼ +3.5 • C from 2000 to 2100) are analysed in counterpoint.

Effect of climate and CO 2 change
The modelling study by Cox et al. (2000) has raised a great concern on the fact that the SOC balance is sensitive to climate via both changes in input and in decomposition. In their study, the drying of the Amazon region leads to forest dieback in the Triffid vegetation model. Forest dieback is accompanied by decreased NPP and subsequently by a large SOC loss of about 170 PgC between 2000 and 2100. This SOC loss in tropical forests increases CO 2 and climate warming by the end of the 21st century. A more detailed study by Jones et al. (2005) using the same climate model and two different soil C models (RothC and Triffid) investigated the opposite role of climate through changed input and changed decomposition. The effect of temperature increase alone causes a decrease of SOC in northern regions during the 21st century, such as Siberia and Canada. The effect of increased NPP is opposite and increased SOC. The two different C models show both a loss of SOC stocks but differ by 26 PgC, due to different climate sensitivity. CO 2 -and climate-induced C input increase would drive an increase in SOC stocks everywhere by 2100, except in the Amazon region. Increase in C input and decomposition would expand the area of regions that will lose SOC to most of the tropics and mid-latitude biomes, except for high latitudes where a small SOC increase is predicted. This illustrates the need to consider future change in the SOC balance as being determined not only by the vulnerability of SOM to warming but also by increased or decreased input to the soil, including changes in litter quality (Fierer et al., 2005).
The current uncertainties on the future global SOC balance can be quantified using data from intercomparisons of many coupled climate-carbon cycle models. Uncertainties related to climate and SOC decay equations can be investigated using the results from C4MIP model intercomparison . Figure 6 presents the SOC changes simulated by 11 models from 2000 and 2099 at both global and regional scales. One common scenario of CO 2 anthropogenic emissions (Special Report on Emission Scenarios (SRES) A2 scenario) was considered. Although the models have been parameterized to be consistent with the current contemporary global land carbon budget, their responses differ markedly. The predicted change in global SOC stocks ranges between −45.8 PgC (HadCM3LC) and 310 PgC (MPI). In most models, soil is projected to gain C during the 21st century since C input driven by enhanced vegetation productivity is predicted to increase more than SOC decomposition. This is true for most models even in the tropics (30 • S-30 • N). However, for 2 of the 11 models (HadCM3LC used by Cox et al. and UMD), soils in tropical and temperate (30 • N-60 • N) regions will lose C leading to global C release from the soil into the atmosphere. Moreover, three other models (BERN-CC, FRCGC and Uvic-2,7) show a decrease in the rate of SOC increase during the last 20 years suggesting that they would simulate SOC release during the 22nd century if atmospheric CO 2 concentration and temperature still increase. This large spread in the prediction of SOC change can be partly explained by differences in model structures (number of SOC pools . . .) and parameterizations (Q10, plant respiration . . .). For example, the use of a single-pool SOC model (HadCM3LC) increases the sensitivity of SOC stocks to C input and temperature changes (Jones et al., 2005). In this study, no consensus emerges among the models for attributing SOC sensitivity to changes in NPP or to change in respiration .
Large uncertainties due to natural vegetation changes and anthropogenic CO 2 emissions were also evidenced by Sitch et al. (2008) using five Dynamic Global Vegetation models or DGVMs (Hyland, Triffid, LPJ, ORCHIDEE and Sheffield), four SRES scenarios of CO 2 emission (A1F1, A2, B1, B2) and a Fig. 6. Future evolution of global and regional soil organic carbon (SOC) stocks from 2000 to 2099 simulated by the 11 C4MIP model ). common climate model. DGVMs are models that simulate changes in vegetation and in its associated biogeochemical and hydrological cycles in response to changes in climatic conditions. As for the C4MIP study, the SOC change in Sitch et al. (2008) over the 21st century differs markedly between DGVMs. The range of model responses to climate change increases with intensive CO 2 emission scenarios. Finally, the global SOC stocks change was predicted to range between a gain of 169 PgC (Sheffield model with A1F1 scenario) and a loss of 58 PgC (LPJ model with A1F1 scenario). Again, major differences between models are localized in the tropics. For the most extreme emissions scenario (A1F1), three of the five DGVMs (Hyland, LPJ, TRIFFID) simulate strong forest dieback and soil carbon loss in both the Amazon and African tropical forests, whereas there is a small increase of SOC stocks for ORCHIDEE and Sheffield (Sitch et al., 2008). Interestingly, even the boreal region remains controversial. Contrarily to the other models, LPJ simulates large reduction in the boreal forest coverage and large high-latitudes losses in soil carbon.

Effect of future agricultural change
Using the OSCAR global carbon model driven by land cover change in the IPCC A2 scenario (Gitz and Ciais, 2004), we estimated potential SOC changes between 2000 and 2099 in response to a land-use change scenario. The results split into four main regions are given in Table 2 and compared to ORCHIDEE-LUC simulations for the 20th century. In all regions except Africa and Latin America, warming-induced increases in NPP dominates over warming-induced increases in decomposition, resulting in a net SOC gain. Thus, it can be seen that land cover change in this particular scenario [SRES A2 with land use from the IMAGE2.0 model (Alcamo, 1994)] is glob-ally reducing the SOC increase due to climate change. In most regions and at global scale, climate-induced NPP increase dominates SOC loss due to land-use change. However, deforestation and warming-increased soil respiration drives the SOC balance to be strongly negative in Africa and Latin American regions, strongly limiting the future C sequestration in the global soil.
Over Europe, more detailed simulations of SOC have been analysed. Smith et al. (2005a,b) using the Roth C model reported that projected climate change to 2080 might have limited impact on cropland NBP due to the balancing effects of increased losses due to faster decomposition and increased inputs due to choice of crops/harvestable fraction and improved technology. They attributed uncertainties in future projections mainly to differences in projected climate  by four GCMs forced by the four IPCC SRES emission scenarios, and differing assumptions about the implementation of technology. Smith et al. (2005a,b) also suggested that croplands could become a net carbon sink under improved technology, but the uncertainty (due to climate projections and interactions between decomposition, NPP and technology) was also large, with a 9% uncertainty due to forecast changes in technology alone.

Potential to increase SOC in agriculture
Agricultural lands occupy about 40-50% of the Earth's land surface. Agricultural practices can make a significant contribution at low cost to increasing soil carbon sinks Smith, 2008).
Soil C sequestration can be achieved by increasing the net flux of C from the atmosphere to the terrestrial biosphere by increasing global C inputs to the soil (via increasing NPP), by storing a larger proportion of the C from NPP in the longer-term C pools in the soil, or by reducing C losses from the soils by Table 2. Partitioning of global and regional SOC changes among land cover change and climate and CO 2 effects using ORCHIDEE-LUC model from 1901 to 2000 and OSCAR model from 2000 to 2099 ORCHIDEE-LUC (1901 OSCAR  Regions Total Climate + CO 2 Land cover Total Climate + CO 2 Land cover slowing decomposition. For soil C sinks, the best options are to increase C stocks in soils that have been depleted in C, that is agricultural soils and degraded soils, or to halt the loss of C from cultivated peatlands . In this section we focus on the potential to increase SOC in agriculture. About 90% of the total mitigation in agriculture arises from sink enhancement and about 10% arises from emission reduction . The most prominent mitigation options in agriculture are the restoration of cultivated organic soils, improved cropland management (including nutrient management, tillage/residue management and water management) and setaside/agro-forestry, improved grazing land management (including grazing intensity, increased productivity, nutrient management, fire management and species introduction, and restoration of degraded lands (using erosion control, organic amendments and nutrient amendments). The management practices needed to deliver this carbon sequestration in agriculture are described in detail in Smith et al. ( , 2008. Many reviews have been published recently discussing options available for soil C sequestration and mitigation potentials (see Smith, 2008 for a review). Most of the estimates for the sequestration potential of activities range from about 0.1 to 0.8 t C ha −1 yr −1 , but some estimates are outside this range (Smith, 2008). When considering soil C sequestration options, it is important also to consider other side effects, including the emission of other greenhouse gases. Smith et al. (2001) suggested that as much as one half of the climate mitigation potential of some C sequestration options could be lost when increased emissions of other greenhouse gases (nitrous oxide; N 2 O and methane; CH 4 ) were included, and Robertson et al. (2000) has shown that some practices that are beneficial for SOC sequestration, may not be beneficial when all greenhouse gases are considered.
The most recent estimate Smith, 2008) of the technical potential for SOC sequestration in agriculture globally is ∼1.3 PgC yr −1 , but this is very unlikely to be realized. Economic potentials for SOC sequestration estimated by  were 0.4, 0.6 and 0.7 PgC yr −1 at carbon prices of 0-20, 0-50 and 0-100 USD t CO 2 -equivalent −1 , respectively. At reasonable C prices, then, global agricultural soil C sequestration seems to be limited to around 0.4-0.7 PgC yr −1 . Even then, there are barriers (e.g. economic, institutional, educational, social) that impeded the economic potential to realize (Trines et al., 2006;Smith and Trines, 2007;Bellassen et al., 2010a). Because forest trees have a technical potential to sequester about 1-2 PgC y −1 (IPCC, 2000), the estimates for C sequestration potential in soils are of the same order. Similarly, the economic potential for C sequestration in forestry is similar to that for soil C sequestration in agriculture (IPCC, 2007).

Uncertainties and future directions
Many sources of uncertainties intervene in SOC change prediction in the context of atmospheric CO 2 increase and climate change, and land use change. These uncertainties include both future evolution of incomplete understanding of several biogeochemical processes (SOC decomposition and stabilization and fire disturbances) and ecosystem management (forest and agriculture). These uncertainties are only partly illustrated by the spread in future SOC stocks changes as predicted by the current process-based models (see Section 5 and Fig. 6). Combination of results by Sitch et al. (2008) and table 2 gives a maximum range of 368 PgC for SOC stocks changes during the 21st century.
First of all, because of incomplete understanding of some important processes (noted below) some potentially important biogeochemical processes are not incorporated in global models. Among the most important are (i) the 'priming effect'; (ii) temperature sensitivity of the different SOC pools and (iii) the role of the mineralogy in SOC stabilization. First, several studies have shown that fresh organic matter input to soil can stimulate SOC mineralization (e.g. Fontaine et al., 2004Fontaine et al., , 2007. This process called 'priming effect' might provide a negative feedback in the carbon cycle-climate system or even it might induce a positive one (Heimann and Reichstein, 2008), if its long-term impact is proven. Indeed, if C inputs stimulate SOC mineralization, enhanced inputs due to increased NPP will also increase SOC mineralization. In this case, the soil C sink forecast by global models would be smaller than currently estimated. Secondly, the vulnerability of the different SOC pools to warming is highly debated. Some authors reported that labile compounds are more sensitive to warming (Liski et al., 1999), that recalcitrant compounds are more sensitive (Fierer et al., 2005) or that both compounds are equally sensitive . At the moment, global C models assume that all pools will respond uniformly to warming. That assumption remains to be demonstrated. Third, several studies observed that most SOC is stabilized through interaction with mineral surfaces (Mikutta et al., 2006;Basile-Doelsch et al., 2009) and that mineral types have a major impact on SOC storage and turnover (e.g. Feller and Beare, 1997;Torn et al., 1997). SOC turnover usually depends on clay content in soil C models. The clay content dependence of the turnover integrates the fact that the higher the surface area the higher the amount of organic molecules that can be adsorbed and protected from microbial enzymes. However, not only the surface area but also the reactivity of this surface will determine the protection of organic molecules (Kleber et al., 2007). Consequently, it may be necessary to take soil mineralogy explicitly into account in global soil C models. These three examples are not exhaustive and illustrate that further breakthrough in soil biogeochemical research may lead to major changes in the architecture itself of soil C models applied at global scales.
Moreover, the effects of agricultural practices on SOC turnover and stability are still poorly understood. For example, even the increased C sequestration under less intensive tillage practices is currently questioned (Baker et al., 2007). Besides the increasing adoption of reduced tillage practices worldwide, recent studies tend to show that the effect of no till adoption on C storage has been overestimated, because too shallow soil depths have been taken into account in most studies. It is now argued that conservation tillage does not accrue SOC stock but only redistribute within the soil profile (Baker et al., 2007;Anger and Eriksen-Hamel, 2005). However, the tillage of previously untilled soils induces a large and rapid loss of soil C, which is poorly described by current models (Gottschalk et al., 2010) and the slower kinetics of C storage under adoption of untilled systems is even less taken into account. The controls exerted by soil structure on SOM mineralization and the feedbacks due to the role of SOM in stabilizing soil structure remain poorly conceptualized and described. Future evolution of cropland management is also a large source of uncertainties. Indeed, the changes in agricultural practices observed since the 1950s will continue and probably accelerate during the next century because of the development of new crop varieties (with improved productivity, smaller root/shoot ratios, etc.), environmental and agricultural policies, adaptation to climate change (irrigation), need for biofuel production, or other factors.
In forest, future management is also difficult to anticipate: on the one hand, more and more forest stands are either put into protection or certified under a sustainable forestry label, both of which limit the amount of wood that can be harvested. On the other hand, the rising demand for wood-energy and the increasing risk of catastrophic event (e.g. fire, pests, storms) with climate change will put an economic incentive to intensify management and shorten rotation length. Recent modelling studies are starting to assess the resulting outcome of these two trends. Pussinen et al. (2009), for example show that if felling is increased to meet the increase in wood demand, soil carbon in European forests will be reduced by about 10 tC ha −1 in 2100 compared to a 'current felling rate' scenario. In their model however, this decrease is offset by an increase in litterfall when the positive effect of climate change on tree growth is taken into account.
Finally, temporal changes in soil water deficit, windspeed, air temperature and humidity will modify the risk frequency and severity of forest fires or insect outbreaks and the consequent rapid loss of carbon from the biosphere (Kurz et al., 2008;Heimann and Reichstein, 2008). Using the Biome-BGC model, Bond-Lamberty et al. (2007) found that the C balance of the Canadian boreal forest was mainly driven by fire disturbances from 1948 and 2005. Thus, the increase in fire disturbances during this period resulted in carbon losses of 6.8 ± 1.0 gC m −2 yr −1 . But for soils, fire may be a global driver of long-term carbon storage because it converts a substantial proportion of plant biomass to black carbon, which remains in the soil for thousands of years (Wardle et al., 2008). Black carbon is therefore often proposed as an important long-term sink of soil carbon (Marris, 2006). Hence, the net effect of fire on soil C balance remains controversial. Less well documented are the impacts of insect outbreaks on the global carbon balance. Climate change will influence insect distribution and abundance (Kurz et al., 2008). For example, it is predicted that outbreaks of mountain pine beetle in Canada will increase in both scale and severity due to the expansion of its habitat and may have an impact on annual CO 2 emissions equivalent to forest fires through reduced NPP and increased heterotrophic respiration (Kurz et al., 2008). Insect outbreak may temporarily increase SOC stocks because of the large transfers of biomass to dead organic matter pools. But, as opposed to fire, no black carbon is produced.

Conclusions
SOC is the largest land biospheric carbon reservoir. Over the past century, climate and land cover change both have significantly changed the regional balance of SOC. Globally, according to simulations with the ORCHIDEE model covering the 20th Century, these two effects of opposite direction have nearly offset each other. In the simulation, tropical regions have gained carbon during the past century because 'CO 2 -and rainfall-induced' NPP increase dominated 'warming increased' decomposition and SOC loss due to conversion of natural ecosystems into agricultural lands. Oppositely, boreal and temperate regions have lost carbon due to clearing of forests into arable lands in the early 20th century. Although mechanization of agriculture in the 1950s has accelerated SOC losses in croplands, development of carbon-sequestering practices over the past decades (no-till, irrigation, reduced use of bare fallows, land conservation program, increased crop residue incorporation, etc.) may have limited SOC losses in the recent years. Nowadays, croplands do not probably represent a strong source of CO 2 to the atmosphere in Europe, China and USA. But regional, systematic inventories of the SOC balance of arable lands are basically lacking.
In the future, climate change may induce SOC release into the atmosphere if temperature and atmospheric CO 2 concentration do not stabilize. According to the results of the terrestrial biosphere model ensembles of Friedlingstein et al. (2006) and Sitch et al. (2008), increased litterfall might not compensate anymore warming-increased decomposition in the future. Because future natural and human-induced land-use changes are not likely to induce large SOC storage, soils are expected to become a net source of CO 2 to the atmosphere, accelerating the atmospheric CO 2 increase and the subsequent warming. The timing and magnitude of projected SOC changes are considerably uncertain. Although most carbon-cycle models predicts that the net SOC balance will still remain positive for the 21st century, several models simulates a fast SOC release during the last decades leading in some case to a negative SOC balance. The most vulnerable regions are consistently the tropical forests because of drought-induced forest dieback. But large high-latitudes C losses could also be expected because of permafrost thawing (not included in the model simulations presented in this review). A better incorporation into models of biogeochemical processes affecting SOC balance together with the identification of the most relevant future anthropogenic CO 2 emissions and land-use scenarios are needed to reduce the uncertainties associated with the predictions of future SOC changes. Besides, the positive influence of some carbon-sequestering agricultural practices on SOC balance may encourage their use for climate change mitigation. More recent estimate of the technical potential for SOC sequestration in agriculture gives a potential global annual sink of ∼1.3 PgC yr −1 , or a total SOC accumulation of ∼27 PgC by 2030. However, this restoration of SOC is unlikely to be realized because of economic, institutional, educational and social barriers.

ORCHIDEE-LUC
The ORCHIDEE model (Krinner et al., 2005) is a dynamic global vegetation model (DGVM) representing key vegetation processes governing terrestrial biogeochemistry and biogeography. In this study, however, we switched off vegetation dynamic simulation module because there are large uncertainties in DGVM modelled vegetation distribution (Krinner et al., 2005;Sitch et al., 2008). ORCHIDEE distinguishes 12 plant functional types with different photosynthetic, phenological and morphological characteristics. Heterotrophic respiration parameterization is taken from CENTURY (Parton et al., 1988). The current version of ORCHIDEE takes into account the effects of land use change on terrestrial C cycle as described in Piao et al. (2009a,b).

ORCHIDEE-STICS
ORCHIDEE-STICS is an evolution of the ORCHIDEE ecosystem model that is coupled to the STICS agronomy model describing crop phenology. More details of the model and applications can be found in de Noblet-Ducoudre et al. (2004), Gervois et al. (2004) and Gervois et al. (2008). ORCHIDEE-STICS was used to simulate changes in the European croplands C balance between 1901 and 2000. This simulation included winter wheat and maize, and ignored summer C3 crops (Gervois et al., 2008). The soil carbon decomposition module of ORCHIDEE-STICS is similar to the CENTURY model equations (Parton et al., 1988). Starting from ancestral farming practice and crop varieties in 1901, the model was driven over Western Europe by rising CO 2 and reconstructed climate at a resolution of 10 km, and by evolving agricultural technology. The technological evolution is reconstructed by Gervois et al. (2008), and includes the use of shorter growth cycle varieties, increased harvest index, increased N-fertilizers and decreased manure applications, maize irrigation, and increased tillage which accelerates SOC decomposition.

ORCHIDEE-FM
ORCHIDEE-FM is an evolution of the ORCHIDEE ecosystem model. It includes a Forest Management Module that explicitly simulates forest stand growth and management (Bellassen et al., 2010b). This is made possible by the disaggregation of forest representation from an average tree to an average stand defined by a population of diameter classes. The model is able to calculate the critical stand self-thinning density above which natural mortality occurs due to competition. Based on the standard forestry principle of avoiding natural mortality, a set of rules are defined to calculate the recurrent intensity and frequency of thinning and forestry operations during the stand lifetime.

OSCAR
OSCAR is a global carbon cycle model including a reduced-form ocean model to quantify the ocean-atmosphere CO2 exchange and a terrestrial carbon cycle model to account for the fluxes between land and atmosphere (Gitz and Ciais, 2004). The parameterization of the biosphere was derived from the spatiallyexplicit CASA-SLAVE model . The terrestrial cycle integrates a detailed land-use module that allow for conversions of biomes, and calculates both the land-use related net CO2 emissions following anthropogenic disturbances as well as the terrestrial uptakes over the remaining undisturbed ecosystems at each time step. The global land cover is based on a simplified vegetation map, which is regionalized into four world regions as defined by the Intergovernmental Panel of Climate Change Third Assessment Report (IPCC, 2000) and shown in table 2: OECD-1990 (North America, Europe, Japan and Australia), Former Soviet Union, ASIA, ALM (Africa, Latin America and Middle East). In each region, six natural biomes are defined, plus three crop types (boreal, temperate, tropical) and cohorts of lands in transition between two biomes. Icy and hot desert were excluded. For each region, a biome has separate biophysical characteristics and it is assigned one surface area in the model.

RothC
RothC is a generic soil carbon model in non-waterlogged soils that allows for the effects of soil texture, temperature, moisture content and plant cover on the SOC turnover process (Jenkinson et al., 1987(Jenkinson et al., , 1991Jenkinson, 1990). Soil carbon is split into four active compartments and a small amount of inert organic matter (IOM). The four active compartments are Decomposable Plant Material (DPM), Resistant Plant Material (RPM), Microbial Biomass (BIO) and Humified Organic Matter. Each compartment decomposes by a first-order process with its own characteristic rate. The IOM compartment is resistant to decomposition. Smith et al. (2005a,b) applied this model over Europe prescribed with input from changing NPP-discounted for harvest-from the LPJ vegetation model (Sitch et al., 2003). The RothC model was initialized and run at a resolution of 10 km from 1900 to 2100 and for four different IPCC future scenarios.

LPJmL
LPJmL where mL stands for managed lands is an evolution of LPJ (Sitch et al., 2003) that has been extended to represent the carbon and water cycles of managed land (Bondeau et al., 2007). The version cited in this study considers 11 crops functional types, eight of them being cultivated in Europe: temperate cereals (wheat like), maize, soybean, rice, temperate roots (sugarbeet like), sunflowers, rapeseed, pulses. The Leff et al. (2004) crop distribution is used to determine the grid cell fraction cover of each of these types within the cropland cover database provided for CarboEurope (Vetter et al., 2008) for the window (15 • W-60 • E; 30 • N-75 • N). The remaining crops are put within the temperate cereals type. For each type and each grid cell, the most appropriate crop variety is determined according to the climatic conditions. The sowing date, the heat unit requirements, and the base temperature, differ for each crop over the European window. For example, spring wheat is grown in higher latitudes, whereas winter wheat is grown over central and southern Europe.

C4MIP
In the context of the Coupled Climate-Carbon Cycle Model Intercomparison Project (C4MIP), seven coupled Ocean-Atmosphere General Circulation Models (OAGCMs) and four models of intermediate complexity performed coupled climate-carbon cycle simulations over the historical period and the 21st century . All models used observed anthropogenic fossil fuel emissions for the historical period (Marland et al., 2005) and the Intergovernmental Panel on Climate Change (IPCC) Special Report on Emissions Scenarios (SRES) A2 emission scenario for the 2000-2100 period. Most models included land-use-associated CO2 emissions provided by Houghton and Hackler (2002) for the historical and by the Integrated Model to Assess the Global Environment (IMAGE)-integrated model for the 21st century (Leemans et al., 1998). Land-use-associated emissions are seen here as an external forcing. More details on the models used in this study are presented in Friedlingstein et al. (2006).