Description: The ocean is a major carbon sink and takes up 25%–30% of the anthropogenically emitted CO2. A state-of-the-art method to quantify this sink are global ocean biogeochemistry models (GOBMs), but their simulated CO2 uptake differs between models and is systematically lower than estimates based on statistical methods using surface ocean pCO2 and interior ocean measurements. Here, we provide an in-depth evaluation of ocean carbon sink estimates from 1980 to 2018 from a GOBM ensemble. As sources of inter-model differences and ensemble-mean biases our study identifies (a) the model setup, such as the length of the spin-up, the starting date of the simulation, and carbon fluxes from rivers and into sediments, (b) the simulated ocean circulation, such as Atlantic Meridional Overturning Circulation and Southern Ocean mode and intermediate water formation, and (c) the simulated oceanic buffer capacity. Our analysis suggests that a late starting date and biases in the ocean circulation cause a too low anthropogenic CO2 uptake across the GOBM ensemble. Surface ocean biogeochemistry biases might also cause simulated anthropogenic fluxes to be too low, but the current setup prevents a robust assessment. For simulations of the ocean carbon sink, we recommend in the short-term to (a) start simulations at a common date before the industrialization and the associated atmospheric CO2 increase, (b) conduct a sufficiently long spin-up such that the GOBMs reach steady-state, and (c) provide key metrics for circulation, biogeochemistry, and the land-ocean interface. In the long-term, we recommend improving the representation of these metrics in the GOBMs.

Description: The products from the Stable Water Isotope Intercomparison Group, Phase 2, are currently used for numerous studies, allowing water isotope model-data comparisons with various isotope-enabled atmospheric general circulation model (AGCMs) outputs. However, the simulations under this framework were performed using different parameterizations and forcings. Therefore, a uniform experimental design with state-of-the-art AGCMs is required to interpret isotope observations rigorously. Here, we evaluate the outputs from three isotope-enabled numerical models nudged by three different reanalysis products and investigate the ability of the isotope-enabled AGCMs to reproduce the spatial and temporal patterns of water isotopic composition observed at the surface and in the atmospheric airborne water. Through correlation analyses at various spatial and temporal scales, we found that the model's performance depends on the model or reanalysis we use, the observations we compare, and the vertical levels we select. Moreover, we employed the stable isotope mass balance method to conduct decomposition analyses on the ratio of isotopic changes in the atmosphere. Our goal was to elucidate the spread in simulated atmospheric column δ18O, which is influenced by factors such as evaporation, precipitation, and horizontal moisture flux. Satisfying the law of conservation of water isotopes, this budget method is expected to explain various fractionation phenomena in atmospheric meteorological and climatic events. It also aims to highlight the spreads in modeled isotope results among different experiments using multiple models and reanalyses, which are primarily dominated by uncertainties in moisture flux and precipitation, respectively.

Description: Significant progress in permafrost carbon science made over the past decades include the identification of vast permafrost carbon stocks, the development of new pan‐Arctic permafrost maps, an increase in terrestrial measurement sites for CO〈jats:sub〉2〈/jats:sub〉 and methane fluxes, and important factors affecting carbon cycling, including vegetation changes, periods of soil freezing and thawing, wildfire, and other disturbance events. Process‐based modeling studies now include key elements of permafrost carbon cycling and advances in statistical modeling and inverse modeling enhance understanding of permafrost region C budgets. By combining existing data syntheses and model outputs, the permafrost region is likely a wetland methane source and small terrestrial ecosystem CO〈jats:sub〉2〈/jats:sub〉 sink with lower net CO〈jats:sub〉2〈/jats:sub〉 uptake toward higher latitudes, excluding wildfire emissions. For 2002–2014, the strongest CO〈jats:sub〉2〈/jats:sub〉 sink was located in western Canada (median: −52 g C m〈jats:sup〉−2〈/jats:sup〉 y〈jats:sup〉−1〈/jats:sup〉) and smallest sinks in Alaska, Canadian tundra, and Siberian tundra (medians: −5 to −9 g C m〈jats:sup〉−2〈/jats:sup〉 y〈jats:sup〉−1〈/jats:sup〉). Eurasian regions had the largest median wetland methane fluxes (16–18 g CH〈jats:sub〉4〈/jats:sub〉 m〈jats:sup〉−2〈/jats:sup〉 y〈jats:sup〉−1〈/jats:sup〉). Quantifying the regional scale carbon balance remains challenging because of high spatial and temporal variability and relatively low density of observations. More accurate permafrost region carbon fluxes require: (a) the development of better maps characterizing wetlands and dynamics of vegetation and disturbances, including abrupt permafrost thaw; (b) the establishment of new year‐round CO〈jats:sub〉2〈/jats:sub〉 and methane flux sites in underrepresented areas; and (c) improved models that better represent important permafrost carbon cycle dynamics, including non‐growing season emissions and disturbance effects.〈/jats:p〉

Description: Current global warming results in rising sea-water temperatures, and the loss of sea ice in arctic and subarctic oceans impacts the community composition of primary producers with cascading effects on the food web and potentially on carbon export rates. This study analyzes metagenomic shotgun and diatom rbcL amplicon-sequencing data from sedimentary ancient DNA (sedaDNA) of the subarctic western Bering Sea that records phyto- and zooplankton community changes over the last glacial–interglacial cycle, including the last interglacial period (Eemian). Our data show that interglacial and glacial plankton communities differ, with distinct Eemian and Holocene plankton communities. The generally warm Holocene period is dominated by pico-sized cyanobacteria and bacteria-feeding heterotrophic protists, while the Eemian period is dominated by eukaryotic pico-sized chlorophytes and Triparmaceae. In contrast, the glacial period is characterized by micro-sized phototrophic protists, including sea-ice associated diatoms in the family Bacillariaceae and co-occurring diatom-feeding crustaceous zooplankton. Our deep-time record of plankton community changes reveals a long-term decrease in phytoplankton cell size coeval with increasing temperatures, and resembling community changes in the currently warming Bering Sea. The phytoplankton community in the warmer-than-present Eemian period is distinct from modern communities and limits the use of the Eemian as an analog for future climate scenarios. However, under enhanced future warming, the expected shift towards the dominance of small-sized phytoplankton and heterotrophic protists might result in an increased productivity, whereas the community’s potential of carbon export will be decreased, thereby weakening the subarctic Bering Sea’s function as an effective carbon sink.