Engineering Papers⌕ Search

Engineering topics

Ciais, Philippe

Publications and source records attributed to Ciais, Philippe.

At least 73 records · Page 4

Global Carbon Budget 2017

Here an accurate assessment of anthropogenic carbon dioxide (CO 2 ) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the global carbon budget – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. CO 2 emissions from fossil fuels and industry ( E FF ) are based on energy statistics and cement production data, respectively, while emissions from land-use change ( E LUC ), mainly deforestation, are based on land-cover change data and bookkeeping models. The global atmospheric CO 2 concentration is measured directly and its rate of growth ( G ATM ) is computed from the annual changes in concentration. The ocean CO 2 sink ( S OCEAN ) and terrestrial CO 2 sink ( S LAND ) are estimated with global process models constrained by observations. The resulting carbon budget imbalance ( B IM ), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1 σ . For the last decade available (2007–2016), E FF was 9.4 ± 0.5 GtC yr –1 , E LUC 1.3 ± 0.7 GtC yr –1 , G ATM 4.7 ± 0.1 GtC yr –1 , S OCEAN 2.4 ± 0.5 GtC yr –1 , and S LAND 3.0 ± 0.8 GtC yr –1 , with a budget imbalance B IM of 0.6 GtC yr –1 indicating overestimated emissions and/or underestimated sinks. For year 2016 alone, the growth in E FF was approximately zero and emissions remained at 9.9 ±\ 0.5 GtC yr –1 . Also for 2016, E LUC was 1.3 ± .7 GtC yr –1 , G ATM was 6.1 ± 0.2 GtC yr –1 , S OCEAN was 2.6 ± 0.5 GtC yr –1 , and S LAND was 2.7 ± 1.0 GtC yr –1 , with a small B IM of –0.3 GtC. G ATM continued to be higher in 2016 compared to the past decade (2007–2016), reflecting in part the high fossil emissions and the small S LAND consistent with El Niño conditions. The global atmospheric CO 2 concentration reached 402.8 ± 0.1 ppm averaged over 2016. For 2017, preliminary data for the first 6–9 months indicate a renewed growth in E FF of +2.0 % (range of 0.8 to 3.0 %) based on national emissions projections for China, USA, and India, and projections of gross domestic product (GDP) corrected for recent changes in the carbon intensity of the economy for the rest of the world. This living data update documents changes in the methods and data sets used in this new global carbon budget compared with previous publications of this data set.

54 ENVIRONMENTAL SCIENCES↗

Missing pieces to modeling the Arctic-Boreal puzzle

NASA has launched the decade-long Arctic-Boreal Vulnerability Experiment (ABoVE). While the initial phases focus on field and airborne data collection, early integration with modeling activities is important to benefit future modeling syntheses. We compiled feedback from ecosystem modeling teams on key data needs, which encompass carbon biogeochemistry, vegetation, permafrost, hydrology, and disturbance dynamics. A suite of variables was identified as part of this activity with a critical requirement that they are collected concurrently and representatively over space and time. Furthermore, individual projects in ABoVE may not capture all these needs, and thus there is both demand and opportunity for the augmentation of field observations, and synthesis of the observations that are collected, to ensure that science questions and integrated modeling activities are successfully implemented.

54 ENVIRONMENTAL SCIENCES↗

ORCHIDEE-PEAT (revision 4596), a model for northern peatland CO 2 , water, and energy fluxes on daily to annual scales

Peatlands store substantial amounts of carbon and are vulnerable to climate change. We present a modified version of the Organising Carbon and Hydrology In Dynamic Ecosystems (ORCHIDEE) land surface model for simulating the hydrology, surface energy, and CO 2 fluxes of peatlands on daily to annual timescales. The model includes a separate soil tile in each 0.5° grid cell, defined from a global peatland map and identified with peat-specific soil hydraulic properties. Runoff from non-peat vegetation within a grid cell containing a fraction of peat is routed to this peat soil tile, which maintains shallow water tables. The water table position separates oxic from anoxic decomposition. The model was evaluated against eddy-covariance (EC) observations from 30 northern peatland sites, with the maximum rate of carboxylation ( V cmax ) being optimized at each site. Regarding short-term day-to-day variations, the model performance was good for gross primary production (GPP) ( r 2 = 0.76; Nash–Sutcliffe modeling efficiency, MEF =0.76) and ecosystem respiration (ER, r 2 = 0.78, MEF = 0.75), with lesser accuracy for latent heat fluxes (LE, r 2 2 exchange (NEE, r 2 = 0.38, MEF = 0.26). Seasonal variations in GPP, ER, NEE, and energy fluxes on monthly scales showed moderate to high r 2 values (0.57–0.86). For spatial across-site gradients of annual mean GPP, ER, NEE, and LE, r 2 values of 0.93, 0.89, 0.27, and 0.71 were achieved, respectively. Water table (WT) variation was not well predicted ( r 2 < 0.1), likely due to the uncertain water input to the peat from surrounding areas. However, the poor performance of WT simulation did not greatly affect predictions of ER and NEE. We found a significant relationship between optimized V cmax and latitude (temperature), which better reflects the spatial gradients of annual NEE than using an average V cmax value

58 GEOSCIENCES↗

Response of Water Use Efficiency to Global Environmental Change Based on Output From Terrestrial Biosphere Models

Here, water use efficiency (WUE), defined as the ratio of gross primary productivity and evapotranspiration at the ecosystem scale, is a critical variable linking the carbon and water cycles. Incorporating a dependency on vapor pressure deficit, apparent underlying WUE (uWUE) provides a better indicator of how terrestrial ecosystems respond to environmental changes than other WUE formulations. Here we used 20th century simulations from four terrestrial biosphere models to develop a novel variance decomposition method. With this method, we attributed variations in apparent uWUE to both the trend and interannual variation of environmental drivers. The secular increase in atmospheric CO 2 explained a clear majority of total variation (66 ± 32%: mean ± one standard deviation), followed by positive trends in nitrogen deposition and climate, as well as a negative trend in land use change. In contrast, interannual variation was mostly driven by interannual climate variability. To analyze the mechanism of the CO 2 effect, we partitioned the apparent uWUE into the transpiration ratio (transpiration over evapotranspiration) and potential uWUE. The relative increase in potential uWUE parallels that of CO 2 , but this direct CO 2 effect was offset by 20 ± 4% by changes in ecosystem structure, that is, leaf area index for different vegetation types. However, the decrease in transpiration due to stomatal closure with rising CO 2 was reduced by 84% by an increase in leaf area index, resulting in small changes in the transpiration ratio. CO 2 concentration thus plays a dominant role in driving apparent uWUE variations over time, but its role differs for the two constituent components: potential uWUE and transpiration.

54 ENVIRONMENTAL SCIENCES↗

Regional contribution to variability and trends of global gross primary productivity

Terrestrial gross primary productivity (GPP) is the largest component of the global carbon cycle and a key process for understanding land ecosystems dynamics. In this study, we used GPP estimates from a combination of eight global biome models participating in the Inter-Sectoral Impact-Model Intercomparison Project phase 2a (ISIMIP2a), the Moderate Resolution Spectroradiometer (MODIS) GPP product, and a data-driven product (Model Tree Ensemble, MTE) to study the spatiotemporal variability of GPP at the regional and global levels. We found the 2000-2010 total global GPP estimated from the model ensemble to be 117±13 Pg C yr-1 (mean ± 1 standard deviation), which was higher than MODIS (112 Pg C yr-1), and close to the MTE (120 Pg C yr-1). The spatial patterns of MODIS, MTE and ISIMIP2a GPP generally agree well, but their temporal trends are different, and the seasonality and inter-annual variability of GPP at the regional and global levels are not completely consistent. For the model ensemble, Tropical Latin America contributes the most to global GPP, Asian regions contribute the most to the global GPP trend, the Northern Hemisphere regions dominate the global GPP seasonal variations, and Oceania is likely the largest contributor to inter-annual variability of global GPP. However, we observed large uncertainties across the eight ISIMIP2a models, which are probably due to the differences in the formulation of underlying photosynthetic processes. The results of this study are useful in understanding the contributions of different regions to global GPP and its spatiotemporal variability, how the model- and observational-based GPP estimates differ from each other in time and space, and the relative strength of the eight models. Our results also highlight the models’ ability to capture the seasonality of GPP that are essential for understanding the inter-annual and seasonal variability of GPP as a major component of the carbon cycle.

Chen, Min↗

Variability and quasi-decadal changes in the methane budget over the period 2000–2012

Following the recent Global Carbon Project (GCP) synthesis of the decadal methane (CH 4 ) budget over 2000–2012 (Saunois et al., 2016), we analyse here the same dataset with a focus on quasi-decadal and inter-annual variability in CH 4 emissions. The GCP dataset integrates results from top-down studies (exploiting atmospheric observations within an atmospheric inverse-modelling framework) and bottom-up models (including process-based models for estimating land surface emissions and atmospheric chemistry), inventories of anthropogenic emissions, and data-driven approaches. The annual global methane emissions from top-down studies, which by construction match the observed methane growth rate within their uncertainties, all show an increase in total methane emissions over the period 2000-2012, but this increase is not linear over the 13 years. Despite differences between individual studies, the mean emission anomaly of the top-down ensemble shows no significant trend in total methane emissions over the period 2000-2006, during the plateau of atmospheric methane mole fractions, and also over the period 2008-2012, during the renewed atmospheric methane increase. However, the top-down ensemble mean produces an emission shift between 2006 and 2008, leading to 22 [16-32]Tg CH 4 yr -1 higher methane emissions over the period 2008-2012 compared to 2002-2006. This emission increase mostly originated from the tropics, with a smaller contribution from mid-latitudes and no significant change from boreal regions. The regional contributions remain uncertain in top-down studies. Tropical South America and South and East Asia seem to contribute the most to the emission increase in the tropics. However, these two regions have only limited atmospheric measurements and remain therefore poorly constrained. The sectorial partitioning of this emission increase between the periods 2002-2006 and 2008-2012 differs from one atmospheric inversion study to another. However, all top-down studies suggest smaller changes in fossil fuel emissions (from oil, gas, and coal industries) compared to the mean of the bottom-up inventories included in this study. This difference is partly driven by a smaller emission change in China from the top-down studies compared to the estimate in the Emission Database for Global Atmospheric Research (EDGARv4.2) inventory, which should be revised to smaller values in a near future. We apply isotopic signatures to the emission changes estimated for individual studies based on five emission sectors and find that for six individual top-down studies (out of eight) the average isotopic signature of the emission changes is not consistent with the observed change in atmospheric 13 CH 4 . However, the partitioning in emission change derived from the ensemble mean is consistent with this isotopic constraint. At the global scale, the top-down ensemble mean suggests that the dominant contribution to the resumed atmospheric CH 4 growth after 2006 comes from microbial sources (more from agriculture and waste sectors than from natural wetlands), with an uncertain but smaller contribution from fossil CH 4 emissions. In addition, a decrease in biomass burning emissions (in agreement with the biomass burning emission databases) makes the balance of sources consistent with atmospheric 13 CH 4 observations. In most of the top-down studies included here, OH concentrations are considered constant over the years (seasonal variations but without any inter-annual variability). As a result, the methane loss (in particular through OH oxidation) varies mainly through the change in methane concentrations and not its oxidants. For these reasons, changes in the methane loss could not be properly investigated in this study, although it may play a significant role in the recent atmospheric methane changes as briefly discussed at the end of the paper.

58 GEOSCIENCES↗

Variability and Quasi-Decadal Changes in the Methane Budget over the Period 2000-2012

Following the recent Global Carbon Project (GCP) synthesis of the decadal methane (CH4) budget over 2000- 2012, we analyse here the same dataset with a focus on quasi-decadal and inter-annual variability in CH4 emissions. The GCP dataset integrates results from top-down studies (exploiting atmospheric observations within an atmospheric inverse-modelling framework) and bottom-up models (including process-based models for estimating land surface emissions and atmospheric chemistry), inventories of anthropogenic emissions, and data-driven approaches. The annual global methane emissions from top-down studies, which by construction match the observed methane growth rate within their uncertainties, all show an increase in total methane emissions over the period 2000-2012, but this increase is not linear over the 13 years. Despite differences between individual studies, the mean emission anomaly of the top-down ensemble shows no significant trend in total methane emissions over the period 2000-2006, during the plateau of atmospheric methane mole fractions, and also over the period 2008-2012, during the renewed atmospheric methane increase. However, the top-down ensemble mean produces an emission shift between 2006 and 2008, leading to 22 [16-32] Tg CH4 yr(exp -1) higher methane emissions over the period 2008-2012 compared to 2002-2006. This emission increase mostly originated from the tropics, with a smaller contribution from mid-latitudes and no significant change from boreal regions. The regional contributions remain uncertain in top-down studies. Tropical South America and South and East Asia seem to contribute the most to the emission increase in the tropics. However, these two regions have only limited atmospheric measurements and remain therefore poorly constrained. The sectorial partitioning of this emission increase between the periods 2002-2006 and 2008-2012 differs from one atmospheric inversion study to another. However, all top-down studies suggest smaller changes in fossil fuel emissions (from oil, gas, and coal industries) compared to the mean of the bottom-up inventories included in this study. This difference is partly driven by a smaller emission change in China from the top-down studies compared to the estimate in the Emission Database for Global Atmospheric Research (EDGARv4.2) inventory, which should be revised to smaller values in a near future. We apply isotopic signatures to the emission changes estimated for individual studies based on five emission sectors and find that for six individual top-down studies (out of eight) the average isotopic signature of the emission changes is not consistent with the observed change in atmospheric 13CH4. However, the partitioning in emission change derived from the ensemble mean is consistent with this isotopic constraint. At the global scale, the top-down ensemble mean suggests that the dominant contribution to the resumed atmospheric CH4 growth after 2006 comes from microbial sources (more from agriculture and waste sectors than from natural wetlands), with an uncertain but smaller contribution from fossil CH4 emissions. In addition, a decrease in biomass burning emissions (in agreement with the biomass burning emission databases) makes the balance of sources consistent with atmospheric 13CH4 observations. In most of the top-down studies included here, OH concentrations are considered constant over the years (seasonal variations but without any inter-annual variability). As a result, the methane loss (in particular through OH oxidation) varies mainly through the change in methane concentrations and not its oxidants. For these reasons, changes in the methane loss could not be properly investigated in this study, although it may play a significant role in the recent atmospheric methane changes as briefly discussed at the end of the paper.

Saunois, Marielle↗

Global wetland contribution to 2000–2012 atmospheric methane growth rate dynamics

Increasing atmospheric methane (CH 4 ) concentrations have contributed to approximately 20% of anthropogenic climate change. Despite the importance of CH 4 as a greenhouse gas, its atmospheric growth rate and dynamics over the past two decades, which include a stabilization period (1999–2006), followed by renewed growth starting in 2007, remain poorly understood. We provide an updated estimate of CH 4 emissions from wetlands, the largest natural global CH 4 source, for 2000–2012 using an ensemble of biogeochemical models constrained with remote sensing surface inundation and inventory-based wetland area data. Between 2000–2012, boreal wetland CH 4 emissions increased by 1.2 Tg yr –1 (–0.2–3.5 Tg yr –1 ), tropical emissions decreased by 0.9 Tg yr –1 (–3.2–1.1 Tg yr –1 ), yet globally, emissions remained unchanged at 184 ± 22 Tg yr –1 . Changing air temperature was responsible for increasing high-latitude emissions whereas declines in low-latitude wetland area decreased tropical emissions; both dynamics are consistent with features of predicted centennial-scale climate change impacts on wetland CH 4 emissions. Despite uncertainties in wetland area mapping, our study shows that global wetland CH 4 emissions have not contributed significantly to the period of renewed atmospheric CH 4 growth, and is consistent with findings from studies that indicate some combination of increasing fossil fuel and agriculture-related CH 4 emissions, and a decrease in the atmospheric oxidative sink.

54 ENVIRONMENTAL SCIENCES↗

Global Wetland Contribution to 2000-2012 Atmospheric Methane Growth Rate Dynamics

Increasing atmospheric methane (CH4) concentrations have contributed to approximately 20% of anthropogenic climate change. Despite the importance of CH4 as a greenhouse gas, its atmospheric growth rate and dynamics over the past two decades, which include a stabilization period (1999-2006), followed by renewed growth starting in 2007, remain poorly understood. We provide an updated estimate of CH4 emissions from wetlands, the largest natural global CH4 source, for 2000-2012 using an ensemble of biogeochemical models constrained with remote sensing surface inundation and inventory-based wetland area data. Between 2000-2012, boreal wetland CH4 emissions increased by 1.2 Tg yr(sup -1) (-0.2-3.5 Tg yr(sup -1), tropical emissions decreased by 0.9 Tg yr(sup -1) (-3.2-1.1 Tg yr(sup -1), yet globally, emissions remained unchanged at 184 +/- 22 Tg yr(sup -1). Changing air temperature was responsible for increasing high-latitude emissions whereas declines in low-latitude wetland area decreased tropical emissions; both dynamics are consistent with features of predicted centennial-scale climate change impacts on wetland CH4 emissions. Despite uncertainties in wetland area mapping, our study shows that global wetland CH4 emissions have not contributed significantly to the period of renewed atmospheric CH4 growth, and is consistent with findings from studies that indicate some combination of increasing fossil fuel and agriculture-related CH4 emissions, and a decrease in the atmospheric oxidative sink.

Poulter, Benjamin↗

Temperature Increase Reduces Global Yields of Major Crops in Four Independent Estimates

Wheat, rice, maize, and soybean provide two-thirds of human caloric intake. Assessing the impact of global temperature increase on production of these crops is therefore critical to maintaining global food supply, but different studies have yielded different results. Here, we investigated the impacts of temperature on yields of the four crops by compiling extensive published results from four analytical methods: global grid-based and local point-based models, statistical regressions, and field-warming experiments. Results from the different methods consistently showed negative temperature impacts on crop yield at the global scale, generally underpinned by similar impacts at country and site scales. Without CO2 fertilization, effective adaptation, and genetic improvement, each degree-Celsius increase in global mean temperature would, on average, reduce global yields of wheat by 6.0%, rice by 3.2%, maize by 7.4%, and soybean by 3.1%. Results are highly heterogeneous across crops and geographical areas, with some positive impact estimates. Multi-method analyses improved the confidence in assessments of future climate impacts on global major crops and suggest crop- and region-specific adaptation strategies to ensure food security for an increasing world population.

Zhao, Chuang↗

Diagnosing phosphorus limitations in natural terrestrial ecosystems in carbon cycle models

Abstract Most of the Earth System Models ( ESMs ) project increases in net primary productivity ( NPP ) and terrestrial carbon (C) storage during the 21st century. Despite empirical evidence that limited availability of phosphorus (P) may limit the response of NPP to increasing atmospheric CO 2 , none of the ESMs used in the previous Intergovernmental Panel on Climate Change assessment accounted for P limitation. We diagnosed from ESM simulations the amount of P need to support increases in carbon uptake by natural ecosystems using two approaches: the demand derived from (1) changes in C stocks and (2) changes in NPP . The C stock‐based additional P demand was estimated to range between −31 and 193 Tg P and between −89 and 262 Tg P for Representative Concentration Pathway ( RCP ) 2.6 and RCP8.5 , respectively, with negative values indicating a P surplus. The NPP ‐based demand, which takes ecosystem P recycling into account, results in a significantly higher P demand of 648–1606 Tg P for RCP2 .6 and 924–2110 Tg P for RCP8 .5. We found that the P demand is sensitive to the turnover of P in decomposing plant material, explaining the large differences between the NPP ‐based demand and C stock‐based demand. The discrepancy between diagnosed P demand and actual P availability (potential P deficit) depends mainly on the assumptions about availability of the different soil P forms. Overall, future P limitation strongly depends on both soil P availability and P recycling on ecosystem scale.

54 ENVIRONMENTAL SCIENCES↗

The Critical Role of the Routing Scheme in Simulating Peak River Discharge in Global Hydrological Models

Global hydrological models (GHMs) have been applied to assess global flood hazards, but their capacity to capture the timing and amplitude of peak river discharge which is crucial in flood simulations has traditionally not been the focus of examination. Here we evaluate to what degree the choice of river routing scheme affects simulations of peak discharge and may help to provide better agreement with observations. To this end we use runoff and discharge simulations of nine GHMs forced by observational climate data (1971-2010) within the ISIMIP2a (Inter-Sectoral Impact Model Intercomparison Project phase 2a) project. The runoff simulations were used as input for the global river routing model CaMa-Flood (Catchment-based Macro-scale Floodplain). The simulated daily discharge was compared to the discharge generated by each GHM using its native river routing scheme. For each GHM both versions of simulated discharge were compared to monthly and daily discharge observations from 1701 GRDC (Global Runoff Data Centre) stations as a benchmark. CaMa-Flood routing shows a general reduction of peak river discharge and a delay of about two to three weeks in its occurrence, likely induced by the buffering capacity of floodplain reservoirs. For a majority of river basins, discharge produced by CaMa-Flood resulted in a better agreement with observations. In particular, maximum daily discharge was adjusted, with a multi-model averaged reduction in bias over about two-thirds of the analysed basin area. The increase in agreement was obtained in both managed and near-natural basins. Overall, this study demonstrates the importance of routing scheme choice in peak discharge simulation, where CaMa-Flood routing accounts for floodplain storage and backwater effects that are not represented in most GHMs. Our study provides important hints that an explicit parameterisation of these processes may be essential in future impact studies.

peak river discharge↗

Global land carbon sink response to temperature and precipitation varies with ENSO phase

Climate variability associated with the El Niño-Southern Oscillation (ENSO) and its consequent impacts on land carbon sink interannual variability have been used as a basis for investigating carbon cycle responses to climate variability more broadly, and to inform the sensitivity of the tropical carbon budget to climate change. Past studies have presented opposing views about whether temperature or precipitation is the primary factor driving the response of the land carbon sink to ENSO. We show that the dominant driver varies with ENSO phase. And whereas tropical temperature explains sink dynamics following El Niño conditions (r TG,P = 0.59, p < 0.01), the post La Niña sink is driven largely by tropical precipitation (r PG,T= -0.46, p = 0.04). This finding points to an ENSO-phase-dependent interplay between water availability and temperature in controlling the carbon uptake response to climate variations in tropical ecosystems. Furthermore, we find that none of a suite of ten contemporary terrestrial biosphere models captures these ENSO-phase-dependent responses, highlighting a key uncertainty in modeling climate impacts on the future of the global land carbon sink.

54 ENVIRONMENTAL SCIENCES↗

Converging Climate Sensitivities of European Forests Between Observed Radial Tree Growth and Vegetation Models

The impacts of climate variability and trends on European forests are unevenly distributed across different bioclimatic zones and species. Extreme climate events are also becoming more frequent and it is unknown how they will affect feed backs of CO2 between forest ecosystems and the atmosphere. An improved understanding of species differences at the regional scale of the response of forest productivity to climate variation and extremes is thus important for forecasting forest dynamics. In this study, we evaluate the climate sensitivity of above ground net primary production (NPP) simulated by two dynamic global vegetation models (DGVM; ORCHIDEE and LPJ-wsl) against tree ring width (TRW) observations from about1000 sites distributed across Europe. In both the model simulations and the TRW observations, forests in northern Europe and the Alps respond positively to warmer spring and summer temperature, and their overall temperature sensitivity is larger than that of the soil-moisture-limited forests in central Europe and Mediterranean regions. Compared with TRW observations, simulated NPP from ORCHIDEE and LPJ-wsl appear to be overly sensitive to climatic factors. Our results indicate that the models lack biological processes that control time lags, such as carbohydrate storage and remobilization, that delay the effects of radial growth dynamics to climate. Our study highlights the need for re-evaluating the physiological controls on the climate sensitivity of NPP simulated by DGVMs. In particular, DGVMs could be further enhanced by a more detailed representation of carbon reserves and allocation that control year-to year variation in plant growth.

Zhang, Zhen↗

Climate mitigation from vegetation biophysical feedbacks during the past three decades

The surface air temperature response to vegetation changes has been studied for the extreme case of land-cover change; yet, it has never been quantified for the slow but persistent increase in leaf area index (LAI) observed over the past 30 years (Earth greening). We isolate the fingerprint of increasing LAI on surface air temperature using a coupled land–atmosphere global climate model prescribed with satellite LAI observations. Furthermore, we found that the global greening has slowed down the rise in global land-surface air temperature by 0.09 ± 0.02 °C since 1982. This net cooling effect is the sum of cooling from increased evapotranspiration (70%), changed atmospheric circulation (44%), decreased shortwave transmissivity (21%), and warming from increased longwave air emissivity (-29%) and decreased albedo (-6%). The global cooling originated from the regions where LAI has increased, including boreal Eurasia, Europe, India, northwest Amazonia, and the Sahel. Increasing LAI did not, but, significantly change surface air temperature in eastern North America and East Asia, where the effects of large-scale atmospheric circulation changes mask local vegetation feedbacks. Overall, the sum of biophysical feedbacks related to the greening of the Earth mitigated 12% of global land-surface warming for the past 30 years.

54 ENVIRONMENTAL SCIENCES↗

Weakening temperature control on the interannual variations of spring carbon uptake across northern lands

Ongoing spring warming allows the growing season to begin earlier, enhancing carbon uptake in northern ecosystems. We use 34 years of atmospheric CO 2 concentration measurements at Barrow, Alaska (BRW, 71° N) to show that the interannual relationship between spring temperature and carbon uptake has recently shifted. Here, we use two indicators: the spring zero-crossing date of atmospheric CO 2 (SZC) and the magnitude of CO 2 drawdown between May and June (SCC). The previously reported strong correlation between SZC, SCC and spring land temperature (ST) was found in the first 17 years of measurements, but disappeared in the last 17 years. As a result, the sensitivity of both SZC and SCC to warming decreased. Simulations with an atmospheric transport model coupled to a terrestrial ecosystem model suggest that the weakened interannual correlation of SZC and SCC with ST in the last 17 years is attributable to the declining temperature response of spring net primary productivity (NPP) rather than to changes in heterotrophic respiration or in atmospheric transport patterns. Reduced chilling during dormancy and emerging light limitation are possible mechanisms that may have contributed to the loss of NPP response to ST. These results thus challenge the ‘warmer spring–bigger sink’ mechanism.

54 ENVIRONMENTAL SCIENCES↗

Global Gridded Crop Model Evaluation: Benchmarking, Skills, Deficiencies and Implications.

Crop models are increasingly used to simulate crop yields at the global scale, but so far there is no general framework on how to assess model performance. Here we evaluate the simulation results of 14 global gridded crop modeling groups that have contributed historic crop yield simulations for maize, wheat, rice and soybean to the Global Gridded Crop Model Intercomparison (GGCMI) of the Agricultural Model Intercomparison and Improvement Project (AgMIP). Simulation results are compared to reference data at global, national and grid cell scales and we evaluate model performance with respect to time series correlation, spatial correlation and mean bias. We find that global gridded crop models (GGCMs) show mixed skill in reproducing time series correlations or spatial patterns at the different spatial scales. Generally, maize, wheat and soybean simulations of many GGCMs are capable of reproducing larger parts of observed temporal variability (time series correlation coefficients (r) of up to 0.888 for maize, 0.673 for wheat and 0.643 for soybean at the global scale) but rice yield variability cannot be well reproduced by most models. Yield variability can be well reproduced for most major producing countries by many GGCMs and for all countries by at least some. A comparison with gridded yield data and a statistical analysis of the effects of weather variability on yield variability shows that the ensemble of GGCMs can explain more of the yield variability than an ensemble of regression models for maize and soybean, but not for wheat and rice. We identify future research needs in global gridded crop modeling and for all individual crop modeling groups. In the absence of a purely observation-based benchmark for model evaluation, we propose that the best performing crop model per crop and region establishes the benchmark for all others, and modelers are encouraged to investigate how crop model performance can be increased. We make our evaluation system accessible to all crop modelers so that other modeling groups can also test their model performance against the reference data and the GGCMI benchmark.

wheat↗

Benchmarking carbon fluxes of the ISIMIP2a biome models

The purpose of this study is to evaluate the eight ISIMIP2a biome models against independent estimates of long-term net carbon fluxes (i.e. Net Biome Productivity, NBP) over terrestrial ecosystems for the recent four decades (1971–2010). Here, we evaluate modeled global NBP against 1) the updated global residual land sink (RLS) plus land use emissions (E LUC ) from the Global Carbon Project (GCP), presented as R + L in this study by Le Quéré et al (2015), and 2) the land CO 2 fluxes from two atmospheric inversion systems: Jena CarboScope s81_v3.8 and CAMS v15r2, referred to as F Jena and F CAMS respectively. The model ensemble-mean NBP (that includes seven models with land-use change) is higher than but within the uncertainty of R + L, while the simulated positive NBP trend over the last 30 yr is lower than that from R + L and from the two inversion systems. ISIMIP2a biome models well capture the interannual variation of global net terrestrial ecosystem carbon fluxes. Tropical NBP represents 31 ± 17% of global total NBP during the past decades, and the year-to-year variation of tropical NBP contributes most of the interannual variation of global NBP. According to the models, increasing Net Primary Productivity (NPP) was the main cause for the generally increasing NBP. Significant global NBP anomalies from the long-term mean between the two phases of El Niño Southern Oscillation (ENSO) events are simulated by all models (p < 0.05), which is consistent with the R + L estimate (p = 0.06), also mainly attributed to NPP anomalies, rather than to changes in heterotrophic respiration (Rh). The global NPP and NBP anomalies during ENSO events are dominated by their anomalies in tropical regions impacted by tropical climate variability. Multiple regressions between R + L, F Jena and F CAMS interannual variations and tropical climate variations reveal a significant negative response of global net terrestrial ecosystem carbon fluxes to tropical mean annual temperature variation, and a non-significant response to tropical annual precipitation variation. According to the models, tropical precipitation is a more important driver, suggesting that some models do not capture the roles of precipitation and temperature changes adequately.

54 ENVIRONMENTAL SCIENCES↗