The greenhouse gas balance of European grasslands

The greenhouse gas (GHG) balance of European grasslands (EU‐28 plus Norway and Switzerland), including CO₂, CH₄ and N₂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 N₂O and CH₄ 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 CO₂ 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 CO₂ fluxes and CO₂ 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 CH₄ and N₂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‐CO₂ equiv. m^?2 year^?1, because the CO₂ sink offsets N₂O and grazing animal CH₄ emissions. However, when considering the farm scale, the GHG balance (NGB) becomes a net GHG source of ?50 g C‐CO₂ 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 CO₂, climate and nitrogen fertilization and the diminishing requirement for grass forage due to the Europe‐wide reduction in livestock numbers.     

Projected changes in mineral soil carbon of European croplands and grasslands, 1990–2080

We present the most comprehensive pan‐European assessment of future changes in cropland and grassland soil organic carbon (SOC) stocks to date, using a dedicated process‐based SOC model and state‐of‐the‐art databases of soil, climate change, land‐use change and technology change. Soil carbon change was calculated using the Rothamsted carbon model on a European 10 × 10′ grid using climate data from four global climate models implementing four Intergovernmental Panel on Climate Change (IPCC) emissions scenarios (SRES). Changes in net primary production (NPP) were calculated by the Lund–Potsdam–Jena model. Land‐use change scenarios, interpreted from the narratives of the IPCC SRES story lines, were used to project changes in cropland and grassland areas. Projections for 1990–2080 are presented for mineral soil only. Climate effects (soil temperature and moisture) will tend to speed decomposition and cause soil carbon stocks to decrease, whereas increases in carbon input because of increasing NPP will slow the loss. Technological improvement may further increase carbon inputs to the soil. Changes in cropland and grassland areas will further affect the total soil carbon stock of European croplands and grasslands. While climate change will be a key driver of change in soil carbon over the 21st Century, changes in technology and land‐use change are estimated to have very significant effects. When incorporating all factors, cropland and grassland soils show a small increase in soil carbon on a per area basis under future climate (1–7 t C ha^?1 for cropland and 3–6 t C ha^?1 for grassland), but when the greatly decreasing area of cropland and grassland are accounted for, total European cropland stocks decline in all scenarios, and grassland stocks decline in all but one scenario. Different trends are seen in different regions. For Europe (the EU25 plus Norway and Switzerland), the cropland SOC stock decreases from 11 Pg in 1990 by 4–6 Pg (39–54%) by 2080, and the grassland SOC stock increases from 6 Pg in 1990 to 1.5 Pg (25%) under the B1 scenario, but decreases to 1–3 Pg (20–44%) under the other scenarios. Uncertainty associated with the land‐use and technology scenarios remains unquantified, but worst‐case quantified uncertainties are 22.5% for croplands and 16% for grasslands, equivalent to potential errors of 2.5 and 1 Pg SOC, respectively. This is equivalent to 42–63% of the predicted SOC stock change for croplands and 33–100% of the predicted SOC stock change for grasslands. Implications for accounting for SOC changes under the Kyoto Protocol are discussed.     

Change in terrestrial ecosystem water‐use efficiency over the last three decades

Defined as the ratio between gross primary productivity (GPP) and evapotranspiration (ET), ecosystem‐scale water‐use efficiency (EWUE) is an indicator of the adjustment of vegetation photosynthesis to water loss. The processes controlling EWUE are complex and reflect both a slow evolution of plants and plant communities as well as fast adjustments of ecosystem functioning to changes of limiting resources. In this study, we investigated EWUE trends from 1982 to 2008 using data‐driven models derived from satellite observations and process‐oriented carbon cycle models. Our findings suggest positive EWUE trends of 0.0056, 0.0007 and 0.0001 g C m^?2 mm^?1 yr^?1 under the single effect of rising CO₂ (‘CO₂’), climate change (‘CLIM’) and nitrogen deposition (‘NDEP’), respectively. Global patterns of EWUE trends under different scenarios suggest that (i) EWUE‐CO₂ shows global increases, (ii) EWUE‐CLIM increases in mainly high latitudes and decreases at middle and low latitudes, (iii) EWUE‐NDEP displays slight increasing trends except in west Siberia, eastern Europe, parts of North America and central Amazonia. The data‐driven MTE model, however, shows a slight decline of EWUE during the same period (?0.0005 g C m^?2 mm^?1 yr^?1), which differs from process‐model (0.0064 g C m^?2 mm^?1 yr^?1) simulations with all drivers taken into account. We attribute this discrepancy to the fact that the nonmodeled physiological effects of elevated CO₂ reducing stomatal conductance and transpiration (TR) in the MTE model. Partial correlation analysis between EWUE and climate drivers shows similar responses to climatic variables with the data‐driven model and the process‐oriented models across different ecosystems. Change in water‐use efficiency defined from transpiration‐based WUE_t (GPP/TR) and inherent water‐use efficiency (IWUE_t, GPP×VPD/TR) in response to rising CO₂, climate change, and nitrogen deposition are also discussed. Our analyses will facilitate mechanistic understanding of the carbon–water interactions over terrestrial ecosystems under global change.     

Effects of land use change and management on the European cropland carbon balance

We model the carbon balance of European croplands between 1901 and 2000 in response to land use and management changes. The process‐based ORCHIDEE‐STICS model is applied here in a spatially explicit framework. We reconstructed land cover changes, together with an idealized history of agro‐technology. These management parameters include the treatment of straw and stubble residues, application of mineral fertilizers, improvement of cultivar species and tillage. The model is integrated for wheat and maize during the period 1901–2000 forced by climate each 1/2‐hour, and by atmospheric CO₂, land cover change and agro‐technology each year. Several tests are performed to identify the most sensitive agro‐technological parameters that control the net biome productivity (NBP) in the 1990s, with NBP equaling for croplands the soil C balance. The current NBP is a small sink of 0.16 t C ha^?1 yr^?1. The value of NBP per unit area reflects past and current management, and to a minor extent the shrinking areas of arable land consecutive to abandonment during the 20th Century. The uncertainty associated with NBP is large, with a 1‐sigma error of 0.18 t C ha^?1 yr^?1 obtained from a qualitative, but comprehensive budget of various error terms. The NBP uncertainty is dominated by unknown historical agro‐technology changes (47%) and model structure (27%), with error in climate forcing playing a minor role. A major improvement to the framework would consist in using a larger number of representative crops. The uncertainty of historical land‐use change derived from three different reconstructions, has a surprisingly small effect on NBP (0.01 t C ha^?1 yr^?1) because cropland area remained stable during the past 20 years in all the tested land use forcing datasets. Regional cross‐validation of modeled NBP against soil C inventory measurements shows that our results are consistent with observations, within the uncertainties of both inventories and model. Our estimation of cropland NBP is however likely to be biased towards a sink, given that inventory data from different regions consistently indicate a small source whereas we model a small sink.     

Reconstruction and attribution of the carbon sink of European forests between 1950 and 2000

European forests are an important carbon sink; however, the relative contributions to this sink of climate, atmospheric CO₂ concentration ([CO₂]), nitrogen deposition and forest management are under debate. We attributed the European carbon sink in forests using ORCHIDEE‐FM, a process‐based vegetation model that differs from earlier versions of ORCHIDEE by its explicit representation of stand growth and idealized forest management. The model was applied on a grid across Europe to simulate changes in the net ecosystem productivity (NEP) of forests with and without changes in climate, [CO₂] and age structure, the three drivers represented in ORCHIDEE‐FM. The model simulates carbon stocks and volume increment that are comparable – root mean square error of 2 m³ ha^?1 yr^?1 and 1.7 kg C m^?2 respectively – with inventory‐derived estimates at country level for 20 European countries. Our simulations estimate a mean European forest NEP of 175 ± 52 g C m^?2 yr^?1 in the 1990s. The model simulation that is most consistent with inventory records provides an upwards trend of forest NEP of 1 ± 0.5 g C m^?2 yr^?2 between 1950 and 2000 across the EU 25. Furthermore, the method used for reconstructing past age structure was found to dominate its contribution to temporal trends in NEP. The potentially large fertilizing effect of nitrogen deposition cannot be told apart, as the model does not explicitly simulate the nitrogen cycle. Among the three drivers that were considered in this study, the fertilizing effect of increasing [CO₂] explains about 61% of the simulated trend, against 26% to changes in climate and 13% only to changes in forest age structure. The major role of [CO₂] at the continental scale is due to its homogeneous impact on net primary productivity (NPP). At the local scale, however, changes in climate and forest age structure often dominate trends in NEP by affecting NPP and heterotrophic respiration.     

A comparison of repeated soil inventory and carbon flux budget to detect soil carbon stock changes after conversion from cropland to grasslands

Assessment of soil carbon (C) stock changes over time is typically based on the application of two methods, namely (i) repeated soil inventory and (ii) determination of the ecosystem C budget or net biome productivity (NBP) by continuous measurement of CO₂ exchange in combination with quantification of other C imports and exports. Here, we applied both methods in parallel to determine C stock changes of two temperate grassland fields previously converted from long‐term cropland. The grasslands differed in management intensity with either intensive management (high fertilization, frequent cutting) or extensive management (no fertilization, less frequent cutting). Soil organic C stocks (0–45 cm depth) were quantified at the beginning (2001) and the end (2006) of a 5 year observational period using the equivalent soil mass approach. For the same period and in both fields, NBP was quantified from net CO₂ fluxes monitored using eddy covariance systems, and measured C import by organic fertilizer and C export by harvest. Both NBP and repeated soil inventories revealed a consistent and significant difference between management systems of 170 ± 48 and 253 ± 182  g C m^?2 a^?1, respectively. For both fields, the inventory method showed a tendency towards higher C loss/smaller C gain than NBP. In the extensive field, a significant C loss was observed by the inventory but not by the NBP approach. Thus neither flux measurements nor repeated soil sampling may be suitable for tracking absolute changes in SOC, but both give similar answers with respect to relative changes.     

Evaluation of terrestrial carbon cycle models for their response to climate variability and to CO2 trends

The purpose of this study was to evaluate 10 process‐based terrestrial biosphere models that were used for the IPCC fifth Assessment Report. The simulated gross primary productivity (GPP) is compared with flux‐tower‐based estimates by Jung et al. [Journal of Geophysical Research 116 (2011) G00J07] (JU11). The net primary productivity (NPP) apparent sensitivity to climate variability and atmospheric CO₂ trends is diagnosed from each model output, using statistical functions. The temperature sensitivity is compared against ecosystem field warming experiments results. The CO₂ sensitivity of NPP is compared to the results from four Free‐Air CO₂ Enrichment (FACE) experiments. The simulated global net biome productivity (NBP) is compared with the residual land sink (RLS) of the global carbon budget from Friedlingstein et al. [Nature Geoscience 3 (2010) 811] (FR10). We found that models produce a higher GPP (133 ± 15 Pg C yr^?1) than JU11 (118 ± 6 Pg C yr^?1). In response to rising atmospheric CO₂ concentration, modeled NPP increases on average by 16% (5–20%) per 100 ppm, a slightly larger apparent sensitivity of NPP to CO₂ than that measured at the FACE experiment locations (13% per 100 ppm). Global NBP differs markedly among individual models, although the mean value of 2.0 ± 0.8 Pg C yr^?1 is remarkably close to the mean value of RLS (2.1 ± 1.2 Pg C yr^?1). The interannual variability in modeled NBP is significantly correlated with that of RLS for the period 1980–2009. Both model‐to‐model and interannual variation in model GPP is larger than that in model NBP due to the strong coupling causing a positive correlation between ecosystem respiration and GPP in the model. The average linear regression slope of global NBP vs. temperature across the 10 models is ?3.0 ± 1.5 Pg C yr^?1 °C^?1, within the uncertainty of what derived from RLS (?3.9 ± 1.1 Pg C yr^?1 °C^?1). However, 9 of 10 models overestimate the regression slope of NBP vs. precipitation, compared with the slope of the observed RLS vs. precipitation. With most models lacking processes that control GPP and NBP in addition to CO₂ and climate, the agreement between modeled and observation‐based GPP and NBP can be fortuitous. Carbon–nitrogen interactions (only separable in one model) significantly influence the simulated response of carbon cycle to temperature and atmospheric CO₂ concentration, suggesting that nutrients limitations should be included in the next generation of terrestrial biosphere models.     

Potential carbon sequestration of European arable soils estimated by modelling a comprehensive set of management practices

Bottom–up estimates from long‐term field experiments and modelling are the most commonly used approaches to estimate the carbon (C) sequestration potential of the agricultural sector. However, when data are required at European level, important margins of uncertainty still exist due to the representativeness of local data at large scale or different assumptions and information utilized for running models. In this context, a pan‐European (EU + Serbia, Bosnia and Herzegovina, Montenegro, Albania, Former Yugoslav Republic of Macedonia and Norway) simulation platform with high spatial resolution and harmonized data sets was developed to provide consistent scenarios in support of possible carbon sequestration policies. Using the CENTURY agroecosystem model, six alternative management practices (AMP) scenarios were assessed as alternatives to the business as usual situation (BAU). These consisted of the conversion of arable land to grassland (and vice versa), straw incorporation, reduced tillage, straw incorporation combined with reduced tillage, ley cropping system and cover crops. The conversion into grassland showed the highest soil organic carbon (SOC) sequestration rates, ranging between 0.4 and 0.8 t C ha^?1 yr^?1, while the opposite extreme scenario (100% of grassland conversion into arable) gave cumulated losses of up to 2 Gt of C by 2100. Among the other practices, ley cropping systems and cover crops gave better performances than straw incorporation and reduced tillage. The allocation of 12 to 28% of the European arable land to different AMP combinations resulted in a potential SOC sequestration of 101–336 Mt CO₂ eq. by 2020 and 549‐2141 Mt CO₂ eq. by 2100. Modelled carbon sequestration rates compared with values from an ad hoc meta‐analysis confirmed the robustness of these estimates.     

Dynamics of soil carbon following destruction of tropical rainforest and the subsequent establishment of Imperata grassland in Indonesian Borneo using stable carbon isotopes

Southeast Asia has the highest rate of tropical rainforest deforestation worldwide, and large deforested areas have been replaced ultimately by the highly invasive grass Imperata cylindrica. However, information on the carbon (C) budget with such land transition is very scarce. This study presents the dynamics of soil C following rainforest destruction and the subsequent establishment of Imperata grassland in the lowland humid tropics of Indonesian Borneo using stable C isotopes. To evaluate the relative contribution of organic matter originating from primary forest (C₃) and grasslands (C₄), we compared soil C stock and natural ¹³C abundance from six sites to a depth of 100 cm using samples with a wide range of soil textures. Twelve years after the first soil sampling in the grasslands, we re‐sampled to examine temporal changes in soil organic matter. The grassland topsoil (0–5 cm) is an active layer with rapid decomposition and incorporation of fresh C (mean residence time: 7.5 year) and a substantial proportion of the stable C pool (37%). The decline in forest‐derived C was slight, even at 5–10 cm depths, and subsoil (20–100 cm depth) forest‐derived C did not change along the forest‐to‐grassland chronosequence. Grassland‐derived C stock increased significantly in the subsurface and subsoils (5–100 cm). Simulation indicated that total soil C stock (0–100 cm) increased by 18.6 Mg ha^?1 from initial primary forest (58.0 Mg ha^?1) to a new equilibrium state of the grassland (76.6 Mg ha^?1) after 30–50 years of grassland establishment. This research indicates that the soil did not function as a CO₂ source when the deforested area was replaced by Imperata grassland on the Ultisols of the Asian humid tropics. Instead, increased soil C stocks offset CO₂ emissions, with the C offset accounting for 6.6–7.4% of the loss of biomass C stock.     

Partitioning of ecosystem respiration of CO2 released during land‐use transition from temperate agricultural grassland to Miscanthus × giganteus

Conversion of large areas of agricultural grassland is inevitable if European and UK domestic production of biomass is to play a significant role in meeting demand. Understanding the impact of these land‐use changes on soil carbon cycling and stocks depends on accurate predictions from well‐parameterized models. Key considerations are cultivation disturbance and the effect of autotrophic root input stimulation on soil carbon decomposition under novel biomass crops. This study presents partitioned parameters from the conversion of semi‐improved grassland to Miscanthus bioenergy production and compares the contribution of autotrophic and heterotrophic respiration to overall ecosystem respiration of CO₂ in the first and second years of establishment. Repeated measures of respiration from within and without root exclusion collars were used to produce time‐series model integrations separating live root inputs from decomposition of grass residues ploughed in with cultivation of the new crop. These parameters were then compared to total ecosystem respiration derived from eddy covariance sensors. Average soil surface respiration was 13.4% higher in the second growing season, increasing from 2.9 to 3.29 g CO₂‐C m^?2 day^?1. Total ecosystem respiration followed a similar trend, increasing from 4.07 to 5.4 g CO₂‐C m^?2 day^?1. Heterotrophic respiration from the root exclusion collars was 32.2% lower in the second growing season at 1.20 g CO₂‐C m^?2 day^?1 compared to the previous year at 1.77 g CO₂‐C m^?2 day^?1. Of the total respiration flux over the two‐year time period, aboveground autotrophic respiration plus litter decomposition contributed 38.46% to total ecosystem respiration while belowground autotrophic respiration and stimulation by live root inputs contributed 46.44% to soil surface respiration. This figure is notably higher than mean figures for nonforest soils derived from the literature and demonstrates the importance of crop‐specific parameterization of respiration models.     

Seasonal responses of terrestrial ecosystem water‐use efficiency to climate change

Ecosystem water‐use efficiency (EWUE) is an indicator of carbon–water interactions and is defined as the ratio of carbon assimilation (GPP) to evapotranspiration (ET). Previous research suggests an increasing long‐term trend in annual EWUE over many regions and is largely attributed to the physiological effects of rising CO₂. The seasonal trends in EWUE, however, have not yet been analyzed. In this study, we investigate seasonal EWUE trends and responses to various drivers during 1982–2008. The seasonal cycle for two variants of EWUE, water‐use efficiency (WUE, GPP/ET), and transpiration‐based WUE (WUE_t, the ratio of GPP and transpiration), is analyzed from 0.5° gridded fields from four process‐based models and satellite‐based products, as well as a network of 63 local flux tower observations. WUE derived from flux tower observations shows moderate seasonal variation for most latitude bands, which is in agreement with satellite‐based products. In contrast, the seasonal EWUE trends are not well captured by the same satellite‐based products. Trend analysis, based on process‐model factorial simulations separating effects of climate, CO₂, and nitrogen deposition (NDEP), further suggests that the seasonal EWUE trends are mainly associated with seasonal trends of climate, whereas CO₂ and NDEP do not show obvious seasonal difference in EWUE trends. About 66% grid cells show positive annual WUE trends, mainly over mid‐ and high northern latitudes. In these regions, spring climate change has amplified the effect of CO₂ in increasing WUE by more than 0.005 gC m⁻² mm⁻¹ yr⁻¹ for 41% pixels. Multiple regression analysis further shows that the increase in springtime WUE in the northern hemisphere is the result of GPP increasing faster than ET because of the higher temperature sensitivity of GPP relative to ET. The partitioning of annual EWUE to seasonal components provides new insight into the relative sensitivities of GPP and ET to climate, CO_2, and NDEP.     

The European carbon balance. Part 3: forests

We present a new synthesis, based on a suite of complementary approaches, of the primary production and carbon sink in forests of the 25 member states of the European Union (EU‐25) during 1990–2005. Upscaled terrestrial observations and model‐based approaches agree within 25% on the mean net primary production (NPP) of forests, i.e. 520±75 g C m^?2 yr^?1 over a forest area of 1.32 × 10⁶ km² to 1.55 × 10⁶ km² (EU‐25). New estimates of the mean long‐term carbon forest sink (net biome production, NBP) of EU‐25 forests amounts 75±20 g C m^?2 yr^?1. The ratio of NBP to NPP is 0.15±0.05. Estimates of the fate of the carbon inputs via NPP in wood harvests, forest fires, losses to lakes and rivers and heterotrophic respiration remain uncertain, which explains the considerable uncertainty of NBP. Inventory‐based assessments and assumptions suggest that 29±15% of the NBP (i.e., 22 g C m^?2 yr^?1) is sequestered in the forest soil, but large uncertainty remains concerning the drivers and future of the soil organic carbon. The remaining 71±15% of the NBP (i.e., 53 g C m^?2 yr^?1) is realized as woody biomass increments. In the EU‐25, the relatively large forest NBP is thought to be the result of a sustained difference between NPP, which increased during the past decades, and carbon losses primarily by harvest and heterotrophic respiration, which increased less over the same period.     

Causes of slowing‐down seasonal CO2 amplitude at Mauna Loa

Changing amplitude of the seasonal cycle of atmospheric CO₂ (SCA) in the northern hemisphere is an emerging carbon cycle property. Mauna Loa (MLO) station (20°N, 156°W), which has the longest continuous northern hemisphere CO₂ record, shows an increasing SCA before the 1980s (p < .01), followed by no significant change thereafter. We analyzed the potential driving factors of SCA slowing‐down, with an ensemble of dynamic global vegetation models (DGVMs) coupled with an atmospheric transport model. We found that slowing‐down of SCA at MLO is primarily explained by response of net biome productivity (NBP) to climate change, and by changes in atmospheric circulations. Through NBP, climate change increases SCA at MLO before the 1980s and decreases it afterwards. The effect of climate change on the slowing‐down of SCA at MLO is mainly exerted by intensified drought stress acting to offset the acceleration driven by CO₂ fertilization. This challenges the view that CO₂ fertilization is the dominant cause of emergent SCA trends at northern sites south of 40°N. The contribution of agricultural intensification on the deceleration of SCA at MLO was elusive according to land–atmosphere CO₂ flux estimated by DGVMs and atmospheric inversions. Our results also show the necessity to adequately account for changing circulation patterns in understanding carbon cycle dynamics observed from atmospheric observations and in using these observations to benchmark DGVMs.     

The European carbon balance. Part 2: croplands

We estimated the long‐term carbon balance [net biome production (NBP)] of European (EU‐25) croplands and its component fluxes, over the last two decades. Net primary production (NPP) estimates, from different data sources ranged between 490 and 846 gC m^?2 yr^?1, and mostly reflect uncertainties in allocation, and in cropland area when using yield statistics. Inventories of soil C change over arable lands may be the most reliable source of information on NBP, but inventories lack full and harmonized coverage of EU‐25. From a compilation of inventories we infer a mean loss of soil C amounting to 17 g m^?2 yr^?1. In addition, three process‐based models, driven by historical climate and evolving agricultural technology, estimate a small sink of 15 g C m^?2 yr^?1 or a small source of 7.6 g C m^?2 yr^?1. Neither the soil C inventory data, nor the process model results support the previous European‐scale NBP estimate by Janssens and colleagues of a large soil C loss of 90 ± 50 gC m^?2 yr^?1. Discrepancy between measured and modeled NBP is caused by erosion which is not inventoried, and the burning of harvest residues which is not modeled. When correcting the inventory NBP for the erosion flux, and the modeled NBP for agricultural fire losses, the discrepancy is reduced, and cropland NBP ranges between ?8.3 ± 13 and ?13 ± 33 g C m^?2 yr^?1 from the mean of the models and inventories, respectively. The mean nitrous oxide (N₂O) flux estimates ranges between 32 and 37 g C Eq m^?2 yr^?1, which nearly doubles the CO₂ losses. European croplands act as small CH₄ sink of 3.3 g C Eq m^?2 yr^?1. Considering ecosystem CO₂, N₂O and CH₄ fluxes provides for the net greenhouse gas balance a net source of 42–47 g C Eq m^?2 yr^?1. Intensifying agriculture in Eastern Europe to the same level Western Europe amounts is expected to result in a near doubling of the N₂O emissions in Eastern Europe. N₂O emissions will then become the main source of concern for the impact of European agriculture on climate.     

Accelerating regrowth of temperate‐maritime forests due to environmental change

To understand how environmental changes have influenced forest productivity, stemwood biomass (B) dynamics were analyzed at 1267 permanent inventory plots, covering a combined 209 ha area of unmanaged temperate‐maritime forest in southwest British Columbia, Canada. Net stemwood production (ΔB) was derived from periodic remeasurements of B collected over a 40‐year measurement period (1959–1998) in stands ranging from 20 to 150 years old. Comparison between the integrated age response of net stemwood production, ΔB(A), and the age response of stemwood biomass, B(A), suggested a 58 ± 11% increase in ΔB between the first 40 years of the chronosequence period (1859–1898) and the measurement period. To estimate extrinsic forcing on ΔB, several different candidate models were developed to remove variation explained by intrinsic factors. All models exhibited temporal bias, with positive trends in (observed minus predicted) residual ΔB ranging between of 0.40 and 0.64% yr^?1. Applying the same methods to stemwood growth (G) indicated residual increases ranging from 0.43 and 0.67% yr^?1. Higher trend estimates corresponded with models that included site index (SI) as a predictor, which may reflect exaggeration of the age‐decline in SI tables. Choosing a model that excluded SI, suggested that ΔB increased by 0.40 ± 0.18% yr^?1, while G increased by 0.43 ± 0.12% yr^?1 over the measurement period. Residual G was significantly correlated with atmospheric carbon dioxide (CO₂), temperature (T), and climate moisture index (CMI). However, models driven with climate and CO₂, alone, could not simultaneously explain long‐term and measurement‐period trends without additional representation of indirect effects, perhaps reflecting compound interest on direct physiological responses to environmental change. Evidence of accelerating forest regrowth highlights the value of permanent inventories to detect and understand systematic changes in forest productivity caused by environmental change.     

The European carbon balance. Part 4: integration of carbon and other trace‐gas fluxes

Overviewing the European carbon (C), greenhouse gas (GHG), and non‐GHG fluxes, gross primary productivity (GPP) is about 9.3 Pg yr^?1, and fossil fuel imports are 1.6 Pg yr^?1. GPP is about 1.25% of solar radiation, containing about 360 × 10¹⁸ J energy – five times the energy content of annual fossil fuel use. Net primary production (NPP) is 50%, terrestrial net biome productivity, NBP, 3%, and the net GHG balance, NGB, 0.3% of GPP. Human harvest uses 20% of NPP or 10% of GPP, or alternatively 1‰ of solar radiation after accounting for the inherent cost of agriculture and forestry, for production of pesticides and fertilizer, the return of organic fertilizer, and for the C equivalent cost of GHG emissions. C equivalents are defined on a global warming potential with a 100‐year time horizon. The equivalent of about 2.4% of the mineral fertilizer input is emitted as N₂O. Agricultural emissions to the atmosphere are about 40% of total methane, 60% of total NO‐N, 70% of total N₂O‐N, and 95% of total NH₃‐N emissions of Europe. European soils are a net C sink (114 Tg yr^?1), but considering the emissions of GHGs, soils are a source of about 26 Tg CO₂ C‐equivalent yr^?1. Forest, grassland and sediment C sinks are offset by GHG emissions from croplands, peatlands and inland waters. Non‐GHGs (NH₃, NOx) interact significantly with the GHG and the C cycle through ammonium nitrate aerosols and dry deposition. Wet deposition of nitrogen (N) supports about 50% of forest timber growth. Land use change is regionally important. The absolute flux values total about 50 Tg C yr^?1. Nevertheless, for the European trace‐gas balance, land‐use intensity is more important than land‐use change. This study shows that emissions of GHGs and non‐GHGs significantly distort the C cycle and eliminate apparent C sinks.     

Net biome production of the Amazon Basin in the 21st century

Global change includes multiple stressors to natural ecosystems ranging from direct climate and land‐use impacts to indirect degradation processes resulting from fire. Humid tropical forests are vulnerable to projected climate change and possible synergistic interactions with deforestation and fire, which may initiate a positive feedback to rising atmospheric CO₂. Here, we present results from a multifactorial impact analysis that combined an ensemble of climate change models with feedbacks from deforestation and accidental fires to quantify changes in Amazon Basin carbon cycling. Using the LPJmL Dynamic Global Vegetation Model, we modelled spatio‐temporal changes in net biome production (NBP); the difference between carbon fluxes from fire, deforestation, soil respiration and net primary production. By 2050, deforestation and fire (with no CO₂ increase or climate change) resulted in carbon losses of 7.4–20.3 Pg C with the range of uncertainty depending on socio‐economic storyline. During the same time period, interactions between climate and land use either compensated for carbon losses due to wetter climate and CO₂ fertilization or exacerbated carbon losses from drought‐induced forest mortality (?20.1 to +4.3 Pg C). By the end of the 21st century, depending on climate projection and the rate of deforestation (including its interaction with fire), carbon stocks either increased (+12.6 Pg C) or decreased (?40.6 Pg C). The synergistic effect of deforestation and fire with climate change contributed up to 26–36 Pg C of the overall decrease in carbon stocks. Agreement between climate projections (n=9), not accounting for deforestation and fire, in 2050 and 2098 was relatively low for the directional change in basin‐wide NBP (19–37%) and aboveground live biomass (13–24%). The largest uncertainty resulted from climate projections, followed by implementation of ecosystem dynamics and deforestation. Our analysis partitions the drivers of tropical ecosystem change and is relevant for guiding mitigation and adaptation policy related to global change.     

CO2 emissions from land‐use change affected more by nitrogen cycle,than by the choice of land‐cover data

The high uncertainty in land‐based CO₂ fluxes estimates is thought to be mainly due to uncertainty in not only quantifying historical changes among forests, croplands, and grassland, but also due to different processes included in calculation methods. Inclusion of a nitrogen (N) cycle in models is fairly recent and strongly affects carbon (C) fluxes. In this study, for the first time, we use a model with C and N dynamics with three distinct historical reconstructions of land‐use and land‐use change (LULUC) to quantify LULUC emissions and uncertainty that includes the integrated effects of not only climate and CO₂ but also N. The modeled global average emissions including N dynamics for the 1980s, 1990s, and 2000–2005 were 1.8 ± 0.2, 1.7 ± 0.2, and 1.4 ± 0.2 GtC yr^?1, respectively, (mean and range across LULUC data sets). The emissions from tropics were 0.8 ± 0.2, 0.8 ± 0.2, and 0.7 ± 0.3 GtC yr^?1, and the non tropics were 1.1 ± 0.5, 0.9 ± 0.2, and 0.7 ± 0.1 GtC yr^?1. Compared to previous studies that did not include N dynamics, modeled net LULUC emissions were higher, particularly in the non tropics. In the model, N limitation reduces regrowth rates of vegetation in temperate areas resulting in higher net emissions. Our results indicate that exclusion of N dynamics leads to an underestimation of LULUC emissions by around 70% in the non tropics, 10% in the tropics, and 40% globally in the 1990s. The differences due to inclusion/exclusion of the N cycle of 0.1 GtC yr^?1 in the tropics, 0.6 GtC yr^?1 in the non tropics, and 0.7 GtC yr^?1 globally (mean across land‐cover data sets) in the 1990s were greater than differences due to the land‐cover data in the non tropics and globally (0.2 GtC yr^?1). While land‐cover information is improving with satellite and inventory data, this study indicates the importance of accounting for different processes, in particular the N cycle.     

Carbon stocks,sequestration, and emissions of wetlands in south eastern Australia

Nontidal wetlands are estimated to contribute significantly to the soil carbon pool across the globe. However, our understanding of the occurrence and variability of carbon storage between wetland types and across regions represents a major impediment to the ability of nations to include wetlands in greenhouse gas inventories and carbon offset initiatives. We performed a large‐scale survey of nontidal wetland soil carbon stocks and accretion rates from the state of Victoria in south‐eastern Australia—a region spanning 237,000 km² and containing >35,000 temperate, alpine, and semi‐arid wetlands. From an analysis of >1,600 samples across 103 wetlands, we found that alpine wetlands had the highest carbon stocks (290 ± 180 Mg C_org ha^?1), while permanent open freshwater wetlands and saline wetlands had the lowest carbon stocks (110 ± 120 and 60 ± 50 Mg C_org ha^?1, respectively). Permanent open freshwater sites sequestered on average three times more carbon per year over the last century than shallow freshwater marshes (2.50 ± 0.44 and 0.79 ± 0.45 Mg C_org ha^?1 year^?1, respectively). Using this data, we estimate that wetlands in Victoria have a soil carbon stock in the upper 1 m of 68 million tons of C_org, with an annual soil carbon sequestration rate of 3 million tons of CO₂ eq. year^?1—equivalent to the annual emissions of about 3% of the state's population. Since European settlement (~1834), drainage and loss of 260,530 ha of wetlands may have released between 20 and 75 million tons CO₂ equivalents (based on 27%–90% of soil carbon converted to CO₂). Overall, we show that despite substantial spatial variability within wetland types, some wetland types differ in their carbon stocks and sequestration rates. The duration of water inundation, plant community composition, and allochthonous carbon inputs likely play an important role in influencing variation in carbon storage.     

Attribution of seasonal leaf area index trends in the northern latitudes with “optimally” integrated ecosystem models

Significant increases in remotely sensed vegetation indices in the northern latitudes since the 1980s have been detected and attributed at annual and growing season scales. However, we presently lack a systematic understanding of how vegetation responds to asymmetric seasonal environmental changes. In this study, we first investigated trends in the seasonal mean leaf area index (LAI) at northern latitudes (north of 30°N) between 1982 and 2009 using three remotely sensed long‐term LAI data sets. The most significant LAI increases occurred in summer (0.009 m² m^?2 year^?1, p < .01), followed by autumn (0.005 m² m^?2 year^?1, p < .01) and spring (0.003 m² m^?2 year^?1, p < .01). We then quantified the contribution of elevating atmospheric CO₂ concentration (eCO₂), climate change, nitrogen deposition, and land cover change to seasonal LAI increases based on factorial simulations from 10 state‐of‐the‐art ecosystem models. Unlike previous studies that used multimodel ensemble mean (MME), we used the Bayesian model averaging (BMA) to optimize the integration of model ensemble. The optimally integrated ensemble LAI changes are significantly closer to the observed seasonal LAI changes than the traditional MME results. The BMA factorial simulations suggest that eCO₂ provides the greatest contribution to increasing LAI trends in all seasons (0.003–0.007 m² m^?2 year^?1), and is the main factor driving asymmetric seasonal LAI trends. Climate change controls the spatial pattern of seasonal LAI trends and dominates the increase in seasonal LAI in the northern high latitudes. The effects of nitrogen deposition and land use change are relatively small in all seasons (around 0.0002 m² m^?2 year^?1 and 0.0001–0.001 m² m^?2 year^?1, respectively). Our analysis of the seasonal LAI responses to the interactions between seasonal changes in environmental factors offers a new perspective on the response of global vegetation to environmental changes.