Temporal and spatial MRI responses to subsecond visual activation

The temporal and spatial characteristics of oxygenation-sensitive MRI responses to very brief visual stimuli (five Hz reversing black and white checkerboard pattern versus darkness) were investigated (nine subjects) by means of serial single-shot gradient-echo echo-planar imaging (2.0 T, TR = 400 ms, mean TE = 54 ms, flip angle 30°). The use of a 0.2-s stimulus and a 90-s control phase resulted in an initial latency phase (about 2 s, no signal change), a positive MRI response (2.5% signal increase peaking at 5 s after stimulus onset), and a post-stimulus undershoot (1% signal decrease peaking at 15 s after stimulus onset) lasting for about 50–60 s. The finding that a subsecond visual stimulus elicits both a strong positive MRI response and a long-lasting undershoot provides further evidence for the neuronal origin of slow signal fluctuations seen in the absence of functional challenge and their utility for mapping functional connectivity. The additional observation that a reduction of the inter-stimulus control phase from 90 s to 9.8 s does not seem to affect the spatial extent of cortical activation in pertinent maps is of major relevance for the design and analysis of “event-related” MRI studies.     

Fast frequency selective MR imaging

With the proposed fast frequency selective MR imaging (FFSMRI) method, we focused on the elimination of all off-resonance components from the image of the observed object. To maintain imaging speed and simultaneously achieve good frequency selectivity, MRI was divided into two steps: signal acquisition and postprocessing. After the preliminary phase in which we determine imaging parameters, MRI takes place; the signal from the same object is successively acquired M times. As a result, we obtain M partial signals in k-space, from which we calculate the image of the observed object in postprocessing phase, after signal acquisition has been completed. With proper selection of parameters, it is possible to exclude from the image a majority of off-resonance components present in the observed object. However, we can decide to keep only a chosen off-resonance component in the image and eliminate all other components, including the on-resonance component and thus producing a different image from the same acquisition. The experiments with Fe(OH)(3) and oil showed that signal-to-noise ratio (SNR) can be improved by about a factor of four. The proposed FFSMRI method is suitable for frequency selective MR imaging and quantitative measurements in dynamic MRI where exclusion of off-resonance components can improve the reliability of measurement.     

Investigation of MR signal modulation due to magnetic fields from neuronal currents in the adult human optic nerve and visual cortex

Neuronal currents produce weak transient magnetic fields, and the hypothesis being investigated here is that the components of these parallel to the B0 field can potentially modulate the MR signal, thus providing a means of direct detection of nerve impulses. A theory for the phase and amplitude changes of the MR signal over time due to an external magnetic field has been developed to predict this modulation. Experimentally, a fast gradient-echo EPI sequence (TR = 158 ms, TE = 32.4 ms) was employed in an attempt to directly detect these neuronal currents in the adult human optic nerve and visual cortex using a 280-mm quadrature head coil at 1.5 T. A symmetrical intravoxel field distribution, which can be plausibly hypothesized for the axonal fields in the optic nerve and visual cortex, would result in phase cancellation within a voxel, and hence, only amplitude changes would be expected. On the other hand, an asymmetrical intravoxel field distribution would produce both phase and amplitude changes. The in vivo magnitude image data sets show a significant nerve firing detection rate of 56%, with zero detection using the phase image data sets. The percentage magnitude signal changes relative to the fully relaxed equilibrium signal fall within a predicted RMS field range of 1.2-2.1 nT in the optic nerve and 0.4-0.6 nT in the visual cortex, according to the hypothesis that the axonal fields create a symmetrical Lorentzian field distribution within the voxel.     

Cochlear reflectivity in transmission-line models and otoacoustic emission characteristic time delays

In transmission-line models of cochlear mechanics, predictions about otoacoustic-emission delays depend on the place- or wave-fixed nature of the emission generation mechanism. In this work, transient evoked otoacoustic emissions (TEOAEs), recorded at different stimulus levels in 10 young subjects, were analyzed using wavelet-based time-frequency analysis to determine the latency of each frequency component of the response. The same wave forms were Fourier analyzed to evaluate the phase-gradient delay as a function of frequency. Interpreting the relation between these two characteristic delays using cochlear models shows that most of the TEOAE response can be attributed to place-fixed reflection mechanisms. The causality principle explains observed correlations between fluctuations of the TEOAE amplitude and phase-gradient delay.     

Sound field variations caused by intense internal waves in a shallow sea with a weak thermocline

Results of an experiment on investigating the fluctuations of a sound signal with a frequency of 295 Hz on a stationary 32-km-long acoustic path in a shallow sea are presented. Hydrological conditions on the path were characterized by the presence of a weak thermocline and a frequent occurrence of intense internal wave trains. The space-time characteristics of these internal waves were measured in detail. Experimental evidence of the repetition of the forms of internal solitons in the variations of the amplitude and phase of an acoustic signal was obtained.     

Robust generation,properties and potential applications of quadratic spatial solitons generated by optical parametric amplification

Quadratic spatial solitary waves are predicted and observed experimentally near degeneracy for Type II optical parametric amplification in bulk KTP, by seeding an intense pump optical field with a control signal at half the pump wave frequency. The self-trapping of light at the two wavelengths has been shown to be insensitive to phase, polarization and magnitude of the control input, and can be exploited in novel schemes for robust all-optical processing.     

Rate dependence of human visual cortical response due to brief stimulation: an event-related fMRI study

The human brain response to a wide range of visual stimulus rates presented over a prolonged time period has been investigated by various neuroimaging techniques. However, to date, no imaging study has been performed to study the dynamic human brain response to various stimulus rates when presented in a short time. This report describes activation in the human brain due to brief visual stimulus presentation (1 s) for stimulus rates varying from 1 to 20 Hz using event-related functional MRI (fMRI). Our results show that the amplitude of the fMRI response increases with the stimulus frequency and plateaus at 6 Hz. This finding differs slightly from the results of previous blocked task paradigm experiments (with a longer time of stimulus presentation), in which the response peaks at approximately 8 Hz and then decreases. Our results are in close agreement with previously published psychophysical studies, suggesting that the fMRI signal in this experiment is indicative of cortical activity related to visual processing.     

Swept-tone transient-evoked otoacoustic emissions

Transient-evoked otoacoustic emissions (TEOAE) are responses generated within the inner ear in response to acoustic stimuli and are indicative of normal cochlear function. They are commonly acquired by averaging post-stimulus acoustic responses recorded near the eardrum in response to brief stimuli such as clicks or tone pips. In this study a new long duration stimulus consisting of a frequency swept tone is introduced for the acquisition of TEOAEs. Like stimulus frequency generated OAEs, swept-tone responses contain embedded OAEs. With swept-tone analysis, OAEs can be recovered by convolving it with a time reversed swept-tone signal resulting in time-compression. In addition, higher order nonlinear OAE responses were removed from the linear TEOAE. The results show comparable phase and time-frequency properties between the click and swept-tone evoked OAEs. Swept-tone acquisition of TEOAEs has beneficial noise properties, improving the signal to noise ratio by 6 dB compared to click evoked responses thus offering testing time savings. Additionally, swept-tone analysis removed synchronized spontaneous OAE activity from the recordings of subjects exhibiting such responses in conventional click TEOAEs. Since swept-tone stimulus consists of a single frequency component at any instantaneous moment, its analysis also provides for direct comparison with stimulus-frequency OAEs and click evoked OAEs.     

Anomalous sea surface reverberation scale model experiments

Low frequency sea surface sound backscattering from approximately 100 Hz to a few kHz observed from the 1960s broadband measurements using explosive charges to the Critical Sea Test measurements conducted in the 1990 s is substantially higher than explained by rough sea surface scattering theory. Alternative theories for explaining this difference range from scattering by bubble plumes/clouds formed by breaking waves to stochastic scattering from fluctuating bubble layers near the sea surface. In each case, theories focus on reverberation in the absence of the large-scale surface wave height fluctuations that are characteristic of a sea that produces bubble clouds and plumes. At shallow grazing angles, shadowing of bubble plumes and clouds caused by surface wave height fluctuations may induce first order changes in the backscattered signal strength. To understand the magnitude of shadowing effects under controlled and repeatable conditions, scale model experiments were performed in a 3 m x 1.5 m x 1.5 m tank at the Technical University of Denmark. The experiments used a 1 MHz transducer as the source and receiver, a computer controlled data acquisition system, a scale model target, and a surface wave generator. The scattered signal strength fluctuations observed at shallow angles are characteristic of the predicted ocean environment. These experiments demonstrate that shadowing has a first order impact on bubble plume and cloud scattering strength and emphasize the usefulness of model scale experiments for studying underwater acoustic events under controlled conditions.     

Phase vs. magnitude information in functional magnetic resonance imaging time series: toward understanding the noise

Although it has been shown that the phase of the MR signal from the brain is particularly prone to variation due to respiration, the overall physiological information contained in phase time series is not well understood. Here, we explore the different physiological processes contributing to the phase time series noise, identify their spatiotemporal characteristics and examine their relationship to BOLD-related and non-BOLD-related physiological noise in the magnitude time series. This was performed by manipulating the contribution of physiological noise to the total signal variance by modulating the TE and voxel volume, and using a short TR in order to adequately sample physiological signal fluctuations. The phase and magnitude signals were compared both before and after removal of signal fluctuations at the primary respiratory and cardiac frequencies with RETROICOR. We found that the temporal phase noise increased with TE at a faster rate than predicted by 1/TSNR as a result of physiological noise. As suggested by previous studies, the primary contributor to phase physiological noise was respiration-related effects which were manifested at a large scale (>1 cm). Notably, RETROICOR removed respiration-related large-scale artifacts and this resulted in considerable improvements in the temporal phase stability (7–90%). Physiological noise in the magnitude time series after RETROICOR consisted of low-frequency BOLD-related fluctuations (<0.13 Hz) localized to gray matter and the vasculature, and fluctuations in the vasculature correlated with slow (<0.1 Hz) variations in respiration volume and cardiac rhythm. Physiological noise in the phase signal after RETROICOR also occurred in frequencies below 0.13 Hz and was consistent with (1) residual large-scale magneto-mechanical effects correlated with slow variations in respiration volume and cardiac rhythm over time, and (2) local scale (<1 cm) effects localized in gray matter and vasculature most likely due to vascular dephasing mediated by a BOLD susceptibility change. While BOLD-related magnitude noise exhibited a TE dependence similar to BOLD, the ‘BOLD-related’ noise in the phase data increased with increasing TE and thus caused the overall phase noise to increase at a faster rate with TE than predicted by 1/TSNR. Interestingly, the spatial specificity of this effect was more evident for the higher resolution phase data, as opposed to the magnitude data, suggesting that at a higher spatial resolution the phase signal may contain more information on physiological processes than the magnitude signal.     

An empirical investigation of motion effects in eMRI of interictal epileptiform spikes

We recently developed a functional neuroimaging technique called encephalographic magnetic resonance imaging (eMRI). Our method acquires rapid single-shot gradient-echo echo-planar MRI (repetition time=47 ms); it attempts to measure an MR signal more directly linked to neuronal electromagnetic activity than existing methods. To increase the likelihood of detecting such an MR signal, we recorded concurrent MRI and scalp electroencephalography (EEG) during fast (20-200 ms), localized, high-amplitude (>50 μV on EEG) cortical discharges in a cohort of focal epilepsy patients. Seen on EEG as interictal spikes, these discharges occur in between seizures and induced easily detectable MR magnitude and phase changes concurrent with the spikes with a lag of milliseconds to tens of milliseconds. Due to the time scale of the responses, localized changes in blood flow or hemoglobin oxygenation are unlikely to cause the MR signal changes that we observed. While the precise underlying mechanisms are unclear, in this study, we empirically investigate one potentially important confounding variable — motion. Head motion in the scanner affects both EEG and MR recording. It can produce brief “spike-like” artifacts on EEG and induce large MR signal changes similar to our interictal spike-related signal changes. In order to explore the possibility that interictal spikes were associated with head motions (although such an association had never been reported), we had previously tracked head position in epilepsy patients during interictal spikes and explicitly demonstrated a lack of associated head motion. However, that study was performed outside the MR scanner, and the root-mean-square error in the head position measurement was 0.7 mm. The large inaccuracy in this measurement therefore did not definitively rule out motion as a possible signal generator. In this study, we instructed healthy subjects to make deliberate brief (<500 ms) head motions inside the MR scanner and imaged these head motions with concurrent EEG and MRI. We compared these artifactual MR and EEG data to genuine interictal spikes. While per-voxel MR and per-electrode EEG time courses for the motion case can mimic the corresponding time courses associated with a genuine interictal spike, head motion can be unambiguously differentiated from interictal spikes via scalp EEG potential maps. Motion induces widespread changes in scalp potential, whereas interictal spikes are localized and have a regional fall-off in amplitude. These findings make bulk head motion an unlikely generator of the large spike-related MR signal changes that we had observed. Further work is required to precisely identify the underlying mechanisms.     

Sources of functional magnetic resonance imaging signal fluctuations in the human brain at rest: a 7 T study

Signal fluctuations in functional magnetic resonance imaging (fMRI) can result from a number of sources that may have a neuronal, physiologic or instrumental origin. To determine the relative contribution of these sources, we recorded physiological (respiration and cardiac) signals simultaneously with fMRI in human volunteers at rest with their eyes closed. State-of-the-art technology was used including high magnetic field (7 T), a multichannel detector array and high-resolution (3 mm³) echo-planar imaging. We investigated the relative contribution of thermal noise and other sources of variance to the observed fMRI signal fluctuations both in the visual cortex and in the whole brain gray matter. The following sources of variance were evaluated separately: low-frequency drifts due to scanner instability, effects correlated with respiratory and cardiac cycles, effects due to variability in the respiratory flow rate and cardiac rate, and other sources, tentatively attributed to spontaneous neuronal activity. We found that low-frequency drifts are the most significant source of fMRI signal fluctuations (3.0% signal change in the visual cortex, TE=32 ms), followed by spontaneous neuronal activity (2.9%), thermal noise (2.1%), effects due to variability in physiological rates (respiration 0.9%, heartbeat 0.9%), and correlated with physiological cycles (0.6%). We suggest the selection and use of four lagged physiological noise regressors as an effective model to explain the variance related to fluctuations in the rates of respiration volume change and cardiac pulsation. Our results also indicate that, compared to the whole brain gray matter, the visual cortex has higher sensitivity to changes in both the rate of respiration and the spontaneous resting-state activity. Under the conditions of this study, spontaneous neuronal activity is one of the major contributors to the measured fMRI signal fluctuations, increasing almost twofold relative to earlier experiments under similar conditions at 3 T.     

Speech processing in the auditory system. I: The representation of speech sounds in the responses of the auditory nerve

In a previous paper the speech evoked spatio-temporal response patterns recorded in large populations of auditory-nerve fibers in the cat were examined [M.I. Miller and M.B. Sachs, J. Acoust. Soc. Am. 74, 502-517 (1983)]. The distribution of the relative phases of synchronized activity emerges as an important response feature reflecting the stimulus spectral parameters. Specifically, each strong low-order harmonic of the stimulus (less than or equal to 1.5-2 kHz) dominates the synchrony of a relatively broad segment of fibers near its corresponding characteristic frequency (CF) location in a pattern which mirrors the underlying traveling wave component. Each such fiber segment can be roughly subdivided into two regions: (1) a region basal to the point of resonance of the harmonic where the fiber PST histograms accumulate only small delays (or phase shifts) relative to each other reflecting the fast speed of propagation of the traveling wave, and (2) a region at or very near the point of resonance where the responses exhibit drastic relative phase shifts owing to the sudden slow down of the traveling wave and the consequent rapid accumulation of phase shifts. These rapid phase shifts thus manifest themselves as steep and localized spatial discontinuities in an otherwise relatively uniform instantaneous pattern of activity across the fiber array, all occurring at the CF locations corresponding to the low-order harmonics of the stimulus.     

Sound fluctuations caused by internal waves in a shallow sea

Estimates are presented for the fluctuations of the parameters of low-frequency sound fields in shallow-water regions of the Barents Sea, in the presence of seasonal internal gravity waves. The objective of the experiments is to reveal the main mechanisms that govern the sound fluctuations and their statistical parameters on paths of moderate lengths (50–60 to 100–120 km). Another objective is to determine the features of the sound interaction with internal waves for the sound speed profile of the summer—autumn type for which the water stratification is most pronounced. As the probing signals, continuous tonal ones produced by bottommoored sources at the frequencies about 100 and 300 Hz are used along with the 1/3-octave noise signals with the central frequency 1000 Hz, which are generated by a source deployed from a vessel. For the signal reception, both fixed bottom-moored hydrophones and a vertical chain of hydrophones are used, the chain also being deployed from the vessel. The water temperature, the salinity, and the thermocline displacements are monitored with standard hydrographic sensors. The following main results are presented: the estimate of the degree of correlation between the sound fluctuations and the parameters of the water layer, the comparison of the fluctuations in the signal amplitude envelope with the data obtained in other regions, and the estimate of the statistical parameters of the signal amplitude fluctuations, including their dependence on the path length. One more result consists in the proof of the wave nature of the interaction of sound and internal waves, which manifests itself in a strong dependence of the sound interaction with internal waves of discrete frequencies on the frequency of the probing signal and on the angle at which these wave beams intersect. An attempt is made to explain the observed phenomena by the synchronism in the interacting sound and gravity waves. The data obtained can be used to analyze and compare the fluctuations of the sound fields in the ocean, especially in shallow-water regions.     

Ku/K波段双频涡旋电磁阵列天线设计

涡旋电磁波具有携带轨道角动量的特性,利用这一特性,采用涡旋电磁波作为信号的载体,可以实现同一时间、同一频段的多路信号传输,极大地提高系统容量和频带利用率。以同轴馈电的半圆型开槽微带天线为单元,设计出了一种能工作在Ku波段和K波段的涡旋电磁阵列天线。使用三维电磁场仿真软件建模并且优化参数,最终得到在中心频率分别为17.1 GHz和19.7 GHz时,阵列天线产生的电磁波携带有轨道角动量。结论表明：该阵列天线能够产生双频涡旋电磁波。     

Investigation of the phase velocities of guided acoustic waves in soft porous layers

A new experimental method for measuring the phase velocities of guided acoustic waves in soft poroelastic or poroviscoelastic plates is proposed. The method is based on the generation of standing waves in the material and on the spatial Fourier transform of the displacement profile of the upper surface. The plate is glued on a rigid substrate so that it has a free upper surface and a nonmoving lower surface. The displacement is measured with a laser Doppler vibrometer along a line corresponding to the direction of propagation of plane surface waves. A continuous sine with varying frequencies was chosen as excitation signal to maximize the precision of the measurements. The spatial Fourier transform provides the wave numbers, and the phase velocities are obtained from the relationship between wave number and frequency. The phase velocities of several guided modes could be measured in a highly porous foam saturated by air. The modes were also studied theoretically and, from the theoretical results, the experimental results, and a fitting procedure, it was possible to determine the frequency behavior of the complex shear modulus and of the complex Poisson ratio from 200 Hz to 1.4 kHz, in a frequency range higher than the traditional methods.     

Whistler mode based explanation for the fast reconnection rate measured in the mit versatile toroidal facility

Despite the widely discussed role of whistler waves in mediating magnetic reconnection (MR), the direct connection between such waves and the MR has not been demonstrated by comparing the characteristic temporal and spatial features of the waves and the MR process. Using the whistler wave dispersion relation, we theoretically predict the experimentally measured rise time (τ(rise)) of a few microseconds for the fast rising MR rate in the Versatile Toroidal Facility at MIT. The rise time is closely given by the inverse of the frequency bandwidth of the whistler waves generated in the evolving current sheet. The wave frequencies lie much above the ion cyclotron frequency, but they are limited to less than 0.1% of the electron cyclotron frequency in the argon plasma. The maximum normalized MR rate R=0.35 measured experimentally is precisely predicted by the angular dispersion of the whistler waves.     

Phase Fluctuations of Radio Waves in the Saturated-Scintillation Regime

We consider the problem of phase fluctuations of radio waves behind a strong phase screen and in an optically thick layer. It is shown that the phase-fluctuation distribution of the received radiation at an observation point located in the saturated-scintillation area behind a turbulent phase screen is almost identical to the normal distribution of phase fluctuations of the wave on the screen. Amplitude and phase fluctuations of the received radiation are uncorrelated both for single-point and space-diversity reception if, in the latter case, the distance between the observation points exceeds the spatial scale of the diffraction component of the complex field of the received signal. Expressions for the mean square and the structural function of phase fluctuations of radio waves behind a turbulent phase screen in the saturated-scintillation regime are obtained. It is shown that the structural function of phase fluctuations in the diffraction component of the scattered field, which is exactly the function that forms saturated scintillations of the received radiation, almost coincides with the structural function of phase fluctuations on the screen. It is also shown that the diffraction effects can be neglected and the geometric-optical approximation should be used when calculating statistical characteristics of phase fluctuations of a plane wave in the saturated-scintillation area in an optically thick layer with large-scale refractive-index irregularities. __________ Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 48, No. 4, pp. 275–282, April 2005     

Magnetic resonance elastography using an air ball-actuator

The purpose of this study was to develop a new technique for a powerful compact MR elastography (MRE) actuator based on a pneumatic ball-vibrator. This is a compact actuator that generates powerful centrifugal force vibrations via high speed revolutions of an internal ball using compressed air. This equipment is easy to handle due to its simple principles and structure. Vibration frequency and centrifugal force are freely adjustable via air pressure changes (air flow volume), and replacement of the internal ball. In order to achieve MRI compatibility, all parts were constructed from non-ferromagnetic materials. Vibration amplitudes (displacements) were measured optically by a laser displacement sensor. From a bench test of displacement, even though the vibration frequency increased, the amount of displacement did not decrease. An essential step in MRE is the generation of mechanical waves within tissue via an actuator, and MRE sequences are synchronized to several phase offsets of vibration. In this system, the phase offset was detected by a four-channel optical-fiber sensor, and it was used as an MRI trigger signal. In an agarose gel phantom experiment, this actuator was used to make an MR elastogram. This study shows that the use of a ball actuator for MRE is feasible.     

Combination of BOLD-fMRI and VEP recordings for spin-echo MRI detection of primary magnetic effects caused by neuronal currents

In the present paper, for the first time, the feasibility to detect primary magnetic field changes caused by neuronal activity in vivo by spin-echo (SE) magnetic resonance imaging (MRI) is investigated. The detection of effects more directly linked to brain activity than secondary hemodynamic–metabolic changes would enable the study of brain function with improved specificity. However, the detection of neuronal currents by MRI is hampered by such accompanying hemodynamic changes. Therefore, SE image acquisition, rather than gradient-echo (GE) image acquisition, was preferred in the present work since the detection of primary neuronal and not blood oxygenation level-dependent (BOLD)-related effects may be facilitated by this approach. First of all, a precise spatiotemporal synchronization of image acquisition with the neuronal event had to be performed to avoid refocusing of the dephasing phenomenon during the course of the SE sequence. At this aim, we propose the combined use of visual evoked potential (VEP) recordings and BOLD-fMRI measurements prior to SE MRI scanning. Moreover, we exemplify by theory and experimentation how the control of artefactual signal changes due to BOLD and movement effects may be further improved by the experimental design. Finally, results from a pilot study using the proposed combination of VEP recordings and MRI techniques are reported, suggesting the feasibility of this method.