Uncertainty quantification of parenchymal tracer distribution using random diffusion and convective velocity fields

Influx and clearance of substances in the brain parenchyma occur by a combination of diffusion and convection, but the relative importance of these mechanisms is unclear. Accurate modeling of tracer distributions in the brain relies on parameters that are partially unknown and with literature values varying by several orders of magnitude. In this work, we rigorously quantified the variability of tracer distribution in the brain resulting from uncertainty in diffusion and convection model parameters.

Using the convection-diffusion-reaction equation, we simulated tracer distribution in the brain parenchyma after intrathecal injection. Several models were tested to assess the uncertainty both in type of diffusion and velocity fields and also the importance of their magnitude. Our results were compared with experimental MRI results of tracer enhancement.

In models of pure diffusion, the expected amount of tracer in the gray matter reached peak value after 15 h, while the white matter did not reach peak within 24 h with high likelihood. Models of the glymphatic system were similar qualitatively to the models of pure diffusion with respect to expected time to peak but displayed less variability. However, the expected time to peak was reduced to 11 h when an additional directionality was prescribed for the glymphatic circulation. In a model including drainage directly from the brain parenchyma, time to peak occured after 6-8 h for the gray matter.

Even when uncertainties are taken into account, we find that diffusion alone is not sufficient to explain transport of tracer deep into the white matter as seen in experimental data. A glymphatic velocity field may increase transport if a large-scale directional structure is included in the glymphatic circulation.