Instability of planetary waves in a high vertical resolution model

Abstract

A high vertical resolution model is used to examine the instability of a baroclinic zonal flow and a finite amplitude topographically forced wave. Two families of unstable modes are found, consisting of zonally propagating most unstable modes, and stationary unstable modes. The former have time scale and spatial structure similar to baroclinic synoptic disturbances, but are localized in space due to interaction with the zonally asymmetric forcing. These modes transport heat efficiently in both the zonal and meridional directions. The second family of stationary unstable modes has characteristics of modes of low frequency variability of the atmosphere. They have time scales of 10 days and longer, and are of planetary scale with an equivalent barotropic vertical structure. The horizontal structure resembles blocking flows. They are maintained by available potential energy of the basic wave, and have large zonal heat fluxes. The results for both families of modes are interpreted in terms of an interaction between forcing and baroclinic instability to create favoured regions for eddy development. Applications to baroclinic planetary waves are also considered.     

Barotropic and baroclinic processes in the transport variability of the Antarctic Circumpolar Current

Synoptic scale variability of the Southern Ocean wind field in the high-frequency range of barotropic Rossby waves results in transport variations of the Antarctic Circumpolar Current (ACC), which are highly coherent with the bottom pressure field all around the Antarctic continent. The coherence pattern, in contrast to the steady state ACC, is steered by the geostrophic f/h contours passing through Drake Passage and circling closely around the continent. At lower frequencies, with interannual and decadal periods, the correlation with the bottom pressure continues, but baroclinic processes gain importance. For periods exceeding a few years, variations of the ACC transport are in geostrophic balance with the pressure field associated with the baroclinic potential energy stored in the stratification, whereas bottom pressure plays a minor role. The low-frequency variability of the ACC transport is correlated with the baroclinic state variable in the entire Southern Ocean, mediated by baroclinic topographic–planetary Rossby waves that are not bound to f/h contours. To clarify the processes of wave dynamics and pattern correlation, we apply a circulation model with simplified physics (the barotropic–baroclinic-interaction model BARBI) and use two types of wind forcing: the National Centers for Environmental Prediction (NCEP) wind field with integrations spanning three decades and an artificial wind field constructed from the first three empirical orthogonal functions of NCEP combined with a temporal variability according to an autoregressive process. Experiments with this Southern Annular Mode type forcing have been performed for 1,800 years. We analyze the spin-up, trends, and variability of the model runs. Particular emphasis is placed on coherence and correlation patterns between the ACC transport, the wind forcing, the bottom pressure field and the pressure associated with the baroclinic potential energy. A stochastic dynamical model is developed that describes the dominant barotropic and baroclinic processes and represents the spectral properties for a wide range of frequencies, from monthly periods to hundreds of years.     

Six circumpolar currents—on the forcing of the Antarctic Circumpolar Current by wind and mixing

The transport of the Antarctic Circumpolar Current (ACC) is influenced by a variety of processes and parameters. A proper implementation of basin geometry, ocean topography and baroclinicity is known to be a fundamental requisite for a realistic simulation of the circulation and transport. Other, more subtle parameters are those of eddy-induced transports and diapycnal mixing of thermohaline tracers or buoyancy, either treated by eddy resolution or by a proper parameterization. Quite a number of realistic numerical simulations of the circulation in the Southern Ocean have recently been published. Many concepts on relations of the ACC transport to model parameters and forcing function are in discussion, however, without much generality and little success. We present a series of numerical simulations of circumpolar flow with a simplified numerical model, ranging from flat-bottom wind-driven flow to baroclinic flow with realistic topography and wind and buoyancy forcing. Analysis of the balances of momentum, vorticity, and baroclinic potential energy enables us to develop a new transport theory, which combines the most important mechanisms driving the circulation of the ACC and determining its zonal transport. The theory is based on the importance of the bottom vertical velocity in generating vorticity and shaping the baroclinic potential energy of the ACC. It explains the breaking of the -constraint by baroclinicity and brings together in one equation the wind and buoyancy forcing of the current. The theory emphasizes the role of Ekman pumping and eddy diffusion of buoyancy to determine the transport. It also demonstrates that eddy viscosity effects are irrelevant in the barotropic vorticity balance and that friction arises via eddy diffusion of density. In this regime, the classical Stommel model of vorticity balance is revived where the bottom friction coefficient is replaced by (with the Gent–McWilliams coefficient and the baroclinic Rossby radius ) and a modified wind curl forcing appears.     

Wave-induced topographic formstress in baroclinic channel flow

Large-scale zonal flow driven across submarine topography establishes standing Rossby waves. In the presence of stratification, the wave pattern can be represented by barotropic and baroclinic Rossby waves of mixed planetary topographic nature, which are locked to the topography. In the balance of momentum, the wave pattern manifests itself as topographic formstress. This wave-induced formstress has the net effect of braking the flow and reducing the zonal transport. Locally, it may lead to acceleration, and the parts induced by the barotropic and baroclinic waves may have opposing effects. This flow regime occurs in the circumpolar flow around Antarctica. The different roles that the wave-induced formstress plays in homogeneous and stratified flows through a zonal channel are analyzed with the BARBI (BARotropic-Baroclinic-Interaction ocean model, Olbers and Eden, J Phys Oceanogr 33:2719–2737, 2003) model. It is used in complete form and in a low-order version to clarify the different regimes. It is shown that the barotropic formstress arises by topographic locking due to viscous friction and the baroclinic one due to eddy-induced density advection. For the sinusoidal topography used in this study, the transport obeys a law in which friction and wave-induced formstress act as additive resistances, and windstress, the effect of Ekman pumping on the density stratification, and the buoyancy forcing (diapycnal mixing of the stratified water column) of the potential energy stored in the stratification act as additive forcing functions. The dependence of the resistance on the system parameters (lateral viscosity ε, lateral diffusivity κ of eddy density advection, Rossby radius λ, and topography height δ) as well as the dependence of transport on the forcing functions are determined. While the current intensity in a channel with homogeneous density decreases from the viscous flat bottom case in an inverse quadratic law ~δ ^–2 with increasing topography height and always depends on ε, a stratified system runs into a saturated state in which the transport becomes independent of δ and ε and is determined by the density diffusivity κ rather than the viscosity: κ/λ ² acts as a vertical eddy viscosity, and the transport is λ ²/κ times the applied forcing. Critical values for the topographic heights in these regimes are identified.     

On the origin of jets in the ocean

We show a mechanism whereby the jets result during the development of β-plumes (i.e., low-frequency Rossby waves that establish gyre circulations) in a model of ocean-basin circulation. The energy originates in baroclinic meanders of circulation at the eastern boundary of the ocean. Eddies are intimately related and occur as a result of the instability of this process. This mechanism does not rely on the existence of the small-scale turbulence to establish zonal flows. Zonal jets can then be amplified by eddies arranged in certain order in the flow. The underlying dynamics include the propagation of linear and nonlinear basin scale Rossby waves. The related barotropic theory for these waves is developed here. We demonstrate the radiative development of jets and β-plumes in a laboratory experiment using a rotating fluid with a paraboloidal free surface. The dynamical fields are measured by the laboratory analog of the satellite altimetry.     

斜压基本气流对东亚夏季风区气旋扰动低频发展影响的数值模拟研究

本文将夏季气候平均的基本气流分解为正压和斜压分量,使用一个线性斜压模式,研究了不同斜压基本气流对热带西北太平洋地区初始气旋性环流扰动低频发展演变的重要作用.其中,控制试验较好地模拟出初始气旋扰动向西北方向传播、在西北太平洋季风槽附近停滞增强、在东亚地区出现经向波列和在南海到海洋大陆地区形成西北—东南向波列等特征.改变斜压分量的敏感性试验结果表明,正压基流不能为西传的初始扰动供给足够的能量;海陆热力差异引起东亚地区的纬向温度梯度和北风垂直切变,是东亚太平洋型经向波列形成和维持的重要因素;当基本气流中的斜压纬向偏差部分线性增大时,扰动的能量会呈e指数迅速增强,提示在气候变化的背景下,基本气流微小的改变可能带来天气或季节内扰动强度的剧烈响应.     

Roles of forced and inertially unstable convection development in the onset process of Indian summer monsoon

The NCEP/NCAR R1 reanalysis data are employed to investigate the impact of forced and inertial instability in the lower troposphere over the Arabian Sea on the onset process of Indian summer monsoon(ISM),and to reveal the important role of zonal advection of zonal geostrophic momentum played in the forced unstable convection.Results show that during the ISM onset the zero absolute vorticity contour(??=0)shifts northward due to the strong cross-equatorial pressure gradient in the lower troposphere over southern Arabian Sea.Thus a region with negative absolute vorticity is generated near the equator in the northern hemisphere,manifesting the evident free inertial instability.When a southerly passes through this region,under the influence of friction a lower convergence that facilitates the convection flourishing at the lower latitudes appears to the north of zero absolute vorticity contour.However,owing to such a traditional inertial instability,the convection is confined near the equator which does not have direct influence on the ISM onset.On the contrary in the region to the north of the zero absolute vorticity contour and to the south of the low pressure center near the surface,although the atmosphere there is inertially stable,the lower westerly jet can develop and bring on the apparent zonal advection of zonal geostrophic momentum.Both theoretical study and diagnosing analysis present that such a zonal advection of geostrophic momentum is closely associated with the zonal asymmetric distribution of meridional land-sea thermal contrast,which induces a convergence center near and further north of the westerly jet in the lower troposphere over the southwestern coast of the Indian Peninsula,providing a favorable lower circulation for the ISM onset.It illustrates that the development of convection over the Arabian Sea in late spring and early summer is not only due to the frictional inertial instability but also strongly affected by the zonal asymmetric distribution of land-sea thermal contrast.Moreover,before the ISM onset due to the eastward development of the South Asian High(SAH)in the upper troposphere,high potential vorticity is transported to the region over the Arabian Sea.Then a local trumpet-shaped stream field is generated to cause the evident upper divergence-pumping effect which favors the ISM onset.When the upper divergence is vertically coupled with the lower convergence resulted from the aforementioned forced unstable convection development near the southwestern coast of Indian Peninsula,the atmospheric baroclinic unstable development is stimulated and the ISM onset is triggered.     

A wind-driven model of the ocean surface layer with wave radiation physics

Ocean Dynamics - Surface windstress transfers energy to the surface mixed layer of the ocean, and this energy partly radiates as internal gravity waves with near-inertial frequencies into the...     

The generation of baroclinic Rossby waves by stationary and and translating currents part 1. Stationary current forcing

Abstract

This paper investigates the generation of linear, baroclinic Rossby waves by an imposed current distribution, in a reduced gravity ocean, both with and without an eastern coast. A zonal current is impulsively applied and maintained along the northern edge of the domain of solution. Using Green's function techniques, analytical solutions are found, and these are evaluated for small times. Numerical solutions are obtained for larger times. The upper layer depth field consists of a transient response, due to the sudden application of the current. Maintenance of the current causes a response which is singular along the line of imposed non-zero h _y. The interior field decays with time (this is shown asymptotically). The parameters used are appropriate for the mid-latitude North Pacific, and the results are relevant to sudden transport changes in the North Pacific Current.     

Causes and underlying dynamic processes of the mid-winter suppression in the North Pacific storm track

Baroclinic wave activity in the North Pacific exhibit peaks in late fall and early spring, and a local minimum in midwinter, when by linear baroclinic instability theory it should attain its maximum. This counterintuitive phenomenon, or"midwinter suppression"(MWM) as called, is investigated with a functional analysis apparatus, multiscale window transform(MWT), and the MWT-based theory of canonical transfer and localized multi-scale energetics analysis, together with a feature tracking technique, using the data from the European Centre for Medium-Range Weather Forecasts ReAnalysis(ERA-40). It is found that the MWM results from a variety of different physical processes, including baroclinic canonical transfer, diabatic effect, energy flux divergence, and frictional dissipation. On one hand, baroclinic canonical transfer and diabatic effect achieve their respective maxima in late fall. More transient available potential energy is produced and then converted to transient kinetic energy, resulting in a stronger storm track in late fall than in midwinter. On the other hand, in early spring, although baroclinic instability and buoyancy conversion are weak, energy flux convergences are substantially strengthened, leading to a net energy inflow into the storm track. Meanwhile, frictional dissipation is greatly reduced in spring; as a result, less transient energy is dissipated in early spring than in midwinter. It is further found that the weakening of baroclinic canonical transfer in midwinter(compared to late fall) is due to the far distance between the storm and the jet stream(located at its southernmost point), which suppresses the interaction between them. Regarding the increase in energy flux convergence in early spring, it appears to originate from the increase(enhancement) in the number(strength) of storms from the upstream into the Pacific.     

Marine algae are ‘taught’ the basics of angular momentum

Advanced modelling studies and high-resolution observations have shown that flows related to instability of the mesoscale (~ 1–10 km scale) may provide both the fertilisation mechanism for nutrient-depleted (oligotrophic) surface waters and a subduction mechanism for the rapid export of phytoplankton biomass to the deep ocean. Here, a detailed multidisciplinary analysis of the data from an example high-resolution observational campaign is presented. The data provide direct observations of the transport of phytoplankton through baroclinic instability. Furthermore, the data confirm that this transport is constrained by the requirement to conserve angular momentum, expressed in a stratified water column as the conservation of potential vorticity. This constraint is clearly seen to produce long thin filaments of phytoplankton populations strained out along isopycnal vorticity annuli associated with mesoscale frontal instabilities.     

初值和海温强迫对延伸期可预报性时空分布的影响

利用全球谱模式T106L19和增长模繁殖法,分别在气候海温和预测海温强迫下,进行了动力延伸集合预报试验.基于方差分析思想,利用集合预报结果,定义和计算了初值影响指数、海温强迫影响指数、潜在可预报性指数以及波动活动指数.通过分析四个指数,揭示了初值和海温强迫对延伸期可预报性时空分布以及潜在可预报性的影响,并探讨了其影响机理.结果表明：初值影响指数分布具有地域和季节的差异,初值的影响在中高纬度地区大于热带地区;相同季节,海温强迫影响指数分布与初值影响指数分布相似;潜在可预报性指数呈带状分布,大值集中在热带地区,且在低纬度地区,高层的潜在可预报性大于低层;初值和海温强迫对延伸期可预报性时空分布的影响,依赖于大气环流形势,初值和海温强迫影响的显著区正是大气长波的活跃区和西风急流区,急流区的强风切变为长波活动提供了斜压不稳定能量,而长波的发展调控着初值和海温强迫的影响,这说明延伸期的可预报性具有明显的流依赖性,大气外强迫的作用也与大气内部的动力过程密切相关.     

Mesoscale eddies and baroclinic instability over the eastern Sakhalin shelf of the Sea of Okhotsk: a model-based analysis

Based on an eddy-permitting numerical model, the mesoscale variability in the East-Sakhalin Current is investigated during the winter-spring period. Analysis of necessary conditions for the development of baroclinic instability showed that the nearshore component of the East-Sakhalin Current is potentially baroclinic unstable in the first half-year. The simulated circulation uncovered a generation of anticyclonic eddies on the eastern Sakhalin shelf. It was established that a spatial scale of these eddies and the first baroclinic Rossby radius of deformation are values of the same order; a lifetime of these eddies varies from 4 to 6 weeks, given the Rossby number varies from 0.05 to 0.2. Analysis of the rate of eddy energy conversion on the eastern Sakhalin shelf showed that the generation of the revealed mesoscale eddies results from, mainly, baroclinic instability, whereas barotropic instability can be both favoring and preventing to the generation of these eddies.     

Impact of the surface temperature and vertical shear of zonal wind on the dynamics of a simple two-layer model of the atmosphere

The paper presents and analyzes, from the point of view of smooth dynamic systems theory, a two-layer baroclinic model of the troposphere in geostrophic approximation. The model describes airflow in β-channel within the tropospheric part of the main Hadley circulation cell. It enables to obtain, after application of the Galerkin method, a fairly simple low-parametric dynamic system describing the phenomena of non-linear interactions, bifurcations and blocking in the atmosphere. This enables to take into consideration such basic factors influencing the atmospheric dynamics like the heat exchange within the surface, orography, vertical variability of zonal wind and hydrostatic stability. Impact of zonal thermal variability of the surface and vertical shear of zonal wind in the troposphere on the orographic bifurcation was investigated and the oscillation character in the dynamic system after Hopf bifurcation of the second kind was analyzed. Additionally, the model dynamics was investigated in conditions including momentum forcing in the upper and lower parts of the troposphere and excluding orographic interaction, as well as in the conditions of thermal interaction between the troposphere and the surface for the vertical shear of zonal wind in both tropospheric layers. Impact of the mean zonal wind in the troposphere on the properties of model dynamics was assessed. It was proved that zonally varied surface temperature and layered mean zonal wind in the atmosphere are the parameters that have basic influence on the model dynamics. They cause numerous bifurcations and strongly influence the periods of oscillations of the model variables. They are often Hopf bifurcations of the second kind during which tropospheric states fairly distant from the ones before the bifurcations are generated. This significantly influences the model predictability.     

Cyclonic eddies and upper thermocline fine-scale structures in the Antarctic Circumpolar Current

Ocean Dynamics - Mesoscale eddies in the open ocean are mostly formed by baroclinic instability, in which the available potential energy from the large-scale slope of the isopycnals is converted...     

The transition from symmetric to baroclinic instability in the Eady model

Here, we explore the transition from symmetric instability to ageostrophic baroclinic instability in the Eady model; an idealised representation of a submesoscale mixed layer front. We revisit the linear stability problem considered by Stone (J Atmos Sci, 23, 390–400, (Stone 1966)), Stone (J Atmos Sci, 27, 721–726, (Stone 1970)), Stone (J Atmos Sci, 29, 419–426, (Stone 1972)) with a particular focus on three-dimensional ‘mixed modes’ (which are neither purely symmetric or baroclinic) and find that these modes can have growth rates within just a few percent of the corresponding two-dimensional growth rate maximum. In addition, we perform very high resolution numerical simulations allowing an exploration of the transition from symmetric to baroclinic instability. Three-dimensional mixed modes represent the largest contribution to the turbulent kinetic energy during the transition period between symmetric and baroclinic instability. In each simulation, we see the development of sharp fronts with associated high rms vertical velocities of up to 30 mm s^?1. Furthermore, we see significant transfer of energy to small scales, demonstrated by time-integrated mixing and energy dissipation by small-scale three-dimensional turbulence totalling about 30 % of the initial kinetic energy in all cases.     

Maintenance of the mean kinetic energy in the global ocean by the barotropic and baroclinic energy routes: the roles of JEBAR and Ekman dynamics

In order to determine the maintenance mechanisms of the currents of the global ocean, this study investigates the budget of the annual mean kinetic energy (KE) in a high-resolution (0.1° × 0.1°) semi-global ocean simulation. The analysis is based on a separation of the mean KE using the barotropic (i.e., depth-averaged) and baroclinic (the residual) components of velocity. The barotropic and baroclinic KEs dominate in higher and lower latitudes, respectively, with their global average being comparable to each other. The working rates of wind forcing on the barotropic and baroclinic circulations in the global ocean are 243 and 747 gigawatts, respectively. This study presents at least three new results for the budget of the barotropic KE. Firstly, an energy diagram is rederived to show that the work of the barotropic component of the horizontal pressure gradient (HPG) is connected to the work related to the joint effect of baroclinicity and bottom relief (JEBAR), and then to the budget of potential energy (PE). Secondly, the model analysis shows that the globally averaged work of the barotropic HPG (which is connected to the work related to JEBAR and then to the budget of the PE) is nearly zero. This indicates that the wind- and buoyancy-induced barotropic circulations in the global ocean are of the same strength with opposite sign. Thirdly, it is found that the work of the wind forcing on the barotropic component of the simulated Antarctic Circumpolar Current (ACC) is canceled by the combined effect, in equal measure, of the work of the barotropic HPG and the work of dissipative processes for mean KE. This result makes a significant contribution to the discussion on the depth-integrated momentum balance of the ACC. The barotropic KE is dissipated by the effects of bottom frictional stress, lateral frictional stress, and the Reynolds stress, of which more than half is attributed to an unexpectedly large contribution from biharmonic horizontal friction. Future studies should pay more attention to the role of biharmonic friction used in high-resolution numerical models.     

Baroclinic energetics and zonal plane distribution of monsoon disturbances

The adiabatic, quasi-geostraphic, 25-layer, numerical, linear model with Ekman boundary layer friction is utilised to perform the baroclinic stability analysis of the mean monsoon zonal wind profile. It is shown thec _i is a function of the resultant wavenumber alone. This relation is able to explain the effects of the lateral walls on the unstable waves.The energetics and zonal plane distribution of the short and long preferred viscous waves are computed. The upward motion of the short wave together with the warm (cold) core lies to the west of the surface trough position above (below) 850 mb. Further, it is shown that the main source of kinetic energy for the wave lies in the middle layer (850–700 mb) which is transported to the lower and upper layers. Computed is found to be in good agreement with observed values.     

Effects of wave-wave and wave-mean flow interactions on the evolution of a baroclinic wave

Abstract

The weakly nonlinear evolution of a free baroclinic wave in the presence of slightly supercritical, vertically sheared zonal flow and a forced stationary wave field that consists of a single zonal scale and an arbitrary number of meridional harmonics is examined within the context of the conventional two-layer model. The presence of the (planetary-scale) stationary wave introduces zonal variations in the supercriticality and is shown to alter the growth rate and asymptotic equilibrium of the (synoptic-scale) baroclinic wave via two distinct mechanisms: The first is due to the direct interaction of the stationary wave with the shorter synoptic wave (wave-wave mechanism), and the second is due to the interaction of the synoptic wave with that portion of the mean field that is corrected by the zonally rectified stationary wave fluxes (wave-mean mechanism). These mechanisms can oppose or augment each other depending on the amplitude and spatial structure of the stationary wave field. If the stationary wave field is confined primarily to the upper (lower) layer and consists of only the gravest cross-stream mode, conditions are favorable (unfavorable) for nonzero equilibrium of the free wave.

In addition to the time dependent heat flux generated by baroclinic growth of the free wave, its interaction with a stationary wave field consisting of two or more meridional harmonics generates time dependent heat fluxes that vary with period of the free wave. However, if the stationary wave field contains several meridional harmonics of sufficiently large amplitude, the free baroclinic wave is destroyed.     

Eddy energy sources and mesoscale eddies in the Sea of Okhotsk

Based on eddy-permitting ocean circulation model outputs, the mesoscale variability is studied in the Sea of Okhotsk. We confirmed that the simulated circulation reproduces the main features of the general circulation in the Sea of Okhotsk. In particular, it reproduced a complex structure of the East-Sakhalin current and the pronounced seasonal variability of this current. We established that the maximum of mean kinetic energy was associated with the East-Sakhalin Current. In order to uncover causes and mechanisms of the mesoscale variability, we studied the budget of eddy kinetic energy (EKE) in the Sea of Okhotsk. Spatial distribution of the EKE showed that intensive mesoscale variability occurs along the western boundary of the Sea of Okhotsk, where the East-Sakhalin Current extends. We revealed a pronounced seasonal variability of EKE with its maximum intensity in winter and its minimum intensity in summer. Analysis of EKE sources and rates of energy conversion revealed a leading role of time-varying (turbulent) wind stress in the generation of mesoscale variability along the western boundary of the Sea of Okhotsk in winter and spring. We established that a contribution of baroclinic instability predominates over that of barotropic instability in the generation of mesoscale variability along the western boundary of the Sea of Okhotsk. To demonstrate the mechanism of baroclinic instability, the simulated circulation was considered along the western boundary of the Sea of Okhotsk from January to April 2005. In April, the mesoscale anticyclonic eddies are observed along the western boundary of the Sea of Okhotsk. The role of the sea ice cover in the intensification of the mesoscale variability in the Sea of Okhotsk was discussed.