How the body controls brain temperature: the temperature shielding effect of cerebral blood flow

Mingming Zhu, Joseph J. H. Ackerman, Alexander L. Sukstanskii, Dmitriy A. Yablonskiy


Normal brain functioning largely depends on maintaining brain temperature. However, the mechanisms protecting brain against a cooler environment are poorly understood. Reported herein is the first detailed measurement of the brain-temperature profile. It is found to be exponential, defined by a characteristic temperature shielding length, with cooler peripheral areas and a warmer brain core approaching body temperature. Direct cerebral blood flow (CBF) measurements with microspheres show that the characteristic temperature shielding length is inversely proportional to the square root of CBF in excellent agreement with a theoretical model. This “temperature shielding effect” quantifies the means by which CBF prevents “extracranial cold” from penetrating deep brain structures. The effect is crucial for research and clinical applications; the relationship between brain, body, and extracranial temperatures can now be quantitatively predicted.

  • brain temperature regulation
  • cerebral metabolism
  • hypothermia
  • bioheat equation

blood flow serves various roles in brain functioning. It delivers nutrients, removes waste products and, importantly, supports brain temperature regulation. Indeed, the temperature of incoming arterial blood is the main determinant of the brain temperature (9). Although this is correct for deep brain structures, it is not clear how the temperature distribution in superficial brain structures depends on the extracranial temperature. How deeply an external “cold assault” penetrates into the functioning brain remains an open question. This issue is especially important because of numerous current attempts to use mild hypothermia for treatment of stroke, multiple sclerosis, and other brain injuries (see, for example, Refs. 6, 1012, 17, 23, 24). Therapeutic results have been conflicting, in large part because quantifying the resultant cooling of human brain in vivo is beyond current technology. Therefore, developing a biophysical framework describing the determinants of heat flow, thus temperature regulation, in human brain would be a substantial achievement. Within such a framework one could predict the brain temperature response, or lack thereof, to external devices designed to impose mild hypothermia. Herein we report experimental measurements of the brain temperature profile that together with theoretical considerations quantitatively establish a general phenomenon: the temperature shielding effect of blood flow, which is responsible for brain protection against external cooling.

Major mechanisms responsible for body temperature regulation in mammals are well known (see, for example, Ref. 28). In our experiments, the body temperature was kept constant by circulating warm water as described in materials and methods. Blood flow acts as a heat exchanger with pipes (blood vessels) penetrating into all brain structures, equilibrating brain and body temperature. As a result, brain temperature has very little dependence on brain metabolism [the maximum effect is less than 1°C (32)] and is primarily defined by the temperature of incoming arterial blood. This dependence, however, is broken near the brain surface where the temperature decreases because of heat exchange with the environment. Computer simulations in adult (18) and neonatal (29) human brain suggest that temperatures substantially different from deep brain values exist only in superficial regions of several millimeters thickness in human brain. A recently proposed theoretical model (27) attributed this phenomenon to the temperature shielding effect of blood flow, which efficiently replenishes tissue heat diffusing down a temperature gradient and thus restricts the latter to a superficial region. This leads to an exponential temperature profile with a characteristic shielding length Δ, which is completely determined by cerebral blood flow (CBF) independent of mammalian species, Math(1) Here K is the tissue thermal conductivity, ρ is the tissue density, and ρb and cb are the density and specific heat of blood [in the present manuscript we use ρ = 1.0 g/cm3, ρb = 1.05 g/cm3, cb = 3.8 J·g−1·°C−1, K = 5.03·10−3 W·cm−1·°C−1, and units of ml (blood)·g (tissue)−1·min−1 for CBF]. A practical consequence, if proven experimentally, is that inducing mild hypothermia in the adult human brain via extracranial selective head cooling will decrease temperature only in a narrow superficial region, not affecting deep brain temperature.

Unfortunately, measurement technology allowing direct validation of the temperature shielding effect in humans is not yet in place. However, because of its dependence on blood flow, the temperature shielding effect, as quantitatively described by Eq. 1, is expected in other mammals. Herein we report direct measurements of brain temperature distribution and blood flow in rats and quantitatively validate Eq. 1.

In normal adult humans, gray matter blood flow is ∼0.67 ml·g−1·min−1 (25), resulting in a characteristic shielding length Δ of 3.6 mm. This shielding length is much smaller than the adult human brain size (radius ∼7 cm). In the case of small animals, like the rat, cortical blood flow can vary from 1.81 ml·g−1·min−1 in the conscious state to 0.58 ml·g−1·min−1 in the α-chloralose-anesthetized state (16). In this case, the characteristic shielding length Δ is in the range 2–4 mm, comparable to the dimensions of the rat brain. Given the specific shape of the rat head (curvature along the direction from nose to tail is much less than the curvature in the transverse plane), the temperature distribution as a function of the distance from the top of the head surface, x, and the head transverse curvature radius R (see inset in Fig. 1) can be approximated as Math(2) where I0 and I1 are modified Bessel functions, Ta is the arterial blood temperature, Tm is the body-brain temperature shift due to metabolic heat generation [under normal physiological conditions Tm in humans is ∼0.3–0.4°C (32)], and h [W·cm−2·°C−1] is an “effective” heat transfer coefficient between brain and environment. This coefficient is determined by the direct heat transfer coefficient (h0) between head surface (skin) and environment and by the presence of the intermediate layers existing between the brain and air: scalp, skull, and CSF (see details in appendix). In this model, the thickness of these insulated layers is assumed negligible; hence x can be considered as a distance from the top of the brain surface along the vertical direction (see inset in Fig. 1). With this convention for x, Eq. 2 describes brain temperature profile only from the top of the brain because side and bottom parts of the brain are insulated from the environment by thick surrounding tissue (muscles). Equation 2 can estimate profiles in these regions of the brain but requires a different convention for x. When Δ≪R and xR, Eq. 2 simplifies to an exponential function Eq. 10 in appendix. A detailed derivation of Eq. 2 is also provided in the appendix. Note that h0 depends on the possible presence of insulation between the head surface and the environment, e.g., hair, fur, clothing, etc., and on environmental conditions such as humidity, wind, etc.

Fig. 1.

Data sets showing difference between brain and body temperature (Tbrain − Tbody) vs. brain depth. At each depth, the temperature was averaged over the final 1 min (30 data points) before the probe was moved to the next depth. Symbols, experimental data; solid lines, curves representing the theoretical model, Eq. 2. Inset shows a transverse magnetic resonance image of a rat head with a cartoon that represents the cylindrical model in transverse plane with a radius R. Parameter x is the depth from the brain surface. The variation in the data confirms temperature shielding effect of blood flow as the characteristic shielding length depends on cerebral blood flow (CBF).


Animal selection.

All surgical procedures were conducted under the guidelines of Washington University Institutional Animal Care and Use Committee. Male Sprague-Dawley rats weighting 260–400 g were divided into two groups. The first group was employed for brain temperature distribution measurement (n = 10) and the second group was employed for absolute CBF measurement employing microsphere tracer techniques (n = 6).

Measurement of brain temperature distribution.

Rats in the first group (n = 10) were initially anesthetized with intraperitoneal injections of a ketamine-xylazine mixture at a dose of 72.9 mg/kg ketamine plus 10.4 mg/kg xylazine. Rat heart rate and O2 saturation values were monitored noninvasively by a veterinary pulse oximeter (model 8500V, Nonin Medical, Plymouth, MN) attached to the animal's right leg. The left femoral artery was exposed and cannulated with PE-50 tubing for mean arterial blood pressure (MABP) measurement by a pressure transducer (BLPR, WPI, Sarasota, FL). The same arterial line was also used for blood sampling. Rats were then endotracheal intubated with a 16-gauge Teflon catheter (Terumo Medical, Somerset, NJ) and mechanically ventilated (model 680; Harvard Apparatus, Holliston, MA). Breathing gas was provided with oxygen-enriched air (30% O2, balanced with nitrogen, Gateway Airgas, St. Louis, MO), which flowed by the gas-in opening of the ventilator at 1 l/min during the experimental session. To stabilize breathing gas temperature, flow-by gas was passed through a heat exchanger in equilibrium with a room temperature water bath before reaching the animal. Rats were secured in a stereotaxic headframe (David Kopf Instruments, Tujunga, CA). To maintain constant body core temperature (Tbody, considered the same as rectal in this paper), the torso of each rat was surrounded by Tygon tubing through which temperature-controlled water circulated (model DC10-B3; Thermo Haake, Karlsruhe, Germany). A T-type Teflon-insulated thermocouple wire (OMEGA, Stamford, CT) was inserted rectally to a depth of 7 cm to monitor Tbody.

Twenty minutes after the initial dose of ketamine-xylazine, 1.2% isoflurane was introduced into the breathing gas. A small incision was then made along the interaural line on the top of the rat's head. Minimal retraction of the scalp exposed a small area of skull around left side of the lambda. A 0.3-mm diameter burr hole was drilled 2–3 mm anterior to and lateral from the lambda [using orientations defined in Paxinos and Watson's The Rat Brain in Stereotaxic Coordinates (19a)]. A small amount of Vaseline was applied around the hole for thermal insulation from the environment. A 36-gauge T-type thermocouple probe (Physitemp Instruments, Clifton, NJ) was fixed to the stereotaxic headframe support and connected to a data-acquisition module (model OMB-DAQ-55; OMEGA, Stamford, CT). The stereotaxic support can move vertically to within an accuracy of 0.1 mm. Initially the tip of the thermocouple probe was placed in the space between skull and brain surface, touching superficial cerebral cortex. The incision was closed to its presurgical position, and a layer of cotton baffle was put around the surgical area for thermal insulation purposes. Nothing further was done for a period ending ∼30 min after administration of 1.2% isoflurane during which temperature and other physiological parameters stabilized. After this stabilization period, the temperature probe was moved vertically downward to the core of brain in 1-mm steps until it reached a point 9 mm deep as measured from the initial surface position. At each step, the probe was left in position at least 2 min to ensure a steady-state brain temperature measurement. Both brain and body-core temperatures were continuously monitored and recorded by using Personal DaqView data-acquisition software (IOtech, Cleveland, OH) with a 0.5 s−1 sampling rate.

Microsphere CBF measurement.

Rats in the second group (n = 6) were anesthetized and ventilated as described above. With each subject the left heart ventricle was catheterized through the right carotid artery with heparinized PE-50 tubing. A pressure transducer was connected to the arterial line to determine the moment of penetration into the left ventricle valve, i.e., when a sudden broadening of blood pressure range was observed. The left femoral artery was catheterized and used to measure arterial blood pH, Pco2 (PaCO2), and Po2 (PaO2) (i-STAT Portable Clinical Analyzer, Heska, Fort Collins, CO). The same femoral arterial line was then connected to a precision syringe pump (model 22, Harvard Apparatus) to withdraw blood at a predetermined rate as a reference for the CBF calculation. Once blood-gas parameters stabilized, a rapid 0.8-ml bolus of a well-mixed solution of microspheres (Samarium labeled, 15-μm diameter, 25 million/ml, BioPAL, Worcester, MA) was injected into the left ventricle. Immediately before the microsphere injection, the syringe pump was switched on to collect femoral arterial blood at a rate of 0.5 ml/min for 2 min. At the end of this 2-min sampling period, each rat was euthanized with intraperitoneal injection of 0.5 ml of Euthasol and its brain was carefully removed. Each brain was immediately dissected into four pieces (left unspecific cortex, left striatum, right unspecific cortex, and right stratum) and each was weighed. These brain tissue sections were then rinsed in sansSaLine sodium-free buffer solution (BioPAL). All tissue and blood samples were dried at 70°C overnight and sent to BioPAL for microsphere counting analysis.

These microsphere measurements provided regional CBF (rCBF) of each dissected brain tissue section. The theory underlying the method and its applications have been discussed elsewhere (5). In short, the absolute rCBF values (ml·g−1·min−1) were calculated based on Math(3) where Os and w are the total microsphere content and weight, respectively, of each tissue sample, Oref is the total microsphere content of the reference blood sample, and Fref is the known reference blood collection rate (0.5 ml/min).

Data collection and statistical analysis.

The difference (Tdiff) between the brain temperature (Tbrain) and Tbody was calculated at each location (depth) in the brain, using the convention Tdiff = Tbody − Tbrain. The value of Tdiff for each depth within the brain was determined by averaging data over a 1-min period (30 individual data points) after the temperature had reached steady state following the new positioning of the probe. Rectal temperature was controlled at 37.0 ± 0.5°C in all experiments. Thermocouple temperature probes were calibrated to within the accuracy of 0.05°C via a factory-certified thermometer (Fisher Scientific, Pittsburgh, PA) before and after each experiment. The temperature of the room in which these experiments were performed was maintained at 23.0 ± 0.3°C. In each experiment, 5 min before brain temperature measurements commenced or microsphere injection, ∼0.1 ml blood was retrieved from femoral artery for blood-gas analysis. Paired Student's t-test was used to analyze results from the same group. Unpaired Student's t-test was used to compare results from different groups. Statistical significance was accepted at the 95% confidence level. Unless otherwise noted, measurements are reported as means ± SD.


The brain-body temperature difference vs. depth into the brain is plotted in Fig. 1. Brain temperature profiles were well described by Eq. 2. For seven rats all the parameters (Tm, h, Δ, and R) were estimated. The derived head radius R = 0.99 ± 0.13 cm (mean ± SD) is quite close to the actual head dimension. Modeling the temperature profiles of the other three rats resulted in greater uncertainty in estimated parameters. In these three cases, R was held fixed at 0.99 cm and only three parameters (Tm, h, and Δ) were estimated. The estimated value (over all 10 subjects) of the effective heat transfer coefficient was h = (4.18 ± 0.06)·10−3 W·cm−2·°C−1 and the estimated value of Tm was 0.19 ± 0.19°C.

Modeling data via Eq. 2 allow estimation of the characteristic length Δ with corresponding CBF determined via Eq. 1. Table 1 lists Δ and calculated results for CBF together with physiological parameters. Data demonstrate relatively diverse interanimal physiological parameters as well as derived CBF values. Statistical analysis of derived CBF values vs. all physiological parameters showed good correlation only with MABP (R2 = 0.90) and PaCO2 (R2 = 0.55) but not with PaO2 (R2 = 0.00) or pH (R2 = 0.31) values. Statistical analysis also showed very low correlation between MABP and PaCO2 (R2 = 0.34), which indicated independency of the two parameters. A linear multiregression analysis of derived CBF against MABP and PaCO2 for all rats gave (R2 = 0.96) Math(4) where CBF0 is the CBF at MABP = 100 mmHg and PaCO2= 40 mmHg; A and B are correlation coefficients for MABP and PaCO2, respectively. The obtained values are CBF0 = 0.94 ± 0.19 ml·g−1·min−1, A = 0.0193 ± 0.0024 ml·g−1·min−1·mmHg−1, and B = 0.0124 ± 0.0043 ml·g−1·min−1·mmHg−1. This equation is written in a form that reflects the deviation from a normal physiological condition of MABP = 100 mmHg and PaCO2 = 40 mmHg. The quality of the correlation between derived CBF and blood flow-related physiological parameters (PaCO2 and MABP) is demonstrated in Fig. 2. Although the correlation of CBF with PaCO2 is well appreciated (25), a substantial component of CBF depends independently on MABP, likely reflective of different depths of isoflurane anesthesia. Lower levels of anesthesia generate higher MABP, which induces higher CBF (e.g., Ref. 22). No correlation of CBF with changes in PaO2 was observed, which is consistent with the known independence of CBF for PaO2 above 50 mmHg (25).

Fig. 2.

Graphical depiction of CBF correlation with both mean arterial blood pressure (MABP) and arterial Pco2 (PaCO2) as determined by linear multiregression analysis, Eq. 4. Abscissa is CBF calculated from regression expression with MABP and PaCO2 values as inputs. Ordinate is CBF obtained either from the theoretical modeling of measured brain temperature distribution (first group rats, ▪) or from microsphere measurement of CBF (second group rats, ○). Continuous lines (on both sides) show 95% prediction limit based on multiregression analysis of group 1 subjects. Dotted line is y = x. Units of CBF for both axes are ml·g−1·min−1.

View this table:
Table 1.

Measured physiological parameters and CBF derived from fitting the measured temperature profiles to the theoretical model

To validate the above described relationship between characteristic shielding length and blood flow we employed the “gold standard” microsphere CBF measurement. Results, including corresponding physiological parameters, are reported in Table 2. The CBF values from right striatum and left striatum were not different. There was a trend toward a difference between CBF from right cortex and left cortex but the difference did not reach statistical significance (P = 0.06). Microsphere data from the right cortex and right striatum were not used in measuring mean CBF because of a well-documented potential bias resulting from occlusion of the right carotid artery. The microsphere-measured mean CBF of the left cortex and left striatum was 0.70 ± 0.13 ml·g−1·min−1 (n = 6).

View this table:
Table 2.

Physiological parameters and CBF from microsphere flow measurements

To test the equivalence of microsphere-measured and temperature distribution-derived CBF, we plot microsphere-measured CBF data on the same graph (Fig. 2) as temperature distribution-derived CBF. The results are entirely consistent: all microsphere-measured CBF data are within 95% prediction limit based on multiregression analysis of temperature distribution-derived CBF data from group 1 subjects.


Despite the differences in brain dimension and in some of the anatomical structures between rat and human, the vascular networks supplying brain have great similarities. For example, both rats and humans possess a circle of Willis, a structure that allows mixing of incoming arterial blood from different arteries before distributing to the entire brain. The blood supply directions for both rat and human are also similar: arterial blood flows from bottom of brain to the top and then venous blood drains from the top (superior sagittal sinus) back to the neck. A similar vascular architecture is characteristic of most mammals and, importantly, provides a dense blood flow distribution throughout the brain. From this point of view, the temperature gradient theory introduced herein is broadly applicable. It should be pointed out, however, that some differences of vascular structures compared between the brain of animals and humans may change the thermal parameters employed in the theoretical model. One such example is a carotid rete, a compact network of intertwined arteries that lies within a venous lake. It has been found in cat, sheep, goat, ox, and pig, but not in human, rabbit, or rat brains (4). The carotid rete allows constant heat exchange between incoming arterial blood and outgoing venous blood, which leads to a lower incoming blood Ta than that of the body core temperature Tbody. However, the carotid rete is located only in the initial region of the brain's arterial supply; it should not change the surface temperature distribution and the theoretical model discussed in our study.

Our measurements confirm the biophysical model's predictions of an exponential temperature distribution near the brain surface and of a crucial dependence of the characteristic temperature shielding length on the blood flow, Eq. 1. Importantly, this agreement of measurement and theory establishes a fundamental restriction on the thickness of the layer of the brain surface that can be affected by extracranial brain cooling. Although external cold devices have been broadly applied in clinical studies to effect physiological protection via brain cooling, outcomes have been mixed. Our results speak directly to this. Figure 3 demonstrates the dependence of temperature shielding length on CBF. In rats the characteristic temperature shielding length is ∼2–4 mm, in a normal adult human it is ∼4 mm, in a term infant it is ∼6 mm, and in a premature (27 wk) infant it can extend to 12 mm.

Fig. 3.

Dependence of the brain temperature characteristic shielding length on the CBF. In general, the greater the latter, the shorter the former, hence the better the temperature shielding of the brain by the blood flow. CBF values corresponding to normal adult humans (25), term and premature (27 wk) infants (31), and nonanesthetized rats (16) are marked with arrows. Characteristic shielding length values are derived from Eq. 1 with physical parameters the same as in Table 1.

Of course, the specific temperature differential between deep brain regions and the brain surface (at x = 0 in Eqs. 2 and 10) also depends on other experimental conditions. For our experimental situation, the temperature differential between deep brain regions and the brain surface in rats is ∼2–3°C (see Fig. 1). This quantity, however, crucially depends on the effective heat transfer coefficient h describing the heat exchange between brain and the environment. If this exchange is very low (formally, h → 0), the temperature on the brain surface coincides with deep brain temperature and there is no differential. This extreme can be readily approximated by insulating the head. The opposite extreme of infinitely high heat exchange, however, cannot be achieved because the presence of the intermediate layers: CSF, skull, and scalp. As a result, the effective heat transfer coefficient h is always smaller than the direct heat transfer coefficient h0 (see appendix, Eqs. 89). The latter is determined by heat radiation, evaporation, etc. This quantity is usually unknown because its evaluation requires knowledge of numerous factors including air humidity, the type and quality of hair or fur (or absence thereof) on the head, etc. However, it can be estimated by means of Eqs. 89 if the temperature distribution in the brain is known. Using the effective heat transfer coefficient h = 4.18·10−3 W·cm−2·°C−1 obtained by fitting the model function (2) to experimental data, and the thermal resistance of the intermediate layers in rats ρ = 40 W−1·cm2·°C (see appendix), we get h0 = 6·10−3 W·cm−2·°C−1. Note that the difference between h and h0 is of the same order as h0. Obviously, for smaller values of the direct transfer coefficient (better insulation), the relative contribution of the intermediate layers to the brain's protection against cold becomes smaller.

In humans, the temperature of the brain surface and, consequently, the temperature differential between deep brain regions and the brain surface is unknown. Assuming the direct heat transfer coefficient h0 in humans has the same value h0 = 6·10−3 W·cm−2·°C−1 and, using Eqs. 89 along with the estimated value for the heat resistance of the intermediate layers for humans ρ = 170 W−1·cm2·°C, we get the effective heat transfer coefficient for human h = 3·10−3 W·cm−2·°C−1. Substituting this value of h into Eq. 11, for Ta = 37°C, Te = 20°C and the typical characteristic shielding length in humans Δ = 3.6 mm, we find the temperature differential between deep brain regions and the brain surface ∼3°C. Note, however, that such a high value of the differential is obtained by using the same value of h0 characteristic to a rat that has been subjected to a surgery (preceded by shaving). For unperturbed head surfaces, the direct heat transfer coefficient can be much smaller and, as a consequence, the temperature differential between deep brain regions and the brain surface would also be smaller. We measured a scalp temperature in five normal human volunteers seating quietly at ambient temperature of 21°C. Temperature on the bare skin surface was 32.1 ± 1.0°C and temperature on the skin surface covered by hair was 34.2 ± 0.6°C. By making use of theoretical equations derived previously (27), we can estimate that the temperature differential between deep brain and brain surface is ∼2.5°C in the area of bare skin and 0.9°C in the area of hair-covered skin. Corresponding effective heat transfer coefficients are h = 1.5·10−3 and 0.8·10−3 W·cm−2·°C−1. As expected, these numbers are smaller compared with those estimated above for surgically operated case.

The temperature distribution in the brains of small animals, like the rats employed in this study, is much more strongly influenced by the environmental temperature than will be the case with large subjects. The characteristic temperature shielding length in rats, 2–4 mm depending on the CBF, is comparable to the rat brain radius, which is ∼5 mm. Hence, the temperature distribution in the rat brain, especially under anesthesia, is substantially inhomogeneous; even deep brain regions are affected by lower environmental temperatures. Because of their larger body surface-to-body volume ratio (e.g., in rats surface-to-volume ratio is ∼10 times higher compared with adult humans), small animals loose heat more readily to the surrounding environment than do larger animals. Countering this effect, small animals have higher CBF, which helps in maintaining temperature equilibration between body and brain.

Because the effectiveness of extracranial cooling depends on the relationship between the brain size and the characteristic shielding length, projecting the hypothermic neuroprotection observed with small animals to humans is problematic. The characteristic temperature shielding length in small animals (like rats and mice) is on the order of the brain size. This allows successful extracranial brain cooling in these species. However, in humans and large animals, the characteristic temperature shielding length is much smaller than the brain size, limiting the potential cooling zone to the superficial brain regions only.

Other methods for brain cooling rely on either cooling the subject's body (see, for example, Refs. 11, 23) or directly cooling the incoming arterial blood (10, 17). In these cases, changes in the brain temperature are caused by changes in the temperature of arterial blood Ta. This approach can be successful; it causes synchronous global changes in the brain temperature compared with a superficial effect achieved through extracranial cooling.

Brain temperature is maintained by CBF, which serves as a body-coupled heat exchanger system penetrating all brain structures. Data presented herein quantitatively validate the manner in which CBF protects against extracranial cooling via the temperature shielding effect. It is often assumed that cerebral metabolism is responsible for maintaining brain temperature. However, according to Eq. 2 cerebral metabolism can only change brain temperature, compared with the temperature of incoming arterial blood, by the quantity Tm. Accordingly, for hematocrit level of 40% and arterial blood oxygen saturation of 100%, the maximum effect of metabolism on brain temperature cannot exceed 0.9°C even if all delivered oxygen is consumed [oxygen extraction fraction (OEF) = 1]. In practice, OEF rarely exceeds 40%, which results in a 0.3–0.4°C increase in brain temperature compared with that of incoming blood (32). Experimental measurements of arterial-venous temperature difference in monkeys (9) and humans (19) agree with this prediction. Small animals normally have higher metabolic rates compared with large animals. But this does not significantly affect brain temperature because more efficient heat removal is provided by higher CBF, the dominant effect in temperature regulation. Thus cerebral metabolic activity does not markedly affect brain temperature, and, indeed, it has been reported that patients in whom cerebral metabolism was reduced to 55% of normal did not show a larger brain-rectal temperature difference (15). A corollary is that body temperature is the primary determinant of deep brain temperature. Our experimental data demonstrate that deep brain temperature in rats under our experimental conditions is slightly (<1°C) lower than the deep body temperature. This can be due to already discussed effects of high heat exchange coefficient for surgically operated rats employed in this study and most importantly due to comparable characteristic shielding length and rat brain size.

Another interesting finding of our study is a linear correlation of CBF with both PaCO2 and MABP. Although the former is well known, the proportionality of CBF and MABP found in Eq. 2 has never been studied extensively. A change of blood pressure in an animal under controlled experimental conditions normally implies a loss of physiological autoregulation. In our study, a modest range of blood pressure was found for different subjects (75–100 mmHg). However, for each rat the blood pressure was well regulated and did not change markedly during the experiment. We speculate that the modest variation of intersubject CBF and MABP arises from different achieved levels of anesthesia despite administering the same concentration of isoflurane. The physiological significance of this correlation is an area for further study.

The presence of temperature gradients, reported in this paper, results in physiologically different conditions for different parts of the brain. Several observations suggest these differences might be important for brain functioning. The temperature of the human body is normally regulated in a very narrow range; sudden deviations, even by 1°C, are a sign of health abnormalities. Interventional changes in body temperature are known to substantially impact physiological processes. In the past, for example, placement of a patient in a hot tub was widely used as an “amplifier” for diagnosis of some neurological disorders (7, 26). It is well established that temperature substantially affects the affinity of hemoglobin for oxygen (blood oxygen saturation level may change by several percent per 1°C) (8), the rate of chemical reactions [the average van't Hoff temperature coefficient in the brain, Q10, is 2.3 (28)], and especially, the activation rate of heat-evoked membrane currents [Q10 up to 25.6 (30)]. It was recently reported that small temperature variations can have profound effects on brain integrative functions by affecting protein geometry, protein assembly, and alterations in protein expression (2, 3). Furthermore, dramatic changes in spine number on mature hippocampal dendrites have been found with alterations in temperature (13), and profound temperature dependence of the rate of vesicle pool depletion and recruitment have also been reported (14).

In summary, these measurements establish an exponential temperature distribution near the brain surface with the characteristic temperature shielding length crucially dependent on the CBF. This temperature shielding effect of blood flow imposes a fundamental restriction on the thickness of the layer of the brain surface that can be affected by the ambient temperature. Our results imply that brain temperature regulation is dominated by CBF in more ways than was thought previously: not only is deep brain temperature mostly dictated by temperature of incoming arterial blood but the temperature profile in the brain is also determined by CBF.


Theory of Temperature Distribution in a Small Animal Brain

Our approach is based on a static bioheat equation originally proposed by Pennes (20) for description of the temperature distribution in organs: Math(5) where the function T = T(r) describes the temperature distribution, q is the rate of metabolic heat generation [W/cm3], and all other notations are the same as in Eqs. 12. All the parameters entering the bioheat equation are assumed to be uniform across the brain. Boundary conditions at the brain's surface supplement the bioheat equation and will be specified below. If surface effects (heat exchange with environment) are ignored, the temperature distribution is homogeneous (32): Math(6) where ΔH° = 470 kJ per mole of O2 is the enthalpy of the net chemical reaction of oxygen and glucose (25), ΔHb = 28 kJ is the energy used to release oxygen from hemoglobin(1), [O2] is oxygen concentration in blood, and OEF is the oxygen extraction fraction. According to the general concept of a physiologically defined baseline state of brain function (achieved when normal subjects rest quietly but awake with eyes closed) (21), the OEF is practically uniform across the brain, despite significant differences in the blood flow and oxygen consumption in different brain regions. Hence, as described in Eq. 6, the uniformity of the OEF in the resting state leads to a uniformity of the baseline temperature (temperature in the resting state) in deep brain regions (32).

If the characteristic temperature shielding length Δ is much smaller than the brain size, the curvature of the brain can be ignored. In this case, the temperature distribution depends only on the distance x from the brain surface, T = T1(x). The boundary condition for the bioheat equation is Math(7) where h is the effective heat transfer coefficient. This quantity is determined by a direct heat transfer coefficient (h0) between skin and air as well as by the presence of the “intermediate” layers between the brain and air: CSF, skull, and scalp. A general expression relating h and h0 has a rather cumbersome structure depending on blood flow in the scalp (Fscalp) and the thicknesses (dj) and thermal conductivities (Kj) of the intermediate layers (27). However, for realistic values of Fscalp, it can be approximated by a simple expression independent of Fscalp: Math(8) where ρ is the thermal resistance of the intermediate layers: Math(9) The thermal conductivities of the intermediate layers are known from literature: (KCSF, Kskull, Kscalp) = (5.82, 11.6, 3.4)·10−3 W/(cm·°C) [cited from Table 1 in (18)]. In rats (dskulldscalp ≈ 1 mm, dCSF < 0.1 mm), the thermal resistance is ρ ≈ 40 W−1·cm2·°C; in humans dskull ≈ 5 mm, dscalp ≈ 3 mm, dCSF ≈ 2 mm (dskull ≈ 5 mm, dscalp ≈ 3 mm, dCSF ≈ 2 mm), ρ ≈ 170 W−1·cm2·°C. However, in both cases, these values of ρ are much smaller than 1/h0, hence the contribution of the intermediate layers to the effective heat transfer coefficient in Eq. 8 is small (see discussion).

The solution to Eq. 5 with the appropriate boundary condition (11) is Math(10) The temperature gradient ΔT between deep brain regions (x ≪ Δ), where temperature is equal to Ta + Tm, and the brain surface (at x = 0) is given by Math(11) As we already mentioned, in humans, blood flow in the gray matter is ∼0.67 ml·g−1·min−1 (25), resulting in the characteristic length Δ of 3.6 mm. This number is much smaller than the human brain size (radius ∼7 cm), hence Eq. 10 adequately describes the temperature distribution in human brain.

The situation is different in a rat brain where blood flow can vary from 1.81 ml·g−1·min−1 in conscious rats to 0.58 ml·g−1·min−1 in α-chloralose-anesthetized rats in the motor cortex area (16). In this case, the characteristic length Δ is in the range 2–4 mm, which is comparable to the size of the rat brain. Hence, correction for the brain curvature should be made. Given the specific shape of the rat head, a cylindrical model can be used to describe the temperature distribution. In the cylindrical model with long-axis parallel to the rat body's principal axis, the cylinder radius R represents the rat head radius in the transverse plane. The temperature distribution in this model depends on the distance from the cylinder's axis r, T = T2(r). The boundary condition at the surface r = R is Math(12) Solution of the bioheat Eq. 5 with the appropriate boundary condition (12) is given by Eq. 2, which is the theoretical model used in this paper.

To further demonstrate that the brain curvature is an important parameter in estimating blood flow in rats, we tested the adequacy of both models (Eq. 10 and Eq. 2) to fit experimental data for brain temperature profiles in rats. The model that takes into account brain curvature, Eq. 2, has consistently higher R2 values (P = 0.009) than the planar model, Eq. 10. This corroborates our proposition that correction for brain curvature should be taken into account in the theory of brain temperature distribution in rats.


This study was supported by NIH Grants RO1-NS41519 and R24-CA83060 (NCI Small Animal Imaging Resource Program).


The authors are grateful to Drs. Marcus Raichle and Sheng-Kwei Song for helpful suggestions and discussions.


  • The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked “advertisement” in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.


  1. 1.
  2. 2.
  3. 3.
  4. 4.
  5. 5.
  6. 6.
  7. 7.
  8. 8.
  9. 9.
  10. 10.
  11. 11.
  12. 12.
  13. 13.
  14. 14.
  15. 15.
  16. 16.
  17. 17.
  18. 18.
  19. 19.
  20. 19a.
  21. 20.
  22. 21.
  23. 22.
  24. 23.
  25. 24.
  26. 25.
  27. 26.
  28. 27.
  29. 28.
  30. 29.
  31. 30.
  32. 31.
  33. 32.
View Abstract