The Marine Biogeochemistry Library (MARBL) is a prognostic ocean biogeochemistry model that simulates marine ecosystem dynamics and the coupled cycles of carbon, nitrogen, phosphorus, iron, silicon, and oxygen. MARBL is a component of the Community Earth System Model (CESM); it supports flexible ecosystem configuration of multiple phytoplankton and zooplankton functional types; it is also portable, designed to interface with multiple ocean circulation models. Here, we present scientific documentation of MARBL, describe its configuration in CESM2 experiments included in the Coupled Model Intercomparison Project version 6 (CMIP6), and evaluate its performance against a number of observational data sets. The model simulates present‐day air‐sea CO₂ flux and many aspects of the carbon cycle in good agreement with observations. However, the simulated integrated uptake of anthropogenic CO₂ is weak, which we link to poor thermocline ventilation, a feature evident in simulated chlorofluorocarbon distributions. This also contributes to larger‐than‐observed oxygen minimum zones. Moreover, radiocarbon distributions show that the simulated circulation in the deep North Pacific is extremely sluggish, yielding extensive oxygen depletion and nutrient trapping at depth. Surface macronutrient biases are generally positive at low latitudes and negative at high latitudes. CESM2 simulates globally integrated net primary production (NPP) of 48 Pg C yr⁻¹ and particulate export flux at 100 m of 7.1 Pg C yr⁻¹. The impacts of climate change include an increase in globally integrated NPP, but substantial declines in the North Atlantic. Particulate export is projected to decline globally, attributable to decreasing export efficiency associated with changes in phytoplankton community composition.

A new set of stratospheric aerosol geoengineering (SAG) model experiments has been performed with Community Earth System Model version 2 (CESM2) with the Whole Atmosphere Community Climate Model (WACCM6) that are based on the Coupled Model Intercomparison Project phase 6 (CMIP6) overshoot scenario (SSP5-34-OS) as a baseline scenario to limit global warming to 1.5 or 2.0 ^∘C above 1850–1900 conditions. The overshoot scenario allows us to applying a peak-shaving scenario that reduces the needed duration and amount of SAG application compared to a high forcing scenario. In addition, a feedback algorithm identifies the needed amount of sulfur dioxide injections in the stratosphere at four pre-defined latitudes, 30^∘ N, 15^∘ N, 15^∘ S, and 30^∘ S, to reach three surface temperature targets: global mean temperature, and interhemispheric and pole-to-Equator temperature gradients. These targets further help to reduce side effects, including overcooling in the tropics, warming of high latitudes, and large shifts in precipitation patterns. These experiments are therefore relevant for investigating the impacts on society and ecosystems. Comparisons to SAG simulations based on a high emission pathway baseline scenario (SSP5-85) are also performed to investigate the dependency of impacts using different injection amounts to offset surface warming by SAG. We find that changes from present-day conditions around 2020 in some variables depend strongly on the defined temperature target (1.5 ^∘C vs. 2.0 ^∘C). These include surface air temperature and related impacts, the Atlantic Meridional Overturning Circulation, which impacts ocean net primary productivity, and changes in ice sheet surface mass balance, which impacts sea level rise. Others, including global precipitation changes and the recovery of the Antarctic ozone hole, depend strongly on the amount of SAG application. Furthermore, land net primary productivity as well as ocean acidification depend mostly on the global atmospheric CO₂ concentration and therefore the baseline scenario. Multi-model comparisons of experiments that include strong mitigation and carbon dioxide removal with some SAG application are proposed to assess the robustness of impacts on societies and ecosystems.

Anthropogenic land cover change (ALCC) is the most important transformation of the Earth system that occurred in the preindustrial Holocene, with implications for carbon, water and sediment cycles, biodiversity and the provision of ecosystem services and regional and global climate. For example, anthropogenic deforestation in preindustrial Eurasia may have led to feedbacks to the climate system: both biogeophysical, regionally amplifying winter cold and summer warm temperatures, and biogeochemical, stabilizing atmospheric CO ${}_2$ concentrations and thus influencing global climate. Quantification of these effects is difficult, however, because scenarios of anthropogenic land cover change over the Holocene vary widely, with increasing disagreement back in time. Because land cover change had such widespread ramifications for the Earth system, it is essential to assess current ALCC scenarios in light of observations and provide guidance on which models are most realistic. Here, we perform a systematic evaluation of two widely-used ALCC scenarios (KK10 and HYDE3.1) in northern and part of central Europe using an independent, pollen-based reconstruction of Holocene land cover (REVEALS). Considering that ALCC in Europe primarily resulted in deforestation, we compare modeled land use with the cover of non-forest vegetation inferred from the pollen data. Though neither land cover change scenario matches the pollen-based reconstructions precisely, KK10 correlates well with REVEALS at the country scale, while HYDE systematically underestimates land use with increasing magnitude with time in the past. Discrepancies between modeled and reconstructed land use are caused by a number of factors, including assumptions of per-capita land use and socio-cultural factors that cannot be predicted on the basis of the characteristics of the physical environment, including dietary preferences, long-distance trade, the location of urban areas and social organization.

Anthropogenic CO₂ emissions are inundating the upper ocean, acidifying the water, and altering the habitat for marine phytoplankton. These changes are thought to be particularly influential for calcifying phytoplankton, namely, coccolithophores. Coccolithophores are widespread and account for a substantial portion of open ocean calcification; changes in their abundance, distribution, or level of calcification could have far‐reaching ecological and biogeochemical impacts. Here, we isolate the effects of increasing CO₂ on coccolithophores using an explicit coccolithophore phytoplankton functional type parameterization in the Community Earth System Model. Coccolithophore growth and calcification are sensitive to changing aqueous CO₂. While holding circulation constant, we demonstrate that increasing CO₂ concentrations cause coccolithophores in most areas to decrease calcium carbonate production relative to growth. However, several oceanic regions show large increases in calcification, such the North Atlantic, Western Pacific, and parts of the Southern Ocean, due to an alleviation of carbon limitation for coccolithophore growth. Global annual calcification is 6% higher under present‐day CO₂ levels relative to preindustrial CO₂ (1.5 compared to 1.4 Pg C/year). However, under 900 μatm CO₂, global annual calcification is 11% lower than under preindustrial CO₂ levels (1.2 Pg C/year). Large portions of the ocean show greatly decreased coccolithophore calcification relative to growth, resulting in significant regional carbon export and air‐sea CO₂ exchange feedbacks. Our study implies that coccolithophores become more abundant but less calcified as CO₂ increases with a tipping point in global calcification (changing from increasing to decreasing calcification relative to preindustrial) at approximately ∼600 μatm CO₂.

Ocean alkalinity is critical to the uptake of atmospheric carbon in surface waters and provides buffering capacity towards the associated acidification. However, unlike dissolved inorganic carbon (DIC), alkalinity is not directly impacted by anthropogenic carbon emissions. Within the context of projections of future ocean carbon uptake and potential ecosystem impacts, especially through Coupled Model Intercomparison Projects (CMIPs), the representation of alkalinity and the main driver of its distribution in the ocean interior, the calcium carbonate cycle, have often been overlooked. Here we track the changes from CMIP5 to CMIP6 with respect to the Earth system model (ESM) representation of alkalinity and the carbonate pump which depletes the surface ocean in alkalinity through biological production of calcium carbonate and releases it at depth through export and dissolution. We report an improvement in the representation of alkalinity in CMIP6 ESMs relative to those in CMIP5, with CMIP6 ESMs simulating lower surface alkalinity concentrations, an increased meridional surface gradient and an enhanced global vertical gradient. This improvement can be explained in part by an increase in calcium carbonate (CaCO₃) production for some ESMs, which redistributes alkalinity at the surface and strengthens its vertical gradient in the water column. We were able to constrain a particulate inorganic carbon (PIC) export estimate of 44–55 Tmol yr⁻¹ at 100 m for the ESMs to match the observed vertical gradient of alkalinity. Reviewing the representation of the CaCO₃ cycle across CMIP5/6, we find a substantial range of parameterizations. While all biogeochemical models currently represent pelagic calcification, they do so implicitly, and they do not represent benthic calcification. In addition, most models simulate marine calcite but not aragonite. In CMIP6, certain model groups have increased the complexity of simulated CaCO₃ production, sinking, dissolution and sedimentation. However, this is insufficient to explain the overall improvement in the alkalinity representation, which is therefore likely a result of marine biogeochemistry model tuning or ad hoc parameterizations. Although modellers aim to balance the global alkalinity budget in ESMs in order to limit drift in ocean carbon uptake under pre-industrial conditions, varying assumptions related to the closure of the budget and/or the alkalinity initialization procedure have the potential to influence projections of future carbon uptake. For instance, in many models, carbonate production, dissolution and burial are independent of the seawater saturation state, and when considered, the range of sensitivities is substantial. As such, the future impact of ocean acidification on the carbonate pump, and in turn ocean carbon uptake, is potentially underestimated in current ESMs and is insufficiently constrained.

As anthropogenic emissions and warming continue to alter Earth's environment, it is essential to highlight future impacts that can be avoided through mitigation. Here we use two ensembles of the Community Earth System Model (CESM) run under the business‐as‐usual scenario, RCP 8.5, and the mitigation scenario, RCP 4.5, to identify avoidable impacts of anthropogenic warming on marine net primary production (NPP). We emphasize the use of ensembles so as to distinguish long‐term, anthropogenic trends in marine productivity from internal variability. Twentieth century globally integrated marine NPP is 55.7 ± 1 Pg C, with much of the variability attributable to certain regions (e.g., the equatorial Pacific). CESM projections indicate that global marine NPP will drop by ∼4% by 2080 if we follow RCP 8.5, but only by 2% under RCP 4.5. The response to warming on a global scale includes compensating regional effects; NPP increases in polar and eastern equatorial Pacific waters but decreases in the Atlantic, western Pacific, and Indian Oceans. The two main phytoplankton groups simulated in CESM show distinct responses: diatoms decrease their NPP, while small phytoplankton NPP increases over the mid‐21st century. Trends in NPP from mid‐21st century to 2080 are significantly different between the two emission scenarios mainly in the Atlantic Ocean basin and therefore impacts here are “avoidable” if we follow RCP 4.5, rather than RCP 8.5. In contrast, changes in NPP on a global scale and in most areas of the Pacific and Indian basins and the Southern Ocean are not distinguishable between forcing scenarios.