Search

In this paper we describe the implementation of the carbon isotopes 13 C and 14 C (radiocarbon) into the marine biogeochemistry model REcoM3. The implementation is tested in long-term equilibrium simulations where REcoM3 is coupled with the ocean general circulation model FESOM2.1, applying a low-resolution configuration and idealized climate forcing. Focusing on the carbon-isotopic composition of dissolved inorganic carbon (δ13 C DIC and Δ14 C DIC), our model results are largely consistent with reconstructions for the pre-anthropogenic period. Our simulations also exhibit discrepancies, e.g. in upwelling regions and the interior of the North Pacific. Some of these differences are due to the limitations of our ocean circulation model setup, which results in a rather shallow meridional overturning circulation. We additionally study the accuracy of two simplified modelling approaches for dissolved inorganic 14 C, which are faster (15 % and about a factor of five, respectively) than the complete consideration of the marine radiocarbon cycle. The accuracy of both simplified approaches is better than 5 %, which should be sufficient for most studies of Δ14 C DIC. [ABSTRACT FROM AUTHOR]

This study describes the coupling of the process-based Model of Early Diagenesis in the Upper Sediment (MEDUSA version 2) to an existing ocean biogeochemistry model consisting of the Finite-volumE Sea ice-Ocean Model (FESOM version 2.1) and the Regulated Ecosystem Model (REcoM version 3). Atmospheric CO2 in the model is a prognostic variable which is determined by the carbonate chemistry in the surface ocean. The model setup and its application to a pre-industrial control climate state is described in detail. In the coupled model 400 PgC are stored in equilibrium in the top 10 cm of the bioturbated sediment, mainly as calcite, but also to 5% as organic matter. Simulated atmospheric CO2 is in equilibrium at 286 ppm in the coupled simulation, which is close to the initially assumed value of the pre-industrial CO2 level. Sediment burial of carbon, alkalinity and nutrients in the coupled simulation is set to be partly compensated by riverine input. The spatial distribution of biological production is altered depending on the location of riverine input and the strength of local nutrient limitation, while the global productivity is not affected substantially. [ABSTRACT FROM AUTHOR]

In this paper we describe the implementation of the carbon isotopes 13C and 14C (radiocarbon) into the marine biogeochemistry model REcoM3. The implementation is tested in long-term equilibrium simulations where REcoM3 is coupled with the ocean general circulation model FESOM2.1, applying a low-resolution configuration and idealized climate forcing. Focusing on the carbon-isotopic composition of dissolved inorganic carbon (δ13CDIC and Δ14CDIC), our model results are largely consistent with reconstructions for the pre-anthropogenic period. Our simulations also exhibit discrepancies, e.g., in upwelling regions and the interior of the North Pacific. Some of these differences are due to the limitations of our ocean circulation model setup which results in a rather shallow meridional overturning circulation. We additionally study the accuracy of two simplified modelling approaches for dissolved inorganic 14C, which are faster (15 % and about a factor of five, respectively) than the complete consideration of the marine radiocarbon cycle. The accuracy of both simplified approaches is better than 5 % which should be sufficient for most studies of Δ14CDIC. [ABSTRACT FROM AUTHOR]

The cycling of carbon in the oceans is affected by feedbacks driven by changes in climate and atmospheric CO2. Understanding these feedbacks is therefore an important prerequisite for projecting future climate. Marine biogeochemistry models are a useful tool but, as with any model, are a simplification and need to be continually improved. In this study, we coupled the Finite-volumE Sea ice–Ocean Model (FESOM2.1) to the Regulated Ecosystem Model version 3 (REcoM3). FESOM2.1 is an update of the Finite-Element Sea ice–Ocean Model (FESOM1.4) and operates on unstructured meshes. Unlike standard structured-mesh ocean models, the mesh flexibility allows for a realistic representation of small-scale dynamics in key regions at an affordable computational cost. Compared to the previous coupled model version of FESOM1.4–REcoM2, the model FESOM2.1–REcoM3 utilizes a new dynamical core, based on a finite-volume discretization instead of finite elements, and retains central parts of the biogeochemistry model. As a new feature, carbonate chemistry, including water vapour correction, is computed by mocsy 2.0. Moreover, REcoM3 has an extended food web that includes macrozooplankton and fast-sinking detritus. Dissolved oxygen is also added as a new tracer. In this study, we assess the ocean and biogeochemical state simulated with FESOM2.1–REcoM3 in a global set-up at relatively low spatial resolution forced with JRA55-do atmospheric reanalysis. The focus is on the recent period (1958–2021) to assess how well the model can be used for present-day and future climate change scenarios on decadal to centennial timescales. A bias in the global ocean–atmosphere preindustrial CO2 flux present in the previous model version (FESOM1.4–REcoM2) could be significantly reduced. In addition, the computational efficiency is 2–3 times higher than that of FESOM1.4–REcoM2. Overall, it is found that FESOM2.1–REcoM3 is a skilful tool for ocean biogeochemical modelling applications. [ABSTRACT FROM AUTHOR]

The radiocarbon analysis of uranium‐thorium‐dated cold‐water corals (CWCs) provides an excellent opportunity for qualitative reconstruction of past ocean circulation and water mass aging. While mid‐depth water mass aging has been studied in the Atlantic Ocean, the evolution of the thermocline is still largely unknown. Here we present a combined 14C and 230Th/U age record obtained from thermocline dwelling CWCs at various sites in the eastern Atlantic Ocean, with intermittently centennial resolution over the last 32 ka. Shallow dwelling CWCs off Angola, located in the South Atlantic, infer a link between the mid‐depth equatorial Atlantic and Southern Ocean. They confirm a 14C drawdown during the Last Glacial Maximum (LGM) and advocate for a consistent Southern Hemisphere radiocarbon aging of upper thermocline waters, as well as strong depth gradients and high variability. Direct comparison with 14C simulations carried out with an Ocean General Circulation Model yield good agreement for Angola. In contrast, the North Atlantic thermocline shows well‐ventilated water with strong variations near the position of today's Azores Front (AF), neither of which are captured by the model. During the Bølling‐Allerød, we confirm the important role of the AF in separating North and South Atlantic thermocline waters and provide further evidence of a 500 year long deep convection interruption within the Younger Dryas (YD). We conclude that the North and South Atlantic thermocline waters were separately acting carbon reservoirs during the LGM and subsequent deglaciation until the modern circulation was established during the YD. Key Points: North Atlantic cold‐water corals trace well‐ventilated thermocline waters near major oceanic fronts since the Last Glacial MaximumAcross the South Atlantic into the Southern Ocean, aged waters with large variability and connectivity are evident during the last glacialThe modern state of radiocarbon ventilation of the thermocline Atlantic is initiated during the Younger Dryas cold reversal [ABSTRACT FROM AUTHOR]

Since it is currently not understood how changes in 14C production rate (Q), and in the carbon cycle, can be combined to explain the reconstructed atmospheric Δ14C record, we discuss possible reasons for this knowledge gap. Reviewing the literature, we exclude that changes in the content of atoms in the atmosphere, which produce cosmogenic 14C after being hit by galactic cosmic rays, might be responsible for parts of the observed differences. When combining Q with carbon cycle changes, one needs to understand the changes in the atmospheric 14C inventory, which are partially counterintuitive. For example, during the Last Glacial Maximum, Δ14C was ∼400‰ higher compared with preindustrial times, but the 14C inventory was 10% smaller. Some pronounced changes in atmospheric Δ14C do not correspond to any significant changes in the atmospheric 14C inventory, since CO2 was changing simultaneously. Using two conceptually different models (BICYCLE‐SE and LSG‐OGCM), we derive hypothetical Qs by forcing the models with identical atmospheric CO2 and Δ14C data. Results are compared with the most recent data‐based estimates of Q derived from cosmogenic isotopes. Millennial‐scale climate change connected to the bipolar seesaw is missing in the applied models, which might explain some, but probably not all, of the apparent model‐data disagreement in Q. Furthermore, Q based on either data from marine sediments or ice cores contains offsets, suggesting an interpretation deficit in the current data‐based approaches. Key Points: No important change found in the level of precursor material, which produces cosmogenic radiocarbon after being hit by galactic cosmic raysTransient changes in the bipolar seesaw are needed in applied carbon cycle models to improve interpretation of the radiocarbon cycleSediment core‐ and ice core‐based radiocarbon production rates differ systematically, suggesting missing processes [ABSTRACT FROM AUTHOR]

Atlantic cod (Gadus morhua) is one of the most commercially important fish species in the North Atlantic. Environmental factors, such as water temperatures, influence growth of individuals over time, thus forming population-specific growth patterns across climatic regions. Here we develop an integrative approach to investigate the role of temperature in shaping geographic differences of cod growth in the Celtic Sea, North Sea, Iceland, and Barents Sea. We combine a physiology-based growth model and 50-years observational temperature data of 0.5 × 0.5° spatial resolution to simulate continuous growth of cod. The model generated weight-at-age data for the period 1959–2007 which we compared to observational data from fishery-independent scientific surveys. In the Celtic and the northern North Sea, simulated growth matches well observational data. We also show that relatively warm temperatures in the Celtic Sea facilitate maximum growth rates; future warming is likely to have a negative impact on growth of these cod stocks. Growth simulations in Icelandic waters and the Barents Sea are less consistent with local observational data. More complex growth patterns in these regions are probably shaped by ontogenetic shifts in temperature regimes, feeding conditions and physiological adaptations. These findings should stimulate further research on critical processes to be considered in population-specific projections of growth of cod and productivity. [ABSTRACT FROM AUTHOR]

The past provides evidence of abrupt climate shifts and changes in the frequency of climate and weather extremes. We explore the nonlinear response to orbital forcing and then consider climate millennial variability down to daily weather events. Orbital changes are translated into regional responses in temperature, where the precessional response is related to nonlinearities and seasonal biases in the system. We question regularities found in climate events by analyzing the distribution of interevent waiting times. Periodicities of about 900 and 1,150 yr are found in ice cores besides the prominent 1,500 yr cycle. However, the variability remains indistinguishable from a random process, suggesting that centennial‐to‐millennial variability is stochastic in nature. New numerical techniques are developed allowing for a high resolution in the dynamically relevant regions like coasts, major upwelling regions, and high latitudes. Using this model, we find a strong sensitivity of the Atlantic meridional overturning circulation depending on where the deglacial meltwater is injected into. Meltwater into the Mississippi and near Labrador hardly affect the large‐scale ocean circulation, whereas subpolar hosing mimicking icebergs yields a quasi shutdown. The same multiscale approach is applied to radiocarbon simulations enabling a dynamical interpretation of marine sediment cores. Finally, abrupt climate events also have counterparts in the recent climate records, revealing a close link between climate variability, the statistics of North Atlantic weather patterns, and extreme events. Plain Language Summary: Predicting the future spread of possible climates, the risk of climate extremes and the risk of rapid transitions is of high socioeconomic relevance. The past provides evidence of abrupt climate change and the frequency of extremes. This allows to separate anthropogenic signals from natural climate variability. Earth system models applied both to past and future scenarios will enhance our ability to detect regime shifts which are necessary to potentially predict climate extremes and transitions. We consider the response of the system to regular orbital forcing and then focus on shorter time scales down to weather. The appearance of precession is linked to nonlinear responses of the climate system to external orbital forcing. Furthermore, we find that centennial‐to‐millennial variability is stochastic in nature. We also discuss recent developments of climate models with superior resolution in typical retrieval regions of paleoclimate records, such as continental margins and coasts. Using this model, we find a strong sensitivity of the Atlantic meridional overturning circulation depending on where the deglacial meltwater is injected into. Meltwater into the Mississippi and near Labrador hardly affect the large‐scale ocean circulation, whereas subpolar hosing related to icebergs yields a quasi shutdown. Our multiscale approach is applied to radiocarbon simulations enabling a dynamical interpretation of marine sediment cores. Finally, we explore a close link between climate variability, the statistics of weather patterns, and extreme events. Key Points: Nonlinear response of the climate system to external orbital forcing yields responses in the precessional bandCentennial and millennial variability are stochastic in natureOcean and radiocarbon modeling with high resolution at the coasts and high latitudes provides a suitable framework for multi‐scale dynamicsThe past provides evidences of abrupt climate changes and weather extremes [ABSTRACT FROM AUTHOR]

The concentration of radiocarbon (14C) differs between ocean and atmosphere. Radiocarbon determinations from samples which obtained their 14C in the marine environment therefore need a marine-specific calibration curve and cannot be calibrated directly against the atmospheric-based IntCal20 curve. This paper presents Marine20, an update to the internationally agreed marine radiocarbon age calibration curve that provides a non-polar global-average marine record of radiocarbon from 0–55 cal kBP and serves as a baseline for regional oceanic variation. Marine20 is intended for calibration of marine radiocarbon samples from non-polar regions; it is not suitable for calibration in polar regions where variability in sea ice extent, ocean upwelling and air-sea gas exchange may have caused larger changes to concentrations of marine radiocarbon. The Marine20 curve is based upon 500 simulations with an ocean/atmosphere/biosphere box-model of the global carbon cycle that has been forced by posterior realizations of our Northern Hemispheric atmospheric IntCal20 14C curve and reconstructed changes in CO2 obtained from ice core data. These forcings enable us to incorporate carbon cycle dynamics and temporal changes in the atmospheric 14C level. The box-model simulations of the global-average marine radiocarbon reservoir age are similar to those of a more complex three-dimensional ocean general circulation model. However, simplicity and speed of the box model allow us to use a Monte Carlo approach to rigorously propagate the uncertainty in both the historic concentration of atmospheric 14C and other key parameters of the carbon cycle through to our final Marine20 calibration curve. This robust propagation of uncertainty is fundamental to providing reliable precision for the radiocarbon age calibration of marine based samples. We make a first step towards deconvolving the contributions of different processes to the total uncertainty; discuss the main differences of Marine20 from the previous age calibration curve Marine13; and identify the limitations of our approach together with key areas for further work. The updated values for ΔR, the regional marine radiocarbon reservoir age corrections required to calibrate against Marine20, can be found at the data base http://calib.org/marine/. [ABSTRACT FROM AUTHOR]

Radiocarbon (14C) ages cannot provide absolutely dated chronologies for archaeological or paleoenvironmental studies directly but must be converted to calendar age equivalents using a calibration curve compensating for fluctuations in atmospheric 14C concentration. Although calibration curves are constructed from independently dated archives, they invariably require revision as new data become available and our understanding of the Earth system improves. In this volume the international 14C calibration curves for both the Northern and Southern Hemispheres, as well as for the ocean surface layer, have been updated to include a wealth of new data and extended to 55,000 cal BP. Based on tree rings, IntCal20 now extends as a fully atmospheric record to ca. 13,900 cal BP. For the older part of the timescale, IntCal20 comprises statistically integrated evidence from floating tree-ring chronologies, lacustrine and marine sediments, speleothems, and corals. We utilized improved evaluation of the timescales and location variable 14C offsets from the atmosphere (reservoir age, dead carbon fraction) for each dataset. New statistical methods have refined the structure of the calibration curves while maintaining a robust treatment of uncertainties in the 14C ages, the calendar ages and other corrections. The inclusion of modeled marine reservoir ages derived from a three-dimensional ocean circulation model has allowed us to apply more appropriate reservoir corrections to the marine 14C data rather than the previous use of constant regional offsets from the atmosphere. Here we provide an overview of the new and revised datasets and the associated methods used for the construction of the IntCal20 curve and explore potential regional offsets for tree-ring data. We discuss the main differences with respect to the previous calibration curve, IntCal13, and some of the implications for archaeology and geosciences ranging from the recent past to the time of the extinction of the Neanderthals. [ABSTRACT FROM AUTHOR]

Beyond ~13.9 cal kBP, the IntCal20 radiocarbon (14C) calibration curve is based upon combining data across a range of different archives including corals and planktic foraminifera. In order to reliably incorporate such marine data into an atmospheric curve, we need to resolve these records into their constituent atmospheric signal and marine reservoir age. We present results of marine reservoir age simulations enabling this resolution, applying the LSG ocean general circulation model forced with various climatic background conditions and with atmospheric radiocarbon changes according to the Hulu Cave speleothem record. Simulating the spatiotemporal evolution of reservoir ages between 54,000 and 10,700 cal BP, we find reservoir ages between 500 and 1400 yr in the low- and mid-latitudes, but also more than 3000 yr in the polar seas. Our results are broadly in agreement with available marine radiocarbon reconstructions, with the caveat that continental margins, marginal seas, or tropical lagoons are not properly resolved in our coarse-resolution model. [ABSTRACT FROM AUTHOR]

Ocean warming may lead to smaller body sizes of marine ectotherms, because metabolic rates increase exponentially with temperature while the capacity of the cardiorespiratory system to match enhanced oxygen demands is limited. Here, we explore the impact of rising sea water temperatures on Atlantic cod ( Gadus morhua), an economically important fish species. We focus on changes in the temperature-dependent growth potential by a transfer function model combining growth observations with climate model ensemble temperatures. Growth potential is expressed in terms of asymptotic body weight and depends on water temperature. We consider changes between the periods 1985-2004 and 2081-2100, assuming that future sea water temperatures will evolve according to climate projections for IPCC AR5 scenario RCP8.5. Our model projects a response of Atlantic cod to future warming, differentiated according to ocean regions, leading to increases of asymptotic weight in the Barents Sea, while weights are projected to decline at the southern margin of the biogeographic range. Southern spawning areas will disappear due to thermal limitation of spawning stages. These projections match the currently observed biogeographic shifts and the temperature- and oxygen-dependent decline in routine aerobic scope at southern distribution limits. [ABSTRACT FROM AUTHOR]

Rapid climate change in the Northeast Atlantic and Arctic poses a threat to some of the world's largest fish populations. Impacts of warming and acidification may become accessible through mechanism-based risk assessments and projections of future habitat suitability. We show that ocean acidification causes a narrowing of embryonic thermal ranges, which identifies the suitability of spawning habitats as a critical life-history bottleneck for two abundant cod species. Embryonic tolerance ranges linked to climate simulations reveal that ever-increasing CO 2 emissions [Representative Concentration Pathway (RCP) 8.5] will deteriorate suitability of present spawning habitat for both Atlantic cod ( Gadus morhua ) and Polar cod ( Boreogadus saida ) by 2100. Moderate warming (RCP4.5) may avert dangerous climate impacts on Atlantic cod but still leaves few spawning areas for the more vulnerable Polar cod, which also loses the benefits of an ice-covered ocean. Emissions following RCP2.6, however, support largely unchanged habitat suitability for both species, suggesting that risks are minimized if warming is held "below 2°C, if not 1.5°C," as pledged by the Paris Agreement.