Prediction of shivering heat production from core and mean skin temperatures

Abstract
Prediction formulae of shivering metabolism (M shiv) are critical to the development of models of thermoregulation for cold exposure, especially when the extrapolation of survival times is required. Many such formulae, however, have been calibrated with data that are limited in their range of core temperatures (T c), seldom involving values of less than 36°C. Certain recent studies of cold-water immersion have reported T c as low as 33.25°C. These data comprise measurements of T c (esophageal) and mean skin temperature ( s), and metabolism from 14 males [mean (SD); age = 28 (5) years; height = 1.78 (0.06) m; body mass = 77.7 (6.9) kg; body fat (BF) = 18.4 (4.5)%] during immersion in water as cold as 8°C for up to 1 h and subsequent self-rewarming via shivering under dry blanketed conditions. The data contain 3343 observations with mean (SD) T c and s of 35.92 (0.93)°C and 23.4 (8.9)°C, respectively, and have been used to re-examine the prediction of M shiv. Rates of changes of these temperatures were not used in the analysis. The best fit of the formulae, which are essentially algebraic constructs with and without setpoints, are those with a quadratic expression involving s. This is consistent with the findings of Benzinger (1969) who demonstrated that the thermosensitivity of skin is parabolic downwards with temperature peaking near a value of 20°C. Formulae that included a multiplicative interaction term between T c and s did not predict as well. The best prediction using 37°C and 33°C as the T c and T s setpoints, respectively, was found with BF as an attenuation factor: M shiv (W · m−2) = [155.5 · (37 − T c) + 47.0 · (33 − s) − 1.57 · (33 − s)2] (%BF)0.5.