Uniform approach to linear and nonlinear interrelation patterns in multivariate time series

Currently, a variety of linear and nonlinear measures is in use to investigate spatiotemporal interrelation patterns of multivariate time series. Whereas the former are by definition insensitive to nonlinear effects, the latter detect both nonlinear and linear interrelation. In the present contribution we employ a uniform surrogate-based approach, which is capable of disentangling interrelations that significantly exceed random effects and interrelations that significantly exceed linear correlation. The bivariate version of the proposed framework is explored using a simple model allowing for separate tuning of coupling and nonlinearity of interrelation. To demonstrate applicability of the approach to multivariate real-world time series we investigate resting state functional magnetic resonance imaging (rsfMRI) data of two healthy subjects as well as intracranial electroencephalograms (iEEG) of two epilepsy patients with focal onset seizures. The main findings are that for our rsfMRI data interrelations can be described by linear cross-correlation. Rejection of the null hypothesis of linear iEEG interrelation occurs predominantly for epileptogenic tissue as well as during epileptic seizures.

Perception of emotion in sounded and imagined music

WE STUDIED THE EMOTIONAL RESPONSES BY MUSICIANS to familiar classical music excerpts both when the music was sounded, and when it was imagined.We used continuous response methodology to record response profiles for the dimensions of valence and arousal simultaneously and then on the single dimension of emotionality. The response profiles were compared using cross-correlation analysis, and an analysis of responses to musical feature turning points, which isolate instances of change in musical features thought to influence valence and arousal responses. We found strong similarity between the use of an emotionality arousal scale across the stimuli, regardless of condition (imagined or sounded). A majority of participants were able to create emotional response profiles while imagining the music, which were similar in timing to the response profiles created while listening to the sounded music.We conclude that similar mechanisms may be involved in the processing of emotion in music when the music is sounded and when imagined.

Transient rhythmic network activity in the somatosensory cortex evoked by distributed input in vitro

The initiation and maintenance of physiological and pathophysiological oscillatory activity depends on the synaptic interactions within neuronal networks. We studied the mechanisms underlying evoked transient network oscillation in acute slices of the adolescent rat somatosensory cortex and modeled its underpinning mechanisms. Oscillations were evoked by brief spatially distributed noisy extracellular stimulation, delivered via bipolar electrodes. Evoked transient network oscillation was detected with multi-neuron patch-clamp recordings under different pharmacological conditions. The observed oscillations are in the frequency range of 2-5 Hz and consist of 4-12 mV large, 40-150 ms wide compound synaptic events with rare overlying action potentials. This evoked transient network oscillation is only weakly expressed in the somatosensory cortex and requires increased [K+]o of 6.25 mM and decreased [Ca2+]o of 1.5 mM and [Mg2+]o of 0.5 mM. A peak in the cross-correlation among membrane potential in layers II/III, IV and V neurons reflects the underlying network-driven basis of the evoked transient network oscillation. The initiation of the evoked transient network oscillation is accompanied by an increased [K+]o and can be prevented by the K+ channel blocker quinidine. In addition, a shift of the chloride reversal potential takes place during stimulation, resulting in a depolarizing type A GABA (GABAA) receptor response. Blockade of alpha-amino-3-hydroxy-5-methyl-4-isoxazole-proprionate (AMPA), N-methyl-D-aspartate (NMDA), or GABA(A) receptors as well as gap junctions prevents evoked transient network oscillation while a reduction of AMPA or GABA(A) receptor desensitization increases its duration and amplitude. The apparent reversal potential of -27 mV of the evoked transient network oscillation, its pharmacological profile, as well as the modeling results suggest a mixed contribution of glutamatergic, excitatory GABAergic, and gap junctional conductances in initiation and maintenance of this oscillatory activity. With these properties, evoked transient network oscillation resembles epileptic afterdischarges more than any other form of physiological or pathophysiological neocortical oscillatory activity.

Theta burst transcranial magnetic stimulation is associated with increased EEG synchronization in the stimulated relative to unstimulated cerebral hemisphere

Theta burst transcranial magnetic stimulation (TBS) may induce behavioural changes that outlast the stimulation period. The neurophysiological basis of these behavioural changes are currently under investigation. Given the evidence that cortical information processing relies on transient synchronization and desynchronization of neuronal assemblies, we set out to test whether TBS is associated with changes of neuronal synchronization as assessed by surface EEG. In four healthy subjects one TBS train of 600 pulses (200 bursts, each burst consisting of 3 pulses at 30 Hz, repeated at intervals of 100 ms) was applied over the right frontal eye field and EEG synchronization was assessed in a time-resolved manner over 60 min by using a non-overlapping moving window. For each time step the linear cross-correlation matrix for six EEG channels of the right and for the six homotopic EEG channels of the left hemisphere were computed and their largest eigenvalues used to assess changes of synchronization. Synchronization was computed for broadband EEG and for the delta, theta, alpha, beta and gamma frequency bands. In all subjects EEG synchronization of the stimulated hemisphere was significantly and persistently increased relative to EEG synchronization of the unstimulated hemisphere. This effect occurred immediately after TBS for the theta, alpha, beta and gamma frequency bands and 10-20 min after TBS for broadband and delta frequency band EEG. Our results demonstrate that TBS is associated with increased neuronal synchronization of the cerebral hemisphere ipsilateral to the stimulation site relative to the unstimulated hemisphere. We speculate that enhanced synchronization interferes with cortical information processing and thus may be a neurophysiological correlate of the impaired behavioural performance detected previously.

Plenoptic camera : theory and experimental results

The purpose of this research was to develop a working physical model of the focused plenoptic camera and develop software that can process the measured image intensity, reconstruct this into a full resolution image, and to develop a depth map from its corresponding rendered image. The plenoptic camera is a specialized imaging system designed to acquire spatial, angular, and depth information in a single intensity measurement. This camera can also computationally refocus an image by adjusting the patch size used to reconstruct the image. The published methods have been vague and conflicting, so the motivation behind this research is to decipher the work that has been done in order to develop a working proof-of-concept model. This thesis outlines the theory behind the plenoptic camera operation and shows how the measured intensity from the image sensor can be turned into a full resolution rendered image with its corresponding depth map. The depth map can be created by a cross-correlation of adjacent sub-images created by the microlenslet array (MLA.) The full resolution image reconstruction can be done by taking a patch from each MLA sub-image and piecing them together like a puzzle. The patch size determines what object plane will be in-focus. This thesis also goes through a very rigorous explanation of the design constraints involved with building a plenoptic camera. Plenoptic camera data from Adobe © was used to help with the development of the algorithms written to create a rendered image and its depth map. Finally, using the algorithms developed from these tests and the knowledge for developing the plenoptic camera, a working experimental system was built, which successfully generated a rendered image and its corresponding depth map.

Parameter-free extraction of the functional network topology from intracranial EEG recordings: a time-resolved study of graph properties in focal onset seizures

Rationale: Focal onset epileptic seizures are due to abnormal interactions between distributed brain areas. By estimating the cross-correlation matrix of multi-site intra-cerebral EEG recordings (iEEG), one can quantify these interactions. To assess the topology of the underlying functional network, the binary connectivity matrix has to be derived from the cross-correlation matrix by use of a threshold. Classically, a unique threshold is used that constrains the topology [1]. Our method aims to set the threshold in a data-driven way by separating genuine from random cross-correlation. We compare our approach to the fixed threshold method and study the dynamics of the functional topology. Methods: We investigate the iEEG of patients suffering from focal onset seizures who underwent evaluation for the possibility of surgery. The equal-time cross-correlation matrices are evaluated using a sliding time window. We then compare 3 approaches assessing the corresponding binary networks. For each time window: * Our parameter-free method derives from the cross-correlation strength matrix (CCS)[2]. It aims at disentangling genuine from random correlations (due to finite length and varying frequency content of the signals). In practice, a threshold is evaluated for each pair of channels independently, in a data-driven way. * The fixed mean degree (FMD) uses a unique threshold on the whole connectivity matrix so as to ensure a user defined mean degree. * The varying mean degree (VMD) uses the mean degree of the CCS network to set a unique threshold for the entire connectivity matrix. * Finally, the connectivity (c), connectedness (given by k, the number of disconnected sub-networks), mean global and local efficiencies (Eg, El, resp.) are computed from FMD, CCS, VMD, and their corresponding random and lattice networks. Results: Compared to FMD and VMD, CCS networks present: *topologies that are different in terms of c, k, Eg and El. *from the pre-ictal to the ictal and then post-ictal period, topological features time courses that are more stable within a period, and more contrasted from one period to the next. For CCS, pre-ictal connectivity is low, increases to a high level during the seizure, then decreases at offset. k shows a ‘‘U-curve’’ underlining the synchronization of all electrodes during the seizure. Eg and El time courses fluctuate between the corresponding random and lattice networks values in a reproducible manner. Conclusions: The definition of a data-driven threshold provides new insights into the topology of the epileptic functional networks.

Campylobacter monitoring in German broiler flocks: an explorative time series analysis

Campylobacter, a major zoonotic pathogen, displays seasonality in poultry and in humans. In order to identify temporal patterns in the prevalence of thermophilic Campylobacter spp. in a voluntary monitoring programme in broiler flocks in Germany and in the reported human incidence, time series methods were used. The data originated between May 2004 and June 2007. By the use of seasonal decomposition, autocorrelation and cross-correlation functions, it could be shown that an annual seasonality is present. However, the peak month differs between sample submission, prevalence in broilers and human incidence. Strikingly, the peak in human campylobacterioses preceded the peak in broiler prevalence in Lower Saxony rather than occurring after it. Significant cross-correlations between monthly temperature and prevalence in broilers as well as between human incidence, monthly temperature, rainfall and wind-force were identified. The results highlight the necessity to quantify the transmission of Campylobacter from broiler to humans and to include climatic factors in order to gain further insight into the epidemiology of this zoonotic disease.

Climatic and human impacts on mountain vegetation at Lauenensee (Bernese Alps, Switzerland) during the last 14,000 years

Lake sediments from Lauenensee (1381 m a.s.l.), a small lake in the Bernese Alps, were analysed to reconstruct the vegetation and fire history. The chronology is based on 11 calibrated radiocarbon dates on terrestrial plant macrofossils suggesting a basal age of 14,200 cal. BP. Pollen and macrofossil data imply that treeline never reached the lake catchment during the Bølling–Allerød interstadial. Treeline north of the Alps was depressed by c. 300 altitudinal meters, if compared with southern locations. We attribute this difference to colder temperatures and to unbuffered cold air excursions from the ice masses in northern Europe. Afforestation started after the Younger Dryas at 11,600 cal. BP. Early-Holocene tree-Betula and Pinus sylvestris forests were replaced by Abies alba forests around 7500 cal. BP. Continuous high-resolution pollen and macrofossil series allow quantitative assessments of vegetation dynamics at 5900–5200 cal. BP (first expansion of Picea abies, decline of Abies alba) and 4100–2900 cal. BP (first collapse of Abies alba). The first signs of human activity became noticeable during the late Neolithic c. 5700–5200 cal. BP. Cross-correlation analysis shows that the expansion of Alnus viridis and the replacement of Abies alba by Picea abies after c. 5500 cal. BP was most likely a consequence of human disturbance. Abies alba responded very sensitively to a combination of fire and grazing disturbance. Our results imply that the current dominance of Picea abies in the upper montane and subalpine belts is a consequence of anthropogenic activities through the millennia.

Short-term feed intake is regulated by macronutrient oxidation in lactating Holstein cows

In addition to plasma metabolites and hormones participating as humoral signals in the control of feed intake, oxidative metabolic processes in peripheral organs also generate signals to terminate feeding. Although the degree of oxidation over longer periods is relatively constant, recent work suggests that the periprandial pattern of fuel oxidation is involved in regulating feeding behavior in the bovine. However, the association between periprandial oxidative metabolism and feed intake of dairy cows has not yet been studied. Therefore, the aim of this study was to elucidate possible associations existing between single feed intake events and whole-body net fat and net carbohydrate oxidation as well as their relation to plasma metabolite concentrations. To this end, 4 late-lactating cows equipped with jugular catheters were kept in respiratory chambers with continuous and simultaneous recording of gas exchange and feed intake. Animals were fed ad libitum (AL) for 24h and then feed restricted (RE) to 50% of the previous AL intake for a further 24h. Blood samples were collected hourly to analyze β-hydroxybutyrate (BHBA), glucose, nonesterified fatty acids (NEFA), insulin, and acylated ghrelin concentrations. Cross-correlation analysis revealed an offset ranging between 30 and 42 min between the maximum of a feed intake event and the lowest level of postprandial net fat oxidation (FOX(net)) and the maximum level of postprandial net carbohydrate oxidation (COX(net)), respectively. During the AL period, FOX(net) did not increase above -0.2g/min, whereas COX(net) did not decrease below 6g/min before the start of the next feed intake event. A strong inverse cross-correlation was obtained between COX(net) and plasma glucose concentration. Direct cross-correlations were observed between COXnet and insulin, between heat production and BHBA, between insulin and glucose, and between BHBA and ghrelin. We found no cross-correlation between FOX(net) and NEFA. During RE, FOX(net) increased with an exponential slope, exceeded the threshold of -0.2g/min as indicated by increasing plasma NEFA concentrations, and approached a maximum rate of 0.1g/min, whereas COX(net) decayed in an exponential manner, approaching a minimal COX(net) rate of about 2.5 g/min in all cows. Our novel findings suggest that, in late-lactating cows, postprandial increases in metabolic oxidative processes seem to signal suppression of feed intake, whereas preprandially an accelerated FOX(net) rate and a decelerated COX(net) rate initiate feed intake.

Possible earthquake trigger for 6th century mass wasting deposit at Lake Ohrid (Macedonia/Albania)

Abstract. Lake Ohrid shared by the Republics of Albania and Macedonia is formed by a tectonically active graben within the south Balkans and suggested to be the oldest lake in Europe. Several studies have shown that the lake provides a valuable record of climatic and environmental changes and a distal tephrostratigraphic record of volcanic eruptions from Italy. Fault structures identified in seismic data demonstrate that sediments have also the potential to record tectonic activity in the region. Here, we provide an example of linking seismic and sedimentological information with tectonic activity and historical documents. Historical documents indicate that a major earthquake destroyed the city of Lychnidus (today: city of Ohrid) in the early 6th century AD. Multichannel seismic profiles, parametric sediment echosounder profiles, and a 10.08m long sediment record from the western part of the lake indicate a 2m thick mass wasting deposit, which is tentatively correlated with this earthquake. The mass wasting deposit is chronologically well constrained, as it directly overlays the AD472/AD 512 tephra. Moreover, radiocarbon dates and cross correlation with other sediment sequences with similar geochemical characteristics of the Holocene indicate that the mass wasting event took place prior to the onset of the Medieval Warm Period, and is attributed it to one of the known earthquakes in the region in the early 6th century AD.

Epileptic seizures as condensed sleep: an analysis of network dynamics from electroencephalogram signals

Both deepening sleep and evolving epileptic seizures are associated with increasing slow-wave activity. Larger-scale functional networks derived from electroencephalogram indicate that in both transitions dramatic changes of communication between brain areas occur. During seizures these changes seem to be 'condensed', because they evolve more rapidly than during deepening sleep. Here we set out to assess quantitatively functional network dynamics derived from electroencephalogram signals during seizures and normal sleep. Functional networks were derived from electroencephalogram signals from wakefulness, light and deep sleep of 12 volunteers, and from pre-seizure, seizure and post-seizure time periods of 10 patients suffering from focal onset pharmaco-resistant epilepsy. Nodes of the functional network represented electrical signals recorded by single electrodes and were linked if there was non-random cross-correlation between the two corresponding electroencephalogram signals. Network dynamics were then characterized by the evolution of global efficiency, which measures ease of information transmission. Global efficiency was compared with relative delta power. Global efficiency significantly decreased both between light and deep sleep, and between pre-seizure, seizure and post-seizure time periods. The decrease of global efficiency was due to a loss of functional links. While global efficiency decreased significantly, relative delta power increased except between the time periods wakefulness and light sleep, and pre-seizure and seizure. Our results demonstrate that both epileptic seizures and deepening sleep are characterized by dramatic fragmentation of larger-scale functional networks, and further support the similarities between sleep and seizures.

Development of automated image analysis tools for verification of radiotherapy field accuracy with an electronic portal imaging device

The successful management of cancer with radiation relies on the accurate deposition of a prescribed dose to a prescribed anatomical volume within the patient. Treatment set-up errors are inevitable because the alignment of field shaping devices with the patient must be repeated daily up to eighty times during the course of a fractionated radiotherapy treatment. With the invention of electronic portal imaging devices (EPIDs), patient's portal images can be visualized daily in real-time after only a small fraction of the radiation dose has been delivered to each treatment field. However, the accuracy of human visual evaluation of low-contrast portal images has been found to be inadequate. The goal of this research is to develop automated image analysis tools to detect both treatment field shape errors and patient anatomy placement errors with an EPID. A moments method has been developed to align treatment field images to compensate for lack of repositioning precision of the image detector. A figure of merit has also been established to verify the shape and rotation of the treatment fields. Following proper alignment of treatment field boundaries, a cross-correlation method has been developed to detect shifts of the patient's anatomy relative to the treatment field boundary. Phantom studies showed that the moments method aligned the radiation fields to within 0.5mm of translation and 0.5$\sp\circ$ of rotation and that the cross-correlation method aligned anatomical structures inside the radiation field to within 1 mm of translation and 1$\sp\circ$ of rotation. A new procedure of generating and using digitally reconstructed radiographs (DRRs) at megavoltage energies as reference images was also investigated. The procedure allowed a direct comparison between a designed treatment portal and the actual patient setup positions detected by an EPID. Phantom studies confirmed the feasibility of the methodology. Both the moments method and the cross-correlation technique were implemented within an experimental radiotherapy picture archival and communication system (RT-PACS) and were used clinically to evaluate the setup variability of two groups of cancer patients treated with and without an alpha-cradle immobilization aid. The tools developed in this project have proven to be very effective and have played an important role in detecting patient alignment errors and field-shape errors in treatment fields formed by a multileaf collimator (MLC). ^

HUMAN PROPHASE CHROMOSOME UNIQUE BAND SEQUENCES; DEFINITION, CHARACTERIZATION, AND APPLICATION (CYTOGENETICS, IMAGE PROCESSING)

Extensive experience with the analysis of human prophase chromosomes and studies into the complexity of prophase GTG-banding patterns have suggested that at least some prophase chromosomal segments can be accurately identified and characterized independently of the morphology of the chromosome as a whole. In this dissertation the feasibility of identifying and analyzing specified prophase chromosome segments was thus investigated as an alternative approach to prophase chromosome analysis based on whole chromosome recognition. Through the use of prophase idiograms at the 850-band-stage (FRANCKE, 1981) and a comparison system based on the calculation of cross-correlation coefficients between idiogram profiles, we have demonstrated that it is possible to divide the 24 human prophase idiograms into a set of 94 unique band sequences. Each unique band sequence has a banding pattern that is recognizable and distinct from any other non-homologous chromosome portion.^ Using chromosomes 11p and 16 thru 22 to demonstrate unique band sequence integrity at the chromosome level, we found that prophase chromosome banding pattern variation can be compensated for and that a set of unique band sequences very similar to those at the idiogram level can be identified on actual chromosomes.^ The use of a unique band sequence approach in prophase chromosome analysis is expected to increase efficiency and sensitivity through more effective use of available banding information. The use of a unique band sequence approach to prophase chromosome analysis is discussed both at the routine level by cytogeneticists and at an image processing level with a semi-automated approach to prophase chromosome analysis. ^

A New Velocity Map for Byrd Glacier, East Antarctica, from Sequential Aster Satellite Imagery

New ice-velocity measurements are obtained for the main trunk of Byrd Glacier, East Antarctica, using recently acquired Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) imagery. The velocities are derived from the application of a cross-correlation technique to sequential images acquired in 2000 and 2001. Images were co-registered and ortho-rectified with the aid of a digital elevation model (DEM) generated from ASTER stereo imagery. This paper outlines the process of DEM generation, image co-registration and correction, and the application of the cross-correlation technique to obtain ice velocities. Comparison of the new velocity map with earlier measurements of velocity from 1978 indicates that the glacier has undergone a substantial deceleration between observations. Portions of the glacier flowing at speeds of similar to 850 m a(-1) in 1978/79 were flowing at similar to 650 m a(-1) in 2000/01. The cause of this change in ice dynamics is not known, but the observation shows that East Antarctic outlet glaciers can undergo substantial changes on relatively short timescales.

Spatial and Temporal Melt Variability at Helheim Glacier, East Greenland, and Its Effect on Ice Dynamics

Understanding the behavior of large outlet glaciers draining the Greenland Ice Sheet is critical for assessing the impact of climate change on sea level rise. The flow of marine-terminating outlet glaciers is partly governed by calving-related processes taking place at the terminus but is also influenced by the drainage of surface runoff to the bed through moulins, cracks, and other pathways. To investigate the extent of the latter effect, we develop a distributed surface-energy-balance model for Helheim Glacier, East Greenland, to calculate surface melt and thereby estimate runoff. The model is driven by data from an automatic weather station operated on the glacier during the summers of 2007 and 2008, and calibrated with independent measurements of ablation. Modeled melt varies over the deployment period by as much as 68% relative to the mean, with melt rates approximately 77% higher on the lower reaches of the glacier trunk than on the upper glacier. We compare melt variations during the summer season to estimates of surface velocity derived from global positioning system surveys. Near the front of the glacier, there is a significant correlation (on >95% levels) between variations in runoff (estimated from surface melt) and variations in velocity, with a 1 day delay in velocity relative to melt. Although the velocity changes are small compared to accelerations previously observed following some calving events, our findings suggest that the flow speed of Helheim Glacier is sensitive to changes in runoff. The response is most significant in the heavily crevassed, fast-moving region near the calving front. The delay in the peak of the cross-correlation function implies a transit time of 12-36 h for surface runoff to reach the bed.