Metabifurcation analysis of a mean field model of the cortex

Mean field models (MFMs) of cortical tissue incorporate salient, average features of neural masses in order to model activity at the population level, thereby linking microscopic physiology to macroscopic observations, e.g., with the electroencephalogram (EEG). One of the common aspects of MFM descriptions is the presence of a high-dimensional parameter space capturing neurobiological attributes deemed relevant to the brain dynamics of interest. We study the physiological parameter space of a MFM of electrocortical activity and discover robust correlations between physiological attributes of the model cortex and its dynamical features. These correlations are revealed by the study of bifurcation plots, which show that the model responses to changes in inhibition belong to two archetypal categories or “families”. After investigating and characterizing them in depth, we discuss their essential differences in terms of four important aspects: power responses with respect to the modeled action of anesthetics, reaction to exogenous stimuli such as thalamic input, and distributions of model parameters and oscillatory repertoires when inhibition is enhanced. Furthermore, while the complexity of sustained periodic orbits differs significantly between families, we are able to show how metamorphoses between the families can be brought about by exogenous stimuli. We here unveil links between measurable physiological attributes of the brain and dynamical patterns that are not accessible by linear methods. They instead emerge when the nonlinear structure of parameter space is partitioned according to bifurcation responses. We call this general method “metabifurcation analysis”. The partitioning cannot be achieved by the investigation of only a small number of parameter sets and is instead the result of an automated bifurcation analysis of a representative sample of 73,454 physiologically admissible parameter sets. Our approach generalizes straightforwardly and is well suited to probing the dynamics of other models with large and complex parameter spaces.

Realistic mean field forward predictions for the integration of co-registered EEG/fMRI

Brain activity can be measured non-invasively with functional imaging techniques. Each pixel in such an image represents a neural mass of about 105 to 107 neurons. Mean field models (MFMs) approximate their activity by averaging out neural variability while retaining salient underlying features, like neurotransmitter kinetics. However, MFMs incorporating the regional variability, realistic geometry and connectivity of cortex have so far appeared intractable. This lack of biological realism has led to a focus on gross temporal features of the EEG. We address these impediments and showcase a "proof of principle" forward prediction of co-registered EEG/fMRI for a full-size human cortex in a realistic head model with anatomical connectivity, see figure 1. MFMs usually assume homogeneous neural masses, isotropic long-range connectivity and simplistic signal expression to allow rapid computation with partial differential equations. But these approximations are insufficient in particular for the high spatial resolution obtained with fMRI, since different cortical areas vary in their architectonic and dynamical properties, have complex connectivity, and can contribute non-trivially to the measured signal. Our code instead supports the local variation of model parameters and freely chosen connectivity for many thousand triangulation nodes spanning a cortical surface extracted from structural MRI. This allows the introduction of realistic anatomical and physiological parameters for cortical areas and their connectivity, including both intra- and inter-area connections. Proper cortical folding and conduction through a realistic head model is then added to obtain accurate signal expression for a comparison to experimental data. To showcase the synergy of these computational developments, we predict simultaneously EEG and fMRI BOLD responses by adding an established model for neurovascular coupling and convolving "Balloon-Windkessel" hemodynamics. We also incorporate regional connectivity extracted from the CoCoMac database [1]. Importantly, these extensions can be easily adapted according to future insights and data. Furthermore, while our own simulation is based on one specific MFM [2], the computational framework is general and can be applied to models favored by the user. Finally, we provide a brief outlook on improving the integration of multi-modal imaging data through iterative fits of a single underlying MFM in this realistic simulation framework.

A high frequency $hp$ boundary element method for scattering by convex polygons

In this paper we propose and analyze a hybrid $hp$ boundary element method for the solution of problems of high frequency acoustic scattering by sound-soft convex polygons, in which the approximation space is enriched with oscillatory basis functions which efficiently capture the high frequency asymptotics of the solution. We demonstrate, both theoretically and via numerical examples, exponential convergence with respect to the order of the polynomials, moreover providing rigorous error estimates for our approximations to the solution and to the far field pattern, in which the dependence on the frequency of all constants is explicit. Importantly, these estimates prove that, to achieve any desired accuracy in the computation of these quantities, it is sufficient to increase the number of degrees of freedom in proportion to the logarithm of the frequency as the frequency increases, in contrast to the at least linear growth required by conventional methods.

Trade-offs and responsiveness of the single-layer urban canopy parametrization in WRF: an offline evaluation using the MOSCEM optimization algorithm and field observations

For an increasing number of applications, mesoscale modelling systems now aim to better represent urban areas. The complexity of processes resolved by urban parametrization schemes varies with the application. The concept of fitness-for-purpose is therefore critical for both the choice of parametrizations and the way in which the scheme should be evaluated. A systematic and objective model response analysis procedure (Multiobjective Shuffled Complex Evolution Metropolis (MOSCEM) algorithm) is used to assess the fitness of the single-layer urban canopy parametrization implemented in the Weather Research and Forecasting (WRF) model. The scheme is evaluated regarding its ability to simulate observed surface energy fluxes and the sensitivity to input parameters. Recent amendments are described, focussing on features which improve its applicability to numerical weather prediction, such as a reduced and physically more meaningful list of input parameters. The study shows a high sensitivity of the scheme to parameters characterizing roof properties in contrast to a low response to road-related ones. Problems in partitioning of energy between turbulent sensible and latent heat fluxes are also emphasized. Some initial guidelines to prioritize efforts to obtain urban land-cover class characteristics in WRF are provided. Copyright © 2010 Royal Meteorological Society and Crown Copyright.

Effect of Rht alleles on the tolerance of wheat grain set to high temperature and drought stress during booting and anthesis

Factorial pot experiments were conducted to compare the responses of GA-sensitive and GA-insensitive reduced height (Rht) alleles in wheat for susceptibility to heat and drought stress during booting and anthesis. Grain set (grains/spikelet) of near isogenic lines (NILs) was assessed following three day transfers to controlled environments imposing day temperatures (t) from 20 to 40°C. Transfers were during booting and/or anthesis and pots maintained at field capacity (FC) or had water withheld. Logistic responses (y = c/1+e-b(t -m)) described declining grain set with increasing t, and t5 was that fitted to give a 5% reduction in grain set. Averaged over NIL, t5 for anthesis at FC was 31.7±0.47°C (S.E.M, 26 d.f.). Drought at anthesis reduced t5 by <2°C. Maintaining FC at booting conferred considerable resistance to high temperatures (t5=33.9°C) but booting was particularly heat susceptible without water (t5 =26.5°C). In one background (cv. Mercia), for NILs varying at the Rht-D1 locus, there was progressive reduction in t5 with dwarfing and reduced gibberellic acid (GA) sensitivity (Rht-D1a, tall, 32.7±0.72; Rht-D1b, semi-dwarf, 29.5±0.85; Rht-D1c, severe dwarf, 24.2±0.72). This trend was not evident for the Rht-B1 locus, or for Rht-D1b in an alternative background (Maris Widgeon). The GA-sensitive severe dwarf Rht12 was more heat tolerant (t5=29.4±0.72) than the similarly statured GA-insensitive Rht-D1c. The GA-sensitive, semi-dwarfing Rht8 conferred greater drought tolerance in one experiment. Despite the effects of Rht-D1 alleles in Mercia on stress tolerance, the inconsistency of the effects over background and locus led to the conclusion that semi-dwarfing with GA-insensitivity did not necessarily increase sensitivity to stress at booting and flowering. In comparison to effects of semi-dwarfing alleles, responses to heat stress are much more dramatically affected by water availability and the precise growth stage at which the stress is experienced by the plants.

Flow structure and near-field dispersion in arrays of building-like obstacles

Dispersion in the near-field region of localised releases in urban areas is difficult to predict because of the strong influence of individual buildings. Effects include upstream dispersion, trapping of material into building wakes and enhanced concentration fluctuations. As a result, concentration patterns are highly variable in time and mean profiles in the near field are strongly non-Gaussian. These aspects of near-field dispersion are documented by analysing data from direct numerical simulations in arrays of building-like obstacles and are related to the underlying flow structure. The mean flow structure around the buildings is found to exert a strong influence over the dispersion of material in the near field. Diverging streamlines around buildings enhance lateral dispersion. Entrainment of material into building wakes in the very near field gives rise to secondary sources, which then affect the subsequent dispersion pattern. High levels of concentration fluctuations are also found in this very near field; the fluctuation intensity is of order 2 to 5.

Use of virulence determinants and seropathotypes to distinguish high- and low-risk Escherichia coli O157 and non-O157 isolates from Europe

The presence of 10 virulence genes was examined using polymerase chain reaction (PCR) in 365 European O157 and non-O157 Escherichia coli isolates associated with verotoxin production. Strain-specific PCR data were analysed using hierarchical clustering. The resulting dendrogram clearly separated O157 from non-O157 strains. The former clustered typical high-risk seropathotype (SPT) A strains from all regions, including Sweden and Spain, which were homogenous by Cramer's V statistic, and strains with less typical O157 features mostly from Hungary. The non-O157 strains divided into a high-risk SPTB harbouring O26, O111 and O103 strains, a group pathogenic to pigs, and a group with few virulence genes other than for verotoxin. The data demonstrate SPT designation and selected PCR separated verotoxigenic E. coli of high and low risk to humans; although more virulence genes or pulsed-field gel electrophoresis will need to be included to separate high-risk strains further for epidemiological tracing.

RMetS special interest group meeting: high resolution data assimilation

Data assimilation (DA) systems are evolving to meet the demands of convection-permitting models in the field of weather forecasting. On 19 April 2013 a special interest group meeting of the Royal Meteorological Society brought together UK researchers looking at different aspects of the data assimilation problem at high resolution, from theory to applications, and researchers creating our future high resolution observational networks. The meeting was chaired by Dr Sarah Dance of the University of Reading and Dr Cristina Charlton-Perez from the MetOffice@Reading. The purpose of the meeting was to help define the current state of high resolution data assimilation in the UK. The workshop assembled three main types of scientists: observational network specialists, operational numerical weather prediction researchers and those developing the fundamental mathematical theory behind data assimilation and the underlying models. These three working areas are intrinsically linked; therefore, a holistic view must be taken when discussing the potential to make advances in high resolution data assimilation.

Field line resonances as a trigger and a tracer for substorm onset

In this paper, we show that periodic auroral arc structures are seen at the location of one particular auroral substorm onset for the 15 min preceding onset, suggesting that field line resonances should be considered a strong candidate for triggering substorm onset. Irrespective of whether this field line resonance is coincidentally or causally linked to this substorm onset, the characteristics of the field line resonance can be used to remote sense the characteristics of the geomagnetic field line that supports substorm onset. In this instance, the eigenfrequency of this resonance is around 12 mHz. Interestingly, however, there is no evidence of this field line resonance in a seven satellite major Time History of Events and Macroscale Interactions during Substorms (THEMIS)-GOES conjunction, ranging from geosynchronous orbit to ~30 RE. However, using space-based cross-phase measurements of the local field line eigenfrequency at the inner THEMIS locations, we find that the local field line eigenfrequency is 6–10 mHz. Hence, we can reliably say that this 12 mHz Field Line Resonance (FLR) must lie inside of THEMIS locations. Our conclusion is that a high-m field line resonance can both represent a strong candidate for a trigger for substorm onset, as first proposed by Samson et al. (1992), and that its characteristics can provide invaluable information as to where substorm onset occurs in the magnetosphere.

Modulation of UK lightning by heliospheric magnetic field polarity

Observational studies have reported solar magnetic modulation of terrestrial lightning on a range of time scales, from days to decades. The proposed mechanism is two-step: lightning rates vary with galactic cosmic ray (GCR) flux incident on Earth, either via changes in atmospheric conductivity and/or direct triggering of lightning. GCR flux is, in turn, primarily controlled by the heliospheric magnetic field (HMF) intensity. Consequently, global changes in lightning rates are expected. This study instead considers HMF polarity, which doesnʼt greatly affect total GCR flux. Opposing HMF polarities are, however, associated with a 40–60% difference in observed UK lightning and thunder rates. As HMF polarity skews the terrestrial magnetosphere from its nominal position, this perturbs local ionospheric potential at high latitudes and local exposure to energetic charged particles from the magnetosphere. We speculate as to the mechanism(s) by which this may, in turn, redistribute the global location and/or intensity of thunderstorm activity.

Comparative Cluster/Double Star observations of the high and low latitude dayside magnetopause

The launch of the Double Star mission has provided the opportunity to monitor events at distinct locations on the dayside magnetopause, in coordination with the quartet of Cluster spacecraft. We present results of two such coordinated studies. In the first, 6 April 2004, both Cluster and the Double Star TC-1 spacecraft were on outbound transits through the dawn-side magnetosphere. Cluster observed northward moving FTEs with +/- polarity, whereas TC-1 saw -/+ polarity FTEs. The strength, motion and occurrence of the FTE signatures changes somewhat according to changes in IMF clock angle. These observations are consistent with ongoing reconnection on the dayside magnetopause, resulting in a series of flux transfer events (FTEs) seen both at Cluster and TC-1. The observed polarity and motion of each FTE signature advocates the existence of an active reconnection region consistently located between the positions of Cluster and TC-1, lying north and south of the reconnection line, respectively. This scenario is supported by the application of a model, designed to track flux tube motion, to conditions appropriate for the prevailing interplanetary conditions. The results from the model confirm the observational evidence that the low-latitude FTE dynamics is sensitive to changes in convected upstream conditions. In particular, changing the interplanetary magnetic field (IMF) clock angle in the model predicts that TC-1 should miss the resulting FTEs more often than Cluster, as is observed. For the second conjunction, on the 4 Jan 2005, the Cluster and TC-1 spacecraft all exited the dusk-side magnetosphere almost simultaneously, with TC-1 lying almost equatorial and Cluster at northern latitudes at about 4 RE from TC-1. The spacecraft traverse the magnetopause during a strong reversal in the IMF from northward to southward and a number of magnetosheath FTE signatures are subsequently observed. One coordinated FTE, studied in detail by Pu et al, [this issue], carries and inflowing energetic electron population and shows a motion and orientation which is similar at all spacecraft and consistent with the predictions of the model for the flux tube dynamics, given a near sub-solar reconnection line. This event can be interpreted either as the passage of two parallel flux tubes arising from adjacent x-line positions, or as a crossing of a single flux tube at different positions.

A comparison between ion characteristics observed by the POLAR and DMSP spacecraft in the high-latitude magnetosphere

We study here the injection and transport of ions in the convection-dominated region of the Earth’s magnetosphere. The total ion counts from the CAMMICE MICS instrument aboard the POLAR spacecraft are used to generate occurrence probability distributions of magnetospheric ion populations. MICS ion spectra are characterised by both the peak in the differential energy flux, and the average energy of ions striking the detector. The former permits a comparison with the Stubbs et al. (2001) survey of He2+ ions of solar wind origin within the magnetosphere. The latter can address the occurrences of various classifications of precipitating particle fluxes observed in the topside ionosphere by DMSP satellites (Newell and Meng, 1992). The peak energy occurrences are consistent with our earlier work, including the dawn-dusk asymmetry with enhanced occurrences on the dawn flank at low energies, switching to the dusk flank at higher energies. The differences in the ion energies observed in these two studies can be explained by drift orbit effects and acceleration processes at the magnetopause, and in the tail current sheet. Near noon at average ion energies of _1 keV, the cusp and open LLBL occur further poleward here than in the Newell and Meng survey, probably due to convectionrelated time-of-flight effects. An important new result is that the pre-noon bias previously observed in the LLBL is most likely due to the component of this population on closed field lines, formed largely by low energy ions drifting earthward from the tail. There is no evidence here of mass and momentum transfer from the solar wind to the LLBL by nonreconnection coupling. At higher energies (_2–20 keV), we observe ions mapping to the auroral oval and can distinguish between the boundary and central plasma sheets. We show that ions at these energies relate to a transition from dawnward to duskward dominated flow, this is evidence of how ion drift orbits in the tail influence the location and behaviour of the plasma populations in the magnetosphere.

Open solar flux estimates from near-Earth measurements of the interplanetary magnetic field: comparison of the first two perihelion passes of the Ulysses spacecraft

Results from all phases of the orbits of the Ulysses spacecraft have shown that the magnitude of the radial component of the heliospheric field is approximately independent of heliographic latitude. This result allows the use of near- Earth observations to compute the total open flux of the Sun. For example, using satellite observations of the interplanetary magnetic field, the average open solar flux was shown to have risen by 29% between 1963 and 1987 and using the aa geomagnetic index it was found to have doubled during the 20th century. It is therefore important to assess fully the accuracy of the result and to check that it applies to all phases of the solar cycle. The first perihelion pass of the Ulysses spacecraft was close to sunspot minimum, and recent data from the second perihelion pass show that the result also holds at solar maximum. The high level of correlation between the open flux derived from the various methods strongly supports the Ulysses discovery that the radial field component is independent of latitude. We show here that the errors introduced into open solar flux estimates by assuming that the heliospheric field’s radial component is independent of latitude are similar for the two passes and are of order 25% for daily values, falling to 5% for averaging timescales of 27 days or greater. We compare here the results of four methods for estimating the open solar flux with results from the first and second perehelion passes by Ulysses. We find that the errors are lowest (1–5% for averages over the entire perehelion passes lasting near 320 days), for near-Earth methods, based on either interplanetary magnetic field observations or the aa geomagnetic activity index. The corresponding errors for the Solanki et al. (2000) model are of the order of 9–15% and for the PFSS method, based on solar magnetograms, are of the order of 13–47%. The model of Solanki et al. is based on the continuity equation of open flux, and uses the sunspot number to quantify the rate of open flux emergence. It predicts that the average open solar flux has been decreasing since 1987, as Correspondence to: M. Lockwood (m.lockwood@rl.ac.uk) is observed in the variation of all the estimates of the open flux. This decline combines with the solar cycle variation to produce an open flux during the second (sunspot maximum) perihelion pass of Ulysses which is only slightly larger than that during the first (sunspot minimum) perihelion pass.

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.

Coherent EISCAT Svalbard Radar spectra from the dayside cusp/cleft and their implications for transient field-aligned currents

Naturally enhanced incoherent scatter spectra from the vicinity of the dayside cusp/cleft, interpreted as being due to plasma turbulence driven by short bursts of intense field-aligned current, are compared with high-resolution narrow-angle auroral images and meridian scanning photometer data. Enhanced spectra have been observed on many occasions in association with nightside aurora, but there has been only one report of such spectra seen in the cusp/cleft region. Narrow-angle images show considerable change in the aurora on timescales shorter than the 10-s radar integration period, which could explain spectra observed with both ion lines simultaneously enhanced. Enhanced radar spectra are generally seen inside or beside regions of 630-nm auroral emission, indicative of sharp F region conductivity gradients, but there appears also to be a correlation with dynamic, small-scale auroral forms of order 100 m and less in width.