Modeling cerebrospinal fluid dynamics across the entire intracranial space through integration of four-dimensional flow and intravoxel incoherent motion magnetic resonance imaging

Bidirectional reciprocal motion of cerebrospinal fluid (CSF) was quantified using four-dimensional (4D) flow magnetic resonance imaging (MRI) and intravoxel incoherent motion (IVIM) MRI. To estimate various CSF motions in the entire intracranial region, we attempted to integrate the flow parameters calculated using the two MRI sequences. To elucidate how CSF dynamics deteriorate in Hakim's disease, an age-dependent chronic hydrocephalus, flow parameters were estimated from the two MRI sequences to assess CSF motion in the entire intracranial region. This study included 127 healthy volunteers aged ≥ 20 years and 44 patients with Hakim's disease. On 4D flow MRI for measuring CSF motion, velocity encoding was set at 5 cm/s. For the IVIM MRI analysis, the diffusion-weighted sequence was set at six b-values (i.e., 0, 50, 100, 250, 500, and 1000 s/mm ), and the biexponential IVIM fitting method was adapted. The relationships between the fraction of incoherent perfusion (f) on IVIM MRI and 4D flow MRI parameters including velocity amplitude (VA), absolute maximum velocity, stroke volume, net flow volume, and reverse flow rate were comprehensively evaluated in seven locations in the ventricles and subarachnoid spaces. Furthermore, we developed a new parameter for fluid oscillation, the Fluid Oscillation Index (FOI), by integrating these two measurements. In addition, we investigated the relationship between the measurements and indices specific to Hakim's disease and the FOIs in the entire intracranial space. The VA on 4D flow MRI was significantly associated with the mean f-values on IVIM MRI. Therefore, we estimated VA that could not be directly measured on 4D flow MRI from the mean f-values on IVIM MRI in the intracranial CSF space, using the following formula; e + 0.25. To quantify fluid oscillation using one integrated parameter with weighting, FOI was calculated as VA × 10 + f × 0.02. In addition, the FOIs at the left foramen of Luschka had the strongest correlations with the Evans index (Pearson's correlation coefficient: 0.78). The other indices related with Hakim's disease were significantly associated with the FOIs at the cerebral aqueduct and bilateral foramina of Luschka. FOI at the cerebral aqueduct was also elevated in healthy controls aged ≥ 60 years. We estimated pulsatile CSF movements in the entire intracranial CSF space in healthy individuals and patients with Hakim's disease using FOI integrating VA from 4D flow MRI and f-values from IVIM MRI.