Introduction

The tourism industry is one of the largest industries in the world and has always affected the economies of many tourist-friendly countries (Jover and Díaz-Parra 2020). The industry can improve relations between countries, cultural exchanges (Yu 2021), and various occupations and help countries to generate more revenue by creating opportunities for entertainment and recreation for tourists (Alipour et al. 2012; Matthew et al. 2021). The latest published statistics show that after France, with an annual attraction of about 89 million foreign tourists, Spain has hosted 82 million foreign tourists annually, which this volume of foreign tourists has set a clear vision for the tourism industry of this country (Argandoña 2010; Manganelli et al. 2016). However, industrial tourism is highly dependent on climate and its changes (Saarinen and Tervo 2006; Wall and Badke 1994). So that, one of the most crucial information needed by tourists is the weather conditions of tourist destinations (Hopkins 2014; Saarinen and Tervo 2006) and proper thermal comfort as the desirability and acceptability of the climate play an important role in deciding for tourist destinations (Lin and Matzarakis 2008). Indeed, it is an undeniable fact that urban areas are growing on a global scale (Afsharzadeh et al. 2021; He 2018), which has affected climate change (increasing levels of hardship) (Gachkar et al. 2021; Mohammad et al. 2021). These led to disruption in urban environments (increased heat waves and heat stress conditions) and the presence of tourists in cities (Bigano et al. 2008; Lopes et al. 2021; Mohammad and Goswami 2019).

Numerous studies have been conducted on thermal comfort in urban microclimates, all of which point to a direct relationship between climatic conditions and the use of outdoor environments (Chen et al. 2020; Deng and Wong 2020; Karimi et al. 2021; Knaus and Haase 2020; Mohammad et al. 2019; Mohammad and Goswami 2021a, 2021b). Thermal comfort in open metropolitan spaces is of extensive significance because of its huge effect on human wellbeing (Falah et al. 2020; Xiong et al. 2015), urban vitality (Afsharzadeh et al. 2021) and the energy consumption of buildings (Ibrahim et al. 2021). It is better to identify the indexes that influence thermal comfort to create the proper conditions (Gachkar et al. 2021; Ng and Cheng 2012). Various factors are influencing the promotion of thermal situation in open urban spaces, which can be mentioned to plants (Shafiee et al. 2020; Simon 2016), the albedo of material (Gachkar et al. 2021; Mohammad et al. 2019, 2021; Zhou et al. 2018), the urban form including direction (Abaas 2020; Athamena 2022; Karimimoshaver and Shahrak 2022), sky view factor (SVF) (Jamei and Rajagopalan 2017; Song and Jeong 2016), and aspect ratio (Vassiliades et al. 2022; Xiaomin et al. 2006), which affects small urban microclimates and the role of water (Chan et al. 2017). For example, the role of the percentage of urban trees on thermal comfort in hot and humid climates of Roorkee taluk in Uttarakhand state of India has been investigated by (Mohammad et al. 2021). Their result shows that when the urban canyon is cover with 50% of the proper trees, the air temperature decreases by 0.97 °C during the day and by 0.71 °C during the night. Another study by Fahed et al. (2020) reveals that the ambient temperature of the roof top decreases with an increasing albedo of the material and reaches a maximum value of 2.5 °C during mid-noon. Also, Ali-Toudert and Mayer (2006) investigated the effect of various orientations of urban forms in summer and winter seasons in a tropical city of Ghardaia and found that the N-S orientation has less exposure to solar radiation in winter than in summer and leads to making better thermal comfort. Another research is conducted by Gachkar et al. (2021) regarding the effects of water on thermal comfort in Urmia and the results illustrated that a proper arrangement of the water pond creates a significant cooling effect for the thermal comfort of pedestrians.

Understanding the relationship between outdoor thermal comfort impact and tourism has been one of the most important research indicators among scientists since the late eighteenth century. So that, Mieczkowski provided a tourism climate comfort index (TCCI) to assess tourist climate relations (Mieczkowski 1985). A few years later, in 2003, Freitas estimated the tourism climate in light of three essential perspectives of aesthetic, physical, and thermal (de Freitas 2003). In connection with the relationship between climate parameters and tourism, various research have been done. Nasrollahi et al. (2017) assessed several climatic variables in the hot and dry district of Iran examined and its relationship with the thermal comfort of tourists and finally set on the appropriate time for visitors in this historical city of Isfahan. In a similar research, Matzarakis et al. (2013) conducted studies on meteorological indicators and their impact on tourism. Zare et al. (2018) have used ISO 7730 to analyze thermal comfort in Kerman based on both subjective and objective procedures. As well, different research have been reported in other parts of planet earth, including the city of Taiwan (Lin and Matzarakis 2008), Italy (Salata et al. 2017), China (Ma et al. 2019), Canada (Scott et al. 2007), and Portugal (Lopes et al. 2011). A survey of studies led on various environments showed that the majority of these studies have been done on an enormous scale and utilizing the information given by different meteorological stations to assess the significant pointers and decide the scope of thermal comfort reasonably for sightseers. However, unlike other previous studies, which focuses the impact of thermal comfort on historical sites, this study considered the range of thermal comfort and comparing them in different urban plazas of Seville and Madrid. This study tries to identify the most appropriate time for tourists to visit each of the historical sites of Seville and Madrid during the day and night.

Method

Selected area

The aim of this study is to investigate the thermal comfort conditions of various cultural sites that were compared during the day and night hours to distinguish the reasonable time to visit each cultural site in Sevilla and Madrid.

Sevilla is situated in the southwest of the Iberian Peninsula (37° 21′ 55″ N, 6° 0′ 30″ W) and as the capital of the Independent People Community of Andalusia (Castaño-Rosa et al. 2018) has had encountered a huge development in recent 20 years (Herrera-Gomez et al. 2017; Jover and Díaz-Parra 2020). It has a total urban area of about 140 km2 with a population of around 700,000 inhabitants and a population density of 5000 citizens/km2. The city experiences Mediterranean climate (Csa) as per Köppen-Geiger climatic classification (Kottek et al. 2006), with mild winters and blazing hot summers. The warmest month of the year was recorded in July in the city of Seville, with maximum and minimum temperatures of 34 °C and 18.9 °C. On the other hand, January is considered as the coolest month with average maximum temperatures of 16 °C and average minimum temperatures of 5.7 °C (Galán et al. 2003), and due to the short distance from the Atlantic Ocean (60 km), sometimes the temperature is changed by the gentle breeze of the sea (Rodrigo 2019). It receives an average total annual rainfall of about 554 mm, mainly within the month of October to April.

The city of Madrid as one of the tourism hubs is also examined in this study. Madrid (40° 25′ 0.3900′′ N and 3° 42′ 13.6440′′ W), the capital of Spain, with a population of over 16.642 million people and 604 square kilometers area is one of the most important tourist cities in Europe (Gómez-Villarino et al. 2021), like other cities (Barcelona and Sevilla), and is belong to CSA zone (Soutullo et al. 2020). CSA (climate-smart agriculture) is an approach related to agricultural systems and its climate change consequence that help people to increase their productivity and incomes with minimum greenhouse gas emissions wherever possible. Rainfall is concentrated in spring and autumn, while cloudless days are more common in summer. Regionally, its climate is slightly influenced by the presence of the central mountain range, which divides the inner plateau of the Iberian Peninsula into two parts and affects the wind direction (mainly NE-SW). The condition of the Manzanares River, which flows from north to south of the city, further contributes to this impact, directing cold air from the mountains into the city (Garcla and Gómez 1996).

The study areas are among the hot spots in Europe that have high summer temperature, making the area prone to heat islands and sometimes heat wave and heat stress conditions (Follos et al. 2020). High humidity and low wind speed conditions are other features of these cities (Peña et al. 2018), and when such conditions match with high temperatures, the recurrence of fatalities because of hotness stress increments drastically (Kalkstein 1991). This in turn highlights the significance of understanding the consequences of planning decisions and the topic of mitigating risk through developed urban planning. Research conducted in each of the surveyed cities has covered three of the most important historical plazas. One of the most important factors in selecting these urban plazas was a distance between the selected places because the probability of visiting different ancient monuments and places with short distances at a certain time increases. A brief description of each of these urban plazas (three in Sevilla (Fig. 1b) and three in Madrid (Fig. 1c)) is given below:

Fig. 1
figure 1

a Location of Sevilla and Madrid in Spain, b locations of the three plaza’s places in Sevilla, and c locations of the three plaza’s places in Madrid

Plaza de España: This urban plaza is the most admirable building in Seville, Spain, after the cathedral of this city, and the reason is the large and extreme scale along with the grandeur and beauty that can be seen in this structure. It is a landmark example of Regionalism Architecture, mixing elements of the Baroque Revival, Renaissance Revival, and Moorish Revival (Neo-Mudéjar) styles of Spanish architecture (Nsang 2012).

Alameda de Hercules: This is the most established public garden in Europe since the 1992 World’s Fair and has since formed into one of the most well-known social spots for local people and vacationers (Muñoz 1998).

Plaza Nueva: This is a public square in the city center of Seville which was established in 1856. The land in this part of the San Fernando Monastery, which was previously built from 1270 to 1840, was later taken over by the local government and turned into a public square (Navarro De Pablos 2017).

Plaza de Colón: This famous square in the city center of Madrid, dedicated to the man who discovered America. It is one of most important open spaces and a strategic location in the city, leading out to prominent streets such as El Paseo de Recoletos, El Paseo de la Castellana, Calle Génova, Calle Serrano, and Calle Goya. It is also home to eminent buildings such as the Torres de Colón, a monument to Christopher Columbus presides over the square from the center of the fountain around which the traffic flows. The square is also home to the Jardines del Descubrimiento, a park that first opened back in 1970, where you can admire various sculptures as well as a monument dedicated to the discovery of America. Underground to this plaza, Fernán Gómez Cultural Centre is situated.

Plaza Mayor: Plaza Mayor is a square with arched passages in the city of Madrid, which is equally popular with tourists and locals, and its symmetrical rectangular buildings have a uniform architectural system. In medieval times, the landscape of the plaza mayor was the only space for a market outside the city walls. In the 1560s, King Philip II asked Juan de Herrera, the architect of Escorial, to turn the market into a real square, and this lasted until 1617 and the reign of King Philip III. The result was a large field measuring about 120 by 90 m with wooden buildings around it. All the buildings around the square were completely burnt down by fire three times in 1631, 1672, and 1790, and this is the last reconstruction we are witnessing today (Basterra et al. 2009).

Plaza de Santa Ana: Santa Ana Square is one of the sights of Spain, located in the center of Madrid. Santa Ana Square dates to the seventeenth century, but its appearance was very different from the current shape of the square. Construction of the current shape of Santa Ana Square was completed in 1880. The square is famous for its monument to Pedro Calderón de la Barca, the famous writer, and Federico García Lorca, the Granadian poet (Del Corral 2008).

Research period

This study was conducted to find the most appropriate time for tourists to visit each urban plaza along with the importance of their thermal comfort in hot a summer month (Karimi et al. 2020; Mohammad et al. 2021). Meteorological data of Seville and Madrid station were collected for the last ten years (from 2010 to 2020). The average monthly air temperature (Ta) and wind speed (Ws) over Sevilla (Fig. 2a) and Madrid (Fig. 2b) were plotted to select the suitable hottest month for our simulation. From this, we selected July as the suitable month for simulation as it is showing maximum air temperature and minimum wind speed over the last 10 year. Also, the Barajas and Sevilla airport station’s meteorological data near the study areas, situated at 12–14 km and 17–20 km, respectively, was used as climatic indicators during the last 10 years (2010 − 2020).

Fig. 2
figure 2

Dispensation of change of monthly meteorological data averaged from 2010 to 2020 (for selecting suitable month), for Sevilla (a) and Madrid (b); dispensation in July 2020 (for choosing the suitable day for simulation) for Sevilla (c) and Madrid (d)

For finding the suitable day in July month for the simulation, we compared the meteorological data in July 2021 and find that 13th–14th July 2021 in Sevilla (Fig. 2c) and 10th–11th July 2021 in Madrid (Fig. 2d) would be possible day for evaluating the tourist’s thermal comfort in the hottest day of the year with minimum wind speed.

Thermal comfort index

Based on the assessments from previous literatures, it is seen that more than 165 thermal comfort indicators have been designed; however, four common indicators of thermal comfort, including standard effective temperature (SET*), physiological equivalent temperature (PET), predicted mean vote (PMV), and universal thermal climate index (UTCI), have been more commonly used (Mohammad et al. 2021). Among these four indexes, PET is more frequently used as it is produced by taking into account the influence of short- and long-wave radiation in outdoor environments on human energy levels (Matzarakis et al. 2007, 2021). It is widely used in urban textures with intricate shadow patterns (Middel et al. 2017; Mohammad et al. 2021). Moreover, as seen from the many previous researches that have been done in hot and humid climates (Binarti et al. 2021; Hegazy and Qurnfulah, 2020; Lian et al. 2020), it is supposed to use the PET index in the present study with hot and humid climatic condition. It is easily calculated by Ray-Man software advanced at the Meteorological Institute of the University of Freiburg based on the instructions of the German Engineering Association (Gachkar et al. 2021).

Simulation

In this research, like most previous studies of outdoor thermal comfort (Gachkar et al. 2021; Karimi et al. 2020; Mohammad et al. 2021; Nasrollahi et al. 2017), Envi-met was used as a reliable software to investigate climatic conditions in open urban spaces (Wai et al. 2020; Xiong et al. 2020). This software is a non-hydrostatic three-dimensional model designed according to the fundamental of thermodynamics, fluid dynamics, and the laws of atmospheric physics and is able to calculate all influential indexes of outdoor thermal comfort (Forouzandeh 2021; Liu et al. 2021; Yilmaz et al. 2021). Six important historical plazas of Sevilla and Madrid along with their urban areas are simulated in this present study (Fig. 3). Since none of them had a time limit for visitors, the thermal comfort conditions of historic sites were evaluated over a 24-h period for finding the best suitable time for vising tourist. Also, given that in the first 2 h, the simulations lack the necessary credits, so the simulations have been started before the hottest day of July, consider for simulation. In addition, the default values of all other indexes in the Envi-met model were considered in the simulations. The details of the simulations for the Sevilla and Madrid are given in Table 1 and Table 2, respectively.

Fig. 3
figure 3

The location of receptors in each case study

Table 1 ENVI-met input data for simulations in each case study in Sevilla
Table 2 ENVI-met input data for simulations in each case study in Madrid

To calculate the physiological equivalent temperature (PET), Rayman software that made by Matzarakis et al. 2007 has been used and needs different indexes included of air temperature, humidity, wind speed, mean radiant temperature, human activity, and type of cover. Heat transfer resistance of 0.5 Clo (suitable for summer dressing habits) and 80 W activity (walking outdoors) was considered (Taleghani and Berardi 2018; Wang et al. 2016). In addition, as shown in Table 3, the default values of the other model parameters are used in the simulations.

Table 3 RayMan input data for calculating thermal comfort

Field study

Measurement

To accredit the simulation results performed by Envi-met software, field studies were performed with the metrological device on the same day of measurement at a suitable height from the ground as seen in priors studies (Nasrollahi et al. 2017, 2021). The thermal sensing votes (TSV) were examined on a 7-point semantic differential scale (cold, cool, slightly cool, neutral, slightly warm, warm, hot) and, to satisfy tourists, along with questionnaires. Different climatic variables collected during the field measurement includes air temperature (Ta), relative humidity (Rh), wind speed (Ws), and global temperature (TGlobe). This collected parameter helps in the calculation of mean radiant temperature (MRT). Certainly, based on different studies (Chen et al. 2014; Gál and Kántor 2020; Lee et al. 2019), MRT is obtained from Eq. 1, where D is global diameter (Tg) and ε is global emissivity. The obtained results are compared with the simulation outputs, which determine the amount of thermal comfort of each classical plazas of Seville and Madrid. Table 4 illustrates the features of the instruments used in this research.

Table 4 Features of the instrument used in this research
$${T}_{MRT}= {\left[{\left({T}_{GLOBE}+273.15\right)}^{4}+\frac{1.1*{10}^{8}*W{s}^{0.6}}{{\varepsilon *D}^{0.4}}*({T}_{GLOBE}-{T}_{A})\right]}^{0.25}-273.15$$
(1)

Questionnaire

This questionnaire has been prepared in accordance with previous studies (Mohammad et al. 2021; Sharmin et al. 2015) and ASHRAE-55 (Heating et al. 2004). It is mainly divided into two parts. The first part is related to the individual characteristics of tourists such as gender, height, their resident, weight, and age of people, whereas the second part is expressed to the feeling of warmth, thermal preference, thermal satisfaction, and length of stay in plazas. The questionnaires were distributed and completed among tourists in the morning (8 a.m. to 1 p.m.) and in the afternoon (4 p.m. to 9 p.m.) time, with maintaining the conditions of health protocols arising due to COVID-19 situation. Each questionnaire took between 5 and 7 min to complete and was designed for both visitors in English and Spanish. An example of a questionnaire is shown in Fig. 4.

Fig. 4
figure 4

Sample questionnaire used in the research

Validation

In this research, to evaluate the efficiency of ENVI-met software, the criteria coefficient of determination (R2) and the Willmott agreement index (d) have been used like other research (Gachkar et al. 2021; Salata et al. 2016). The Willmott agreement index (d) is calculated using Eq. (2), where Pi and Oi are the predicted and observed values, respectively, n is the total number of observations, and \(\overline{{Q}_{i}}\) is the mean of the observed values. The performance evaluation of three parameters, namely Ta, TMRT, and Ws, which have a significant effect on heat perception of individuals, has been considered in this study. Table 5 reports the mean values of R2 and d using the evaluation performed for each case study.

Table 5  Quantitative assessment of the ENVI-met efficiency based on the evaluation of R2 and d in Sevilla and Madrid
$$d=1-\frac{{\sum }_{i-1}^{n}{\left({P}_{i}-{Q}_{i}\right)}^{2}}{\sum {\left(\left|{P}_{i}-\overline{{Q}_{i}}\right|\right)+\left(\left|{P}_{i}-\overline{{Q}_{i}}\right|\right)}^{2}}$$
(2)

According to conducted results in the plazas of Seville, the correlation ratio (R2) for Ta was between 0.82 and 0.89, MRT was between 0.79 and 0.84, and Ws was between 0.77 and 0.81, while d varies from 0.67 to 0.81 for the different plazas of Seville. As well, in the Madrid case studies good correlation were observed in the different plazas. The correlation ratio (R2) of Ta, MRT, and Ws was between 0.8 and 0.9 in the three selected plazas, whereas d is from 0.68 to 0.79. The gained results is comparable with the outcome of the prior reported studies (Gachkar et al. 2021; Karimi et al. 2020; Lee et al. 2016). It can be accounted for obtaining results are advised satisfactory, and the extent of this difference can be related to elements like limited simulation accuracy in the ENVI-met model (Acero and Herranz-Pascual 2015), incompatibility of native plant profiles with software plants, and variation in MRT (Krüger et al. 2011).

Results

Questionnaire results

Questionnaires were distributed in different sections as shown in Fig. 4 in the selected plazas of Seville and Madrid to assess the feeling of thermal comfort. An overall of 180 questionnaires were collected in 7 days in two shifts in the morning and evening, and each questionnaire lasted about 5 to 7 min to ensure the stability of outside weather conditions. The highest number of forms was entered in Plaza de Espana (92) and the lowest in Plaza Nueva (42). While, in Madrid the completed questionnaires were the lowest and highest in Plaza de Santa Ana (36) and Plaza de Colón (87), respectively. Overall, participants completing the questionnaire in the city of Seville were 55.36% male and 44.64% female. While in the plazas of Madrid, women (58.35%) participated more than men (41.65%) in completing the questionnaire. Domestic tourists were 9.06% from Seville, 8.33% from Malaga, and 7.97% from Cadiz. Most foreign tourists included 5.07% from France, 2.75% from Switzerland, 3.5% from Morocco, 8.05% from Netherland, and 4.5% from Italy. While polls in plazas of Madrid showed that domestic tourists were less than foreign tourists and more than 80% of those polled were from France, Switzerland, Italy, Netherlands, Morocco, and Turkey. In terms of age, the highest frequency is related to tourists between 41 and 50 years and then to the age range of 31 to 40 years in the city of Seville, while the highest frequency in Madrid was related to tourists between 51 to 60 and 41 to 50 years. About 67.86% and 54.86% of tourists in Seville and Madrid surveyed the weather conditions to determine their travel destination, and about 60.85% and 42.65% of them were informed of the reasonable weather situation to visit the ancient sites of Seville and Madrid, respectively.

Thermal comfort sensation of tourists

Given that the questionnaires were evaluated and completed in one of the warm months of the year, the thermal comfort sensation of travelers was more disposed to the hotness range. The explanation of the thermal comfort sensation of the interviewees is illustrated in Fig. 5. With increasing attendance time at each of the assessed plazas in Madrid and Seville, the increase in the feeling of warmth led to a decrease in the level of comfort of tourists and their dissatisfaction with the thermal situation and their priority for cooler conditions developed. While when the visitors are in shady places, their thermal sensation has been increased. In addition, it is necessary to mention that there is no relation between tourist’s emotional states (happiness and sadness) and their thermal comfort sensation. The importance of the relationship between tourist’s sense of thermal comfort and various factors is shown in Tables 6 and 7 for Sevilla and Madrid, respectively.

Fig. 5
figure 5

Percentage distributions of thermal sensation votes (TSV) in all the six sites of Sevilla and Madrid

Table 6  Gamma correlation for thermal sensation and other parameters in the plazas of the city of Sevilla
Table 7  Gamma correlation for thermal sensation and other parameters in the plazas of the city of Madrid

Thermal comfort satisfaction of tourists

Based on the assessment of the results of the questionnaire conducted in the city of Seville, about 38.11% was satisfied, 21.09% were slightly satisfied, and 9.09% were very satisfied with the thermal conditions as most of whom were elderly. While surveys conducted in Madrid showed that about 28.19% were satisfied, 21.09% were slightly satisfied, and 9.31% were very satisfied with the thermal conditions. This is due to high heat stress, and tourists’ satisfaction may be in line with their expectations of this hot season. Figure 6 shows the percentages of satisfaction of tourists in each of the selected plazas of Seville and Madrid.

Fig. 6
figure 6

Percentage distributions of satisfaction in the plazas of the city of Sevilla and Madrid

Percentage distributions of satisfaction at each plaza of Sevilla and Madrid

Evaluations in the city of Seville show that visitors to the Plaza Nueva and Plaza de Espana were more satisfied than Alameda de Hercules (Fig. 7). Plaza Nueva has an almost smaller outdoor space than other squares, and, the shaded trees, also the semi-open porch of Plaza de Espana, got the highest level of tourist satisfaction. However, Alameda de Hercules got the least unlevel of thermal satisfaction among every one of the plazas regardless of having many shade trees. On the other hand, the vegetation cover of Plaza de Colón plays a similar role in overall thermal comfort satisfaction. The visitors to the Plaza de Santa Ana and Plaza de Mayor were less satisfied than the Plaza de Colón. At the same time, the Plaza de Mayor and Plaza de Santa Ana have the lowest thermal satisfaction among all urban plazas in Madrid and Seville due to the lack of shade trees and suitable porches or the vast open space. The psychological effects of tourists expecting from the capital’s squares play a considerable role in terms of thermal comfort. In general, shadings of trees and buildings, radiations, and site morphology in these plazas can affect the perceived satisfaction of tourists.

Fig. 7
figure 7

The percentage of satisfaction distribution in each urban plaza in Sevilla and Madrid

Thermal comfort zone of tourists in Sevilla and Madrid

Climatic indicators required to calculate PET were measured when completing the questionnaires and as input parameters in Rayman have been considered. In addition to climatic indicators, the 80-W metabolic rate for a 35-year-old man with a height of 1.75 cm and a weight of 75 was added to the climatic data, and the PET value was calculated for each of the designated points (Höppe 1999). A metabolic rate of 80 watts is a normal rate for someone walking at 1.2 m per second (Mayer and Höppe 1987). According to previous research (Kántor et al. 2012; Mohammad et al. 2021) which stated that metabolism, clothing insulation, height, weight, and age of individuals do not have a significant influence on the calculation of PET. To make an autonomous index of individual conduct, the amount of clothing and activity in PET calculations are kept constant (Höppe 1999). It is fundamental to compute the mean thermal sensation votes (MTSV) for every 1 °C PET interval to comprehend the connection among PET and TSV (De Dear and Brager 2002). In this research, the relevance between PET and MTSV with Eq. 4 (for Seville) and Eq. 5 (for Madrid) is defined in accordance with prior researches (ASHRAE 2001; Gachkar et al. 2021; Nakano and Tanabe 2004; Nasrollahi et al. 2021), with thermal comfort domain between − 0.5 and + 0.5 is acceptable. In this climate, the neutral PET of Sevilla was 28.44 °C when MTSV was equal to zero (Eq. (4)). While in Madrid, substituting 0 in Eq. (5), the neutral PET was 27.18 °C, and finally, the outdoor thermal comfort in Sevilla ranged 28.42 °C to 30.87 °C PET, while the PET range is between 24.5 and 29.82 °C in Madrid (Fig. 8).

Fig. 8
figure 8

Correlation between the mean thermal sensations votes (MTSVs) and PET in Sevilla and Madrid

$$MTSV=0.203 PET-5.794$$
(4)
$$MTSV=0.188 PET-5.113$$
(5)

Simulation results

Figure 9 shows the average Ta, Rh, Ws, and MRT changes at Z = 1.7 m, and the map of Ta, Rh, Ws, and MRT at 12:00 p.m. is shown in Fig. 10 for the different plazas of Seville. The evaluation results of Sevilla plazas showed that the difference of Ta in several sites was less than the various other factors, so that the maximum difference between the lowest temperature in Plaza de España and Alameda de Hercules was 0.37 °C that is due to its subsidence in comparison with other plazas. On the other hand, Rh shows another trend from the Ta. Alameda de Hercules and Plaza Nueva because of high vaporization have the most elevated measure of Rh than Plaza de España. The MRT index depends on the shade of trees (both the tree canopy and height), subsidence of some plazas, and their walls, which in some hours of the day due to lack of suitable vegetation covers and some of them have experienced temperatures up to 75.2 °C. The shadow of the trees played a more important role in reducing the MRT than the shadows of the walls. So that, due to the presence of suitable vegetation in Plaza Nueva (both tree canopy and height) compared to Alameda de Hercules and Plaza de España, MRT was lower than the other two plazas, while the average MRT difference between Plaza de España and Alameda de Hercules has reached 2.32 °C due to lack of proper coverage. In connection with the effect of Ws on each case study, it can be stated that due to the considerable size of Plaza de España and Alameda de Hercules in comparison with Plaza Nueva, the Ws values have recorded higher than Plaza Nueva. However, the relatively taller walls and shallower depths of the Plaza de España have a lower Ws than the Alameda de Hercules. The most suitable OTC was obtained in Plaza Nueva, where relatively good vegetation and small space compared to other fields were the most important factors of superiority over other plazas. Generally none of the plazas were within the PET range and only some of them were in the thermal comfort zone during some hours of the night.

Fig. 9
figure 9

Variation of Ta, Rh, Ws, and MRT in each receptor of the three selected plazas of Seville

Fig. 10
figure 10

Ta, Rh, Ws, and MRT maps at 12:00 for each selected plazas of Seville

The average Ta, Rh, Ws, and MRT changes at Z = 1.7 m and the map of Ta, Rh, Ws, and MRT at 12:00 p.m. is illustrated in Figs. 11 and 12, respectively, for the different plazas of Madrid. Analysis of Madrid squares showed that the difference of Ta in different places due to being at the same level (no difference in height) is even less than the temperature difference in Sevilla plazas, so that the maximum difference reaches up to 0.23 °C. On the other hand, Plaza de colon has the highest Rh content compared to Plaza de Santa Ana and Plaza de Mayor (32.36%), due to its relatively good vegetation compared to other plazas and high evaporation. The least daytime MRT was seen under the plane trees in Plaza de Colón due to more shading coverage (48.46 °C). In addition, the maximum values in MRT (28.17 °C) during the night occurred in this scenario because of forestalling the active radiation to the sky. Conversely, Plaza de Mayor had the biggest degree of radiation during the daytime among all plazas because of the less shadow power, where the MRT was 5.42 °C and 7.22 °C higher than in the Plaza de colon and Plaza de Santa Ana individually. Nonetheless, the greatest decrease in MRT during the night occurred in Plaza de Mayor because of diminished measure of seen long-wave radiation. The most suitable ventilation in Plaza de Colon has been recorded due to its location in an open space compared to other options, while the placement of other plazas in a closed space has led to a decrease in wind speed. The most suitable ventilation among the surveyed plazas is assigned to Plaza de Colón, where due to its location in a more open space than the other two plazas even despite having vegetation cover, the wind speed has increased significantly compared to other plazas. Lack of proper ventilation, small space and limited, and having the most intense thermal stress during sunny hours among other plazas are the most reasons of creating the worst thermal comfort condition in Madrid (Figs. 13 and 14).

Fig. 11
figure 11

Variation of Ta, Rh, Ws, and MRT in each receptor of the three selected plazas of Madrid

Fig. 12
figure 12

Ta, Rh, Ws, and MRT maps at 12:00 for each selected plazas of Madrid

Fig. 13
figure 13

Variation PET between selected plazas of Seville and Madrid in different times

Fig. 14
figure 14

A comparison of thermal comfort (i.e., PET) between selected plazas of Seville and Madrid in different times

Comparison of thermal conditions in different urban sites

To compare the thermal comfort situations in various urban areas, climatic indicators including Ta, Ws, Rh, and MRT were determined in the planned places of each urban plazas, and their information was utilized to estimate PET’s thermal comfort index in Rayman software. The thermal comfort situation in the simulated microclimates of each urban plaza was not quite the same as different destinations during the hours of the day. For example, according to estimates made in the city of Seville, the amount of PET is less than in Madrid, so the difference between the highest value in Sevilla plazas (Plaza de Espana) and the highest amount of pet in Madrid plazas (Plaza de Colón) reaches 4.61 °C.

The thermal comfort schedule has outlined the most reasonable opportunity to visit each popular urban plaza in a day with the mean hottest month of the year (Table 8). The thermal situation of the considered urban plazas was divided into 2-h intervals between 1 a.m. and 12 noon. Each color addresses a 3 °C scope of PET’s thermal index. Albeit, none of the urban plazas is consistent with the PET adequate thermal ranges from 8 a.m. to 9 p.m.; however, these types of decisions are considerably comparative, and destinations with better thermal comfort conditions are positioned appropriately to visit. For example, thermal comfort conditions in the periods 9–11 a.m. show that in the surveyed plazas in Madrid, Plaza de Santa Ana is in the priority of visitors, while Plaza de Mayor and Colon are in the next priorities. On the other hand, in the surveyed plazas in the city of Seville, Plaza Nueva provides the most suitable conditions for thermal comfort for visitors during this period. Alameda de Hercules of Sevilla and Plaza de Santa Ana of Madrid were a better option to visit during the peak hours (11 a.m. to 3 p.m.) of heat than other sites. Semi-open spaces and having suitable trees in Plaza Nueva and Plaza de Mayor created more favorable thermal comfort conditions than other urban plazas and considered of the most important reasons for choosing visitors in the evening.

Table 8  Thermal comfort calendar of the historical sites in average during the hottest month

Comparison of simulation and questionnaire results

In the simulation models, two popular urban plazas with higher thermal comfort conditions (Plaza Nueva and Plaza de Colon) reported higher satisfaction levels in the questionnaires (Table 9). However, three local urban areas with lower thermal comfort conditions (Plaza de Mayor, Plaza de Santa Ana, and Alameda Hercules) got lower fulfillment levels. On the other hand, based on the output of the simulation, Plaza de Espana had the most exceedingly terrible thermal comfort conditions during the day in contrast with other popular destinations, and after Plaza de Mayor and Alameda de Hercules, Plaza de Espana had a higher level of dissatisfaction than different locales. This may be because of individuals’ assumptions from Plaza de Espana Palace. Thus, their satisfaction level with thermal comfort in this urban square diminished. In several investigations, it is obvious that the mental assumptions for individuals influence their understanding (Nasrollahi et al. 2017, 2021; Zeng and Dong 2015). Similarly, Cohen et al. (2013) and Gachkar et al. (2021) researched the expectation of thermal comfort of people in paths, stations, and squares in historical garden, respectively. In both research, psychological effects were observed on their perception of thermal comfort, although with the experience of better thermal comfort in the urban plazas, they anticipated other thermal conditions.

Table 9  Thermal satisfaction of tourist in historical place in hottest month

Discussion

Tourism is so important in the economic and social development of any country, and nowadays economists have called it as invisible exports, because tourism goods and services are offered without the cost of transporting goods for export by attracting foreign tourists inside. Nature tourism is one of the branches of tourism based on the preservation and sustainability of the environment and guarantees the sustainable development of the tourism industry in relation to the needs of tourists and the host community. However, it is highly dependent on climate change and always as one of the most important indices necessary by tourists and desirable thermal comfort conditions and plays an important role in increasing the presence of tourists in urban spaces. Current research analyzed the thermal comfort conditions by using meteorological data irregardless of microclimates created in urban plazas. The current research is conducted to assess thermal comfort conditions in the microclimates of the metropolitan squares of Sevilla and Madrid in a constantly with extreme hotness stress. The thermal comfort sensation of tourists at the urban plazas of Sevilla ranges from 28.42 to 30.87 °C of PET, while the pet range in urban plazas of Madrid is varied between 24.5 and 29.82 °C. Notwithstanding the extreme heat stress, about 38.11% and 28.09% of tourists in Sevilla and Madrid, respectively, were satisfied with the thermal conditions, while 13.14% and 11.65% were dissatisfied. This shows the mental adaptability of travelers with the current situation in spite of the great hotness stress. Albeit, the respondents are from various environments and nations, yet their thermal comfort inclination is near the comfort zones of warm climates (subtropical, hot, humid, and tropical). Based on the simulation results, approximately none of the historic urban areas are in their comfort zone during the day, except from 9 to 11 am in Plaza Colon. These well-known urban plazas don’t have a period impediment to visit and can be visited during the night hours. The thermal comfort situation analogation of each urban plazas is the most important factor regarding choosing the suitable visiting time range.

The overall satisfaction of the tourists with the thermal comfort conditions due to the psychological effects is evaluated in this present study in hottest heat stress month of July. The best time to visit the Plaza de Santa Ana is the early morning hours, while Plaza de Mayor and Plaza de Colon are in the next priorities. On the other hand, in the surveyed plazas in the city of Seville, Plaza Nueva provides the most suitable conditions for thermal comfort for visitors during this period. In the peak hours of heat (11 a.m. to 3 p.m.), the priority goes Alameda de Hercules of Sevilla and Plaza Santa de Ana of Madrid in terms of thermal comfort. Semi-open spaces and having suitable trees in Plaza Nueva and Plaza de mayor provided more favorable thermal comfort conditions than other urban plazas and considered of the most important reasons for choosing visitors in the evening.

It is deduced in this research that urban plazas with higher thermal pressure in the simulation cycle have a higher amount of thermal dissatisfaction in the survey. Albeit both results have a bit different from each other, there are similarities in outcomes. The thermal comfort conditions are predicted through simulation. Indeed, the differences in simulation and field study results, which accounts for the psychological effects, could not be simulated.

Conclusion

This study examined the factors related to the thermal comfort in open space, considering that the tourists can choose the most appropriate time to visit each of the historic sites of Seville and Madrid during the day and night. According to field measurements, questionnaires, and analyses managed in this study, the acceptable range of outdoor thermal comfort for visitors in the well-known urban plazas of Sevilla varies from 28.42 to 30.87 °C in the hot summer, while the PET range in Madrid is between 24.5 and 29.82 °C. Despite the high heat stress condition, about 38.11% and 28.09% of tourists in Sevilla and Madrid, respectively, were satisfied with the thermal conditions, while around 13.14% and 11.65% were dissatisfied. This indicates the psychological compatibility of tourists with the conditions despite the high heat stress. According to the Envi-met different scenario evaluated results, almost none of the surveyed sites were in the range of thermal comfort. However, Plaza de Santa Ana is the best place for visitors in the early morning hours in Madrid and Plaza Nueva in Sevilla at the same time. In addition, during the peak heat hours, Alameda de Hercules of Sevilla and Plaza de Santa Ana of Madrid in terms of thermal comfort are the most suitable historical places for visitors. Concerning semi-open spaces and relatively suitable vegetation, Plaza Nueva and Plaza de Mayor bring more favorable conditions in terms of thermal comfort for visitors in the evening. The result of the present study witnessed that the urban plazas with relatively high thermal stresses in the simulation process have a higher rate of thermal dissatisfaction in the questionnaire. However, there are effective psychological aspects that differentiate the results of simulation and questionnaires because the effects of psychology could not be simulated.