Barotropic unstable modes in zonal and meridional channel on the beta-plane

Abstract

The linear stability of a non-divergent barotropic parallel shear flow in a zonal and a non-zonal channel on the β plane was examined numerically. When the channel is non-zonal, the governing equation is slightly modified from the Orr-Sommerfeld equation. Numerical solutions were obtained by solving the discretized linear perturbation equation as an eigenvalue problem of a matrix. When the channel is zonal and lateral viscosity is neglected the problem is reduced to the ordinary barotropic instability problem described by Kuo's (1949) equation. The discrepancy between the stability properties of westward and eastward flows, which have been indicated by earlier studies, was reconfirmed. It has also been suggested that the unstable modes are closely related to the continuous modes discretized by a finite differential approximation. When the channel is non-zonal, the properties of unstable modes were quite different from those of the zonal problem in that: (1) The phase speed of the unstable modes can exceed the maximum value of the basic flow speed; (2) The unstable modes are not accompanied by their conjugate mode; and (3) The basic flow without an inflection point can be unstable. The dispersion relation and the spatial structure of the unstable modes suggested that, irrespective of the orientation of the channel, they have close relation to the neutral modes (Rossby channel modes) which are the solutions in the absence of a basic shear flow. The features mentioned above are not dependent on whether or not the flow velocity at the boundary is zero.     

Three‐dimensional non‐geostrophic disturbances in a baroclinic zonal flow

Abstract

A study is made to determine the stability properties of a baroclinic zonal current on which small amplitude three‐dimensional non‐geostrophic disturbances are superimposed. The flow is assumed to be bounded to the north and south by rigid vertical walls and the Rossby number Ro is taken to be small compared to unity. It is then shown that if the perturbation quantities are expanded in power series in Ro the leading or zero order terms in the series correspond to the quasi‐geostrophic solution obtained by Eady (1949) and that the higher order terms represent the “non‐geostrophic” effects neglected by the latter.

It is shown that to the second order in Ro the non‐geostrophic effects decrease the growth rates of those disturbances which are found to be unstable according to Eady's analysis but do not alter their speed of propagation. The results indicate, on the other hand, that to the same order of approximation the stable waves travel at a speed which is different from that given by Eady's solution. The modification of the perturbation wave structure by the non‐geostrophic effects is also investigated. It is found in particular that to the first order in Ro the latter produce a northward tilt with height in the ridge (or trough) lines of the meridional and vertical particle velocity fields away from the lateral boundaries.     

The nonlinear development of disturbances in a zonal shear flow

Abstract

A study is made of the nonlinear stability of a weakly supercritical zonal shear flow in the β-plane approximation. The dynamics of initially small disturbances are examined. The main nonlinear effects are associated with the rearrangement of the critical layer. It is shown that as the wave grows in amplitude, linear regimes of the critical layer (viscous and nonstationary) change over to a nonlinear regime while the exponential law of disturbance growth becomes a power-law.     

Inertial waves in a fluid partially filling a cylindrical cavity during spin-up from rest

Abstract

Inertial waves are excited in a fluid contained in a slightly tilted rotating cylindrical cavity while the fluid is spinning up from rest. The surface of the fluid is free. Since the perturbation frequency is equal to the rotation speed resonance occurs at a critical height to radius aspect ratio of the fluid. Detailed study of a particular inertial wave shows that in solid body rotation this “eigenratio” agrees with predictions from linear inviscid theory to within 0.5%. Measured time dependence of the eigenratio during spin-up from rest is a function of the tilt amplitude and agrees favorably with predictions from a numerical study. Mean flow associated with the inertial wave becomes unstable during spin-up and in the steady state. A boundary for the unstable region is found experimentally.     

The instability of barotropic circular vortices

Abstract

The linear, normal mode instability of barotropic circular vortices with zero circulation is examined in the f-plane quasigeostrophic equations. Equivalents of Rayleigh's and Fjortoft's criteria and the semicircle theorem for parallel shear flow are given, and the energy equation shows the instability to be barotropic. A new result is that the fastest growing perturbation is often an internal instability, having a finite vertical scale, but may also be an external instability, having no vertical structure. For parallel shear flow the fastest growing perturbation is always an external instability; this is Squire's theorem. Whether the fastest growing perturbation is internal or external depends upon the profile: for mean flow streamfunction profiles which monotonically decrease with radius, the instability is internal for less steep profiles with a broad velocity extremum and external for steep profiles with a narrow velocity extremum. Finite amplitude, numerical model calculations show that this linear instability analysis is not valid very far into the finite amplitude range, and that a barotropic vortex, whose fastest growing perturbation is internal, is vertically fragmented by the instability.     

Topographic instability and multiple equilibria on an f-plane

Abstract

The flow of a two-layer flow in a rotating channel on an f-plane over topography with sinusoidal variation of height in a direction parallel to the flow is investigated. When the two layers flow in opposite directions a resonance is found when the topographic scale matches the free mode of the system. We examine the stability of the forced mode in the vicinity of this resonance by means of a perturbation expansion of the topographic height. Both subresonant and super-resonant instabilities are found and their equilibration is examined. For small values of the dissipation multiple equilibria are found. The topographic drag releases potential energy even when the flow is baroclinically stable.     

Nonlocal modons on the beta-plane

Abstract

The vortex pair known as a modon is a classical solitary wave in the sense that it decays exponentially with distance from the center of the wave whenever the modon's phase speed of the wave is outside the linear range. In contrast, when ?1 < c < 0, the modon “far field” is oscillatory so that the modon is “nonlocal” in the sense that it has nonzero amplitude even at arbitrarily far distances from the vortex maximum. However, Tribbia and Verkley have independently noted that the oscillatory far field may be very weak for some parameter ranges.     

On the stability of plane flow with horizontal shear to three-Dimensional nondivergent disturbances

Abstract

Arnold's (1965a) method is used to investigate the stability of a stationary, nonparallel, plane flow, with horizontal shear, to three-dimensional nondivergent disturbances in a Boussinesq fluid. It is shown that, if the fluid is statically stable, the Rayleigh condition is not sufficient to insure inertial stability to all disturbance modes. For channel flow it is possible to establish the sufficiency condition for stability to some of these modes.     

Macrodynamics of α2 dynamos

Abstract

Two distributions of the α-effect in a sphere are considered. The inviscid limit is approached both by direct numerical solution and by solution of a simpler nonlinear eigenvalue problem deriving from asymptotic boundary layer analysis for the case of stress-free boundaries. The inviscid limit in both cases is dominated by the need to satisfy the Taylor constraint which states that the integral of the Lorentz force over cylindrical (geostrophic) contours in a homogeneous fluid must tend to zero. For a small supercritical range in α, this condition can only be met by magnetic fields which vanish as the viscosity goes to zero. In this range, the agreement of the two approaches is excellent. In a portion of this range, the method of finite amplitude perturbation expansion is useful, and serves as a guide for understanding the numerical results. For larger α, evidence from the nonlinear eigenvalue problem suggests both that the Taylor state exists, and that the transition from small to large amplitude can require a finite amplitude (oscillatory) instability in accord with the findings of Soward and Jones (1983). However, solutions of the full equations have not been found which are independent of viscosity at larger values of α.     

Experimental study of upstream influence in the two‐dimensional flow of a stratified fluid over an obstacle†

An experimental investigation was made of the upstream influence in front of two‐dimensional obstacles when they were towed in a linearly stratified fluid. The experiments were performed in a plexiglas channel 30.5 feet long, 2 feet high and 14 inches wide filled with a linearly stratified salt solution. Velocity measurements and flow visualization were obtained by neutrally buoyant liquid droplets and dye lines. Density measurements were made by a salinity probe.

The existence of unattenuated upstream influence in front of an obstacle was quantitatively documented for the first time. It occurred in the form of multiple unattenuated horizontal jets when there was a separated open wake behind the obstacle. These jets were identified to be the super‐position of “columnar disturbance modes”. The total number of columnar modes was determined solely by the Froude number of the flow and was equal to the number of lee‐wave modes excited. The drag due to upstream columnar modes was estimated and found to be lower than the drag due to the lee wave modes:     

The nonlinear stabilization of a zonal shear flow instability

Abstract

The development of initially small perturbations in a weakly supercritical zonal shear flow on a β-plane is studied. Two different scenarios of evolution are possible. If the supercriticality is sufficiently small, the growth of a perturbation is stopped in the viscous critical layer regime; for this case the evolution equation (corrected by the inclusion of a quintic nonlinearity) is derived. At greater supercriticality the nonlinearity cannot stop the growth of the perturbation in a linear (viscous or unsteady) critical layer regime, and the evolution is more complicated. Transition to a nonlinear critical layer regime leads to a reduction in the growth rate and to a slowing (but not a stopping) of the increase in amplitude, A. These are connected to the formation of a plateau (S=constant) of width L=O(A ^½) in the profile of absolute vorticity, S. Careful analysis reveals that the growth in amplitude ceases only when the whole instability domain (where the slope of unperturbed S-profile is positive) becomes covered again by the plateau.     

On steady and travelling waves over bottom topography in the model of homogeneous flow in a beta-plane channel

Abstract

A spectral low-order model is proposed in order to investigate some effects of bottom corrugation on the dynamics of forced and free Rossby waves. The analysis of the interaction between the waves and the topographic modes in the linear version of the model shows that the natural frequencies lie between the corresponding Rossby wave frequencies for a flat bottom and those applying in the “topographic limit” when the beta-effect is zero. There is a possibility of standing or eastward-travelling free waves when the integrated topograhic effect exceeds the planetary beta-effect.

The nonlinear interactions between forced waves in the presence of topography and the beta-effect give rise to a steady dynamical mode correlated to the topographic mode. The periodic solution that includes this steady wave is stable when the forcing field moves to the West with relatively large phase speed. The energy of this solution may be transferred to the steady zonal shear flow if the spatial scale of this zonal mode exceeds the scale of the directly forced large-scale dynamical mode.     

Investigating hydrological regimes and processes in a set of catchments with temporary waters in Mediterranean Europe

Abstract

Seven catchments of diverse size in Mediterranean Europe were investigated in order to understand the main aspects of their hydrological functioning. The methods included the analysis of daily and monthly precipitation, monthly potential evapotranspiration rates, flow duration curves, rainfall—runoff relationships and catchment internal data for the smaller and more instrumented catchments. The results showed that the catchments were less “dry” than initially considered. Only one of them was really semi-arid throughout the year. All the remaining catchments showed wet seasons when precipitation exceeded potential evapotrans-piration, allowing aquifer recharge, “wet” runoff generation mechanisms and relevant baseflow contribution. Nevertheless, local infiltration excess (Hortonian) overland flow was inferred during summer storms in some catchments and urban overland flow in some others. The roles of karstic groundwater, human disturbance and low winter temperatures were identified as having an important impact on the hydrological regime in some of the catchments.     

On the inadequacies of Long's model for steady two-dimensional stratified flows

Nonlinear analysis of two-dimensional steady flows with density stratification in the presence of gravity is considered. Inadequacies of Long's model for steady stratified flow over topography are explored. These include occurrence of closed streamline regions and waves propagating upstream. The usual requirements in Long's model of constant dynamic pressure and constant vertical density gradient in the upstream condition are believed to be the cause of these inadequacies. In this article, we consider a relaxation of these requirements, and also provide a systematic framework to accomplish this. As illustrations of this generalized formulation, exact solutions are given for the following two special flow configurations: the stratified flow over a barrier in an infinite channel; the stratified flow due to a line sink in an infinite channel. These solutions exhibit again closed-streamline regions as well as waves propagating upstream. The persistence of these inadequacies in the generalized Long's model appears to indicate that they are not quite consequences of the assumptions of constant dynamic pressure and constant vertical density gradient in Long's model, contrary to previous belief.

On the other hand, solutions admitted by the generalized Long's model show that departures from Long's model become small as the flow becomes more and more supercritical. They provide a nonlinear mechanism for the generation of columnar disturbances upstream of the obstacle and lead in subcritical flows to qualitatively different streamline topological patterns involving saddle points, which may describe the lee-wave-breaking process in subcritical flows and could serve as seats of turbulence in real flows. The occurrences of upstream disturbances in the presence of lee-wave-breaking activity described by the present solution are in accord with the experiments of Long (Long, R.R., “Some aspects of the flow of stratified fluids, Part 3. Continuous density gradients”, Tellus 7, 341--357 (1955)) and Davis (Davis, R.E., “The two-dimensional flow of a stratified fluid over an obstacle”, J. Fluid Mech. 36, 127–143 ()).     

The effects of vertical shear and stratification on stationary rossby waves

Abstract

The generation of stationary Rossby waves by sources of potential vorticity in a westerly flow is examined here in the context of a two-layer, quasi-geostrophic, β-plane model. The response in each layer consists of a combination of a barotropic Rossby wave disturbance that extends far downstream of the source, and a baroclinic disturbance which is evanescent or wave-like in character, depending on the shear and degree of stratification. Contributions from each of these modes in each layer are strongly dependent on the basic flows in each layer; the degree of stratification; and the depths of the two layers. The lower layer response is dominated by an evanescent baroclinic mode when the upper layer westerlies are much larger than those in the lower layer. In this case, weak stationary Rossby waves of large wavelengths are confined to the upper layer and the disturbance in the lower layer is confined to the source region.

Increasing the upper layer flow (with the lower layer flow fixed) increases the Rossby wavelength and decreases the amplitude. Decreasing the lower layer flow (with the upper layer flow fixed) decreases the wavelength and increases the amplitude. Stratification increases the contribution from the barotropic wave-like mode and causes the response to be confined to the lower layer.

The finite amplitude response to westerly flow over two sources of potential vorticity is also considered. In this case stationary Rossby waves induced by both sources interact to reinforce or diminish the downstream wave pattern depending on the separation distance of the sources relative to the Rossby wavelength. For fixed separation distance, enhancement of the downstreatm Rossby waves will only occur for a narrow range of flow variables and stratification.     

Quantification of the stability of river flow regimes

Abstract

An important characteristic of a river flow regime type is the time of year when high and low flows are likely to occur. How likely is it, however, to observe an identified seasonal pattern each individual year? Stability is an often neglected property of a flow regime, though shifts in the seasonal behaviour of flows affect both environmental and economic activities. An approach to characterize objectively the stability of a flow regime type, based on the concept of entropy, is presented. The stabilities of river flow maxima and minima are studied separately to investigate their respective contributions to the stability character of a particular regime type. A quantitative “instability index” permits a study of the development of a flow regime's stability in time, especially important in the context of a possible climate change. The method is presented using the example of a quantitative flow regime classification developed for Scandinavia and western Europe.     

Solitary Rossby Waves in Zonally Varying Jet Flows

The dynamics of solitary Rossby waves (SRWs) embedded in a meridionally sheared, zonally varying background flow are examined using a non-divergent barotropic model centered on a midlatitude g -plane. The zonally varying background flow, which is produced by an external potential vorticity (PV) forcing, yields a modified Korteweg-de Vries (K-dV) equation that governs the spatial-temporal evolution of a disturbance field that contains both Rossby wave packets and SRWs. The modified K-dV equation differs from the classical equation in that the zonally varying background flow, which varies on the same scale as the disturbance field, directly affects the disturbance linear translation speed and linear growth characteristics. In the limit of a locally parallel background flow, equations governing the amplitude and propagation characteristics of SRWs are derived analytically. These equations show, for example, that a sufficiently large (small) translation speed and/or a sufficiently weak (strong) background zonal shear favor transmission (reflection) of the SRW through (from) the jet. Conservation equations are derived showing that time changes in the domain averaged amplitude ("mass") or squared amplitude ("momentum") are due to zonal variation in both the linear, long-wave phase speed and linear growth; dispersion and nonlinearity do not affect the "mass" or "momentum". Provided (1) the background PV forcing is sufficiently small, or (2) the background PV forcing is meridionally symmetric and the disturbance is a SRW, the dynamics of the disturbance field is Hamiltonian and mass and energy are thus conserved. Numerical solutions of the K-dV equation show that the zonally varying background flow yields three general classes of behavior: reflection, transmission, or trapping. Within each class there exists SRWs and Rossby wave packets. SRWs that become trapped within the zonally localized jet region may exhibit the following behaviors: (1) an oscillatory decay to a steady state at the jet center, (2) the creation of additional SRWs within the jet region, or (3) a steady-state wherein the solution has a smoothed step-like structure located downstream along the jet axis.     

A channel head stability criterion for first‐order catchments

The present analysis derives a stability criterion for long‐term equilibrium channel heads. The concept of finite perturbation analysis is presented, during which the surface is subjected to perturbations of a finite amplitude and resulting changes in flow path structure and slope are computed. Based on these quantities the analysis predicts whether the perturbed location is going to erode, be filled in or remain steady. The channel head is defined geometrically as the focus point of converging flow lines at the bottom of hollows. It is demonstrated that stability at the channel head grows out of the competition between the rate of flow path convergence and the degree of profile concavity. Analytical functions are derived to compute channel head‐contributing area and ‐slope, flow path convergence and profile concavity as a function of perturbation depth, distance from the crest and the initial slope. In a numerical model these quantities point to the long‐term equilibrium channel head position, which is shown to depend also on the width to length ratio of hollows. It is also demonstrated that the equilibrium channel head position is sensitive to the base‐level lowering/non‐dimensional slope length ratio and to the slope of the initial topography. Morphometrical measurements both in the field and on simulated topographies were used to test the theoretical predictions.     

Hydromagnetic waves in rapidly rotating spherical shells generated by poloidal decay modes

Abstract

This paper presents the first attempt to examine the stability of a poloidal magnetic field in a rapidly rotating spherical shell of electrically conducting fluid. We find that a steady axisymmetric poloidal magnetic field loses its stability to a non-axisymmetric perturbation when the Elsasser number A based on the maximum strength of the field exceeds a value about 20. Comparing this with observed fields, we find that, for any reasonable estimates of the appropriate parameters in planetary interiors, our theory predicts that all planetary poloidal fields are stable, with the possible exception of Jupiter. The present study therefore provides strong support for the physical relevance of magnetic stability analysis to planetary dynamos. We find that the fluid motions driven by magnetic instabilities are characterized by a nearly two-dimensional columnar structure attempting to satisfy the Proudman-Taylor theorm. This suggests that the most rapidly growing perturbation arranges itself in such a way that the geostrophic condition is satisfied to leading order. A particularly interesting feature is that, for the most unstable mode, contours of the non-axisymmetric azimuthal flow are closely aligned with the basic axisymmetric poloidal magnetic field lines. As a result, the amplitude of the azimuthal component of the instability is smaller than or comparable with that of the poloidal component, in contrast with the instabilities generated by toroidal decay modes (Zhang and Fearn, 1994). It is shown, by examining the same system with and without fluid inertia, that fluid inertia plays a secondary role when the magnetic Taylor number T_m ? 10⁵. We find that the direction of propagation of hydromagnetic waves driven by the instability is influenced strongly by the size of the inner core.