Transpolar voltage and polar cap flux during the substorm cycle and steady convection events

Transpolar voltages observed during traversals of the polar cap by the Defense Meteorological Satellite Program (DMSP) F-13 spacecraft during 2001 are analyzed using the expanding-contracting polar cap model of ionospheric convection. Each of the 10,216 passes is classified by its substorm phase or as a steady convection event (SCE) by inspection of the AE indices. For all phases, we detect a contribution to the transpolar voltage by reconnection in both the dayside magnetopause and in the crosstail current sheet. Detection of the IMF influence is 97% certain during quiet intervals and >99% certain during substorm/SCE growth phases but falls to 75% in substorm expansion phases: It is only 27% during SCEs. Detection of the influence of the nightside voltage is only 19% certain during growth phases, rising during expansion phases to a peak of 96% in recovery phases: During SCEs, it is >99%. The voltage during SCEs is dominated by the nightside, not the dayside, reconnection. On average, substorm expansion phases halt the growth phase rise in polar cap flux rather than reversing it. The main destruction of the excess open flux takes place during the 6- to 10-hour interval after the recovery phase (as seen in AE) and at a rate which is relatively independent of polar cap flux because the NENL has by then retreated to the far tail. The best estimate of the voltage associated with viscous-like transfer of closed field lines into the tail is around 10 kV.

An overview of Earth’s global electric circuit and atmospheric conductivity

The Earth’s global atmospheric electric circuit depends on the upper and lower atmospheric boundaries formed by the ionosphere and the planetary surface. Thunderstorms and electrified rain clouds drive a DC current (∼1 kA) around the circuit, with the current carried by molecular cluster ions; lightning phenomena drive the AC global circuit. The Earth’s near-surface conductivity ranges from 10−7 S m−1 (for poorly conducting rocks) to 10−2 S m−1 (for clay or wet limestone), with a mean value of 3.2 S m−1 for the ocean. Air conductivity inside a thundercloud, and in fair weather regions, depends on location (especially geomagnetic latitude), aerosol pollution and height, and varies from ∼10−14 S m−1 just above the surface to 10−7 S m−1 in the ionosphere at ∼80 km altitude. Ionospheric conductivity is a tensor quantity due to the geomagnetic field, and is determined by parameters such as electron density and electron–neutral particle collision frequency. In the current source regions, point discharge (coronal) currents play an important role below electrified clouds; the solar wind-magnetosphere dynamo and the unipolar dynamo due to the terrestrial rotating dipole moment also apply atmospheric potential differences. Detailed measurements made near the Earth’s surface show that Ohm’s law relates the vertical electric field and current density to air conductivity. Stratospheric balloon measurements launched from Antarctica confirm that the downward current density is ∼1 pA m−2 under fair weather conditions. Fortuitously, a Solar Energetic Particle (SEP) event arrived at Earth during one such balloon flight, changing the observed atmospheric conductivity and electric fields markedly. Recent modelling considers lightning discharge effects on the ionosphere’s electric potential (∼+250 kV with respect to the Earth’s surface) and hence on the fair weather potential gradient (typically ∼130 V m−1 close to the Earth’s surface. We conclude that cloud-to-ground (CG) lightning discharges make only a small contribution to the ionospheric potential, and that sprites (namely, upward lightning above energetic thunderstorms) only affect the global circuit in a miniscule way. We also investigate the effects of mesoscale convective systems on the global circuit.

An optical study of multiple NEIAL events driven by low energy electron precipitation

Optical data are compared with EISCAT radar observations of multiple Naturally Enhanced Ion-Acoustic Line (NEIAL) events in the dayside cusp. This study uses narrow field of view cameras to observe small-scale, short-lived auroral features. Using multiple-wavelength optical observations, a direct link between NEIAL occurrences and low energy (about 100 eV) optical emissions is shown. This is consistent with the Langmuir wave decay interpretation of NEIALs being driven by streams of low-energy electrons. Modelling work connected with this study shows that, for the measured ionospheric conditions and precipitation characteristics, growth of unstable Langmuir (electron plasma) waves can occur, which decay into ion-acoustic wave modes. The link with low energy optical emissions shown here, will enable future studies of the shape, extent, lifetime, grouping and motions of NEIALs.

Dicopper(II) complexes of a new di-para-xylyldioxatetraazamacrocycle and cascade species with dicarboxylate anions: thermodynamics and structural properties

The new dioxatetraazamacrocycle (L-1) was synthesized by a 2 + 2 condensation and characterized. Stability constants of its copper(II) complexes were determined by spectrophotometry in DMSO at 298.2 K in 0. 10 mol dm(-3) KClO4. Mainly dinuclear complexes are formed and the presence of mononuclear species is dependent on the counterion (Cl- or ClO4-). The association constants of the dinuclear copper(II) complexes with dicarboxylate anions [oxalate (oxa(2-)), malonate (mal(2-)), succinate (suc(2-)), and glutarate (glu(2-))] were also determined by spectrophotometry at 298.2 K in DMSO, and it was found that values decrease with an increase of the alkyl chain between the carboxylate groups. X-Band EPR spectra of the dicopper(II) complexes and of their cascade species in frozen DMSO exhibit dipole-dipole coupling, and their simulation, together with their UV-vis spectra, showed that the copper centres of the complexes in solution had square pyramidal geometries though with different distortions. From the experimental data, it was also possible to predict the Cu...Cu distances, the minimum being found at 6.4 angstrom for the (Cu2LCl4)-Cl-1 complex and then successively this distance slightly increases when the chloride anions are replaced by dicarboxylate anions, from 6.6 angstrom for oxa(2-) to 7.8 for glu(2-). The crystal structures of the dinuclear copper cascade species with oxa(2-) and suc(2-) were determined and showed one anion bridging both copper centres and Cu...Cu distances of 5.485(7) angstrom and 6.442(8) angstrom, respectively.

Potassium selective calix[4]semitubes

A new class of ionophore consisting of two calix[4]arene units linked through the lower rim by two ethylene chains, in combination with propyl ether and phenolic functional groups, has been developed. These calix[4]semitube molecules exhibit remarkable selectivity and fast complexation kinetics for potassium over all Group 1 metal cations. Molecular modelling studies, using structural models derived from crystallographic data, suggest the potassium cation is complexed by a horizontal, side-on route and not through the calix[4]arene annulus. The length of the bridging alkylene chain between the respective calix[4]arenes of the semitube structure dictates the strength and selectivity of alkali metal cation binding.

Molecular characterisation of the gut microflora of healthy and inflammatory bowel disease cats using fluorescence in situ hybridisation with special reference to Desulfovibrio spp

Inflammatory bowel disease (IBD) is a common cause of chronic large bowel diarrhoea in cats. Although the aetiology of IBD is unknown, an immune-mediated response to a luminal antigen is thought to be involved. As knowledge concerning the colonic microflora of cats is limited and requires further investigation, the purpose of this study was to determine the presence of specific bacterial groups in normal and IBD cats, and the potential role they play in the health of the host. Total bacterial populations, Bacteroides spp., Bifidobacterium spp., Clostridium histolyticum subgp., Lactobacillus-Enterococcus subgp. and Desulfovibrio spp. were enumerated in 34 healthy cats and 11 IBD cats using fluorescence in situ hybridisation. The study is one of the first to show the presence of Desulfovibrio in cats. Total bacteria, Bifidobacterium spp. and Bacteroides spp. counts were all significantly higher in healthy cats when compared with IBD cats, whereas Desulfovibrio spp. (producers of toxic sulphides) numbers were found to be significantly higher in colitic cats. The information obtained from this study suggests that modulation of bacterial flora by increasing bifidobacteria and decreasing Desulfovibrio spp. may be beneficial to cats with IBD. Dietary intervention may be an important aspect of their treatment.

Twentieth century secular decrease in the atmospheric potential gradient

Current flowing in the global atmospheric electrical circuit (AEC) substantially decreased during the twentieth century. Fair-weather potential gradient (PG) observations in Scotland and Shetland show a previously unreported annual decline from 1920 to 1980, when the measurements ceased. A 25% reduction in PG occurred in Scotland 1920–50, with the maximum decline during the winter months. This is quantitatively explained by a decrease in cosmic rays (CR) increasing the thunderstorm-electrosphere coupling resistance, reducing the ionospheric potential VI. Independent measurements of VI also suggest a reduction of 27% from 1920–50. The secular decrease will influence fair weather atmospheric electrical parameters, including ion concentrations and aerosol electrification. Between 1920–50, the PG showed a negative correlation with global temperature, despite the positive correlation found recently between surface temperature and VI. The 1980s stabilisation in VI may arise from compensation of the continuing CR-induced decline by increases in global temperature and convective electrification.

A potential lag between the open solar magnetic source flux and solar EUV and X-ray emissions as measured by the Earth's ionosphere during total solar eclipses

Measurements of the ionospheric E-region during total solar eclipses have been used to provide information about the evolution of the solar magnetic field and EUV and X-ray emissions from the solar corona and chromosphere. By measuring levels of ionisation during an eclipse and comparing these measurements with an estimate of the unperturbed ionisation levels (such as those made during a control day, where available) it is possible to estimate the percentage of ionising radiation being emitted by the solar corona and chromosphere. Previously unpublished data from the two eclipses presented here are particularly valuable as they provide information that supplements the data published to date. The eclipse of 23 October 1976 over Australia provides information in a data gap that would otherwise have spanned the years 1966 to 1991. The eclipse of 4 December 2002 over Southern Africa is important as it extends the published sequence of measurements. Comparing measurements from eclipses between 1932 and 2002 with the solar magnetic source flux reveals that changes in the solar EUV and X-ray flux lag the open source flux measurements by approximately 1.5 years. We suggest that this unexpected result comes about from changes to the relative size of the limb corona between eclipses, with the lag representing the time taken to populate the coronal field with plasma hot enough to emit the EUV and X-rays ionising our atmosphere.

Electromagnetic atmosphere-plasma coupling: the global atmospheric electric circuit

A description is given of the global atmospheric electric circuit operating between the Earth’s surface and the ionosphere. Attention is drawn to the huge range of horizontal and vertical spatial scales, ranging from 10−9 m to 1012 m, concerned with the many important processes at work. A similarly enormous range of time scales is involved from 10−6 s to 109 s, in the physical effects and different phenomena that need to be considered. The current flowing in the global circuit is generated by disturbed weather such as thunderstorms and electrified rain/shower clouds, mostly occurring over the Earth’s land surface. The profile of electrical conductivity up through the atmosphere, determined mainly by galactic cosmic ray ionization, is a crucial parameter of the circuit. Model simulation results on the variation of the ionospheric potential, ∼250 kV positive with respect to the Earth’s potential, following lightning discharges and sprites are summarized. Experimental results comparing global circuit variations with the neutron rate recorded at Climax, Colorado, are then discussed. Within the return (load) part of the circuit in the fair weather regions remote from the generators, charge layers exist on the upper and lower edges of extensive layer clouds; new experimental evidence for these charge layers is also reviewed. Finally, some directions for future research in the subject are suggested.

Retrieval of mesospheric electron densities using an optimal estimation inverse method

We present a new method to determine mesospheric electron densities from partially reflected medium frequency radar pulses. The technique uses an optimal estimation inverse method and retrieves both an electron density profile and a gradient electron density profile. As well as accounting for the absorption of the two magnetoionic modes formed by ionospheric birefringence of each radar pulse, the forward model of the retrieval parameterises possible Fresnel scatter of each mode by fine electronic structure, phase changes of each mode due to Faraday rotation and the dependence of the amplitudes of the backscattered modes upon pulse width. Validation results indicate that known profiles can be retrieved and that χ2 tests upon retrieval parameters satisfy validity criteria. Application to measurements shows that retrieved electron density profiles are consistent with accepted ideas about seasonal variability of electron densities and their dependence upon nitric oxide production and transport.

A linear model of gravity wave drag for hydrostatic sheared flow over elliptical mountains

An analytical model of orographic gravity wave drag due to sheared flow past elliptical mountains is developed. The model extends the domain of applicability of the well-known Phillips model to wind profiles that vary relatively slowly in the vertical, so that they may be treated using a WKB approximation. The model illustrates how linear processes associated with wind profile shear and curvature affect the drag force exerted by the airflow on mountains, and how it is crucial to extend the WKB approximation to second order in the small perturbation parameter for these effects to be taken into account. For the simplest wind profiles, the normalized drag depends only on the Richardson number, Ri, of the flow at the surface and on the aspect ratio, γ, of the mountain. For a linear wind profile, the drag decreases as Ri decreases, and this variation is faster when the wind is across the mountain than when it is along the mountain. For a wind that rotates with height maintaining its magnitude, the drag generally increases as Ri decreases, by an amount depending on γ and on the incidence angle. The results from WKB theory are compared with exact linear results and also with results from a non-hydrostatic nonlinear numerical model, showing in general encouraging agreement, down to values of Ri of order one.

Effect of shear rupture on aggregate scale formation in sea ice

A discrete element model is used to study shear rupture of sea ice under convergent wind stresses. The model includes compressive, tensile, and shear rupture of viscous elastic joints connecting floes that move under the action of the wind stresses. The adopted shear rupture is governed by Coulomb’s criterion. The ice pack is a 400 km long square domain consisting of 4 km size floes. In the standard case with tensile strength 10 times smaller than the compressive strength, under uniaxial compression the failure regime is mainly shear rupture with the most probable scenario corresponding to that with the minimum failure work. The orientation of cracks delineating formed aggregates is bimodal with the peaks around the angles given by the wing crack theory determining diamond-shaped blocks. The ice block (floe aggregate) size decreases as the wind stress gradient increases since the elastic strain energy grows faster leading to a higher speed of crack propagation. As the tensile strength grows, shear rupture becomes harder to attain and compressive failure becomes equally important leading to elongation of blocks perpendicular to the compression direction and the blocks grow larger. In the standard case, as the wind stress confinement ratio increases the failure mode changes at a confinement ratio within 0.2–0.4, which corresponds to the analytical critical confinement ratio of 0.32. Below this value, the cracks are bimodal delineating diamond shape aggregates, while above this value failure becomes isotropic and is determined by small-scale stress anomalies due to irregularities in floe shape.

The correlation of ULF waves and auroral intensity before, during and after substorm expansion phase onset

We present case studies of the evolution of magnetic wave amplitudes and auroral intensity through the late growth phase and the expansion phase of the substorm cycle. We present strong evidence that substorm-related auroral enhancements are clearly and demonstrably linked to ULF wave amplitudes observed at the same location. In most cases, we find that the highest correlations are observed when the magnetometer time series is advanced in time, indicating that the ULF wave amplitudes start to grow before measured auroral intensities, though interestingly this is not always the case. Further we discuss the four possible reasons that may be able to explain both the timing and the high correlations between these two phenomena, including: a simple coincidence, an artifact of instrumental effects, the response of the ionosphere to magnetic waves and auroral particle precipitation, and finally that ULF waves and auroral particle precipitation are physically linked. We discount coincidence and instrumental effects since in the studies presented here they are unlikely or in general will contribute negligible effects, and we find that the ionospheric response to waves and precipitation can explain some, but not all of the results contained within this paper. Specifically, ionospheric response to substorm waves and auroral precipitation cannot explain that the result that previous studies have shown, that onset of ULF wave activity and the onset of auroral particle precipitation occur at the same time and in the same location. This leaves the possibility that ULF waves and auroral particles are physically linked.

Ultralow-frequency modulation of whistler-mode wave growth

Measurements from ground-based magnetometers and riometers at auroral latitudes have demonstrated that energetic (~30-300keV) electron precipitation can be modulated in the presence of magnetic field oscillations at ultra-low frequencies. It has previously been proposed that an ultra-low frequency (ULF) wave would modulate field and plasma properties near the equatorial plane, thus modifying the growth rates of whistler-mode waves. In turn, the resulting whistler-mode waves would mediate the pitch-angle scattering of electrons resulting in ionospheric precipitation. In this paper, we investigate this hypothesis by quantifying the changes to the linear growth rate expected due to a slow change in the local magnetic field strength for parameters typical of the equatorial region around 6.6RE radial distance. To constrain our study, we determine the largest possible ULF wave amplitudes from measurements of the magnetic field at geosynchronous orbit. Using nearly ten years of observations from two satellites, we demonstrate that the variation in magnetic field strength due to oscillations at 2mHz does not exceed ±10% of the background field. Modifications to the plasma density and temperature anisotropy are estimated using idealised models. For low temperature anisotropy, there is little change in the whistler-mode growth rates even for the largest ULF wave amplitude. Only for large temperature anisotropies can whistler-mode growth rates be modulated sufficiently to account for the changes in electron precipitation measured by riometers at auroral latitudes.

Optical characterization of the growth and spatial structure of a substorm onset arc

We present a detailed case study of the characteristics of auroral forms that constitute the first ionospheric signatures of substorm expansion phase onset. Analysis of the optical frequency and along-arc (azimuthal) wave number spectra provides the strongest constraint to date on the potential mechanisms and instabilities in the near-Earth magnetosphere that accompany auroral onset and which precede poleward arc expansion and auroral breakup. We evaluate the frequency and growth rates of the auroral forms as a function of azimuthal wave number to determine whether these wave characteristics are consistent with current models of the substorm onset mechanism. We find that the frequency, spatial scales, and growth rates of the auroral forms are most consistent with the cross-field current instability or a ballooning instability, most likely triggered close to the inner edge of the ion plasma sheet. This result is supportive of a near-Earth plasma sheet initiation of the substorm expansion phase. We also present evidence that the frequency and phase characteristics of the auroral undulations may be generated via resonant processes operating along the geomagnetic field. Our observations provide the most powerful constraint to date on the ionospheric manifestation of the physical processes operating during the first few minutes around auroral substorm onset.