Cardiovascular Stress and Characteristics of Cold-Induced Vasodilation in Women and Men during Cold-Water Immersion: A Randomized Control Study.

BACKGROUND: Cold-induced vasodilation (CIVD) is a phenomenon that refers to a paradoxical increase in finger temperature that sometimes occurs during cold exposure. The aim of this study was to compare CIVD responses between women and men, during exposure to different environmental conditions.
METHODS: Seven men and seven women participated in a matched controlled study consisting of a familiarization protocol followed by three experimental sessions (cool (10.8 °C WBGT), thermoneutral (17.2 °C WBGT), and hot (27.2 °C WBGT)). In each session, participants were asked to immerse their left hand and foot in warm water (35 ± 1 °C) for five minutes. Thereafter, the left hand and foot were immersed in cold water (8 ± 1 °C) for 40 min. After that, the left hand and foot were removed from the water and participants remained seated for five minutes.
RESULTS: For a matched thermal stress, women experienced an elevated cardiovascular strain (heart rate and in some cases mean arterial pressure) and higher frequency of CIVD reactions (men: 31 vs. women: 60) in comparison to their male counterparts.
CONCLUSIONS: The present study demonstrated that women experienced elevated cardiovascular strain and higher frequency of CIVD reactions, particularly in the toes, compared to their male counterparts during cold-water immersion.

1. Introduction

Cold-induced vasodilation (CIVD) is a phenomenon that has been investigated since 1930 [1,2,3,4,5,6,7] and refers to a paradoxical increase in finger and toe temperature that sometimes occurs during cold exposure [1,2,8,9,10]. For this reason, it has been hypothesized that CIVD may act as a protective mechanism to reduce the risk of cold injury in the limbs [5,11,12,13]. In this light, some authors support the notion that CIVD reactions are driven by central mechanisms caused by an increase in deep body temperature [3,4,14,15] while others suggest that it may emerge as a response to peripheral stimuli in which arteriovenous anastomoses play a key role [8,16,17]. Despite the plethora of studies conducted on CIVD, the triggering mechanism of this phenomenon remains unclear, and its appearance has been described as nonsystematic and highly heterogenous [12,18]. Indeed, CIVD is reportedly affected by many factors which may lead to an unpredictable occurrence of the phenomenon. For instance, intra-individual differences such as acclimatization or adaptation to cold [19], alcohol consumption [20], diet [21], exposure to hypoxic environments [22], menstrual cycle, mental stress [23], physical fitness, and smoking [24], as well as interindividual differences such as age [25], sex [10,26,27,28], ethnicity [29], and long-term adaptation to cold have been previously included in the long list of potential factors involved in CIVD [1,8,10]. The majority of these factors have been repeatedly investigated uncovering the importance of inter- and intra-individual variability for the occurrence and magnitude of CIVD [1,8,10], but there is still limited understanding on this matter.

Differences between sexes on the CIVD response have been explored in only a handful of studies [10,26,27,28]. These studies investigated finger skin temperature but did not evaluate, toe skin temperature, blood flow in the fingers or toes, important contributing factors such as clothing, as well as potential underlying mechanisms of cutaneous vasomotion. On the whole, our knowledge on the potential impacts of sex differences on CIVD is limited and this may have important implications for workers and how they cope with exposure to cold environments. This knowledge gap is important to address because most modern workplace practices involve females in the same work tasks as men, and female workers can be found doing jobs in the cold that were previously predominantly occupied by men, such as butchering [30], fishing [31], and performing army-related tasks [32,33].

The aim of the present study was to investigate and compare CIVD and other cardiovascular responses between women and men, during exposure to different environmental conditions, while evaluating an array of physiological and perceptual parameters that may have a contributing role in the occurrence of the CIVD phenomenon.

2. Materials and Methods

2.1. Ethical Approval and Participants

Participants were requested to provide a medical history, used to exclude patients with Raynaud’s syndrome/phenomenon and those under prescription medication for hypertension or other drugs that could affect vasomotion. The minimum required sample size for investigating “repeated measures, within-between factors” was calculated using the results of a previous repeated-measures study (18) which assessed CIVD reactions in men who were in three different conditions, more specifically in thermoneutral, mildly hypothermic, and mildly hyperthermic conditions. Specifically, we calculated the effect size (f) [34] for the comparisons among fingers and among participants, for each condition, reported in the previously published study regarding minimum and maximum finger temperature (Tf) [13]. To ensure high statistical power, our calculation of the minimum required sample size was based on the lowest effect size identified. This was an effect size (f) equal to 0.5 (equating to an effect size (d) of 1.0) for the comparison of maximum Tf in the mild hypothermic (8.3 ± 1.7 °C) versus the thermoneutral (11 ± 3.2 °C) condition [13]. Sample size calculations were conducted using G*Power 3.1.9.4 [35], while setting statistical power and α error probabilities at 0.90 and 0.05, respectively. Based on this calculation, 6 participants per group (12 participants in total) would provide enough power to detect differences between men and women in our study. Accordingly, seven male and seven female subjects matched for age, body mass index, and body surface area (i.e., p > 0.05 between groups) participated in the study. The anthropometric characteristics of all 14 participants are shown in Table 1.

Table 1

Anthropometric characteristics (mean ± SD) of the participants in the two groups.

Group	n	Age(Years)	Weight(kg)	Height(m)	BMI(kg/m2)	BSA(m2)
Men	7	34.0 ± 12.2	81.1 ± 17.7	1.76 ± 0.5	26.0 ± 5.0	2.0 ± 0.2
Women	7	32.1 ± 8.5	75.5 ± 15.9	1.72 ± 0.5	25.4 ± 5.1	1.9 ± 1.8

Note: BMI: body mass index; BSA: body surface area.

The experimental protocol (ClinicalTrials.gov ID: NCT04215939) conformed to the standards set by the Declaration of Helsinki and was approved by the Bioethics Review Board of the University of Thessaly Department of Physical Education and Sport Science (protocol no.: 1320). Four out of fourteen participants (two men and two women) were smokers and were asked to refrain from smoking at least 10 h prior to each experimental session. The menstrual cycle was not controlled since there is recent evidence demonstrating that women’s heat dissipation [36] and finger blood flow [37] do not differ between menstrual phases. All volunteers were given a full explanation of all the procedures and signed a written informed consent prior to participating in the study.

2.2. Experimental Protocol

The study included a familiarization protocol followed by three experimental sessions. During the familiarization session, participants were informed about all data collection procedures/equipment and underwent anthropometric and body composition (dual-energy X-ray absorptiometry) assessments. Each of the three experimental sessions was performed inside a 32.5 m3 environmental chamber (HGX22G/190-4, Bock GmbH, Frickenhausen,,Germany) under different conditions: cool (wet-bulb globe temperature: 10.8 °C; air temperature: 16 ± 1°C; relative humidity: 45 ± 5%; air velocity: 0.2 m/s; solar radiation: 0 W/m2), thermoneutral (wet-bulb globe temperature: 17.2 °C; air temperature: 23 ± 1 °C; relative humidity: 45 ± 5%; air velocity: 0.2 m/s; solar radiation: 0 W/m2), and hot (wet-bulb globe temperature: 27.2 °C; air temperature: 34 ± 1 °C; 45 ± 5% relative humidity; air velocity: 0.2 m/s; solar radiation: 0 W/m2) environments. The overall thermal stress experienced by the participants in our study was expressed by means of wet-bulb globe temperature which was found to be the most efficacious thermal-stress indicator for assessing the physiological strain [38,39,40]. The present study is a randomized control trial. The experimental sessions were administered in a random order, based on a random allocation algorithm implemented in Excel Spreadsheets (Microsoft Office, Microsoft, Washington, DC, USA).

For each session, participants arrived at the lab at the same time of the day. They had been asked to refrain from caffeine for at least two hours, from food for at least three hours, and from alcohol and exercise for at least 12 h. Upon arrival, participants wore standard clothing consisting of a long-sleeve (100% cotton) shirt and a pair of pants (100% cotton) to avoid the influence of different fabric properties on physiological data in each experimental session (Figure 1). The room temperature in the lab where the participants changed clothes before entering the environmental chamber was thermoneutral (wet-bulb globe temperature: 16.3 °C; air temperature: 22 °C; relative humidity: 45%; air velocity: 0.2 m/s; solar radiation: 0 W/m2). The estimated clothing insulation was 0.67 clo for men (underwear: 0.04; long-sleeve shirt: 0.29; pants: 0.34) and 0.69 clo for women (underwear: 0.04; bra: 0.02; long-sleeve shirt: 0.29; pants: 0.34) [41,42]. Thereafter, participants entered the environmental chamber (with the trial order [temperature] counterbalance, but randomly allocated to each participant) and assumed a seated position for 10–15 min to complete instrumentation. Once all sensors were placed, participants were requested to relax seated for 20 min with their hands supported at the level of the heart (Figure 1). After the baseline period, participants were asked to simultaneously immerse their left hand (up to the ulnar styloid process) and left foot (up to the lateral malleolus) in warm water (35 ± 1 °C) for five minutes to ensure similar starting local tissue temperature [2]. Thereafter, the left hand and foot were immersed in cold water (8 ± 1 °C) for 40 min (Figure 1 and Figure 2). After the end of the cold-water immersion period, the left hand and foot were removed from the water and participants remained seated for an additional five minutes to monitor the recovery phase. The total duration of the experiment was 70 min: 20 min baseline period, five minutes warm-water immersion, 40 min cold-water immersion, and five minutes recovery.

Figure 1

A male participant during the data collection. Note: the right sleeve of the shirt was cut at elbow height to avoid constricting blood flow. The blood pressure cuff was comfortably placed at the upper arm without constricting blood flow except when assessing arterial blood pressure. Photograph used with permission from the participant featured.

Figure 2

Representative data for two cold-induced vasodilation (CIVD) reactions (red color) in the index finger of a male participant. Onset time (green color) corresponds to the time needed to reach minimum temperature (Tmin) prior to the onset of CIVD. Amplitude represents the difference between Tmin and the highest temperature during the CIVD (Tmax). Peak time corresponds to the time interval between Tmin and Tmax.

The water tanks used for the warm immersion were 0.58 m3 for the hand and 0.85 m3 for the foot. The water tanks used for the cold immersion were 0.91 m3. Water temperature was controlled within 1 °C using an air-to-water heat pump (model KF120-B, Dongguan, China) and the water was continuously circulating to ensure uniform exposure. Throughout the protocol, both hands were positioned at the level of the heart. During the water immersion periods, the left foot was placed in a water tank on the floor, and the right foot rested on the floor (Figure 1). Water-permeable mesh fabrics were used to line the bottom of the water tanks and the floor to ensure that the immersed hand and foot made no contact with the bottom of the water tanks, and that the non-immersed foot made no contact with the floor of the environmental chamber (Figure 1).

2.3. Measurements

Gastrointestinal (Tgi) and weighted mean skin (Tsk) temperatures (both assessed in 1 min intervals), as well as finger temperature (Tf), toe temperature (Tt), skin blood flow (SkBF), sweat rate, and heart rate (HR) were continuously monitored. Raw data from these variables—assessed in 0.05 s (for SkBF), 1 s (for HR), 8 s (for Tf and Tt), and 10 s (for sweat rate) intervals—were used to provide 1 min averages. In addition, 20 s averages were calculated for Tf, Tt, and SkBF and were used to identify CIVD reactions. In terms of the remaining measurements, arterial blood pressure was assessed via automatic auscultation at the non-immersed hand every 10 min during baseline, at the end of the warm-water immersion, every 10 min during cold-water immersion, and at the end of the recovery period. Thermal comfort and thermal sensation were recorded every 10 min during baseline and every five minutes during cold-water immersion. Pain sensation was assessed at the start of cold-water immersion and then every five minutes. Tactile sensitivity of the immersed limbs was tested at the end of the baseline period, at the end of the warm-water immersion, as well as at minutes 5 and 40 of the cold-water immersion (for the latter two, the hand and foot were removed from the water for ~30 s).

The Tgi was monitored using telemetric capsules (BodyCap, Caen, France) as an indicator of central body core temperature. Although the time following ingestion does not significantly influence the validity of telemetry capsule measurements of core temperature in the absence of fluid consumption [43], the capsules were always ingested at the same time before the start of the protocol for each participant. Skin temperature from four sites (arm, chest, leg, and thigh) was recorded using wireless thermistors (iButtons type DS1921H, Maxim/Dallas Semiconductor Corp., Sunnyvale, CA, USA) and was used to calculate Tsk according to Ramanathan [Tsk = 0.3 × (chest + arm) + 0.2 × (thigh + leg)] [44].

Using surgical tape (3M Transpore Tape, 3M Canada, Inc., London, ON, Canada), six ceramic chip skin thermistors (MA-100, Thermometrics, Northridge, CA, USA) were attached on the lower part of the pad on the 2nd finger (i.e., index) of both hands and on the lower part of the pad on the 1st, 3rd, and 5th toe of the immersed foot, as well as on the 1st toe of the non-immersed foot. Data were recorded using a data logger (Smartreader 8 Plus, ACR, Vancouver, BC, Canada).

The SkBF was monitored with a laser Doppler flowmeter (PF4000 LDPM, Perimed, Stockholm, Sweden and PF5010 LDPM, Perimed, Stockholm, Sweden) at the pad of the 2nd finger in each hand as well as at the distal edge of the 1st toe of each foot. The probe (PR 407 small straight probe, Perimed) in the non-immersed 2nd finger was held in place with a plastic mini holder (diameter: 5 mm; PH 07-5, Perimed), which was fixed to the skin using double-sided adhesive strips (PF 105-3, Perimed) without constricting the finger. All other probes (413 Integrating Probe, Perimed, Stockholm, Sweden) were held in place with a plastic holder (PH 13, Perimed, Stockholm, Sweden). The SkBF data were expressed in perfusion units (PU).

Sweat rate was measured at the forehead and at the belly of the gastrocnemius muscle using a ventilated capsule system that utilizes a digital mass flow meter sensor (SFM4100, Sensirion, Staefa, Switzerland). Heart rate was monitored using a wireless heart rate system (Polar Team2, Polar Electro Oy, Kempele, Finland). Arterial blood pressure was assessed at the non-immersed hand via automatic auscultation (Omron Healthcare, M6 comfort, Kyoto, Japan). Mean arterial pressure was calculated as follows [45]: where, MAP is the mean arterial blood pressure (mmHg), DBP is the diastolic blood pressure (mmHg) and SBP is the systolic blood pressure (mmHg).

In addition, thermal comfort (1 = comfortable to 5 = extremely uncomfortable) and thermal sensation (−3 = cold to +3 = hot) were recorded using standard scales [46]. Pain intensity (1 = no pain to 10 = worst pain imaginable), as well as pain distress (in Likert (1 = no pain to 10 = unbearable pain) and visual analog scale (left edge = no pain; right edge = unbearable pain)) were recorded using standardized scales [2]. Tactile sensitivity of the tip of the middle finger and the second toe of the immersed limbs was assessed using Semmes–Weinstein monofilaments and a digital esthesiometer, based on previous methods [2,47]. For both measurements, lower values indicate increased tactile sensitivity.

2.4. Statistical Analysis

The 20 s average values for Tf and Tt were used to identify CIVD reactions based on previous criteria (Figure 2) [2,8] as follows:

Number of waves (n): the minimum temperature increase to define a CIVD response was set at 1 °C.

Minimum temperature (Tmin): the lowest temperature just before the onset of CIVD in °C.

Maximum temperature (Tmax): the highest temperature during the CIVD in °C.

Onset time (Tonset): the time from immersion to Tmin in minutes.

Peak time (Tpeak): the time to reach the maximum temperature (time at Tmax minus time at Tmin) in minutes.

Average temperature (Tavg): the average temperature during the cold-water immersion period without the first five minutes of the cold-water immersion in °C.

Temperature amplitude (ΔT): the difference between Tmin and Tmax in °C.

Cardiovascular strain (systolic, diastolic, and mean arterial blood pressure, as well as HR) during the first minute of cold-water immersion was calculated by means of delta (Δ) values expressed as the difference between the last minute of the warm-water immersion and the 1st minute of the cold-water immersion.

Chi-square tests were used to compare the frequency of CIVD reactions across (1) the three different environments, (2) the fingers, (3) the toes, as well as (4) the limbs (i.e., the immersed hand vs. the foot). Paired sample t-tests along with effect sizes (d) were used to examine potential differences in the physiological strain (Tf, Tt, SkBF, Tsk, Tgi, arm/chest/leg/thigh skin temperature, perceptual scales (thermal comfort, thermal sensation, pain), systolic blood pressure, diastolic blood pressure, and HR) experienced between women and men. The comparisons were performed separately for each phase of the protocol (baseline, warm immersion, cold immersion, and recovery). Additional paired-sample t-tests with effect sizes (d) were conducted to examine potential differences in the cardiovascular strain during the first minute of cold-water immersion between women and men. The level of statistical significance in the t-tests was adjusted for multiple comparisons using the Bonferroni correction, resulting in an adjusted level alpha value of 0.008. The magnitude of Cohen’s (d) effect sizes was determined as follows: d (0.01) = very small; d (0.2) = small; d (0.5) = medium; d (0.8) = large; d (1.2) = very large; and d (2.0) = huge [48]. Effect sizes (d) were computed with Excel Spreadsheets (Microsoft Office, Microsoft, Washington, DC, USA) and all other statistical analyses were conducted with SPSS v27.0 (IBM, Armonk, NY, USA).

3. Results

Under identical experimental conditions that led to a similar heat strain between women and men (i.e., no clinically meaningful differences in Tsk, Tgi, or sweat rate; see subsection “other physiological responses”), our study demonstrated that the studied women experienced increased cardiovascular strain and more frequent CIVD reactions, especially during exposure to the cool environment. Detailed information is presented in the following subsections.

3.1. Frequency of CIVD

The total number of CIVD reactions in fingers and toes was 91 (men: 31, women: 60) across all three environments: 14 in the cool, 22 in the thermoneutral, and 55 in the hot environment (p < 0.001; Figure 3). Men showed no CIVD in the foot during the cool and neutral exposures. Chi-square analyses revealed that women had more CIVD reactions in the foot in comparison to their male counterparts (p < 0.001). No statistically significant differences between sexes in CIVD occurrence were detected in the hand or in both limbs together across all environments (Figure 3).

Figure 3

Differences (mean ± SD) between men (M) and women (F) in the occurrence of cold-induced vasodilation reactions (CIVD; top bars) and cardiovascular strain during the first minute of cold-water immersion (bottom bars) across cool (blue bars), thermoneutral (yellow bars), and hot (orange bars) environments. Cardiovascular strain during the first minute of cold-water immersion was expressed as the delta difference from the last minute of warm-water immersion and the first minute of cold-water immersion.

For CIVD reactions detected in the toes, two participants experienced CIVDs in the cool trial (men: zero, women: two; Table 2), one in the thermoneutral (men: zero, women: one; Table 3), and ten in the hot (men: four, women: six; Table 4) environment. When considering all participants simultaneously, chi-square analysis showed that a significantly higher number of participants demonstrated CIVD reactions in the toes during exposure in the hot environment (n = 10) compared to the thermoneutral environment (n = 1; p < 0.001) and the cool environment (n = 2; p < 0.001). Overall, we found significant differences in the number of participants demonstrating CIVDs in the toes across groups and environments (i.e., participants with CIVD in the toes × group × environment; p < 0.001), with CIVDs being more prevalent for women in the heat compared to the thermoneutral and cool environments. Men demonstrated no CIVD reactions in the toes in cool and thermoneutral environments.

Table 2

Heat map presenting the number of participants that demonstrated a CIVD reaction for each minute of the cold immersion during exposure in the cool environment (16 ± 1 °C). Colors indicate frequencies. Gray color indicates lack of CIVD.

		Time (min)
		0	1	2	3	4	5	6	7	8	9	10	11	12	13	14	15	16	17	18	19	20	21	22	23	24	25	26	27	28	29	30	31	32	33	34	35	36	37	38	39	40
Women	Second finger											1	1	2	2	2	3	4	4	4	4	4	4	4	4	3	3	3	3	3	2	2	1	1	1	1	1	1	1	1	1
	First toe														1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1
	Third toe																					1	1	1	1	1
	Fifth toe																				1	1	1	1	1	1
	First toe *
Men	Second finger									2	2	2	2	3	4	4	4	4	4	5	5	5	5	5	5	5	3	3	3	3	3	3	3	2	2	2	2	2	2	2	2	2
	First toe
	Third toe
	Fifth toe
	First toe *

Notes: * = non-immersed foot (all other fingers and toes refer to the immersed limbs).

Table 3

Heat map presenting the number of participants that demonstrated a CIVD reaction for each minute of the cold immersion during exposure in the thermoneutral environment (23 ± 1 °C). Colors indicate frequencies. Gray color indicates lack of CIVD.

		Time (min)
		0	1	2	3	4	5	6	7	8	9	10	11	12	13	14	15	16	17	18	19	20	21	22	23	24	25	26	27	28	29	30	31	32	33	34	35	36	37	38	39	40
Women	2nd finger									1	1	2	4	5	5	5	5	5	5	4	3	1	1	2	2	2	2	2	2	4	4	4	3	3	3	3	2	2	2	1	1	1
	1st toe										1	1	1	1	1	1	1	1	1	1																			1	1	1	1
	3rd toe											1	1	1	1	1																		1	1	1	1	1	1	1	1
	5th toe											1	1	1	1	1	1	1	1	1	1				1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1
	1st toe *
Men	2nd finger							1	1	1	2	2	3	3	3	3	3	3	3	3	3	3	3	3	3	3	3	3	2	2	3	2	2	2	3	3	3	3	3	2
	1st toe
	3rd toe
	5th toe
	1st toe *

Notes: * = non-immersed foot (all other fingers and toes refer to the immersed limbs).

Table 4

Heat map presenting the number of participants that demonstrated a CIVD reaction for each minute of the cold immersion during exposure in the hot environment (34 ± 1 °C). Colors indicate frequencies. Gray color indicates lack of CIVD.

		Time (min)
		0	1	2	3	4	5	6	7	8	9	10	11	12	13	14	15	16	17	18	19	20	21	22	23	24	25	26	27	28	29	30	31	32	33	34	35	36	37	38	39	40
Women	2nd finger								1	1	1	3	4	4	5	5	5	5	5	5	5	5	5	5	5	5	5	5	5	5	5	5	5	5	3	3	3	3	3	2	2	2
	1st toe									1	2	2	2	2	2	3	4	4	4	4	3	3	4	4	4	3	4	2	3	3	3	3	4	4	4	4	4	3	3	3	2	2
	3rd toe							1	1	3	3	3	4	4	5	5	5	5	4	4	4	3	2	3	2	3	3	3	3	4	4	4	4	4	3	3	3	2	1	1	1	1
	5th toe									2	3	3	4	5	5	5	5	5	5	4	4	4	3	2	2	2	2	2	2	2	2	2	3	2	2	2	2	1	1	1	1
	1st toe *						1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1			1	1	1	1	1	1
Men	2nd finger							1	1	1	1	1	1	2	2	3	3	3	3	2	2	3	3	3	3	3	3	3	4	4	4	4	4	4	4	4	4	3	2	2	2	2
	1st toe													2	2	2	2	2	2	2	2	2	2	2	1	2	2	2	2	2	2	2	1	1	1	1	1	1	1	1
	3rd toe													1	2	3	3	3	3	2	2	2	3	3	2	2	2	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1
	5th toe												1	2	2	2	2	2	2	2	2	2	2	3	1	1	1	1	1	2	2	2	2	2	1	1	1	1	1
	1st toe *				1	2	2	2	2	2	1	1	1	1	1	1	1	1	1	1	2	1	1	1	1	1	1	1	1	1	1	1	1

Notes: * = non-immersed foot (all other fingers and toes refer to the immersed limbs).

For the CIVDs detected in the fingers, four women revealed CIVD reactions in the cool environment, five in the thermoneutral environment, and five in the hot environment, revealing no statistically significant differences across the three environments (p = 0.458). In the male group, five participants revealed CIVDs in the fingers in the cool environment, three in the thermoneutral environment, and four in the heat, revealing no statistically significant differences across the three environments (p = 0.710). When considering all participants simultaneously, chi-square analysis showed that there were no differences in the number of participants who demonstrated CIVD reactions in the finger during exposure in the hot environment (n = 9) in comparison to the thermoneutral (n = 8; p < 0.001) and cool (n = 9; p < 0.001) environments. Overall, we found no differences in the number of participants demonstrating CIVDs in the fingers between sexes and across environments (p = 0.757).

3.2. Characteristics of CIVD

The characteristics of CIVD waves across groups and environments are shown in Appendix A (Table A1) as well as in Figure 4, and Appendix A Figure A1 and Figure A2. For the toes, as the male participants showed no CIVD in the cool or the thermoneutral environments, we were only able to compare the CIVD characteristics between groups for the CIVD waves observed in the hot environment. Thus, for the hot environment we found that the first toe demonstrated small effect sizes for the minimum (d = 0.35) and maximum (d = 0.39) temperatures as well as for the onset time (d = 0.39) between men and women (Table A1) while there were large effect sizes for the peak time (d = 0.94), the average temperature (d = 0.87), and the amplitude (d = 0.85) with women reaching faster peak time, higher average temperature, and lower amplitude. The third toe demonstrated medium effect sizes for the Tmin (d = 0.90) and very small for Tmax (d = 0.09), while there were no substantial effect sizes in the Tonset and Tpeak (very small to small effect sizes; d = 0.26 to 0.10) between the two groups. Women demonstrated lower Tavg (large effect size; d = 0.81) and higher amplitude (medium effect size; d = 0.46) compared to men. Finally, for the fifth toe, Tonset was earlier for women compared to men, demonstrating a large effect size (d = 0.84). At the same time, women demonstrated higher (by ~2 °C) toe Tavg compared to their male counterparts (medium effect size; d = 0.66).

Table A1

Characteristics of CIVD waves across groups and environments.

			n	Minimum (°C)	Maximum (°C)	Onset Time (mm:ss)	Peak Time (mm:ss)	Average (°C)	Amplitude (°C)
1st toe immersed leg	Wave 1
	Cool	men	0
		women	2	8.7 ± 0.1	10.9 ± 0.1	17:23 ± 13:45	14:00 ± 04:15	10.0 ± 0.6	2.2 ± 0.0
	Neutral	men	0
		women	1	9.7	12.1	09:23	07:00	10.9	2.4
	Hot	men	3	10.3 ± 1.4	14.3 ± 3.6	16:10 ± 07:08	09:00 ± 00:20	11.3 ± 1.3	4.0 ± 2.6
		women	5	10.8 ± 1.7	13.2 ± 2.1	13:43 ± 05:41	06:16 ± 03:33	12.4 ± 1.3	2.4 ± 1.4
	Wave 2
	Cool	men	0
		women	0
	Neutral	men	0
		women	1	10.7	12.0	27:44	06:00	10.9	1.3
	Hot	men	0
		women	3	11.1 ± 0.7	14.4 ± 2.1	20:56 ± 03:27	05:00 ± 03:06	13.4 ± 0.7	3.3 ± 2.5
3rd toe immersed leg	Wave 1
	Cool	men	0
		women	1	9.4	10.9	19:41	01:40	10.0	1.4
	Neutral	men	0
		women	1	10.0	14.2	09:43	06:40	12.0	4.2
	Hot	men	3	10.3 ± 0.9	12.9 ± 1.0	12:50 ± 01:02	05:00 ± 02:36	11.7 ± 1.1	2.6 ± 1.8
		women	6	9.2 ± 1.4	12.6 ± 3.2	11:31 ± 05:55	04:43 ± 02:57	10.3 ± 1.9	3.4 ± 1.9
	Wave 2
	Cool	men	0
		women	0
	Neutral	men	0
		women	1	11.4	14.7	31:43	06:00	12.0	3.3
	Hot	men	1	8.5	20.0	21:24	09:40	12.9	11.6
		women	3	10.3 ± 1.1	13.3 ± 1.6	24:57 ± 02:51	06:00 ± 01:46	11.9 ± 1.0	3.0 ± 1.2
	Wave 3
	Cool	men	0
		women	0
	Neutral	men	0
		women	0
	Hot	men	0
		women	1	9.7	10.8	35:07	02:00	11.2	1.0
5th toe immersed leg	Wave 1
	Cool	men	0
		women	1	10.0	11.7	19:01	03:40	11.2	1.7
	Neutral	men	0
		women	1	9.0	14.4	09:43	05:40	11.6	5.4
	Hot	men	4	10.6 ± 2.5	14.6 ± 4.3	17:08 ± 06:32	06:00 ± 02:42	11.4 ± 1.8	3.9 ± 2.7
		women	6	11.7 ± 4.8	14.9 ± 5.4	11:46 ± 06:20	05:48 ± 04:27	13.3 ± 3.5	3.2 ± 1.5
	Wave 2
	Cool	men	0
		women	0
	Neutral	men	0
		women	1	9.0	12.5	22:43	04:40	11.3	3.5
	Hot	men	2	12.1 ± 0.1	16.2 ± 4.5	30:23 ± 06:35	05:10 ± 04:00	12.7 ± 1.8	4.2 ± 4.3
		women	3	12.8 ± 3.9	15.0 ± 5.5	21:30 ± 04:49	02:20 ± 00:00	14.2 ± 5.1	2.2 ± 1.7
	Wave 3
	Cool	men	0
		women	0
	Neutral	men	0
		women	1	8.6	15.5	30:43	07:00	11.6	6.9
	Hot	men	0
		women	1	14.0	15.8	31:21	02:00	13.5	1.9
2nd finger immersed hand	Wave 1
	Cool	men	5	9.5 ± 0.9	12.7 ± 2.3	11:10 ± 06:00	08:07 ± 10:40	10.6 ± 1.2	3.2 ± 1.4
		women	4	8.4 ± 0.9	11.6 ± 1.7	13:29 ± 02:36	09:10 ± 04:08	10.1 ± 1.5	3.2 ± 1.6
	Neutral	men	3	10.0 ± 0.5	12.5 ± 1.3	10:57 ± 02:33	10:40 ± 01:09	12.4 ± 1.8	5.2 ± 1.1
		women	5	9.1 ± 1.2	12.5 ± 2.1	09:57 ± 01:38	07:36 ± 00:59	10.8 ± 1.6	3.4 ± 1.4
	Hot	men	4	9.0 ± 0.8	15.0 ± 2.6	10:53 ± 03:25	08:50 ± 05:12	12.5 ± 1.8	6.0 ± 3.2
		women	5	11.1 ± 1.4	18.7 ± 1.9	08:50 ± 01:49	08:33 ± 04:54	15.1 ± 2.3	7.6 ± 1.7
	Wave 2
	Cool	men	1	10.2	12.9	15:10	01:40	11.27	2.7
		women	0
	Neutral	men	3	10.9 ± 2.2	14.7 ± 2.2	30:37 ± 04:24	07:20 ± 04:11	12.4 ± 1.8	3.8 ± 1.2
		women	4	9.9 ± 1.9	11.9 ± 1.8	25:25 ± 03:19	05:05 ± 03:06	10.4 ± 1.6	2.0 ± 0.9
	Hot	men	2	10.4 ± 1.4	13.3 ± 1.7	24:04 ± 03:46	04:50 ± 01:11	11.5 ± 1.3	2.9 ± 0.3
		women	3	11.1 ± 0.7	14.4 ± 2.1	20:56 ± 03:27	05:40 ± 00:20	13.4 ± 0.7	3.3 ± 2.5

Notes: Minimum = minimum value prior the CIVD. Maximum = maximum value during the CIVD. Onset time = time of the minimum value − start time of the cold immersion. Peak time = time of the minimum value − time of the maximum value. Average = mean values of the last 35 min of the cold immersion. Amplitude = maximum − minimum temperature.

Figure 4

Skin temperature and blood flow (mean ± SD) in the fingers and toes during exposure to the cool environment (16 ± 1 °C) in the two groups. The first 20 min (00:00 to 00:20) indicate data collected during the baseline phase, the next five min (00:20 to 00:25) indicate responses during the warm-water immersion, the next 40 min (00:30 to 00:70) indicate responses during the cold-water immersion, and the final five min (00:70 to 00:75) show responses during the recovery phase. Finger and toe temperatures are indicated with continuous yellow lines in men and with dashed yellow lines in women. Finger and toe skin blood flow data are indicated with continuous purple lines in men and with dashed purple lines in women.

Figure A1

Skin temperature and blood flow (mean ± SD) in the fingers and toes during exposure to the thermoneutral environment (23 ± 1 °C) in the two groups. The first 20 min (00:00 to 00:20) indicate data collected during the baseline phase, the next five min (00:20 to 00:25) indicate responses during the warm-water immersion, the next 40 min (00:30 to 00:70) indicate responses during the cold-water immersion, and the final five min (00:70 to 00:75) show responses during the recovery phase. Finger and toe temperatures are indicated with continuous yellow lines in men and with dashed yellow lines in women. Finger and toe skin blood flow data are indicated with continuous purple lines in men and with dashed purple lines in women.

Figure A2

Skin temperature and blood flow (mean ± SD) in the fingers and toes during exposure to the hot environment (34 ± 1 °C) in the two groups. The first 20 min (00:00 to 00:20) indicate data collected during the baseline phase, the next five min (00:20 to 00:25) indicate responses during the warm-water immersion, the next 40 min (00:30 to 00:70) indicate responses during the cold-water immersion, and the final five min (00:70 to 00:75) show responses during the recovery phase. Finger and toe temperatures are indicated with continuous yellow lines in men and with dashed yellow lines in women. Finger and toe skin blood flow data are indicated with continuous purple lines in men and with dashed purple lines in women.

For the second finger, in the cool environment women demonstrated lower Tmin (d = 1.74) and Tmax (d = 0.56), delayed Tonset and Tpeak (very small to small effect sizes; d = 0.02 to 0.40), as well as negligible differences in Tavg and the amplitude in comparison to men. In the thermoneutral environment, women demonstrated lower Tmin, Tmax, Tavg, and amplitude (d = 0.82 to 1.43), as well as earlier Tonset and Tpeak (medium to huge effect sizes; d = 0.51 to 2.94) compared to men. In the hot environment, women demonstrated higher Tmin, Tmax, Tavg, and amplitude (medium to very large effect sizes; d = 0.65 to 1.74) as well as earlier Tonset (d = 0.78) compared to men, while Tpeak was similar between sexes (d = 0.06).

3.3. Cardiovascular Strain

Women demonstrated higher cardiovascular strain during the entire protocol across all environments in comparison to their male counterparts. Specifically, women presented higher HR values (medium to large effect sizes; d = 0.58 to 1.81; Table A2, Figure A3, Figure A4 and Figure A5) in comparison to men. Moreover, they also demonstrated higher diastolic blood pressure (very small to large effect sizes; d = 0.11 to 0.99) and lower systolic blood pressure (very small to large effect sizes; d = 0.05 to 0.73) compared to men. In the cool environment, during the first minute of cold-water immersion, women experienced higher cardiovascular strain (as expressed by delta differences from the last minute of warm-water immersion) in comparison to men (Figure 3).

Table A2

Physiological and perceptual parameters across groups and environments. Effect sizes indicate between group comparisons.

			Baseline	Warm Immersion	Cold Immersion	Recovery	Effect Size (d)
							Base	WI	CI	Rec
Gastrointestinal temperature (°C)	Cool	men	37.18 ± 0.17	37.24 ± 0.23	37.01 ± 0.32	36.94 ± 0.44	0.86	0.33	0.24	0.10
		women	37.30 ± 0.10	37.30 ± 0.11	37.07 ± 0.17	36.90 ± 0.15
	Neutral	men	37.15 ± 0.25	37.16 ± 0.27	37.00 ± 0.32	36.83 ± 0.55	0.32	0.36	0.88	0.73
		women	37.24 ± 0.28	37.28 ± 0.35	37.26 ± 0.27	37.17 ± 0.35
	Hot	men	37.03 ± 0.10	37.07 ± 0.16	37.14 ± 0.19	37.01 ± 0.26	1.40	1.01	0.76	0.19
		women	37.25 ± 0.20	37.31 ± 0.27	37.28 ± 0.18	37.05 ± 0.07
Mean skin temperature (°C)	Cool	men	30.34 ± 0.86	30.06 ± 0.96	29.55 ± 0.91	29.42 ± 1.06	0.23	0.10	0.01	0.18
		women	30.51 ± 0.61	30.15 ± 0.70	29.56 ± 0.82	29.59 ± 0.91
	Neutral	men	32.27 ± 0.68	32.20 ± 0.56	31.92 ± 0.53	31.84 ± 0.79	0.32	0.41	0.33	0.25
		women	32.48 ± 0.64	32.43 ± 0.56	32.09 ± 0.53	32.00 ± 0.43
	Hot	men	34.80 ± 0.97	34.90 ± 1.03	35.05 ± 1.06	35.08 ± 1.29	0.43	0.54	0.53	0.68
		women	35.13 ± 0.47	35.35 ± 0.52	35.47 ± 0.36	35.72 ± 0.31
Arm temperature (°C)	Cool	men	31.35 ± 1.01	31.09 ± 1.11	30.67 ± 1.08	30.29 ± 1.08	0.54	0.52	0.60	0.02
		women	30.85 ± 0.85	30.60 ± 0.72	30.06 ± 0.87	30.30 ± 0.78
	Neutral	men	33.20 ± 0.77	33.31 ± 0.64	33.04 ± 0.72	33.10 ± 1.03	0.31	0.39	0.19	0.21
		women	32.95 ± 0.85	33.04 ± 0.75	32.92 ± 0.54	32.93 ± 0.42
	Hot	men	35.58 ± 0.51	35.79 ± 0.51	35.97 ± 0.51	36.07 ± 0.50	0.33	0.33	0.53	0.44
		women	35.38 ± 0.71	35.60 ± 0.67	35.71 ± 0.49	35.88 ± 0.30
Chest temperature (°C)	Cool	men	30.67 ± 1.27	30.57 ± 1.38	30.23 ± 1.47	30.72 ± 1.82	1.49	1.26	1.32	1.12
		women	32.30 ± 0.88	32.11 ± 1.04	31.89 ± 1.01	32.32 ± 0.90
	Neutral	men	32.49 ± 0.74 0.009	32.48 ± 0.6 0.006	32.37 ± 0.56	32.50 ± 0.81	1.62	2.11	1.71	1.37
		women	33.60 ± 0.64	33.73 ± 0.58	33.44 ± 0.68	33.54 ± 0.70
	Hot	men	34.39 ± 2.05	34.36 ± 2.10	34.44 ± 2.25	34.29 ± 3.00	0.78	0.84	0.80	0.90
		women	35.56 ± 0.47	35.64 ± 0.45	35.73 ± 0.37	36.22 ± 0.36
Leg temperature (°C)	Cool	men	29.17 ± 1.07	28.59 ± 1.09	27.79 ± 1.20	27.10 ± 1.35	0.40	0.34	0.48	0.65
		women	28.71 ± 1.62	28.13 ± 1.58	27.02 ± 1.95	25.92 ± 2.20
	Neutral	men	31.20 ± 0.88	30.91 ± 0.83	30.28 ± 0.75	29.62 ± 0.73	0.24	0.28	0.42	0.47
		women	30.99 ± 0.81	30.67 ± 0.89	29.93 ± 0.95	29.23 ± 0.90
	Hot	men	34.43 ± 0.53	34.52 ± 0.68	34069 ± 0.63	34.87 ± 0.40	0.16	0.35	0.30	0.36
		women	34.36 ± 0.39	34.73 ± 0.55	34.84 ± 0.35	34.75 ± 0.24
Thigh temperature (°C)	Cool	men	29.47 ± 1.03	29.21 ± 1.16	28.63 ± 1.19	28.47 ± 1.45	0.33	0.59	0.65	0.31
		women	29.10 ± 1.25	28.52 ± 1.19	27.86 ± 1.17	28.10 ± 0.88
	Neutral	men	31.61 ± 0.72	31.43 ± 0.62	31.20 ± 0.70	31.19 ± 0.83	0.01	0.11	0.31	0.15
		women	31.60 ± 0.75	31.35 ± 0.72	30.98 ± 0.73	31.06 ± 0.75
	Hot	men	34.63 ± 1.25	34.76 ± 1.23	34.90 ± 1.42	34.98 ± 1.83	0.27	0.41	0.41	0.54
		women	34.89 ± 0.47	35.16 ± 0.58	35.34 ± 0.53	35.70 ± 0.50
1st toe immersed leg (°C)	Cool	men	21.48 ± 1.39	31.51 ± 2.48	10.10 ± 0.60	11.63 ± 0.87	0.02	0.40	0.49	0.59
		women	21.44 ± 2.27	30.34 ± 3.38	10.48 ± 0.93	12.73 ± 2.47
	Neutral	men	26.30 ± 3.24	33.08 ± 1.02	10.80 ± 1.17	13.92 ± 2.09	0.42	0.42	0.51	0.20
		women	28.05 ± 4.83	33.76 ± 2.05	10.32 ± 0.62	14.33 ± 1.93
	Hot	men	34.40 ± 0.83	35.47 ± 0.77	13.03 ± 2.59	23.10 ± 5.26	0.89	0.61	0.39	0.69
		women	35.01 ± 0.55	35.96 ± 0.86	12.19 ± 1.59	20.24 ± 2.65
3rd toe immersed leg (°C)	Cool	men	19.55 ± 1.55	32.17 ± 2.42	9.91 ± 0.79	12.05 ± 0.83	0.45	0.56	0.21	0.27
		women	20.46 ± 2.338	30.57 ± 3.24	10.10 ± 1.05	12.46 ± 1.92
	Neutral	men	25.15 ± 3.22	32.74 ± 1.11	10.62 ± 1.05	13.54 ± 1.35	0.41	0.93	0.59	0.68
		women	26.90 ± 5.07	33.87 ± 1.33	9.93 ± 1.30	14.81 ± 2.29
	Hot	men	34.43 ± 0.97	35.22 ± 0.56	12.40 ± 1.08	23.05 ± 4.78	0.34	0.38	0.67	0.34
		women	34.72 ± 0.74	35.47 ± 0.73	11.29 ± 2.05	21.43 ± 4.66
5th toe immersed leg (°C)	Cool	men	20.10 ± 1.71	31.66 ± 1.72	10.97 ± 1.88	12.19 ± 1.21	0.62	0.53	0.10	0.49
		women	21.74 ± 3.32	30.13 ± 3.69	10.81 ± 1.22	13.23 ± 2.79
	Neutral	men	25.54 ± 3.15	32.02 ± 1.23	10.91 ± 0.73	13.30 ± 1.81	0.41	1.38	0.57	0.95
		women	27.21 ± 4.76	33.68 ± 1.18	10.39 ± 1.07	15.40 ± 2.56
	Hot	men	34.02 ± 1.25	35.06 ± 1.44	13.28 ± 2.24	21.87 ± 5.01	0.49	0.39	0.01	0.31
		women	34.58 ± 1.00	35.54 ± 1.00	13.24 ± 4.26	20.46 ± 4.12
1st toe non immersed leg (°C)	Cool	men	21.82 ± 1.46	21.10 ± 1.06	18.31 ± 0.92	16.93 ± 0.91	0.03	0.17	0.13	0.20
		women	21.89 ± 2.33	21.38 ± 2.14	18.47 ± 1.38	17.12 ± 0.99
	Neutral	men	26.45 ± 3.85	26.67 ± 3.74	24.34 ± 2.97	24.33 ± 3.81	0.34	0.16	0.34	0.81
		women	27.79 ± 4.00	27.24 ± 3.59	23.52 ± 1.70	22.09 ± 0.89
	Hot	men	34.44 ± 0.91	34.80 ± 0.62	34.08 ± 0.94	34.48 ± 0.47	0.88	0.85	0.27	1.17
		women	35.10 ± 0.55	35.30 ± 0.56	34.32 ± 0.77	33.03 ± 1.69
2nd finger immersed hand (°C)	Cool	men	20.77 ± 4.30	30.74 ± 2.38	11.00 ± 1.18	14.25 ± 2.96	0.13	0.55	0.03	0.14
		women	21.35 ± 4.80	29.27 ± 3.00	11.04 ± 0.98	14.64 ± 2.48
	Neutral	men	29.37 ± 3.48	33.99 ± 1.24	11.87 ± 1.48	18.15 ± 5.32	0.70	0.36	0.45	0.25
		women	31.49 ± 2.43	34.51 ± 1.64	11.18 ± 1.59	17.13 ± 2.19
	Hot	men	35.46 ± 0.59	35.17 ± 0.44	13.02 ± 2.46	28.50 ± 4.83	1.16	1.95	0.69	0.02
		women	36.01 ± 0.40	36.01 ± 0.43	14.80 ± 2.68	28.41 ± 3.21
2nd finger immersed hand (°C)	Cool	men	21.04 ± 3.99	19.73 ± 4.02	17.36 ± 1.41	16.52 ± 0.53	0.12	0.24	0.15	0.07
		women	20.54 ± 4.19	18.81 ± 3.50	17.13 ± 1.59	16.57 ± 0.72
	Neutral	men	28.65 ± 3.17	27.85 ± 2.89	24.84 ± 2.09	24.07 ± 2.34	1.07	0.62	0.18	0.69
		women	31.32 ± 1.58	29.64 ± 2.92	25.25 ± 2.30	22.92 ± 0.30
	Hot	men	35.50 ± 0.59	35.51 ± 0.51	35.21 ± 0.59	35.39 ± 0.17	0.05	0.18	0.40	1.35
		women	35.53 ± 0.57	35.39 ± 0.74	34.93 ± 0.79	34.58 ± 0.83
SkBF 1st toe immersed leg (PU)	Cool	men	8.47 ± 5.29	39.54 ± 31.54	20.88 ± 11.09	38.74 ± 7.97	0.44	0.42	0.41	0.53
		women	17.76 ± 29.25	28.59 ± 19.62	28.05 ± 22.40	56.91 ± 47.73
	Neutral	men	52.14 ± 47.67	79.53 ± 28.06	20.68 ± 13.98	33.26 ± 15.27	0.37	0.43	0.71	0.90
		women	80.81 ± 98.23	102.68 ± 70.38	42.56 ± 41.34	60.48 ± 39.83
	Hot	men	262.48 ± 72.15	282.51 ± 67.05	82.22 ± 45.71	157.16 ± 108.53	0.12	0.92	0.01	0.29
		women	253.28 ± 78.81	210.34 ± 88.42	81.59 ± 57.26	131.52 ± 63.37
SkBF 1st toe non immersed leg (PU)	Cool	men	11.81 ± 5.15	14.02 ± 10.42	14.84 ± 8.88	12.50 ± 2.42	0.52	0.52	0.43	0.25
		women	28.11 ± 44.44	50.32 ± 98.45	21.51 ± 20.34	11.29 ± 6.38
	Neutral	men	63.75 ± 43.87	48.65 ± 38.23	28.82 ± 26.19	38.20 ± 33.55	0.03	0.07	0.50	1.08
		women	62.39 ± 56.62	45.74 ± 43.06	18.504.03	12.44 ± 1.33
	Hot	men	221.14 ± 75.13	233.12 ± 62.66	193.68 ± 63.98	210.73 ± 58.43	0.18	0.39	0.33	3.15
		women	203.14 ± 120.04	197.30 ± 112.78	160.92 ± 123.61	67.76 ± 26.62
SkBF 2nd finger immersed hand (PU)	Cool	men	89.75 ± 121.99	105.82 ± 110.72	103.67 ± 64.71	108.13 ± 13.86	0.14	0.13	0.27	0.33
		women	72.60 ± 116.16	91.30 ± 116.72	85.49 ± 71.23	139.03 ± 82.32
	Neutral	men	224.96 ± 126.19	202.69 ± 118.12	100.06 ± 62.93	178.11 ± 94.11	0.63	0.70	0.63	0.06
		women	293.38 ± 87.65	275.99 ± 89.71	137.65 ± 56.96	173.18 ± 83.76
	Hot	men	329.08 ± 75.45	301.39 ± 91.38	151.01 ± 66.33	290.01 ± 101.76	1.04	0.50	0.75	0.86
		women	392.06 ± 39.75	339.34 ± 56.45	199.13 ± 62.69	354.43 ± 28.27
SkBF 2nd finger non immersed hand (PU)	Cool	men	46.71 ± 54.65	40.38 ± 63.60	19.44 ± 8.97	15.34 ± 7.40	0.21	0.10	0.23	0.22
		women	59.57 ± 69.78	34.40 ± 50.08	2.75 ± 18.24	17.48 ± 11.46
	Neutral	men	148.45 ± 73.52	95.70 ± 54.05	44.06 ± 24.55	42.13 ± 21.85	0.76	0.55	0.24	0.63
		women	205 ± 0.37 ± 76.72	137.69 ± 94.24	54.56 ± 56.13	28.05 ± 22.58
	Hot	men	248.85 ± 78.47	253.97 ± 75.61	209.77 ± 68.63	246.72 ± 31.68	0.02	0.29	0.22	0.38
		women	247.17 ± 56.27	233.98 ± 61.27	210.44 ± 47.45	224.17 ± 76.81
Sweat rate forehead (mg/cm2/min)	Cool	men	1.20 ± 0.36	1.42 ± 0.30	1.51 ± 0.40	1.31 ± 0.61	0.27	0.28	0.11	0.57
		women	1.07 ± 0.58	1.28 ± 0.64	1.45 ± 0.76	1.74 ± 0.89
	Neutral	men	0.93 ± 0.19	1.14 ± 0.25	1.23 ± 0.35	1.01 ± 0.36	0.66	0.76	0.84	2.47
		women	1.75 ± 1.74	2.00 ± 1.56	2.15 ± 1.50	3.00 ± 1.08
	Hot	men	1.39 ± 0.71	1.57 ± 0.85	1.67 ± 0.82	1.64 ± 1.03	0.13	0.30	0.41	0.59
		women	1.49 ± 0.83	1.84 ± 0.91	2.03 ± 0.89	2.28 ± 1.13
Sweat rate leg (mg/cm2/min)	Cool	men	0.56 ± 0.38	0.61 ± 0.48	0.74 ± 0.46	0.78 ± 0.57	0.58	0.23	0.06	0.02
		women	0.80 ± 0.44	0.70 ± 0.28	0.76 ± 0.30	0.77 ± 0.33
	Neutral	men	0.75 ± 0.33	0.90 ± 0.47	0.89 ± 0.50	0.85 ± 0.66	0.66	0.23	0.56	0.26
		women	1.04 ± 1.30	0.76 ± 0.79	0.64 ± 0.55	0.71 ± 0.38
	Hot	men	1.34 ± 0.29	1.40 ± 0.42	1.55 ± 0.71	1.64 ± 0.95	0.15	0.30	0.02	0.91
		women	1.26 ± 0.71	1.53 ± 1.08	1.20 ± 0.54	0.99 ± 0.36
Heart rate (beats/min)	Cool	men	64.89 ± 4.62	64.31 ± 3.57	71.17 ± 6.33	68.46 ± 9.52	1.43	1.05	0.90	1.29
		women	75.68 ± 9.60	71.78 ± 9.45	77.76 ± 8.21	80.90 ± 9.74
	Neutral	men	65.46 ± 2.87	65.11 ± 3.04	68.63 ± 4.08	64.07 ± 5.26	1.81	1.55	0.73	0.58
		women	78.68 ± 9.90	76.15 ± 9.57	73.11 ± 7.61	68.98 ± 10.76
	Hot	men	72.48 ± 9.43	72.34 ± 7.84	72.53 ± 6.41	68.02 ± 4.22	0.87	0.94	0.62	0.66
		women	81.20 ± 10.48	81.52 ± 11.44	78.75 ± 12.63	75.38 ± 15.12
Systolic blood pressure (mmHg)	Cool	men	108.48 ± 8.17	109.71 ± 5.96	122.20 ± 7.48	127.00 ± 11.75	0.17	0.07	0.84	0.10
		women	110.17 ± 11.35	109.14 ± 9.04	131.95 ± 14.60	125.50 ± 18.86
	Neutral	men	106.45 ± 6.07	107.67 ± 11.34	126.31 ± 13.98	121.00 ± 8.25	0.44	0.70	0.42	0.73
		women	103.17 ± 8.49	100.83 ± 7.88	121.17 ± 10.45	113.75 ± 11.27
	Hot	men	102.83 ± 10.64	104.17 ± 8.73	111.03 ± 12.26	102.75 ± 12.50	0.05	0.25	0.46	0.22
		women	102.39 ± 7.16	101.50 ± 12.37	106.40 ± 7.34	105.00 ± 7.16
Diastolic blood Pressure (mmHg)	Cool	men	71.86 ± 8.58	71.43 ± 4.83 <0.001	86.52 ± 7.71	86.75 ± 12.50	0.66	0.94	1.08	0.28
		women	77.17 ± 7.59	75.86 ± 4.63	95.38 ± 8.64	92.50 ± 26.06
	Neutral	men	70.06 ± 8.16	68.33 ± 10.67	86.78 ± 9.83	78.28 ± 10.37	0.42	0.44	0.11	0.18
		women	73.00 ± 5.56	72.33 ± 7.09	87.75 ± 6.96	80.00 ± 9.20
	Hot	men	67.78 ± 9.22	66.00 ± 7.80	74.88 ± 10.11	67.25 ± 10.69	0.67	0.47	0.36	0.99
		women	72.83 ± 5.23	69.00 ± 4.60	77.53 ± 2.91	78.25 ± 11.41
Pain index hand (Likert)	Cool	men			7.09 ± 1.89				0.09
		women			7.32 ± 2.99
	Neutral	men			6.86 ± 2.16				0.15
		women			6.45 ± 3.32
	Hot	men			5.35 ± 1.96				0.59
		women			3.98 ± 2.62
Pain index foot (Likert)	Cool	men			7.58 ± 1.05				0.19
		women			7.15 ± 3.02
	Neutral	men			6.97 ± 1.98				0.21
		women			6.38 ± 3.31
	Hot	men			5.37 ± 1.91				0.47
		women			4.31 ± 2.62
Pain index hand (VAS)	Cool	men			6.84 ± 2.04				0.01
		women			6.81 ± 3.05
	Neutral	men			6.46 ± 2.33				0.07
		women			6.27 ± 306
	Hot	men			5.08 ± 2.21				0.34
		women			4.25 ± 2.65
Pain index foot (VAS)	Cool	men			7.10 ± 1.25				0.20
		women			6.63 ± 3.04
	Neutral	men			6.54 ± 1.91				0.08
		women			6.32 ± 3.15
	Hot	men			4.90 ± 1.78				0.24
		women			4.38 ± 2.53
Distress index hand (Likert)	Cool	men			7.20 ± 1.89				0.06
		women			7.37 ± 3.48
	Neutral	men			6.83 ± 2.44				0.03
		women			6.93 ± 3.61
	Hot	men			5.51 ± 1.95				0.52
		women			4.23 ± 2.88
Distress index foot (Likert)	Cool	men			7.54 ± 0.95				0.11
		women			7.26 ± 3.46
	Neutral	men			6.66 ± 1.86				0.09
		women			6.91 ± 3.62
	Hot	men			5.15 ± 2.04				0.26
		women			4.51 ± 2.77
Tactile sensitivity hand (monofil.)	Cool	men	0.62 ± 0.94	0.35 ± 0.73	26.14 ± 60.70		0.55	0.17	0.26
		women	1.42 ± 1.80	0.25 ± 0.25	51.55 ± 122.21
	Neutral	men	0.07 ± 0.00	0.13 ± 0.13	29.34 ± 59.73		0.69	0.47	0.67
		women	0.45 ± 0.77	0.39 ± 0.79	1.04 ± 0.87
	Hot	men	0.39 ± 0.79	0.07 ± 0.00	25.56 ± 60.99		0.36	0.58	0.50
		women	0.71 ± 1.00	0.39 ± 0.79	3.64 ± 7.41
Tactile sensitivity foot (monofil.)	Cool	men	0.72 ± 0.89	1.72 ± 0.73	77.01 ± 124.08		0.66	0.01	0.01
		women	1.51 ± 1.43	1.73 ± 1.35	76.25 ± 124.78
	Neutral	men	0.82 ± 0.92	1.15 ± 0.94	76.19 ± 124.66		0.00	0.22	0.20
		women	0.82 ± 0.92	1.42 ± 1.55	51.73 ± 121.63
	Hot	men	1.41 ± 0.92	0.88 ± 0.88	1.63 ± 0.48		0.91	0.51	0.58
		women	100.52 ± 154.52	1.36 ± 1.00	51.38 ± 121.80
Tactile sensitivity hand (esthesiometer)	Cool	men	2.57 ± 0.74	3.46 ± 0.92	4.37 ± 1.96		0.37	0.40	0.74
		women	3.01 ± 1.52	2.99 ± 1.41	2.97 ± 1.83
	Neutral	men	3.11 ± 1.27	3.32 ± 1.11	3.76 ± 1.55		0.10	0.49	0.47
		women	2.99 ± 1.08	2.68 ± 1.44	2.94 ± 1.93
	Hot	men	2.80 ± 0.74	2.46 ± 0.81	3.02 ± 0.55		0.63	0.19	0.12
		women	2.17 ± 1.20	2.26 ± 1.25	2.92 ± 1.08
Tactile sensitivity Foot (esthesiometer)	Cool	men	8.28 ± 2.52	6.82 ± 2.08	7.75 ± 2.31 0.005		0.70	0.29	0.80
		women	6.32 ± 3.08	6.15 ± 2.51	6.03 ± 2.00
	Neutral	men	5.59 ± 1.05	5.02 ± 2.56	4.93 ± 2.39		0.63	0.15	0.53
		women	6.97 ± 2.91	5.43 ± 2.91	6.58 ± 365
	Hot	men	7.30 ± 1.39	6.40 ± 2.05	6.13 ± 2.41		0.36	0.56	0.01
		women	6.43 ± 3.18	4.86 ± 3.31	6.15 ± 2.86
Thermal comfort	Cool	men	2.31 ± 0.94	2.36 ± 1.11	2.97 ± 0.96	2.13 ± 0.63	0.29	0.06	0.13	0.93
		women	2.05 ± 0.84	2.29 ± 1.25	3.13 ± 1.45	3.38 ± 1.80
	Neutral	men	1.33 ± 0.30	1.33 ± 0.61	1.47 ± 0.43	1.88 ± 1.11	0.18	0.32	0.63	0.49
		women	1.25 ± 0.61	1.17 ± 0.41	1.88 ± 0.82	2.38 ± 0.95
	Hot	men	1.72 ± 0.49	2.00 ± 0.89	1.58 ± 0.57	1.38 ± 0.25	0.10	0.16	0.64	0.00
		women	1.78 ± 0.64	2.00 ± 0.89	1.29 ± 0.29	1.38 ± 0.25
Thermal sensation	Cool	men	1.81 ± 0.74	1.43 ± 0.98	1.84 ± 0.57	1.50 ± 1.91	0.17	0.42	0.98	0.65
		women	1.67 ± 0.90	1.86 ± 1.07	2.50 ± 0.78	2.50 ± 1.00
	Neutral	men	0.56 ± 0.46	0.17 ± 0.98	0.19 ± 0.38	1.00 ± 1.15	0.82	0.16	0.80	0.62
		women	0.22 ± 0.35	0.00 ± 1.10	0.67 ± 0.75	1.75 ± 1.26
	Hot	men	1.72 ± 0.49	1.67 ± 0.82	1.01 ± 0.74	1.50 ± 0.58	0.14	0.58	0.44	0.26
		women	1.78 ± 0.27	0.00 ± 1.10	1.40 ± 1.00	1.25 ± 1.26
Key: Superscript numbers indicate p values for statistically significant differences at p < 0.008 from paired-samples t tests with a Bonferroni correction.
Note: Base = baseline; WI = warm-water immersion; CI = cold-water immersion; Rec = recovery; SkBF = skin blood flow; VAS = visual analogue scale.
Effect size (d) interpretation:			very small0.010.19	small0.20.49	medium0.50.79	large0.81.19	very large1.21.99			huge≥2.0
higher values in men
lower values in men

Figure A3

Physiological parameters (mean ± SD) during exposure to the cool environment (16 ± 1 °C) in the two groups. The first 20 min (00:00 to 00:20) indicate data collected during the baseline phase, the next five min (00:20 to 00:25) indicate responses during the warm-water immersion, the next 40 min (00:30 to 00:70) indicate responses during the cold-water immersion, and the final five min (00:70 to 00:75) show responses during the recovery phase. Continuous yellow lines indicate results for men participants and with dashed yellow lines present results for the women individuals.

Figure A4

Physiological parameters (mean ± SD) during exposure to the thermoneutral environment (23 ± 1 °C) in the two groups. The first 20 min (00:00 to 00:20) indicate data collected during the baseline phase, the next five min (00:20 to 00:25) indicate responses during the warm-water immersion, the next 40 min (00:30 to 00:70) indicate responses during the cold-water immersion, and the final five min (00:70 to 00:75) show responses during the recovery phase. Continuous yellow lines indicate results for men participants and with dashed yellow lines present results for the women individuals.

Figure A5

Physiological parameters (mean ± SD) during exposure to the hot environment (34 ± 1 °C) in the two groups. The first 20 min (00:00 to 00:20) indicate data collected during the baseline phase, the next five min (00:20 to 00:25) indicate responses during the warm-water immersion, the next 40 min (00:30 to 00:70) indicate responses during the cold-water immersion, and the final five min (00:70 to 00:75) show responses during the recovery phase. Continuous yellow lines indicate results for men participants and with dashed yellow lines present results for the women individuals.

3.4. Other Physiological Responses

The physiological parameters monitored during the protocol are shown in Appendix A (Table A2; Figure A3, Figure A4 and Figure A5). In most cases, the differences observed between groups did not reach the Bonferroni adjusted threshold of p < 0.008 for post hoc comparisons, yet a number of moderate or large effect sizes were detected. Specifically, the Tgi was lower for men compared to women in the cool environment during baseline (d = 0.86), in the thermoneutral environment during cold-water immersion (d = 0.88) and recovery (d = 0.73), as well as in the hot environment during baseline (d = 1.40), warm-water immersion (d = 1.01), and cold-water immersion (d = 0.76).

No substantial differences were identified between the Tsk of men and women (very small to medium effect sizes; d = 0.01 to 0.68), with the female group demonstrating a slightly increased Tsk (Table A2). Nevertheless, it is important to note that across all environments and phases, women had higher chest skin temperature (large to very large effect sizes; d = 0.80 to 2.11) and lower limb (gastrocnemius) temperature (small to medium effect sizes; d = 0.01 to 0.65) compared to men (Table A2). On the other hand, in most cases, similar Tt and SkBF in the immersed foot were observed between men and women across all environments (d = 0.02 to 0.69; Table A2).

No marked and/or meaningful differences were observed between men and women in the sweat rate of the forehead and the leg (gastrocnemius), except for the thermoneutral environment where women showed markedly higher forehead sweat rate (medium to huge effect sizes, d = 0.66 to 2.47) and lower leg sweat rate for the hot environment during recovery (d = 0.91).

3.5. Perceptual Data

Pain and distress were similar between men and women across all environments (very small to medium effect sizes; d = 0.01 to 0.59; Table A2). Compared to men, the tactile sensitivity of women tended to be slightly decreased when assessed using the Semmes–Weinstein monofilaments (very small to medium effect sizes; d = 0.01 to 0.69) and slightly increased when assessed using a digital esthesiometer (very small to large effect sizes; d = 0.01 to 0.80; Table A2). This difference suggests that adoption of both methods provides a more robust assessment of tactile sensitivity in relation to CIVD. During the cold-water immersion and the recovery, women reported being more thermally uncomfortable compared to men in the cool and the thermoneutral environments (small to large effect sizes; d = 0.13 to 0.93), and this was reversed in the heat (no to medium effect sizes; d = 0.00 to 0.64). Finally, in terms of thermal sensation, men reported feeling slightly warmer than women in the cool environment (small to large effect sizes; d = 0.17 to 0.98), and these results were reversed in the hot environment (small to medium effect sizes; d = 0.14 to 0.58; Table A2).

4. Discussion

For a matched thermal stress, women in the present study experienced an elevated cardiovascular strain and higher frequency of CIVD reactions in comparison to their male counterparts. Specifically, the female group had higher HR and twice as many CIVD reactions compared to their matched male counterparts during exposure across all environments, while they had increased mean arterial pressure in the cool environment and almost the same mean arterial pressure in the other two environments (i.e., thermoneutral and hot). When both fingers and toes were considered, we found no statistically significant effects of sex in the overall CIVD occurrence. However, females demonstrated a statistically higher occurrence of CIVD was in the toes across all environments accompanied by an anticipated increase in the local SkBF of the immersed foot, which may be triggered as a response to the lower Tt.

The few previous studies [10,26,27,28] evaluating sex differences in the occurrence of CIVD indicate no systematic difference between women and men. However, these studies investigated the occurrence of CIVD responses only in the fingers of their participants, and their results are in line with our findings of no substantial sexual differences in CIVD reactions in the fingers. In contrast, we show for the first time that women experienced twice as many CIVD reactions compared to men in the toes. This may be explained through two theoretical frameworks. First, previous studies state that there are potential differences in vasomotion between men and women due to hormonal status [49] or due to an increase in sympathetic outflow to the cutaneous circulation in women [23]. Furthermore, other studies indicated that there are no sex differences in cutaneous vasodilation [50], and that the sex differences are probably due to epidermis thickness [51] or the size of the hand and foot rather than due to the hormonal status [52]. Thus, the second, and probably more important, theory is that these differences may be a result of the smaller surface area which characterizes the immersed feet of women [52]. As a rule, the average woman (531 cm2) has approximately an 18% smaller foot compared to the average man (650 cm2) [53], which means that female feet may cool more quickly when exposed to cold stimuli, a notion that was previously proposed by others [52]. This may, therefore, affect local mechanisms that have been reported to induce CIVD [8], resulting in more frequent CIVD reactions.

Modern workplace practices involve females in the same work tasks as men, and female workers can be currently found doing jobs requiring frequent immersion of their limbs to cold water, such as butchering [30], fishing [31], and performing army-related tasks [32,33]. In the coming decades, more women are expected to enter the world of male-dominated work; therefore, more studies should be directed toward identifying workplace strategies for achieving gender equity. In this light, the literature suggests that sex should be considered an important contributing factor capable of modifying the physiological heat strain experienced by workers [54,55,56]. Yet, the majority of relevant studies investigate the impact of occupational heat-stress, not cold exposure. Only a handful of studies on army-related tasks or occupations such as fishing or butchering have assessed the risk for cold injuries in the extremities during exposure to cold [57,58,59]. Some of these studies found that women have increased risk for cold injuries compared to men [60,61] but none of them mention protective footwear for occupations such as fishing, butchering, and especially the disparate cold injury risks for males and females. These sex differences were further confirmed in the present study, indicating that foot clothing and shoe industries should contemplate implementing new technologies that will consider sex differences in an attempt to better protect worker health and well-being. In addition to this, during exposure to cool environments, female workers may be placed in jobs that do not require frequent exposure of their limbs to cold water. It is important to note that, although this is the first study to examine sex differences in cardiovascular stress and CIVD responses in feet and hands during cold-water immersion, it involved monitoring a small group of Caucasian individuals that were not habituated to frequent exposure of their limbs to cold-water immersion. Future studies should consider investigating how inter- and intra- individual factors such as age, ethnicity, and body composition may contribute to the occurrence of the CIVD responses in the feet. Finally, it is worth noting that the lack of heterogeneity and sex differences in CIVD reactions in the hands is likely due to the single measurement site, contrary to the toes where we used three sites and reported significant sex differences. Nevertheless, it is important to note that all previous studies on the topic reported no effects of sex on CIVD in the fingers [10,26,27,28].

5. Conclusions

The present investigation demonstrated that the studied women experienced elevated cardiovascular strain and higher frequency of CIVD reactions, particularly in the toes, compared to their male counterparts, especially when exposed to cool environments.