A steady state of the disc dynamo

Abstract

We discuss the steady states of the αω-dynamo in a thin disc which arise due to α-quenching. Two asymptotic regimes are considered, one for the dynamo numberD near the generation thresholdD ₀, and the other for |D| ? 1. Asymptotic solutions for |D—D ₀| ? |D ₀| have a rather universal character provided only that the bifurcation is supercritical. For |D| ? 1 the asymptotic solution crucially depends on whether or not the mean helicity α, as a function ofB, has a positive root (hereB is the mean magnetic field). When such a root exists, the field value in the major portion of the disc is O(l), while near the disc surface thin boundary layers appear where the field rapidly decreases to zero (if the disc is surrounded by vacuum). Otherwise, when α = O(|B|^?s) for |B| → ∞, we demonstrate that |B| = O(|D|^1/s ) and the solution is free of boundary layers. The results obtained here admit direct comparison with observations of magnetic fields in spiral galaxies, so that an appropriate model of nonlinear galactic dynamos hopefully could be specified.     

Time domain scattering of acoustic plane waves by vertical faults

Acoustic plane wave scattering at a vertical fault structure represents the simplest two-dimensional model of geophysical exploration that can be investigated by analytical techniques. The exact and complete solution, in the time domain, for the scattering of the pressure field of an acoustic plane wave normally incident on a vertical fault structure is determined adapting previous results given for the frequency domain. The wave form of the pressure field of the incident plane wave is expressed by a causal time function that decays exponentially with time at every point above the fault (z<0). The zero-order term of the scattered pressure field has been computed above the fault. This zero-order term consists of an inverse Fourier transform which reduces to a closed expression forx=0, and contains an integral of a Hankel function forx#0. The high frequency part of the inverse Fourier transform forx#0 is computed employing asymptotic expressions for the Hankel function. The integral of the asymptotic expression of the Hankel function reduces to: (i) a Fresnel integral which contains a plane wave term for |x||z|; and (ii) a stationary point plane wave term plus an upper limit term for |x|=O(|z|). For the latter case the plane wave term cancels, leaving a cylindrical wave emanated from the edge of the fault. The wave front is well defined in shape, in phase and in amplitude. The amplitude of the scattered field is discontinuous atx=0, presents a jump and is well defined for |x| small and is rather smooth for |x| large.     

Echo intensity data as a directional antenna for observing processes above sloping ocean bottoms

Relative ‘echo intensity’ data (dI) from a bottom-mounted four-beam 300 kHz acoustic Doppler current profiler (ADCP) are used to infer propagation of vigorous processes above a continental slope. The 3- to 60-m horizontal beam spread and the 2-Hz sampling allow the distinction of different arrival times t _i , i = 1,..., 4, at different distances in the acoustic beams from sharp changes in dI-content associated with frontal non-linear and turbulent bores or ‘waves’. The changes in dI are partially due to variations in amounts of resuspended material carried by the near-bottom turbulence and partially due to the fast variations in density stratification (‘stratified turbulence’), as inferred from 1-Hz sampled thermistor string data above the ADCP. Such bores are observed to pass the mooring up to 80 m above the bottom, having typical propagation speeds c = 0.15–0.5 m s⁻¹, as determined from dI(t _i ). Particle speeds in the immediate environment of a bore amount to |u|_env=c ± 0.05 m s⁻¹, the equality being a necessary condition for kinematic instability, whilst the maximum particle speeds amount |u|_max = 1.2–2c. The dI-determined directions of up-, down- and alongslope processes are all to within ±10° of the ADCP’s beam-spread averaged current (particle velocity) data.     

Dependence of geomagnetic activity during magnetic storms on the solar wind parameters for different types of streams

The dependence of the maximal values of the |Dst| and AE geomagnetic indices observed during magnetic storms on the value of the interplanetary electric field (E _y ) was studied based on the catalog of the large-scale solar wind types created using the OMNI database for 1976–2000 [Yermolaev et al., 2009]. An analysis was performed for eight categories of magnetic storms caused by different types of solar wind streams: corotating interaction regions (CIR, 86 storms); magnetic clouds (MC, 43); Sheath before MCs (Sh_MC, 8); Ejecta (95); Sheath (Sh_E, 56); all ICME events (MC + Ejecta, 138); all compression regions Sheaths before MCs and Ejecta (Sh_MC + Sh_E, 64); and an indeterminate type of storm (IND, 75). It was shown that the |Dst| index value increases with increasing electric field E _y for all eight types of streams. When electric fields are strong (E _y > 11 mV m⁻¹), the |Dst| index value becomes saturated within magnetic clouds MCs and possibly within all ICMEs (MC + Ejecta). The AE index value during magnetic storms is independent of the electric field value E _y for almost all streams except magnetic clouds MCs and possibly the compressed (Sheath) region before them (Sh_MC). The AE index linearly increases within MC at small values of the electric field (E _y < 11 mV m⁻¹) and decrease when these fields are strong (E _y > 11 mV m⁻¹). Since the dynamic pressure (Pd) and IMF fluctuations (σB) correlate with the E _y value in all solar wind types, both geomagnetic indices (|Dst| and AE) do not show an additional dependence on Pd and IMF δB. The nonlinear relationship between the intensities of the |Dst| and AE indices and the electric field E _y component, observed within MCs and possibly all ICMEs during strong electric fields E _y , agrees with modeling the magnetospheric-ionospheric current system of zone 1 under the conditions of the polar cap potential saturation.     

Finite amplitude holmboe waves

Abstract

We investigate the evolution of a parallel shear flow which has embedded within it a thin, symmetrically positioned layer of stable density stratification. The primary instability of this flow may deliver either Kelvin-Helmholtz waves or Holmboe waves, depending on the strength of the stratification. In this paper we describe a sequence of numerical simulations which reveal for the first time the behavior of the Holmboe wave at finite amplitude and clarify its structural relationship to the Kelvin-Helmholtz wave.

The flows investigated have initial profiles of horizontal velocity and Brunt-Vaisala frequency given in nondimensional form by U = tanhζ and N ²=J sech² RCζ, respectively, in which ζ is a nondimensional vertical coordinate, J is the value of the gradient Richardson number N ²/(dU/dζ)² at ζ=0, and R = 3. Linear stability theory predicts that the flow will develop Holmboe instability when J exceeds some critical value J_c' and Kelvin-Helmholtz instability when J is less than J_c; J_c being approximately equal to 0.25 when R=3. We simulate the evolution of flows with J=0.9, J=0.45, and J = 0.22, and find that the first two simulations yield Holmboe waves while the third yields a Kelvin-Helmholtz wave, as predicted.

The Holmboe wave is a superposition of two oppositely propagating disturbances, a right-going mode whose energy is concentrated in the region above the centre of the shear layer, and a left-going mode whose energy is concentrated below the centre of the shear layer. The horizontal speed of the modes varies periodically, and the variations are most pronounced at low values of J. If J ζ J_c' the minimum horizontal speed of the modes vanishes and the modes become phase-locked, whereupon they roll up to form a Kelvin-Helmholtz wave as predicted by Holmboe (1962). When J is moderately greater than J_c' the Holmboe wave ejects long, thin plumes of fluid into the regions above and below the shear layer, as has often been observed in laboratory experiments, and we examine in detail the mechanism by which this occurs.     

An antidynamo theorem for partly symmetric flows

Abstract

It is shown that magnetic fields generated by flows v _r,(r,t)e_r+v_T where v_T is an arbitrary toroidal component (e_r˙v_T≡V≡v_T≡0), cannot be maintained indefinitely against ohmic dissipation. The poloidal field variable max |r ² B _r| is shown to decay strictly monotonically with an undetermined decay rate. A bound on the growth of the toroidal field norm ∥T∥₁ is established solely dependent on the rate of conversion of poloidal to toroidal field, so that when the poloidal field is negligible then ∥T∥₁ decays strictly monotonically. The main application of these results is to models of stellar evolution based on axisymmetric differential rotation and spherically symmetric contraction. This symmetric velocity theorem overlaps with two already known theorems, namely the toroidal velocity theorem where v _r≡0 and the radial velocity theorem where v_T≡0. The new theorem does not entirely include the already established ones, principal differences being in the rates of decay and the field variables for which the decay is proven (see Table 1).     

Nonlinear dynamics of boussinesq convection in a deep rotating spherical shell-i

Abstract

Convection in a rotating spherical shell has wide application for understanding the dynamics of the atmospheres and interiors of many celestial bodies. In this paper we review linear results for convection in a shell of finite depth at substantial but not asymptotically large Taylor numbers, present nonlinear multimode calculations for similar conditions, and discuss the model and results in the context of the problem of solar convection and differential rotation. Detailed nonlinear calculations are presented for Taylor number T = 10⁵, Prandtl number P = 1, and Rayleigh number R between 1 |MX 10⁴ and 4 |MX 10⁴ (which is between about 4 and 16 times critical) for a shell of depth 20% of the outer radius. Sixteen longitudinal wave numbers are usually included (all even wave numbers m between 0 and 30) the amplitudes of which are computed on a staggered grid in the meridian plane.

The kinetic energy spectrum shows a peak in the wave number range m = 12–18 at R = 10⁴, which straddles the critical wave number m = 14 predicted by linear theory. These are modes which peak near the equator. The spectrum shows a second strong peak at m = 0, which represents the differential rotation driven by the peak convective modes. As R is increased, the amplitude of low wave numbers increases relative to high wave numbers as convection fills in in high and middle latitudes, and as the longitudinal scale of equatorial convection grows. By R = 3 |MX 10⁴, m = 8 is the peak convective mode. There is a clear minimum in the total kinetic energy at middle latitudes relative to low and high, well into the nonlinear regime, representing the continued dominance of equatorial and polar modes found in the linear case. The kinetic energy spectrum for m > 0 is maintained primarily by buoyancy work in each mode, but with substantial nonlinear transfer of kinetic energy from the peak modes to both lower and higher wave numbers.

For R = 1 to 2 |MX 10⁴, the differential rotation takes the form of an equatorial acceleration, with angular velocity generally decreasing with latitude away from the equator (as on the sun) and decreasing inwards. By R = 4 |MX 10⁴, this equatorial profile has completely reversed, with angular velocity increasing with depth and latitude. Also, a polar vortex which has positive rotation relative to the reference frame (no evidence of which has been seen on the sun) builds up as soon as polar modes become important. Meridional circulation is quite weak relative to differential rotation at R = 10⁴, but grows relative to it as R is increased. This circulation takes the farm of a single cell of large latitudinal extent in equatorial regions, with upward flow near the equator, together with a series of narrower cells in high latitudes. It is maintained primarily by axisymmetric buoyancy forces. The differential rotation is maintained at all R primarily by Reynolds stresses, rather than meridional circulation. Angular momentum transport toward the equator for R = 1–2 |MX 10⁴ maintains the equatorial acceleration while radially inward transport maintains the opposite profile at R = 4 |MX 10⁴.

The total heat flux out the top of the convective shell always shows two peaks for the range of R studied, one at the equator and the other near the poles (no significant variation with latitude is seen on the sun), while heat flux in at the bottom shows only a polar peak at large R. The meridional circulation and convective cells transport heat toward the equator to maintain this difference.

The helicity of the convection plus the differential rotation produced by it suggest the system may be capable of driving a field reversing dynamo, but the toroidal field may migrate with lime in each cycle toward the poles and equator, rather than just toward the equator as apparently occurs on the sun.

We finally outline additions to the physics of the model to make it more realistic for solar application.     

Patterns of thermal convection in slowly rotating,weakly magnetic fluid shells

The effects of rotation and a toroidal magnetic field on the preferred pattern of small amplitude convection in spherical fluid shells are considered. The convective motions are described in terms of associated Legendre functions P_l^|m| (cos θ). For a given pair of Prandtl number P and magnetic Prandtl number $\text{P}$_m the physically realized solution is represented either by m = 0 or |m| = l depending on the ratio of the rotation rate Λ to the magnetic field amplitude H. The case of m = 0 is preferred if this ratio ranges below a critical value, which is a function of the shell thickness, and |m| = l otherwise.     

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.     

Forces on spheres moving horizontally in a rotating stratified fluid

Abstract

Measurements have been made of the net horizontal force F acting on a sphere moving with horizontal velocity U (Reynolds numbers in the range 10²-10⁴) through a stratified fluid rotating about a vertical axis with uniform angular velocity Ω. In both homogeneous and stratified rotating fluids with small Rossby number R(R = U/Ωa ? 1 where a is the radius of the sphere) the force F is of magnitude 2ΩρUV (where ρ is the density of the fluid and V is the volume of the sphere). In a homogeneous fluid the relative directions of F and U were found to depend on the quantity F = 8Ωa ²/UD (where D is the depth of the fluid in which the object is placed (Mason, 1975)). In a rotating stratified fluid the relative directions of F and U are found to depend on the inverse Froude number k(k = Na/U where N ² = (g/δ)?ρ/?z) provided D > 4aΩ/N. In a homogeneous fluid with F ? 1 the force F is mainly in the U direction (a drag force due to inertial wave radiation) and is ～ ?0.4 |MX 2ΩρUV For F ? 1 a “Taylor column” occurs and the force, in correspondence with theoretical expectations, is ～ - 2Ω |MX UρV In a rotating stratified fluid with N ～2Ω and k ? 1 the force F is mainly in the U direction but is roughly one half of that occurring in the homogeneous situation with F ? 1 (tentatively explained as due to the evanescence of inertia-gravity disturbances). In a rotating stratified fluid with k ? 1 the flow should have no vertical motion (as with F ? 1) and again in correspondence with theoretical expectations the drag is ～ ?2 Ω |MX UρV. In a non-rotating stratified fluid the drag coefficient C _D(C _D = F U/½?ρU ²) was measured in the range k = 0.1 to 10 and had a maximum value ～ 1.2 for k ～ 3.     

Farlie-Gumbel-Morgenstern bivariate densities: Are they applicable in hydrology?

Certain bivariate densities constructed from marginals have recently been suggested as models of hydrologic variates such as rainfall intensity and depth. It is pointed out that (i) these densities belong to the families of the Farlie-Gumbel-Morgenstern densities and the Farlie polynomial densities, which have been extensively studied in the statistical literature, and that (ii) these densities have a limited potential applicability in hydrology since they can model only weakly associated variates, whose product-moment correlationR is within the range |R|1/3, under the first family of densities, and |R|1/2 under the second family.     

Farlie-Gumbel-Morgenstern bivariate densities: Are they applicable in hydrology?

Certain bivariate densities constructed from marginals have recently been suggested as models of hydrologic variates such as rainfall intensity and depth. It is pointed out that (i) these densities belong to the families of the Farlie-Gumbel-Morgenstern densities and the Farlie polynomial densities, which have been extensively studied in the statistical literature, and that (ii) these densities have a limited potential applicability in hydrology since they can model only weakly associated variates, whose product-moment correlationR is within the range |R|1/3, under the first family of densities, and |R|1/2 under the second family.     

On the frontal condition for finite amplitude α2Ω-dynamo wave trains in stellar shells

Abstract

In a previous paper, Bassom et al. (Proc. R. Soc. Lond. A, 455, 1443–1481, 1999) (BKS) investigated finite amplitude αΩ-dynamo wave trains in a thin turbulent, differentially rotating convective stellar shell; nonlinearity arose from α-quenching. There asymptotic solutions were developed based upon the small aspect ratio ε of the shell. Specifically, as a consequence of a prescribed latitudinally dependent α-effect and zonal shear flow, the wave trains have smooth amplitude modulation but are terminated abruptly across a front at some high latitude θ_F. Generally, the linear WKB-solution ahead of the front is characterised by the vanishing of the complex group velocity at a nearby point θ_f; this is essentially the Dee–Langer criterion, which determines both the wave frequency and front location.

Recently, Griffiths et al. (Geophys. Astrophys. Fluid Dynam. 94, 85–133, 2001) (GBSK) obtained solutions to the α²Ω-extension of the model by application of the Dee—Langer criterion. Its justification depends on the linear solution in a narrow layer ahead of the front on the short O(θ_f—θ_F) length scale; here conventional WKB-theory, used to describe the solution elsewhere, is inadequate because of mode coalescence. This becomes a highly sensitive issue, when considering the transition from the linear solution, which occurs when the dynamo number D takes its critical value D _c corresponding to the onset of kinematic dynamo action, to the fully nonlinear solutions, for which the Dee—Langer criterion pertains.

In this paper we investigate the nature of the narrow layer for α²Ω-dynamos in the limit of relatively small but finite α-effect Reynolds numbers R _α, explicitly ε^½ ? R ² _α ? 1. Though there is a multiplicity of solutions, our results show that the space occupied by the corresponding wave train is generally maximised by a solution with θ_f—θ_F small; such solutions are preferred as evinced by numerical simulations. This feature justifies the application by GBSK of the Dee—Langer criterion for all D down to the minimum D _min that the condition admits. Significantly, the frontal solutions are subcritical in the sense that |D _min| ≤ |D _c|; equality occurs as the α-effect Reynolds number tends to zero. We demonstrate that the critical linear solution is not connected by any parameter track to the preferred nonlinear solution associated with D _min. By implication, a complicated bifurcation sequence is required to make the connection between the linear and nonlinear states. This feature is in stark contrast to the corresponding results for αΩ-dynamos obtained by BKS valid in the limit R _{2 α} ? ε^½, which, though exhibiting a weak subcriticality, showed that the connection follows a clearly identifiable nonbifurcating track.     

Development of a robust runoff-prediction model by fusing the Rational Equation and a modified SCS-CN method

Abstract

The objective of this study is to develop a Modified Rational Equation (MoRE) that combines the advantages of the Rational Equation (e.g. simplicity and global acceptance) and those of the standard US Department of Agriculture (USDA) Soil Conservation Service (SCS) curve number (CN) method (e.g. easy parameterization and extensive verification across the world). Herein, the hypothesis is that the MoRE is more accurate, consistent and robust than the SCS-CN method and its improved versions in predicting runoff in watersheds with limited data. The MoRE was designed to have a simple structure that is described by four intrinsic parameters: CN, permanent wilting point, field capacity and saturation soil moisture, and does not include initial abstraction as a variable. An evaluation of 77 USDA small agricultural watersheds indicated that CN of the MoRE has different physical meanings from CN of the SCS-CN method. The MoRE (mean Nash-Sutcliffe coefficient, E > 0.73) performed better than the SCS-CN (mean E < 0.32) and the four improved models (mean E < 0.56) in reproducing the runoff of the study watersheds. Performance of all six models varied greatly between watersheds, as well as between events, but was independent of watershed drainage area. However, the model performances tend to be better for watersheds and/or events with a runoff-to-rainfall ratio of between 0.1 and 0.3 than for those with a ratio outside this range. The MoRE has the most consistent and robust performance.

Editor D. Koutsoyiannis; Associate editor I. Nalbantis

Citation Wang, X., Liu, T., and Yang, W., 2012. Development of a robust runoff-prediction model by fusing the rational equation and a modified SCS-CN method. Hydrological Sciences Journal, 57 (6), 1118–1140.     

Relativistic magnetohydrodynamic winds of finite temperature

Abstract

The singular differential equations for finite temperature relativistic magnetohydrodynamic (MHD) winds integrate to two algebraic equations when the source magnetic field is a monopole. This simplification enables an extensive characterization of the asymptotic wind solutions in terms of source parameters. We will consider only the critical solutions-those that pass smoothly through both an intermediate (Alfvenic) and a fast MHD critical point and expand to zero pressure at infinite radial distance from the source. Because the constants of motion must be specified to extremely high accuracy, the critical solutions cannot be found analytically. Synopsis of many numerical solutions reveals a uniform parametric characterization of the asymptotic wind in terms of one combination of source parameters, Z, the mean source particle energy divided by mc²[sgrave]^½, where [sgrave] is a generalization of Michel's (1969) cold relativistic wind strength parameter. Cool winds, with Z<1, behave asymptotically much as Michel's cold wind minimum torque solution; Z1 hot winds have quite different, but simply characterized, asymptotic solutions. Thus, the strength of magnetized relativistic outflows can depend critically upon the temperature of the source.     

Topological invariants of field lines rooted to planes

Abstract

Two open curves with fixed endpoints on a boundary surface can be topologically linked. However, the Gauss linkage integral applies only to closed curves and cannot measure their linkage. Here we employ the concept of relative helicity in order to define a linkage for open curves. For a magnetic field consisting of closed field lines, the magnetic helicity integral can be expressed as the sum of Gauss linkage integrals over pairs of lines. Relative helicity extends the helicity integral to volumes where field lines may cross the boundary surface. By analogy, linkages can be defined for open lines by requiring that their sum equal the relative helicity.

With this definition, the linkage of two lines which extend between two parallel planes simply equals the number of turns the lines take about each other. We obtain this result by first defining a gauge-invariant, one-dimensional helicity density, i.e. the relative helicity of an infinitesimally thin plane slab. This quantity has a physical interpretation in terms of the rate at which field lines lines wind about each other in the direction normal to the plane. A different method is employed for lines with both endpoints on one plane; this method expresses linkages in terms of a certain Gauss linkage integral plus a correction term. In general, the linkage number of two curves can be put in the form L=r + n, |r|≦1J2, where r depends only on the positions of the endpoints, and n is an integer which reflects the order of braiding of the curves.

Given fixed endpoints, the linkage numbers of a magnetic field are ideal magneto-hydrodynamic invariants. These numbers may be useful in the analysis of magnetic structures not bounded by magnetic surfaces, for example solar coronal fields rooted in the photosphere. Unfortunately, the set of linkage numbers for a field does not uniquely determine the field line topology. We briefly discuss the problem of providing a complete and economical classification of field topologies, using concepts from the theory of braid equivalence classes.     

Paradoxes of geoelectrics

The problem of energy paradoxes revealed in geoelectrics are discussed. The experimental facts illustrating the anomalous energy characteristics of the magnetotelluric (MT) field are presented. An attempt is made to interpret these anomalies from the standpoint of directional analysis. Two three-layer models corresponding to the situation |Q| > 1 and \(\widetilde {{S_z}} < 0\) are found by the numerical modeling. The possibility of accounting for the observed paradoxes within the resonance model “heterogeneous plane wave—layered medium” is discussed.     

Bare soil evaporation under high evaporation demand: a proposed modification to the FAO-56 model

Abstract

This study evaluates the evaporation component of the FAO-56 model under high evaporation demand. To perform this, two data sets were used as field evaluation, and a second model was used for comparison (a model based on the square root of time, SRT). The results show that although FAO-56, the field data and the SRT model present similar cumulative evaporation over the study period (approximately one month), when the data are analysed daily, FAO-56 overestimated evaporation at the beginning of the process and underestimated it at the end. A correction for FAO-56 is proposed to amend the mismatch between FAO-56 and the field-measured data under high evaporation conditions. Consequently, the parameters used by the FAO-56 evaporation component are discussed.

Citation Torres, E. A. & Calera, A. (2010) Bare soil evaporation under high evaporation demand: a proposed modification to the FAO-56 model. Hydrol. Sci. J. 55(3), 303–315.     

Nonlinear development of phase-mixed alfvén waves

Abstract

We derive an equation governing the nonlinear propagation of a linearly polarized Alfvén wave in a two-dimensional, anisotropic, slightly compressible, highly magnetized, viscous plasma, where nonlinearities arise from the interaction of the Alfvén wave with fast and slow magnetoacoustic waves. The phase mixing of such a wave has been suggested as a mechanism for heating the outer solar atmosphere (Heyvaerts and Priest, 1983).

We find that cubic wave damping dominates shear linear dissipation whenever the Alfvén wave velocity amplitude δv_y exceeds a few times ten metres per second. In the nonlinear regime, phase-mixed waves are marginally stable, while non-phase-mixed waves of wavenumber k_a are damped over a timescale k_uRe ₀|δ v_y/v_A |^?2, Re ₀ being the Reynolds number corresponding to the Braginskij viscosity coefficient η₀ and v_A the Alfvén speed. Dissipation is most effective where β = (v_s /v_A) ² ≈ 1, v_s being the speed of sound.