South Atlantic circulation in a world ocean model

The circulation in the South Atlantic Ocean has been simulated within a global ocean general circulation model. Preliminary analysis of the modelled ocean circulation in the region indicates a rather close agreement of the simulated upper ocean flows with conventional notions of the large-scale geostrophic currents in the region. The modelled South Atlantic Ocean witnesses the return flow and export of North Atlantic Deep Water (NADW) at its northern boundary, the inflow of a rather barotropic Antarctic Circumpolar Current (ACC) through the Drake Passage, and the inflow of warm saline Agulhas water around the Cape of Good Hope. The Agulhas leakage amounts to 8.7 Sv, within recent estimates of the mass transport shed westward at the Agulhas retroflection. Topographic steering of the ACC dominates the structure of flow in the circumpolar ocean. The Benguela Current is seen to be fed by a mixture of saline Indian Ocean water (originating from the Agulhas Current) and fresher Subantarctic surface water (originating in the ACC). The Benguela Current is seen to modify its flow and fate with depth; near the surface it flows north-westwards bifurcating most of its transport northward into the North Atlantic Ocean (for ultimate replacement of North Atlantic surface waters lost to the NADW conveyor). Deeper in the water column, more of the Benguela Current is destined to return with the Brazil Current, though northward flows are still generated where the Benguela Current extension encounters the coast of South America. At intermediate levels, these northward currents trace the flow of Antarctic Intermediate Water (AAIW) equatorward, though even more AAIW is seen to recirculate poleward in the subtropical gyre. In spite of the model’s rather coarse resolution, some subtle features of the Brazil-Malvinas Confluence are simulated rather well, including the latitude at which the two currents meet. Conceptual diagrams of the recirculation and interocean exchange of thermocline, intermediate and deep waters are constructed from an analysis of flows bound between isothermal and isobaric surfaces. This analysis shows how the return path of NADW is partitioned between a cold water route through the Drake Passage (6.5 Sv), a warm water route involving the Agulhas Current sheeding thermocline water westward (2.5 Sv), and a recirculation of intermediate water originating in the Indian Ocean (1.6 Sv).     

Dynamic model of mesoscale eddies. Eddy parameterization for coarse resolution ocean circulation models

In the framework of the eddy dynamic model developed in two previous papers (Dubovikov, M.S., Dynamical model of mesoscale eddies, Geophys. Astophys. Fluid Dyn., 2003, 97, 311–358; Canuto, V.M. and Dubovikov, M.S., Modeling mesoscale eddies, Ocean Modelling, 2004, 8, 1–30 referred as I–II), we compute the contribution of unresolved mesoscale eddies to the large-scale dynamic equations of the ocean. In isopycnal coordinates, in addition to the bolus velocity discussed in I–II, the mesoscale contribution to the large scale momentum equation is derived. Its form is quite different from the traditional down-gradient parameterization. The model solutions in isopycnal coordinates are transformed to level coordinates to parameterize the eddy contributions to the corresponding large scale density and momentum equations. In the former, the contributions due to the eddy induced velocity and to the residual density flux across mean isopycnals (so called Σ-term) are derived, both contributions being shown to be of the same order. As for the large scale momentum equation, as well as in isopycnal coordinates, the eddy contribution has a form which is quite different from the down-gradient expression.     

Equilibration mechanisms in an adjoint ocean general circulation model

We examine the equilibrated and time-evolving adjoint solutions of an ocean general circulation model. Adjoint models calculate the sensitivity of a diagnostic, (here, the strength of the meridional overturning) to all forcing fields in a single integration. The time evolution of the sensitivity patterns demonstrates the validity of the adjoint modeling approach over climatological time scales in coarse-resolution ocean models. Our objective is to identify the principle adjustment mechanisms through which the meridional overturning strength adapts to perturbations in wind and buoyancy forcing. The adjoint approach is shown to be a valuable alternative to traditional perturbation methods in highlighting the processes and time scales important to ocean and climate modeling.     

A source-Sink laboratory model of the ocean circulation

Abstract

A two gyre circulation and inertial western boundary currents have been observed in a sloping bottom laboratory model of a barotropic ocean circulation. Water of viscosity v is contained in a rotating (angular velocity ω), square basin of side L (30 cm) with a flat top and a bottom slope (tan θ) such that the depth (H) varies from 12 to 15 cm. The flow is driven by a distributed source and sink at the upper surface, a plate drilled with 342 holes. The hole distribution and size is arranged so that the average imposed vertical velocity, w = w ₀ sin (2πy′/30), models the Ekman divergence from a two gyre zonal wind stress. Fluid flow is observed with the thymol blue technique over the ranges of Rossby numbers (w ₀/2ωL tan θ) from 1.44 × 10^?3 to 1.41 × 10^?2 and Ekman numbers (v/2ωH ²) from 2.13 × 10^?5 to 2.10 × 10^?3. At the largest Rossby numbers the flow pattern changes markedly, but the non-uniformity of the imposed vertical velocity also penetrates deep into the fluid in this regime.     

开边界海盆尺度环流模式模拟北太平洋CFCs分布

利用一个开边界北太平洋海盆尺度环流模式模拟了北太平洋CFCs的吸收和分布.模式计算了CFC_11海气界面通量及其在海洋中的输运和再分布,考察了开边界对模式模拟结果的影响以及开边界模式中不同沿等密度面扩散系数的影响.打开南边界可以很好地改善赤道以南海区的CFC_11模拟.增加沿等密度面扩散系数使得137°E和165°E断面模拟的高CFC_11水向南、向下输送的强度和影响范围都有明显增大,能很好地改善西北太平洋的模拟结果.西北太平洋是CFC_11的一个重要的汇区,此外在赤道东太平洋上升流区还存在一个小的汇区.     

Processes affecting dissolved iron across the Subtropical North Atlantic: a model study

Ocean Dynamics - Trace metal measurements in recent years have revealed a complex distribution of dissolved iron (dFe) in the ocean that models still struggle to reproduce. The GEOTRACES section...     

Optimisation of a parallel ocean general circulation model

This paper presents the development of a general-purpose parallel ocean circulation model, for use on a wide range of computer platforms, from traditional scalar machines to workstation clusters and massively parallel processors. Parallelism is provided, as a modular option, via high-level message-passing routines, thus hiding the technical intricacies from the user. An initial implementation highlights that the parallel efficiency of the model is adversely affected by a number of factors, for which optimisations are discussed and implemented. The resulting ocean code is portable and, in particular, allows science to be achieved on local workstations that could otherwise only be undertaken on state-of-the-art supercomputers.     

Impact of partial steps and momentum advection schemes in a global ocean circulation model at eddy-permitting resolution

Series of sensitivity tests were performed with a z-coordinate, global eddy-permitting (1/4°) ocean/sea-ice model (the ORCA-R025 model configuration developed for the DRAKKAR project) to carefully evaluate the impact of recent state-of-the-art numerical schemes on model solutions. The combination of an energy–enstrophy conserving (EEN) scheme for momentum advection with a partial step (PS) representation of the bottom topography yields significant improvements in the mean circulation. Well known biases in the representation of western boundary currents, such as in the Atlantic the detachment of the Gulf Stream, the path of the North Atlantic Current, the location of the Confluence, and the strength of the Zapiola Eddy in the south Atlantic, are partly corrected. Similar improvements are found in the Pacific, Indian, and Southern Oceans, and characteristics of the mean flow are generally much closer to observations. Comparisons with other state-of-the-art models show that the ORCA-R025 configuration generally performs better at similar resolution. In addition, the model solution is often comparable to solutions obtained at 1/6 or 1/10° resolution in some aspects concerning mean flow patterns and distribution of eddy kinetic energy. Although the reasons for these improvements are not analyzed in detail in this paper, evidence is shown that the combination of EEN with PS reduces numerical noise near the bottom, which is likely to affect current–topography interactions in a systematic way. We conclude that significant corrections of the mean biases presently seen in general circulation model solutions at eddy-permitting resolution can still be expected from the development of numerical methods, which represent an alternative to increasing resolution.     

A non-breaking-wave-generated turbulence mixing scheme for a global ocean general circulation model

Ocean Dynamics - A novel vertical mixing scheme to describe the influence of the non-breaking surface waves in ocean general circulation models is proposed based on the second-order turbulence...     

Eddy-permitting simulations of the sub-polar North Atlantic: impact of the model bias on water mass properties and circulation

Some previous studies demonstrated that model bias has a strong impact on the quality of long-term prognostic model simulations of the sub-polar North Atlantic Ocean. Relatively strong bias of water mass characteristics is observed in both eddy-permitting and eddy-resolving simulations, suggesting that an increase of model resolution does not reduce significantly the model bias. This study is an attempt to quantify the impact of model bias on the simulated water mass and circulation characteristics in an eddy-permitting model of the sub-polar ocean. This is done through comparison of eddy-permitting prognostic model simulations with the results from two other runs in which the bias is constrained by using spectral nudging. In the first run, the temperature and salinity are nudged towards climatology in the whole column. In the second run, the spectral nudging is applied in the surface 30 m layer and at depths below 560 m only. The biases of the model characteristics of the unconstrained run are similar to those reported in previous eddy-permitting and eddy-resolving studies. The salinity in the surface and intermediate waters of the Labrador Sea waters increases with respect to the climatology, which reduces the stability of the water column. The deep convection in the unconstrained run is artificially intensified and the transport in the sub-polar gyre stronger than in the observations. In particular, the transport of relatively salty and warm Irminger waters into the Labrador Sea is unrealistically high. While the water mass temperature and salinity in the run with spectral nudging in the whole column are closest to the observations, the depth of the winter convection is underestimated in the model. The water mass characteristics and water transport in the run with spectral nudging in the surface and deep layers only are close to observations and at the same time represent well the deep convection in terms of its intensity and position. The source of the bias in the prognostic model run is discussed.     

Simulating the three-dimensional circulation and hydrography of Halifax Harbour using a multi-nested coastal ocean circulation model

Halifax Harbour is located on the Atlantic coast of Nova Scotia, Canada. It is one of the world’s largest, ice-free natural harbours and of great economic importance to the region. A good understanding of the physical processes controlling tides, flooding, transport and dispersion, and hydrographic variability is required for pollution control and sustainable development of the Harbour. For the first time, a multi-nested, finite difference coastal ocean circulation model is used to reconstruct the three-dimensional circulation and hydrography of the Harbour and its variability on timescales of hours to months for 2006. The model is driven by tides, wind and sea level pressure, air-sea fluxes of heat, and terrestrial buoyancy fluxes associated with river and sewage discharge. The predictive skill of the model is assessed by comparing the model simulations with independent observations of sea level from coastal tide gauges and currents from moored instruments. The simulated hydrography is also compared against a new monthly climatology created from all available temperature and salinity observations made in the Harbour over the last century. It is shown that the model can reproduce accurately the main features of the observed tides and storm surge, seasonal mean circulation and hydrography, and wind driven variations. The model is next used to examine the main physical processes controlling the circulation and hydrography of the Harbour. It is shown that non-linear interaction between tidal currents and complex topography occurs over the Narrows. The overall circulation can be characterized as a two-layer estuarine circulation with seaward flow in the thin upper layer and landward flow in the broad lower layer. An important component of this estuarine circulation is a relatively strong, vertically sheared jet situated over a narrow sill connecting the inner Harbour to the deep and relatively quiescent Bedford Basin. Local wind driven variability is strongest in winter as expected but it is also shown that a significant part of the temperature and salinity variability is driven by physical processes occurring on the adjacent inner continental shelf, especially during storm and coastal upwelling events.     

Process study of coastal circulation over the inner Scotian Shelf using a nested-grid ocean circulation model, with a special emphasis on the storm-induced circulation during tropical storm Alberto in 2006

This study examines main physical processes affecting the three-dimensional (3D) circulation and hydrographic distributions over the inner Scotian Shelf (ISS) in June and July 2006 using a nested-grid coastal ocean circulation modeling system known as the NCOPS-LB. The nested-grid system has five relocatable downscaling submodels, with the outermost submodel of a coarse horizontal resolution of (1/12)° for simulating storm surges and barotropic shelf waves over the Eastern Canadian shelf and the innermost submodel of a fine resolution of ～180 m for simulating the 3D coastal circulation and hydrography over Lunenburg Bay of Nova Scotia in the default setup. The NCOPS-LB is driven by meteorological and astronomical forcing and used to study the storm-induced circulation over the ISS during tropical storm Alberto. Model results demonstrate that the coastal circulation and hydrographic distributions over the ISS are affected significantly by tides, local wind forcing, and remotely generated coastal waves during the study period.     

Dissolution rates of calcium carbonate in the deep ocean; an in-situ experiment in the North Atlantic Ocean

An in-situ water circulator (ISWAC) was developed to allow accurate measurement of dissolution rates of various carbonate particles with minimal stagnation and mechanical weight loss. Three ISWAC packages were deployed at 3600, 4800 and 5518 m for 79 days in the Sargasso Sea (Northwest Atlantic). Weight losses for different particles during 79 days at 5518 m were as follows: pteropods and synthetic aragonite, 72.8%; reagent calcite, 57.5%; foraminifera, 23–36%; coccoliths, 11.3–24%, and diatoms, 12%. These weight losses are 2.5–7.5 times higher than those reported in earlier in-situ experiments. Normalization of weight losses with respect to BET specific surface area for different calcite particles yielded specific dissolution rates that differed by more than 2 orders of magnitude. Bleached biogenic particles dissolve significantly faster than non-bleached although their surface area is identical. We suggest that the BET surface area does not represent the reactive surface area available for dissolution, especially in biogenic calcite particles. Coatings, probably of organic matter, may reduce the reactive surface area and thus retard dissolution rates.The existence of a chemical lysocline in the Northwest Atlantic was confirmed. However, it seems that different kinds of particles have different lysoclines. The origin of the lysocline cannot be attributed to water flow or to the thermodynamic transition from supersaturation to undersaturation (Ω = 1). It seems to be a kinetic phenomenon.A simple model comparing the complete dissolution time and the residence time of a particle on the sediment-water interface suggests that coccoliths can be preserved in the sediments of the deep Northwest Atlantic below the CCD, in good agreement with SEM observations.     

The impact of freshwater discharges on the ocean circulation in the Skagerrak/northern North Sea area Part I: model validation

We assess the performance of an eddy-recognizing numerical ocean model in simulating the pattern and variability of the hydrography in the Skagerrak/northern North Sea area. The model we use is a version of the widely used Princeton ocean model employing a terrain-following vertical coordinate. Results from a series of five multi-year simulations of the mesoscale response are described. The simulations differ in their representation of the lateral freshwater supply to the model ocean of which the first is a reference simulation. The next four are variations in which the river discharges and/or the Baltic outflow are given more realistic representations. For validation, we have used in situ hydrographic data. A novelty is that we use the concepts of freshwater height and potential energy anomaly as objective validation tools. We find that, in general, the model faithfully reproduces many of the observed hydrographic features including their mean patterns and their variance. Not surprisingly, we find that the Baltic outflow is by far the most significant freshwater source in terms of its influence on the hydrography in the area, a result corroborating earlier findings. The best validation is obtained when all freshwater supply is made as realistic as possible, in particular the Baltic outflow. We also find that the large scale cyclonic circulation and the location of fronts are robust characteristics of the Skagerrak/northern North Sea circulation given the impact changes in the freshwater input has on the hydrography. Finally, we find that a further exploration of the impact of the lateral open boundary forcing, e.g., the input of Atlantic water, is needed.     

A numerical study of the string function using a primitive equation ocean model

Abstract

We use results from a primitive-equation ocean numerical model (SCRUM) to test a theoretical 'string function' formulation put forward by Tyler and Käse in another article in this issue. The string function acts as a stream function for the large-scale potential energy flow under the combined beta and topographic effects. The model results verify that large-scale anomalies propagate along the string function contours with a speed correctly given by the cross-string gradient. For anomalies having a scale similar to the Rossby radius, material rates of change in the layer mass following the string velocity are balanced by material rates of change in relative vorticity following the flow velocity. It is shown that large-amplitude anomalies can be generated when wind stress is resonant with the string function configuration.