A Field Study on Thermal Comfort in Multi-Storey Residential Buildings in the Karst Area of Guilin

: It is important to consider reducing energy use while improving occupants’ indoor thermal comfort. The actual thermal comfort needs and demands should be considered to determine the indoor thermal environment design. In previous studies, research has not been carried out on thermal comfort in karst areas. Thus, a long-term field investigation was carried out on multi-storey residential buildings in the karst area of Guilin city centre during summer (from August 2019 to September 2019) and winter (from December 2019 to January 2020). In this study, the indoor thermal environments of three categories of dwellings were analysed. A total of 77 residential buildings with 144 households were randomly selected, and 223 occupants from 18 to 80 years old participated. A total of 414 effective questionnaires were collected from the subjects. The results show that there was an obvious conflict between the predicted mean vote (PMV) and the thermal sensation vote (TSV). The neutrality temperatures calculated by the regression method were 24.2 ◦ C in summer and 16.2 ◦ C in winter. The thermal comfort range was observed at operative temperatures of 20.9–27.5 ◦ C in summer and 12.2–20.1 ◦ C in winter. The desired thermal sensation for people in the Guilin karst area was not always reflected in the thermal neutrality range. A preference for warmness was identified in the survey.


Introduction
Thermal comfort research has been a hot topic in recent years [1].It is important to consider reducing energy use while improving occupants' indoor thermal comfort.Due to economic growth, there has been a continued and growing demand for the improvement of indoor thermal environments and, consequently, the growth of energy demand for both heating and cooling [2].In residential buildings, varieties of thermal sensation and comfort requirements significantly impact energy consumption [3].There would be a potential waste of energy to maintain the indoor thermal comfort thresholds using the thermal comfort standards directly [4] because of the different thermal comfort characteristics in different outdoor climates and areas [5].The indoor thermal environment should be carefully analysed and controlled with consideration for the actual thermal comfort threshold needs and requirements.
Thermal comfort standardisation and modelling have been widely used.The predicted mean vote (PMV) model is now most commonly used to predict and evaluate environmental thermal comfort [6].The international standards of ISO 7730:2005 [7], ASHRAE 55:2017 [8], and EN16798:2019 [9] adopted the PMV model to evaluate indoor thermal environments.The physiological and physical parameters, including air temperature (T a ), mean radiant temperature (T r ), air velocity (v a ), relative humidity (HR), activity level (met), and clothing insulation (clo), were regarded as the main influential variables for thermal comfort in the PMV model.In this model, the neutral temperature is defined as the optimum condition where habitants do not want a warmer or colder environment.
However, many previous studies have found discrepancies between the PMV and the actual thermal sensation vote (TSV).Habitants have a wider acceptable temperature range in naturally ventilated buildings than that in air-conditioned buildings because of the dynamic response of thermal sensations to the outdoor microclimate [10,11].There is a potential interaction between microclimates and human habitual adaptive behaviour, which adjusts to a comfortable indoor environment according to the occupants' thermal expectations [12].The physical, physiological, and psychological parameters, as well as thermal preference, should be considered in the reaction description of the human body to the variable indoor thermal environment [13].The authors of plenty of field studies have claimed that occupants' thermal comfort changes with block type, the indoor environment, and living habits [14][15][16][17][18][19].Many studies have also found that differences exist in occupants' thermal sensations in different areas and local microclimates [20][21][22][23].Although previous field studies have been conducted on diverse building types in various locations and climates, each study has had unique sample characteristics.
Guilin (109 • 36 50"-111 • 29 30" E, 24 • 15 23"-26 • 23 30" N) is located on areas of alluvial plains between the limestone towers in western China.The tower-shaped karsts of Guilin are the most beautiful landforms in China and the world [24].Along both sides of the Lijiang River in the city centre of Guilin, many limestone towers of its well-known fenglin karst are lined up in an orderly formation and covered with verdant bush.In the city centre of Guilin, the residential buildings are mainly multi-storey houses that are no taller than seven storeys due to a height restriction policy of the Guilin government.Therefore, a unique architectural microclimate environment that potentially influences the occupants' thermal sensations was formed amongst the residential buildings, fenglin karst, plains, and rivers.However, most field studies performed on residential thermal comfort in China have focused on the developed plain cities in southern and central China and there has been no research done on thermal comfort in the karst areas.There are no generally recommended acceptable comfort range and specific thermal comfort prediction models for the existing residential buildings in undeveloped areas in western China, particularly in Guilin City.
In this study, we researched the possible correlations between the occupants' thermal sensations to on-site environmental monitoring and in situ measurements of multi-storey residential buildings under natural ventilation in the Guilin karst area in order to analyse the effect of the fenglin karst and river factors on the thermal comfort results and to develop a rigorous database for designing an optimal indoor thermal environment.

Climate
Our study took place in the Guilin karst area.Guilin is located in a hot summer and cold winter climate (HSCW) zone in China.The Guilin area is characterised by a subtropical monsoon climate, which is dominated by the East Asian monsoon and characterised by two distinct seasons of a cold and rainy winter and a hot and rainy summer.The annual mean air temperature of Guilin City is 18.9 • C. The hottest months are July and August, with an average outdoor temperature of 28.2 • C, while the coldest month is January, the mean temperature of which was 8.1 • C from 1951 to 2020.The average annual rainy days and the average annual rainfall are 172.5 days and 1886 mm, respectively.The rainy season lasts from April to August and accounts for 77.4% of the total rainfall.The maximum mean monthly outdoor relative humidity is 81% in June, while the minimum relative humidity is 66% in December.Figure 1 shows the meteorological data of a typical year in Guilin City.
Sustainability 2021, 13, x FOR PEER REVIEW minimum relative humidity is 66% in December.Figure 1 shows the meteorologica of a typical year in Guilin City.

Location and Description of the Buildings Surveyed
To investigate neutral thermal comfort level and thermal preference, this stud ployed a questionnaire survey, in situ measurements, and environmental monitor 77 residential buildings in 20 residential sub-districts in Guilin city centre.In this s the 20 residential sub-districts were divided into 3 categories: 4 residential sub-di near fenglin karst (NF), 8 residential sub-districts near river (NR), and 8 resid sub-districts neither near fenglin karst nor near river (NN).The dwellings investiga the 3 categories of residential sub-districts were abbreviated to dwelling NF, dw NR, and dwelling NN, respectively.Among the buildings surveyed in the s brick-concrete structure and frame structure accounted for 42.5% and 55.5%, respect Only 3.9% of buildings used external wall insulation; however, buildings sun-shading accounted for about 76.9%.
The surveys were conducted over 61 days during summer (from 1 August 2 30 September 2019) and 42 days during winter (from 20 December 2019 to 31 Ja 2020).In total, 144 households and their dwellings were surveyed, 96 during su (22 for NF, 35 for NR, and 39 for NN) and 48 during winter (12 for NF, 17 for NR, a for NN).Indoor thermal environmental parameters were monitored concurrently the questionnaire surveys were carried out with occupants sampled.One or two pants in a dwelling were chosen.In summer, 145 occupants from 20 to 80 years o participated in the study, 62.8% female and 37.2% male.The distribution of the p pants in summer was as follows: 38 occupants from the dwelling NF (near fenglin k 45 occupants from the dwelling NR (near river), and 62 occupants from the dwellin (neither near fenglin karst nor near river).In winter, 78 occupants from 18 to 80 ye age participated in the study, 61.5% female and 38.5% male.The distribution of th ticipants in winter was as follows: 20 occupants from the dwelling NF, 25 occu from the dwelling NR, and 33 occupants from the dwelling NN.Table 1 shows the acteristics of the participants in our study.The clothing insulation was calculated on the standard ASHRAE 55 and individual survey.
The thermal environment for the main space for daily life (bedroom and room) of each household was chosen as the research environment for both indoor ical parameters and questionnaire.The environment measurement and question survey were in the natural ventilation state.The air conditioning or heating equip was not turned on.

Location and Description of the Buildings Surveyed
To investigate neutral thermal comfort level and thermal preference, this study employed a questionnaire survey, in situ measurements, and environmental monitoring of 77 residential buildings in 20 residential sub-districts in Guilin city centre.In this study, the 20 residential sub-districts were divided into 3 categories: 4 residential sub-districts near fenglin karst (NF), 8 residential sub-districts near river (NR), and 8 residential subdistricts neither near fenglin karst nor near river (NN).The dwellings investigated in the 3 categories of residential sub-districts were abbreviated to dwelling NF, dwelling NR, and dwelling NN, respectively.Among the buildings surveyed in the study, brick-concrete structure and frame structure accounted for 42.5% and 55.5%, respectively.Only 3.9% of buildings used external wall insulation; however, buildings with sun-shading accounted for about 76.9%.
The surveys were conducted over 61 days during summer (from 1 August 2019 to 30 September 2019) and 42 days during winter (from 20 December 2019 to 31 January 2020).In total, 144 households and their dwellings were surveyed, 96 during summer (22 for NF, 35 for NR, and 39 for NN) and 48 during winter (12 for NF, 17 for NR, and 19 for NN).Indoor thermal environmental parameters were monitored concurrently when the questionnaire surveys were carried out with occupants sampled.One or two occupants in a dwelling were chosen.In summer, 145 occupants from 20 to 80 years of age participated in the study, 62.8% female and 37.2% male.The distribution of the participants in summer was as follows: 38 occupants from the dwelling NF (near fenglin karst), 45 occupants from the dwelling NR (near river), and 62 occupants from the dwelling NN (neither near fenglin karst nor near river).In winter, 78 occupants from 18 to 80 years of age participated in the study, 61.5% female and 38.5% male.The distribution of the participants in winter was as follows: 20 occupants from the dwelling NF, 25 occupants from the dwelling NR, and 33 occupants from the dwelling NN.Table 1 shows the characteristics of the participants in our study.The clothing insulation was calculated based on the standard ASHRAE 55 and individual survey.
The thermal environment for the main space for daily life (bedroom and living room) of each household was chosen as the research environment for both indoor physical parameters and questionnaire.The environment measurement and questionnaire survey were in the natural ventilation state.The air conditioning or heating equipment was not turned on.

Environmental Parameters Measurements
In the field study, the outdoor environmental conditions, including the outdoor air temperature (t a ) and relative humidity (RH), were monitored with a thermometer of a HOBO MX2300 series data logger (resolution 0.04 • C, RH 0.05%) to investigate the neutral thermal comfort for the climate of Guilin.The indoor environmental conditions were recorded using a thermometer of an AZ8829S data logger (resolution 0.1 • C, RH 0.1%) and a hot-wire anemometer of ST733.The globe temperature (t g ) was recorded with a 150 mm diameter globe thermometer (resolution 0.1 • C).Table 2 shows the variables monitored and information about the devices used, based on the conditions of instruments for measuring physical quantities in the standard of ISO 7726:1998 [25].
The assessment and analysis processes of the indoor thermal environment were according to the standards of ASHRAE 55 and GB 50785-2012 [26].Data were measured at one central point if the indoor area was less than 16 m 2 .If the indoor area was between 16 m 2 and 30 m 2 , data were measured at the two dividing points that divided the diagonal line of the space into three equal parts.Data were measured at three points when the indoor area was more than 30 m 2 .All instruments were checked, as required, before testing.They were placed far away from heat sources.The height of 1.1 m above the floor was selected at each point [27], and the measurement intervals were 1 min.

Thermal Comfort Surveys
The horizontal and longitudinal surveys were conducted simultaneously in the dwellings NF, NR, and NN.The participants were asked to complete the longitudinal questionnaire at least twice.The interval between two questionnaires should be more than 20 min.Participants were asked to sit indoors for at least 15 min before the survey and then completed questionnaires.A total of 242 effective responses were collected over the summer survey period, 64 responses in the dwelling NF, 78 responses in the dwelling NR, and 100 responses in the dwelling NN, while a total of 172 effective responses were collected over the winter survey period, 40 responses in the dwelling NF, 59 responses in the dwelling NR, and 73 responses in the dwelling NN.
The content of the questionnaire was explained by the surveyor face-to-face.If the participant (such as the elderly) was unable to write, the survey was conducted through an interview.The questionnaire included: (1) Basic information of buildings and participants, e.g., surroundings environment, carpet area, insulation strategies, sun-shading, gender, age, clothing, height and weight, time living in Guilin, annual family income; (2) Thermal subjective sensation vote: questionnaire options with 2 to 7 scales are listed in Table 3.
This study respects the fundamental principles established in the Declaration of Helsinki and was approved by the ethics committee of the Guilin University of Technology.

Thermal Environmental Conditions
Table 4 shows the average, minimum, and maximum values and their standard deviation of the indoor and outdoor variables recorded during the field study.As shown in Table 4, the indoor air temperature (t a,in ) varied between 26.5 • C and 33.3 • C with the mean value of 30.0 ± 1.7 • C during the summer survey and between 10.4 • C and 22.8 • C with the mean value of 16.1 ± 2.6 • C during the winter.The average indoor relative humidity (RH in ) was about 46.0% and 63.6% during summer and winter, respectively.The maximum outdoor air temperature (t a,out ) was 33.3 ± 1.9 • C during summer, while it was about 20.8 ± 3.1 • C during the cold period.
The variations of the indoor thermal environment in the dwellings NF, NR, and NN are shown in Figure 2. In summer, regarding the t a,in and the globe temperature (t g ), similar values were monitored in the dwellings NF, NR, and NN, although the t a,in fluctuation was larger in the dwelling NN.The average RH in in the dwelling NF was the lowest at about 30.5%, which was approximately 20% lower than that in the dwellings NR and NN.In winter, the average t a,in and the average t g in the dwelling NF were the lowest, and those in the dwelling NR were the highest.No significant diversity was observed in the average RH in in the spaces.The most unstable indoor thermal environment occurred in the dwelling NN.

Distribution of Thermal Sensation Votes (TSV)
In the surveys, the thermal sensation votes (TSV) were applied to respond and evaluate how comfortable the subjects found the temperature in their home according to the widely used ASHRAE thermal comfort scale (as shown in Table 3).
The distribution of TSV is shown in Figure 3.The majority of the participants reported a hot thermal sensation accounting for 43.7% in summer, and a neutral thermal sensation accounting for 49.1% in winter during the interviews.It is generally the case that people voting within the central three categories of thermal sensation (−1, 0, and +1) are considered comfortable [28].In the survey, more than 81% of the votes in winter were within the comfort range of TSV (−1, 0, and +1), while only 53% of the votes of comfort were in summer.stainability 2021, 13, x FOR PEER REVIEW that people voting within the central three categories of thermal sensation (−1, 0 are considered comfortable [28].In the survey, more than 81% of the votes in win within the comfort range of TSV (−1, 0, and +1), while only 53% of the votes of were in summer.

Predicted Mean Votes (PMV) and TSV
The mathematical equations for PMV were determined by the principle o body heating balance [6,7] (Equations ( 1) and ( 2)).We computed the PMV value code edited by the statistical programming language VB and validated it by th input-outputs.For the dwellings in the Guilin karst area, the average TSV a values were 1.42 and 1.60, respectively, in summer and −0.50 and −0.31, respec winter, as shown in Figure 4. where M is the metabolic rate (W/m 2 ); W is useful work (W/m 2 ); Pa is partial pr water vapour (mmHg); ta is air temperature (°C); fcl is dressing area coefficie mined by the clothing thermal resistance Icl (fcl = 1 + 0.3Icl); tcl is surface tempe clothing (°C); tmrt is mean radiant temperature (°C); and hc is convective heat tra efficient (W/m 2 °C).The violin plots in Figure 5 show the probability density of the data for b and TSV for the range of possible values on the vote scale (from −3 to +3) in th ings NF, NR, and NN.In summer, percentages of the TSV in the interval [−1, + dwellings NF, NR, and NN were 35%, 35%, and 29%, respectively.In winter, t centages were 87%, 90%, and 74%, respectively.The classical PMV mo der-predicted the true comfort of participants due to the significant differences

Predicted Mean Votes (PMV) and TSV
The mathematical equations for PMV were determined by the principle of human body heating balance [6,7] (Equations ( 1) and ( 2)).We computed the PMV value with the code edited by the statistical programming language VB and validated it by the known input-outputs.For the dwellings in the Guilin karst area, the average TSV and PMV values were 1.42 and 1.60, respectively, in summer and −0.50 and −0.31, respectively, in winter, as shown in Figure 4.
where L is the body thermal load (W/m 2 ).
where M is the metabolic rate (W/m 2 ); W is useful work (W/m 2 ); P a is partial pressure of water vapour (mmHg); t a is air temperature ( • C); f cl is dressing area coefficient determined by the clothing thermal resistance I cl (f cl = 1 + 0.3I cl ); t cl is surface temperature of clothing ( • C); t mrt is mean radiant temperature ( • C); and h c is convective heat transfer coefficient (W/m 2 • C).The violin plots in Figure 5 show the probability density of the data for both PMV and TSV for the range of possible values on the vote scale (from −3 to +3) in the dwellings NF, NR, and NN.In summer, percentages of the TSV in the interval [−1, +1] in the dwellings NF, NR, and NN were 35%, 35%, and 29%, respectively.In winter, their percentages were 87%, 90%, and 74%, respectively.The classical PMV model under-predicted the true comfort of participants due to the significant differences between TSV and PMV.Additionally, the TSV showed a much narrower range.

Regression Method and Comparison between TSV and PMV
As is confirmed by references [29,30], the mean radiant temperature (tmrt) and air temperature (ta) influence the uncertainty in thermal comfort evaluation.It non-uniform that the radiant temperature distributed on the globe thermometers [31 order to reduce the uncertainty, the operative temperature (top), which is a synth temperature comprehensively being considered as the effects of ta and tmrt on the the

Regression Method and Comparison between TSV and PMV
As is confirmed by references [29,30], the mean radiant temperature (tmrt) and the air temperature (ta) influence the uncertainty in thermal comfort evaluation.It was non-uniform that the radiant temperature distributed on the globe thermometers [31].In order to reduce the uncertainty, the operative temperature (top), which is a synthetic temperature comprehensively being considered as the effects of ta and tmrt on the thermal perception of the human body, was recommended for use [29].The top is calculated as in Equation ( 3) [32]: where top is operative temperature (°C); ta is air temperature (°C); tmrt is mean radiant temperature (°C); and A is the average weight factor of the ta and tmrt.

. Regression Method and Comparison between TSV and PMV
As is confirmed by references [29,30], the mean radiant temperature (t mrt ) and the air temperature (t a ) influence the uncertainty in thermal comfort evaluation.It was nonuniform that the radiant temperature distributed on the globe thermometers [31].In order to reduce the uncertainty, the operative temperature (t op ), which is a synthetic temperature comprehensively being considered as the effects of t a and t mrt on the thermal perception of the human body, was recommended for use [29].The t op is calculated as in Equation ( 3) [32]: where t op is operative temperature ( • C); t a is air temperature ( • C); t mrt is mean radiant temperature ( • C); and A is the average weight factor of the t a and t mrt .
Under the natural convection condition, the t mrt was calculated as in Equation ( 4) [33]: where t mrt is mean radiant temperature ( • C); t g is globe temperature ( • C); and t a is air temperature ( • C).
To compare the TSV and the PMV, the whole range of t op was used to analyse the variation of the data.A scatter diagram of t op and TSV and PMV for dwellings surveyed in the Guilin karst area was planned and equations of linear regression were established, as shown in Figure 6.The intervals of 1 • C were defined for the t op .In summer, the linear regression equation for the mean TSV was TSV = 0.2576t op − 6.2295 (R 2 = 0.8911) and for the PMV was PMV = 0.3929t op − 10.106 (R 2 = 0.8292), as shown in Figure 6a.Additionally, as t op increased at the interval of t op < 28.7 • C, the difference between PMV and TSV decreased gradually.While t op increased at the interval of t op > 28.7 • C, the difference between PMV and TSV increased.In winter, the linear regression equations were TSV = 0.2143t op − 3.4618 (R 2 = 0.8858) and PMV = 0.201t op − 3.496 (R 2 = 0.5202) for the TSV and PMV, respectively, as shown in Figure 6b.The two regression lines of TSV and PMV were almost parallel.
To further investigate the effect of the fenglin karst and river factors on the results, the data were analysed independently for TSV and PMV in the dwellings NF, NR, and NN, respectively, as shown in Figure 7.In summer, six regression lines were found to intersect at around the point value of t op = 29 • C. At the interval of t op < 29 • C, the three regression lines of TSV in dwellings NF, NR, and NN were higher than the PMV, while at the interval of t op > 29 • C, the TSV was overestimated in the PMV values.In winter, the TSV always was underestimated in the PMV values.Similarly, Yang [34] found that the PMV for the elderly was lower than the TSV in winter and was higher than the TSV in summer under a naturally ventilated condition.Alike results were also confirmed in surveys for naturally ventilated residential and office buildings by Dhaka [35].
The regression and Griffiths' methods are the two main methods to identify thermal comfort temperature.In this study, occupants' acceptance of the indoor thermal environment and their neutral estimates were calculated by the linear regression method, as shown in Figures 6 and 7.All the regression equations passed the goodness of fit (R2 > 0.5).When the regression equation equals zero, the calculated temperature can be considered as the neutral temperature (t n ) of participants.Their accepted temperature (t ac ) range of 80% is −0.85 < TSV < +0.85 [36].
Table 5 shows the actual and predicted t n and t ac in the dwellings NF, NR, and NN during summer and winter.The actual t n was lower than the predicted t n in the dwellings NF, NR, and NN, respectively.Yao [37] and Yu [38] proposed similar results in different buildings and different climates by their previous research.In summer, there was no significant difference of the t n between the dwellings NF, NR, and NN.In winter, the t n in the dwelling NF was the lowest in all of the dwellings.The range of actual t ac (80%) calculated by the regression method in summer and winter was from 20.9 to 27.5 • C and from 12.2 to 20.1 • C, respectively, as shown in Table 5.

Preferred Temperature
The occupants' preferred temperature (t p ) was explored to be compared with their actual t n .For the purpose, the "probit" analysis was conducted for the thermal preferred votes (TPV) and their corresponding t op .Probit regression analysis was first presented by Webb [39].It is generally applied to illustrate the probability of an objective for data with binary response variables.Based on the investigation, the overall TPV were divided into two parts as preferences for cooler (assigned "−1") and preferences for warmer (assigned "+1") surroundings.The result of cumulation of TPV was set as "1" according to the Probit analysis.The t p was identified as the intersection point of the curve with "cooler" preference and the curve with "warmer" preference [40], as shown in Figure 8.
The t p in summer and winter were 25.8 • C and 18.6 • C, respectively.The result displayed an obvious difference between the t p and the actual t n .In accordance with the analysis, t p was around 1.6 • C higher than the actual t n in summer and 2.4 • C in winter.The analysis result clarified that occupants have a tendency towards warmer indoor environments for their t p in Guilin.

Preferred Temperature
The occupants' preferred temperature (tp) was explored to be compared with their actual tn.For the purpose, the "probit" analysis was conducted for the thermal preferred votes (TPV) and their corresponding top.Probit regression analysis was first presented by Webb [39].It is generally applied to illustrate the probability of an objective for data with binary response variables.Based on the investigation, the overall TPV were divided into two parts as preferences for cooler (assigned "−1") and preferences for warmer (assigned "+1") surroundings.The result of cumulation of TPV was set as "1" according to the Probit analysis.The tp was identified as the intersection point of the curve with "cooler" preference and the curve with "warmer" preference [40], as shown in Figure 8.
The tp in summer and winter were 25.8 °C and 18.6 °C, respectively.The result displayed an obvious difference between the tp and the actual tn.In accordance with the analysis, tp was around 1.6 °C higher than the actual tn in summer and 2.4 °C in winter.The analysis result clarified that occupants have a tendency towards warmer indoor environments for their tp in Guilin.

Thermal Comfort Characteristics Comparison with Previous Studies
Table 6 summarises the discrepancy of tn, tac (80%), and TSV equations in the previous field studies [4,[41][42][43][44][45][46][47][48][49] which were conducted in naturally ventilated buildings of various climates and locations in China.As the results of comparison among the SC/C (severe cold/cold), HSCW (hot summer and cold winter), and HSWW (hot summer and warm winter) zones, occupants' tn and the maximum of tac generally increased when the outdoor atmosphere became warmer.In the HSCW zone, both tn and tac were higher in

Thermal Comfort Characteristics Comparison with Previous Studies
Table 6 summarises the discrepancy of t n , t ac (80%), and TSV equations in the previous field studies [4,[41][42][43][44][45][46][47][48][49] which were conducted in naturally ventilated buildings of various climates and locations in China.As the results of comparison among the SC/C (severe cold/cold), HSCW (hot summer and cold winter), and HSWW (hot summer and warm winter) zones, occupants' t n and the maximum of t ac generally increased when the outdoor atmosphere became warmer.In the HSCW zone, both t n and t ac were higher in other cities than those in the Guilin karst area in summer due to the potential impact of the fenglin karst and the unique microclimate of the surroundings in Guilin.The results also presented that the occupants' t ac ranges of indoor thermal environment under natural convection in a HSCW climate zone were different.Hence, the design t ac range of different locations in the HSCW zone should be varied, which was in contradiction with the Chinese standard [26].

Conclusions
This paper analyses the data of a thermal comfort field study conducted in three types of naturally ventilated dwellings NF, NR, and NN in the Guilin karst area, China, during summer (61 days) and winter (42 days).The results are derived from environmental measurements, indoor and outdoor, taken in 144 dwellings and thermal comfort questionnaires collected for 414 responses of 223 participants, which could be extrapolated to karst areas in the hot summer and cold winter zone of China.The main conclusions are follows: 1.
In summer, an indoor thermal environment in dwelling NF was on average 1 • C cooler compared with that in dwelling NR and 0.5 • C cooler than that in dwelling NN.In winter, average indoor temperature in dwelling NR was 2.7 • C warmer than that in dwelling NF and 1 • C warmer than that in dwelling NN.During both summer and winter, the fluctuations of indoor temperature and relative humidity in dwelling NN were the largest; 2.
According to the TSVs, most occupants at about 81% were comfortable in winter, however, only 53% of the votes were comfortable in summer.The average TSV and PMV values were 1.42 and 1.60, respectively, in summer and −0.50 and −0.31, respectively, in winter.The actual thermal comfort of occupants could not be precisely predicted by the classical PMV model, in most cases, as shown by the significant differences between TSV and PMV; 3.
Although the most desired sensation was "neutral", the desired temperature was not always the thermal neutral temperature for occupants.As seen in the results of this study on thermal comfort, a negative relationship between thermal sensation votes and thermal preference votes suggested that occupants preferred warmness in the Guilin karst area of the HSCW zone during both summer and winter; 4.
The range of actual accepted temperature (80%) in summer and in winter was 20.9-27.5

Figure 1 .
Figure 1.The meteorological data of typical year in Guilin City.

Figure 1 .
Figure 1.The meteorological data of typical year in Guilin City.

Figure 2 .
Figure 2. Variations of indoor thermal environment in the dwellings NF, NR, and NN.(a) Summer; (b) Winter.

Figure 3 .
Figure 3.The percentage of occupants' thermal sensation votes in Guilin karst area.

Figure 3 .
Figure 3.The percentage of occupants' thermal sensation votes in Guilin karst area.

Figure 5 .
Figure 5. Density plots for TSV and PMV of the dwellings NF, NR, and NN.(a) Summer; (b) Winter.

Figure 6 .
Figure 6.Relationship between participants' TSV/PMV and operative temperature in the survey dwellings in Guilin karst area.(a) Summer; (b) Winter.

Table 1 .
The characteristics of the participants in our study.

Table 2 .
Monitoring instruments used in the survey.

Table 4 .
Thermal physical variables during the field study.Note: t a,in -Indoor air temperature, RH in -Indoor relative humidity, v a,in -Indoor air speed, t g -Globe temperature, t op -Operative temperature, t a,out -Outdoor air temperature, RH out -Outdoor relative humidity, S.D.-Standard deviation.

Table 4 .
Thermal physical variables during the field study.Indoor air temperature, RHin-Indoor relative humidity, va,in-Indoor air speed, tg-Globe temperature, top-Operative temperature, ta,out-Outdoor air temperature, RHout-Outdoor relative humidity, S.D.-Standard deviation.

Table 5 .
Neutral temperatures calculated by the regression method.

Table
Thermal comfort characteristics in different locations.
• C and 12.2-20.1 • C, respectively.The actual thermal T n of 24.2 • C and 16.2 • C in summer and in winter, respectively, was lower than the predicted thermal T n of 25.7 • C and 17.4 • C in summer and in winter, respectively.