The modelled occurrence of non-thermal plasma in the ionospheric F-region and the possible consequences for ion outflows into the magnetosphere

A global, time-dependent, three-dimensional, coupled ionosphere-thermosphere model is used to predict the spatial distribution of non-thermal plasma in the F-layer. It is shown that, even for steady-state conditions with Kp as low as 3, the difference between the ion and neutral velocities often exceeds the neutral thermal speed by a factor, D', which can be as large as 4. Theoretically, highly non-Maxwellian, and probably toroidal, ion velocity distributions are expected when D' exceeds about 1.5. The lack of response of the neutral winds to sunward ion drifts in the dawn sector of the auroral oval cause this to be the region most likely to contain toroidal distributions. The maximum in D' is found in the throat region of the convection pattern, where the strong neutral winds of the afternoon sector meet the eastward ion flows of the morning sector. These predictions are of interest, not only to radar scientists searching for non-thermal ionospheric plasma, but also as one possible explanation of the initial heating and upward flows of ions in the cleft ion fountain and nightside auroral oval, both of which are a major source of plasma for the magnetosphere.

Low-energy ion outflows from the ionosphere during a major polar cap expansion — evidence for equatorward motion of inverted-V structures

Data from the Dynamics Explorer 1 satellite and the EISCAT and Sondrestrom incoherent scatter radars, have allowed a study of low-energy ion outflows from the ionosphere into the magnetosphere during a rapid expansion of the polar cap. From the combined radar data, a 200kV increase in cross-cap potential is estimated. The upflowing ions show “X” signatures in the pitch angle-time spectrograms in the expanding midnight sector of the auroral oval. These signatures reveal low-energy (below about 60eV), light-ion beams sandwiched between two regions of ion conics and are associated with inverted-V electron precipitation. The lack of mass dispersion of the poleward edge of the event, despite great differences in the times of flight, reflects the equatorward expansion of the acceleration regions at velocities similar to those of the antisunward convection. In addition, a transient burst of upflow of 0+ is observed within the cap, possibly due to enhanced Joule heating during the event.

Upwelling O+ ion source characteristics

Recent observations from the Dynamics Explorer 1 (DE-1) spacecraft have shown that the dayside auroral zone is an important source of very low-energy superthermal O^+ ions for the polar magnetosphere. When observed at 2000- to 5000-km altitude, the core of the O^+ distribution exhibits transverse heating to energies on the order of 10 eV, significant upward heat flux, and subsonic upward flow at significant flux levels exceeding 10^8 cm^{-2}s^{-1}. The term "upwelling ions" has been adopted to label these flows, which stand out in sharp contrast to the light ion polar wind flows observed in the same altitude range in the polar cap and subauroral magnetosphere. We have chosen a typical upwelling ion event for detailed study, correlating retarding ion mass spectrometer observations of the low-energy plasma with energetic ion observations and local electromagnetic field observations. The upwelling ion signature is colocated with the magnetospheric cleft as marked by precipitating energetic magnetosheath ions. The apparent ionospheric heating is clearly linked with the magnetic field signatures of strong field-aligned currents in the vicinity of the dayside polar cap boundary. Electric field and ion plasma measurements indicate that a very strong and localized convection channel or jet exists coincident with the other signatures of this event. These observations indicate that transverse ion heating to temperatures on the order of 10^5 K in the 2000- to 5000-km ionosphere is an important factor in producing heavy ion outflows into the polar magnetosphere. This result contrasts with recent suggestions that electron heating to temperatures of order 10^4 K is the most important parameter with regard to O^+ outflow.

A survey of simultaneous observations of the high-latitude ionosphere and interplanetary magnetic field with EISCAT and AMPTE-UKS

This paper surveys the results of simultaneous observations by the EISCAT incoherent scatter radar and the AMPTE-UKS satellite, made during three periods in September and October 1984, when AMPTE-UKS was in the solar wind on the dayside of the Earth and the UK-POLAR EISCAT experiment was measuring ionospheric parameters at invariant latitudes 70.8–75.0°. A total of 42 h of EISCAT convection velocity data, with 2.5 min resolution, were obtained, together with 28 h of simultaneous 5 s resolution AMPTE-UKS observations of the solar wind and interplanetary magnetic field (IMF). The general features of the AMPTE-UKS data are described in Section 2 and those of the EISCAT data are described in Sections 3 and 4. The main subjects discussed are the form of the plasma convection patterns and their dependence on all three components of the IMF (Section 5), the ionospheric response to abrupt changes in the IMF (Section 6), in particular a sharp ‘southward turning’ of the IMF on 27 October 1984, and a crossing of an IMF sector boundary. Section 7 describes ‘short lived rapid flow burst’, which are believed to be related to flux transfer events at the magnetopause.

The geomagnetic mass spectrometer— mass and energy dispersions of ionospheric ion flows into the magnetosphere

NASA's Dynamics Explorer (DE) mission was designed to study the coupling between the Earth's magnetosphere, ionosphere and neutral thermosphere1. One area of major interest is the outflow of ionospheric plasma into the magnetosphere, the scale and significance of which is only now becoming apparent with the advent of mass-resolving, low-energy ion detectors. Here we compare observations of ion flows in the polar magnetosphere, made by the retarding ion mass spectrometer (RIMS)2 on DE1, with those made simultaneously in the topside ionosphere by the ion drift meter (IDM)3 on the lower-altitude DE2 spacecraft. The results show the dayside auroral ionosphere to be a significant and highly persistent source of plasma for the magnetosphere. The upwelling ionospheric ions are spatially dispersed, according to both their energy and mass, by the combined actions of the geomagnetic field and the dawn-to-dusk convection electric field, in an effect analogous to the operation of an ion mass spectrometer.

The cleft ion fountain: a two-dimensional kinetic model

The transport of ionospheric ions from a source in the polar cleft ionosphere through the polar magnetosphere is investigated using a two-dimensional, kinetic, trajectory-based code. The transport model includes the effects of gravitation, longitudinal magnetic gradient force, convection electric fields, and parallel electric fields. Individual ion trajectories as well as distribution functions and resulting bulk parameters of density, parallel average energy, and parallel flux for a presumed cleft ionosphere source distribution are presented for various conditions to illustrate parametrically the dependences on source energies, convection electric field strengths, ion masses, and parallel electric field strengths. The essential features of the model are consistent with the concept of a cleft-based ion fountain supplying ionospheric ions to the polar magnetosphere, and the resulting plasma distributions and parameters are in general agreement with recent low-energy ion measurements from the DE 1 satellite.

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.

A diabatically-generated potential vorticity structure near the extratropical tropopause in three simulated extratropical cyclones

The structure of near-tropopause potential vorticity (PV) acts as a primary control on the evolution of extratropical cyclones. Diabatic processes such as the latent heating found in ascending moist warm conveyor belts modify PV. A dipole in diabatically-generated PV (hereafter diabatic PV) straddling the extratropical tropopause, with the positive pole above the negative pole, was diagnosed in a recently published analysis of a simulated extratropical cyclone. This PV dipole has the potential to significantly modify the propagation of Rossby waves and the growth of baroclinically-unstable waves. This previous analysis was based on a single case study simulated with 12-km horizontal grid spacing and parameterized convection. Here, the dipole is investigated in three additional cold-season extratropical cyclones simulated in both convection-parameterizing and convection-permitting model configurations. A diabatic PV dipole across the extratropical tropopause is diagnosed in all three cases. The amplitude of the dipole saturates approximately 36 hours from the time diabatic PV is accumulated. The node elevation of the dipole varies between 2-4 PVU in the three cases, and the amplitude of the system-averaged dipole varies between 0.2-0.4 PVU. The amplitude of the negative pole is similar in the convection-parameterizing and convection-permitting simulations. The positive pole, which is generated by long-wave radiative cooling, is weak in the convection-permitting simulations due to the small domain size which limits the accumulation of diabatic tendencies within the interior of the domain. The possible correspondence between the diabatic PV dipole and the extratropical tropopause inversion layer is discussed.

A process-based analysis of ocean heat uptake in an AOGCM with an eddy-permitting ocean component

About 90% of the anthropogenic increase in heat stored in the climate system is found the oceans. Therefore it is relevant to understand the details of ocean heat uptake. Here we present a detailed, process-based analysis of ocean heat uptake (OHU) processes in HiGEM1.2, an atmosphere-ocean general circulation model (AOGCM) with an eddy-permitting ocean component of 1/3 degree resolution. Similarly to various other models, HiGEM1.2 shows that the global heat budget is dominated by a downward advection of heat compensated by upward isopycnal diffusion. Only in the upper tropical ocean do we find the classical balance between downward diapycnal diffusion and upward advection of heat. The upward isopycnal diffusion of heat is located mostly in the Southern Ocean, which thus dominates the global heat budget. We compare the responses to a 4xCO2 forcing and an enhancement of the windstress forcing in the Southern Ocean. This highlights the importance of regional processes for the global ocean heat uptake. These are mainly surface fluxes and convection in the high latitudes, and advection in the Southern Ocean mid-latitudes. Changes in diffusion are less important. In line with the CMIP5 models, HiGEM1.2 shows a band of strong OHU in the mid-latitude Southern Ocean in the 4xCO2 run, which is mostly advective. By contrast, in the high-latitude Southern Ocean regions it is the suppression of convection that leads to OHU. In the enhanced windstress run, convection is strengthened at high Southern latitudes, leading to heat loss, while the magnitude of the OHU in the Southern mid-latitudes is very similar to the 4xCO2 results. Remarkably, there is only very small global OHU in the enhanced windstress run. The wind stress forcing just leads to a redistribution of heat. We relate the ocean changes at high southern latitudes to the effect of climate change on the Antarctic Circumpolar Current (ACC). It weakens in the 4xCO2 run and strengthens in the wind stress run. The weakening is due to a narrowing of the ACC, caused by an expansion of the Weddell Gyre, and a flattening of the isopycnals, which are explained by a combination of the wind stress forcing and increased precipitation.

The impact of monsoon intraseasonal variability on renewable power generation in India

India is increasingly investing in renewable technology to meet rising energy demands, with hydropower and other renewables comprising one-third of current installed capacity. Installed wind-power is projected to increase 5-fold by 2035 (to nearly 100GW) under the International Energy Agency’s New Policies scenario. However, renewable electricity generation is dependent upon the prevailing meteorology, which is strongly influenced by monsoon variability. Prosperity and widespread electrification are increasing the demand for air conditioning, especially during the warm summer. This study uses multi-decadal observations and meteorological reanalysis data to assess the impact of intraseasonal monsoon variability on the balance of electricity supply from wind-power and temperature-related demand in India. Active monsoon phases are characterised by vigorous convection and heavy rainfall over central India. This results in lower temperatures giving lower cooling energy demand, while strong westerly winds yield high wind-power output. In contrast, monsoon breaks are characterised by suppressed precipitation, with higher temperatures and hence greater demand for cooling, and lower wind-power output across much of India. The opposing relationship between wind-power supply and cooling demand during active phases (low demand, high supply) and breaks (high demand, low supply) suggests that monsoon variability will tend to exacerbate fluctuations in the so-called demand-net-wind (i.e., electrical demand that must be supplied from non-wind sources). This study may have important implications for the design of power systems and for investment decisions in conventional schedulable generation facilities (such as coal and gas) that are used to maintain the supply/demand balance. In particular, if it is assumed (as is common) that the generated wind-power operates as a price-taker (i.e., wind farm operators always wish to sell their power, irrespective of price) then investors in conventional facilities will face additional weather-volatility through the monsoonal impact on the length and frequency of production periods (i.e. their load-duration curves).

The DYMECS project: a statistical approach for the evaluation of convective storms in high-resolution NWP models

A new frontier in weather forecasting is emerging by operational forecast models now being run at convection-permitting resolutions at many national weather services. However, this is not a panacea; significant systematic errors remain in the character of convective storms and rainfall distributions. The DYMECS project (Dynamical and Microphysical Evolution of Convective Storms) is taking a fundamentally new approach to evaluate and improve such models: rather than relying on a limited number of cases, which may not be representative, we have gathered a large database of 3D storm structures on 40 convective days using the Chilbolton radar in southern England. We have related these structures to storm life-cycles derived by tracking features in the rainfall from the UK radar network, and compared them statistically to storm structures in the Met Office model, which we ran at horizontal grid length between 1.5 km and 100 m, including simulations with different subgrid mixing length. We also evaluated the scale and intensity of convective updrafts using a new radar technique. We find that the horizontal size of simulated convective storms and the updrafts within them is much too large at 1.5-km resolution, such that the convective mass flux of individual updrafts can be too large by an order of magnitude. The scale of precipitation cores and updrafts decreases steadily with decreasing grid lengths, as does the typical storm lifetime. The 200-m grid-length simulation with standard mixing length performs best over all diagnostics, although a greater mixing length improves the representation of deep convective storms.

Potential vorticity in warm conveyor belt outflow

Warm conveyor belts (WCBs) are the main ascending air masses within extratropical cyclones. They often exhibit strong condensation and precipitation, associated with ascent on large scales or embedded convection. Most of the air outflows into the upper troposphere as part of a ridge. Such ridges are an integral part of Rossby waves propagating along the tropopause and are identified with a negative potential vorticity (PV) anomaly and associated anticyclonic circulation. It has been argued that diabatic modification of PV in WCBs has an important influence on the extent of the ridge, propagation of Rossby waves and weather impacts downstream. Following the coherent ensemble of trajectories defining a WCB, PV is expected to increase with time while below the level of maximum latent heating and then decrease as trajectories ascend above the heating maximum. In models, it is found that the net change is approximately zero, so that the average PV of the WCB outflow is almost equal to the PV of its inflow. Here, the conditions necessary for this evolution are explored analytically using constraints arising from the conservation of circulation. It is argued that the net PV change is insensitive to the details of diabatic processes and the PV maximum midway along a WCB depends primarily on the net diabatic transport of mass from the inflow to the outflow layer. The main effect of diabatic processes within a WCB is to raise the isentropic level of the outflow, rather than to modify PV.

On the microphysical effects of observed cloud edge charging

Liquid layer clouds are abundant globally. Lacking strong convection, they do not become electrified by the usual thunderstorm mechanisms of collisional electrification between hydrometeors of different phases. Instead, the background global circuit current flow in fair weather is largely unaffected by the layer cloud’s presence, and, if the layer cloud is extensive horizontally, the vertical fair weather conduction current passes through the cloud. A consequence of the vertical current flow is that, at the cloud-air boundary where there is a conductivity transition and droplets form or evaporate, droplet charging occurs. Charge can affect both droplet evaporation and droplet-droplet collisions. Using new radiosonde instrumentation, the charge observed at layer cloud edges is evaluated for both these microphysical droplet processes. This shows that the charging is more likely to affect collision processes than activation, for small droplets. Enhancing the collection efficiency of small droplets modifies their evolution and propagates through the size distribution to shorten the autoconversion timescale to rain drops, and the cloud radiative properties. Because the conduction current density is influenced by both external (e.g. solar modulation of high energy particles) and internal (e.g. ENSO) factors, current flow leading to layer cloud edge charging provides a possible route for expressing solar influences on the climate system and a teleconnection mechanism for communicating internal climate variability.

Vertical structure and physical processes of the Madden-Julian oscillation: synthesis and summary

The "Vertical structure and physical processes of the Madden-Julian oscillation (MJO)" project comprises three experiments, designed to evaluate comprehensively the heating, moistening and momentum associated with tropical convection in general circulation models (GCMs). We consider here only those GCMs that performed all experiments. Some models display relatively higher or lower MJO fidelity in both initialized hindcasts and climate simulations, while others show considerable variations in fidelity between experiments. Fidelity in hindcasts and climate simulations are not meaningfully correlated. The analysis of each experiment led to the development of process-oriented diagnostics, some of which distinguished between GCMs with higher or lower fidelity in that experiment. We select the most discriminating diagnostics and apply them to data from all experiments, where possible, to determine if correlations with MJO fidelity hold across scales and GCM states. While normalized gross moist stability had a small but statistically significant correlation with MJO fidelity in climate simulations, we find no link with fidelity in medium-range hindcasts. Similarly, there is no association between timestep-to-timestep rainfall variability, identified from short hindcasts, and fidelity in medium-range hindcasts or climate simulations. Two metrics that relate precipitation to free-tropospheric moisture--the relative humidity for extreme daily precipitation, and variations in the height and amplitude of moistening with rain rate--successfully distinguish between higher- and lower-fidelity GCMs in hindcasts and climate simulations. To improve the MJO, developers should focus on relationships between convection and both total moisture and its rate of change. We conclude by offering recommendations for further experiments.

Vertical structure and physical processes of the Madden-Julian oscillation: exploring key model physics in climate simulations

Aimed at reducing deficiencies in representing the Madden-Julian oscillation (MJO) in general circulation models (GCMs), a global model evaluation project on vertical structure and physical processes of the MJO was coordinated. In this paper, results from the climate simulation component of this project are reported. It is shown that the MJO remains a great challenge in these latest generation GCMs. The systematic eastward propagation of the MJO is only well simulated in about one-fourth of the total participating models. The observed vertical westward tilt with altitude of the MJO is well simulated in good MJO models, but not in the poor ones. Damped Kelvin wave responses to the east of convection in the lower troposphere could be responsible for the missing MJO preconditioning process in these poor MJO models. Several process-oriented diagnostics were conducted to discriminate key processes for realistic MJO simulations. While large-scale rainfall partition and low-level mean zonal winds over the Indo-Pacific in a model are not found to be closely associated with its MJO skill, two metrics, including the low-level relative humidity difference between high and low rain events and seasonal mean gross moist stability, exhibit statistically significant correlations with the MJO performance. It is further indicated that increased cloud-radiative feedback tends to be associated with reduced amplitude of intraseasonal variability, which is incompatible with the radiative instability theory previously proposed for the MJO. Results in this study confirm that inclusion of air-sea interaction can lead to significant improvement in simulating the MJO.