Comparison of compartmental analytical Blood‐Oxygen‐Level‐Dependent functional Magnetic Resonance Imaging models against Monte Carlo simulations performed over cortical micro‐angiograms

Blood oxygen level-dependent functional magnetic resonance imaging (BOLD fMRI) arises from a physiological and physical cascade of events taking place at the level of the cortical microvasculature which constitutes a medium with complex geometry. Several analytical models of the BOLD contrast have been developed, but these have not been compared directly against detailed bottom-up modeling methods. Using a 3D modeling method based on experimentally measured images of mice microvasculature and Monte Carlo simulations, we quantified the accuracy of two analytical models to predict the amplitude of the BOLD response from 1.5 to 7 T, for different echo time (TE) and for both gradient echo and spin echo acquisition protocols. We also showed that accounting for the tridimensional structure of the microvasculature results in more accurate prediction of the BOLD amplitude, even if the values for SO₂ were averaged across individual vascular compartments. A secondary finding is that modeling the venous compartment as two individual compartments results in more accurate prediction of the BOLD amplitude compared with standard homogenous venous modeling, arising from the bimodal distribution of venous SO₂ across the microvasculature in our data.