scholarly journals A probabilistic assessment of calcium carbonate export and dissolution in the modern ocean

2015 ◽  
Vol 12 (24) ◽  
pp. 20223-20282 ◽  
Author(s):  
G. Battaglia ◽  
M. Steinacher ◽  
F. Joos

Abstract. The marine cycle of calcium carbonate (CaCO3) is an important element of the carbon cycle and co-governs the distribution of carbon and alkalinity within the ocean. However, CaCO3 fluxes and mechanisms governing CaCO3 dissolution are highly uncertain. We present an observationally-constrained, probabilistic assessment of the global and regional CaCO3 budgets. Parameters governing pelagic CaCO3 export fluxes and dissolution rates are sampled using a Latin-Hypercube scheme to construct a 1000 member ensemble with the Bern3D ocean model. Ensemble results are constrained by comparing simulated and observation-based fields of excess dissolved calcium carbonate (TA*). The minerals calcite and aragonite are modelled explicitly and ocean–sediment fluxes are considered. For local dissolution rates either a strong, a weak or no dependency on CaCO3 saturation is assumed. Median (68 % confidence interval) global CaCO3 export is 0.82 (0.67–0.98) Gt PIC yr−1, within the lower half of previously published estimates (0.4–1.8 Gt PIC yr−1). The spatial pattern of CaCO3 export is broadly consistent with earlier assessments. Export is large in the Southern Ocean, the tropical Indo–Pacific, the northern Pacific and relatively small in the Atlantic. Dissolution within the 200 to 1500 m depth range (0.33; 0.26–0.40 Gt PIC yr−1) is substantially lower than inferred from the TA*-CFC age method (1 ± 0.5 Gt PIC yr−1). The latter estimate is likely biased high as the TA*-CFC method neglects transport. The constrained results are robust across a range of diapycnal mixing coefficients and, thus, ocean circulation strengths. Modelled ocean circulation and transport time scales for the different setups were further evaluated with CFC11 and radiocarbon observations. Parameters and mechanisms governing dissolution are hardly constrained by either the TA* data or the current compilation of CaCO3 flux measurements such that model realisations with and without saturation-dependent dissolution achieve skill. We suggest to apply saturation-independent dissolution rates in Earth System Models to minimise computational costs.

2016 ◽  
Vol 13 (9) ◽  
pp. 2823-2848 ◽  
Author(s):  
Gianna Battaglia ◽  
Marco Steinacher ◽  
Fortunat Joos

Abstract. The marine cycle of calcium carbonate (CaCO3) is an important element of the carbon cycle and co-governs the distribution of carbon and alkalinity within the ocean. However, CaCO3 export fluxes and mechanisms governing CaCO3 dissolution are highly uncertain. We present an observationally constrained, probabilistic assessment of the global and regional CaCO3 budgets. Parameters governing pelagic CaCO3 export fluxes and dissolution rates are sampled using a Monte Carlo scheme to construct a 1000-member ensemble with the Bern3D ocean model. Ensemble results are constrained by comparing simulated and observation-based fields of excess dissolved calcium carbonate (TA*). The minerals calcite and aragonite are modelled explicitly and ocean–sediment fluxes are considered. For local dissolution rates, either a strong or a weak dependency on CaCO3 saturation is assumed. In addition, there is the option to have saturation-independent dissolution above the saturation horizon. The median (and 68 % confidence interval) of the constrained model ensemble for global biogenic CaCO3 export is 0.90 (0.72–1.05) Gt C yr−1, that is within the lower half of previously published estimates (0.4–1.8 Gt C yr−1). The spatial pattern of CaCO3 export is broadly consistent with earlier assessments. Export is large in the Southern Ocean, the tropical Indo–Pacific, the northern Pacific and relatively small in the Atlantic. The constrained results are robust across a range of diapycnal mixing coefficients and, thus, ocean circulation strengths. Modelled ocean circulation and transport timescales for the different set-ups were further evaluated with CFC11 and radiocarbon observations. Parameters and mechanisms governing dissolution are hardly constrained by either the TA* data or the current compilation of CaCO3 flux measurements such that model realisations with and without saturation-dependent dissolution achieve skill. We suggest applying saturation-independent dissolution rates in Earth system models to minimise computational costs.


2007 ◽  
Vol 20 (13) ◽  
pp. 2978-2993 ◽  
Author(s):  
Tommy G. Jensen

Abstract Composites of Florida State University winds (1970–99) for four different climate scenarios are used to force an Indian Ocean model. In addition to the mean climatology, the cases include La Niña, El Niño, and the Indian Ocean dipole (IOD). The differences in upper-ocean water mass exchanges between the Arabian Sea and the Bay of Bengal are investigated and show that, during El Niño and IOD years, the average clockwise Indian Ocean circulation is intensified, while it is weakened during La Niña years. As a consequence, high-salinity water export from the Arabian Sea into the Bay of Bengal is enhanced during El Niño and IOD years, while transport of low-salinity waters from the Bay of Bengal into the Arabian Sea is enhanced during La Niña years. This provides a venue for interannual salinity variations in the northern Indian Ocean.


2012 ◽  
Vol 25 (6) ◽  
pp. 1814-1826 ◽  
Author(s):  
Dimitrios Giannakis ◽  
Andrew J. Majda

Abstract An information-theoretic framework is developed to assess the predictive skill and model error in imperfect climate models for long-range forecasting. Here, of key importance is a climate equilibrium consistency test for detecting false predictive skill, as well as an analogous criterion describing model error during relaxation to equilibrium. Climate equilibrium consistency enforces the requirement that long-range forecasting models should reproduce the climatology of prediction observables with high fidelity. If a model meets both climate consistency and the analogous criterion describing model error during relaxation to equilibrium, then relative entropy can be used as an unbiased superensemble measure of the model’s skill in long-range coarse-grained forecasts. As an application, the authors investigate the error in modeling regime transitions in a 1.5-layer ocean model as a Markov process and identify models that are strongly persistent but their predictive skill is false. The general techniques developed here are also useful for estimating predictive skill with model error for Markov models of low-frequency atmospheric regimes.


2013 ◽  
Vol 2013 ◽  
pp. 1-18 ◽  
Author(s):  
Yanyun Liu ◽  
Lian Xie ◽  
John M. Morrison ◽  
Daniel Kamykowski

The regional impact of global climate change on the ocean circulation around the Galápagos Archipelago is studied using the Hybrid Coordinate Ocean Model (HYCOM) configured for a four-level nested domain system. The modeling system is validated and calibrated using daily atmospheric forcing derived from the NCEP/NCAR reanalysis dataset from 1951 to 2007. The potential impact of future anthropogenic global warming (AGW) in the Galápagos region is examined using the calibrated HYCOM with forcing derived from the IPCC-AR4 climate model. Results show that although the oceanic variability in the entire Galápagos region is significantly affected by global climate change, the degree of such effects is inhomogeneous across the region. The upwelling region to the west of the Isabella Island shows relatively slower warming trends compared to the eastern Galápagos region. Diagnostic analysis suggests that the variability in the western Galápagos upwelling region is affected mainly by equatorial undercurrent (EUC) and Panama currents, while the central/east Galápagos is predominantly affected by both Peru and EUC currents. The inhomogeneous responses in different regions of the Galápagos Archipelago to future AGW can be explained by the incoherent changes of the various current systems in the Galápagos region as a result of global climate change.


2005 ◽  
Vol 35 (1) ◽  
pp. 13-32 ◽  
Author(s):  
A. Birol Kara ◽  
Alan J. Wallcraft ◽  
Harley E. Hurlburt

Abstract A 1/25° × 1/25° cos(lat) (longitude × latitude) (≈3.2-km resolution) eddy-resolving Hybrid Coordinate Ocean Model (HYCOM) is introduced for the Black Sea and used to examine the effects of ocean turbidity on upper-ocean circulation features including sea surface height and mixed layer depth (MLD) on annual mean climatological time scales. The model is a primitive equation model with a K-profile parameterization (KPP) mixed layer submodel. It uses a hybrid vertical coordinate that combines the advantages of isopycnal, σ, and z-level coordinates in optimally simulating coastal and open-ocean circulation features. This model approach is applied to the Black Sea for the first time. HYCOM uses a newly developed time-varying solar penetration scheme that treats attenuation as a continuous quantity. This scheme includes two bands of solar radiation penetration, one that is needed in the top 10 m of the water column and another that penetrates to greater depths depending on the turbidity. Thus, it is suitable for any ocean general circulation model that has fine vertical resolution near the surface. With this scheme, the optical depth–dependent attenuation of subsurface heating in HYCOM is given by monthly mean fields for the attenuation of photosynthetically active radiation (kPAR) during 1997–2001. These satellite-based climatological kPAR fields are derived from Sea-Viewing Wide Field-of-View Sensor (SeaWiFS) data for the spectral diffuse attenuation coefficient at 490 nm (k490) and have been processed to have the smoothly varying and continuous coverage necessary for use in the Black Sea model applications. HYCOM simulations are driven by two sets of high-frequency climatological forcing, but no assimilation of ocean data is then used to demonstrate the importance of including spatial and temporal varying attenuation depths for the annual mean prediction of upper-ocean quantities in the Black Sea, which is very turbid (kPAR > 0.15 m−1, in general). Results are reported from three model simulations driven by each atmospheric forcing set using different values for the kPAR. A constant solar-attenuation optical depth of ≈17 m (clear water assumption), as opposed to using spatially and temporally varying attenuation depths, changes the surface circulation, especially in the eastern Black Sea. Unrealistic sub–mixed layer heating in the former results in weaker stratification at the base of the mixed layer and a deeper MLD than observed. As a result, the deep MLD off Sinop (at around 42.5°N, 35.5°E) weakens the surface currents regardless of the atmospheric forcing used in the model simulations. Using the SeaWiFS-based monthly turbidity climatology gives a shallower MLD with much stronger stratification at the base and much better agreement with observations. Because of the high Black Sea turbidity, the simulation with all solar radiation absorbed at the surface case gives results similar to the simulations using turbidity from SeaWiFS in the annual means, the aspect of the results investigated in this paper.


2006 ◽  
Vol 3 (3) ◽  
pp. 637-669 ◽  
Author(s):  
S. Natale ◽  
R. Sorgente ◽  
S. Gaberšek ◽  
A. Ribotti ◽  
A. Olita

Abstract. Ocean forecasts over the Central Mediterranean, produced by a near real time regional scale system, have been evaluated in order to assess their predictability. The ocean circulation model has been forced at the surface by a medium, high or very high resolution atmospheric forcing. The simulated ocean parameters have been compared with satellite data and they were found to be generally in good agreement. High and very high resolution atmospheric forcings have been able to form noticeable, although short-lived, surface current structures, due to their ability to detect transient atmospheric disturbances. The existence of the current structures has not been directly assessed due to lack of measurements. The ocean model in the slave mode was not able to develop dynamics different from the driving coarse resolution model which provides the boundary conditions.


MAUSAM ◽  
2021 ◽  
Vol 48 (4) ◽  
pp. 657-668
Author(s):  
XIAOMING LIU ◽  
JOHN M. MORRISON ◽  
LIAN XIE

Two sets of atmospheric forcing from NCEP/NCAR 40-year reanalysis project, one based on monthly averaged climatological data and the other on 1982-83 monthly averaged data, are used to derive the global Miami Isopycnic Coordinate Ocean Model (MICOM). These two runs are referred to as the climatological experiments and 1982-83 El Nino experiments. Sensitivity tests of tropical Pacific SST to different bulk parameterizations of air-sea heat and momentum fluxes are carried out in the two experiments. Primary results show that constant transfer coefficients                          (1.2 × 10-3) for heat flux greatly overestimate the tropical Pacific SST, whereas the Liu-Katsaros-Businger (Liu et al. 1979) method can significantly improve the SST simulation especially under very low-wind speed conditions. On the other hand, Large and Pond (1982) formulation of the drag coefficient made little difference on the tropical Pacific SST simulation although it might modify the surface ocean circulation. The SST seasonal cycle and interannual variability of tropical Pacific SST are also examined in this study. Since SST is the most important oceanic parameter that provides the link between the atmosphere and the ocean, this evaluation of different parameterization schemes may facilitate future studies on coupling ocean-atmospheric numeric models.    


2020 ◽  
Vol 13 (11) ◽  
pp. 5465-5483
Author(s):  
Clément Bricaud ◽  
Julien Le Sommer ◽  
Gurvan Madec ◽  
Christophe Calone ◽  
Julie Deshayes ◽  
...  

Abstract. Ocean biogeochemical models are key tools for both scientific and operational applications. Nevertheless the cost of these models is often expensive because of the large number of biogeochemical tracers. This has motivated the development of multi-grid approaches where ocean dynamics and tracer transport are computed on grids of different spatial resolution. However, existing multi-grid approaches to tracer transport in ocean modelling do not allow the computation of ocean dynamics and tracer transport simultaneously. This paper describes a new multi-grid approach developed for accelerating the computation of passive tracer transport in the Nucleus for European Modelling of the Ocean (NEMO) ocean circulation model. In practice, passive tracer transport is computed at runtime on a grid with coarser spatial resolution than the hydrodynamics, which reduces the CPU cost of computing the evolution of tracers. We describe the multi-grid algorithm, its practical implementation in the NEMO ocean model, and discuss its performance on the basis of a series of sensitivity experiments with global ocean model configurations. Our experiments confirm that the spatial resolution of hydrodynamical fields can be coarsened by a factor of 3 in both horizontal directions without significantly affecting the resolved passive tracer fields. Overall, the proposed algorithm yields a reduction by a factor of 7 of the overhead associated with running a full biogeochemical model like PISCES (with 24 passive tracers). Propositions for further reducing this cost without affecting the resolved solution are discussed.


2012 ◽  
Vol 58 (212) ◽  
pp. 1227-1244 ◽  
Author(s):  
Carl V. Gladish ◽  
David M. Holland ◽  
Paul R. Holland ◽  
Stephen F. Price

AbstractA numerical model for an interacting ice shelf and ocean is presented in which the ice- shelf base exhibits a channelized morphology similar to that observed beneath Petermann Gletscher’s (Greenland) floating ice shelf. Channels are initiated by irregularities in the ice along the grounding line and then enlarged by ocean melting. To a first approximation, spatially variable basal melting seaward of the grounding line acts as a steel-rule die or a stencil, imparting a channelized form to the ice base as it passes by. Ocean circulation in the region of high melt is inertial in the along-channel direction and geostrophically balanced in the transverse direction. Melt rates depend on the wavelength of imposed variations in ice thickness where it enters the shelf, with shorter wavelengths reducing overall melting. Petermann Gletscher’s narrow basal channels may therefore act to preserve the ice shelf against excessive melting. Overall melting in the model increases for a warming of the subsurface water. The same sensitivity holds for very slight cooling, but for cooling of a few tenths of a degree a reorganization of the spatial pattern of melting leads, surprisingly, to catastrophic thinning of the ice shelf 12 km from the grounding line. Subglacial discharge of fresh water along the grounding line increases overall melting. The eventual steady state depends on when discharge is initiated in the transient history of the ice, showing that multiple steady states of the coupled system exist in general.


2021 ◽  
Author(s):  
Chen Zhao ◽  
Rupert Gladstone ◽  
Ben Galton-Fenzi ◽  
David Gwyther

<p>The ocean-driven basal melting has important implications for the stability of ice shelves in Antarctic, which largely affects the ice sheet mass balance, ocean circulation, and subsequently global sea level rise. Due to the limited observations in the ice shelf cavities, the couple ice sheet ocean models have been playing a critical role in examining the processes governing basal melting. In this study we use the Framework for Ice Sheet-Ocean Coupling (FISOC) to couple the Elmer/Ice full-stokes ice sheet model and the Regional Ocean Modeling System (ROMS) ocean model to model ice shelf/ocean interactions for an idealised three-dimensional domain. Experiments followed the coupled ice sheet–ocean experiments under the first phase of the Marine Ice Sheet–Ocean Model Intercomparison Project (MISOMIP1). A periodic pattern in the simulated mean basal melting rates is found to be highly consistent with the maximum barotropic stream function and also the grounding line retreat row by row,  which is likely to be related with the gyre break down near the grounding line caused by some non-physical instability events from the ocean bottom. Sensitivity tests are carried out, showing that this periodic pattern is not sensitive to the choice of couple time intervals and horizontal eddy viscosities but sensitive to vertical resolution in the ocean model, the chosen critical water column thickness in the wet-dry scheme, and the tracer properties for the nudging dry cells at the ice-ocean interface boundary. Further simulations are necessary to better explain the mechanism involved in the couple ice-ocean system, which is very significant for its application on the realistic ice-ocean systems in polar regions.</p>


Sign in / Sign up

Export Citation Format

Share Document