Simple applications of McWilliams' “A note on a consistent quasigeostrophic model in a multiply connected domain”

We have solved several linear quasigeostrophic problems in singly and multiply connected domains to illustrate the use and effects of correct boundary conditions which balance circulation changes with forcing along each boundary. The results with these boundary conditions differ markedly from those obtained by setting the streamfunction to zero on the walls either when the fluid is forced by a time-dependent wind stress curl or when the fluid oscillates in a free basin mode in a domain which contains interior islands.     

On the interactions between planetary geostrophy and mesoscale eddies

Multiscale asymptotics are used to derive three systems of equations connecting the planetary geostrophic (PG) equations for gyre-scale flow to a quasigeostrophic (QG) equation set for mesoscale eddies. Pedlosky (1984), following similar analysis, found eddy buoyancy fluxes to have only a small effect on the large-scale flow; however, numerical simulations disagree. While the impact of eddies is relatively small in most regions, in keeping with Pedlosky’s result, eddies have a significant effect on the mean flow in the vicinity of strong, narrow currents.First, the multiple-scales analysis of Pedlosky is reviewed and amplified. Novel results of this analysis include new multiple-scales models connecting large-scale PG equations to sets of QG eddy equations. However, only introducing anisotropic scaling of the large-scale coordinates allows us to derive a model with strong two-way coupling between the QG eddies and the PG mean flow. This finding reconciles the analysis with simulations, viz. that strong two-way coupling is observed in the vicinity of anisotropic features of the mean flow like boundary currents and jets. The relevant coupling terms are shown to be eddy buoyancy fluxes. Using the Gent-McWilliams parameterization to approximate these fluxes allows solution of the PG equations with closed tracer fluxes in a closed domain, which is not possible without mesoscale eddy (or other small-scale) effects. The boundary layer width is comparable to an eddy mixing length when the typical eddy velocity is taken to be the long Rossby wave phase speed, which is the same result found by Fox-Kemper and Ferrari (2009) in a reduced gravity layer.     

A description of numerical antarctic circumpolar currents

A sequence of numerical calculations has been made for the equilibrium balances of eddies and mean currents in open and partially blocked, periodic channels. The physical model employed is a two-layer, quasigeostrophic, wind-driven one, with important bottom friction and weak lateral friction. The resolved eddies provide the interior fluxes of momentum and potential vorticity which allow the mean state to be a balanced one. The set of calculations does not provide a parameter study as such, but does provide examples of the influences of alternative physical processes and geometrical constraints. These alternatives include the presence or absence of a partial barrier across the channel, the length of the channel, the addition of a transient component to the wind-driving, and the addition of a topographic sill across the channel gap. Particular attention is focused upon the steadily driven general circulation of a β-plane channel, because of the structural simplicity of the solution. The results may be broadly summarized as follows. The eddies are generated by a baroclinic instability of the mean flow. They act to intensity the upper layer mean jet and mean cross-jet potential vorticity gradient (through eddy horizontal Reynolds stress and relative vorticity flux divergence, respectively) and to transfer downwards mean zonal momentum, energy, and potential vorticity gradient (through eddy interfacial pressure drag, vertical pressure work, and vortex stretching flux divergence, respectively). In the case of a zonally uniform channel, the meridional heat flux is found not to conform closely to previously proposed parameterizations. The presence of a partial meridional barrier and a topographic obstacle are found to strongly influence the equilibrium solution, while neither a change in the basin length nor the presence of a transient wind component appear to importantly alter the solution.     

A coupled surface boundary-layer-quasigeostrophic model

A three-dimensional model of the mesoscale surface boundary layer of the open ocean is developed through scale analysis of the primitive equations with mixing included. A set of surface boundary-layer equations appropriate for a broad range of oceanic and atmospheric scales is thereby derived. The essential basis of the model is a coupling between quasigeostrophic dynamics away from the boundary layer and arbitrary mixing models within the mixed layer. The coupling consists of advection of the boundary layer by the horizontal and vertical components of the interior quasigeostrophic flow and forcing of the interior by the boundary layer in the form of divergence within the boundary layer which leads to vortex stretching/compression in the interior. The divergence is generalized for mesoscale wind-driven flows and includes nonlinear interaction between the directly wind-driven boundary-layer flow and the interior flow in the form of interior relative vorticity advection by the wind-driven flow. The nature of the equations leads us to apply a numerical algorithm to their solution. This algorithm is calibrated through application to idealized problems to determine the temporal and spatial grid requirements. The model is initialized with a realistic ocean flow having the properties of the Gulf Stream.     

Influence of bottom topography roughness on the jet and inertial recirculation of a mid-latitude gyre

Several numerical experiments are conducted to examine the influence of mesoscale, bottom topography roughness on the inertial circulation of a wind-driven, mid-latitude ocean gyre. The ocean model is based on the quasi-geostrophic formulation, and is eddy-resolving as it features high vertical and horizontal resolutions (six layers and a 10 km grid). An antisymmetrical double-gyre wind stress curl forces the baroclinic modes and generates a strong surface jet. In the case of a flat bottom, inertia and inverse energy cascade force the barotropic mode, and the resulting circulation features strong, barotropic, inertial gyres. The sea-floor roughness inhibits the inertial circulation in the deep layers; the barotropic component of the flow is then forced by eddy-topography interactions, and its energy concentrates at the scales of the topography. As a result, the baroclinicity of the flow is intesified: the barotropic mode is reduced with regard to the baroclinic modes, and the bottom flow (constrained by the mesoscale sea-floor roughness) is decoupled from the surface flow (forced by the gyre-scale wind). Rectified, mesoscale bottom circulation induces an interfacial form stress at the thermocline, which enhances horizontal shear instability and opposes the eastward penetration of the jet. The mean jet is consequently shortened, but the instantaneous jet remains very turbulent, with meanders of large meridional extent. The sea-floor roughness modifies the energy pathways, and the eddies have an even more important role in the establishment of the mean circulation: below the thermocline, rectification processes are dominant, and eddies transfer energy toward permanent mesoscale circulations strongly correlated with topography, whereas above the thermocline mean flow and eddy generation are influenced by the mean bottom circulation through interfacial stress. The topography modifies the vorticity of the barotropic and highest baroclinic modes. Vorticity accumulates at the small topographic scales, and the vorticity content of the highest modes, which is very weak in the flat-bottom case, increases significantly. Few changes occur in surface-intensified modes. In the deep layers of the model, the inverse correlation between relative vorticity and topography at small scales ensures the homogenization of the potential vorticity, which mainly retains the largest scales of the bottom flow and the scale of β.     

Impact of ocean model resolution on CCSM climate simulations

The current literature provides compelling evidence suggesting that an eddy-resolving (as opposed to eddy-permitting or eddy-parameterized) ocean component model will significantly impact the simulation of the large-scale climate, although this has not been fully tested to date in multi-decadal global coupled climate simulations. The purpose of this paper is to examine how resolved ocean fronts and eddies impact the simulation of large-scale climate. The model used for this study is the NCAR Community Climate System Model version 3.5 (CCSM3.5)—the forerunner to CCSM4. Two experiments are reported here. The control experiment is a 155-year present-day climate simulation using a 0.5° atmosphere component (zonal resolution 0.625 meridional resolution 0.5°; land surface component at the same resolution) coupled to ocean and sea-ice components with zonal resolution of 1.2° and meridional resolution varying from 0.27° at the equator to 0.54° in the mid-latitudes. The second simulation uses the same atmospheric and land-surface models coupled to eddy-resolving 0.1° ocean and sea-ice component models. The simulations are compared in terms of how the representation of smaller scale features in the time mean ocean circulation and ocean eddies impact the mean and variable climate. In terms of the global mean surface temperature, the enhanced ocean resolution leads to a ubiquitous surface warming with a global mean surface temperature increase of about 0.2?°C relative to the control. The warming is largest in the Arctic and regions of strong ocean fronts and ocean eddy activity (i.e., Southern Ocean, western boundary currents). The Arctic warming is associated with significant losses of sea-ice in the high-resolution simulation. The sea surface temperature gradients in the North Atlantic, in particular, are better resolved in the high-resolution model leading to significantly sharper temperature gradients and associated large-scale shifts in the rainfall. In the extra-tropics, the interannual temperature variability is increased with the resolved eddies, and a notable increases in the amplitude of the El Ni?o and the Southern Oscillation is also detected. Changes in global temperature anomaly teleconnections and local air-sea feedbacks are also documented and show large changes in ocean–atmosphere coupling. In particular, local air-sea feedbacks are significantly modified by the increased ocean resolution. In the high-resolution simulation in the extra-tropics there is compelling evidence of stronger forcing of the atmosphere by SST variability arising from ocean dynamics. This coupling is very weak or absent in the low-resolution model.     

轴对称加热条件下圆盘实验中波的转变过程的初步分析

对于一控制在中心降冷、边缘加热的旋转圆盘内的流体,增加或减少其温差,可引起流体中波数的转变。本文对四波向三波及三波向四波的转变过程作了较详细的分析,发现在四波向三波转变时,平均经圈环流和西风强度均发生迅速的变化。维持四波时,平均经圈环流为靠近热源处下沉和靠近冷源处上升的反环流。当转变过程发生时,经圈环流转变成正环流,转换完成后则恢复成反环流。在经圈环流变化的同时,西风强度也发生由弱变强而后再由强变弱的转变。三波向四波转变时,经圈环流及西风强度没有上述变化,只有强度的不同。 根据热量输送的计算结果,波数不同,它所产生的热量的涡动输送也不一样,三波时热量的涡动输送较强,四波时热量的涡动输送较弱。最后,我们联系热量输送的特点对上述结果进行了初步的讨论。     

On deep mean flow generation mechanisms and the abyssal circulation of numerical model gyres

Mesoscale resolution ocean general circulation model (EGCM) experiments have been carried out under a variety of different model physical assumptions, and the different model systems often produce very different deep mean flow fields. The flat bottom, rectangular basin experiments exhibit two distinct types of deep mean flow, which are here called “corotating” and “counterrotating”. Counterrotating deep flow, in which two adjacent deep gyres, with circulation of opposite senses, underlie the upper ocean eastward jet and its recirculation, has been found only in models with adiabetic two-layer model physics. None of the more complex model systems exhibit counterrotating deep flows; this type of flow is apparently restricted to a particular range of forcing/dissipation parameter space and/or particular model physical assumptions.Since the deep flow in these EGCM systems is generally weak, geostrophic dynamics provides the basic deep flow interior balance and the mean vertical velocity field, through the lower layer vorticity equation, largely determines the deep interior flow. The dynamical constraints on the mean vertical velocity field introduced by different model physical equations are reviewed and the adiabatic quasi-geostrophic (QG) two-layer model system is shown to be strongly constrained in several respects. In particular, the idea that eddy and mean heat flux divergence (or “layer thickness flux divergence”) drive the mean vertical velocity does not generalize to more complicated dynamical systems in which there is the possibility of altering the mean vertical density profile and/or in which the horizontal flow can be divergent. As a consequence of the constraints, there can be no basin net vorticity input to the lower layer via vortex stretching in the QG system.Because of the adiabatic QG constraints and the particular parametric regime in which the published adiabatic QG EGCM experiments exist, a very plausible explanation can be found for the existence of the deep cyclonic circulation of the model subtropical gyre. It is this cyclonic circulation that causes these deep flows to differ so dramatically from those of the more physically complex model systems. Because all the published adiabatic QG experiments that have non-trivial deep flows exhibit the counterrotating behavior, and because available ocean data do not support the existence of such a gyre in the North Atlantic, it seems important to thoroughly understand the reasons for the existence or absence of the deep cyclonic circulations. If they are an invitable feature of adiabatic QG systems, these models may need to be treated with caution as tools for understanding the mean ocean circulation.     

黄海典型气旋与反气旋式海洋涡旋特征及影响机制模拟研究

为探讨黄海海洋涡旋的三维结构特征、能量输送与转换及影响机制,对黄海海域典型台风海洋气旋与近海海湾反气旋式涡旋个例进行数值模拟和时空诊断分析。采用FVCOM(Finite Volume Community Ocean Model)区域海洋数值模式精细化描述台风海洋涡旋与近海海洋中小尺度涡旋系统。对涡旋能量传输特征模拟显示,气旋式和反气旋式海洋涡旋中,非对称强流区动能能量下传比涡旋中心部位的强度更强,维持时间更长,下传深度更深。反气旋式海洋涡旋因Ekman流动形成的向中心辐合作用,造成此类差异更显著。气旋涡的动能主要来源于台风的近海面风应力动能和海洋涡旋有效位能的转换,反气旋涡旋区域风动力偏弱,其动能强度维持在低位,其涡旋增强伴随着有效位能的增加。环境因子影响机制从风浪,底摩擦和地形三方面讨论。结果显示:耦合波浪模块后,台风强风应力和风浪的综合作用扩大台风海洋涡旋尺度,并增强涡旋环流强度,同时对相邻的反气旋涡有压缩和减弱作用。风浪效应对台风海洋涡旋有正贡献。强台风过程表层环流响应台风应力而浅水地形和底摩擦强烈影响涡旋下层,造成台风海洋涡旋结构在垂直方向上偏移,并影响到下层环流速度减小,流向与表层相反。在海洋气旋涡和反气旋涡的显著辐散区,其混合层下方有温盐要素的涌升对应,辐合区有温盐要素的下沉对应;同时海底地形的升降也造成温盐强迫上升与下降,其强度与地形起伏尺度成正比,较环流系统作用更强。     

Poleward heat flux by an ocean gyre

A calculation to determine the poleward heat flux by an oceanic subtropical gyre is described. The circulation is given by Stommel's (1948) barotropic wind-driven model, and the temperature field is treated as a passive tracer. Air-sea heat exchange is parameterized as linearly proportional to the difference between air and sea temperatures. The problem so formulated reduces to the solution of the advection-diffusion (heat) equation for the temperature. The results show how ocean circulation, with a narrow western boundary layer, transports significantly more heat than does a symmetric circulation. A more intense circulation transports more heat but the increase levels off beyond a certain point. A wide ocean basin transports more of its heat by conduction than does a narrow basin. Approximate analytical solutions for both weak and strong circulations yield results that agree semiquantitatively with the numerical results.     

Energy analysis of open regions of turbulent flows — mean eddy energetics of a numerical ocean circulation experiment

The concepts involved in the interpretation of energy budgets in subregions of a turbulent flow are examined in order to determine the processes responsible for the production, transport, and dissipation of energy throughout a dynamically inhomogeneous circulation. An interpretation of the effects of Reynolds stress—mean flow interaction work for open regions is presented in terms of the change in the total mean kinetic energy. In an arbitrary volume of fluid the changes in kinetic energy of the mean flow and the mean kinetic energy of the eddy flow are not generally equal and opposite, so this process is not generally responsible for a conversion of energy between the two forms. These ideas are then applied to a regional kinetic energy analysis of the mesoscale resolution general ocean circulation numerical experiment of Robinson et al. (1977). The spatial structure of the various terms in the equation for the mean eddy kinetic energy is examined. The issues involved in selection of a set of analysis regions are discussed and explored via examination of budgets over different subregions of this flow. Thereby a relatively simple picture of the regional energetics emerges. Mean eddy kinetic energy is produced by conversion of kinetic energy of the mean flow in the net over the recirculation and near field of the northern boundary current system and roughly half of this energy is lost to each of mean eddy pressure work transport and diffusion work. Budgets over subregions of this net source region are much more complex. The interior eddy field is driven by pressure work influx, while the southwestern region has eddy buoyancy work conversion of mean potential energy as its energy source. At every depth level the eddy field draws its kinetic energy from the mean flow, when averaged over the horizontal extent of the basin or over the recirculation and near field.     

A high resolution numerical study of Gulf of Mexico fronts and eddies

Summary The Gulf of Mexico (GOM) circulation is simulated using the DieCAST ocean model, with a horizontal resolution of 1/12° and 20 vertical layers. The results compare well with observations of both large and small scale features, including Loop Current frontal occlusions associated with frontal eddies. The simulation is carried out without any data assimilation. The frontal eddies tend to be spaced at about 90° intervals around the Loop Current, leading to a Loop Current head shaped like a square with rounded corners. The pattern rotates as the eddies circle the Loop, and frontal eddies elongate as they squeeze through the Florida Strait. Major warm core eddies separate regularly from the Loop Current and propagate to the western GOM. Old warm core eddies in the western Gulf dissipate through bottom drag effects, which also generate cyclonic parasitic eddies. Newly arrived warm core eddies merge with old ones in the western GOM. Recently separated elongated Loop Current eddies can rotate and reattach temporarily to the Loop Current. The barotropic flow component develops eddies between the main separated warm core eddy and the Loop Current due to eastward dispersion, as the main eddy itself propagates westward into the Gulf.With 10 Figures     

Transports and budgets of volume,heat, and salt from a global eddy-resolving ocean model

The results from an integration of a global ocean circulation model have been condensed into an analysis of the volume, heat, and salt transports among the major ocean basins. Transports are also broken down between the model's Ekman, thermocline, and deep layers. Overall, the model does well. Horizontal exchanges of mass, heat, and salt between ocean basins have reasonable values; and the volume of North Atlantic Deep Water (NADW) transport is in general agreement with what limited observations exist. On a global basis the zonally integrated meridional heat transport is poleward at all latitudes except for the latitude band 30°S to 45°S. This anomalous transport is most likely a signature of the model's inability to form Antarctic Intermediate (AAIW) and Antarctic bottom water (AABW) properly. Eddy heat transport is strong at the equator where its convergence heats the equatorial Pacific about twice as much as it heats the equatorial Atlantic. The greater heating in the Pacific suggests that mesoscale eddies may be a vital mechanism for warming and maintaining an upwelling portion of the global conveyor-belt circulation. The model's fresh water transport compares well with observations. However, in the Atlantic there is an excessive southward transport of fresh water due to the absence of the Mediterranean outflow and weak northward flow of AAIW. Eddies in the mid-latitudes act to redistribute heat and salt down the mean gradients. Residual fluxes calculated from a sum of the computed advective (including eddies), forced, and stored fluxes of heat and salt represent transport mostly due to vertical sub-grid scale mixing processes. Perhaps the model's greatest weakness is the lack of strong AAIW and AABW circulation cells. Accurate thermohaline forcing in the North Atlantic (based on numerous hydrographic observations) helps the model adequately produce NADW. In contrast, the southern ocean is an area of sparse observation. Better thermohaline observations in this area may be needed if models such as this are to produce the deep convection that will achieve more accurate simulations of the global 3-dimensional circulation.     

Eddy fluxes and sensitivity of the water cycle to spatial resolution in idealized regional aquaplanet model simulations

A multi-scale moisture budget analysis is used to identify the mechanisms responsible for the sensitivity of the water cycle to spatial resolution using idealized regional aquaplanet simulations. In the higher resolution simulations, moisture transport by eddy fluxes dry the boundary layer enhancing evaporation and precipitation. This effect of eddies, which is underestimated by the physics parameterizations in the low-resolution simulations, is found to be responsible for the sensitivity of the water cycle both directly, and through its upscale effects on the transport of mean moisture by the mean circulation. Correlations among moisture transport by eddies at adjacent ranges of scales provides a potential for reducing this sensitivity by representing the unresolved eddies by their marginally resolved counterparts.     

盆地环流型及冷湖的形成和消散的数值研究

本文建立了一个三维原始方程数值模式，模拟三维盆地中的环流，重点讨论在热力强迫环流型的转变过程中能量场的动量场的配置情况及环流演变过程。模拟结果表明，夜间地面强烈的辐射冷却过程引起盆地各坡面冷迳流的发展，盆地中成为冷空气湖，逆温层结建立。日出后对流边界层从大气低层向上发展，逆温从底部随之消蚀，仅在上层存在了夜间逆温残余，即所谓稳定核，上坡风在对流边界层发展，而稳定核内仍维持夜间下坡风环流。     

Spatial and temporal variability of major ocean currents and mesoscale eddies

A global picture of ocean current variability may be obtained by analyzing surface drift currents in terms of their mean and eddy kinetic energies. High values for both quantities are found in the western boundary currents and in the equatorial current system; low values are found in the interior of major gyres. However, nowhere are eddy energies less than 200 cm² s^–2, indicating that, even in the least energetic parts of the oceans, surface speeds of ~20 cm s^–1 prevail. Recent experimental studies also support the widespread occurrence of mesoscale mid-oceanic eddies. Another type of eddy is abundant in the vicinity of boundary currents: examples include Gulf Stream Rings, the Great Whirl of the Somali Current, and disturbances of the predominantly zonal equatorial flow manifested by large-scale meandering about the equator. Recent numerical models using low-viscosity and high-resolution computational grids also reveal the ubiquitous existence of mesoscale structures. The importance of eddies is that they seem to be energetic enough and sufficiently widespread so as to play some part — not yet understood — in the circulation of the world ocean. Speculative analogies to the atmosphere suggest that the mesoscale ocean eddies are the storms and weather systems of the sea. We need global statistics on their distribution, their occurrence in various oceanic regions, their dimensions, and their lifetimes. The prospect of even a single global oceanic weather map, comparable to those obtained daily for the atmosphere, is hopeless in terms ofin situ oceanographic observations. Remote sensing may provide a partial solution.In the past, sea-surface temperature observations by satellites have revealed cold and warm eddies shed by western boundary currents. Satellite observations, moreover, have shown the thermal effects of continental shelf waves and areas with pronounced upwelling. These phenomena are characterized by strong temperature gradients and relatively large differences in surface elevation. Mid-ocean eddies are far more subtle and difficult to observe since they possess smaller differences of temperature and sea-surface topography. Three representative examples, including recent experimental results in western boundary currents, the equatorial region, and a typical mid-ocean region, are discussed in detail. Typical signals of temperature and sea-level topography, as well as typical temporal and spatial scales of the observed phenomena, are given.     

Evolving Perspectives on Abrupt Seasonal Changes of the General Circulation

Professor Duzheng YE(Tu-cheng YEH) was decades ahead of his time in proposing a model experiment to investigate whether abrupt seasonal changes of the general circulation can arise through circulation feedbacks alone, unrelated to underlying inhomogeneities at the lower boundary. Here, we introduce Professor YEH's ideas during the 1950 s and 1960 s on the general circulation and summarize the results and suggestions of Yeh et al.(1959) on abrupt seasonal changes. We then review recent advances in understanding abrupt seasonal changes arising from model experiments like those proposed by Yeh et al.(1959). The model experiments show that circulation feedbacks can indeed give rise to abrupt seasonal transitions.In these transitions, large-scale eddies that originate in midlatitudes and interact with the zonal mean flow and meridional overturning circulations in the tropics play central roles.     

Modeled Sverdrup flow in the North Atlantic from 11 different wind stress climatologies

In studies of large-scale ocean dynamics, often quoted values of Sverdrup transport are computed using the Hellerman–Rosenstein wind stress climatology. The Sverdrup solution varies, however, depending on the wind set used. We examine the differences in the large-scale upper ocean response to different surface momentum forcing fields for the North Atlantic Ocean by comparing the different Sverdrup interior/Munk western boundary layer solutions produced by a 1/16° linear numerical ocean model forced by 11 different wind stress climatologies. Significant differences in the results underscore the importance of careful selection of a wind set for Sverdrup transport calculation and for driving nonlinear models. This high-resolution modeling approach to solving the linear wind-driven ocean circulation problem is a convenient way to discern details of the Sverdrup flow and Munk western boundary layers in areas of complicated geometry such as the Caribbean and Bahamas. In addition, the linear solutions from a large number of wind sets provide a well-understood baseline oceanic response to wind stress forcing and thus, (1) insight into the dynamics of observed circulation features, by themselves and in conjunction with nonlinear models, and (2) insight into nonlinear model sensitivity to the choice of wind-forcing product.The wind stress products are evaluated and insight into the linear dynamics of specific ocean features is obtained by examining wind stress curl patterns in relation to the corresponding high-resolution linear solutions in conjunction with observational knowledge of the ocean circulation. In the Sverdrup/Munk solutions, the Gulf Stream pathway consists of two branches. One separates from the coast at the observed separation point, but penetrates due east in an unrealistic manner. The other, which overshoots the separation point at Cape Hatteras and continues to flow northward along the continental boundary, is required to balance the Sverdrup interior transport. A similar depiction of the Gulf Stream is commonly seen in the mean flow of nonlinear, eddy-resolving basin-scale models of the North Atlantic Ocean. An O(1) change from linear dynamics is required for realistic simulation of the Gulf Stream pathway. Nine of the eleven Sverdrup solutions have a C-shaped subtropical gyre, similar to what is seen in dynamic height contours derived from observations. Three mechanisms are identified that can contribute to this pattern in the Sverdrup transport contours. Along 27°N, several wind sets drive realistic total western boundary current transport (within 10% of observed) when a 14 Sv global thermohaline contribution is added (COADS, ECMWF 10 m re-analysis and operational, Hellerman–Rosenstein and National Centers for Environmental Prediction (NCEP) surface stress re-analysis), a few drive transport that is substantially too high (ECMWF 1000 mb re-analysis and operational and Isemer–Hasse) and Fleet Numerical Meteorology and Oceanography Center (FNMOC) surface stresses give linear transport that is slightly weaker than observed. However, higher order dynamics are required to explain the partitioning of this transport between the Florida Straits and just east of the Bahamas (minimal in the linear solutions vs. 5 Sv observed east of the Bahamas). Part of the Azores Current transport is explained by Sverdrup dynamics. So are the basic path of the North Atlantic Current (NAC) and the circulation features within the Intra-Americas Sea (IAS), when a linear rendition of the northward upper ocean return flow of the global thermohaline circulation is added in the form of a Munk western boundary layer.     

An intermediate model for large-scale ocean circulation studies

The design and purposes of an intermediate model are discussed along with fundamentals of the model and results of numerical experiments. The main purposes of the model are reconstructions of the schemes of the ocean large-scale circulation and paleocirculation. For these problems numerical effectiveness is the key factor. A novel feature is the parameterization of the side wall Ekman boundary layers which was introduced to enable the use of geostrophy for calculating baroclinic velocity. This approach of connecting the side frictional layers with a non-viscous interior is not model-specific and can be used in any model employing geostrophy in the interior. The method can facilitate the no-flux and no-slip boundary conditions at the side walls in such models. Preliminary numerical experiments with simple basin geometry and idealized forcing aimed at a comparison with primitive equation and planetary geostrophic models are carried out. A direct comparison with the Geophysical Fluid Dynamics Laboratory (GFDL) primitive equation model was performed for a quantitative test of the proposed model. The results of the experiments are discussed in the context of the applicability of intermediate models for studying ocean climate dynamics.     

Impact of 1/8° to 1/64° resolution on Gulf Stream model–data comparisons in basin-scale subtropical Atlantic Ocean models

We investigate the impact of 1/8°, 1/16°, 1/32°, and 1/64° ocean model resolution on model–data comparisons for the Gulf Stream system mainly between the Florida Straits and the Grand Banks. This includes mean flow and variability, the Gulf Stream pathway, the associated nonlinear recirculation gyres, the large-scale C-shape of the subtropical gyre and the abyssal circulation. A nonlinear isopycnal, free surface model covering the Atlantic from 9°N to 47°N or 51°N, including the Caribbean and Gulf of Mexico, and a similar 1/16° global model are used. The models are forced by winds and by a global thermohaline component via ports in the model boundaries. When calculated using realistic wind forcing and Atlantic model boundaries, linear simulations with Munk western boundary layers and a Sverdrup interior show two unrealistic mean Gulf Stream pathways between Cape Hatteras and the Grand Banks, one proceeding due east from Cape Hatteras and a second one continuing northward along the western boundary until forced eastward by the regional northern boundary. The northern pathway is augmented when a linear version of the upper ocean global thermohaline contribution to the Gulf Stream is added as a Munk western boundary layer. A major change is required to obtain a realistic pathway in nonlinear models. Resolution of 1/8° is eddy-resolving but mainly gives a wiggly version of the linear model Gulf Stream pathway and weak abyssal flows except for the deep western boundary current (DWBC) forced by ports in the model boundaries. All of the higher resolution simulations show major improvement over the linear and 1/8° nonlinear simulations. Additional major improvement is seen with the increase from 1/16° to 1/32° resolution and modest improvement with a further increase to 1/64°. The improvements include (1) realistic separation of the Gulf Stream from the coast at Cape Hatteras and a realistic Gulf Stream pathway between Cape Hatteras and the Grand Banks based on comparisons with Gulf Stream pathways from satellite IR and from GEOSAT and TOPEX/Poseidon altimetry (but 1/32° resolution was required for robust results), (2) realistic eastern and western nonlinear recirculation gyres (which contribute to the large-scale C-shape of the subtropical gyre) based on comparisons with mean surface dynamic height from the generalized digital environmental model (GDEM) oceanic climatology and from the pattern and amplitude of sea surface height (SSH) variability surrounding the eastern gyre as seen in TOPEX/Poseidon altimetry, (3) realistic upper ocean and DWBC transports based on several types of measurements, (4) patterns and amplitude of SSH variability which are generally realistic compared to TOPEX/Poseidon altimetry, but which vary from simulation to simulation for specific features and which are most realistic overall in the 1/64° simulation, (5) a basin wide explosion in the number and strength of mesoscale eddies (with warm core rings (WCRs) north of the Gulf Stream, the regional eddy features best observed by satellite IR), (6) realistic statistics for WCRs north of the Gulf Stream based on comparison to IR analyses (low at 1/16° resolution and most realistic at 1/64° resolution for mean population and rings generated/year; realistic ring diameters at all resolutions), and (7) realistic patterns and amplitude of abyssal eddy kinetic energy (EKE) in comparison to historical measurements from current meters.