Long duration stimuli and nonlinearities in the neural–haemodynamic coupling

Recent studies have shown that the haemodynamic responses to brief (<2 secs) stimuli can be well characterised as a linear convolution of neural activity with a suitable haemodynamic impulse response. In this paper, we show that the linear convolution model cannot predict measurements of blood flow responses to stimuli of longer duration (>2 secs), regardless of the impulse response function chosen. Modifying the linear convolution scheme to a nonlinear convolution scheme was found to provide a good prediction of the observed data. Whereas several studies have found a nonlinear coupling between stimulus input and blood flow responses, the current modelling scheme uses neural activity as an input, and thus implies nonlinearity in the coupling between neural activity and blood flow responses. Neural activity was assessed by current source density analysis of depth-resolved evoked field potentials, while blood flow responses were measured using laser Doppler flowmetry. All measurements were made in rat whisker barrel cortex after electrical stimulation of the whisker pad for 1 to 16 secs at 5 Hz and 1.2 mA (individual pulse width 0.3 ms).

The resting-state neurovascular coupling relationship: rapid changes in spontaneous neural activity in the somatosensory cortex are associated with haemodynamic fluctuations that resemble stimulus-evoked haemodynamics

Although promise exists for patterns of resting-state blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) brain connectivity to be used as biomarkers of early brain pathology, a full understanding of the nature of the relationship between neural activity and spontaneous fMRI BOLD fluctuations is required before such data can be correctly interpreted. To investigate this issue, we combined electrophysiological recordings of rapid changes in multi-laminar local field potentials from the somatosensory cortex of anaesthetized rats with concurrent two-dimensional optical imaging spectroscopy measurements of resting-state haemodynamics that underlie fluctuations in the BOLD fMRI signal. After neural ‘events’ were identified, their time points served to indicate the start of an epoch in the accompanying haemodynamic fluctuations. Multiple epochs for both neural ‘events’ and the accompanying haemodynamic fluctuations were averaged. We found that the averaged epochs of resting-state haemodynamic fluctuations taken after neural ‘events’ closely resembled the temporal profile of stimulus-evoked cortical haemodynamics. Furthermore, we were able to demonstrate that averaged epochs of resting-state haemodynamic fluctuations resembling the temporal profile of stimulus-evoked haemodynamics could also be found after peaks in neural activity filtered into specific electroencephalographic frequency bands (theta, alpha, beta, and gamma). This technique allows investigation of resting-state neurovascular coupling using methodologies that are directly comparable to that developed for investigating stimulus-evoked neurovascular responses.

Interhemispheric coupling, the West Antarctic Ice Sheet and warm Antarctic interglacials

Ice core evidence indicates that even though atmospheric CO2 concentrations did not exceed ~300 ppm at any point during the last 800 000 years, East Antarctica was at least ~3–4 °C warmer than preindustrial (CO2~280 ppm) in each of the last four interglacials. During the previous three interglacials, this anomalous warming was short lived (~3000 years) and apparently occurred before the completion of Northern Hemisphere deglaciation. Hereafter, we refer to these periods as "Warmer than Present Transients" (WPTs). We present a series of experiments to investigate the impact of deglacial meltwater on the Atlantic Meridional Overturning Circulation (AMOC) and Antarctic temperature. It is well known that a slowed AMOC would increase southern sea surface temperature (SST) through the bipolar seesaw and observational data suggests that the AMOC remained weak throughout the terminations preceding WPTs, strengthening rapidly at a time which coincides closely with peak Antarctic temperature. We present two 800 kyr transient simulations using the Intermediate Complexity model GENIE-1 which demonstrate that meltwater forcing generates transient southern warming that is consistent with the timing of WPTs, but is not sufficient (in this single parameterisation) to reproduce the magnitude of observed warmth. In order to investigate model and boundary condition uncertainty, we present three ensembles of transient GENIE-1 simulations across Termination II (135 000 to 124 000 BP) and three snapshot HadCM3 simulations at 130 000 BP. Only with consideration of the possible feedback of West Antarctic Ice Sheet (WAIS) retreat does it become possible to simulate the magnitude of observed warming.

Thermodynamics/dynamics coupling in weakly compressible turbulent stratified fluids

In traditional and geophysical fluid dynamics, it is common to describe stratified turbulent fluid flows with low Mach number and small relative density variations by means of the incompressible Boussinesq approximation. Although such an approximation is often interpreted as decoupling the thermodynamics from the dynamics, this paper reviews recent results and derive new ones that show that the reality is actually more subtle and complex when diabatic effects and a nonlinear equation of state are retained. Such an analysis reveals indeed: (1) that the compressible work of expansion/contraction remains of comparable importance as the mechanical energy conversions in contrast to what is usually assumed; (2) in a Boussinesq fluid, compressible effects occur in the guise of changes in gravitational potential energy due to density changes. This makes it possible to construct a fully consistent description of the thermodynamics of incompressible fluids for an arbitrary nonlinear equation of state; (3) rigorous methods based on using the available potential energy and potential enthalpy budgets can be used to quantify the work of expansion/contraction B in steady and transient flows, which reveals that B is predominantly controlled by molecular diffusive effects, and act as a significant sink of kinetic energy.

Eddy-driven exchange between the open ocean and a sub-ice shelf cavity

The exchange between the open ocean and sub-ice shelf cavities is important to both water mass transformations and ice shelf melting. Here we use a high-resolution (500 m) numerical model to investigate to which degree eddies produced by frontal instability at the edge of a polynya are capable of transporting dense High Salinity Shelf Water (HSSW) underneath an ice shelf. The applied surface buoyancy flux and ice shelf geometry is based on Ronne Ice Shelf in the southern Weddell Sea, an area of intense wintertime sea ice production where a flow of HSSW into the cavity has been observed. Results show that eddies are able to enter the cavity at the southwestern corner of the polynya where an anticyclonic rim current intersects the ice shelf front. The size and time scale of simulated eddies are in agreement with observations close to the Ronne Ice Front. The properties and strength of the inflow are sensitive to the prescribed total ice production, flushing the ice shelf cavity at a rate of 0.2–0.4 × 106 m3 s−1 depending on polynya size and magnitude of surface buoyancy flux. Eddy-driven HSSW transport into the cavity is reduced by about 50% if the model grid resolution is decreased to 2-5 km and eddies are not properly resolved.

The integrated WRF/urban modelling system: development, evaluation, and applications to urban environmental problems

To bridge the gaps between traditional mesoscale modelling and microscale modelling, the National Center for Atmospheric Research, in collaboration with other agencies and research groups, has developed an integrated urban modelling system coupled to the weather research and forecasting (WRF) model as a community tool to address urban environmental issues. The core of this WRF/urban modelling system consists of the following: (1) three methods with different degrees of freedom to parameterize urban surface processes, ranging from a simple bulk parameterization to a sophisticated multi-layer urban canopy model with an indoor–outdoor exchange sub-model that directly interacts with the atmospheric boundary layer, (2) coupling to fine-scale computational fluid dynamic Reynolds-averaged Navier–Stokes and Large-Eddy simulation models for transport and dispersion (T&D) applications, (3) procedures to incorporate high-resolution urban land use, building morphology, and anthropogenic heating data using the National Urban Database and Access Portal Tool (NUDAPT), and (4) an urbanized high-resolution land data assimilation system. This paper provides an overview of this modelling system; addresses the daunting challenges of initializing the coupled WRF/urban model and of specifying the potentially vast number of parameters required to execute the WRF/urban model; explores the model sensitivity to these urban parameters; and evaluates the ability of WRF/urban to capture urban heat islands, complex boundary-layer structures aloft, and urban plume T&D for several major metropolitan regions. Recent applications of this modelling system illustrate its promising utility, as a regional climate-modelling tool, to investigate impacts of future urbanization on regional meteorological conditions and on air quality under future climate change scenarios. Copyright © 2010 Royal Meteorological Society

Expression of the Rap1 guanine nucleotide exchange factor, MR-GEF, is altered in individuals with bipolar disorder

In the rodent forebrain GABAergic neurons are generated from progenitor cells that express the transcription factors Dlx1 and Dlx2. The Rap-1 guanine nucleotide exchange factor, MR-GEF, is turned on by many of these developing GABAergic neurons. Expression of both Dlx1/2 and MR-GEF is retained in both adult mouse and human forebrain where, in human, decreased Dlx1 expression has been associated with psychosis. Using in situ hybridization studies we show that MR-GEF expression is significantly down-regulated in the forebrain of Dlx1/2 double mutant mice suggesting that MR-GEF and Dlx1/2 form part of a common signalling pathway during GABAergic neuronal development. We therefore compared MR-GEF expression by in situ hybridization in individuals with major psychiatric disorders (schizophrenia, bipolar disorder, major depression) and control individuals. We observed a significant positive correlation between layers II and IV of the dorso-lateral prefrontal cortex (DLPFC) in the percentage of MR-GEF expressing neurons in individuals with bipolar disorder, but not in individuals with schizophrenia, major depressive disorder or in controls. Since MR-GEF encodes a Rap1 GEF able to activate G-protein signalling, we suggest that changes in MR-GEF expression could potentially influence neurotransmission.

Expression of the guanine nucleotide exchange factor, mr-gef, is regulated during the differentiation of specific subsets of telencephalic neurons

We have performed a screen combining subtractive hybridization with PCR to isolate genes that are regulated when neuroepithelial (NE) cells differentiate into neurons. From this screen, we have isolated a number of known genes that have not previously been associated with neurogenesis, together with several novel genes. Here we report that one of these genes, encoding a guanine nucleotide exchange factor (GEF), is regulated during the differentiation of distinct neuronal populations. We have cloned both rat and mouse GEF genes and shown that they are orthologs of the human gene, MR-GEF, which encodes a GEF that specifically activates the small GTPase, Rap1. We have therefore named the rat gene rat mr-gef (rmr-gef) and the mouse gene mouse mr-gef (mmr-gef). Here, we will collectively refer to these two rodent genes as mr-gef. Expression studies show that mr-gef is expressed by young neurons of the developing rodent CNS but not by progenitor cells in the ventricular zone (VZ). The expression pattern of mr-gef during early telencephalic neurogenesis is strikingly similar to that of GABA and the LIM homeobox gene Lhx6, a transcription factor expressed by GABAergic interneurons generated in the ventral telencephalon, some of which migrate into the cortex during development. These observations suggest that mr-gef encodes a protein that is part of a signaling pathway involved in telencephalic neurogenesis; particularly in the development of GABAergic interneurons.

Market transparency and collusion: on the UK agricultural tractor registration exchange

We review the decision by the European Commission in the case of the UK Agricultural Registration Exchange. We propose a theoretical model, offering a basis for some of the intuitive arguments used by the Commission on the anti-competitive role of information exchange in the case of price and non price collusion. Market transparency on non price data is shown to be a collusion facilitating device which may achieve stability in otherwise unstable cartels.

Coupling ensemble weather predictions based on TIGGE database with Grid-Xinanjiang model for flood forecast

The incorporation of numerical weather predictions (NWP) into a flood forecasting system can increase forecast lead times from a few hours to a few days. A single NWP forecast from a single forecast centre, however, is insufficient as it involves considerable non-predictable uncertainties and lead to a high number of false alarms. The availability of global ensemble numerical weather prediction systems through the THORPEX Interactive Grand Global Ensemble' (TIGGE) offers a new opportunity for flood forecast. The Grid-Xinanjiang distributed hydrological model, which is based on the Xinanjiang model theory and the topographical information of each grid cell extracted from the Digital Elevation Model (DEM), is coupled with ensemble weather predictions based on the TIGGE database (CMC, CMA, ECWMF, UKMO, NCEP) for flood forecast. This paper presents a case study using the coupled flood forecasting model on the Xixian catchment (a drainage area of 8826 km2) located in Henan province, China. A probabilistic discharge is provided as the end product of flood forecast. Results show that the association of the Grid-Xinanjiang model and the TIGGE database gives a promising tool for an early warning of flood events several days ahead.

Evaluating forecasts from SETAR models of exchange rates

We consider the forecasting performance of two SETAR exchange rate models proposed by Kräger and Kugler [J. Int. Money Fin. 12 (1993) 195]. Assuming that the models are good approximations to the data generating process, we show that whether the non-linearities inherent in the data can be exploited to forecast better than a random walk depends on both how forecast accuracy is assessed and on the ‘state of nature’. Evaluation based on traditional measures, such as (root) mean squared forecast errors, may mask the superiority of the non-linear models. Generalized impulse response functions are also calculated as a means of portraying the asymmetric response to shocks implied by such models.

Quantile forecasts of daily exchange rate returns from forecasts of realized volatility

Quantile forecasts are central to risk management decisions because of the widespread use of Value-at-Risk. A quantile forecast is the product of two factors: the model used to forecast volatility, and the method of computing quantiles from the volatility forecasts. In this paper we calculate and evaluate quantile forecasts of the daily exchange rate returns of five currencies. The forecasting models that have been used in recent analyses of the predictability of daily realized volatility permit a comparison of the predictive power of different measures of intraday variation and intraday returns in forecasting exchange rate variability. The methods of computing quantile forecasts include making distributional assumptions for future daily returns as well as using the empirical distribution of predicted standardized returns with both rolling and recursive samples. Our main findings are that the Heterogenous Autoregressive model provides more accurate volatility and quantile forecasts for currencies which experience shifts in volatility, such as the Canadian dollar, and that the use of the empirical distribution to calculate quantiles can improve forecasts when there are shifts

The role of air–sea coupling in the simulation of the Madden–Julian oscillation in the Hadley Centre model

The role of air–sea coupling in the simulation of the Madden–Julian oscillation (MJO) is explored using two configurations of the Hadley Centre atmospheric model (AGCM), GA3.0, which differ only in F, a parameter controlling convective entrainment and detrainment. Increasing F considerably improves deficient MJO-like variability in the Indian and Pacific Oceans, but variability in and propagation through the Maritime Continent remains weak. By coupling GA3.0 in the tropical Indo-Pacific to a boundary-layer ocean model, KPP, and employing climatological temperature corrections, well resolved air–sea interactions are simulated with limited alterations to the mean state. At default F, when GA3.0 has a poor MJO, coupling produces a stronger MJO with some eastward propagation, although both aspects remain deficient. These results agree with previous sensitivity studies using AGCMs with poor variability. At higher F, coupling does not affect MJO amplitude but enhances propagation through the Maritime Continent, resulting in an MJO that resembles observations. A sensitivity experiment with coupling in only the Indian Ocean reverses these improvements, suggesting coupling in the Maritime Continent and West Pacific is critical for propagation. We hypothesise that for AGCMs with a poor MJO, coupling provides a “crutch” to artificially augment MJO-like activity through high-frequency SST anomalies. In related experiments, we employ the KPP framework to analyse the impact of air–sea interactions in the fully coupled GA3.0, which at default F shows a similar MJO to uncoupled GA3.0. This is due to compensating effects: an improvement from coupling and a degradation from mean-state errors. Future studies on the role of coupling should carefully separate these effects.

Medieval foreign exchange: a time series analysis

This chapter applies rigorous statistical analysis to existing datasets of medieval exchange rates quoted in merchants’ letters sent from Barcelona, Bruges and Venice between 1380 and 1310, which survive in the archive of Francesco di Marco Datini of Prato. First, it tests the exchange rates for stationarity. Second, it uses regression analysis to examine the seasonality of exchange rates at the three financial centres and compares them against contemporary descriptions by the merchant Giovanni di Antonio da Uzzano. Third, it tests for structural breaks in the exchange rate series.

Protease-activated receptor 2 (PAR2) protein and transient receptor potential vanilloid 4 (TRPV4) protein coupling is required for sustained inflammatory signaling

G protein-coupled receptors of nociceptive neurons can sensitize transient receptor potential (TRP) ion channels, which amplify neurogenic inflammation and pain. Protease-activated receptor 2 (PAR(2)), a receptor for inflammatory proteases, is a major mediator of neurogenic inflammation and pain. We investigated the signaling mechanisms by which PAR(2) regulates TRPV4 and determined the importance of tyrosine phosphorylation in this process. Human TRPV4 was expressed in HEK293 cells under control of a tetracycline-inducible promoter, allowing controlled and graded channel expression. In cells lacking TRPV4, the PAR(2) agonist stimulated a transient increase in [Ca(2+)](i). TRPV4 expression led to a markedly sustained increase in [Ca(2+)](i). Removal of extracellular Ca(2+) and treatment with the TRPV4 antagonists Ruthenium Red or HC067047 prevented the sustained response. Inhibitors of phospholipase A(2) and cytochrome P450 epoxygenase attenuated the sustained response, suggesting that PAR(2) generates arachidonic acid-derived lipid mediators, such as 5',6'-EET, that activate TRPV4. Src inhibitor 1 suppressed PAR(2)-induced activation of TRPV4, indicating the importance of tyrosine phosphorylation. The TRPV4 tyrosine mutants Y110F, Y805F, and Y110F/Y805F were expressed normally at the cell surface. However, PAR(2) was unable to activate TRPV4 with the Y110F mutation. TRPV4 antagonism suppressed PAR(2) signaling to primary nociceptive neurons, and TRPV4 deletion attenuated PAR(2)-stimulated neurogenic inflammation. Thus, PAR(2) activation generates a signal that induces sustained activation of TRPV4, which requires a key tyrosine residue (TRPV4-Tyr-110). This mechanism partly mediates the proinflammatory actions of PAR(2).