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.

Ensemble projections for wine production in the Douro Valley of Portugal

Wine production is largely governed by atmospheric conditions, such as air temperature and precipitation, together with soil management and viticultural/enological practices. Therefore, anthropogenic climate change is likely to have important impacts on the winemaking sector worldwide. An important winemaking region is the Portuguese Douro Valley, which is known by its world-famous Port Wine. The identification of robust relationships between atmospheric factors and wine parameters is of great relevance for the region. A multivariate linear regression analysis of a long wine production series (1932–2010) reveals that high rainfall and cool temperatures during budburst, shoot and inflorescence development (February-March) and warm temperatures during flowering and berry development (May) are generally favourable to high production. The probabilities of occurrence of three production categories (low, normal and high) are also modelled using multinomial logistic regression. Results show that both statistical models are valuable tools for predicting the production in a given year with a lead time of 3–4 months prior to harvest. These statistical models are applied to an ensemble of 16 regional climate model experiments following the SRES A1B scenario to estimate possible future changes. Wine production is projected to increase by about 10 % by the end of the 21st century, while the occurrence of high production years is expected to increase from 25 % to over 60 %. Nevertheless, further model development will be needed to include other aspects that may shape production in the future. In particular, the rising heat stress and/or changes in ripening conditions could limit the projected production increase in future decades.

Quantifying the climatological relationship between extratropical cyclone intensity and atmospheric precursors

We introduce a novel technique in which linear regression analysis is applied to clusters of tracked cyclones to statistically assess the factors controlling cyclone development. We illustrate this technique by evaluating the differences between cyclones forming in the west and east North Atlantic (herein termed west and east Atlantic cyclones). Enhanced cyclone intensity 2 days after genesis is found to be associated with deeper upper-level troughs upstream of the cyclone center at the genesis time in both west and east Atlantic cyclones. However, whilst west Atlantic cyclones are also enhanced by the presence of strong fronts, east Atlantic cyclones are not. Instead, east Atlantic cyclones exhibit an enhancement when diabatically generated midlevel potential vorticity is present (with the enhancement being of approximately equal magnitude to that associated with the potential vorticity in the upper-level trough). This is consistent with the paradigm of latent heat release in the warm conveyor belt region playing an important role in the development of east Atlantic cyclones.

The stability of a two-dimensional vorticity filament under uniform strain

The quantitative effects of uniform strain and background rotation on the stability of a strip of constant vorticity (a simple shear layer) are examined. The thickness of the strip decreases in time under the strain, so it is necessary to formulate the linear stability analysis for a time-dependent basic flow. The results show that even a strain rate γ (scaled with the vorticity of the strip) as small as 0.25 suppresses the conventional Rayleigh shear instability mechanism, in the sense that the r.m.s. wave steepness cannot amplify by more than a certain factor, and must eventually decay. For γ < 0.25 the amplification factor increases as γ decreases; however, it is only 3 when γ e 0.065. Numerical simulations confirm the predictions of linear theory at small steepness and predict a threshold value necessary for the formation of coherent vortices. The results help to explain the impression from numerous simulations of two-dimensional turbulence reported in the literature that filaments of vorticity infrequently roll up into vortices. The stabilization effect may be expected to extend to two- and three-dimensional quasi-geostrophic flows.

Influence of dietary protein intake and glycemic index on the association between TCF7L2 HapA and weight gain

BACKGROUND: Genetic polymorphisms of transcription factor 7-like 2 (TCF7L2) have been associated with type 2 diabetes and BMI. OBJECTIVE: The objective was to investigate whether TCF7L2 HapA is associated with weight development and whether such an association is modulated by protein intake or by the glycemic index (GI). DESIGN: The investigation was based on prospective data from 5 cohort studies nested within the European Prospective Investigation into Cancer and Nutrition. Weight change was followed up for a mean (±SD) of 6.8 ± 2.5 y. TCF7L2 rs7903146 and rs10885406 were successfully genotyped in 11,069 individuals and used to derive HapA. Multiple logistic and linear regression analysis was applied to test for the main effect of HapA and its interaction with dietary protein or GI. Analyses from the cohorts were combined by random-effects meta-analysis. RESULTS: HapA was associated neither with baseline BMI (0.03 ± 0.07 BMI units per allele; P = 0.6) nor with annual weight change (8.8 ± 11.7 g/y per allele; P = 0.5). However, a previously shown positive association between intake of protein, particularly of animal origin, and subsequent weight change in this population proved to be attenuated by TCF7L2 HapA (P-interaction = 0.01). We showed that weight gain becomes independent of protein intake with an increasing number of HapA alleles. Substitution of protein with either fat or carbohydrates showed the same effects. No interaction with GI was observed. CONCLUSION: TCF7L2 HapA attenuates the positive association between animal protein intake and long-term body weight change in middle-aged Europeans but does not interact with the GI of the diet.

Association of low adiponectin levels with the metabolic syndrome--the Chennai Urban Rural Epidemiology Study (CURES-4)

The aim of the study was to assess the relation of adiponectin levels with the metabolic syndrome in Asian Indians, a high-risk group for diabetes and premature coronary artery disease. The study was conducted on 100 (50 men and 50 women) type 2 diabetic subjects and 100 age and sex matched subjects with normal glucose tolerance selected from the Chennai Urban Rural Epidemiology Study, an ongoing population study in Chennai in southern India. Metabolic syndrome was defined using modified Adult Treatment Panel III (ATPIII) guidelines. Adiponectin values were significantly lower in diabetic subjects (men: 5.2 vs 8.3 microg/mL, P=.00l; women: 7.6 vs 11.1 microg/mL, P<.00l) and those with the metabolic syndrome (men: 5.0 vs 6.8 microg/mL, P=.01; women: 6.5 vs 9.9 microg/mL, P=.001) compared with those without. Linear regression analysis revealed adiponectin to be associated with body mass index (P<.05), waist circumference (P<.01), fasting plasma glucose (P=.001), glycated hemoglobin (P<.001), triglycerides (P<.00l), high-density lipoprotein (HDL) cholesterol (P<.001), cholesterol/HDL ratio (P<.00l), and insulin resistance measured by homeostasis assessment model (P<.00l). Factor analysis identified 2 factors: factor 1, negatively loaded with adiponectin and HDL cholesterol and positively loaded with triglycerides, waist circumference, and insulin resistance measured by homeostasis assessment model; and factor 2, with a positive loading of waist circumference and systolic and diastolic blood pressure. Logistic regression analysis revealed adiponectin to be negatively associated with metabolic syndrome (odds ratio [OR], 0.365; P<.001) even after adjusting for age (OR, 0.344; P<.00l), sex (OR, 0.293; P<.001), and body mass index (OR, 0.292; P<.00l). Lower adiponectin levels are associated with the metabolic syndrome per se and several of its components, particularly, diabetes, insulin resistance, and dyslipidemia in this urban south Indian population.

Multiple stationary solutions of an irradiated slab

A mathematical model describing the heat budget of an irradiated medium is introduced. The one-dimensional form of the equations and boundary conditions are presented and analysed. Heat transport at one face of the slab occurs by absorption (and reflection) of an incoming beam of short-wave radiation with a fraction of this radiation penetrating into the body of the slab, a diffusive heat flux in the slab and a prescribed incoming heat flux term. The other face of the slab is immersed in its own melt and is considered to be a free surface. Here, temperature continuity is prescribed and evolution of the surface is determined by a Stefan condition. These boundary conditions are flexible enough to describe a range of situations such as a laser shining on an opaque medium, or the natural environment of polar sea ice or lake ice. A two-stream radiation model is used which replaces the simple Beer’s law of radiation attenuation frequently used for semi-infinite domains. The stationary solutions of the governing equations are sought and it is found that there exists two possible stationary solutions for a given set of boundary conditions and a range of parameter choices. It is found that the existence of two stationary solutions is a direct result of the model of radiation absorption, due to its effect on the albedo of the medium. A linear stability analysis and numerical calculations indicate that where two stationary solutions exist, the solution corresponding to a larger thickness is always stable and the solution corresponding to a smaller thickness is unstable. Numerical simulations reveal that when there are two solutions, if the slab is thinner than the smaller stationary thickness it will melt completely, whereas if the slab is thicker than the smaller stationary thickness it will evolve toward the larger stationary thickness. These results indicate that other mechanisms (e.g. wave-induced agglomeration of crystals) are necessary to grow a slab from zero initial thickness in the parameter regime that yields two stationary solutions.

A mathematical model of crystallization in an emulsion

A mathematical model incorporating many of the important processes at work in the crystallization of emulsions is presented. The model describes nucleation within the discontinuous domain of an emulsion, precipitation in the continuous domain, transport of monomers between the two domains, and formation and subsequent growth of crystals in both domains. The model is formulated as an autonomous system of nonlinear, coupled ordinary differential equations. The description of nucleation and precipitation is based upon the Becker–Döring equations of classical nucleation theory. A particular feature of the model is that the number of particles of all species present is explicitly conserved; this differs from work that employs Arrhenius descriptions of nucleation rate. Since the model includes many physical effects, it is analyzed in stages so that the role of each process may be understood. When precipitation occurs in the continuous domain, the concentration of monomers falls below the equilibrium concentration at the surface of the drops of the discontinuous domain. This leads to a transport of monomers from the drops into the continuous domain that are then incorporated into crystals and nuclei. Since the formation of crystals is irreversible and their subsequent growth inevitable, crystals forming in the continuous domain effectively act as a sink for monomers “sucking” monomers from the drops. In this case, numerical calculations are presented which are consistent with experimental observations. In the case in which critical crystal formation does not occur, the stationary solution is found and a linear stability analysis is performed. Bifurcation diagrams describing the loci of stationary solutions, which may be multiple, are numerically calculated.

Variation in genes related to hepatic lipid metabolism and changes in waist circumference and body weight.

We analysed single nucleotide polymorphisms (SNPs) tagging the genetic variability of six candidate genes (ATF6, FABP1, LPIN2, LPIN3, MLXIPL and MTTP) involved in the regulation of hepatic lipid metabolism, an important regulatory site of energy balance for associations with body mass index (BMI) and changes in weight and waist circumference. We also investigated effect modification by sex and dietary intake. Data of 6,287 individuals participating in the European prospective investigation into cancer and nutrition were included in the analyses. Data on weight and waist circumference were followed up for 6.9 ± 2.5 years. Association of 69 tagSNPs with baseline BMI and annual changes in weight as well as waist circumference were investigated using linear regression analysis. Interactions with sex, GI and intake of carbohydrates, fat as well as saturated, monounsaturated and polyunsaturated fatty acids were examined by including multiplicative SNP-covariate terms into the regression model. Neither baseline BMI nor annual weight or waist circumference changes were significantly associated with variation in the selected genes in the entire study population after correction for multiple testing. One SNP (rs1164) in LPIN2 appeared to be significantly interacting with sex (p = 0.0003) and was associated with greater annual weight gain in men (56.8 ± 23.7 g/year per allele, p = 0.02) than in women (-25.5 ± 19.8 g/year per allele, p = 0.2). With respect to gene-nutrient interaction, we could not detect any significant interactions when accounting for multiple testing. Therefore, out of our six candidate genes, LPIN2 may be considered as a candidate for further studies.

Novel single trial movement classification based on temporal dynamics of EEG

Various complex oscillatory processes are involved in the generation of the motor command. The temporal dynamics of these processes were studied for movement detection from single trial electroencephalogram (EEG). Autocorrelation analysis was performed on the EEG signals to find robust markers of movement detection. The evolution of the autocorrelation function was characterised via the relaxation time of the autocorrelation by exponential curve fitting. It was observed that the decay constant of the exponential curve increased during movement, indicating that the autocorrelation function decays slowly during motor execution. Significant differences were observed between movement and no moment tasks. Additionally, a linear discriminant analysis (LDA) classifier was used to identify movement trials with a peak accuracy of 74%.

Contrasting interannual and multidecadal NAO variability

Decadal and longer timescale variability in the winter North Atlantic Oscillation (NAO) has considerable impact on regional climate, yet it remains unclear what fraction of this variability is potentially predictable. This study takes a new approach to this question by demonstrating clear physical differences between NAO variability on interannual-decadal (<30 year) and multidecadal (>30 year) timescales. It is shown that on the shorter timescale the NAO is dominated by variations in the latitude of the North Atlantic jet and storm track, whereas on the longer timescale it represents changes in their strength instead. NAO variability on the two timescales is associated with different dynamical behaviour in terms of eddy-mean flow interaction, Rossby wave breaking and blocking. The two timescales also exhibit different regional impacts on temperature and precipitation and different relationships to sea surface temperatures. These results are derived from linear regression analysis of the Twentieth Century and NCEP-NCAR reanalyses and of a high-resolution HiGEM General Circulation Model control simulation, with additional analysis of a long sea level pressure reconstruction. Evidence is presented for an influence of the ocean circulation on the longer timescale variability of the NAO, which is particularly clear in the model data. As well as providing new evidence of potential predictability, these findings are shown to have implications for the reconstruction and interpretation of long climate records.

Plasma free fatty acids do not provide the link between obesity and insulin resistance or β-cell dysfunction: results of the Reading, Imperial, Surrey, Cambridge, Kings (RISCK) study

Aims To investigate the relationship between adiposity and plasma free fatty acid levels and the influence of total plasma free fatty acid level on insulin sensitivity and β-cell function. Methods An insulin sensitivity index, acute insulin response to glucose and a disposition index, derived from i.v. glucose tolerance minimal model analysis and total fasting plasma free fatty acid levels were available for 533 participants in the Reading, Imperial, Surrey, Cambridge, Kings study. Bivariate correlations were made between insulin sensitivity index, acute insulin response to glucose and disposition index and both adiposity measures (BMI, waist circumference and body fat mass) and total plasma free fatty acid levels. Multivariate linear regression analysis was performed, controlling for age, sex, ethnicity and adiposity. Results After adjustment, all adiposity measures were inversely associated with insulin sensitivity index (BMI: β = −0.357; waist circumference: β = −0.380; body fat mass: β = −0.375) and disposition index (BMI: β = −0.215; waist circumference: β = −0.248; body fat mass: β = −0.221) and positively associated with acute insulin response to glucose [BMI: β = 0.200; waist circumference: β = 0.195; body fat mass β = 0.209 (P values <0.001)]. Adiposity explained 13, 4 and 5% of the variation in insulin sensitivity index, acute insulin response to glucose and disposition index, respectively. After adjustment, no adiposity measure was associated with free fatty acid level, but total plasma free fatty acid level was inversely associated with insulin sensitivity index (β = −0.133), acute insulin response to glucose (β = −0.148) and disposition index [β = −0.218 (P values <0.01)]. Plasma free fatty acid concentration accounted for 1.5, 2 and 4% of the variation in insulin sensitivity index, acute insulin response to glucose and disposition index, respectively. Conclusions Plasma free fatty acid levels have a modest negative association with insulin sensitivity, β-cell secretion and disposition index but no association with adiposity measures. It is unlikely that plasma free fatty acids are the primary mediators of obesity-related insulin resistance or β-cell dysfunction.

An event-based approach to understanding decadal fluctuations in the Atlantic meridional overturning circulation

Many previous studies have shown that unforced climate model simulations exhibit decadal-scale fluctuations in the Atlantic meridional overturning circulation (AMOC), and that this variability can have impacts on surface climate fields. However, the robustness of these surface fingerprints across different models is less clear. Furthermore, with the potential for coupled feedbacks that may amplify or damp the response, it is not known whether the associated climate signals are linearly related to the strength of the AMOC changes, or if the fluctuation events exhibit nonlinear behaviour with respect to their strength or polarity. To explore these questions, we introduce an objective and flexible method for identifying the largest natural AMOC fluctuation events in multicentennial/multimillennial simulations of a variety of coupled climate models. The characteristics of the events are explored, including their magnitude, meridional coherence and spatial structure, as well as links with ocean heat transport and the horizontal circulation. The surface fingerprints in ocean temperature and salinity are examined, and compared with the results of linear regression analysis. It is found that the regressions generally provide a good indication of the surface changes associated with the largest AMOC events. However, there are some exceptions, including a nonlinear change in the atmospheric pressure signal, particularly at high latitudes, in HadCM3. Some asymmetries are also found between the changes associated with positive and negative AMOC events in the same model. Composite analysis suggests that there are signals that are robust across the largest AMOC events in each model, which provides reassurance that the surface changes associated with one particular event will be similar to those expected from regression analysis. However, large differences are found between the AMOC fingerprints in different models, which may hinder the prediction and attribution of such events in reality.

Motor imagery-induced EEG patterns in individuals with spinal cord injury and their impact on brain–computer interface accuracy

Objective. Assimilating the diagnosis complete spinal cord injury (SCI) takes time and is not easy, as patients know that there is no ‘cure’ at the present time. Brain–computer interfaces (BCIs) can facilitate daily living. However, inter-subject variability demands measurements with potential user groups and an understanding of how they differ to healthy users BCIs are more commonly tested with. Thus, a three-class motor imagery (MI) screening (left hand, right hand, feet) was performed with a group of 10 able-bodied and 16 complete spinal-cord-injured people (paraplegics, tetraplegics) with the objective of determining what differences were present between the user groups and how they would impact upon the ability of these user groups to interact with a BCI. Approach. Electrophysiological differences between patient groups and healthy users are measured in terms of sensorimotor rhythm deflections from baseline during MI, electroencephalogram microstate scalp maps and strengths of inter-channel phase synchronization. Additionally, using a common spatial pattern algorithm and a linear discriminant analysis classifier, the classification accuracy was calculated and compared between groups. Main results. It is seen that both patient groups (tetraplegic and paraplegic) have some significant differences in event-related desynchronization strengths, exhibit significant increases in synchronization and reach significantly lower accuracies (mean (M) = 66.1%) than the group of healthy subjects (M = 85.1%). Significance. The results demonstrate significant differences in electrophysiological correlates of motor control between healthy individuals and those individuals who stand to benefit most from BCI technology (individuals with SCI). They highlight the difficulty in directly translating results from healthy subjects to participants with SCI and the challenges that, therefore, arise in providing BCIs to such individuals

Motor imagery induced EEG patterns in spinal cord injury patients and their impact on Brain-Computer Interface accuracy

OBJECTIVE: Assimilating the diagnosis complete spinal cord injury (SCI) takes time and is not easy, as patients know that there is no 'cure' at the present time. Brain-computer interfaces (BCIs) can facilitate daily living. However, inter-subject variability demands measurements with potential user groups and an understanding of how they differ to healthy users BCIs are more commonly tested with. Thus, a three-class motor imagery (MI) screening (left hand, right hand, feet) was performed with a group of 10 able-bodied and 16 complete spinal-cord-injured people (paraplegics, tetraplegics) with the objective of determining what differences were present between the user groups and how they would impact upon the ability of these user groups to interact with a BCI. APPROACH: Electrophysiological differences between patient groups and healthy users are measured in terms of sensorimotor rhythm deflections from baseline during MI, electroencephalogram microstate scalp maps and strengths of inter-channel phase synchronization. Additionally, using a common spatial pattern algorithm and a linear discriminant analysis classifier, the classification accuracy was calculated and compared between groups. MAIN RESULTS: It is seen that both patient groups (tetraplegic and paraplegic) have some significant differences in event-related desynchronization strengths, exhibit significant increases in synchronization and reach significantly lower accuracies (mean (M) = 66.1%) than the group of healthy subjects (M = 85.1%). SIGNIFICANCE: The results demonstrate significant differences in electrophysiological correlates of motor control between healthy individuals and those individuals who stand to benefit most from BCI technology (individuals with SCI). They highlight the difficulty in directly translating results from healthy subjects to participants with SCI and the challenges that, therefore, arise in providing BCIs to such individuals.