Spontaneous tangential discontinuities and the optical analogy for static magnetic fields. VI. Topology of current sheets

Abstract

The electric surface current in a tangential discontinuity in a force-free magnetic field is conserved. The direction of the current is halfway between the direction of the continuous fields on either side of the surface of discontinuity. Hence the current sheets, i.e. the surface of tangential discontinuity, have a topology that is distinct from the lines of force of the field. The precise nature of the topology of the current sheet depends upon the form of the winding patterns in the field. Hence, invariant winding patterns and random winding patterns are treated separately. Current sheets may have edges, at the junction of two or more topological separatrices. The current lines may, in special cases, be closed on themselves. The lines of force that lie on either side of a current sheet somewhere pass off the sheet across a junction onto another sheet. In most cases the current sheets extending along a field make an irregular honeycomb.

The honeycomb pattern varies along the field if the winding pattern of the field varies. The surface current density in a tangential discontinuity declines inversely, or faster, with distance from its region of origin. The edges of weaker tangential discontinuities (originating in more distant regions) are bounded by the stronger tangential discontinuities (of nearby origin).

An examination of the force-free field equations in a small neighborhood of the line of intersection of two tangential discontinuities shows that the lines of force twist around to cross the line of intersection at right angles. If the angle between the tangential discontinuities exceeds π/2, there is also the possibilitity that the lines twist around so as to come tangent to the line of intersection as they cross it.     

Tangential discontinuities and the optical analogy for stationary fields. VII. Effects of resistivity

Abstract

The relatively large resistivity in the solar photosphere and chromosphere softens the ideal tangential discontinuities of magnetostatic equilibrium into continuous transitions in field direction over scales of 0.1–10 km. This softening is communicated upward at the Alfvén speed into the active solar corona. The degree of softening is a vital part of the theory of magnetic heat input to the active X-ray corona, because the very low resistivity of the coronal gas provides effective dissipation only if the current sheets are reduced to a thickness of 10^?2km.

A close examination of the problem shows that the Alfvén transit time up into the corona is large compared to the characteristic time of 1 sec in which the coronal tangential discontinutities are formed. It also shows that the principal effect of the resistivity is to create a thin surface layer of fluid on adjacent flux bundles, which causes a general drift of the flux but does not directly broaden the current sheets higher up in the field. In fact the motions of the surface layers do not extend upwards beyond the first winding pattern at each end of a coronal loop.

It appears that the photospheric and chromospheric resistivity is without striking consequences for magnetic heating in the corona.     

Global aspects of solar wind–ionosphere interactions

Recent observations have quantified the auroral wind O⁺ outflow in response to magnetospheric inputs to the ionosphere, notably Poynting energy flux and precipitating electron density. For moderate to high activity periods, ionospheric O⁺ is observed to become a significant or dominant component of plasma pressure in the inner plasma sheet and ring current regions. Using a global circulation model of magnetospheric fields and its imposed ionospheric boundary conditions, we evaluate the global ionospheric plasma response to local magnetospheric conditions imposed by the simulation and evaluate magnetospheric circulation of solar wind H⁺, polar wind H⁺, and auroral wind O⁺. We launch and track the motions of millions of test particles in the global fields, launched at randomly distributed positions and times. Each particle is launched with a flux weighting and perpendicular and parallel energies randomly selected from defined thermal ranges appropriate to the launch point. One sequence is driven by a two-hour period of southward interplanetary magnetic field for average solar wind intensity. A second is driven by a 2-h period of enhanced solar wind dynamic pressure for average interplanetary field. We find that the simulated ionospheric O⁺ becomes a significant plasma pressure component in the inner plasma sheet and outer ring current region, particularly when the solar wind is intense or its magnetic field is southward directed. We infer that the reported empirical scalings of auroral wind O⁺ outflows are consistent with a substantial pressure contribution to the inner plasma sheet and plasma source surrounding the ring current. This result violates the common assumption that the ionospheric load is entirely confined to the F layer, and shows that the ionosphere is often an important dynamic element throughout the magnetosphere during moderate to large solar wind disturbances.     

Studying ice sheet stability using the method of separation of variables

Abstract

Stability of two-dimensional stationary continental and marine ice sheets is studied using perturbations of ice sheet surface elevation and the margin position about a steady-state solution. Ice accumulation rate at the upper surface is specified as a function of elevation and span. Ice is considered as a Newtonian fluid. Linearisation and separation of variables yield a non-homogeneous eigenvalue problem. In case when the margin moves, a necessary condition for the existence of a solution is very restrictive on the functions of the bed profile and accumulation rate. The method of separation of variables is applicable when the margin is assumed to be stationary. Singularity of the perturbation at the margin is treated using the method of matched asymptotic expansions. Numerical experiments show that dependence of accumulation rate on elevation contributes strongly towards ice sheet instability. Bed slope, ice thickness at the grounding line of a marine ice sheet and equilibrium line inclination are the main parameters determining the ice sheet reaction to a surface perturbation.     

Pore-fluid pressures and frictional heating on a fault surface

This study considers the effects of heat transfer and fluid flow on the thernal, hydrologic, and mechanical response of a fault surface during seismic failure. Numerical modeling techniques are used to account for the coupling of the thermal, fluid-pressure, and stress fields. Results indicate that during an earthquake the failure surface is heated to a tempeature required for the thermal expansion of pore fluids to balance the rate of fluid loss due to flow and the fluid-volume changes due to pore dilatation. Once this condition is established, the pore fluids pressurize and the shear strength decreases rapidly to a value sufficient to maintain the thermal pressurization of pore fluids at near-lithostatic values. If the initial fluid pressure is hydrostatic, the final temperature attained on the failure surface will increase with depth, because a greater pressure increase can occur before a near-lithostatic pressure is reached. The rate at which thermal pressurization proceeds depends primarily on the hydraulic characteristics of the surrounding porous medium, the coefficient of friction on the fault surface, and the slip velocity. If either the permeability exceeds 10^–15 m² or the porous medium compressibility exceeds 10^–8 Pa^–1, then frictional melting may occur on the fault surface before thermal pressurization becomes significant. If the coefficient of friction is less than 10^–1 and if the slip velocity is less than 10^–2 msec^–1, then it is doubtful that either thermal pressurization or frictional melting on the fault surface could cause a reduction in the dynamic shear strength of a fault during an earthquake event.     

Magnetohydrodynamic simulation of a solar flare: 2. Flare model and simulation using active-region magnetic maps

The results of a three-dimensional MHD simulation and data obtained using specialized spacecraft made it possible to construct an electrodynamic model of solar flares. A flare results from explosive magnetic reconnection in a current sheet above an active region, and electrons accelerated in field-aligned currents cause hard X rays on the solar surface. In this review, we considered works where the boundary and initial conditions on the photosphere were specified directly from the magnetic maps, obtained by SOHO MDI in the preflare state, in order to simulate the formation of a current sheet. A numerical solution of the complete set of MHD equations, performed using the new-generation PERESVET program, demonstrated the formation of several current sheets before a series of flares. A comparison of the observed relativistic proton spectra and the simulated proton acceleration along a magnetic field singular line made it possible to estimate the magnetic reconnection rate during a flare (∼10⁷ cm s⁻¹). Great flares (of the X class) originate after an increase in the active region magnetic flux up to 10²² Mx.     

The kinematics of hexagonal magnetoconvection

Abstract

Calculations are presented for the evolution of a magnetic field which is subject to the effect of three-dimensional motions in a convecting layer of highly conducting fluid with hexagonal symmetry. The back reaction of the field on the motions via the Lorentz force is neglected. We consider cases where the imposed field is either vertical or horizontal. In the former case, flux accumulates at cell centres, with subsidiary concentrations at the vertices of the pattern. In the latter, topological asymmetries between up- and down-moving fluid regions generate positive flux at the base of the layer and negative flux at the top, though the system is actually an amplifier rather than a self-excited dynamo. Spiral field lines form in the interiors of the cells, and the phenomenon of “flux expulsion” found in two-dimensional solutions is somewhat altered when the imposed field is horizontal. Applications for stellar magnetic fields include a possible mechanism for burying flux at the base of a convection zone.     

Thermal convection in differentially rotating systems

Abstract

The onset of convection in a cylindrical fluid annulus is analyzed in the case when the cylindrical walls are rotating differentially, a temperature gradient in the radial direction is applied, and the centrifugal force dominates over gravity. The small gap approximation is used and no-slip conditions on the cylindrical walls are assumed. It is found that over a considerable range of the parameter space either convection rolls aligned with the axis of rotation or rolls in the perpendicular (azimuthal) direction are preferred. It is shown that by a suitable redefinition of parameters, results for finite amplitude Taylor vortices and for convection rolls in the presence of shear can be applied to the present problem. Weakly nonlinear results for transverse rolls in a Couette flow indicate the possibility of subcritical bifurcation for Prandtl numbers P less than 0.82. Heat and momentum transports are derived as functions of P and the problem of interaction between transverse and longitudinal rolls is considered. The relevance of the analysis for problems of convection in planetary and stellar atmospheres is briefly discussed.     

Instability of a geostrophic front and its energetics

Abstract

The instability of a current with a geostrophic surface density front is investigated by means of a reduced gravity model having a velocity profile with nearly uniform potential vorticity. It is shown that currents are unstable when the mean potential vorticity decreases toward the surface front at the critical point of the frontal trapped waves investigated by Paldor (1983). This instability is identical with that demonstrated by Killworth (1983) in the longwave limit.

The cross-stream component of mass flux and the rates of energy conversions among the five energy forms defined by Orlanski (1968) are also calculated. The main results are as follows, (a) The mass flux toward the surface front is positive near the front and negative around the critical point. The positive mass flux near the front does not vanish at the position of the undisturbed surface front, so that the mean position of the front moves outward and the region of the strong current spreads. (b) The potential energy of the mean flow integrated over the fluid is released through the work done by the force of the pressure gradient of the mean flow on the fluid, and is converted into the kinetic energy of the mean flow. (c) In the critical layer, the mean flow is rapidly accelerated with the growth of the unstable wave. This acceleration is caused by the rapid phase shift of the unstable wave in the critical layer.     

Two-phase hydrodynamic and sediment transport modeling of wave-generated sheet flow

This numerical investigation was carried out to advance mechanistic understanding of sediment transport under sheet flow conditions. An Euler–Euler coupled two-phase flow model was developed to simulate fluid–sediment oscillatory sheet flow. Since the concentration of sediment particles is high in such flows, the kinematics of the fluid and sediment phases are strongly coupled. This model includes interaction forces, intergranular stresses and turbulent stress closure. Each phase was modeled via the Reynolds-Averaged Navier–Stokes equations, with interphase momentum conservation accounting for the interaction between the phases. The generation and transformation of turbulence was modeled using the two-equation k–ε$k''–''\epsilon$ turbulence model. Concentration and sediment flux profiles were compared with experimental data for sheet flow conditions considering both symmetric and asymmetric oscillatory flows. Sediment and fluid velocity variations, concentration profiles, sediment flux and turbulence parameters of wave-generated sheet flow were studied numerically with a focus on sediment transport characteristics. In all applications, the model predictions compared well with the experimental data. Unlike previous investigations in which the flow is driven by a horizontal pressure gradient, the present model solves the Navier–Stokes equations under propagating waves. The model’s ability to predict sediment transport under oscillatory sheet flow conditions underscores its potential for understanding the evolution of beach morphology.     

Thermomechanical balances of ice sheet flows

Abstract

The flow of large natural ice masses under gravity is described by the mass, momentum, and energy balances of an incompressible, homogeneous, heat conducting, non-linearly viscous fluid in which the shear response includes a strongly temperature-dependent rate factor. Dimensionless analysis and co-ordinate stretching reflecting the long aspect ratio show that series expansions in a small parameter which determines the surface slope magnitude are uniformly valid even when temperature variation induces a strongly non-uniform mechanical response. The normalised energy balance shows that both horizontal and vertical advection are significant in thin and thick grounded sheets and in floating shelves, and that viscous dissipation can be significant in basal regions of a grounded sheet, and hence there is strong thermomechanical coupling. Moreover, though a thermal basal boundary layer may arise in a thick sheet, it would only give rise to significantly enhanced temperature and strain-rate gradients in extreme circumstances. The leading order relations for steady plane flow of a grounded sheet are reduced to a parabolic system for the temperature and two velocity components, which incorporates the unknown surface slope in coefficients and boundary conditions. This provides a useful starting point for numerical solution of the thermomechanically coupled problem. A fixed domain mapping is presented as an attractive alternative formulation when the bed topography is close to planar.     

Distribution of magnetic energy in αΩ-dynamos,I: The method

Abstract

In this paper a method for solving the equation for the mean magnetic energy <BB> of a solar type dynamo with an axisymmetric convection zone geometry is developed and the main features of the method are described. This method is referred to as the finite magnetic energy method since it is based on the idea that the real magnetic field B of the dynamo remains finite only if <BB> remains finite. Ensemble averaging is used, which implies that fields of all spatial scales are included, small-scale as well as large-scale fields. The method yields an energy balance for the mean energy density ε ≡ B ²/8π of the dynamo, from which the relative energy production rates by the different dynamo processes can be inferred. An estimate for the r.m.s. field strength at the surface and at the base of the convection zone can be found by comparing the magnetic energy density and the outgoing flux at the surface with the observed values. We neglect resistive effects and present arguments indicating that this is a fair assumption for the solar convection zone. The model considerations and examples presented indicate that (1) the energy loss at the solar surface is almost instantaneous; (2) the convection in the convection zone takes place in the form of giant cells; (3) the r.m.s. field strength at the base of the solar convection zone is no more than a few hundred gauss; (4) the turbulent diffusion coefficient within the bulk of the convection zone is about 10¹⁴cm²s^?1, which is an order of magnitude larger than usually adopted in solar mean field models.     

Landscape preservation under Fennoscandian ice sheets determined from in situ produced Be and Al

Some areas within ice sheet boundaries retain pre-existing landforms and thus either remained as ice free islands (nunataks) during glaciation, or were preserved under ice. Differentiating between these alternatives has significant implications for paleoenvironment, ice sheet surface elevation, and ice volume reconstructions. In the northern Swedish mountains, in situ cosmogenic ¹⁰Be and ²⁶Al concentrations from glacial erratics on relict surfaces as well as glacially eroded bedrock adjacent to these surfaces, provide consistent last deglaciation exposure ages (∼8-13 kyr), confirming ice sheet overriding as opposed to ice free conditions. However, these ages contrast with exposure ages of 34-61 kyr on bedrock surfaces in these same relict areas, demonstrating that relict areas were preserved with little erosion through multiple glacial cycles. Based on the difference in radioactive decay between ²⁶Al and ¹⁰Be, the measured nuclide concentration in one of these bedrock surfaces suggests that it remained largely unmodified for a minimum period of 845₋₄₁₈⁺⁴⁶¹ kyr. These results indicate that relict areas need to be accounted for as frozen bed patches in basal boundary conditions for ice sheet models, and in landscape development models. Subglacial preservation also implies that source areas for glacial sediments in ocean cores are considerably smaller than the total area covered by ice sheets. These relict areas also have significance as potential long-term subglacial biologic refugia.     

Effect of a uniform magnetic field on nonlinear magnetocenvection in a rotating fluid spherical shell

Abstract

Dynamic interaction between magnetic field and fluid motion is studied through a numerical experiment of nonlinear three-dimensional magnetoconvection in a rapidly rotating spherical fluid shell to which a uniform magnetic field parallel to its spin axis is applied. The fluid shell is heated by internal heat sources to maintain thermal convection. The mean value of the magnetic Reynolds number in the fluid shell is 22.4 and 10 pairs of axially aligned vortex rolls are stably developed. We found that confinement of magnetic flux into anti-cyclonic vortex rolls was crucial on an abrupt change of the mode of magnetoconvection which occurred at Δ = 1 ～ 2, where A is the Elsasser number. After the mode change, the fluid shell can store a large amount of magnetic flux in itself by changing its convection style, and the magnetostrophic balance among the Coriolis, Lorentz and pressure forces is established. Furthermore, the toroidal/poloidal ratio of the induced magnetic energy becomes less than unity, and the magnetized anti-cyclones are enlarged due to the effect of the magnetic force. Using these key ideas, we investigated the causes of the mode change of magnetoconvection. Considering relatively large magnetic Reynolds number and a rapid rotation rate of this model, we believe that these basic ideas used to interpret the present numerical experiment can be applied to the dynamics in the Earth's and other planetary cores.     

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.     

Sausage mode instability of thin current sheets as a cause of magnetospheric substorms

Observations have shown that, prior to substorm explosions, thin current sheets are formed in the plasma sheet of the Earth’s magnetotail. This provokes the question, to what extent current-sheet thinning and substorm onsets are physically, maybe even causally, related. To answer this question, one has to understand the plasma stability of thin current sheets. Kinetic effects must be taken into account since particle scales are reached in the course of tail current-sheet thinning. We present the results of theoretical investigations of the stability of thin current sheets and about the most unstable mode of their decay. Our conclusions are based upon a non-local linear dispersion analysis of a cross-magnetic field instability of Harris-type current sheets. We found that a sausage-mode bulk current instability starts after a sheet has thinned down to the ion inertial length. We also present the results of three-dimensional electromagnetic PIC-code simulations carried out for mass ratios up to M_i/m_e = 64. They verify the linearly predicted properties of the sausage mode decay of thin current sheets in the parameter range of interest.     

Tangential discontinuities and the optical analogy for stationary fields. v. formal integration of the force-free field equations

Abstract

This paper demonstrates the appearance of tangential discontinuities in deformed force-free fields by direct integration of the field equation ? x B = αB. To keep the mathematics tractable the initial field is chosen to be a layer of linear force-free field B_x = + B ₀cosqz, B_y = — B ₀sinqz, B_z = 0, anchored at the distant cylindrical surface ? = (x ² + y ²)^1/2 = R and deformed by application of a local pressure maximum of scale l centered on the origin x = y = 0. In the limit of large R/l the deformed field remains linear, with α = q[1 + O(l ²/R ²)]. The field equations can be integrated over ? = R showing a discontinuity extending along the lines of force crossing the pessure maximum. On the other hand, examination of the continuous solutions to the field equations shows that specification of the normal component on the enclosing boundary ? = R completely determines the connectivity throughout the region, in a form unlike the straight across connections of the initial field. The field can escape this restriction only by developing internal discontinuities.

Casting the field equation in a form that the connectivity can be specified explicitly, reduces the field equation to the eikonal equation, describing the optical analogy, treated in papers II and III of this series. This demonstrates the ubiquitous nature of the tangential discontinuity in a force-free field subject to any local deformation.     

Large aspect ratio cells in two-dimensional thermal convection

Numerical experiments have been carried out on two-dimensional thermal convection, in a Boussinesq fluid with infinite Prandtl number, at high Rayleigh numbers. With stress free boundary conditions and fixed heat flux on upper and lower boundaries, convection cells develop with aspect ratios (width/depth) λ? 5, if heat is supplied either entirely from within or entirely from below the fluid layer. The preferred aspect ratio is affected by the lateral boundary conditions. If the temperature, rather than the heat flux, is fixed on the upper boundary the cells haveλ ≈ 1. At Rayleigh numbers of 2.4 × 10⁵ and greater, small sinking sheets are superimposed on the large aspect ratio cells, though they do not disrupt the circulation. Similar two-scale flows have been proposed for convection in the earth's mantle. The existence of two scales of flow in two-dimensional numerical experiments when the viscosity is constant will allow a variety of geophysically important effects to be investigated.     

鄱阳湖夏季水热通量特征及环境要素影响分析

气候变化加速了全球水文循环过程,然而,气候变化如何影响水体蒸发及其水热通量交换仍然不清楚.基于涡度相关系统观测鄱阳湖水体水热通量过程,在小时和日尺度分析了水热通量的变化规律及其主要影响因子.研究表明,潜热通量日变化波动剧烈,大部分为正值,变化范围在-50~580 W/m2之间.而感热通量数值较小,变化范围在-50~50 W/m2之间.8月份潜热通量和感热通量均呈波动下降趋势,均值分别为167.4和15.9 W/m2.8月份日平均潜热通量和感热通量之和大于净辐射,这是由于这一时段储存在水体中的热量释放并补充潜热通量和感热通量.小时尺度上潜热通量日变化在相位上与净辐射无显著相关性,而与风速显著相关.在日尺度变化趋势上,8月份日平均潜热通量仍主要受到风速和水温的影响,感热通量则主要受到风速和饱和水汽压差的影响.