High-latitude particle precipitation and its relationship to magnetospheric source regions

Two central issues in magnetospheric research are understanding the mapping of the low-altitude ionosphere to the distant regions of the magnetsphere, and understanding the relationship between the small-scale features detected in the various regions of the ionosphere and the global properties of the magnetosphere. The high-latitude ionosphere, through its magnetic connection to the outer magnetosphere, provides an important view of magnetospheric boundaries and the physical processes occurring there. All physical manifestations of this magnetic connectivity (waves, particle precipitation, etc.), however, have non-zero propagation times during which they are convected by the large-scale magnetospheric electric field, with phenomena undergoing different convection distances depending on their propagation times. Identification of the ionospheric signatures of magnetospheric regions and phenomena, therefore, can be difficult. Considerable progress has recently been made in identifying these convection signatures in data from low- and high-altitude satellites. This work has allowed us to learn much about issues such as: the rates of magnetic reconnection, both at the dayside magnetopause and in the magnetotail; particle transport across the open magnetopause; and particle acceleration at the magnetopause and the magnetotail current sheets.

Low and middle altitude cusp particle signatures for general magnetopause reconnection rate variations: 1. Theory

We present predictions of the signatures of magnetosheath particle precipitation (in the regions classified as open low-latitude boundary layer, cusp, mantle and polar cap) for periods when the interplanetary magnetic field has a southward component. These are made using the “pulsating cusp” model of the effects of time-varying magnetic reconnection at the dayside magnetopause. Predictions are made for both low-altitude satellites in the topside ionosphere and for midaltitude spacecraft in the magnetosphere. Low-altitude cusp signatures, which show a continuous ion dispersion signature, reveal "quasi-steady reconnection" (one limit of the pulsating cusp model), which persists for a period of at least 10 min. We estimate that “quasi-steady” in this context corresponds to fluctuations in the reconnection rate of a factor of 2 or less. The other limit of the pulsating cusp model explains the instantaneous jumps in the precipitating ion spectrum that have been observed at low altitudes. Such jumps are produced by isolated pulses of reconnection: that is, they are separated by intervals when the reconnection rate is zero. These also generate convecting patches on the magnetopause in which the field lines thread the boundary via a rotational discontinuity separated by more extensive regions of tangential discontinuity. Predictions of the corresponding ion precipitation signatures seen by midaltitude spacecraft are presented. We resolve the apparent contradiction between estimates of the width of the injection region from midaltitude data and the concept of continuous entry of solar wind plasma along open field lines. In addition, we reevaluate the use of pitch angle-energy dispersion to estimate the injection distance.

Dependence of convective flows and particle precipitation in the high-latitude dayside ionosphere on theXandYcomponents of the interplanetary magnetic field

The asymmetries in the convective flows, current systems, and particle precipitation in the high-latitude dayside ionosphere which are related to the equatorial plane components of the interplanetary magnetic field (IMF) are discussed in relation to the results of several recent observational studies. It is argued that all of the effects reported to date which are ascribed to the y component of the IMF can be understood, at least qualitatively, in terms of a simple theoretical picture in which the effects result from the stresses exerted on the magnetosphere consequent on the interconnection of terrestrial and interplanetary fields. In particular, relaxation under the action of these stresses allows, in effect, a partial penetration of the IMF into the magnetospheric cavity, such that the sense of the expected asymmetry effects on closed field lines can be understood, to zeroth order, in terms of the “dipole plus uniform field” model. In particular, in response to IMF By, the dayside cusp should be displaced in longitude about noon in the same sense as By in the northern hemisphere, and in the opposite sense to By in the southern hemisphere, while simultaneously the auroral oval as a whole should be shifted in the dawn-dusk direction in the opposite sense with respect to By. These expected displacements are found to be consistent with recently published observations. Similar considerations lead to the suggestion that the auroral oval may also undergo displacements in the noon-midnight direction which are associated with the x component of the IMF. We show that a previously published study of the position of the auroral oval contains strong initial evidence for the existence of this effect. However, recent results on variations in the latitude of the cusp are more ambiguous. This topic therefore requires further study before definitive conclusions can be drawn.

The equivalent-weights particle filter in a high-dimensional system

In general, particle filters need large numbers of model runs in order to avoid filter degeneracy in high-dimensional systems. The recently proposed, fully nonlinear equivalent-weights particle filter overcomes this requirement by replacing the standard model transition density with two different proposal transition densities. The first proposal density is used to relax all particles towards the high-probability regions of state space as defined by the observations. The crucial second proposal density is then used to ensure that the majority of particles have equivalent weights at observation time. Here, the performance of the scheme in a high, 65 500 dimensional, simplified ocean model is explored. The success of the equivalent-weights particle filter in matching the true model state is shown using the mean of just 32 particles in twin experiments. It is of particular significance that this remains true even as the number and spatial variability of the observations are changed. The results from rank histograms are less easy to interpret and can be influenced considerably by the parameter values used. This article also explores the sensitivity of the performance of the scheme to the chosen parameter values and the effect of using different model error parameters in the truth compared with the ensemble model runs.

The effect of the equivalent-weights particle filter on dynamical balance in a primitive equation model

The disadvantage of the majority of data assimilation schemes is the assumption that the conditional probability density function of the state of the system given the observations [posterior probability density function (PDF)] is distributed either locally or globally as a Gaussian. The advantage, however, is that through various different mechanisms they ensure initial conditions that are predominantly in linear balance and therefore spurious gravity wave generation is suppressed. The equivalent-weights particle filter is a data assimilation scheme that allows for a representation of a potentially multimodal posterior PDF. It does this via proposal densities that lead to extra terms being added to the model equations and means the advantage of the traditional data assimilation schemes, in generating predominantly balanced initial conditions, is no longer guaranteed. This paper looks in detail at the impact the equivalent-weights particle filter has on dynamical balance and gravity wave generation in a primitive equation model. The primary conclusions are that (i) provided the model error covariance matrix imposes geostrophic balance, then each additional term required by the equivalent-weights particle filter is also geostrophically balanced; (ii) the relaxation term required to ensure the particles are in the locality of the observations has little effect on gravity waves and actually induces a reduction in gravity wave energy if sufficiently large; and (iii) the equivalent-weights term, which leads to the particles having equivalent significance in the posterior PDF, produces a change in gravity wave energy comparable to the stochastic model error. Thus, the scheme does not produce significant spurious gravity wave energy and so has potential for application in real high-dimensional geophysical applications.

Quantifying particle size and turbulent scale dependence of dust uplift in the Sahara using aircraft measurements

The first size-resolved airborne measurements of dust fluxes and the first dust flux measurements from the central Sahara are presented and compared with a parameterization by Kok (2011a). High-frequency measurements of dust size distribution were obtained from 0.16 to 300 µm diameter, and eddy covariance fluxes were derived. This is more than an order of magnitude larger size range than previous flux estimates. Links to surface emission are provided by analysis of particle drift velocities. Number flux is described by a −2 power law between 1 and 144 µm diameter, significantly larger than the 12 µm upper limit suggested by Kok (2011a). For small particles, the deviation from a power law varies with terrain type and the large size cutoff is correlated with atmospheric vertical turbulent kinetic energy, suggesting control by vertical transport rather than emission processes. The measured mass flux mode is in the range 30–100 µm. The turbulent scales important for dust flux are from 0.1 km to 1–10 km. The upper scale increases during the morning as boundary layer depth and eddy size increase. All locations where large dust fluxes were measured had large topographical variations. These features are often linked with highly erodible surface features, such as wadis or dunes. We also hypothesize that upslope flow and flow separation over such features enhance the dust flux by transporting large particles out of the saltation layer. The tendency to locate surface flux measurements in open, flat terrain means these favored dust sources have been neglected in previous studies.

Energetic particle influence on the Earth's atmosphere

This manuscript gives an up-to-date and comprehensive overview of the effects of energetic particle precipitation (EPP) onto the whole atmosphere, from the lower thermosphere/mesosphere through the stratosphere and troposphere, to the surface. The paper summarizes the different sources and energies of particles, principally galactic cosmic rays (GCRs), solar energetic particles (SEPs) and energetic electron precipitation (EEP). All the proposed mechanisms by which EPP can affect the atmosphere are discussed, including chemical changes in the upper atmosphere and lower thermosphere, chemistry-dynamics feedbacks, the global electric circuit and cloud formation. The role of energetic particles in Earth’s atmosphere is a multi-disciplinary problem that requires expertise from a range of scientific backgrounds. To assist with this synergy, summary tables are provided, which are intended to evaluate the level of current knowledge of the effects of energetic particles on processes in the entire atmosphere.

Enhancing the recovery of tiger nut (Cyperus esculentus) oil by mechanical pressing: moisture content, particle size, high pressure and enzymatic pre-treatment effects

Tiger nut (Cyperus esculentus) tuber contains oil that is high in monounsaturated fatty acids, and this oil makes up about 23% of the tuber. The study aimed at evaluating the impact of several factors and enzymatic pre-treatment on the recovery of pressed tiger nut oil. Smaller particles were more favourable for pressing. High pressure pre-treatment did not increase oil recovery but enzymatic treatment did. The highest yield obtained by enzymatic treatment prior to mechanical extraction was 33 % on a dry defatted basis, which represents a recovery of 90 % of the oil. Tiger nut oil consists mainly of oleic acid; its acid and peroxide values reflect the high stability of the oil.

Twin experiments with the equivalent weights particle filter and HadCM3

This paper investigates the use of a particle filter for data assimilation with a full scale coupled ocean–atmosphere general circulation model. Synthetic twin experiments are performed to assess the performance of the equivalent weights filter in such a high-dimensional system. Artificial 2-dimensional sea surface temperature fields are used as observational data every day. Results are presented for different values of the free parameters in the method. Measures of the performance of the filter are root mean square errors, trajectories of individual variables in the model and rank histograms. Filter degeneracy is not observed and the performance of the filter is shown to depend on the ability to keep maximum spread in the ensemble.

Aspects of particle filtering in high-dimensional spaces

Nonlinear data assimilation is high on the agenda in all fields of the geosciences as with ever increasing model resolution and inclusion of more physical (biological etc.) processes, and more complex observation operators the data-assimilation problem becomes more and more nonlinear. The suitability of particle filters to solve the nonlinear data assimilation problem in high-dimensional geophysical problems will be discussed. Several existing and new schemes will be presented and it is shown that at least one of them, the Equivalent-Weights Particle Filter, does indeed beat the curse of dimensionality and provides a way forward to solve the problem of nonlinear data assimilation in high-dimensional systems.

Implicit equal-weights particle filter

Filter degeneracy is the main obstacle for the implementation of particle filter in non-linear high-dimensional models. A new scheme, the implicit equal-weights particle filter (IEWPF), is introduced. In this scheme samples are drawn implicitly from proposal densities with a different covariance for each particle, such that all particle weights are equal by construction. We test and explore the properties of the new scheme using a 1,000-dimensional simple linear model, and the 1,000-dimensional non-linear Lorenz96 model, and compare the performance of the scheme to a Local Ensemble Kalman Filter. The experiments show that the new scheme can easily be implemented in high-dimensional systems and is never degenerate, with good convergence properties in both systems.

Single-particle tracking uncovers dynamics of glutamate-induced retrograde transport of NF-κB p65 in living neurons

Retrograde transport of NF-κB from the synapse to the nucleus in neurons is mediated by the dynein/dynactin motor complex and can be triggered by synaptic activation. The calibre of axons is highly variable ranging down to 100 nm, aggravating the investigation of transport processes in neurites of living neurons using conventional light microscopy. In this study we quantified for the first time the transport of the NF-κB subunit p65 using high-density single-particle tracking in combination with photoactivatable fluorescent proteins in living mouse hippocampal neurons. We detected an increase of the mean diffusion coefficient (Dmean) in neurites from 0.12 ± 0.05 µm2/s to 0.61 ± 0.03 µm2/s after stimulation with glutamate. We further observed that the relative amount of retrogradely transported p65 molecules is increased after stimulation. Glutamate treatment resulted in an increase of the mean retrograde velocity from 10.9 ± 1.9 to 15 ± 4.9 µm/s, whereas a velocity increase from 9 ± 1.3 to 14 ± 3 µm/s was observed for anterogradely transported p65. This study demonstrates for the first time that glutamate stimulation leads to an increased mobility of single NF-κB p65 molecules in neurites of living hippocampal neurons.

Coupling an aerosol box model with one-dimensional flow: a tool for understanding observations of new particle formation events

Field observations of new particle formation and the subsequent particle growth are typically only possible at a fixed measurement location, and hence do not follow the temporal evolution of an air parcel in a Lagrangian sense. Standard analysis for determining formation and growth rates requires that the time-dependent formation rate and growth rate of the particles are spatially invariant; air parcel advection means that the observed temporal evolution of the particle size distribution at a fixed measurement location may not represent the true evolution if there are spatial variations in the formation and growth rates. Here we present a zero-dimensional aerosol box model coupled with one-dimensional atmospheric flow to describe the impact of advection on the evolution of simulated new particle formation events. Wind speed, particle formation rates and growth rates are input parameters that can vary as a function of time and location, using wind speed to connect location to time. The output simulates measurements at a fixed location; formation and growth rates of the particle mode can then be calculated from the simulated observations at a stationary point for different scenarios and be compared with the ‘true’ input parameters. Hence, we can investigate how spatial variations in the formation and growth rates of new particles would appear in observations of particle number size distributions at a fixed measurement site. We show that the particle size distribution and growth rate at a fixed location is dependent on the formation and growth parameters upwind, even if local conditions do not vary. We also show that different input parameters used may result in very similar simulated measurements. Erroneous interpretation of observations in terms of particle formation and growth rates, and the time span and areal extent of new particle formation, is possible if the spatial effects are not accounted for.