scholarly journals A global scavenging and circulation ocean model of thorium-230 and protactinium-231 with improved particle dynamics (NEMO–ProThorP 0.1)

2018 ◽  
Vol 11 (9) ◽  
pp. 3537-3556 ◽  
Author(s):  
Marco van Hulten ◽  
Jean-Claude Dutay ◽  
Matthieu Roy-Barman

Abstract. In this paper we set forth a 3-D ocean model of the radioactive trace isotopes 230Th and 231Pa. The interest arises from the fact that these isotopes are extensively used for investigating particle transport in the ocean and reconstructing past ocean circulation. The tracers are reversibly scavenged by biogenic and lithogenic particles. Our simulations of 230Th and 231Pa are based on the NEMO–PISCES ocean biogeochemistry general circulation model, which includes biogenic particles, namely small and big particulate organic carbon, calcium carbonate and biogenic silica. Small and big lithogenic particles from dust deposition are included in our model as well. Their distributions generally compare well with the small and big lithogenic particle concentrations from recent observations from the GEOTRACES programme, except for boundary nepheloid layers for which, as of today, there are no non-trivial prognostic models available on a global scale. Our simulations reproduce 230Th and 231Pa dissolved concentrations: they compare well with recent GEOTRACES observations in many parts of the ocean. Particulate 230Th and 231Pa concentrations are significantly improved compared to previous studies, but they are still too low because of missing particles from nepheloid layers. Our simulation reproduces the main characteristics of the 231Pa∕230Th ratio observed in the sediments and supports a moderate affinity of 231Pa to biogenic silica as suggested by recent observations relative to 230Th. Future model development may further improve understanding, especially when this will include a more complete representation of all particles, including different size classes, manganese hydroxides and nepheloid layers. This can be done based on our model as its source code is readily available.

2017 ◽  
Author(s):  
Marco van Hulten ◽  
Jean-Claude Dutay ◽  
Matthieu Roy-Barman

Abstract. In this paper, we set forth a 3-D ocean model of the radioactive trace isotopes 230Th and 231Pa. The interest arises from the fact that these isotopes are extensively used for investigating particle transport in the ocean and reconstructing past ocean circulation. The tracers are reversibly scavenged by biogenic and lithogenic particles. Our simulations of 230Th and 231Pa are based on the NEMO–PISCES ocean biogeochemistry general circulation model, which includes biogenic particles, namely small and big particulate organic carbon, calcium carbonate and biogenic silica. Small and big lithogenic particles from dust deposition are included in our model as well. Their distributions generally compare well with the small and big lithogenic particle concentrations from recent observations from the GEOTRACES programme, except for boundary nepheloid layers for which, as up to today, there are no non-trivial, prognostic models available on a global scale. Our simulations reproduce 230Th and 231Pa dissolved concentrations: they compare well with recent GEOTRACES observations in many parts of the ocean. Particulate 230Th and 231Pa concentrations are significantly improved compared to previous studies, but they are still too low because of missing particles from nepheloid layers. Our simulation reproduces the main characteristics of the 231Pa / 230Th ratio observed in the sediments, and supports a moderate affinity of 231Pa to biogenic silica as suggested by recent observations, relative to 230Th. Future model development may further improve understanding, especially when this will include a more complete representation of all particles, including different size classes, manganese hydroxides and nepheloid layers. This can be done based on our model, as its source code is readily available.


Ocean Science ◽  
2009 ◽  
Vol 5 (4) ◽  
pp. 461-473 ◽  
Author(s):  
P. Oddo ◽  
M. Adani ◽  
N. Pinardi ◽  
C. Fratianni ◽  
M. Tonani ◽  
...  

Abstract. A new numerical general circulation ocean model for the Mediterranean Sea has been implemented nested within an Atlantic general circulation model within the framework of the Marine Environment and Security for the European Area project (MERSEA, Desaubies, 2006). A 4-year twin experiment was carried out from January 2004 to December 2007 with two different models to evaluate the impact on the Mediterranean Sea circulation of open lateral boundary conditions in the Atlantic Ocean. One model considers a closed lateral boundary in a large Atlantic box and the other is nested in the same box in a global ocean circulation model. Impact was observed comparing the two simulations with independent observations: ARGO for temperature and salinity profiles and tide gauges and along-track satellite observations for the sea surface height. The improvement in the nested Atlantic-Mediterranean model with respect to the closed one is particularly evident in the salinity characteristics of the Modified Atlantic Water and in the Mediterranean sea level seasonal variability.


2005 ◽  
Vol 35 (6) ◽  
pp. 933-948 ◽  
Author(s):  
Henk A. Dijkstra ◽  
Wilbert Weijer

Abstract A study of the stability of the global ocean circulation is performed within a coarse-resolution general circulation model. Using techniques of numerical bifurcation theory, steady states of the global ocean circulation are explicitly calculated as parameters are varied. Under a freshwater flux forcing that is diagnosed from a reference circulation with Levitus surface salinity fields, the global ocean circulation has no multiple equilibria. It is shown how this unique-state regime transforms into a regime with multiple equilibria as the pattern of the freshwater flux is changed in the northern North Atlantic Ocean. In the multiple-equilibria regime, there are two branches of stable steady solutions: one with a strong northern overturning in the Atlantic and one with hardly any northern overturning. Along the unstable branch that connects both stable solution branches (here for the first time computed for a global ocean model), the strength of the southern sinking in the South Atlantic changes substantially. The existence of the multiple-equilibria regime critically depends on the spatial pattern of the freshwater flux field and explains the hysteresis behavior as found in many previous modeling studies.


2007 ◽  
Vol 37 (2) ◽  
pp. 203-213 ◽  
Author(s):  
Manfred Wenzel ◽  
Jens Schröter

Abstract The mass budget of the ocean in the period 1993–2003 is studied with a general circulation model. The model has a free surface and conserves mass rather than volume; that is, freshwater is exchanged with the atmosphere via precipitation and evaporation and inflow from land is taken into account. The mass is redistributed by the ocean circulation. Furthermore, the ocean’s volume changes by steric expansion with changing temperature and salinity. To estimate the mass changes, the ocean model is constrained by sea level measurements from the Ocean Topography Experiment (TOPEX)/Poseidon mission as well as by hydrographic data. The modeled ocean mass change within the years 2002–03 compares favorably to measurements from the Gravity Recovery and Climate Experiment (GRACE), and the evolution of the global mean sea level for the period 1993–2003 with annual and interannual variations can be reproduced to a 0.15-cm rms difference. Its trend has been measured as 3.37 mm yr−1 while the constrained model gives 3.34 mm yr−1 considering only the area covered by measurements (3.25 mm yr−1 for the total ocean). A steric rise of 2.50 mm yr−1 is estimated in this period, as is a gain in the ocean mass that is equivalent to an eustatic rise of 0.74 mm yr−1. The amplitude and phase (day of maximum value since 1 January) of the superimposed eustatic annual cycle are also estimated to be 4.6 mm and 278°, respectively. The corresponding values for the semiannual cycle are 0.42 mm and 120°. The trends in the eustatic sea level are not equally distributed. In the Atlantic Ocean (80°S–67°N) the eustatic sea level rises by 1.8 mm yr−1 and in the Indian Ocean (80°S–30°N) it rises by 1.4 mm yr−1, but it falls by −0.20 mm yr−1 in the Pacific Ocean (80°S–67°N). The latter is mainly caused by a loss of mass through transport divergence in the Pacific sector of the Antarctic Circumpolar Current (−0.42 Sv; Sv ≡ 109 kg s−1) that is not balanced by the net surface water supply. The consequence of this uneven eustatic rise is a shift of the oceanic center of mass toward the Atlantic Ocean and to the north.


2009 ◽  
Vol 6 (2) ◽  
pp. 1093-1127 ◽  
Author(s):  
P. Oddo ◽  
M. Adani ◽  
N. Pinardi ◽  
C. Fratianni ◽  
M. Tonani ◽  
...  

Abstract. A new numerical general circulation ocean model for the Mediterranean Sea has been implemented nested within an Atlantic general circulation model within the framework of the Marine Environment and Security for the European Area project (MERSEA, Desaubies, 2006). A 4-year twin experiment was carried out from January 2004 to December 2007 with two different models to evaluate the impact on the Mediterranean Sea circulation of open lateral boundary conditions in the Atlantic Ocean. One model considers a closed lateral boundary in a large Atlantic box and the other is nested in the same box in a global ocean circulation model. Impact was observed comparing the two simulations with independent observations: ARGO for temperature and salinity profiles and tide gauges and along-track satellite observations for the sea surface height. The improvement in the nested Atlantic-Mediterranean model with respect to the closed one is particularly evident in the salinity characteristics of the Modified Atlantic Water and in the Mediterranean sea level seasonal variability.


2012 ◽  
Vol 53 (60) ◽  
pp. 59-69 ◽  
Author(s):  
Patrick Heimbach ◽  
Martin Losch

AbstractWe investigate the sensitivity of sub-ice-shelf melt rates under Pine Island Ice Shelf, West Antarctica, to changes in the oceanic state using an adjoint ocean model that is capable of representing the flow in sub-ice-shelf cavities. The adjoint code is based on algorithmic differentiation (AD) of the Massachusetts Institute of Technology’s ocean general circulation model (MITgcm). The adjoint model was extended by adding into the AD process the corresponding sub-ice-shelf cavity code, which implements a three-equation thermodynamic melt-rate parameterization to infer heat and freshwater fluxes at the ice-shelf/ocean boundary. The inferred sensitivities reveal dominant timescales of 30–60 days over which the shelf exit is connected to the deep interior via advective processes. They exhibit rich three-dimensional time-evolving patterns that can be understood in terms of a combination of the buoyancy forcing by inflowing water masses, the cavity geometry and the effect of rotation and topography in steering the flow in the presence of prominent features in the bedrock bathymetry. Dominant sensitivity pathways are found over a sill, as well as ‘shadow regions’ of very low sensitivities. To the extent that these transient patterns are robust they carry important information for decision-making in observation deployment and monitoring.


2007 ◽  
Vol 37 (2) ◽  
pp. 277-295 ◽  
Author(s):  
A. J. Meijers ◽  
N. L. Bindoff ◽  
J. L. Roberts

Abstract The large-scale volume, heat, and freshwater ocean transports in the Southern Hemisphere are investigated using time-averaged output from a seasonless, high-resolution general circulation model. The ocean circulation is realistic, and property transports are comparable to observations. The Antarctic Circumpolar Current (ACC) carries 144 Sv (Sv ≡ 106 m3 s−1) of water eastward across Drake Passage, increasing to 155 Sv south of Australia because of the Indonesian Throughflow (ITF). There is a clear Indo-Pacific gyre around Australia exchanging −10 Sv, 0.9 PW of heat, and 0.2 Sv of freshwater through the ITF, and there is a 9-Sv leakage from the Tasman Sea to the Indian Ocean. The transport of heat and freshwater by eddies is localized to the upper 1000 m of the water column and specific regions, such as western boundary currents, confluences, and the subantarctic front (SAF). Eddy transport of heat and freshwater is negligible in gyre interiors and south of the SAF but is vital across the northern edge of the ACC, in particular at the Agulhas Retroflection where eddies accomplish almost 100% of the net ocean heat and 60% of the southward freshwater transport. The eddy transport is almost zero across the latitude of Drake Passage while in a quasi-Lagrangian frame eddy transports are significant across the ACC but surprisingly are still smaller than the mean transport of heat. Mean and eddy property transport divergences are found to be strongly compensating in areas of high eddy activity. This is caused by increased baroclinic instability in strong mean flows, which induces an opposing eddy transport. This relationship is observed to be stronger in the case of horizontal heat transport than in corresponding horizontal freshwater transports.


1997 ◽  
Vol 25 ◽  
pp. 111-115 ◽  
Author(s):  
Achim Stössel

This paper investigates the long-term impact of sea ice on global climate using a global sea-ice–ocean general circulation model (OGCM). The sea-ice component involves state-of-the-art dynamics; the ocean component consists of a 3.5° × 3.5° × 11 layer primitive-equation model. Depending on the physical description of sea ice, significant changes are detected in the convective activity, in the hydrographic properties and in the thermohaline circulation of the ocean model. Most of these changes originate in the Southern Ocean, emphasizing the crucial role of sea ice in this marginally stably stratified region of the world's oceans. Specifically, if the effect of brine release is neglected, the deep layers of the Southern Ocean warm up considerably; this is associated with a weakening of the Southern Hemisphere overturning cell. The removal of the commonly used “salinity enhancement” leads to a similar effect. The deep-ocean salinity is almost unaffected in both experiments. Introducing explicit new-ice thickness growth in partially ice-covered gridcells leads to a substantial increase in convective activity, especially in the Southern Ocean, with a concomitant significant cooling and salinification of the deep ocean. Possible mechanisms for the resulting interactions between sea-ice processes and deep-ocean characteristics are suggested.


2016 ◽  
Vol 12 (5) ◽  
pp. 1181-1198 ◽  
Author(s):  
Daniel J. Lunt ◽  
Alex Farnsworth ◽  
Claire Loptson ◽  
Gavin L. Foster ◽  
Paul Markwick ◽  
...  

Abstract. During the period from approximately 150 to 35 million years ago, the Cretaceous–Paleocene–Eocene (CPE), the Earth was in a “greenhouse” state with little or no ice at either pole. It was also a period of considerable global change, from the warmest periods of the mid-Cretaceous, to the threshold of icehouse conditions at the end of the Eocene. However, the relative contribution of palaeogeographic change, solar change, and carbon cycle change to these climatic variations is unknown. Here, making use of recent advances in computing power, and a set of unique palaeogeographic maps, we carry out an ensemble of 19 General Circulation Model simulations covering this period, one simulation per stratigraphic stage. By maintaining atmospheric CO2 concentration constant across the simulations, we are able to identify the contribution from palaeogeographic and solar forcing to global change across the CPE, and explore the underlying mechanisms. We find that global mean surface temperature is remarkably constant across the simulations, resulting from a cancellation of opposing trends from solar and palaeogeographic change. However, there are significant modelled variations on a regional scale. The stratigraphic stage–stage transitions which exhibit greatest climatic change are associated with transitions in the mode of ocean circulation, themselves often associated with changes in ocean gateways, and amplified by feedbacks related to emissivity and planetary albedo. We also find some control on global mean temperature from continental area and global mean orography. Our results have important implications for the interpretation of single-site palaeo proxy records. In particular, our results allow the non-CO2 (i.e. palaeogeographic and solar constant) components of proxy records to be removed, leaving a more global component associated with carbon cycle change. This “adjustment factor” is used to adjust sea surface temperatures, as the deep ocean is not fully equilibrated in the model. The adjustment factor is illustrated for seven key sites in the CPE, and applied to proxy data from Falkland Plateau, and we provide data so that similar adjustments can be made to any site and for any time period within the CPE. Ultimately, this will enable isolation of the CO2-forced climate signal to be extracted from multiple proxy records from around the globe, allowing an evaluation of the regional signals and extent of polar amplification in response to CO2 changes during the CPE. Finally, regions where the adjustment factor is constant throughout the CPE could indicate places where future proxies could be targeted in order to reconstruct the purest CO2-induced temperature change, where the complicating contributions of other processes are minimised. Therefore, combined with other considerations, this work could provide useful information for supporting targets for drilling localities and outcrop studies.


Sign in / Sign up

Export Citation Format

Share Document