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.

The TransCom 3 experiment was begun to explore the estimation of carbon sources and sinks via the inversion of simulated tracer transport. We build upon previous TransCom work by presenting the seasonal inverse results which provide estimates of carbon flux for 11 land and 11 ocean regions using 12 atmospheric transport models. The monthly fluxes represent the mean seasonal cycle for the 1992 to 1996 time period. The spread among the model results is larger than the average of their estimated flux uncertainty in the northern extratropics and vice versa in the tropical regions. In the northern land regions, the model spread is largest during the growing season. Compared to a seasonally balanced biosphere prior flux generated by the CASA model, we find significant changes to the carbon exchange in the European region with greater growing season net uptake which persists into the fall months. Both Boreal North America and Boreal Asia show lessened net uptake at the onset of the growing season with Boreal Asia also exhibiting greater peak growing season net uptake. Temperate Asia shows a dramatic springward shift in the peak timing of growing season net uptake relative to the neutral CASA flux while Temperate North America exhibits a broad flattening of the seasonal cycle. In most of the ocean regions, the inverse fluxes exhibit much greater seasonality than that implied by the ΔpCO₂ derived fluxes though this may be due, in part, to misallocation of adjacent land flux. In the Southern Ocean, the austral spring and fall exhibits much less carbon uptake than implied by ΔpCO₂ derived fluxes. Sensitivity testing indicates that the inverse estimates are not overly influenced by the prior flux choices. Considerable agreement exists between the model mean, annual mean results of this study and that of the previously published TransCom annual mean inversion. The differences that do exist are in poorly constrained regions and tend to exhibit compensatory fluxes in order to match the global mass constraint. The differences between the estimated fluxes and the prior model over the northern land regions could be due to the prior model respiration response to temperature. Significant phase differences, such as that in the Temperate Asia region, may be due to the limited observations for that region. Finally, differences in the boreal land regions between the prior model and the estimated fluxes may be a reflection of the timing of spring thaw and an imbalance in respiration versus photosynthesis

 Inverse estimation of carbon dioxide (CO₂) sources and sinks uses atmospheric CO₂ observations, mostly made near the Earth's surface. However, transport models used in such studies lack perfect representation of atmospheric dynamics and thus often fail to produce unbiased forward simulations. The error is generally larger for observations over the land than those over the remote/marine locations. The range of this error is estimated by using multiple transport models (16 are used here). We have estimated the remaining differences in CO₂ fluxes due to the use of ocean-only versus all-sites (i.e., over ocean and land) observations of CO₂ in a time-independent inverse modeling framework. The fluxes estimated using the ocean-only networks are more robust compared to those obtained using all-sites networks. This makes the global, hemispheric, and regional flux determination less dependent on the selection of transport model and observation network.

Information about regional carbon sources and sinks can be derived from variations in observed atmospheric CO2 concentrations via inverse modelling with atmospheric tracer transport models. A consensus has not yet been reached regarding the size and distribution of regional carbon fluxes obtained using this approach, partly owing to the use of several different atmospheric transport models. Here we report estimates of surface-atmosphere CO2 fluxes from an intercomparison of atmospheric CO2 inversion models (the TransCom 3 project), which includes 16 transport models and model variants. We find an uptake of CO2 in the southern extratropical ocean less than that estimated from ocean measurements, a result that is not sensitive to transport models or methodological approaches. We also find a northern land carbon sink that is distributed relatively evenly among the continents of the Northern Hemisphere, but these results show some sensitivity to transport differences among models, especially in how they respond to seasonal terrestrial exchange of CO2. Overall, carbon fluxes integrated over latitudinal zones are strongly constrained by observations in the middle to high latitudes. Further significant constraints to our understanding of regional carbon fluxes will therefore require improvements in transport models and expansion of the CO2 observation network within the tropics.

The relative skill of 21 regional and global biogeochemical models was assessed in terms of how well the models reproduced observed net primary productivity (NPP) and environmental variables such as nitrate concentration (NO3), mixed layer depth (MLD), euphotic layer depth (Zeu), and sea ice concentration, by comparing results against a newly updated, quality-controlled in situ NPP database for the Arctic Ocean (1959-2011). The models broadly captured the spatial features of integrated NPP (iNPP) on a pan-Arctic scale. Most models underestimated iNPP by varying degrees in spite of overestimating surface NO3, MLD, and Zeu throughout the regions. Among the models, iNPP exhibited little difference over sea ice condition (ice-free versus ice-influenced) and bottom depth (shelf versus deep ocean). The models performed relatively well for the most recent decade and toward the end of Arctic summer. In the Barents and Greenland Seas, regional model skill of surface NO3 was best associated with how well MLD was reproduced. Regionally, iNPP was relatively well simulated in the Beaufort Sea and the central Arctic Basin, where in situ NPP is low and nutrients are mostly depleted. Models performed less well at simulating iNPP in the Greenland and Chukchi Seas, despite the higher model skill in MLD and sea ice concentration, respectively. iNPP model skill was constrained by different factors in different Arctic Ocean regions. Our study suggests that better parameterization of biological and ecological microbial rates (phytoplankton growth and zooplankton grazing) are needed for improved Arctic Ocean biogeochemical modeling.

Projections of climate change impacts on marine ecosystems have revealed long-term declines in global marine animal biomass and unevenly distributed impacts on fisheries. Here we apply an enhanced suite of global marine ecosystem models from the Fisheries and Marine Ecosystem Model Intercomparison Project (Fish-MIP), forced by new-generation Earth system model outputs from Phase 6 of the Coupled Model Intercomparison Project (CMIP6), to provide insights into how projected climate change will affect future ocean ecosystems. Compared with the previous generation CMIP5-forced Fish-MIP ensemble, the new ensemble ecosystem simulations show a greater decline in mean global ocean animal biomass under both strong-mitigation and high-emissions scenarios due to elevated warming, despite greater uncertainty in net primary production in the high-emissions scenario. Regional shifts in the direction of biomass changes highlight the continued and urgent need to reduce uncertainty in the projected responses of marine ecosystems to climate change to help support adaptation planning.