Fire is an important global ecological process that influences the distribution of biomes, with consequences for carbon, water, and energy budgets. Therefore it is impossible to appropriately model ...the history and future of the terrestrial ecosystems and the climate system without including fire. This study incorporates the process-based prognostic fire module SPITFIRE into the global vegetation model ORCHIDEE, which was then used to simulate burned area over the 20th century. Special attention was paid to the evaluation of other fire regime indicators such as seasonality, fire size and fire length, next to burned area. For 2001–2006, the simulated global spatial extent of fire agrees well with that given by satellite-derived burned area data sets (L3JRC, GLOBCARBON, GFED3.1), and 76–92% of the global burned area is simulated as collocated between the model and observation, depending on which data set is used for comparison. The simulated global mean annual burned area is 346 Mha yr−1, which falls within the range of 287–384 Mha yr−1 as given by the three observation data sets; and is close to the 344 Mha yr−1 by the GFED3.1 data when crop fires are excluded. The simulated long-term trend and variation of burned area agree best with the observation data in regions where fire is mainly driven by climate variation, such as boreal Russia (1930–2009), along with Canada and US Alaska (1950–2009). At the global scale, the simulated decadal fire variation over the 20th century is only in moderate agreement with the historical reconstruction, possibly because of the uncertainties of past estimates, and because land-use change fires and fire suppression are not explicitly included in the model. Over the globe, the size of large fires (the 95th quantile fire size) is underestimated by the model for the regions of high fire frequency, compared with fire patch data as reconstructed from MODIS 500 m burned area data. Two case studies of fire size distribution in Canada and US Alaska, and southern Africa indicate that both number and size of large fires are underestimated, which could be related with short fire patch length and low daily fire size. Future efforts should be directed towards building consistent spatial observation data sets for key parameters of the model in order to constrain the model error at each key step of the fire modelling.
Future climate warming is expected to enhance plant growth in temperate ecosystems and to increase carbon sequestration. But although severe regional heatwaves may become more frequent in a changing ...climate, their impact on terrestrial carbon cycling is unclear. Here we report measurements of ecosystem carbon dioxide fluxes, remotely sensed radiation absorbed by plants, and country-level crop yields taken during the European heatwave in 2003. We use a terrestrial biosphere simulation model to assess continental-scale changes in primary productivity during 2003, and their consequences for the net carbon balance. We estimate a 30 per cent reduction in gross primary productivity over Europe, which resulted in a strong anomalous net source of carbon dioxide (0.5 Pg C yr-1) to the atmosphere and reversed the effect of four years of net ecosystem carbon sequestration. Our results suggest that productivity reduction in eastern and western Europe can be explained by rainfall deficit and extreme summer heat, respectively. We also find that ecosystem respiration decreased together with gross primary productivity, rather than accelerating with the temperature rise. Model results, corroborated by historical records of crop yields, suggest that such a reduction in Europe's primary productivity is unprecedented during the last century. An increase in future drought events could turn temperate ecosystems into carbon sources, contributing to positive carbon-climate feedbacks already anticipated in the tropics and at high latitudes.
Celotno besedilo
Dostopno za:
DOBA, IJS, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
The European CARBOEUROPE/FLUXNET monitoring sites, spatial remote sensing observations via the EOS-MODIS sensor and ecosystem modelling provide independent and complementary views on the effect of ...the 2003 heatwave on the European biosphere's productivity and carbon balance. In our analysis, these data streams consistently demonstrate a strong negative anomaly of the primary productivity during the summer of 2003. FLUXNET eddy-covariance data indicate that the drop in productivity was not primarily caused by high temperatures ('heat stress') but rather by limitation of water (drought stress) and that, contrary to the classical expectation about a heat wave, not only gross primary productivity but also ecosystem respiration declined by up to more than to 80 gC m⁻² month⁻¹. Anomalies of carbon and water fluxes were strongly correlated. While there are large between-site differences in water-use efficiency (WUE, 1-6 kg C kg⁻¹ H₂O) here defined as gross carbon uptake divided by evapotranspiration (WUE=GPP/ET), the year-to-year changes in WUE were small (<1 g kg⁻¹) and quite similar for most sites (i.e. WUE decreased during the year of the heatwave). Remote sensing data from MODIS and AVHRR both indicate a strong negative anomaly of the fraction of absorbed photosynthetically active radiation in summer 2003, at more than five standard deviations of the previous years. The spatial differentiation of this anomaly follows climatic and land-use patterns: Largest anomalies occur in the centre of the meteorological anomaly (central Western Europe) and in areas dominated by crops or grassland. A preliminary model intercomparison along a gradient from data-oriented models to process-oriented models indicates that all approaches are similarly describing the spatial pattern of ecosystem sensitivity to the climatic 2003 event with major exceptions in the Alps and parts of Eastern Europe, but differed with respect to their interannual variability.
The global carbon budget 1959–2011 Le Quéré, C; Andres, R. J; Boden, T ...
Earth system science data,
05/2013, Letnik:
5, Številka:
1
Journal Article
Recenzirano
Odprti dostop
Accurate assessments of anthropogenic carbon dioxide (CO sub(2)) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global ...carbon cycle, support the climate policy process, and project future climate change. Present-day analysis requires the combination of a range of data, algorithms, statistics and model estimates and their interpretation by a broad scientific community. Here we describe datasets and a methodology developed by the global carbon cycle science community to quantify all major components of the global carbon budget, including their uncertainties. We discuss changes compared to previous estimates, consistency within and among components, and methodology and data limitations. CO sub(2) emissions from fossil fuel combustion and cement production (E sub(FF)) are based on energy statistics, while emissions from Land-Use Change (E sub(LUC)), including deforestation, are based on combined evidence from land cover change data, fire activity in regions undergoing deforestation, and models. The global atmospheric CO sub(2) concentration is measured directly and its rate of growth (G sub(ATM)) is computed from the concentration. The mean ocean CO sub(2) sink (S sub(OCEAN)) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. Finally, the global residual terrestrial CO sub(2) sink (S sub(LAND)) is estimated by the difference of the other terms. For the last decade available (2002-2011), E sub(FF) was 8.3 plus or minus 0.4 PgC yr super(-1), E sub(LUC) 1.0 plus or minus 0.5 PgC yr super(-1), G sub(ATM) 4.3 plus or minus 0.1PgC yr super(-1), S sub(OCEAN) 2.5 plus or minus 0.5 PgC yr super(-1), and S sub(LAND) 2.6 plus or minus 0.8 PgC yr super(-1). For year 2011 alone, E sub(FF) was 9.5 plus or minus 0.5 PgC yr super(-1), 3.0 percent above 2010, reflecting a continued trend in these emissions; E sub(LUC) was 0.9 plus or minus 0.5 PgC yr super(-1), approximately constant throughout the decade; G sub(ATM) was 3.6 plus or minus 0.2 PgC yr super(-1), S sub(OCEAN) was 2.7 plus or minus 0.5 PgC yr super(-1), and S sub(LAND) was 4.1 plus or minus 0.9 PgC yr super(-1). G sub(ATM) was low in 2011 compared to the 2002-2011 average because of a high uptake by the land probably in response to natural climate variability associated to La Nina conditions in the Pacific Ocean. The global atmospheric CO sub(2) concentration reached 391.31 plus or minus 0.13 ppm at the end of year 2011. We estimate that E sub(FF) will have increased by 2.6% (1.9-3.5%) in 2012 based on projections of gross world product and recent changes in the carbon intensity of the economy. All uncertainties are reported as plus or minus 1 sigma (68% confidence assuming Gaussian error distributions that the real value lies within the given interval), reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. This paper is intended to provide a baseline to keep track of annual carbon budgets in the future. All data presented here can be downloaded from the Carbon Dioxide Information Analysis Center (doi:10.3334/CDIAC/GCP_V2013).
This study describes how management of grasslands is included in the Organizing Carbon and Hydrology in Dynamic Ecosystems (ORCHIDEE) process-based ecosystem model designed for large-scale ...applications, and how management affects modeled grassland-atmosphere CO sub(2) fluxes. The new model, ORCHIDEE-GM (grassland management) is enabled with a management module inspired from a grassland model (PaSim, version 5.0), with two grassland management practices being considered, cutting and grazing. The evaluation of the results from ORCHIDEE compared with those of ORCHIDEE-GM at 11 European sites, equipped with eddy covariance and biometric measurements, shows that ORCHIDEE-GM can realistically capture the cut-induced seasonal variation in biometric variables (LAI: leaf area index; AGB: aboveground biomass) and in CO sub(2) fluxes (GPP: gross primary productivity; TER: total ecosystem respiration; and NEE: net ecosystem exchange). However, improvements at grazing sites are only marginal in ORCHIDEE-GM due to the difficulty in accounting for continuous grazing disturbance and its induced complex animal-vegetation interactions. Both NEE and GPP on monthly to annual timescales can be better simulated in ORCHIDEE-GM than in ORCHIDEE without management. For annual CO sub(2) fluxes, the NEE bias and RMSE (root mean square error) in ORCHIDEE-GM are reduced by 53% and 20%, respectively, compared to ORCHIDEE. ORCHIDEE-GM is capable of modeling the net carbon balance (NBP) of managed temperate grasslands (37 plus or minus 30 gC m super(-2) yr super(-1) (P < 0.01) over the 11 sites) because the management module contains provisions to simulate the carbon fluxes of forage yield, herbage consumption, animal respiration and methane emissions.
The response of crops to changing climate and atmospheric CO2 concentration (CO2) could have large effects on food production, and impact carbon, water, and energy fluxes, causing feedbacks to the ...climate. To simulate the response of temperate crops to changing climate and CO2, which accounts for the specific phenology of crops mediated by management practice, we describe here the development of a process-oriented terrestrial biogeochemical model named ORCHIDEE-CROP (v0), which integrates a generic crop phenology and harvest module, and a very simple parameterization of nitrogen fertilization, into the land surface model (LSM) ORCHIDEEv196, in order to simulate biophysical and biochemical interactions in croplands, as well as plant productivity and harvested yield. The model is applicable for a range of temperate crops, but is tested here using maize and winter wheat, with the phenological parameterizations of two European varieties originating from the STICS agronomical model. We evaluate the ORCHIDEE-CROP (v0) model against eddy covariance and biometric measurements at seven winter wheat and maize sites in Europe. The specific ecosystem variables used in the evaluation are CO2 fluxes (net ecosystem exchange, NEE), latent heat, and sensible heat fluxes. Additional measurements of leaf area index (LAI) and aboveground biomass and yield are used as well. Evaluation results revealed that ORCHIDEE-CROP (v0) reproduced the observed timing of crop development stages and the amplitude of the LAI changes. This is in contrast to ORCHIDEEv196 where, by default, crops have the same phenology as grass. A halving of the root mean square error for LAI from 2.38 ± 0.77 to 1.08 ± 0.34 m2 m−2 was obtained when ORCHIDEEv196 and ORCHIDEE-CROP (v0) were compared across the seven study sites. Improved crop phenology and carbon allocation led to a good match between modeled and observed aboveground biomass (with a normalized root mean squared error (NRMSE) of 11.0–54.2 %), crop yield, daily carbon and energy fluxes (with a NRMSE of ∼ 9.0–20.1 and ∼ 9.4–22.3 % for NEE), and sensible and latent heat fluxes. The simulated yields for winter wheat and maize from ORCHIDEE-CROP (v0) showed a good match with the simulated results from STICS for three sites with available crop yield observations, where the average NRMSE was ∼ 8.8 %. The model data misfit for energy fluxes were within the uncertainties of the measurements, which themselves showed an incomplete energy balance closure within the range 80.6–86.3 %. The remaining discrepancies between the modeled and observed LAI and other variables at specific sites were partly attributable to unrealistic representations of management events by the model. ORCHIDEE-CROP (v0) has the ability to capture the spatial gradients of carbon and energy-related variables, such as gross primary productivity, NEE, and sensible and latent heat fluxes across the sites in Europe, which is an important requirement for future spatially explicit simulations. Further improvement of the model, with an explicit parameterization of nutritional dynamics and management, is expected to improve its predictive ability to simulate croplands in an Earth system model.
Soil moisture availability is important in regulating photosynthesis and controlling land surface-climate feedbacks at both the local and global scale. Recently, global remote-sensing datasets for ...soil moisture have become available. In this paper we assess the possibility of using remotely sensed soil moisture – AMSR-E (LPRM) – to similate soil moisture dynamics of the process-based vegetation model ORCHIDEE by evaluating the correspondence between these two products using both correlation and autocorrelation analyses. We find that the soil moisture product of AMSR-E (LPRM) and the simulated soil moisture in ORCHIDEE correlate well in space and time, in particular when considering the root zone soil moisture of ORCHIDEE. However, the root zone soil moisture in ORCHIDEE has on average a higher temporal autocorrelation relative to AMSR-E (LPRM) and in situ measurements. This may be due to the different vertical depth of the two products – AMSR-E (LPRM) at the 2–5 cm surface depth and ORCHIDEE at the root zone (max. 2 m) depth – to uncertainty in precipitation forcing in ORCHIDEE, and to the fact that the structure of ORCHIDEE consists of a single-layer deep soil, which does not allow simulation of the proper cascade of time scales that characterize soil drying after each rain event. We conclude that assimilating soil moisture, using AMSR-E (LPRM) in a land surface model like ORCHIDEE with an improved hydrological model of more than one soil layer, may significantly improve the soil moisture dynamics, which could lead to improved CO2 and energy flux predictions.
This paper describes a new model for the calculation of daily, high-resolution (up to 1 km) fire emissions, developed in the framework of the APIFLAME (Analysis and Prediction of the Impact of Fires ...on Air quality ModEling) project. The methodology relies on the classical approach, multiplying the burned area by the fuel load consumed and the emission factors specific to the vegetation burned. Emissions can be calculated on any user-specified domain, horizontal grid, and list of trace gases and aerosols, providing input information on the burned area (location, extent), and emission factors of the targeted species are available. The applicability to high spatial resolutions and the flexibility to different input data (including vegetation classifications) and domains are the main strength of the proposed algorithm. The modification of the default values and databases proposed does not require any change in the core of the model. The code may be used for the calculation of global or regional inventories. However, it has been developed and tested more specifically for Europe and the Mediterranean area. A regional analysis of fire activity and the resulting emissions in this region is provided. The burning season extends from June to October in most regions, with generally small but frequent fires in eastern Europe, western Russia, Ukraine and Turkey, and large events in the Mediterranean area. The resulting emissions represent a significant fraction of the total yearly emissions (on average amounting to ~ 30% of anthropogenic emissions for PM2.5, ~ 20% for CO). The uncertainty regarding the daily carbon emissions is estimated at ~ 100% based on an ensemble analysis. Considering the large uncertainties regarding emission factors, the potential error on the emissions for the various pollutants is even larger. Comparisons with other widely used emission inventories show good correlations but discrepancies of a factor of 2–4 in the amplitude of the emissions, our results being generally on the higher end.
Climate change is leading to a disproportionately large warming in the high northern latitudes, but the magnitude and sign of the future carbon balance of the Arctic are highly uncertain. Using 40 ...terrestrial biosphere models for the Alaskan Arctic from four recent model intercomparison projects – NACP (North American Carbon Program) site and regional syntheses, TRENDY (Trends in net land atmosphere carbon exchanges), and WETCHIMP (Wetland and Wetland CH4 Inter-comparison of Models Project) – we provide a baseline of terrestrial carbon cycle uncertainty, defined as the multi-model standard deviation (σ) for each quantity that follows. Mean annual absolute uncertainty was largest for soil carbon (14.0 ± 9.2 kg C m−2), then gross primary production (GPP) (0.22 ± 0.50 kg C m−2 yr−1), ecosystem respiration (Re) (0.23 ± 0.38 kg C m−2 yr−1), net primary production (NPP) (0.14 ± 0.33 kg C m−2 yr−1), autotrophic respiration (Ra) (0.09 ± 0.20 kg C m−2 yr−1), heterotrophic respiration (Rh) (0.14 ± 0.20 kg C m−2 yr−1), net ecosystem exchange (NEE) (−0.01 ± 0.19 kg C m−2 yr−1), and CH4 flux (2.52 ± 4.02 g CH4 m−2 yr−1). There were no consistent spatial patterns in the larger Alaskan Arctic and boreal regional carbon stocks and fluxes, with some models showing NEE for Alaska as a strong carbon sink, others as a strong carbon source, while still others as carbon neutral. Finally, AmeriFlux data are used at two sites in the Alaskan Arctic to evaluate the regional patterns; observed seasonal NEE was captured within multi-model uncertainty. This assessment of carbon cycle uncertainties may be used as a baseline for the improvement of experimental and modeling activities, as well as a reference for future trajectories in carbon cycling with climate change in the Alaskan Arctic and larger boreal region.