Abstract. We developed a complex Earth system model by coupling terrestrial and oceanic carbon cycle components into the Bergen Climate Model. For this study, we have generated two model simulations (one with climate change inclusions and the other without) to study the large scale climate and carbon cycle variability as well as its feedback for the period 1850–2100. The simulations are performed based on historical and future IPCC CO₂ emission scenarios. Globally, a pronounced positive climate-carbon cycle feedback is simulated by the terrestrial carbon cycle model, but smaller signals are shown by the oceanic counterpart. Over land, the regional climate-carbon cycle feedback is highlighted by increased soil respiration, which exceeds the enhanced production due to the atmospheric CO₂ fertilization effect, in the equatorial and northern hemisphere mid-latitude regions. For the ocean, our analysis indicates that there are substantial temporal and spatial variations in climate impact on the air-sea CO₂ fluxes. This implies feedback mechanisms act inhomogeneously in different ocean regions. In the North Atlantic subpolar gyre, the simulated future cooling of SST improves the CO₂ gas solubility in seawater and, hence, reduces the strength of positive climate carbon cycle feedback in this region. In most ocean regions, the changes in the Revelle factor is dominated by changes in surface pCO₂, and not by the warming of SST. Therefore, the solubility-associated positive feedback is more prominent than the buffer capacity feedback. In our climate change simulation, the retreat of Southern Ocean sea ice due to melting allows an additional ~20 Pg C uptake as compared to the simulation without climate change.