LOSCAR: Long-term Ocean-atmosphere-Sediment CArbon cycle Reservoir Model v2.0.4

The LOSCAR model is designed to efficiently compute the partitioning of carbon between ocean, atmosphere, and sediments on time scales ranging from centuries to millions of years. While a variety of computationally inexpensive carbon cycle models are already available, many are missing a critical sediment component, which is indispensable for long-term integrations. One of LOSCAR's strengths is the coupling of ocean-atmosphere routines to a computationally efficient sediment module. This allows, for instance, adequate computation of CaCO<sub>3</sub> dissolution, calcite compensation, and long-term carbon cycle fluxes, including weathering of carbonate and silicate rocks. The ocean component includes various biogeochemical tracers such as total carbon, alkalinity, phosphate, oxygen, and stable carbon isotopes. LOSCAR's configuration of ocean geometry is flexible and allows for easy switching between modern and paleo-versions. We have previously published applications of the model tackling future projections of ocean chemistry and weathering, <i>p</i>CO<sub>2</sub> sensitivity to carbon cycle perturbations throughout the Cenozoic, and carbon/calcium cycling during the Paleocene-Eocene Thermal Maximum. The focus of the present contribution is the detailed description of the model including numerical architecture, processes and parameterizations, tuning, and examples of input and output. Typical CPU integration times of LOSCAR are of order seconds for several thousand model years on current standard desktop machines. The LOSCAR source code in C can be obtained from the author by sending a request to [email protected]

CO2 perturbation experiments: similarities and differences between dissolved inorganic carbon and total alkalinity manipulations

Increasing atmospheric carbon dioxide (CO2) through human activities and invasion of anthropogenic CO2 into the surface ocean alters the seawater carbonate chemistry, increasing CO2 and bicarbonate (HCO3−) at the expense of carbonate ion (CO32−) concentrations. This redistribution in the dissolved inorganic carbon (DIC) pool decreases pH and carbonate saturation state (Ω). Several components of the carbonate system are considered potential key variables influencing for instance calcium carbonate precipitation in marine calcifiers such as coccolithophores, foraminifera, corals, mollusks and echinoderms. Unravelling the sensitivities of marine organisms and ecosystems to CO2 induced ocean acidification (OA) requires well-controlled experimental setups and accurate carbonate system manipulations. Here we describe and analyse the chemical changes involved in the two basic approaches for carbonate chemistry manipulation, i.e. changing DIC at constant total alkalinity (TA) and changing TA at constant DIC. Furthermore, we briefly introduce several methods to experimentally manipulate DIC and TA. Finally, we examine responses obtained with both approaches using published results for the coccolithophore Emiliania huxleyi. We conclude that under most experimental conditions in the context of ocean acidification DIC and TA manipulations yield similar changes in all parameters of the carbonate system, which implies direct comparability of data obtained with the two basic approaches for CO2 perturbation

Early detection of ocean acidification effects on marine calcification

Ocean acidification is likely to impact calcification rates in many pelagic organisms, which may in turn cause significant changes in marine ecosystem structure. We examine effects of changes in marine CaCO3 production on total alkalinity (TA) in the ocean using the global biogeochemical ocean model HAMOCC. We test a variety of future calcification scenarios because experimental studies with different organisms have revealed a wide range of calcification sensitivities to CaCO3 saturation state. The model integrations start at a preindustrial steady state in the year 1800 and run until the year 2300 forced with anthropogenic CO2 emissions. Calculated trends in TA are evaluated taking into account the natural variability in ocean carbonate chemistry, as derived from repeat hydrographic transects. We conclude that the data currently available does not allow discerning significant trends in TA due to changes in pelagic calcification caused by ocean acidification. Given different calcification scenarios, our model calculations indicate that the TA increase over time will start being detectable by the year 2040, increasing by 5-30 mu mol kg(-1) compared to the present-day values. In a scenario of extreme reductions in calcification, large TA changes relative to preindustrial conditions would have occurred at present, which we consider very unlikely. However, the time interval of reliable TA observations is too short to disregard this scenario. The largest increase in surface ocean TA is predicted for the tropical and subtropical regions. In order to monitor and quantify possible early signs of acidification effects, we suggest to specifically target those regions during future ocean chemistry surveys

Ongoing transients in carbonate compensation

Uptake of anthropogenic CO2 is acidifying the oceans. Over the next 2000 years, this will modify the dissolution and preservation of sedimentary carbonate. By coupling new formulas for the positions of the calcite saturation horizon, zsat, the compensation depth, zcc, and the snowline, zsnow, to a biogeochemical model of the oceanic carbonate system, we evaluate how these horizons will change with ongoing ocean acidification. Our model is an extended Havardton-Bear-type box model, which includes novel kinetic descriptions for carbonate dissolution above, between, and below these critical depths. In the preindustrial ocean, zsat and zcc are at 3939 and 4750 m, respectively. When forced with the IS92a CO2 emission scenario, the model forecasts (1) that zsat will rise rapidly (“runaway” conditions) so that all deep water becomes undersaturated, (2) that zcc will also rise and over 1000 years will pass before it will be stabilized by the dissolution of previously deposited CaCO3, and (3) that zsnow will respond slowly to acidification, rising by ∼1150 m during a 2000 year timeframe. A further simplified model that equates the compensation and saturation depths produces quantitatively different results. Finally, additional feedbacks due to acidification on calcification and increased atmospheric CO2 on organic matter productivity strongly affect the positions of the compensation horizons and their dynamics.

Inorganic carbon concentrating mechanisms in free-living and symbiotic dinoflagellates and chromerids

Photosynthetic dinoflagellates are ecologically and biogeochemically important in marine and freshwater environments. However, surprisingly little is known of how this group acquires inorganic carbon or how these diverse processes evolved. Consequently, how CO2 availability ultimately influences the success of dinoflagellates over space and time remains poorly resolved compared to other microalgal groups. Here we review the evidence. Photosynthetic core dinoflagellates have a Form II RuBisCO (replaced by Form IB or Form ID in derived dinoflagellates). The in vitro kinetics of the Form II RuBisCO from dinoflagellates are largely unknown, but dinoflagellates with Form II (and other) RuBisCOs have inorganic carbon concentrating mechanisms (CCMs), as indicated by in vivo internal inorganic C accumulation and affinity for external inorganic C. However, the location of the membrane(s) at which the essential active transport component(s) of the CCM occur(s) is (are) unresolved; isolation and characterization of functionally competent chloroplasts would help in this respect. Endosymbiotic Symbiodiniaceae (in Foraminifera, Acantharia, Radiolaria, Ciliata, Porifera, Acoela, Cnidaria, and Mollusca) obtain inorganic C by transport from seawater through host tissue. In corals this transport apparently provides an inorganic C concentration around the photobiont that obviates the need for photobiont CCM. This is not the case for tridacnid bivalves, medusae, or, possibly, Foraminifera. Overcoming these long-standing knowledge gaps relies on technical advances (e.g., the in vitro kinetics of Form II RuBisCO) that can functionally track the fate of inorganic C forms

Marine biogeochemical responses to the North Atlantic Oscillation in a coupled climate model

In this study a coupled ocean-atmosphere model containing interactive marine biogeochemistry is used to analyze interannual, lagged, and decadal marine biogeochemical responses to the North Atlantic Oscillation (NAO), the dominant mode of North Atlantic atmospheric variability. The coupled model adequately reproduces present-day climatologies and NAO atmospheric variability. It is shown that marine biogeochemical responses to the NAO are governed by different mechanisms according to the time scale considered. On interannual time scales, local changes in vertical mixing, caused by modifications in air-sea heat, freshwater, and momentum fluxes, are most relevant in influencing phytoplankton growth through light and nutrient limitation mechanisms. At subpolar latitudes, deeper mixing occurring during positive NAO winters causes a slight decrease in late winter chlorophyll concentration due to light limitation and a 10%–20% increase in spring chlorophyll concentration due to higher nutrient availability. The lagged response of physical and biogeochemical properties to a high NAO winter shows some memory in the following 2 years. In particular, subsurface nutrient anomalies generated by local changes in mixing near the American coast are advected along the North Atlantic Current, where they are suggested to affect downstream chlorophyll concentration with 1 year lag. On decadal time scales, local and remote mechanisms act contemporaneously in shaping the decadal biogeochemical response to the NAO. The slow circulation adjustment, in response to NAO wind stress curl anomalies, causes a basin redistribution of heat, freshwater, and biogeochemical properties which, in turn, modifies the spatial structure of the subpolar chlorophyll bloom

Constraints on Earth system functioning at the Paleocene-Eocene Thermal Maximum from the marine silicon cycle

The Paleocene‐Eocene Thermal Maximum (PETM, ca. 56 Ma) is marked by a negative carbon isotope excursion (CIE) and increased global temperatures. The CIE is thought to result from the release of 13C‐depleted carbon, although the source(s) of carbon and triggers for its release, its rate of release, and the mechanisms by which the Earth system recovered are all debated. Many of the proposed mechanisms for the onset and recovery phases of the PETM make testable predictions about the marine silica cycle, making silicon isotope records a promising tool to address open questions about the PETM. We analyzed silicon isotope ratios (δ30Si) in radiolarian tests and sponge spicules from the Western North Atlantic (ODP Site 1051) across the PETM. Radiolarian δ30Si decreases by 0.6‰ from a background of 1‰ coeval with the CIE, while sponge δ30Si remains consistent at 0.2‰. Using a box model to test the Si cycle response to various scenarios, we find the data are best explained by a weak silicate weathering feedback, implying the recovery was mostly driven by nondiatom organic carbon burial, the other major long‐term carbon sink. We find no resolvable evidence for a volcanic trigger for carbon release, or for a change in regional oceanography. Better understanding of radiolarian Si isotope fractionation and more Si isotope records spanning the PETM are needed to confirm the global validity of these conclusions, but they highlight how the coupling between the silica and carbon cycles can be exploited to yield insight into the functioning of the Earth system