A Database of Static Thermal Insulation and Evaporative Resistance Values of Dutch Firefighter Clothing Items and Ensembles

Simple Summary Rescue services personnel may be exposed to a wide range of harsh environments that require the use of protective gear. In order to achieve the best performance, the optimal protection for specific incident scenarios should be selected. In order to be able to specify the protective clothing ensemble that best matches the protection needs, activity and environmental conditions, while keeping up the performance and reduce impact of thermal stress, we need to know the clothing properties that affect heat transfer via the protective layers. The thermal insulation and evaporative resistance values of firefighter clothing items and ensembles of the whole protective system are not easily available. The work presented in this paper fills this gap. Collected data allow for validation and selection of prediction models for exposure evaluation, and with this, it contributes to the best performing protective ensemble choice for specific incident scenarios. The paper also presents the relationships for summing individual clothing items’ insulation of the firefighter protective clothing system and therefore reduces the need for separate testing of all possible clothing ensemble configurations. The wide literature background allows for extrapolation of the results to other branches where protective clothing is used. Abstract The rescue operations’ environment can impair firefighters’ performance and increase the risk of injuries, e.g., burns and hyperthermia. The bulk and carried weight of heavy protection contributes to lower physical performance, higher metabolic load and internal body heat production. For recommending optimal protection for the tasks and incident scenarios, knowledge of clothing thermal properties is needed. However, detailed data on firefighter protective clothing systems are not available. The aim of the study was to provide scientific background and a dataset that would allow for validation of thermo-physiological models for task-specific conditions of rescue work. Thermal insulation of 37 single items and their variations and 25 realistic protective clothing ensembles were measured on a thermal manikin. Twelve (12) ensembles that evenly covered the whole insulation range were selected for evaporative resistance testing. The equations for summing up individual item’s insulation to ensemble insulation and calculating clothing area factor were derived from the dataset. The database of a firefighter clothing system was created. In addition, the local and regional thermal properties of the clothing ensembles were provided for use in future validation of advanced thermo-physiological models for rescue worker exposure predictions and for designing decision aid tools.


Introduction
During rescue operations, the environment can impair firefighters' performance and increase the risk of injuries such as burns [1,2] and hyperthermia [3][4][5]. The firefighter ensembles offer high protection levels. However, increased protection also increases the bulk and carried weight that means lower physical performance, higher metabolic load and internal body heat production [6][7][8]. Reduced mobility, flexibility, and dexterity may also restrict the ability to successfully accomplish a task [9][10][11][12]. Therefore, an objective aim of firefighter organizations is being able to select optimal protection for the tasks and incident scenarios.
In the Netherlands, the fire services are currently switching to a new station wear, called operational uniform (OU) for low-risk tasks. It is very important to look at how firefighters can in the future be optimally protected within the entire risk spectrum of both low-risk and high-risk tasks. A current limitation is an insufficient database and knowledge about the thermal properties of firefighter protective clothing items and systems. There is currently no coordination of (protective) clothing packages in use for operational tasks. At every alarm, the protective clothing for firefighting is put on, with or without the current station uniform, or in the future, on the operational uniform. For volunteer firefighters, often, only the protective clothing for firefighting is available, and this is worn over their daily clothing. Risk profiles, protection requirements and wearing comfort play an important role in order to arrive at the protection concept that is best suited to the tasks. It is desirable that such a protection concept leaves space for integration of future developments.
With regard to functionality and protection, the minimum level has been established and formalized by means of (inter)national standards (NEN-EN and ISO) [13][14][15][16][17][18][19][20][21][22]. However, most of the requirements are related to materials and material packages and not for the complete items or clothing ensembles. The current clothing systems in use are combined where possible and should ultimately lead to a fully integrated protective clothing system, of which digital support can be a part, for example, mobile applications, such as ClimApp [23,24] and decision aids [25,26]. In the future, possibly, the clothing system can be combined with the necessary sensors and smart solutions to acquire real-time input from the incident site to allow for the optimal selection of gear. Firefighters have to be adequately protected for the tasks they have to perform, whereby each task may require specific protection. An aim is that the clothing can be built up by modules into an integral protection concept. Modularity in clothing systems may also support sustainable development with easier inspection, care and maintenance of the gear, allowing for replacement of worn out items and longer lasting of the whole protective clothing system.
Clothing insulation and evaporative resistance are the basic clothing parameters that affect human heat exchange with the environment. They are used as a common behavioural thermoregulatory measures and are means to maintain thermal comfort in a wide range of temperatures [27]. Commonly, these properties are not evenly distributed over the whole body surface due to the material choice, clothing design, layering, fit, etc. Physiological parameters, such as heat generation in muscles of different body parts during different work tasks, body internal heat distribution and areas for heat loss [28,29], protection criteria of vital regions, sweating patterns [30][31][32], but also freedom of motion, may set the specific demands on the need of regional evaporative resistance [33] and clothing insulation [34]. Clothing parameters for body regions can be measured, and the change due to walking and wind can be estimated [35][36][37][38]. Local values can be used in advanced physiological models for exposure evaluation [26,[39][40][41][42][43][44], for occupational health and safety purposes, as a feedback for clothing manufacturers for improving their products, and/or for industries for selecting the clothing provided to their employees. In addition, some publications provide information on how to count for fit, garment size and/or airgaps [37,43,45,46]. The methods to measure or estimate clothing items' or their combinations' insulation [47][48][49][50] and evaporative resistance [51,52] are standardized, or detailed suggestions for testing and data treatment are available [53]. ISO 9920 [50] contains databases on clothing items and ensembles based on earlier studies [54][55][56], and presents algorithms to sum individual items, calculate effect of wind and motion, etc. Despite the extensive content of this standard [50], it has a limited coverage of different clothing styles other than Western office and workwear, while recent publications cover that gap [57][58][59]. Another problem is that the majority of the standard's database was compiled mainly around 1980s, and both the design and the materials have changed. This has now been overcome by many recent studies on modern office and casual clothing that have been measured using the latest technology and providing detailed information on different aspects of clothing thermal performance [36][37][38][39]42,43,[60][61][62][63]. At the same time, when modern office and casual clothing are measured in a systematic way, data on industrial workwear are not as easily available. Some exceptions are studies on military [64][65][66][67][68][69], ambulance [70,71], cold protective [72][73][74][75][76][77], and some agricultural [78,79] clothing systems. Although, the studied clothing systems cover a wide variety of clothing configurations, it may be suspected that specialized clothing, e.g., the ones used by firefighters, may behave differently due to the protective layers and their combinations. Considering any future developments in clothing science and advanced support systems [23,25,26], and to select the best suited protection for any incident scenario in high stress situations quickly, more detailed information on thermal properties of clothing that is used by rescue workers is needed. As it is quite impossible to measure all possible clothing combinations used in the field, it is important to create an individual clothing items' wardrobe database based on manikin measurements. Development and/or validation of specific clothing insulation and evaporative resistance summation methods that match the selected clothing system would allow for evaluation of human exposure and to select the proper protection over a wide range of environments and working conditions [23,26,40,78,[80][81][82][83][84][85].
This work aimed to create a detailed firefighter clothing database of static thermal insulation values of the individual clothing items, some PPE items, and their combinations that are the elements of Dutch firefighter station wear (operational uniform) and incident scenario-based turnout gear. An aim was also to report the static thermal insulation of realistic clothing ensembles of station wear and turnout gear and evaporative resistance of selected clothing ensembles. As the collected material is quite extensive, local insulation and evaporative resistance values of the clothing ensembles for various body regions that may be useful for advanced thermo-physiological modelling are included in Appendices A and B (Tables A1 and A2, respectively). The items' and ensembles' thermal properties are ranked and visualized in Appendix C ( Figures A1-A3).

The Thermal Manikin Tore
In this study, all items were tested individually or in combinations on the thermal manikin Tore at the Thermal Environment Laboratory, Lund University. Tore is a maleshaped thermal manikin [86,87]. Tore is made of orthopaedic foam with a metal frame inside to support the body parts and for joints. Tore is divided into 17 individually controlled zones: head, chest, back, stomach, buttocks, left and right upper arm, left and right lower arm, left and right hand, left and right thigh, left and right calf, and left and right foot. The weight of the manikin is 32 kg. The surface temperature of the manikin's zones was controlled at 34 • C. Thermal manikin Tore and the used climatic chamber have been used in multiple joint studies and interlaboratory trials, and comparative data with the other manikins and test environments are well available over the wide range of clothing types, postures, walking and wind velocities, and for the variety of purposes, e.g., for thermal insulation [58,63,72,88], evaporative resistance [89], for thermo-physiological modelling [90], and as a part of evaluation of occupational incident [91], cooling systems [92,93], etc.

Manikin Calibration
Offset calibration of the manikin's surface temperature sensors was carried out in the homogenous conditions (34 • C) in a warm chamber with the same Pico equipment that was

Measurements of Clothing Surface Area Factor (f cl )
Photographic method based on 2 photos was used to estimate clothing area factor (f cl ) of the individual garments and the ensembles [52,58]. The frontal and the side photos were used.

Thermal Insulation Measurements
The tests for total thermal insulation (I T ) were carried out on the manikin in a standing position according to ISO 15831 [49] following ISO 9920 [50] recommendations (low air velocity). The basic (intrinsic) insulation of each garment item (I clu ) and clothing ensemble (I cl ) was calculated. The mean air velocity was measured with Swema Air 3000d logger with omnidirectional SWA03 sensor 7 Air velocity in the chamber stayed on average at 0.18 ± 0.07 m/s. It was originally measured as an average of 15 points in air flow 50 cm before the manikin and then checked in the control position at the end of each individual test. Air layer insulation (I a ) was measured with a nude manikin in the same conditions. It has been shown that the calculation method (parallel or serial) may strongly affect the result [76,[94][95][96][97]. When testing highly insulative and less permeable clothing, the manikin heat loss from some zones, especially in the areas with many overlapping layers, e.g., buttocks, abdomen, may drop close to zero, and local temperature may rise above the set temperature. Therefore, the principles of global method [76,97] that are insensitive for such technical effects and are also more reliable for correlating with human exposures, were utilized in the calculations.
The air temperature in the climate chamber was set to 20 • C for individual items and light clothing combinations and to 10 • C for heavy protective combinations. Air temperature variation during the stable state was below 0.1 • C. Ambient air temperature was continuously monitored using three sensors (PT 100, Pico Technology Ltd., St. Neots, UK) positioned adjacent to the ankles, the mid-trunk, and the head (vertical heights of 0.1, 1.1, and 1.7 m from the soles of the manikin). The mean radiant temperature was equal to the air temperature. The temperatures and heat losses were recorded at ten-second intervals. The last 10 min of stable state was used for thermal insulation calculation.

Evaporative Resistance Measurements
The total evaporative resistance (R et ) of selected clothing ensembles was measured according to ASTM F2370-16 [51] on the manikin in standing position under so-called isothermal conditions set to 34 • C. Both the mass loss and heat loss methods were used for further analysis, but as comparable local evaporative resistances of body regions could only be based on heat loss, only the results from the heat loss method are presented in this paper. Average air velocity during these tests was 0.40 ± 0.13 m/s. Data for the whole measuring period (commonly up to 70 min) were saved. Based on the weight change, temperature and heat loss curves, the stable periods, meaning constant evaporative heat loss, were selected for evaporative resistance calculation. Relative humidity was measured with the Humidity Sensor Evaluation Kit EK-H3 with pin SHT75 sensors (Sensirion AG, Stäfa, Switzerland). Relative humidity in the chamber was set to 40%. During the stable state of each measurement, the variation in the relative humidity was less than 2%. The corrections according to Wang et al. [98] were applied to the measured evaporative resistance values. Any further corrections, e.g., for skin material and moist skin insulation, etc. [52,53,[99][100][101], were not used. The further comparisons between different available corrections [53,102] will be presented in a separate paper. Finally, clothing evaporative resistance (R ecl ) and permeability indices (i m and i m,cl ) were calculated.

Calculations
Total insulation of an individual zone (I T,i , m 2 K/W): where T s,i -mean surface temperature of a manikin zone i ( • C); T o -operative temperature ( • C); A i -area of a manikin zone i (m 2 ); H c,i -dry heat loss from a manikin zone i (W). An individual zone insulation as a contribution to the whole body insulation (I T,i,r , m 2 K/W): where A-area of the whole manikin surface (m 2 ). Total insulation for a body zone group (I T,z,i , m 2 K/W): where ∑ A n z,i -area of all manikin zones included in the group (m 2 ); ∑ H n c,z,i -dry heat loss from all manikin zones included in the group (W).
Total insulation for the whole body (I T , m 2 K/W): Basic (intrinsic) clothing insulation of an item (I clu , m 2 K/W) or an ensemble (I cl , m 2 K/W): where I a -air layer insulation (measured on the nude manikin, m 2 K/W); f cl -clothing area factor (n.d.).
Water vapour pressure at manikin surface temperature (Pa): where T s -mean surface temperature of a manikin ( • C); RH sk -relative humidity at manikin (textile) skin temperature (%). Manikin surface temperature ( • C): Water vapour pressure at textile skin temperature (Pa): where T sk -mean (textile) skin temperature of a manikin ( • C), and where mean RH sk over the various conditions irrespective of the clothing combination was 96%.
where T sk,i -mean (textile) skin temperature of a manikin zone i ( • C). Correction of the manikin skin temperature [98,100] (depends on "skin" properties): where Q i -heat flux from a manikin zone i (W/m 2 ). Water vapour pressure in the air (Pa): where T a -air temperature ( • C); RH a -relative air humidity (%). Total evaporative resistance of a zone (m 2 Pa/W): where p sk,i -water vapour pressure at (textile) skin of a manikin zone i (Pa); H i -heat loss from a manikin zone i (W). Heat flux (W/m 2 ): Uncorrected evaporative resistance (m 2 Pa/W): where H e,i -assumed evaporative heat loss from a manikin zone i (W).
The assumption here is that conductive, convective and radiative heat exchange is close to zero, and the test is carried out at the isothermal conditions (T a = T s = T sk ), i.e., the only heat exchange pathway is evaporation.

Clothing
In total, 37 individual items or their variations (see Table 1 for images and description) and 25 clothing ensembles were measured including 150 dry measurements for thermal insulation and 36 wet measurements for evaporative resistance of 12 selected ensembles (see Table 2 for description and Table 3 for images). Each item and combination was tested at least twice. If the difference of two measurements was above 4%, then the measuring settings were checked, and additional measurements were carried out. In addition, the tested items and clothing combinations were photographed from front and side for documentation purposes and for clothing area factor (f cl ) calculations. All the items were weighed.        Table 2. Clothing ensembles and their properties * (for images of the clothing ensembles see Table 3). Table 1 Table 3. Images of the tested clothing ensembles (for properties and composition, see Table 2).           Table 3. Images of the tested clothing ensembles (for properties and composition, see Table 2).

Results and Discussion
The specifications and the insulation measurement results for the individual items (IT, Iclu) are shown in Table 1. Figure A1 in Appendix C visualizes the individual items' basic insulation ranked by the item group and/or body coverage region, e.g., feet, hands, head, underwear, shirts, jackets and their variations, etc. Table 2 presents the thermal insulation of the tested clothing ensembles (IT, Icl), the evaporative resistances (Ret,raw, Ret, Ret,cl) and permeability indices (im, im,cl) of the selected clothing ensembles. In addition, the measured clothing area factors (fcl) for individual items and the clothing ensembles are given in Tables 1 and 2, respectively. Table 3 presents the pictures of the clothing ensembles. Table A1 of Appendix A shows the thermal insulation and Table A2 of Appendix B the evaporative resistance values of the individual body zones and body regions.

Relationship between Measured and Summed Basic Clothing Insulation (Icl)
If the standard method for summing up the individual clothing items' insulation [50] fit practically perfectly for the ambulance clothing system [50,70], then the firefighter clothing system tested in the present study differed clearly in this respect, creating two groups. For the operational uniform (OU) ensembles, the standard summation method overestimated the actual insulation (Figure 1a), and for the incident clothing systems, it underestimated the insulation (Figure 1b). It may be expected that the basic clothing insulation range between 1.4 and 2 clo can be predicted well; however, there were no ensembles in the range from 1.6 to 1.9 clo. When these OU and incident clothing groups were separated, then for clothing ensembles that covered incident scenarios, it was enough to adjust the intercept of the standard equation to acquire a reasonable prediction (Equation (20), Figure 1b).

Results and Discussion
The specifications and the insulation measurement results for the individual items (IT, Iclu) are shown in Table 1. Figure A1 in Appendix C visualizes the individual items' basic insulation ranked by the item group and/or body coverage region, e.g., feet, hands, head, underwear, shirts, jackets and their variations, etc. Table 2 presents the thermal insulation of the tested clothing ensembles (IT, Icl), the evaporative resistances (Ret,raw, Ret, Ret,cl) and permeability indices (im, im,cl) of the selected clothing ensembles. In addition, the measured clothing area factors (fcl) for individual items and the clothing ensembles are given in Tables 1 and 2, respectively. Table 3 presents the pictures of the clothing ensembles. Table A1 of Appendix A shows the thermal insulation and Table A2 of Appendix B the evaporative resistance values of the individual body zones and body regions.

Relationship between Measured and Summed Basic Clothing Insulation (Icl)
If the standard method for summing up the individual clothing items' insulation [50] fit practically perfectly for the ambulance clothing system [50,70], then the firefighter clothing system tested in the present study differed clearly in this respect, creating two groups. For the operational uniform (OU) ensembles, the standard summation method overestimated the actual insulation (Figure 1a), and for the incident clothing systems, it underestimated the insulation (Figure 1b). It may be expected that the basic clothing insulation range between 1.4 and 2 clo can be predicted well; however, there were no ensembles in the range from 1.6 to 1.9 clo. When these OU and incident clothing groups were separated, then for clothing ensembles that covered incident scenarios, it was enough to adjust the intercept of the standard equation to acquire a reasonable prediction (Equation (20), Figure 1b).

Results and Discussion
The specifications and the insulation measurement results for the individual items (IT, Iclu) are shown in Table 1. Figure A1 in Appendix C visualizes the individual items' basic insulation ranked by the item group and/or body coverage region, e.g., feet, hands, head, underwear, shirts, jackets and their variations, etc. Table 2 presents the thermal insulation of the tested clothing ensembles (IT, Icl), the evaporative resistances (Ret,raw, Ret, Ret,cl) and permeability indices (im, im,cl) of the selected clothing ensembles. In addition, the measured clothing area factors (fcl) for individual items and the clothing ensembles are given in Tables 1 and 2, respectively. Table 3 presents the pictures of the clothing ensembles. Table A1 of Appendix A shows the thermal insulation and Table A2 of Appendix B the evaporative resistance values of the individual body zones and body regions.

Relationship between Measured and Summed Basic Clothing Insulation (Icl)
If the standard method for summing up the individual clothing items' insulation [50] fit practically perfectly for the ambulance clothing system [50,70], then the firefighter clothing system tested in the present study differed clearly in this respect, creating two groups. For the operational uniform (OU) ensembles, the standard summation method overestimated the actual insulation (Figure 1a), and for the incident clothing systems, it underestimated the insulation (Figure 1b). It may be expected that the basic clothing insulation range between 1.4 and 2 clo can be predicted well; however, there were no ensembles in the range from 1.6 to 1.9 clo. When these OU and incident clothing groups were separated, then for clothing ensembles that covered incident scenarios, it was enough to adjust the intercept of the standard equation to acquire a reasonable prediction (Equation (20), Figure 1b).

Results and Discussion
The specifications and the insulation measurement results for the individual items (I T , I clu ) are shown in Table 1. Figure A1 in Appendix C visualizes the individual items' basic insulation ranked by the item group and/or body coverage region, e.g., feet, hands, head, underwear, shirts, jackets and their variations, etc. Table 2 presents the thermal insulation of the tested clothing ensembles (I T , I cl ), the evaporative resistances (R et,raw , R et , R et,cl ) and permeability indices (i m , i m,cl ) of the selected clothing ensembles. In addition, the measured clothing area factors (f cl ) for individual items and the clothing ensembles are given in Tables 1 and 2, respectively. Table 3 presents the pictures of the clothing ensembles. Table A1 of Appendix A shows the thermal insulation and Table A2 of Appendix B the evaporative resistance values of the individual body zones and body regions.

Relationship between Measured and Summed Basic Clothing Insulation (I cl )
If the standard method for summing up the individual clothing items' insulation [50] fit practically perfectly for the ambulance clothing system [50,70], then the firefighter clothing system tested in the present study differed clearly in this respect, creating two groups. For the operational uniform (OU) ensembles, the standard summation method overestimated the actual insulation (Figure 1a), and for the incident clothing systems, it underestimated the insulation (Figure 1b). It may be expected that the basic clothing insulation range between 1.4 and 2 clo can be predicted well; however, there were no ensembles in the range from 1.6 to 1.9 clo. When these OU and incident clothing groups were separated, then for clothing ensembles that covered incident scenarios, it was enough to adjust the intercept of the standard equation to acquire a reasonable prediction (Equation (20), Figure 1b).
When the OU ensembles are above 1.3-1.4 clo, the prediction (summation) may be considered as a reasonable match ( Figure 1); then, for lighter clothing ensembles, the overestimation by summation is bigger. For OU ensembles, a better matching relationship had to be developed (R 2 = 0.993, Equation (21), Figure 1a): When the OU ensembles are above 1.3-1.4 clo, the prediction (summation) may be considered as a reasonable match ( Figure 1); then, for lighter clothing ensembles, the overestimation by summation is bigger. For OU ensembles, a better matching relationship had to be developed (R 2 = 0.993, Equation (21), Figure 1a):

Clothing Surface Area Factor (fcl)
As in an earlier study on ambulance clothing [71], for firefighter clothing systems, a wide variety of available equations for fcl calculation were utilized. The best estimations were not the ones derived specifically for a dedicated clothing system based on a layerby-layer approach [71] for a winter clothing system [103] or estimation equations based on a wider database [50], but instead according to the recent estimation for a modern Western clothing database [63] (R 2 = 0.978, Figure 2), possibly indicating similar trends in fashion: 1.01 1.599 .
The specific equation for clothing area factor estimation (R 2 = 0.978; Icl in clo) for this clothing system was 1.02 0.2314 .
Any further tests with other firefighter clothing ensembles may need to prove if the equations above also fit in general with those. Those combinations with the respiratory protective system (SCBA) did not specifically stand out from the general variation, neither for clothing area factor estimation nor for the basic clothing insulation summation method. Figure 1. Separate relationships between measured and summed basic clothing insulation (I cl , clo) for (a) operational uniform and (b) incident clothing ensembles. Icl (sum, ISO 9920) represents summation according to ISO 9920, and Icl (sum, new) modified the summation for these datasets.

Clothing Surface Area Factor (f cl )
As in an earlier study on ambulance clothing [71], for firefighter clothing systems, a wide variety of available equations for f cl calculation were utilized. The best estimations were not the ones derived specifically for a dedicated clothing system based on a layer-bylayer approach [71] for a winter clothing system [103] or estimation equations based on a wider database [50], but instead according to the recent estimation for a modern Western clothing database [63] (R 2 = 0.978, Figure 2), possibly indicating similar trends in fashion: The specific equation for clothing area factor estimation (R 2 = 0.978; I cl in clo) for this clothing system was f cl = 1.02 + 0.2314 × I cl .
Any further tests with other firefighter clothing ensembles may need to prove if the equations above also fit in general with those. Those combinations with the respiratory protective system (SCBA) did not specifically stand out from the general variation, neither for clothing area factor estimation nor for the basic clothing insulation summation method.

The Use of Collected Data
The outcome of the individual items and garment ensembles measurements allows us to evaluate the effects of (a) variation in individual items and their combinations, (b) ways of dressing, and (c) design solutions. Gathering specific basic data on firefighter clothing systems allowed us to adjust the clothing summation method for improved prediction accuracy when selecting and recommending the clothing combinations. The basic data on clothing evaporative resistance give us improved prediction accuracy when selecting and recommending the clothing combinations. In combination with input from Figure 2. Relationship between the measured and the estimated clothing area factors for the firefighter clothing system by various available equations. f cl (or fcl)-clothing area factor; IT-total clothing insulation; Icl-basic clothing insulation.

The Use of Collected Data
The outcome of the individual items and garment ensembles measurements allows us to evaluate the effects of (a) variation in individual items and their combinations, (b) ways of dressing, and (c) design solutions. Gathering specific basic data on firefighter clothing systems allowed us to adjust the clothing summation method for improved prediction accuracy when selecting and recommending the clothing combinations. The basic data on clothing evaporative resistance give us improved prediction accuracy when selecting and recommending the clothing combinations. In combination with input from other relevant research in this area [26,37,42,43,[104][105][106][107], the collected data allow us to model the exposure to a variety of incident scenarios in the presently available clothing systems, and most importantly, validate those available models for rescue services incident scenarios and exposure conditions if the experimental data with the same or similar clothing systems will be available.
As insulation and evaporative resistance values of all individual body parts and body areas are also available (see Appendices A and B), there is a possibility to (a) use these datasets in advanced thermo-physiological models to obtain more detailed predictions; (b) find weak areas in the clothing systems and suggest improvements in clothing design; (c) perform more accurate exposure predictions for improved purchase support. However, experimental data have to be collected for proper model selection and validation. Knowing the clothing properties will greatly assist in decision making in relation to extreme weather events and indoor exposures [23,24]. In the same way, exposure to indoor (office) environments can be estimated according to ISO 7730 [82]. This standard introduces PMV (predicted mean vote) and PPD (predicted percentage of dissatisfied) indexes that relate the environmental parameters to thermal responses of the exposed population. However, the tools are not easily available for much more extreme exposure conditions that the firefighters face.

Conclusions and Future Perspectives
The work of this paper aimed to evaluate the thermal properties of the presently used firefighter clothing systems. Some complementary evaluations should be performed further on, for example: measuring body motion (walking) effects on manikins and the effects of higher air velocity and the combinations of these factors to estimate ventilation effects in the clothing systems. A separate line of research on firefighter clothing systems is the continued information collection and material testing for the effects of heat radiation levels.
Additionally, some common ways of wearing presently used garments should be tested, e.g., front zip open/closed/partly closed, Velcro or rubber bands at wrists or ankles tightened or loose, in order to collect more information on clothing performance for future design solution selection. In addition, the effects and usefulness of available cooling solutions [108] should be evaluated for firefighter use, taking into account logistical challenges. Further, testing of specific design solutions/ideas, including the effects of smart functions, must be considered in order to support the development of standards that are able to cover evaluation and certification needs of the new technical solutions so that they can safely be taken into use by the rescue services [109,110].
As many simple, but also advanced, thermo-physiological models do not support exposure to extreme environmental conditions or highly protective clothing (including impermeable layers), after checking the capacities of existing advanced models [40,44], the needs for improving the prediction models must be settled. Such validation requires considerable experimental work with human subjects. However, the use of such validated models allows us to evaluate clothing design solutions before and under prototyping phases and to diminish the need of later testing on human subjects and minimize costs on changes. At the same time, testing the protective clothing systems for compatibility and ergonomics is still needed [111,112].
The basic data on the present clothing system were collected, and the database was created. This allows us to proceed with: • Selection and recommendations of clothing combinations for the best performance and comfort under specific exposure conditions; • Evaluation of the design aspects of the protective clothing and protective clothing systems; • Modelling and prediction model validation, and the comparison of their deviation from human responses in realistic conditions with human subjects; • Supporting the development of the firefighter's modular protective clothing system of the future; • Developing digital decision support tools for rescue services' purposes; • Creating educational and demonstration materials for optimal protection selection. Appendix A Table A1. Clothing ensembles' local and regional clothing insulation values (I T,i , I T,z,i , m 2 K/W).

Appendix C. Visualisation and Ranking of Measured Insulation and Evaporative Resistance Values
Biology 2022, 11, x FOR PEER REVIEW 24 of 30 Figure A1. Basic clothing insulation of the tested items (Iclu, clo). As the firefighter coverall had over 2 times higher basic insulation, the y-axis scale was not extended to that level, but the value (1.52 clo) is given above the top of this column. Alternating colour shade groups differentiate between product groups and/or body regions covered by these items, e.g., feet (socks, shoes), hands (gloves) head, SCBA element variations, etc.

Figure A1.
Basic clothing insulation of the tested items (I clu , clo). As the firefighter coverall had over 2 times higher basic insulation, the y-axis scale was not extended to that level, but the value (1.52 clo) is given above the top of this column. Alternating colour shade groups differentiate between product groups and/or body regions covered by these items, e.g., feet (socks, shoes), hands (gloves) head, SCBA element variations, etc.   . Measured (Icl (meas)) and summed (Icl (sum, ISO 9920)) basic clothing insulation of all ensembles (Icl, clo) with their relative differences (Diff(%)). For clothing combination descriptions see Tables 1-3. Figure A3. Uncorrected total evaporative resistances (Ret,raw), and according to Wang et al. [98], corrected total (Ret) and clothing (Recl) evaporative resistances of the tested clothing ensembles (m 2 Pa/W). For clothing combination descriptions, see Tables 1-3. Figure A3. Uncorrected total evaporative resistances (Ret,raw), and according to Wang et al. [98], corrected total (Ret) and clothing (Recl) evaporative resistances of the tested clothing ensembles (m 2 Pa/W). For clothing combination descriptions, see Tables 1-3.