Abstract
The blood oxygen level dependent (BOLD) signal evoked by brief neural stimulation, the hemodynamic response function (HRF), is a critical feature of neurovascular coupling. The HRF is directly related to local transient changes in oxygen supplied by cerebral blood flow (CBF) and oxygen demand, the cerebral metabolic rate of oxygen (CMRO2). Previous efforts to explain the HRF have relied upon the hypothesis that CBF produces a non-linear venous dilation within the cortical parenchyma. Instead, the observed dynamics correspond to prompt arterial dilation without venous volume change. This work develops an alternative biomechanical model for the BOLD response based on the hypothesis that prompt upstream dilation creates an arterial flow impulse amenable to linear description. This flow model is coupled to a continuum description of oxygen transport. Measurements using high-resolution fMRI demonstrate the efficacy of the model. The model predicts substantial spatial variations of the oxygen saturation along the length of capillaries and veins, and fits the varied gamut of measured HRFs by the combined effects of corresponding CBF and CMRO2 responses. Three interesting relationships among the hemodynamic parameters are predicted. First, there is an offset linear correlation with approximately unity slope between CBF and CMRO2 responses. Second, the HRF undershoot is strongly correlated to the corresponding CBF undershoot. Third, late-time-CMRO2 response can contribute to a slow recovery to baseline, lengthening the HRF undershoot. The model provides a powerful mathematical framework to understand the dynamics of neurovascular and neurometabolic responses that form the BOLD HRF.