Specific features of development of two-dimensional convection in a rotating compressible conducting fluid

An axisymmetric model of convection in a rotating cylinder in an external uniform magnetic field has been considered. In the considered model, the meridional circulation is created by a nonuniform rotation of the lower boundary relative to the other boundaries. In the considered model, the time of formation of the stationary regime in the magnetic field considerably increases if the vertical density (compressibility) inhomogeneity is taken into account for Ekman numbers of E = E _M = 3 × 10⁻³. This example shows that the compressibility of a medium should be taken into account in the convection and dynamics of the magnetic field when the magnetohydrodynamics of the Earth is analyzed.     

Flow past spheres and spheroids in the presence of a toroidal magnetic field

Summary Stewartson [1]²) has considered the inviscid flow past a sphere in the presence of a uniform magnetic field andMurray andLudford [2] have investigated a similar problem in which the magnetic field originates from an axially symmetric dipole field situated at the centre of the sphere. In connection with the study of earth's magnetic field, the toroidal part of this field plays a dominant part. This gives rise to the importance of studying the effect of a toroidal magnetic field on flows past different bodies of revolution; specially past spheres and spheroids. In the present note inviscid flows past a sphere, and a spheroid, are considered, for the case of a toroidal magnetic field originating in the fluid. In the case of the sphere the field inside the sphere consists of an electric dipole directed along the axis of symmetry together with a uniform electric field which produces a uniform current along the axis. In the case of the spheroid, the field inside it is due to an electric dipole and quadrupole directed along the axis of symmetry, together with a uniform electric field which produces a uniform current along this axis.     

An analytical model for ice-edge upwelling

Abstract

This paper presents an analytical, two-dimensional model of the wind-induced homogeneous circulation near the edge of an ice pack floating on the ocean surface. It is shown that a vertical shear layer arises under the ice edge, by which the wind-driven geostrophic motion in the open ocean is matched to the flow region underneath the ice. As in coastal upwelling models, this shear layer consists of a thin E ^1/2-layer inside a thicker E ^1/4-layer (E being the Ekman number). Under certain conditions the shear layer produces a vertical mass flux from the bottom to the surface Ekman layer. Near the surface this upwelling flux is concentrated in the narrow E ^1/2-layer. Comparison with observations of upwelling at the edge of a polar ice pack shows good agreement.     

α2-Dynamos and taylor's constraint

Abstract

An idealised α²ω-dynamo is considered in which the α-effect is prescribed. The additional ω-effect results from a geostrophic motion whose magnitude is determined indirectly by the Lorentz forces and Ekman suction at the boundary. As the strength of the α-effect is increased, a critical value α? _c is reached at which dynamo activity sets in; α? _c is determined by the solution of the kinematic α²-dynamo problem. In the neighbourhood of the critical value of α? the magnetic field is weak of order E ^1/4(μηρω)^½ due to the control of Ekman suction; E(?1) is the Ekman number. At certain values of α?, viscosity independent solutions are found satisfying Taylor's constraint. They are identified by the bifurcation of a nonlinear eigenvalue problem. Dimensional arguments indicate that following this second bifurcation the magnetic field is strong of order (μηρω)^½. The nature of the transition between the kinematic linear theory and the Taylor state is investigated for various distributions of the α-effect. The character of the transition is found to be strongly model dependent.     

Atmosphere and Earth's rotation

The theoretical aspects of the transfer of angular momentum between atmosphere and Earth are treated with particular emphasis on analytical solutions. This is made possible by the consequent usage of spherical harmonics of low degree and by the development of large-scale atmospheric dynamics in terms of orthogonal wave modes as solutions of Laplace's tidal equations.An outline of the theory of atmospheric ultralong planetary waves is given leading to analytical expressions for the meridional and height structure of such waves. The properties of the atmospheric boundary layer, where the exchange of atmospheric angular momentum with the solid Earth takes place, are briefly reviewed. The characteristic coupling time is the Ekman spin-down time of about one week.The axial component of the atmospheric angular momentum (AAM), consisting of a pressure loading component and a zonal wind component, can be described by only two spherical functions of latitude : the zonal harmonicP ₂ ⁰ (), responsible for pressure loading, and the spherical functionP ₁ ¹ () simulating supperrotation of the zonal wind. All other wind and pressure components merely redistributeAAM internally such that their contributions toAAM disappear if averaged over the globe. It is shown that both spherical harmonics belong to the meridional structure functions of the gravest symmetric Rossby-Haurwitz wave (0, –1)^*. This wave describes retrograde rotation of the atmosphere within the tropics (the tropical easterlies), while the gravest symmetric external wave mode (0, –2) is responsible for the westerlies at midlatitudes. Applying appropriate lower boundary conditions and assuming that secular angular momentum exchange between solid Earth and atmosphere disappears, the sum of both waves leads to an analytical solution of the zonal mean flow which roughly simulates the observed zonal wind structure as a function of latitude and height. This formalism is used as a basis for a quantitative discussion of the seasonal variations of theAAM within the troposphere and middle atmosphere.Atmospheric excitation of polar motion is due to pressure loading configurations, which contain the antisymmetric functionP ₂ ¹ () exp(i) of zonal wavenumberm=1, while the winds must have a superrotation component in a coordinate system with the polar axis within the equator. The Rossby-Haurwitz wave (1, –3)^* can simulate well the atmospheric excitation of the observed polar motion of all periods from the Chandler wobble down to normal modes with periods of about 10 days. Its superrotation component disappears so that only pressure loading contributes to polar motion.The solar gravitational semidiurnal tidal force acting on the thermally driven atmospheric solar semidiurnal tidal wave can accelerate the rotation rat of the Earth by about 0.2 ms per century. It is speculated that the viscous-like friction of the geomagnetic field at the boundary between magnetosphere and solar wind may be responsible for the westward drift of the dipole component of the internal geomagnetic field. Electromagnetic or mechanical coupling between outer core and mantle may then contribute to a decrease of the Earth's rotation rate.     

Rotating flow over long shallow ridges

Abstract

The flow of a rotating homogeneous, incompressible fluid past a long ridge is investigated. An analysis is presented for flows in which E ? 1, Ro ～ E^½, H/D ～ E⁰, h/D ～ E^½ and cosα ～ E⁰ where E is the Ekman number, Ro the Rossby number, H/D the fluid depth to ridge width ratio, h/D the ridge height to ridge width ratio and α the angle between the free stream flow and a line perpendicular to the ridge axis. The analysis includes effects of the nonlinear inertial terms. Particular examples of a ridge of triangular cross section and a sinusoidal topography are investigated in some detail. Experiments are presented for a triangular ridge which are in good agreement with the theory.     

On the spin-up of a stably stratified,electrically conducting fluid. Part 1: Spin-up controlled by the hartmann layer

Abstract

The linear spin-up of a stably stratified, electrically conducting fluid within an electrically insulating cylindrical container in the presence of an applied axial magnetic field is analyzed for those cases in which electric currents generated within the steady Hartmann boundary layer control the fluid interior. It is shown how to obtain the known spin-up times for a homogeneous, nonconducting fluid (τ = E ^-½), a stably stratified, nonconducting fluid (τ = (σS/E, E ^?1) and a homogeneous conducting fluid (τ = α^?1 E ^-½) from the present formulation where τ = v/ωt, E = v/ωL ², σS = vN²/κω² and 2α² = σB²/pω. The problem is solved in the parameter range E?α²?1, α²/E?σS using the Laplace transform and two new spin-up times are obtained. Combined into one expression, they are τ = (1 + δ)α^?1E^-½ where δ = σμv. The spin-up mechanism is investigated and it is found that, in contrast to the homogeneous, conducting case, torsional Alfvén waves may be instrumental in the spin-up of a stratified conducting fluid. The effects of viscous and ohmic diffusion on the torsional Alfvén wave fronts are studied and the following regimes are identified: 0 < δ ?E/α², spin-up by meridional circulation of electric current with no Alfvén waves; E^-½/α ? δ ? 1, spin-up by meridional circulation of electric current with transient Alfvén waves; α/E^½ ? δ ? α²/E, spin-up by meridional circulation of current with weak Alfvén waves; 1 ? δ ? α/E^½, spin-up by strong Alfvén waves; α^½/E ? δ ? α²/E, spin-up by viscous diffusion with transient Alfvén waves; α/E ? δ < ∞, spin-up by viscous diffusion with no Alfvén waves.     

Experiments in ocean circulation modelling

Abstract

In a laboratory model ocean, fluid in a rotating tank of varying depth is subjected to “wind-stress”, For a certain range of the parameters, Ekman number E and Rossby number R, a homogeneous fluid displays steady, westward intensified flow. For the same range of E and R, a two-layer fluid can have baroclinic instabilities. The parameter range for the various kinds of instabilities is mapped in a regime diagram. The northward transport in the western boundary current is measured as it varies with Rossby number for both homogeneous and two-layer fluid.     

Coastal upwelling,circulation and heat balance around New Caledonia’s barrier reef

An outstanding characteristic of New Caledonia upwelling is that most events appear limited to the southern half of the western barrier reef. This north–south difference cannot be explained by alongshore variability of the projected wind stress and no strong evidence for alternative explanations has been proposed. A major objective of this paper is to provide the first dynamical analysis of New Caledonia upwelling and its regional environment, based on numerical simulations. Coastal upwelling around New Caledonia is shown to be modulated by a system of geostrophic currents interacting with the island mass. Upwelling velocities are weaker than expected from the two-dimensional Ekman theory, as Ekman divergence is balanced by “coastal geostrophic convergence”. The cooling effect of upwelling is also attenuated by alongshore transport of warm water by the Alis current, reminiscent of the Leeuwin current off Western Australia. Nevertheless, coastal upwelling can locally modify the large-scale surface water heat budget, dominated by meridional advection warming and surface cooling. The upwelled waters appear to be mostly of western origin and are transported below the surface by the Subtropical Counter Current before upwelling off New Caledonia. This appears in sharp contrast with the eastern barrier reef where the general warming by meridional advection of tropical surface waters is accentuated by the vigorous western boundary type Vauban current.     

On the Mechanism of the Secular Tidal Acceleration of the Solid Inner Core and the Viscosity of the Liquid Core

As is known, the secular deceleration of the Earth's diurnal rotation is explained mainly by the tidal friction in the ocean. Below we consider this mechanism in some detail, taking into account also elastic deformations of the mantle under the action of ocean loading and the interaction between the tide-generating body, ocean tidal wave, liquid outer core, and solid inner core. It is shown that elastic displacements of the core-mantle boundary under the action of ocean loading are of about the same amplitude and phase as the elastic loading displacements of the Earth's outer surface. As a result, side by side with the mechanism of secular deceleration of diurnal rotation of the mantle, there are also (1) the opposite mechanism of secular acceleration of diurnal rotation of the outer liquid core and of the solid inner core and (2) the mechanism of excitation of differential rotation in the liquid core. Taking these effects into account, we compare theoretical and modern observed data on the eastward drift of the solid inner core. It is shown that the best agreement may be obtained if the turbulent viscosity of the liquid core is about 2 × 10 ³ Poise     

Thermal and magnetic instabilities in a rapidly rotating fluid sphere

Abstract

A linear analysis is used to study the stability of a rapidly rotating, electrically-conducting, self-gravitating fluid sphere of radius r ₀, containing a uniform distribution of heat sources and under the influence of an azimuthal magnetic field whose strength is proportional to the distance from the rotation axis. The Lorentz force is of a magnitude comparable with that of the Coriolis force and so convective motions are fully three-dimensional, filling the entire sphere. We are primarily interested in the limit where the ratio q of the thermal diffusivity κ to the magnetic diffusivity η is much smaller than unity since this is possibly of the greatest geophysical relevance.

Thermal convection sets in when the temperature gradient exceeds some critical value as measured by the modified Rayleigh number R_c. The critical temperature gradient is smallest (R_c reaches a minimum) when the magnetic field strength parameter Λ ? 1. [R_c and Λ are defined in (2.3).] The instability takes the form of a very slow wave with frequency of order κ/r ² ₀ and its direction of propagation changes from eastward to westward as Λ increases through Λ _c ? 4.

When the fluid is sufficiently stably stratified and when Λ > Λ_m ? 22 a new mode of instability sets in. It is magnetically driven but requires some stratification before the energy stored in the magnetic field can be released. The instability takes the form of an eastward propagating wave with azimuthal wavenumber m = 1.     

On the steady mass transport induced by wave motions in a viscous rotating fluid

Abstract

The steady second-order motion induced by a first-order wave motion in a homogeneous, viscous and rotating fluid is examined. If the wave motion produces a steady Ekman layer suction by non-linear interactions, this suction must induce a steady component of interior, relative vorticity parallel to the axis of rotation in order to conserve mass. A boundary value problem for the determination of the induced, steady interior mass transport velocity is presented. The mass transport induced by a Kelvin wave is examined as an illustration and possible application of the theory.     

A study on three-dimensional structures of the ionospheric dynamo currents induced by the neutral winds simulated by the Kyushu-GCM

Three-dimensional structures of the ionospheric dynamo currents are examined using the neutral winds in a general circulation model of the middle atmosphere at Kyushu University. A quasi-three-dimensional ionospheric dynamo model is constructed assuming an infinite parallel conductivity in the ionosphere. This model is able to simulate both the equatorial electrojet and the global Sq current system successfully. The simulated results reveal that the equatorial electrojet is confined in quite narrow latitudes around the equator accompanied with meridional current circulations and satisfies a non-divergent structure mainly within the E region. A vertically stratified double layered structure is seen in the east–west current density near the focus latitude of the global Sq current system. It is shown that the stratified structure mainly consists of the east–west Hall current associated with the eastward wind of zonal wavenumbers 1 and 2 in the lower altitudes and the westward wind of zonal wavenumber 2 in the upper altitudes. The day-to-day variation of the neutral winds can significantly vary the induced ionospheric dynamo current system, which is recognized as changes of the focus latitude and/or the maximum value of the equatorial electrojet.     

The structure of separated flow past a circular cylinder in a rotating frame

Abstract

This paper examines the detailed E ^1/4-layer structure of separated flow past a circular cylinder in a low-Rossby-number rotating fluid as the Ekman number E tends to zero. This structure is based on an initial proposal by Page (1987) but with some modifications in response to further evidence, outlined both in this paper and elsewhere, on the behaviour of E ^1/4-layer flows in this context. Numerical calculations for flow in an E ^1/4 shear layer along the separated free streamline are described and the mass flux from this layer is then used to calculate the higher-order flow within the separation bubble. The flow structure is found to have two forms, depending on the value of the O(1) parameter λ, and these are compared with results from published “Navier-Stokes” type calculations for the flow at small but finite values of E.     

Asymptotic treatment of the stability of a rotating layer of fluid with rigid boundaries

Abstract

Chandrasekhar (1961) has summarized the stability results of Bénard convection in a rotating fluid for the cases where the boundary surfaces are both rigid and free, and for both exchange of stabilities and overstability. His analysis provides very accurate results for a limited range of Taylor number J. Bisshopp and Niiler (1965) presented an asymptotic analysis of the rigid boundary problem for exchange of stabilities which is valid for very large Taylor number. The present paper makes use of modern rotating fluid theory to develop an approximate scheme for evaluating the Rayleigh number and other parameters and variables. Known asymptotic results for the free boundary problem at large J are used and an expansion in powers of E^1/6 (the Ekman number, E = 2J ^?½) yields a sequence of equations and appropriate boundary conditions for the rigid boundary problem. After the algorithm for the calculation is developed, results are given for the problem to second order in the expansion parameter for the case of exchange of stabilities and to first order in the expansion parameters for the overstable case. Ekman boundary layers are important in the development as one might expect. However, an additional, diffusive boundary layer of thickness E^? is necessary to provide the details of the temperature field. This boundary layer is the thermal response in the vertical direction to the horizontal spacing of the cells which is also order E^?. The horizontal spacing of the cells is essentially a series of contiguous, Stewartson (1957) layers of thickness E^?.     

Comment on rotating hydraulics

Abstract

Numerical solutions of the axisymmetric flows during the relatively early phase of spin-up from rest of a stratified fluid in a cylinder are presented. Detailed results are given for a cylinder of aspect ratio of O(l) and for a minute Ekman number, showing axisymmetric spin-up for three values of the stratification parameter. As the stratification increases, the meridional circulation is confined to a region closer to the Ekman layers. An axisymmetric shear wave propagates radially inward from the sidewall, but, unlike the strictly vertical front for a homogeneous fluid, the interface which separates rotating from nonrotating fluid is bow-shaped. For a stratified fluid, the axial vorticity distribution is nonuniform both in the vertical and in the radial directions. With increasing stratification, diffusive vorticity production near the sidewall is more pronounced. Axisymmetric flows in the early phase of spin-up of a stratified fluid are controlled by both the inviscid dynamic effect and the viscous diffusion effect. At a location close to the Ekman layers, the inviscid effect outweighs the viscous effect, in much the same way as in a homogeneous fluid. However, at a location close to mid-depth, the viscous diffusion effect, enhanced by substantial flow gradients in that region, is dominant. This points to the necessity of including the direct effect of viscous diffusion in the interior in formulating an analytical model of stratified spin-up problems.     

Evolution of non-linear α 2-dynamos and taylor's constraint

A key non-linear mechanism in a strong-field geodynamo is that a finite amplitude magnetic field drives a flow through the Lorentz force in the momentum equation and this flow feeds back on the field-generation process in the magnetic induction equation, equilibrating the field. We make use of a simpler non-linear?α?²-dynamo to investigate this mechanism in a rapidly rotating fluid spherical shell. Neglecting inertia, we use a pseudo-spectral time-stepping procedure to solve the induction equation and the momentum equation with no-slip velocity boundary conditions for a finitely conducting inner core and an insulating mantle. We present calculations for Ekman numbers (E) in the range 2.5× 10^?3 to 5.0× 10^?5, for?α?=α ₀cos?θ?sin?π?(r?r_i ) (which vanishes on both inner and outer boundaries). Solutions are steady except at lower E and higher values of?α?₀. Then they are periodic with a reversing field and a characteristic rapid increase then equally rapid decrease in magnetic energy. We have investigated the mechanism for this and shown the influence of Taylor's constraint. We comment on the application of our findings to numerical hydrodynamic dynamos.     

Strong Potential Magnetic Field Influence on Slightly Differentially Rotating Spherical Shell

An electrically conducting viscous fluid-filled spherical shell is permeated by an axisymmetric strong potential magnetic field with large Elssaser number ² 1. We describe analytically the steady flow driven by a slightly faster rotation of the conducting inner boundary of the shell. The main flow is controlled by Ekman-Hartmann boundary layers with a small thickness /, where ² is the Ekman number. Asymptotics based on small ^–1 1 reveal the nature of a free shear layer O((/)^1/2) and a super-rotation that allows a part of the fluid to rotate faster than the inner sphere. The free shear is following an imposed field line that is tangent to the inner or outer sphere. Meridional flux is concentrated in the shear and boundary layers. Fluid tends to rotate with the inner sphere and to expel azimuthal magnetic field from an -region restricted by the free shear in the spherical shell. For an imposed axial uniform magnetic field, this -region is outside the cylinder tangent to the inner sphere and rotates with the outer sphere. Weak differential rotation O(/) is inside the cylinder, while almost all difference in rotation rates between spheres is accommodated in the thin O((/)^1/2) free shear. For an imposed dipole magnet, the region has a shape of a lobe touching the outer equator. Inside a super-rotation exists; this is the common case for such when the source of the imposed field is inside.     

Simultaneous observations at different altitudes of ionospheric backscatter in the eastward electrojet

A common feature of evening near-range ionospheric backscatter in the CUTLASS Iceland radar field of view is two parallel, approximately L-shell-aligned regions of westward flow which are attributed to irregularities in the auroral eastward electrojet region of the ionosphere. These backscatter channels are separated by approximately 100–200 km in range. The orientation of the CUTLASS Iceland radar beams and the zonally aligned nature of the flow allows an approximate determination of flow angle to be made without the necessity of bistatic measurements. The two flow channels have different azimuthal variations in flow velocity and spectral width. The nearer of the two regions has two distinct spectral signatures. The eastern beams detect spectra with velocities which saturate at or near the ion-acoustic speed, and have low spectral widths (less than 100ms^–1), while the western beams detect lower velocities and higher spectral widths (above 200ms^–1). The more distant of the two channels has only one spectral signature with velocities above the ionacoustic speed and high spectral widths. The spectral characteristics of the backscatter are consistent with E-region scatter in the nearer channel and upper-E-region or F-region scatter in the further channel. Temporal variations in the characteristics of both channels support current theories of E-region turbulent heating and previous observations of velocity-dependent backscatter cross-section. In future, observations of this nature will provide a powerful tool for the investigation of simultaneous E- and F-region irregularity generation under similar (nearly co-located or magnetically conjugate) electric field conditions.     

Models of solar differential rotation

Abstract

Models of a differentially rotating compressible convection zone are calculated, considering the inertial forces in the poloidal components of the equations of motion. Two driving mechanisms have been considered: latitude dependent heat transport and anisotropic viscosity. In the former case a meridional circulation is induced initially which in turn generates differential rotation, whereas in the latter case differential rotation is directly driven by the anisotropic viscosity, and the meridional circulation is a secondary effect.

In the case of anisotropic viscosity the choice of boundary conditions has a big influence on the results: depending on whether or not the conditions of vanishing pressure perturbation are imposed at the bottom of the convection zone, one obtains differential rotation with a fast (≥ 10 ms^?1) or a slow (～ 1 ms^?1) circulation. In the latter case the rotation law is mainly a function of radius and the rotation rate increases inwards if the viscosity is larger in radial direction than in the horizontal directions.

The models with latitude dependent heat transport exhibit a strong dependence on the Prandtl number. For values of the Prandtl number less than 0.2 the pole-equator temperature difference and the surface velocity of the meridional circulation are compatible with observations. For sufficiently small values of the Prandtl number the convection zone becomes globally unstable like a layer of fluid for which the critical Rayleigh number is exceeded.