Ocean carbon uptake and storage simulated by the Community Earth System Model, version 1–Biogeochemistry [CESM1(BGC)], is described and compared to observations. Fully coupled and ocean-ice configurations are examined; both capture many aspects of the spatial structure and seasonality of surface carbon fields. Nearly ubiquitous negative biases in surface alkalinity result from the prescribed carbonate dissolution profile. The modeled sea–air CO₂ fluxes match observationally based estimates over much of the ocean; significant deviations appear in the Southern Ocean. Surface ocean pCO₂ is biased high in the subantarctic and low in the sea ice zone. Formation of the water masses dominating anthropogenic CO₂ (C_ant) uptake in the Southern Hemisphere is weak in the model, leading to significant negative biases in C_ant and chlorofluorocarbon (CFC) storage at intermediate depths. Column inventories of C_ant appear too high, by contrast, in the North Atlantic. In spite of the positive bias, this marks an improvement over prior versions of the model, which underestimated North Atlantic uptake. The change in behavior is attributable to a new parameterization of density-driven overflows. CESM1(BGC) provides a relatively robust representation of the ocean–carbon cycle response to climate variability. Statistical metrics of modeled interannual variability in sea–air CO₂ fluxes compare reasonably well to observationally based estimates. The carbon cycle response to key modes of climate variability is basically similar in the coupled and forced ocean-ice models; however, the two differ in regional detail and in the strength of teleconnections.