Modeling electrocortical activity through improved local approximations of integral neural field equations

Neural field models of firing rate activity typically take the form of integral equations with space-dependent axonal delays. Under natural assumptions on the synaptic connectivity we show how one can derive an equivalent partial differential equation (PDE) model that properly treats the axonal delay terms of the integral formulation. Our analysis avoids the so-called long-wavelength approximation that has previously been used to formulate PDE models for neural activity in two spatial dimensions. Direct numerical simulations of this PDE model show instabilities of the homogeneous steady state that are in full agreement with a Turing instability analysis of the original integral model. We discuss the benefits of such a local model and its usefulness in modeling electrocortical activity. In particular, we are able to treat “patchy” connections, whereby a homogeneous and isotropic system is modulated in a spatially periodic fashion. In this case the emergence of a “lattice-directed” traveling wave predicted by a linear instability analysis is confirmed by the numerical simulation of an appropriate set of coupled PDEs.

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.

On nonlinear symmetric stability and the nonlinear saturation of symmetric instability

A nonlinear symmetric stability theorem is derived in the context of the f-plane Boussinesq equations, recovering an earlier result of Xu within a more general framework. The theorem applies to symmetric disturbances to a baroclinic basic flow, the disturbances having arbitrary structure and magnitude. The criteria for nonlinear stability are virtually identical to those for linear stability. As in Xu, the nonlinear stability theorem can be used to obtain rigorous upper bounds on the saturation amplitude of symmetric instabilities. In a simple example, the bounds are found to compare favorably with heuristic parcel-based estimates in both the hydrostatic and non-hydrostatic limits.

Rossby number expansions, slaving principles, and balance dynamics

We consider the problem of constructing balance dynamics for rapidly rotating fluid systems. It is argued that the conventional Rossby number expansion—namely expanding all variables in a series in Rossby number—is secular for all but the simplest flows. In particular, the higher-order terms in the expansion grow exponentially on average, and for moderate values of the Rossby number the expansion is, at best, useful only for times of the order of the doubling times of the instabilities of the underlying quasi-geostrophic dynamics. Similar arguments apply in a wide class of problems involving a small parameter and sufficiently complex zeroth-order dynamics. A modified procedure is proposed which involves expanding only the fast modes of the system; this is equivalent to an asymptotic approximation of the slaving relation that relates the fast modes to the slow modes. The procedure is systematic and thus capable, at least in principle, of being carried to any order—unlike procedures based on truncations. We apply the procedure to construct higher-order balance approximations of the shallow-water equations. At the lowest order quasi-geostrophy emerges. At the next order the system incorporates gradient-wind balance, although the balance relations themselves involve only linear inversions and hence are easily applied. There is a large class of reduced systems associated with various choices for the slow variables, but the simplest ones appear to be those based on potential vorticity.

Application of the direct Liapunov method to the problem of symmetric stability in the atmosphere

The problem of symmetric stability is examined within the context of the direct Liapunov method. The sufficient conditions for stability derived by Fjørtoft are shown to imply finite-amplitude, normed stability. This finite-amplitude stability theorem is then used to obtain rigorous upper bounds on the saturation amplitude of disturbances to symmetrically unstable flows.By employing a virial functional, the necessary conditions for instability implied by the stability theorem are shown to be in fact sufficient for instability. The results of Ooyama are improved upon insofar as a tight two-sided (upper and lower) estimate is obtained of the growth rate of (modal or nonmodal) symmetric instabilities.The case of moist adiabatic systems is also considered.

An exact local conservation theorem for finite-amplitude disturbances to non-parallel shear flows, with remarks on Hamiltonian structure and on Arnol'd's stability theorems

Disturbances of arbitrary amplitude are superposed on a basic flow which is assumed to be steady and either (a) two-dimensional, homogeneous, and incompressible (rotating or non-rotating) or (b) stably stratified and quasi-geostrophic. Flow over shallow topography is allowed in either case. The basic flow, as well as the disturbance, is assumed to be subject neither to external forcing nor to dissipative processes like viscosity. An exact, local ‘wave-activity conservation theorem’ is derived in which the density A and flux F are second-order ‘wave properties’ or ‘disturbance properties’, meaning that they are O(a2) in magnitude as disturbance amplitude a [rightward arrow] 0, and that they are evaluable correct to O(a2) from linear theory, to O(a3) from second-order theory, and so on to higher orders in a. For a disturbance in the form of a single, slowly varying, non-stationary Rossby wavetrain, $\overline{F}/\overline{A}$ reduces approximately to the Rossby-wave group velocity, where (${}^{-}$) is an appropriate averaging operator. F and A have the formal appearance of Eulerian quantities, but generally involve a multivalued function the correct branch of which requires a certain amount of Lagrangian information for its determination. It is shown that, in a certain sense, the construction of conservable, quasi-Eulerian wave properties like A is unique and that the multivaluedness is inescapable in general. The connection with the concepts of pseudoenergy (quasi-energy), pseudomomentum (quasi-momentum), and ‘Eliassen-Palm wave activity’ is noted. The relationship of this and similar conservation theorems to dynamical fundamentals and to Arnol'd's nonlinear stability theorems is discussed in the light of recent advances in Hamiltonian dynamics. These show where such conservation theorems come from and how to construct them in other cases. An elementary proof of the Hamiltonian structure of two-dimensional Eulerian vortex dynamics is put on record, with explicit attention to the boundary conditions. The connection between Arnol'd's second stability theorem and the suppression of shear and self-tuning resonant instabilities by boundary constraints is discussed, and a finite-amplitude counterpart to Rayleigh's inflection-point theorem noted

Idealised simulations of sting-jet cyclones

An idealised modelling study of sting-jet cyclones is presented. Sting jets are descending mesoscale jets that occur in some extratropical cyclones and produce localised regions of strong low-level winds in the frontal fracture region. Moist baroclinic lifecycle (LC1) simulations are performed with modifications to produce cyclones resembling observed sting-jet cyclones. A sting jet exists in the idealised control cyclone with similar characteristics to the sting jet in a simulation of windstorm Gudrun (a confirmed sting-jet case). Unlike in windstorm Gudrun, a low-level layer of strong moist static stability prohibits the descent of the strong winds from above the boundary layer to the surface in the idealised case. Conditional symmetric instability (CSI) exists in the cloud head and dissipates as the sting jet leaves the cloud head and descends. The descending, initially moist, sting-jet trajectories consistently have negative or near-zero saturated moist potential vorticity but moist static stability and inertial stability, consistent with CSI release; the moist static stability becomes negative during the period of most rapid descent, by which time the air is relatively dry implying conditional instability release is unlikely. Sensitivity experiments show that the existence of the sting jet is robust to changes in the initial state, and that the initial tropospheric static stability significantly impacts the descent rate of the sting jet. Inertial and conditional instability are probably being released in the experiment with the weakest initial static stability. This suggests that sting jets can arise through the release of all three instabilities associated with negative saturated moist potential vorticity. While evaporative cooling occurs along the sting-jet trajectories, a sensitivity experiment with evaporation effects turned off shows no significant change to the wind strength or descent rate of the sting jet implying that instability release is the dominant sting-jet driving mechanism.

Distinguishing the cold conveyor belt and sting jet air streams in an intense extratropical cyclone

Strong winds equatorwards and rearwards of a cyclone core have often been associated with two phenomena, the cold conveyor belt (CCB) jet and sting jets. Here, detailed observations of the mesoscale structure in this region of an intense cyclone are analysed. The {\it in-situ} and dropsonde observations were obtained during two research flights through the cyclone during the DIAMET (DIAbatic influences on Mesoscale structures in ExTratropical storms) field campaign. A numerical weather prediction model is used to link the strong wind regions with three types of ``air streams'', or coherent ensembles of trajectories: two types are identified with the CCB, hooking around the cyclone center, while the third is identified with a sting jet, descending from the cloud head to the west of the cyclone. Chemical tracer observations show for the first time that the CCB and sting jet air streams are distinct air masses even when the associated low-level wind maxima are not spatially distinct. In the model, the CCB experiences slow latent heating through weak resolved ascent and convection, while the sting jet experiences weak cooling associated with microphysics during its subsaturated descent. Diagnosis of mesoscale instabilities in the model shows that the CCB passes through largely stable regions, while the sting jet spends relatively long periods in locations characterized by conditional symmetric instability (CSI). The relation of CSI to the observed mesoscale structure of the bent-back front and its possible role in the cloud banding is discussed.

Mean climate and variability of the atmosphere and ocean on an aquaplanet

Numerical experiments are described that pertain to the climate of a coupled atmosphere–ocean–ice system in the absence of land, driven by modern-day orbital and CO2 forcing. Millennial time-scale simulations yield a mean state in which ice caps reach down to 55° of latitude and both the atmosphere and ocean comprise eastward- and westward-flowing zonal jets, whose structure is set by their respective baroclinic instabilities. Despite the zonality of the ocean, it is remarkably efficient at transporting heat meridionally through the agency of Ekman transport and eddy-driven subduction. Indeed the partition of heat transport between the atmosphere and ocean is much the same as the present climate, with the ocean dominating in the Tropics and the atmosphere in the mid–high latitudes. Variability of the system is dominated by the coupling of annular modes in the atmosphere and ocean. Stochastic variability inherent to the atmospheric jets drives variability in the ocean. Zonal flows in the ocean exhibit decadal variability, which, remarkably, feeds back to the atmosphere, coloring the spectrum of annular variability. A simple stochastic model can capture the essence of the process. Finally, it is briefly reviewed how the aquaplanet can provide information about the processes that set the partition of heat transport and the climate of Earth.

Computer simulation study of equilibrium properties of magnetic dipolar fluids

Langevin dynamics simulations are used to investigate the equilibrium magnetization properties and structure of magnetic dipolar fluids. The influence of using different boundary conditions are systematically studied. Simulation results on the initial susceptibility and magnetization curves are compared with theoretical predictions. The effect of particle aggregation is discussed in detail by performing a cluster analysis of the microstructure.

Increases in plasma sheet temperature with solar wind driving during substorm growth phases

During the substorm growth phase, magnetic reconnection extracts ~10^15 J from the solar wind through magnetic reconnection at the magnetopause, which is then stored in the magnetotail lobes. Plasma sheet pressure then increases to balance magnetic flux density increases in the lobes. We examine plasma sheet pressure, density and temperature during substorm growth phases using nine years of Cluster data (>316,000 data points). We show that plasma sheet pressure and temperature are higher during growth phases with higher solar wind driving whereas the density is approximately constant. We also show a weak correlation between plasma sheet temperature before onset and the minimum SuperMAG SML auroral index in the subsequent substorm. We discuss how energization of the plasma sheet before onset may result from thermodynamically adiabatic processes; how hotter plasma sheets may result in magnetotail instabilities and how this relates to the onset and size of the subsequent substorm expansion phase.

Ground-based and satellite observations of high-latitude auroral activity in the dusk sector of the auroral oval

On 7 December 2000, during 13:30-15:30 UT the MIRACLE all-sky camera at Ny Alesund observed auroras at high-latitudes (MLAT similar to 76) simultaneously when the Cluster spacecraft were skimming the magnetopause in the same MLT sector (at similar to 16:00-18:00 MLT). The location of the auroras (near the ionospheric convection reversal boundary) and the clear correlation between their dynamics and IMF variations suggests their close relationship with R1 currents. Consequently, we can assume that the Cluster spacecraft were making observations in the magnetospheric region associated with the auroras, although exact magnetic conjugacy between the ground-based and satellite observations did not exist. The solar wind variations appeared to control both the behaviour of the auroras and the magnetopause dynamics. Auroral structures were observed at Ny Alesund especially during periods of negative IMF B-Z. In addition, the Cluster spacecraft experienced periodic (T similar to 4 - 6 min) encounters between magnetospheric and magnetosheath plasmas. These undulations of the boundary can be interpreted as a consequence of tailward propagating magnetopause surface waves. Simultaneous dusk sector ground-based observations show weak, but discernible magnetic pulsations (Pc 5) and occasionally periodic variations (T - 2 - 3 min) in the high-latitude auroras. In the dusk sector, Pc 5 activity was stronger and had characteristics that were consistent with a field line resonance type of activity. When IMF BZ stayed positive for a longer period, the auroras were dimmer and the spacecraft stayed at the outer edge of the magnetopause where they observed electromagnetic pulsations with T similar to 1 min. We find these observations interesting especially from the viewpoint of previously presented studies relating poleward-moving high-latitude auroras with pulsation activity and MHD waves propagating at the magnetospheric boundary layers.

On the cause of a magnetospheric flux transfer event

We present a detailed investigation of a magnetospheric flux transfer event (FTE) seen by the Active Magnetospheric Tracer Explorer (AMPTE) UKS and IRM satellites around 1046 UT on October 28, 1984. This event has been discussed many times previously in the literature and has been cited as support for a variety of theories of FTE formation. We make use of a model developed to reproduce ion precipitations seen in the cusp ionosphere. The analysis confirms that the FTE is well explained as a brief excursion into an open low-latitude boundary layer (LLBL), as predicted by two theories of magnetospheric FTEs: namely, that they are bulges in the open LLBL due to reconnection rate enhancements or that they are indentations of the magnetopause by magnetosheath pressure increases (but in the presence of ongoing steady reconnection). The indentation of the inner edge of the open LLBL that these two models seek to explain is found to be shallow for this event. The ion model reproduces the continuous evolution of the ion distribution function between the sheath-like population at the event center and the surrounding magnetospheric populations; it also provides an explanation of the high-pressure core of the event as comprising field lines that were reconnected considerably earlier than those that are draped over it to give the event boundary layer. The magnetopause transition parameter is used to isolate a field rotation on the boundaries of the core, which is subjected to the tangential stress balance test. The test identifies this to be a convecting structure, which is neither a rotational discontinuity (RD) nor a contact discontinuity, but could possibly be a slow shock. In addition, evidence for ion reflection off a weak RD on the magnetospheric side of this structure is found. The event structure is consistent in many ways with features predicted for the open LLBL by analytic MHD theories and by MHD and hybrid simulations. The de Hoffman-Teller velocity of the structure is significantly different from that of the magnetosheath flow, indicating that it is not an indentation caused by a high-pressure pulse in the sheath but is consistent with the motion of newly opened field lines (different from the sheath flow because of the magnetic tension force) deduced from the best fit to the ion data. However, we cannot here rule out the possibility that the sheath flow pattern has changed in the long interval between the two satellites observing the FTE and subsequently emerging into the magnetosheath; thus this test is not conclusive in this particular case. Analysis of the fitted elapsed time since reconnection shows that the core of the event was reconnected in one pulse and the event boundary layer was reconnected in a subsequent pulse. Between these two pulses is a period of very low (but nonzero) reconnection rate, which lasts about 14 mins. Thus the analysis supports, but does not definitively verify, the concept that the FTE is a partial passage into an open LLBL caused by a traveling bulge in that layer produced by a pulse in reconnection rate.

Inferring the structure of the solar corona and inner heliosphere during the Maunder minimum using global thermodynamic magnetohydrodynamic simulations

Observations of the Sun’s corona during the space era have led to a picture of relatively constant, but cyclically varying solar output and structure. Longer-term, more indirect measurements, such as from 10Be, coupled by other albeit less reliable contemporaneous reports, however, suggest periods of significant departure from this standard. The Maunder Minimum was one such epoch where: (1) sunspots effectively disappeared for long intervals during a 70 yr period; (2) eclipse observations suggested the distinct lack of a visible K-corona but possible appearance of the F-corona; (3) reports of aurora were notably reduced; and (4) cosmic ray intensities at Earth were inferred to be substantially higher. Using a global thermodynamic MHD model, we have constructed a range of possible coronal configurations for the Maunder Minimum period and compared their predictions with these limited observational constraints. We conclude that the most likely state of the corona during—at least—the later portion of the Maunder Minimum was not merely that of the 2008/2009 solar minimum, as has been suggested recently, but rather a state devoid of any large-scale structure, driven by a photospheric field composed of only ephemeral regions, and likely substantially reduced in strength. Moreover, we suggest that the Sun evolved from a 2008/2009-like configuration at the start of the Maunder Minimum toward an ephemeral-only configuration by the end of it, supporting a prediction that we may be on the cusp of a new grand solar minimum.

Statistical characterisation of the growth and spatial scales of the substorm onset arc

We present the first multi-event study of the spatial and temporal structuring of the aurora to provide statistical evidence of the near-Earth plasma instability which causes the substorm onset arc. Using data from ground-based auroral imagers, we study repeatable signatures of along-arc auroral beads, which are thought to represent the ionospheric projection of magnetospheric instability in the near-Earth plasma sheet. We show that the growth and spatial scales of these wave-like fluctuations are similar across multiple events, indicating that each sudden auroral brightening has a common explanation. We find statistically that growth rates for auroral beads peak at low wavenumber with the most unstable spatial scales mapping to an azimuthal wavelength λ≈1700 − 2500 km in the equatorial magnetosphere at around 9-12 RE. We compare growth rates and spatial scales with a range of theoretical predictions of magnetotail instabilities, including the cross-field current instability and the shear-flow ballooning instability. We conclude that, although the cross-field current instability can generate similar magnitude of growth rates, the range of unstable wavenumbers indicates that the shear-flow ballooning instability is the most likely explanation for our observations.