Dynamical control of the mesosphere by orographic and non-orographic gravity wave drag during the extended northern winters of 2006 and 2009

A version of the Canadian Middle Atmosphere Model (CMAM) that is nudged toward reanalysis data up to 1 hPa is used to examine the impacts of parameterized orographic and non-orographic gravity wave drag (OGWD and NGWD) on the zonal-mean circulation of the mesosphere during the extended northern winters of 2006 and 2009 when there were two large stratospheric sudden warmings. The simulations are compared to Aura Microwave Limb Sounder (MLS) observations of mesospheric temperature, carbon monoxide (CO) and derived zonal winds. The control simulation, which uses both OGWD and NGWD, is shown to be in good agreement with MLS. The impacts of OGWD and NGWD are assessed using simulations in which those sources of wave drag are removed. In the absence of OGWD the mesospheric zonal winds in the months preceding the warmings are too strong, causing increased mesospheric NGWD, which drives excessive downwelling, resulting in overly large lower mesospheric values of CO prior to the warming. NGWD is found to be most important following the warmings when the underlying westerlies are too weak to allow much vertical propagation of the orographic gravity waves to the mesosphere. NGWD is primarily responsible for driving the circulation that results in the descent of CO from the thermosphere following the warmings. Zonal mean mesospheric winds and temperatures in all simulations are shown to be strongly constrained by (i.e. slaved to) the stratosphere. Finally, it is demonstrated that the responses to OGWD and NGWD are non-additive due to their dependence and influence on the background winds and temperatures.

Urban surface energy balance models: model characteristics and methodology for a comparison study

Many urban surface energy balance models now exist. These vary in complexity from simple schemes that represent the city as a concrete slab, to those which incorporate detailed representations of momentum and energy fluxes distributed within the atmospheric boundary layer. While many of these schemes have been evaluated against observations, with some models even compared with the same data sets, such evaluations have not been undertaken in a controlled manner to enable direct comparison. For other types of climate model, for instance the Project for Intercomparison of Land-Surface Parameterization Schemes (PILPS) experiments (Henderson-Sellers et al., 1993), such controlled comparisons have been shown to provide important insights into both the mechanics of the models and the physics of the real world. This paper describes the progress that has been made to date on a systematic and controlled comparison of urban surface schemes. The models to be considered, and their key attributes, are described, along with the methodology to be used for the evaluation.

The science case for an orbital mission to Uranus: Exploring the origins and evolution of ice giant planets

Giant planets helped to shape the conditions we see in the Solar System today and they account for more than 99% of the mass of the Sun’s planetary system. They can be subdivided into the Ice Giants (Uranus and Neptune) and the Gas Giants (Jupiter and Saturn), which differ from each other in a number of fundamental ways. Uranus, in particular is the most challenging to our understanding of planetary formation and evolution, with its large obliquity, low self-luminosity, highly asymmetrical internal field, and puzzling internal structure. Uranus also has a rich planetary system consisting of a system of inner natural satellites and complex ring system, five major natural icy satellites, a system of irregular moons with varied dynamical histories, and a highly asymmetrical magnetosphere. Voyager 2 is the only spacecraft to have explored Uranus, with a flyby in 1986, and no mission is currently planned to this enigmatic system. However, a mission to the uranian system would open a new window on the origin and evolution of the Solar System and would provide crucial information on a wide variety of physicochemical processes in our Solar System. These have clear implications for understanding exoplanetary systems. In this paper we describe the science case for an orbital mission to Uranus with an atmospheric entry probe to sample the composition and atmospheric physics in Uranus’ atmosphere. The characteristics of such an orbiter and a strawman scientific payload are described and we discuss the technical challenges for such a mission. This paper is based on a white paper submitted to the European Space Agency’s call for science themes for its large-class mission programme in 2013.

Energy and pitch-angle dispersions of LLBL/cusp ions seen at middle altitudes: predictions by the open magnetosphere model

Numerical simulations are presented of the ion distribution functions seen by middle-altitude spacecraft in the low-latitude boundary layer (LLBL) and cusp regions when reconnection is, or has recently been, taking place at the equatorial magnetopause. From the evolution of the distribution function with time elapsed since the field line was opened, both the observed energy/observation-time and pitch-angle/energy dispersions are well reproduced. Distribution functions showing a mixture of magnetosheath and magnetospheric ions, often thought to be a signature of the LLBL, are found on newly opened field lines as a natural consequence of the magnetopause effects on the ions and their flight times. In addition, it is shown that the extent of the source region of the magnetosheath ions that are detected by a satellite is a function of the sensitivity of the ion instrument . If the instrument one-count level is high (and/or solar-wind densities are low), the cusp ion precipitation detected comes from a localised region of the mid-latitude magnetopause (around the magnetic cusp), even though the reconnection takes place at the equatorial magnetopause. However, if the instrument sensitivity is high enough, then ions injected from a large segment of the dayside magnetosphere (in the relevant hemisphere) will be detected in the cusp. Ion precipitation classed as LLBL is shown to arise from the low-latitude magnetopause, irrespective of the instrument sensitivity. Adoption of threshold flux definitions has the same effect as instrument sensitivity in artificially restricting the apparent source region.

Variability of dayside high latitude convection associated with a sequence of auroral transients

10 second resolution ionospheric convection data covering the invariant latitude range from 71° to 76°, obtained by using the EISCAT UHF and VHF radars, are combined with optical data from Ny Ålesund during a sequence of auroral transients in the post-noon sector (∼ 15 MLT). Satellite observations of polar cap convection patterns suggest negative BZ and BY components of the interplanetary magnetic field. Burst-like enhancements of westward (sunward) post-noon convection were accompanied by eastward moving auroral forms at higher latitudes, above the convection reversal boundary. In this case the background convection was weak, whereas the integrated potential drop across the radar field-of-view associated with the westward flow bursts was typically ∼ 20-35 kV. The auroral phenomenon consists of a series of similar events with a mean repetition period of 8 min. A close correlation between the auroral activity and convection enhancements in the cleft ionosphere is demonstrated.

Measuring ion temperatures and studying the ion energy balance in the high-latitude ionosphere

Data are presented for a nighttime ion heating event observed by the EISCAT radar on 16 December 1988. In the experiment, the aspect angle between the radar beam and the geomagnetic field was fixed at 54.7°, which avoids any ambiguity in derived ion temperature caused by anisotropy in the ion velocity distribution function. The data were analyzed with an algorithm which takes account of the non-Maxwellian line-of-sight ion velocity distribution. During the heating event, the derived spectral distortion parameter (D∗) indicated that the distribution function was highly distorted from a Maxwellian form when the ion drift increased to 4 km s−1. The true three-dimensional ion temperature was used in the simplified ion balance equation to compute the ion mass during the heating event. The ion composition was found to change from predominantly O4 to mainly molecular ions. A theoretical analysis of the ion composition, using the MSIS86 model and published values of the chemical rate coefficients, accounts for the order-of-magnitude increase in the atomic/molecular ion ratio during the event, but does not successfully explain the very high proportion of molecular ions that was observed.

Analysis of incoherent scatter radar data from non-thermal F-region plasma

A procedure is presented for fitting incoherent scatter radar data from non-thermal F-region ionospheric plasma, using theoretical spectra previously predicted. It is found that values of the shape distortion factor D∗, associated with deviations of the ion velocity distribution from a Maxwellian distribution, and ion temperatures can be deduced (the results being independent of the path of iteration) if the angle between the line-of-sight and the geomagnetic field is larger than about 15–20°. The procedure can be used with one or both of two sets of assumptions. These concern the validity of the adopted model for the line-of-sight ion velocity distribution in the one case or for the full three-dimensional ion velocity distribution function in the other. The distribution function employed was developed to describe the line-of-sight velocity distribution for large aspect angles, but both experimental data and Monte Carlo simulations indicate that the form of the field-perpendicular distribution can also describe the distribution at more general aspect angles. The assumption of this form for the line-of-sight velocity distribution at a general aspect angle enables rigorous derivation of values of the one-dimensional, line-of-sight ion temperature. With some additional assumptions (principally that the field-parallel distribution is always Maxwellian and there is a simple relationship between the ion temperature anisotropy and the distortion of the field-perpendicular distribution from a Maxwellian), fits to data for large aspect angles enable determination of line-of-sight temperatures at all aspect angles and hence, of the average ion temperature and the ion temperature anisotropy. For small aspect angles, the analysis is restricted to the determination of the line-of-sight ion temperature because the theoretical spectrum is insensitive to non-thermal effects when the plasma is viewed along directions almost parallel to the magnetic field. This limitation is expected to apply to any realistic model of the ion velocity distribution function and its consequences are discussed. Fit strategies which allow for mixed ion composition are also considered. Examples of fits to data from various EISCAT observing programmes are presented.

Comparisons between EISCAT observations and model calculations of the high latitude ionosphere

Calculations using a numerical model of the convection dominated high latitude ionosphere are compared with observations made by EISCAT as part of the UK-POLAR Special Programme. The data used were for 24–25 October 1984, which was characterized by an unusually steady IMF, with Bz < 0 and By > 0; in the calculations it was assumed that a steady IMF implies steady convection conditions. Using the electric field models of Heppner and Maynard (1983) appropriate to By > 0 and precipitation data taken from Spiroet al. (1982), we calculated the velocities and electron densities appropriate to the EISCAT observations. Many of the general features of the velocity data were reproduced by the model. In particular, the phasing of the change from eastward to westward flow in the vicinity of the Harang discontinuity, flows near the dayside throat and a region of slow flow at higher latitudes near dusk were well reproduced. In the afternoon sector modelled velocity values were significantly less than those observed. Electron density calculations showed good agreement with EISCAT observations near the F-peak, but compared poorly with observations near 211 km. In both cases, the greatest disagreement occurred in the early part of the observations, where the convection pattern was poorly known and showed some evidence of long term temporal change. Possible causes for the disagreement between observations and calculations are discussed and shown to raise interesting and, as yet, unresolved questions concerning the interpretation of the data. For the data set used, the late afternoon dip in electron density observed near the F-peak and interpreted as the signature of the mid-latitude trough is well reproduced by the calculations. Calculations indicate that it does not arise from long residence times of plasma on the nightside, but is the signature of a gap between two major ionization sources, viz. photoionization and particle precipitation.

Flow in the high latitude ionosphere: measurements at 15s resolution made using the EISCAT ‘Polar’ experiment

We present a first overview of flows in the high latitude ionosphere observed at 15 s resolution using the U.K.-Polar EISCAT experiment. Data are described from experiments conducted on two days, 27 October 1984 and 29 August 1985, which together span the local times between about 0200 and 2130MLT and cover five different regions of ionospheric flow. With increasing local time, these are: the dawn auroral zone flow cell, the dayside region of low background flows equatorward of the flow cells, the dusk auroral zone flow cell, the boundary region between the dusk auroral zone and the polar cap, and the evening polar cap. Flows in both the equatorward and poleward portions of the auroral zone cells appear to be relatively smooth, while in the central region of high speed flow considerable variations are generally present. These have the form of irregular fluctuations on a wide range of time scales in the early morning dawn cell, and impulsive wave-like variations with periods of a few minutes in the afternoon dusk cell. In the dayside region between the flow cells, the ionosphere is often essentially stagnant for long intervals, but low amplitude ULF waves with a period of about 5 min can also occur and persist for many cycles. These conditions are punctuated at one to two hour intervals by sudden ‘flow burst’ events with impulsively generated damped wave trains. Initial burst flows are generally directed poleward and can peak at line-of-sight speeds in excess of 1 km s^{−1} after perhaps 45 s. Flows in the polar cap are reasonably smooth on time scales of a few minutes and show no evidence for the presence of ULF waves. Under most, but not all, of the above conditions, the beam-swinging algorithm used to determine background vector flows should produce meaningful results. Comparison of these flow data with simultaneous plasma and magnetic field measurements in the solar wind, made by the AMPTE IRM and UKS spacecraft, emphasizes the strong control exerted on high latitude flows by the north-south component of the IMF.

Model predictions of the occurrence of non-Maxwellian plasmas, and analysis of their effects on EISCAT data

The recent identification of non-thermal plasmas using EISCAT data has been made possible by their occurrence during large, short-lived flow bursts. For steady, yet rapid, ion convection the only available signature is the shape of the spectrum, which is unreliable because it is open to distortion by noise and sampling uncertainty and can be mimicked by other phenomena. Nevertheless, spectral shape does give an indication of the presence of non-thermal plasma, and the characteristic shape has been observed for long periods (of the order of an hour or more) in some experiments. To evaluate this type of event properly one needs to compare it to what would be expected theoretically. Predictions have been made using the coupled thermosphere-ionosphere model developed at University College London and the University of Sheffield to show where and when non-Maxwellian plasmas would be expected in the auroral zone. Geometrical and other factors then govern whether these are detectable by radar. The results are applicable to any incoherent scatter radar in this area, but the work presented here concentrates on predictions with regard to experiments on the EISCAT facility.

EISCAT observations of plasma convection and the high-latitude, winter F-region during substorm activity

A 24 h period of observations by the EISCAT radar and other ground-based instrumentation is used to study the role of plasma convection in determining the morphology of the high-latitude F-region during winter. It is suggested that, in the afternoon sector of the polar convection pattern, rapid zonal (westward) flows caused low F-region electron densities due to an extension of the mid-latitude trough far into the sunlit hemisphere. Low densities on the dawn side prior to 0600 UT may also have been due to a trough-like feature. Although the generation mechanism is unclear, the trough may be the fossil remnant of a substorm. Around midnight, high F-region densities were seen, probably due to plasma flow emerging from the cap through soft particle precipitation in the auroral oval. Two substorms occurred at times when the radar was south of the auroral oval. Both caused enhanced convection speeds, a swing to equatorward flow, enhanced E-region densities and a depleted F-region. The first was seen as a Westward Travelling Surge, and the swing to purely southward flow which followed the surge front did not return to westward flows until 80–110 min later. The Harang discontinuity was observed co-rotating eastwards between the substorms, 65 ± 30 min before the separatrix between the dawn and dusk convection cells.

Field-perpendicular and field-aligned plasma flows observed by EISCAT during a prolonged period of northward IMF

The effect of a prolonged period of strongly northward Interplanetary Magnetic Field (IMF) on the high-latitude F-region is studied using data from the EISCAT Common Programme Zero mode of operation on 11–12 August 1982. The analysis of the raw autocorrelation functions is kept to the directly derived parameters Ne, Te, Ti and velocity, and limits are defined for the errors introduced by assumptions about ion composition and by changes in the transmitted power and system constant. Simple data-cleaning criteria are employed to eliminate problems due to coherent signals and large background noise levels. The observed variations in plasma densities, temperatures and velocities are interpreted in terms of supporting data from ISEE-3 and local riometers and magnetometers. Both field-aligned and field-perpendicular plasma flows at Tromsø showed effects of the northward IMF: convection was slow and irregular and field-aligned flow profiles were characteristic of steady-state polar wind outflow with flux of order 1012 m−2 s−1. This period followed a strongly southward IMF which had triggered a substorm. The substorm gave enhanced convection, with a swing to equatorward flow and large (5 × 1012 m−2 s−1), steady-state field-aligned fluxes, leading to the possibility of O+ escape into the magnetosphere. The apparent influence of the IMF over both field-perpendicular and field-aligned flows is explained in terms of the cross-cap potential difference and the location of the auroral oval.

Discovery of 157W and 161Os

The nuclides 157W and 161Os have been discovered in reactions of 58Ni ion beams with a 106Cd target. The 161Os α -decay energy and half-life were 6890±12 keV and 640±60 μs. The daughter 157W nuclei β -decayed with a half-life of 275±40 ms, populating both low-lying α-decaying states in 157Ta, which is consistent with a 7/2− ground state in 157W. Fine structure observed in the α decay of 161Os places the lowest excited state in 157W with Iπ=9/2− at 318±30 keV. The branching ratio of View the MathML source indicates that 161Os also has a 7/2− ground state. Shell-model calculations analysing the effects of monopole shifts and a tensor force on the relative energies of 2f7/2 and 1h9/2 neutron states in N=83 isotones are presented.

Further considerations of cosmic ray modulation of infra-red radiation in the atmosphere

Understanding effects of ionisation in the lower atmosphere is a new interdisciplinary area, crossing the traditionally distinct scientific boundaries between astro-particle and atmospheric physics and also requiring understanding of both heliospheric and magnetospheric influences on cosmic rays. Following the paper of Erlykin et al. (2014) we develop further the interpretation of our observed changes in long-wave (LW) radiation, Aplin and Lockwood (2013) by taking account of both cosmic ray ionisation yields and atmospheric radiative transfer. To demonstrate this, we show that the thermal structure of the whole atmosphere needs to be considered along with the vertical profile of ionisation. Allowing for, in particular, ionisation by all components of a cosmic ray shower and not just by the muons, reveals that the effect we have detected is certainly not inconsistent with laboratory observations of the LW absorption cross section. The analysis presented here, although very different from that of Erlykin et al., does come to the same conclusion that the events detected by AL were not caused by individual cosmic ray primaries – not because it is impossible on energetic grounds, but because events of the required energy are too infrequent for the 12 h_1 rate at which they were seen by the AL experiment. The present paper numerically models the effect of three different scenario changes to the primary GCR spectrum which all reproduce the required magnitude of the effect observed by AL. However, they cannot solely explain the observed delay in the peak effect which, if confirmed, would appear to open up a whole new and interesting area in the study of water oligomers and their effects on LW radiation. We argue that a technical artefact in the AL experiment is highly unlikely and that our initial observations merit both a wide-ranging follow-up experiment and more rigorous, self-consistent, three-dimensional radiative transfer modelling