Search

Soil carbon (C) responses to environmental change represent a major source of uncertainty in the global C cycle. Feedbacks between soil C stocks and climate drivers could impact atmospheric CO2 levels, further altering the climate. Here, we assessed the reliability of Earth system model (ESM) predictions of soil C change using the Coupled Model Intercomparison Project phases 5 and 6 (CMIP5 and CMIP6). ESMs predicted global soil C gains under the high emission scenario, with soils taking up 43.9 Pg (95% CI: 9.2–78.5 Pg) C on average during the 21st century. The variation in global soil C change declined significantly from CMIP5 (with average of 48.4 Pg [95% CI: 2.0–94.9 Pg] C) to CMIP6 models (with average of 39.3 Pg [95% CI: 23.9–54.7 Pg] C). For some models, a small C increase in all biomes contributed to this convergence. For other models, offsetting responses between cold and warm biomes contributed to convergence. Although soil C predictions appeared to converge in CMIP6, the dominant processes driving soil C change at global or biome scales differed among models and in many cases between earlier and later versions of the same model. Random Forest models, for soil carbon dynamics, accounted for more than 63% variation of the global soil C change predicted by CMIP5 ESMs, but only 36% for CMIP6 models. Although most CMIP6 models apparently agree on increased soil C storage during the 21st century, this consensus obscures substantial model disagreement on the mechanisms underlying soil C response, calling into question the reliability of model predictions.

Abstract: In November 2021, the Artificial Intelligence for Earth System Predictability (AI4ESP) workshop was held, which involved hundreds of researchers from dozens of institutions. There were 17 sessions held at the workshop, including one on ecohydrology. The ecohydrology session included various breakout rooms that addressed specific topics, including 1) soils and belowground areas; 2) watersheds; 3) hydrology; 4) ecophysiology and plant hydraulics; 5) ecology; 6) extremes, disturbance and fire, and land-use and land-cover change; and 7) uncertainty quantification methods and techniques. In this paper, we investigate and report on the potential application of artificial intelligence and machine learning in ecohydrology, highlight outcomes of the ecohydrology session at the AI4ESP workshop, and provide visionary perspectives for future research in this area.

Global estimates of the land carbon sink are often based on simulations by terrestrial biosphere models (TBMs). The use of a large number of models that differ in their underlying hypotheses, structure and parameters is one way to assess the uncertainty in the historical land carbon sink. Here we show that the atmospheric forcing datasets used to drive these TBMs represent a significant source of uncertainty that is currently not systematically accounted for in land carbon cycle evaluations. We present results from three TBMs each forced with three different historical atmospheric forcing reconstructions over the period 1850-2015. We perform an analysis of variance to quantify the relative uncertainty in carbon fluxes arising from the models themselves, atmospheric forcing, and model-forcing interactions. We find that atmospheric forcing in this set of simulations plays a dominant role on uncertainties in global gross primary productivity (GPP) (75% of variability) and autotrophic respiration (90%), and a significant but reduced role on net primary productivity and heterotrophic respiration (30%). Atmospheric forcing is the dominant driver (52%) of variability for the net ecosystem exchange flux, defined as the difference between GPP and respiration (both autotrophic and heterotrophic respiration). In contrast, for wildfire-driven carbon emissions model uncertainties dominate and, as a result, model uncertainties dominate for net ecosystem productivity. At regional scales, the contribution of atmospheric forcing to uncertainty shows a very heterogeneous pattern and is smaller on average than at the global scale. We find that this difference in the relative importance of forcing uncertainty between global and regional scales is related to large differences in regional model flux estimates, which partially offset each other when integrated globally, while the flux differences driven by forcing are mainly consistent across the world and therefore add up to a larger fractional contribution to global uncertainty.

As the effects of anthropogenic climate change become more severe, several approaches for deliberate climate intervention to reduce or stabilize Earth's surface temperature have been proposed. Solar radiation modification (SRM) is one potential approach to partially counteract anthropogenic warming by reflecting a small proportion of the incoming solar radiation to increase Earth's albedo. While climate science research has focused on the predicted climate effects of SRM, almost no studies have investigated the impacts that SRM would have on ecological systems. The impacts and risks posed by SRM would vary by implementation scenario, anthropogenic climate effects, geographic region, and by ecosystem, community, population, and organism. Complex interactions among Earth's climate system and living systems would further affect SRM impacts and risks. We focus here on stratospheric aerosol intervention (SAI), a well-studied and relatively feasible SRM scheme that is likely to have a large impact on Earth's surface temperature. We outline current gaps in knowledge about both helpful and harmful predicted effects of SAI on ecological systems. Desired ecological outcomes might also inform development of future SAI implementation scenarios. In addition to filling these knowledge gaps, increased collaboration between ecologists and climate scientists would identify a common set of SAI research goals and improve the communication about potential SAI impacts and risks with the public. Without this collaboration, forecasts of SAI impacts will overlook potential effects on biodiversity and ecosystem services for humanity.

Large stocks of soil organic carbon (SOC) have accumulated in the Northern Hemisphere permafrost region, but their current amounts and future fate remain uncertain. By analyzing dataset combining >2700 soil profiles with environmental variables in a geospatial framework, we generated spatially explicit estimates of permafrost-region SOC stocks, quantified spatial heterogeneity, and identified key environmental predictors. We estimated that Pg C are stored in the top 3 m of permafrost region soils. The greatest uncertainties occurred in circumpolar toe-slope positions and in flat areas of the Tibetan region. We found that soil wetness index and elevation are the dominant topographic controllers and surface air temperature (circumpolar region) and precipitation (Tibetan region) are significant climatic controllers of SOC stocks. Our results provide first high-resolution geospatial assessment of permafrost region SOC stocks and their relationships with environmental factors, which are crucial for modeling the response of permafrost affected soils to changing climate.

In an era of rapid global change, our ability to understand and predict Earth's natural systems is lagging behind our ability to monitor and measure changes in the biosphere. Bottlenecks to informing models with observations have reduced our capacity to fully exploit the growing volume and variety of available data. Here, we take a critical look at the information infrastructure that connects ecosystem modeling and measurement efforts, and propose a roadmap to community cyberinfrastructure development that can reduce the divisions between empirical research and modeling and accelerate the pace of discovery. A new era of data-model integration requires investment in accessible, scalable, and transparent tools that integrate the expertise of the whole community, including both modelers and empiricists. This roadmap focuses on five key opportunities for community tools: the underlying foundations of community cyberinfrastructure; data ingest; calibration of models to data; model-data benchmarking; and data assimilation and ecological forecasting. This community-driven approach is a key to meeting the pressing needs of science and society in the 21st century.

Various approaches of differing mathematical complexities are being applied for spatial prediction of soil properties. Regression kriging is a widely used hybrid approach of spatial variation that combines correlation between soil properties and environmental factors with spatial autocorrelation between soil observations. In this study, we compared four machine learning approaches (gradient boosting machine, multinarrative adaptive regression spline, random forest, and support vector machine) with regression kriging to predict the spatial variation of surface (0-30 cm) soil organic carbon (SOC) stocks at 250-m spatial resolution across the northern circumpolar permafrost region. We combined 2,374 soil profile observations (calibration datasets) with georeferenced datasets of environmental factors (climate, topography, land cover, bedrock geology, and soil types) to predict the spatial variation of surface SOC stocks. We evaluated the prediction accuracy at randomly selected sites (validation datasets) across the study area. We found that different techniques inferred different numbers of environmental factors and their relative importance for prediction of SOC stocks. Regression kriging produced lower prediction errors in comparison to multinarrative adaptive regression spline and support vector machine, and comparable prediction accuracy to gradient boosting machine and random forest. However, the ensemble median prediction of SOC stocks obtained from all four machine learning techniques showed highest prediction accuracy. Although the use of different approaches in spatial prediction of soil properties will depend on the availability of soil and environmental datasets and computational resources, we conclude that the ensemble median prediction obtained from multiple machine learning approaches provides greater spatial details and produces the highest prediction accuracy. Thus an ensemble prediction approach can be a better choice than any single prediction technique for predicting the spatial variation of SOC stocks.

The Community Land Model (CLM) is the land component of the Community Earth System Model (CESM) and is used in several global and regional modeling systems. In this paper, we introduce model developments included in CLM version 5 (CLM5), which is the default land component for CESM2. We assess an ensemble of simulations, including prescribed and prognostic vegetation state, multiple forcing data sets, and CLM4, CLM4.5, and CLM5, against a range of metrics including from the International Land Model Benchmarking (ILAMBv2) package. CLM5 includes new and updated processes and parameterizations: (1) dynamic land units, (2) updated parameterizations and structure for hydrology and snow (spatially explicit soil depth, dry surface layer, revised groundwater scheme, revised canopy interception and canopy snow processes, updated fresh snow density, simple firn model, and Model for Scale Adaptive River Transport), (3) plant hydraulics and hydraulic redistribution, (4) revised nitrogen cycling (flexible leaf stoichiometry, leaf N optimization for photosynthesis, and carbon costs for plant nitrogen uptake), (5) global crop model with six crop types and time-evolving irrigated areas and fertilization rates, (6) updated urban building energy, (7) carbon isotopes, and (8) updated stomatal physiology. New optional features include demographically structured dynamic vegetation model (Functionally Assembled Terrestrial Ecosystem Simulator), ozone damage to plants, and fire trace gas emissions coupling to the atmosphere. Conclusive establishment of improvement or degradation of individual variables or metrics is challenged by forcing uncertainty, parametric uncertainty, and model structural complexity, but the multivariate metrics presented here suggest a general broad improvement from CLM4 to CLM5.

Over the past several decades, the land modeling community has recognized the importance of nutrient regulation on the global terrestrial carbon cycle. Implementations of nutrient limitation in land models are diverse, varying from applying simple empirical down-regulation of potential gross primary productivity under nutrient deficit conditions to more mechanistic treatments. In this study, we introduce a new approach to model multinutrient (nitrogen [N] and phosphorus [P]) limitations in the Energy Exascale Earth System Model (E3SM) Land Model version 1 (ELMv1-ECA). The development is grounded on (1) advances in representing multiple-consumer, multiple-nutrient competition; (2) a generic dynamic allocation scheme based on water, N, P, and light availability; (3) flexible plant CNP stoichiometry; (4) prognostic treatment of N and P constraints on several carbon cycle processes; and (5) global data sets of plant physiological traits. Through benchmarking the model against best knowledge of global plant and soil carbon pools and fluxes, we show that our implementation of nutrient constraints on the present-day carbon cycle is robust at the global scale. Compared with predecessor versions, ELMv1-ECA better predicts global-scale gross primary productivity, ecosystem respiration, leaf area index, vegetation biomass, soil carbon stocks, evapotranspiration, N2O emissions, and NO3- leaching. Factorial experiments indicate that representing the phosphorus cycle improves modeled carbon fluxes, while considering dynamic allocation improves modeled carbon stock density. We also highlight the value of using the International Land Model Benchmarking (ILAMB) package to evaluate and document performance during model development.

Earth system models are complex and represent a large number of processes, resulting in a persistent spread across climate projections for a given future scenario. Owing to different model performances against observations and the lack of independence among models, there is now evidence that giving equal weight to each available model projection is suboptimal. This Perspective discusses newly developed tools that facilitate a more rapid and comprehensive evaluation of model simulations with observations, process-based emergent constraints that are a promising way to focus evaluation on the observations most relevant to climate projections, and advanced methods for model weighting. These approaches are needed to distil the most credible information on regional climate changes, impacts, and risks for stakeholders and policy-makers.

Climate change is expected to increase the frequency of flooding events and, thus, the risks of flood-related mortality and infrastructure damage. Global-scale assessments of future flooding from Earth system models based only on precipitation changes neglect important processes that occur within the land surface, particularly plant physiological responses to rising CO2. Higher CO2 can reduce stomatal conductance and transpiration, which may lead to increased soil moisture and runoff in some regions, promoting flooding even without changes in precipitation. Here we assess the relative impacts of plant physiological and radiative greenhouse effects on changes in daily runoff intensity over tropical continents using the Community Earth System Model. We find that extreme percentile rates increase significantly more than mean runoff in response to higher CO2. Plant physiological effects have a small impact on precipitation intensity but are a dominant driver of runoff intensification, contributing to one half of the 99th and one third of the 99.9th percentile runoff intensity changes.

The increasing complexity of Earth system models has inspired efforts to quantitatively assess model fidelity through rigorous comparison with best available measurements and observational data products. Earth system models exhibit a high degree of spread in predictions of land biogeochemistry, biogeophysics, and hydrology, which are sensitive to forcing from other model components. Based on insights from prior land model evaluation studies and community workshops, the authors developed an open source model benchmarking software package that generates graphical diagnostics and scores model performance in support of the International Land Model Benchmarking (ILAMB) project. Employing a suite of in situ, remote sensing, and reanalysis data sets, the ILAMB package performs comprehensive model assessment across a wide range of land variables and generates a hierarchical set of web pages containing statistical analyses and figures designed to provide the user insights into strengths and weaknesses of multiple models or model versions. Described here is the benchmarking philosophy and mathematical methodology embodied in the most recent implementation of the ILAMB package. Comparison methods unique to a few specific data sets are presented, and guidelines for configuring an ILAMB analysis and interpreting resulting model performance scores are discussed. ILAMB is being adopted by modeling teams and centers during model development and for model intercomparison projects, and community engagement is sought for extending evaluation metrics and adding new observational data sets to the benchmarking framework.

Understanding how anthropogenic CO2 emissions will influence future precipitation is critical for sustainably managing ecosystems, particularly for drought-sensitive tropical forests. Although tropical precipitation change remains uncertain, nearly all models from the Coupled Model Intercomparison Project Phase 5 predict a strengthening zonal precipitation asymmetry by 2100, with relative increases over Asian and African tropical forests and decreases over South American forests. Here we show that the plant physiological response to increasing CO2 is a primary mechanism responsible for this pattern. Applying a simulation design in the Community Earth System Model in which CO2 increases are isolated over individual continents, we demonstrate that different circulation, moisture and stability changes arise over each continent due to declines in stomatal conductance and transpiration. The sum of local atmospheric responses over individual continents explains the pan-tropical precipitation asymmetry. Our analysis suggests that South American forests may be more vulnerable to rising CO2 than Asian or African forests.

Using the Community Earth System Model, we explore the role of human land use and land cover change (LULCC) in modifying the terrestrial carbon budget in simulations forced by Representative Concentration Pathway 8.5, extended to year 2300. Overall, conversion of land (e.g., from forest to croplands via deforestation) results in a model-estimated, cumulative carbon loss of 490 Pg C between 1850 and 2300, larger than the 230 Pg C loss of carbon caused by climate change over this same interval. The LULCC carbon loss is a combination of a direct loss at the time of conversion and an indirect loss from the reduction of potential terrestrial carbon sinks. Approximately 40% of the carbon loss associated with LULCC in the simulations arises from direct human modification of the land surface; the remaining 60% is an indirect consequence of the loss of potential natural carbon sinks. Because of the multicentury carbon cycle legacy of current land use decisions, a globally averaged amplification factor of 2.6 must be applied to 2015 land use carbon losses to adjust for indirect effects. This estimate is 30% higher when considering the carbon cycle evolution after 2100. Most of the terrestrial uptake of anthropogenic carbon in the model occurs from the influence of rising atmospheric CO2 on photosynthesis in trees, and thus, model-projected carbon feedbacks are especially sensitive to deforestation.

This report presents results from the DOE-sponsored workshop titled, Advancing X-Cutting Ideas for Computational Climate Science Workshop,'' known as AXICCS, held on September 12--13, 2016 in Rockville, MD. The workshop brought together experts in climate science, computational climate science, computer science, and mathematics to discuss interesting but unsolved science questions regarding climate modeling and simulation, promoted collaboration among the diverse scientists in attendance, and brainstormed about possible tools and capabilities that could be developed to help address them. Emerged from discussions at the workshop were several research opportunities that the group felt could advance climate science significantly. These include (1) process-resolving models to provide insight into important processes and features of interest and inform the development of advanced physical parameterizations, (2) a community effort to develop and provide integrated model credibility, (3) including, organizing, and managing increasingly connected model components that increase model fidelity yet complexity, and (4) treating Earth system models as one interconnected organism without numerical or data based boundaries that limit interactions. The group also identified several cross-cutting advances in mathematics, computer science, and computational science that would be needed to enable one or more of these big ideas. It is critical to address the need for organized, verified, and optimized software, which enables the models to grow and continue to provide solutions in which the community can have confidence. Effectively utilizing the newest computer hardware enables simulation efficiency and the ability to handle output from increasingly complex and detailed models. This will be accomplished through hierarchical multiscale algorithms in tandem with new strategies for data handling, analysis, and storage. These big ideas and cross-cutting technologies for enabling breakthrough climate simulation advancements also need the "glue" of outreach and learning across the scientific domains to be successful. The workshop identified several strategies to allow productive, continuous engagement across those who have a broad knowledge of the various angles of the problem. Specific ideas to foster education and tools to make material progress were discussed. Examples include follow-on cross-cutting meetings that enable unstructured discussions of the types this workshop fostered. A concerted effort to recruit undergraduate and graduate students from all relevant domains and provide them experience, training, and networking across their immediate expertise is needed. This will broaden and expand their exposure to the future needs and solutions, and provide a pipeline of scientists with a diversity of knowledge and know-how. Providing real-world experience with subject matter experts from multiple angles may also motivate the students to attack these problems and even come up with the missing solutions.

Rising atmospheric CO2 will make Earth warmer, and many studies have inferred that this warming will cause droughts to become more widespread and severe. However, rising atmospheric CO2 also modifies stomatal conductance and plant water use, processes that are often are overlooked in impact analysis. We find that plant physiological responses to CO2 reduce predictions of future drought stress, and that this reduction is captured by using plant-centric rather than atmosphere-centric metrics from Earth system models (ESMs). The atmosphere-centric Palmer Drought Severity Index predicts future increases in drought stress for more than 70% of global land area. This area drops to 37% with the use of precipitation minus evapotranspiration (P-E), a measure that represents the water flux available to downstream ecosystems and humans. The two metrics yield consistent estimates of increasing stress in regions where precipitation decreases are more robust (southern North America, northeastern South America, and southern Europe). The metrics produce diverging estimates elsewhere, with P-E predicting decreasing stress across temperate Asia and central Africa. The differing sensitivity of drought metrics to radiative and physiological aspects of increasing CO2 partly explains the divergent estimates of future drought reported in recent studies. Further, use of ESM output in offline models may double-count plant feedbacks on relative humidity and other surface variables, leading to overestimates of future stress. The use of drought metrics that account for the response of plant transpiration to changing CO2, including direct use of P-E and soil moisture from ESMs, is needed to reduce uncertainties in future assessment.

Coordinated experimental design and implementation has become a cornerstone of global climate modelling. Model Intercomparison Projects (MIPs) enable systematic and robust analysis of results across many models, by reducing the influence of ad hoc differences in model set-up or experimental boundary conditions. As it enters its 6th phase, the Coupled Model Intercomparison Project (CMIP6) has grown significantly in scope with the design and documentation of individual simulations delegated to individual climate science communities. The Coupled Climate-Carbon Cycle Model Intercomparison Project (C4MIP) takes responsibility for design, documentation, and analysis of carbon cycle feedbacks and interactions in climate simulations. These feedbacks are potentially large and play a leading-order contribution in determining the atmospheric composition in response to human emissions of CO2 and in the setting of emissions targets to stabilize climate or avoid dangerous climate change. For over a decade, C4MIP has coordinated coupled climate-carbon cycle simulations, and in this paper we describe the C4MIP simulations that will be formally part of CMIP6. While the climate-carbon cycle community has created this experimental design, the simulations also fit within the wider CMIP activity, conform to some common standards including documentation and diagnostic requests, and are designed to complement the CMIP core experiments known as the Diagnostic, Evaluation and Characterization of Klima (DECK). C4MIP has three key strands of scientific motivation and the requested simulations are designed to satisfy their needs: (1) pre-industrial and historical simulations (formally part of the common set of CMIP6 experiments) to enable model evaluation, (2) idealized coupled and partially coupled simulations with 1% per year increases in CO2 to enable diagnosis of feedback strength and its components, (3) future scenario simulations to project how the Earth system will respond to anthropogenic activity over the 21st century and beyond. This paper documents in detail these simulations, explains their rationale and planned analysis, and describes how to set up and run the simulations. Particular attention is paid to boundary conditions, input data, and requested output diagnostics. It is important that modelling groups participating in C4MIP adhere as closely as possible to this experimental design.

Organic macromolecules constitute a high percentage of remote sea spray. They enter the atmosphere through adsorption onto bubbles followed by bursting at the ocean surface, and go on to influence the chemistry of the fine mode aerosol. We present a global estimate of mixed-layer macromolecular distributions, driven by offline marine systems model output. The approach permits estimation of oceanic concentrations and bubble film surface coverages for several classes of organic compound. Mixed layer levels are computed from the output of a global ocean ecodynamics model by relating the macromolecules to standard biogeochemical tracers. Steady state is assumed for labile forms, and for longer-lived components we rely on ratios to existing transported variables. Adsorption is then represented through conventional Langmuir isotherms, with equilibria deduced from laboratory analogs. Open water concentrations locally exceed one micromolar carbon for the total of proteins, polysaccharides and refractory heteropolycondensates. The shorter-lived lipids remain confined to regions of strong biological activity. Results are evaluated against available measurements for all compound types, and agreement is generally well within an order of magnitude. Global distributions are further estimated for both fractional coverage of bubble films at the air–water interface and the two-dimensional concentration excess. Overall, we show that macromolecular mapping provides a novel tool for the comprehension of oceanic surfactant patterns. These results may prove useful in planning field experiments and assessing the potential response of surface chemical behaviors to global change.

Version 1 of the Community Earth SystemModel, in the configurationwhere its full carbon cycle is enabled, is introduced and documented. In this configuration, the terrestrial biogeochemical model, which includes carbon-nitrogen dynamics and is present in earlier model versions, is coupled to an ocean biogeochemical model and atmospheric CO2 tracers. The authors provide a description of the model, detail how preindustrial-control and twentieth-century experiments were initialized and forced, and examine the behavior of the carbon cycle in those experiments. They examine how sea- and land-to-air CO2 fluxes contribute to the increase of atmospheric CO2 in the twentieth century, analyze how atmospheric CO2 and its surface fluxes vary on interannual time scales, including how they respond to ENSO, and describe the seasonal cycle of atmospheric CO2 and its surface fluxes. While the model broadly reproduces observed aspects of the carbon cycle, there are several notable biases, including having too large of an increase in atmospheric CO2 over the twentieth century and too small of a seasonal cycle of atmospheric CO2 in the Northern Hemisphere. The biases are related to a weak response of the carbon cycle to climatic variations on interannual and seasonal time scales and to twentieth-century anthropogenic forcings, including rising CO2, land-use change, and atmospheric deposition of nitrogen.

Changes in atmospheric CO2 variability during the twenty-first century may provide insight about ecosystem responses to climate change and have implications for the design of carbon monitoring programs. This paper describes changes in the three-dimensional structure of atmospheric CO2 for several representative concentration pathways (RCPs 4.5 and 8.5) using the Community Earth System Model–Biogeochemistry (CESM1-BGC). CO2 simulated for the historical period was first compared to surface, aircraft, and column observations. In a second step, the evolution of spatial and temporal gradients during the twenty-first century was examined. The mean annual cycle in atmospheric CO2 was underestimated for the historical period throughout the Northern Hemisphere, suggesting that the growing season net flux in the Community Land Model (the land component of CESM) was too weak. Consistent with weak summer drawdown in Northern Hemisphere high latitudes, simulated CO2 showed correspondingly weak north–south and vertical gradients during the summer. In the simulations of the twenty-first century, CESM predicted increases in the mean annual cycle of atmospheric CO2 and larger horizontal gradients. Not only did the mean north–south gradient increase due to fossil fuel emissions, but east–west contrasts in CO2 also strengthened because of changing patterns in fossil fuel emissions and terrestrial carbon exchange. In the RCP8.5 simulation, where CO2 increased to 1150 ppm by 2100, the CESM predicted increases in interannual variability in the Northern Hemisphere midlatitudes of up to 60% relative to present variability for time series filtered with a 2–10-yr bandpass. Such an increase in variability may impact detection of changing surface fluxes from atmospheric observations.

With representation of the global carbon cycle becoming increasingly complex in climate models, it is important to develop ways to quantitatively evaluate model performance against in situ and remote sensing observations. Here we present a systematic framework, the Carbon-LAnd Model Intercomparison Project (C-LAMP), for assessing terrestrial biogeochemistry models coupled to climate models using observations that span a wide range of temporal and spatial scales. As an example of the value of such comparisons, we used this framework to evaluate two biogeochemistry models that are integrated within the Community Climate System Model (CCSM) -Carnegie-Ames-Stanford Approach′ (CASA′) and carbon-nitrogen (CN). Both models underestimated the magnitude of net carbon uptake during the growing season in temperate and boreal forest ecosystems, based on comparison with atmospheric CO 2 measurements and eddy covariance measurements of net ecosystem exchange. Comparison with MODerate Resolution Imaging Spectroradiometer (MODIS) measurements show that this low bias in model fluxes was caused, at least in part, by 1-3 month delays in the timing of maximum leaf area. In the tropics, the models overestimated carbon storage in woody biomass based on comparison with datasets from the Amazon. Reducing this model bias will probably weaken the sensitivity of terrestrial carbon fluxes to both atmospheric CO2 and climate. Global carbon sinks during the 1990s differed by a factor of two (2.4PgCyr-1 for CASA′ vs. 1.2PgCyr-1 for CN), with fluxes from both models compatible with the atmospheric budget given uncertainties in other terms. The models captured some of the timing of interannual global terrestrial carbon exchange during 1988-2004 based on comparison with atmospheric inversion results from TRANSCOM (r =0.66 for CASA′ and r =0.73 for CN). Adding (CASA′) or improving (CN) the representation of deforestation fires may further increase agreement with the atmospheric record. Information from C-LAMP has enhanced model performance within CCSM and serves as a benchmark for future development. We propose that an open source, community-wide platform for model-data intercomparison is needed to speed model development and to strengthen ties between modeling and measurement communities. Important next steps include the design and analysis of land use change simulations (in both uncoupled and coupled modes), and the entrainment of additional ecological and earth system observations. Model results from C-LAMP are publicly available on the Earth System Grid. © 2009 Blackwell Publishing Ltd.