Feasibility of spiral fMRI based on an LTI gradient model

Short-term gradient imperfections in high-resolution EPI lead to Fuzzy Ripple artifacts

PURPOSE

High-resolution fMRI is a rapidly growing research field focused on capturing functional signal changes across cortical layers. However, the data acquisition is limited by low spatial frequency EPI artifacts; termed here as Fuzzy Ripples. These artifacts limit the practical applicability of acquisition protocols with higher spatial resolution, faster acquisition speed, and they challenge imaging in inferior regions of the brain.

METHODS

We characterize Fuzzy Ripple artifacts across commonly used sequences and distinguish them from conventional EPI Nyquist ghosts and off-resonance effects. To investigate their origin, we employ dual-polarity readouts.

RESULTS

Our findings indicate that Fuzzy Ripples are primarily caused by readout-specific imperfections in k-space trajectories, which can be exacerbated by short-term eddy current, and by inductive coupling between third-order shims and readout gradients. We also find that these artifacts can be mitigated through complex-valued averaging of dual-polarity EPI or by disconnecting the third-order shim coils.

CONCLUSION

The proposed mitigation strategies allow overcoming current limitations in layer-fMRI protocols: Achieving resolutions beyond 0.8 mm is feasible, and even at 3T, we achieved 0.53 mm voxel functional connectivity mapping. Sub-millimeter sampling acceleration can be increased to allow sub-second TRs and laminar whole brain protocols with up to GRAPPA 8. Sub-millimeter fMRI is achievable in lower brain areas, including the cerebellum.

Fuzzy ripple artifact in high resolution fMRI: identification, cause, and mitigation

Purpose

High resolution fMRI is a rapidly growing research field focused on capturing functional signal changes across cortical layers. However, the data acquisition is limited by low spatial frequency EPI artifacts; termed here as Fuzzy Ripples. These artifacts limit the practical applicability of acquisition protocols with higher spatial resolution, faster acquisition speed, and they challenge imaging in lower brain areas.

Methods

We characterize Fuzzy Ripple artifacts across commonly used sequences and distinguish them from conventional EPI Nyquist ghosts, off-resonance effects, and GRAPPA artifacts. To investigate their origin, we employ dual polarity readouts.

Results

Our findings indicate that Fuzzy Ripples are primarily caused by readout-specific imperfections in k-space trajectories, which can be exacerbated by inductive coupling between third-order shims and readout gradients. We also find that these artifacts can be mitigated through complex-valued averaging of dual polarity EPI or by disconnecting the third-order shim coils.

Conclusion

The proposed mitigation strategies allow overcoming current limitations in layer-fMRI protocols: (1)Achieving resolutions beyond 0.8mm is feasible, and even at 3T, we achieved 0.53mm voxel functional connectivity mapping.(2)Sub-millimeter sampling acceleration can be increased to allow sub-second TRs and laminar whole brain protocols with up to GRAPPA 8.(3)Sub-millimeter fMRI is achievable in lower brain areas, including the cerebellum.

Beat phenomena of oscillating readouts

PURPOSE

To demonstrate slowly varying, erroneous magnetic field gradients for oscillating readouts due to the mechanically resonant behavior of gradient systems.

METHODS

Projections of a static phantom were acquired using a one-dimensional (1D) EPI sequence with varying EPI frequencies ranging from 1121 to 1580 Hz on clinical 3T systems (30 mT/m, 200 T/m/s). Phase due to static B0 inhomogeneities was eliminated by a complex division of two separate scans with different polarities of the EPI readout. The temporal evolution of phase was evaluated and related to the mechanical resonances of the gradient systems derived from the gradient modulation transfer function. Additionally, the impact of temporally varying mechanical resonance effects on EPI was evaluated using an echo-planar spectroscopic imaging sequence.

RESULTS

A beat phenomenon resulting in a slowly varying phase was observed. Its temporal frequency was given by the difference between the EPI frequency and the mechanical resonance frequency of the activated gradient axis. The maximum erroneous, oscillating phase during phase encoding was ±0.5 rad for an EPI frequency of 1281 Hz. Echo-planar spectroscopic imaging images showed the resulting time-dependent stretching/compression of the FOV.

CONCLUSION

Oscillating readouts such as those used in EPI can result in low-frequency, erroneous phase contributions, which are explained by the beat phenomenon. Therefore, EPI phase-correction approaches may need to include beat effects for accurate image reconstruction.

Fast measurement of the gradient system transfer function at 7 T

PURPOSE

In this work, a new method to determine the gradient system transfer function (GSTF) with high frequency resolution and high SNR is presented, using fast and simple phantom measurements. The GSTF is an effective instrument for hardware characterization and calibration, which can be used to correct for gradient distortions, or enhance gradient fidelity.

METHODS

The thin-slice approach for phantom-based measurements of the GSTF is expanded by adding excitations that are shifted after the application of the probing gradient, to capture long-lasting field fluctuations with high SNR. A physics-informed regularization procedure is implemented to derive high-quality transfer functions from a small number of measurements. The resulting GSTFs are evaluated by means of gradient time-course estimation and pre-emphasis of a trapezoidal test gradient on a 7T scanner.

RESULTS

The GSTFs determined with the proposed method capture sharp mechanical resonances with a high level of detail. The measured trapezoidal gradient progressions are authentically reproduced by the GSTF estimations on all three axes. The GSTF-based pre-emphasis considerably improves the gradient fidelity in the plateau phase of the test gradient and almost completely eliminates lingering field oscillations.

CONCLUSION

The presented approach allows fast and simple characterization of gradient field fluctuations caused by long-living eddy current and vibration effects, which become more pronounced at ultrahigh field strengths.

Calibration of concomitant field offsets using phase contrast MRI for asymmetric gradient coils

PURPOSE

Asymmetric gradient coils introduce zeroth- and first-order concomitant field terms, in addition to higher-order terms common to both asymmetric and symmetric gradients. Salient to compensation strategies is the accurate calibration of the concomitant field spatial offset parameters for asymmetric coils. A method that allows for one-time calibration of the offset parameters is described.

THEORY AND METHODS

A modified phase contrast pulse sequence with single-sided bipolar flow encoding is proposed to calibrate the offsets for asymmetric, transverse gradient coils. By fitting the measured phase offsets to different gradient amplitudes, the spatial offsets were calculated by fitting the phase variation. This was used for calibrating real-time pre-emphasis compensation of the zeroth- and first-order concomitant fields.

RESULTS

Image quality improvement with the proposed corrections was demonstrated in phantom and healthy volunteers with non-Cartesian and Cartesian trajectory acquisitions. Concomitant field compensation using the calibrated offsets resulted in a residual phase error <3% at the highest gradient amplitude and demonstrated substantial reduction of image blur and slice position/selection artifacts.

CONCLUSIONS

The proposed implementation provides an accurate method for calibrating spatial offsets that can be used for real-time concomitant field compensation of zeroth and first-order terms, substantially reducing artifacts without retrospective correction or sequence specific waveform modifications.