Ground magnetic characteristics of the storm-time ring current Asymmetry

Using the minute data of the H component of geomagnetic field from the 20°E magnetic meridian chain and the 30°N magnetic latitudinal chain, the temporal evolution characteristics of the equatorial ring current during the storm on November 7-10, 2004 are studied. It is indicated that the UT-MLT and UT-MLAT graphics extremely exhibit the local time distribution, latitudinal variation and temporal evo- lution of the H component. The results show: (1) The UT-MLT contour clearly shows the increasing of the H component mostly around noon during the initial phase, representing the geomagnetic effect from the magnetopause current system. During the main phase, most negative values of the H com- ponent appear around the dusk-side, indicating the dawn-dusk asymmetric distribution of the ring cur- rent. (2) The contour of UT-MLAT suggests the latitudinal variation of the H component decreasing with the enhancement of the latitudes during geomagnetic storm, which is in good agreement with the Dst index. The latitudinal variations provide a new sight for describing the temporal characteristics of the intensity of the storm-time ring current. (3) Both the contours of UT-MLT and UT-MLAT are useful to monitor the space-time distribution of the equatorial ring current.     

Analysis of the heliospheric current sheet fine structure: Single or multiple current sheets

When studying the heliospheric current sheet (HCS) local structure, it is not unusual to find wide HCS crossings. In this paper, we present one crossing that appears to have a complex internal structure composed of three parallel sheets and several possible HCS crossings that are consecutive and are on the order of minutes. Depending on their origin, different scenarios can explain multiple current sheets such as complex structures of helmet streamer at the Corona flowing into the solar wind, local waviness in the HCS structure, local oscillations of the HCS, and inverted magnetic fields or planar magnetic structures (PMS) close to the HCS. Distinguishing among these scenarios using just one observational point is very difficult because all of them are 3D structures. Nevertheless, we think that nearly parallel sheets are more likely in the first and in the last scenarios, i.e. multiple helmet streamer structure and PMS. In order to make the distinction between them, we have studied the possible reversal in the Q_e·B sign for every event. Our results suggest that the existence of not-wide HCS composed of multiple parallel sheets cannot be rejected.     

Concerning the generation of geomagnetic giant pulsations by drift-bounce resonance ring current instabilities

Giant pulsations are nearly monochromatic ULF-pulsations of the Earth’s magnetic field with periods of about 100 s and amplitudes of up to 40 nT. For one such event ground-magnetic observations as well as simultaneous GEOS-2 magnetic and electric field data and proton flux measurements made in the geostationary orbit have been analysed. The observations of the electromagnetic field indicate the excitation of an odd-mode type fundamental field line oscillation. A clear correlation between variations of the proton flux in the energy range 30–90 keV with the giant pulsation event observed at the ground is found. Furthermore, the proton phase space density exhibits a bump-on-the-tail signature at about 60 keV. Assuming a drift-bounce resonance instability as a possible generation mechanism, the azimuthal wave number of the pulsation wave field may be determined using a generalized resonance condition. The value determined in this way, m = −21±4, is in accord with the value m = −27±6 determined from ground-magnetic measurements. A more detailed examination of the observed ring current plasma distribution function f shows that odd-mode type eigenoscillations are expected for the case ∂f/∂W ≥ 0, much as observed. This result is different from previous theoretical studies as we not only consider local gradients of the distribution function in real space, but also in velocity space. It is therefore concluded that the observed giant pulsation is the result of a drift-bounce resonance instability of the ring current plasma coupling to an odd-mode fundamental standing wave. The generation of the bump-on-the-tail distribution causing ≥f/≥W ≥ 0 can be explained due to velocity dispersion of protons injected into the ring current. Both this velocity dispersion and the necessary substorm activity causing the injection of protons into the nightside magnetosphere are observed.     

Effects of chemical reactions on density-dependent fluid flow: on the numerical formulation and the development of instabilities

A three-dimensional, reactive numerical flow model is developed that couples chemical reactions with density-dependent mass transport and fluid flow. The model includes equilibrium reactions for the aqueous species, kinetic reactions between the solid and aqueous phases, and full coupling of porosity and permeability changes that result from precipitation and dissolution reactions in porous media. A one-step, global implicit approach is used to solve the coupled flow, transport and reaction equations with a fully implicit upstream-weighted control volume discretization. The Newton–Raphson method is applied to the discretized non-linear equations and a block ILU-preconditioned CGSTAB method is used to solve the resulting Jacobian matrix equations. This approach permits the solution of the complete set of governing equations for both concentration and pressure simultaneously affected by chemical and physical processes. A series of chemical transport simulations are conducted to investigate coupled processes of reactive chemical transport and density-dependent flow and their subsequent impact on the development of preferential flow paths in porous media. The coupled effects of the processes driving flow and the chemical reactions occurring during solute transport is studied using a carbonate system in fully saturated porous media. Results demonstrate that instability development is sensitive to the initial perturbation caused by density differences between the solute plume and the ambient groundwater. If the initial perturbation is large, then it acts as a “trigger” in the flow system that causes instabilities to develop in a planar reaction front. When permeability changes occur due to dissolution reactions occurring in the porous media, a reactive feedback loop is created by calcite dissolution and the mixed convective transport of the system. Although the feedback loop does not have a significant impact on plume shape, complex concentration distributions develop as a result of the instabilities generated in the flow system.     

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.     

Multiple current sheets in a double auroral oval observed from the MAGION-2 and MAGION-3 satellites

A case is described of multiple current sheets crossed by the MAGION-2 satellite in the near-midnight quieting auroral oval. The data were obtained by the magnetometer experiment onboard. Results show during a quieting period after a preceding substorm, or during an early growth phase of the next substorm, two double-sheet current bands, POLE and EQUB, located at respectively the polar and equatorial borders of the auroral oval separated by about 500 km in latitude. This is consistent with the double-oval structure during recovery introduced by Elphinstone et al. (1995). Within the POLE, the magnetic field data show simultaneous existence of several narrow parallel bipolar current sheets within the upward current branch (at 69.5–70.3° invariant latitude) with an adjacent downward current branch at its polar side at (70.5–71.3°). The EQUB was similarly stratified and located at 61.2–63.5° invariant latitude. The narrow current sheets were separated on average by about 35 km and 15 km, respectively, within the POLE and EQUB. A similar case of double-oval current bands with small-scale structuring of their upward current branches during a quieting period is found in the data from the MAGION-3 satellite. These observations contribute to the double-oval structure of the late recovery phase, and add a small-scale structuring of the upward currents producing the auroral arcs in the double- oval pattern, at least for the cases presented here. Other observations of multiple auroral current sheets and theories of auroral arc multiplicity are briefly discussed. It is suggested that multiple X-lines in the distant tail, and/or leakage of energetic particles and FA currents from a series of plasmoids formed during preceding magnetic activity, could be one cause of highly stratified upward FA currents at the polar edge of the quieting double auroral oval.     

Numerical MHD simulations of the appearance of a series of current sheets above the active region AR 0365

The first attempt at numerical MHD simulations of the appearance of several current sheets above an active region before a series of elementary flares is described. Energy accumulates in the field of each sheet that can be released during one of the flares. The computations started three days before the appearance of a series of flares, i.e., before the emergence of a new magnetic flux in the active region. The initial (potential) magnetic field was calculated by solving the Laplace equation with an oblique derivative. The boundary conditions on the photosphere were specified from maps of the measured magnetic field in the active region for various instants of time. The Peresvet program solving the full system of MHD equations with dissipative terms was used in the computations. An absolutely implicit scheme conservative relative to the magnetic flux was used. The problem of properly choosing the size of the computational domain and finding the positions of singular magnetic field lines is discussed.     

Relative current effect on short wave growth

Ocean Dynamics - Short waves growth is characterized by nonlinear and dynamic processes that couple ocean and atmosphere. Ocean surface currents can have a strong impact on short wave steepness and...     

Asymmetry in the form of solar spots

We have analyzed the geometric characteristics of sunspots. The form of sunspots has been studied by sunspot image normalization to obtain the average profile of spots and the profile relative to the position of cores. The deviation of the sunspot form from the axisymmetric configuration has been studied. We have found that the spots of leading and trailing polarities have a drop shape. The cores of leading and trailing sunspots are shifted toward the western and eastern edges of the photosphere–penumbra boundary, respectively. The strength of the magnetic field of the cores of leading spots in the eastern hemisphere exceeds the field strength in the western hemisphere. We considered the tilt of the form of sunspots as a function of size. The form of spots of a large area (S > 1000 ppm of solar hemisphere) is elongated along the magnetic axis of the bipole of a group of sunspots.     

山体遮挡对滇池风生流的影响初探

用二维风生流数值模型模拟滇池湖流运动。滇池在主导风向西南风作用下,假定湖面风场是均匀的,数值模拟的湖流流态与实测湖流结果相差很大。而考虑山体遮挡影响,根据实测湖流期间现有的风情资料,在湖面上构造一非均匀风场,数值模拟结果与实测值基本一致。山体遮挡对滇池风生流的影响是不容忽视的。建议进一步进行湖流和湖面风向、风速监测,并建立过山气流数学模型,深入研究山体遮挡对湖泊风生流的影响。     

The effect of wind on the residual current velocities in the inlets of Venice lagoon

The residual flow in the inlets of Venice lagoon subject to Bora and Sirocco winds has been studied. Current velocities have been monitored since 2001 using Acoustic Doppler Current Profilers (ADCP) installed on the beds of the inlets that connect the lagoon to the Adriatic Sea; these inlets are Lido, Malamocco and Chioggia. Wind velocity data have also been continuously measured at an oceanographic platform 14.8 km offshore from the lagoon; these data were subsequently decomposed into Principal Components, which are associated with Bora and Sirocco wind directions. Analyses show that the inflow in Lido inlet is strongly related to the Bora wind. The outflow in Chioggia inlet occurs during Bora events but shows a slightly weaker correlation with the wind speed, while Malamocco inlet shows little or no influence of Bora winds on flow patterns. A net residual inflow through Lido and Malamocco inlet was found, while outflow prevails in Chioggia inlet. During Bora events, the average residual inflow increased three-fold in Lido inlet, whereas the outflow in Chioggia inlet doubled. The current velocities in Lido and Chioggia inlets are best described by an exponential function of wind velocity with exponents of −0.1187 and −0.0924, respectively. The response to Sirocco events was evident mainly in Chioggia inlet. Specifically, there was a slow down of the outflow in linear proportion to wind speed. In excess of 10 m/s a complete current reversal was observed. Lido and Malamocco inlets showed little or no response to Sirocco winds, except during rare cases when wind speeds exceeded 15 m/s.     

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.     

Magnetic instabilities in rapidly rotating spherical geometries. III. The effect of differential rotation

Abstract

We investigate the influence of differential rotation on magnetic instabilities for an electrically conducting fluid in the presence of a toroidal basic state of magnetic field B ₀ = B_MB₀(r, θ)1 _φ and flow U₀ = U_MU₀ (r, θ)1_φ, [(r, θ, φ) are spherical polar coordinates]. The fluid is confined in a rapidly rotating, electrically insulating, rigid spherical container. In the first instance the influence of differential rotation on established magnetic instabilities is studied. These can belong to either the ideal or the resistive class, both of which have been the subject of extensive research in parts I and II of this series. It was found there, that in the absence of differential rotation, ideal modes (driven by gradients of B ₀) become unstable for A_c ? 200 whereas resistive instabilities (generated by magnetic reconnection processes near critical levels, i.e. zeros of B₀) require A_c ? 50. Here, Λ is the Elsasser number, a measure of the magnetic field strength and Λ_c is its critical value at marginal stability. Both types of instability can be stabilised by adding differential rotation into the system. For the resistive modes the exact form of the differential rotation is not important whereas for the ideal modes only a rotation rate which increases outward from the rotation axis has a stabilising effect. We found that in all cases which we investigated Λ_c increased rapidly and the modes disappeared when R_m ≈ O(Λ_C), where the magnetic Reynolds number R_m is a measure of the strength of differential rotation. The main emphasis, however, is on instabilities which are driven by unstable gradients of the differential rotation itself, i.e. an otherwise stable fluid system is destabilised by a suitable differential rotation once the magnetic Reynolds number exceeds a certain critical value (R_m )_c. Earlier work in the cylindrical geometry has shown that the differential rotation can generate an instability if R_m ) ?O(Λ). Those results, obtained for a fixed value of Λ = 100 are extended in two ways: to a spherical geometry and to an analysis of the effect of the magnetic field strength Λ on these modes of instability. Calculations confirm that modes driven by unstable gradients of the differential rotation can exist in a sphere and they are in good agreement with the local analysis and the predictions inferred from the cylindrical geometry. For Λ = O(100), the critical value of the magnetic Reynolds number (R_m )_c Λ 100, depending on the choice of flow U₀ . Modes corresponding to azimuthal wavenumber m = 1 are the most unstable ones. Although the magnetic field B ₀ is itself a stable one, the field strength plays an important role for this instability. For all modes investigated, both for cylindrical and spherical geometries, (R_m )_c reaches a minimum value for 50 ≈ Λ ≈ 100. If Λ is increased, (R_m )_c ∝ Λ, whereas a decrease of Λ leads to a rapid increase of (R_m )_c, i.e. a stabilisation of the system. No instability was found for Λ ≈ 10 — 30. Optimum conditions for instability driven by unstable gradients of the differential rotation are therefore achieved for ≈ Λ 50 — 100, R_m ? 100. These values lead to the conclusion that the instabilities can play an important role in the dynamics of the Earth's core.     

The effect of magnetic substorms on near-ground atmospheric current

Ionosphere-magnetosphere disturbances at high latitudes, e.g. magnetic substorms, are accompanied by energetic particle precipitation and strong variations of the ionospheric electric fields and currents. These might reasonably be expected to modify the local atmospheric electric circuit. We have analysed air-earth vertical currents (AECs) measured by a long wire antenna at Esrange, northern Sweden during 35 geomagnetic substorms. Using superposed epoch analysis we compare the air-earth current variations during the 3 h before and after the time of the magnetic X-component minimum with those for corresponding local times on 35 days without substorms. After elimination of the average daily variation we can conclude that the effect of substorms on AEC is small but distinguishable. It is speculated that the AEC increases observed during about 2 h prior to the geomagnetic X-component minimum, are due to enhancement of the ionospheric electric field. During the subsequent 2 h of the substorm recovery phase, the difference between substorm and quiet atmospheric currents decreases. The amplitude of this substorm variation of AEC is estimated to be less than 50% of the amplitude of the diurnal variation in AEC during the same time interval. The statistical significance of this result was confirmed using the Van der Waerden X-test. This method was further used to show that the average air-earth current and its fluctuations increase during late expansion and early recovery phases of substorms.     

Modeling the partial ring current effect in a disturbed magnetosphere

The mathematical model of the magnetic field of the partial ring current has been proposed. This current is considered as a pair of spatial current circuits in the Northern and Southern hemispheres, either of which includes two ring zones, in the geomagnetic equator plane and ionosphere, and two zones of the field-aligned current, flowing along the geomagnetic dipole field lines and joining ring fragments of the circuit. The model parameters are: colatitude of the eastward electrojet, longitudinal shift relative to the Sun-Earth axis, circuit half-angle, and the total current flowing in the circuit. The Biot-Savart-Laplace law has been used to calculate the magnetic field of the current circuit. The magnetic field of the partial ring current has been calculated under the conditions typical of a strong magnetic storm. The technique for calculating the partial ring current intensity, using the Asym-H geomagnetic index, has been proposed.     

The North-South Asymmetry in the Green Corona

The north-south asymmetry of the Fe XIV 530.3 nm coronal emission line (the green corona) over cycle 22 was investigated. The green corona line brightness was dominant in the southern hemisphere during cycle 22 (A = –0.07), except for short periods of the ascending phase of the activity cycle. The asymmetry of the semi-annual mean over the period 1940 – 1996 was also studied. The asymmetry, during these years, reached its maximum in 1962 – 1966, and then decreased. Important periodicities of the asymmetry in cycle 22, e.g., 158 and 350 days, 2.39 years were found. Similar periodicities were also detected in the years 1940 – 1996. An FFT analysis was used to detect these periodicities.     

Modelling the effect of waves,weathering and beach development on shore platform development

A mathematical model was used to study shore platform development. Mechanical wave erosion was dependent on such variables as tidal range, wave height and period, breaker height and depth, breaker type, surf zone width and bottom roughness, submarine gradient, rock resistance and the elevational frequency of wave action within the intertidal zone. Also included were the effects of sand and pebble accumulation, cliff height and debris mobility, and downwearing associated with tidal wetting and drying. The occurrence, location and thickness of beaches often depended on initially quite minor variations in platform morphology, but owing to their abrasive or protective effect on underlying rock surfaces, they were able to produce marked differences in platform morphology. Generalizations are difficult, but the model suggests that platform gradient increases with tidal range. Platform width also increases with tidal range with slow downwearing but it decreases with fast downwearing. Platform gradient decreases and width increases with wave energy, and decreasing rock resistance and platform roughness. With low tidal range, platform gradient is generally lower and platform width greater with beaches of fine sand than with gravel, but the relationship is more variable with a high tidal range. Platform width increases and platform gradient decreases with the rate of downwearing on bare surfaces, particularly in low tidal range environments, but the pattern is less clear on beach‐covered platforms. Platforms with large amounts of beach sediment tend to be narrower and steeper than bare platform surfaces. Platform gradient increases and platform width decreases with increasing cliff height and with decreasing cliff debris mobility. Copyright © 2005 John Wiley & Sons, Ltd.     

Shear flow instabilities in the Earth’s magnetotail

Shear flow instability is studied in the Earth’s magnetotail by treating plasma as compressible. A dispersion relation is derived from the linearized MHD equations using the oscillating boundary conditions at the inner central plasma sheet/outer central plasma sheet (OCPS) interface and OCPS/plasma-sheet boundary layer (PSBL) interface, whereas the surface-mode boundary condition is used at the PSBL/lobe interface. The growth rates and the real frequencies are obtained numerically for near-Earth (\midX\mid\sim10-15 R_E) and far-Earth (\midX\mid\sim100 R_E) magnetotail parameters. The periods and wavelengths of excited modes depend sensitively on the value of plasma-sheet half thickness, L, which is taken as L=5 R_E for quiet time and L=1 R_E for disturbed time. The plasma-sheet region is found to be stable for constant plasma flows unless M_A3>1.25, where M_A3 is the Alfvén Mach number in PSBL. For near-Earth magnetotail, the excited oscillations have periods of 2–20 min (quiet time) and 0.5-4 min (disturbed time) with typical transverse wavelengths of 2–30 R_E and 0.5-6.5 R_E, respectively; whereas for distant magnetotail, the analysis predicts the oscillation periods of \sim8-80 min for quiet periods and 2–16 min for disturbed periods.