Cross helicity and related dynamo

The turbulent cross helicity is directly related to the coupling coefficients for the mean vorticity in the electromotive force and for the mean magnetic-field strain in the Reynolds stress tensor. This suggests that the cross-helicity effects are important in the cases where global inhomogeneous flow and magnetic-field structures are present. Since such large-scale structures are ubiquitous in geo/astrophysical phenomena, the cross-helicity effect is expected to play an important role in geo/astrophysical flows. In the presence of turbulent cross helicity, the mean vortical motion contributes to the turbulent electromotive force. Magnetic-field generation due to this effect is called the cross-helicity dynamo. Several features of the cross-helicity dynamo are introduced. Alignment of the mean electric-current density J with the mean vorticity Ω , as well as the alignment between the mean magnetic field B and velocity U , is supposed to be one of the characteristic features of the dynamo. Unlike the case in the helicity or α effect, where J is aligned with B in the turbulent electromotive force, we in general have a finite mean-field Lorentz force J ?×? B in the cross-helicity dynamo. This gives a distinguished feature of the cross-helicity effect. By considering the effects of cross helicity in the momentum equation, we see several interesting consequences of the effect. Turbulent cross helicity coupled with the mean magnetic shear reduces the effect of turbulent or eddy viscosity. Flow induction is an important consequence of this effect. One key issue in the cross-helicity dynamo is to examine how and how much cross helicity can be present in turbulence. On the basis of the cross-helicity transport equation, its production mechanisms are discussed. Some recent developments in numerical validation of the basic notion of the cross-helicity dynamo are also presented.     

Lateral friction in the oceans as an effect of potential vorticity mixing

Abstract

Applying a mixing-length calculation to potential vorticity rather than to momentum a new type of lateral friction appears in the oceanic mass transport equations. This friction is evaluated for the special case of horizontally homogeneous, quasi-geostrophic turbulence. The main effect is a westward force arising from the so-called β-term. This produces an additional southward interior transport and a strengthening of the western boundary current. A turbulent exchange coefficient K^H = 10⁸ em²s^?1 is sufficient to give a Gulf Stream transport twice that obtained by the classical Sverdrup model.     

Turbulent transport of magnetic fields. V. Distribution of magnetic energy in a simple α2-dynamo

Abstract

In this paper we analyse the stationary mean energy density tensor T_ij = B_iB_j for the x ²-sphere. This model is one of the simplest possible turbulent dynamos, originally due to Krause and Steenbeck (1967): a conducting sphere of radius R with homogeneous, isotropic and stationary turbulent convection, no differential rotation and negligible resistivity. The stationary solution of the (linear) equation for T_ij is found analytically. Only T_rr , T _θθ and T _φφ are unequal to zero, and we present their dependence on the radial distance r.

The stationary solution depends on two coefficients describing the turbulent state: the diffusion coefficient β≈?u²?τ_c/3 and the vorticity coefficient γ ≈ ?|?×u|²?τ_c/3 where u(r, t) is the turbulent velocity and _c its correlation time. But the solution is independent of the dynamo coefficient α≈??u·?×u?τ_c/3 although α does occur in the equation for T_ij . This result confirms earlier conclusions that helicity is not required for magnetic field generation. In the stationary state, magnetic energy is generated by the vorticity and transported to the boundary, where it escapes at the same rate. The solution presented contains one free parameter that is connected with the distribution of B over spatial scales at the boundary, about which T_ij gives no information. We regard this investigation as a first step towards the analysis of more complicated, solar-type dynamos.     

Superhelicity,helicity and potential vorticity

Abstract

‘‘Helicity'’ density H ≡ u · ω and other pseudo-scalar fields such as P ≡ ω · Vlnρ (which is related to Ertel potential vorticity) are useful quantities in theoretical fluid dynamics and magneto-fluid dynamics. Here u denotes the Eulerian flow velocity relative to the chosen frame of reference, ω ≡ V × u is the corresponding relative vorticity and ρ the mass density of the fluid. A general expression is readily obtained for ?H/?t (where t denotes time) in terms of P and the ‘‘superhelicity'’ density S ≡ ω · V × ω which, in fluids of low viscosity, has its highest values in boundary layers. One need for such a relationship became evident during an attempt to interpret the findings of laboratory experiments on thermal convection in rotating fluids in containers of various geometrical shapes and topological characteristics.

In electrodynamics an analogous expression can be found relating the time rate of change of ‘‘magnetic helicity'’ A · B to ‘‘magnetic superhelicity'’ B · ? × B (where B · ? × A is the magnetic field) and a scalar quantity analogous to P which involves non-Ohmic contributions to the relationship between the electric current density and the electric field.     

The nonlinear spin-up of a stratified ocean

Abstract

The adjustment of a nonlinear, quasigeostrophic, stratified ocean to an impulsively applied wind stress is investigated under the assumption that barotropic advection of vortex tube length is the most important nonlinearity. The present study complements the steady state theories which have recently appeared, and extends earlier, dissipationless, linear models.

In terms of Sverdrup transport, the equation for baroclinic evolution is a forced advection-diffusion equation. Solutions of this equation subject to a “tilted disk” Ekman divergence are obtained analytically for the case of no diffusion and numerically otherwise. The similarity between the present equation and that of a forced barotropic fluid with bottom topography is shown.

Barotropic flow, which is assumed to mature instantly, can reverse the tendency for westward propagation, and thus produce regions of closed geostrophic contours. Inside these regions, dissipation, or equivalently the eddy field, plays a central role. We assume that eddy mixing effects a lateral, down-gradient diffusion of potential vorticity; hence, within the closed geostrophic contours, our model approaches a state of uniform potential vorticity. The solutions also extend the steady-state theories, which require weak diffusion, by demonstrating that homogenization occurs for moderately strong diffusion.

The evoiution of potential vorticity and the thermocline are examined, and it is shown that the adjustment time of the model is governed by dissipation, rather than baroclinic wave propagation as in linear theories. If dissipation is weak, spin-up of a nonlinear ocean may take several times that predicted by linear models, which agrees with analyses of eddy-resolving general circulation models. The inclusion of a western boundary current may accelerate this process, although dissipation will still play a central role.     

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 Reynolds-averaged turbulence modelling approach to the maintenance of the Venus superrotation

Abstract

A maintenance mechanism of an approximately linear velocity profile of the Venus zonal flow or superrotation is explored, with the aid of a Reynolds-averaged turbulence modelling approach. The basic framework is similar to that of Gierasch (Meridional circulation and maintenance of the Venus atmospheric rotation. J. Atmos. Sci. 1975, 32, 1038–1044) in the sense that the mechanism is examined under a given meridional circulation. The profile mimicking the observations of the flow is initially assumed, and its maintenance mechanism in the presence of turbulence effects is investigated from a viewpoint of the suppression of energy cascade. In the present work, the turbulent viscosity is regarded as an indicator of the intensity of the cascade. A novelty of this formalism is the use of the isotropic turbulent viscosity based on a non-local time scale linked to a large-scale flow structure. The mechanism is first discussed qualitatively. On the basis of these discussions, the two-dimensional numerical simulation of the proposed model is performed, with an initially assumed superrotation, and the fast zonal flow is shown to be maintained, compared with the turbulent viscosity lacking the non-local time scale. The relationship of the present model with the current general circulation model simulation is discussed in light of a crucial role of the vertical viscosity.     

Numerical computation of the tidally induced Lagrangian residual current in a model bay

With a depth-averaged numerical model, the tidally induced Lagrangian residual current in a model bay was studied. To correctly reflect the long-term mass transport, it is appropriate to use the Lagrangian residual velocity (LRV) rather than the Eulerian residual velocity (ERV) or the Eulerian residual transport velocity (ETV) to describe the residual current. The parameter κ, which is defined as the ratio of the typical tidal amplitude at the open boundary to the mean water depth, is considered to be the indicator of the nonlinear effect in the system. It is found that the feasibility of making the mass transport velocity (MTV) approximate the LRV is strongly dependent on κ. The error between the MTV and the LRV tends to increase with a growing κ. An additional error will come from the various initial tidal phases due to the Lagrangian drift velocity (LDV) when κ is no longer small. According to the residual vorticity equation based on the MTV, the Coriolis effect is found to influence the residual vorticity mainly through the curl of the tidal stress. A significant difference in the flow pattern indicates that the LRV is sensitive to the bottom friction in different forms.     

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

Abstract

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

Deep circulation driven by strong vertical mixing in the Timor Basin

The importance of deep mixing in driving the deep part of the overturning circulation has been a long debated question at the global scale. Our observations provide an illustration of this process at the Timor Basin scale of ～1000 km. Long-term averaged moored velocity data at the Timor western sill suggest that a deep circulation is present in the Timor Basin. An inflow transport of ～0.15 Sv is observed between 1600 m and the bottom at 1890 m. Since the basin is closed on its eastern side below 1250 m depth, a return flow must be generated above 1600 m with a ～0.15 Sv outflow. The vertical turbulent diffusivity is inferred from a heat and transport balance at the basin scale and from Thorpe scale analysis. Basin averaged vertical diffusivity is as large as 1 × 10^?3 m² s^?1. Observations are compared with regional low-resolution numerical simulations, and the deep observed circulation is only recovered when a strong vertical diffusivity resulting from the parameterization of internal tidal mixing is considered. Furthermore, the deep vertical mixing appears to be strongly dependent on the choice of the internal tide mixing parameterization and also on the prescribed value of the mixing efficiency.     

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.     

Sill effects on physical dynamics in eastern Long Island Sound

This study investigates how Mattituck Sill influences circulation patterns and physical dynamics in eastern Long Island Sound, a major estuary on the U.S. east coast. Observations show there is pronounced across-estuary transport in the area and suggest there may be subtidal anticyclonic flow around the sill. Model runs, with and without sill bathymetry, exhibit this across-estuary transport and anticyclonic circulation. Comparison between these runs indicates that the sill intensifies the anticyclonic circulation. This study finds the sill does not exert internal hydraulic control during neap, mean, or spring tidal conditions. Nevertheless, along-estuary exchange is reduced over the sill and across-estuary fluxes are increased. The Connecticut River plume enters close to the estuary mouth. The sill deflects more of the plume waters towards the mouth, causing less freshwater to take the long looping route through the estuary. The subtidal circulation balance around the sill indicates a barotropic balance between the tidal advection of tidal vorticity and friction. The subtidal vorticity balance indicates the net effect of tidal advection of relative vorticity is balanced with frictional curl associated with lateral speed gradients and vorticity dissipation. Previously developed scalings based on the circulation balance (Nature 290:549–555, 1981), frictional vorticity generation mechanisms (Deep-Sea Res 28:195–212, 1981), and tidal diffusion of potential vorticity (J Phys Oceanogr 29:821–827, 1999) are applicable to Mattituck Sill and predict circulation with a similar magnitudes to model results.     

Diagnosing vertical motion in the Equatorial Atlantic

Estimating the vertical velocity (w) in the oceanic upper-layers is a key issue for understanding the cold tongue development in the Eastern Equatorial Atlantic. In this methodological paper, we develop an expanded and general formulation of the vertical velocity equation based on the primitive equation (PE) system, in order to gain new insight into the physical processes responsible for the Equatorial and Angola upwellings. This approach is more accurate for describing the real ocean than simpler considerations based on just the wind-driven patterns of surface layer divergence. The w-sources/forcings are derived from the PE w-equation and diagnosed from a realistic ocean simulation of the Equatorial Atlantic. Sources of w are numerous and express the high complexity of terms related to the turbulent momentum flux, to the circulation and to the mass fields, some of them depending explicitly on w and others not. The equatorial upwelling is found to be mainly induced by the (i) the zonal turbulent momentum flux, (ii) the curl of turbulent momentum flux and (iii) the imbalance between the circulation and the pressure fields. The Angola upwelling in the eastern part of the basin is controlled by strong curl of turbulent momentum flux. A strong cross-regulation is evidenced between the w-forcings independent of w and dependent on w, which suggests an equatorial balanced-dynamics. The w-forcing depending on w represents the negative feedback of the ocean to the w-forcing independent of w: in the equatorial band, this adjustment is led by non-linear processes and by vortex stretching outside.     

Turbulent “Dynamo” in the Passive Scalar Problem

A turbulent magnetic dynamo can be considered as the evolution of a vector field in a turbulent fluid flow. The problem of evolution of scalar fields (e.g., number density of small particles) in a turbulent fluid flow is similar to the turbulent magnetic dynamo. The dynamo instability results in generation of magnetic field. The most important effect which can cause a generation of mean magnetic field in a turbulent fluid flow is the -effect: = – (1/3) u · ( × u), where u is the turbulent velocity field with the correlation time . A similar instability in the passive scalar problem results in formation of large-scale inhomogeneous structures in a spatial distribution of particles due to the -effect: = u_p ( · u_p), where u _p is the random velocity field of the particles which they acquire in a turbulent fluid velocity field. The effect is caused by inertia of particles which results in divergent velocity field of the particles. This results in additional turbulent nondiffusive flux of particles. The mean-field dynamics of inertial particles are studied by considering the stability of the equilibrium solution of the derived evolution equation for the mean number density of the particles in the limit of large Péclet numbers. The resulting equation is reduced to an eigenvalue problem for a Schrödinger equation with a variable mass, and a modified Rayleigh-Ritz variational method is used to estimate the lowest eigenvalue (corresponding to the growth rate of the instability). This estimate is in good agreement with obtained numerical solution of the Schrödinger equation. Similar effects arise during turbulent transport of gaseous admixtures (or light noninertial particles) in a low-Mach-number compressible fluid flow. The discussed effects are important in planetary and atmospheric physics (cloud formation, pollutant dynamics, preferential concentration of particles in protoplanetary disks and also planetesimals in them).     

An explanation for salinity- and SPM-induced vertical countergradient buoyancy fluxes

Measurements of turbulent fluctuations of velocity, salinity, and suspended particulate matter (SPM) are presented. The data show persistent countergradient buoyancy fluxes. These countergradient fluxes are controlled by the ratio of vertical turbulent kinetic energy (VKE) and available potential energy (APE) terms in the buoyancy flux equation. The onset of countergradient fluxes is found to approximately coincide with larger APE than VKE. It is shown here that the ratio of VKE to APE can be written as the square of a vertical Froude number. This number signifies the onset of the dynamical significance of buoyancy in the transport of mass. That is when motions driven by buoyancy begin to actively determine the vertical turbulent transport of mass. Spectral and quadrant analyses show that the occurrence of countergradient fluxes coincides with a change in the relative importance of turbulent energetic structures and buoyancy-driven motions in the transport of mass. Furthermore, these analyses show that with increasing salinity-induced Richardson number (Ri), countergradient contributions expand to the larger scales of motions and the relative importance of outward and inward interactions increases. At the smaller scales, at moderate Ri, the countergradient buoyancy fluxes are physically associated with an asymmetry in transport of fluid parcels by energetic turbulent motions. At the large scales, at large Ri, the countergradient buoyancy fluxes are physically associated with convective motions induced by buoyancy of incompletely dispersed fluid parcels which have been transported by energetic motions in the past. Moreover, these convective motions induce restratification and enhanced settling of SPM. The latter is generally the result of salinity-induced convective motions, but SPM-induced buoyancy is also found to play a role.     

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.     

Inertial effects in an oceanic circulation model

Abstract

The flow in a mechanically driven thin barotropic rotating fluid system is analysed. The linear theory of Baker and Robinson (1969) is modified and extended into the non-linear regime.

An internal parameter, the “local Rossby number”, is indicative of the onset of nonlinear effects. If this parameter is 0(1) then inertial effects are as important as Coriolis accelerations in the interior of the transport-turning western boundary layer and both of its Ekman layers. The inertial effects in the Ekman layers, ignored in previous explorations of non-linear wind driven oceanic circulation, are retained here and calculated using an approximation of the Oseen type. The circulation problem is reduced to a system of scalar equations in only two independent variables; the system is valid for non-small local Rossby number provided only that the approximate total vorticity is positive.

To complete the solution for small Rossby number a boundary condition for the inertially induced transport is needed. It is found by examining the dynamics controlling this additional transport from the western boundary layer as the transport recirculates through the rest of the ocean basin. The strong constraint of total recirculation within the western boundary layer (zero net inertial transport) is derived.

The calculated primary inertial effects are in agreement with the observations of the laboratory model of Baker and Robinson (1969).

The analysis indicates the extent to which three-dimensional non-linear circulation can be reduced to a two dimensional problem.     

Interactions of isolated vortices I: Modons colliding with modons

Abstract

Collisions between isolated(i.e. localized) dipolar vortex states, called modons, are examined in various numerical solutions of the quasigeostrophic, equivalent barotropic equation. For a range of parameters, the collision interactions are soliton-like in that the vorticity maxima are displaced (phase-shifted) with only small speed changes and little excitation of internal degrees of freedom. For other parameters, new “inelastic” effects are observed, including speed changes due to vorticity rearrangement, vorticity filamentation, modon “capture” or “fusion” in an overtaking collision, and the “fission” of a modon into its component vorticity monopoles in a head-on collision.     

State analysis using the Local Ensemble Transform Kalman Filter (LETKF) and the three-layer circulation structure of the Luzon Strait and the South China Sea

A new circulation model of the western North Pacific Ocean based on the parallelized version of the Princeton Ocean Model and incorporating the Local Ensemble Transform Kalman Filter (LETKF) data assimilation scheme has been developed. The new model assimilates satellite data and is tested for the period January 1 to April 3, 2012 initialized from a 24-year simulation to estimate the ocean state focusing in the South China Sea (SCS). Model results are compared against estimates based on the optimum interpolation (OI) assimilation scheme and are validated against independent Argo float and transport data to assess model skills. LETKF provides improved estimates of the western North Pacific Ocean state including transports through various straits in the SCS. In the Luzon Strait, the model confirms, for the first time, the three-layer transport structure previously deduced in the literature from sparse observations: westward in the upper and lower layers and eastward in the middle layer. This structure is shown to be robust, and the related dynamics are analyzed using the results of a long-term (18 years) unassimilated North Pacific Ocean model. Potential vorticity and mass conservations suggest a basin-wide cyclonic circulation in the upper layer of the SCS (z?>??570 m), an anticyclonic circulation in the middle layer (?570 m?≥?z?>??2,000 m), and, in the abyssal basin (<?2,000 m), the circulation is cyclonic in the north and anticyclonic in the south. The cyclone–anticyclone abyssal circulation is confirmed and explained using a deep-layer reduced-gravity model as being caused by overflow over the deep sill of the Luzon Strait, coupled with intense, localized upwelling west of the strait.