{'fraction_non_alphanumeric': 0.03714956850329367, 'fraction_numerical': 0.02407700556605219, 'mean_word_length': 4.995254859941833, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 5, 'https://': 2, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Weather suitability for outdoor tourism in three European regions in first decades of the twenty-first century

Outdoor tourism and recreational activities strongly depend on actual meteorological conditions. Traditionally, in three studied regions, the peak of tourists’ streams concentrates in summer months. In the present study, we assess suitability of weather conditions for various forms of outdoor tourism in different regions of Serbia, Poland and Ukraine. Additionally, how the location of the station differentiates temporal patterns of weather suitability will be discussed. To analyse the suitability of weather conditions for various forms of outdoor recreation, we have chosen 23 meteorological stations of the national weather networks which represent different tourism areas and destinations. For each weather station, daily data for the period 2000–2017 of air temperature, relative humidity, total cloud cover and wind speed (at 10 m above ground) for 12 UTC as well as the daily maximum and minimum temperature, precipitation totals and snow cover depth were applied. Suitability of climate for outdoor recreation and tourism is assessed by the Weather Suitability Index (WSI) based on Błażejczyk’s bio-thermal weather classification. The results of research show that passive forms of recreation (sun and air bathing) are preferred mostly in months from May till August or September. For the active forms of recreation, weather in summer months is very oppressive, especially in the resorts located in the south (Serbia, southern Ukraine). Active forms of recreation are preferred there in autumn, winter and spring months. Electronic supplementary material The online version of this article (10.1007/s00484-020-01984-z) contains supplementary material, which is available to authorized users.


Introduction
During the last decades, tourism industry has increased all over the world. According to the United Nations World Tourism Organization (UNWTO), international arrivals have raised from 436 million in 1990 to 1400 million in 2018 (UNWTO, 2019). The most frequently visited places are European countries (710 million in 2018).
Tourism is an important branch of national economies. In Poland, in the years 2006-2016 a number of hotels and other tourism objects has increased to 72%. It is about 2% of the European market. The number of domestic tourists has increased to 57% and international tourists to about 48% (Łysoń 2019). The statistical data for the whole Ukraine cover the period 2000-2013 (2000-2010 based on data from the Ministry of Infrastructure of Ukraine, http://mtu.gov.ua, since 2011-from to the State Statistics Service of Ukraine (2019). From 2014, statistical records do not include the temporarily occupied territories of Crimea, Donetsk, Sevastopol and Luhansk. Thus, in the years 2000-2013 about 55% increase in tourists (61% for foreign and 52% for domestic tourists) was noted, and the last years brought a decrease in tourism activity. In Serbia, the number of tourists has increased in the last 17 years (with decreasing and stagnation phases in the years [2009][2010][2011][2012][2013] and in 2018, 9.3 million overnight stays were registered which is 12.1% more than in 2017. The share of foreign tourists in the total number of overnight stays in 2008 was 19.1%, while in 2018 increased to 39.2% (Statistical Yearbook of the Republic of Serbia 2017).
The time series of tourists' nights contains very pronounced seasonal fluctuations in three studied regions, with the highest values in the summer season and in January. There are several factors that influence tourism activity. One of the main reasons of increasing tourism is the growing richness of population and transportation availability.
There is seen a great seasonality in tourism over the year. In Poland and in Serbia, the greatest number of tourists is registered in the summer months. In 2018, about 34% hotel services of domestic and foreign tourists concentrate in the months June-August and 53% in the period of May-September (Łysoń 2019). The similar seasonal pattern of tourists' traffic is also observed in Serbia (Statistical Yearbook of the Republic of Serbia 2017). Such concentration of tourism activity has different reasons (e.g. organisation of school year, holiday habits), and one of them is climate conditions assessed as appropriated for recreation (Andrade et al. 2007, Hamilton 2007, Mansfeld et al. 2007, Scott et al. 2007, Vrtačnik Garbas 2007. There are hundreds of research and papers considering the influence of various climate elements on the different kinds of tourism. In climate-tourism research, the most frequently used are physiological equivalent temperature (PET, Höppe 1999), perceived temperature (PT, Jendritzky 1990, VDI 2008, apparent temperature (AT, Steadman 1984), heat index (Rothfusz 1990), effective temperature (ET, Baranowska and Gabryl 1981), physiological subjective temperature (PST, Błażejczyk 2007), predicted mean vote (PMV, Fanger 1970), and Universal Thermal Climate Index (UTCI, . The most frequently used indices were listed and discussed by Epstein and Moran (2006),  and de Freitas and Grigorieva (2017). While tourism integrally corresponds to the change of the place of stay, special adaptation indices are also used.
Climate information provided to tourists and tourism industry should include not only its general features (mean monthly and annual values of particular climate elements) but also a detailed information regarded to thermal comfort as well as aesthetic and physical weather factors (de Freitas et al. 2008). All of these factors have to be included in an integrated assessment of tourism climate (Matzarakis 2007). These requirements of climate-tourism research represent weather approach which considers the complex of meteorological elements (i.e. the actual weather) influencing the human organism. Several proposals of weather classification have been developed so far (Błażejczyk and Kunert 2011). More recently, there are proposals from Błażejczyk (2007) and Andrade et al. (2007) that are based on human heat balance considerations. There are also the in use DWD (German Weather Service) synoptic classification of weather (Bissoli and Dittmann 2001;Bucher 1991) and Spatial Synoptic Classification (SSC, Kalkstein et al. 1996, Sheridan 2002. However, DWD and SSC classifications are dedicated for general assessment of bioclimate and not for climatetourism research. An integral assessment of weather and climate information for tourism is possible by analysis of the Climate-Tourism-Information-Schemes (CTIS) (Matzarakis 2007), Climate Index for Tourism (CIT, de Freitas et al. 2008) and Weather Suitability Index (WSI, Błażejczyk 2007). All base their data on daily meteorological information and human heat balance measures. However, Mieczkowski's (1985) Tourism Climate Index (TCI) bases on the mean monthly climate data. It is used (after some modifications) for the general assessment of climate usefulness for tourism in continental and global scales (e.g. Amelung and Nicholls 2014;Kovács and Unger 2014).
Temporal and spatial variability is one of the essential features of climate experienced during tourism activity. Tourists' organisms must adapt to the changing climate stimuli. In general, an increase in differences of climate stimuli intensifies the magnitude of adaptation processes in an organism (Jendritzky and de Dear 2008;Koppe and Jendritzky 2005). There are few attempts of quantitative assessments of bioclimatic contrasts between different sites. The concept of Bioclimatic Distance (BD, Mateeva and Filipov 2003) bases on the comparison of clothing necessary for keeping the heat balance of the human body. The Acclimatization Thermal Strain Index (ATSI) of de Freitas and Grigorieva (2009) is based on heat loss by respiration. Błażejczyk and Vinogradova (2014) proposed the Adaptation Strain Index (ASI) which includes three measures of adaptation processes of an organism to different ambient conditions: water loss, insulation predicted and Universal Thermal Climate Index.
There are many studies reporting bioclimatic potential for recreation and tourism in Poland and in Serbia. General characteristics of whole Poland were done in several publications of Błażejczyk and co-authors (Błażejczyk et al. 2015. There are also researches assessing bioclimatic conditions for tourism in different regions of Poland. As an example of research using WSI approach, we can list papers of Radzka and Dragańska (2015), Miszuk (2008), and Mąkosza (2013). However, Mąkosza andMichalska (2010, 2011) have applied human heat balancebased indices (Subjective Temperature Index, STI and Physiological Strain, PhS) to assess bioclimatic potential for tourism and recreation in western Poland. Recently, several papers have been published in Serbia considering heat budget indices in the context of tourism and recreation. Pecelj et al. (2013) analysed the Weather Suitability Index (WSI), physiological subjective temperature (PST) and physiological strain index (PhS) for different physical activities in tourism context, and Basarin et al. (2018) used physiological equivalent temperature (PET) to study bioclimatic changes in Novi Sad. Bioclimatic conditions of Zlatibor, Belgrade and Loznica were analyzed referring to the Universal Thermal Climate Index (UTCI) by Pecelj et al. (2017) and Pecelj et al. (2018). Stojićević et al. (2016) studied the bioclimate condition of Banja Koviljaca Spa using PET.
To the best of our knowledge, for Ukraine, there are several researches addressed to the influence of climate on tourism activity. For instance, a current state of bioclimatic conditions and their usefulness for tourism in Odesa (the coast of the Black sea) was considered by Katerusha and Matzarakis (2015) by applying the CTIS approach. Seasonal and multiannual changes of PET in Kyiv in the period 1961-2015 was assessed by Shevchenko et al. (2020). General information regarding climate resources can be found in tourism geography and tourism economy publications where authors pay attention to the importance of climate and weather in recreational balneology and spa therapy (Beydyk 2001;Hladkyi and Mirzodaieva 2018).
There are also some attempts to compare bioclimatic conditions in various parts of Europe,-e.g. Nicholls and Amelung (2008) assessed tourism potential using the TCI for Nordic countries and Kovács and Unger (2014) for the central European region. Błażejczyk and Matzarakis (2008) used PST, WSI and CTIS to evaluate climate suitability of Helsinki, Athens, Paris and Cracow. Błażejczyk and Vinogradova (2014) studied the adaptation strain of tourists travelling from central and northern Europe to the Mediteranean. Several papers comparing bioclimatic conditions in various sites in Europe were also published by Błażejczyk and collaborators Błażejczyk et al. 2015).
The aim of the present study is to assess suitability of weather conditions for various forms of outdoor tourism and recreation in different regions of Serbia, Poland and Ukraine and to compare how location of the station differentiates temporal patterns of weather suitability.  Fig. 1). When selecting the weather stations, we have taken into consideration their location due to local landscape which significantly modify regional features of climate . All the stations are located outside the urban areas and represent regional features of climate. The influence of local factors (as slope exposure, land use, urbanization level, etc.) can be neglected (Table 1).

Materials and methods
In general, the stations represent two classes of climate based on the Köppen-Geiger classification (Peel et al. 2007): -Temperate climates (C-temperature of the hottest month >10°C and the coldest one from 0 to 18°C), without dry season (f) and with hot (a), warm (b) or cold (c) summer, -Cold climates (D-temperature of the hottest month >10°C and the coldest one ≤0°C) without dry season (f) and with hot (a), warm (b) or cold (c) summer For each weather station, daily data of air temperature, relative humidity, total cloud cover and wind speed at 10 m above ground for 12 UTC was applied. We have also used daily data of maximum and minimum air temperature, precipitation totals and snow cover depth. The observational and homogenised data provided by the national weather services of Poland, Serbia and Ukraine cover the period 2000-2017 (the exception is Yalta, where due to temporal occupation of Crimea, data are available only till 2014). For the entire period, data were 100% complete. Suitability of climate for outdoor tourism and recreation is assessed by the weather suitability index (WSI). It bases on daily meteorological information. WSI is assessed based on bio-thermal classification of weather proposed by Błażejczyk (Błażejczyk 2007, Błażejczyk and Matzarakis 2007, Błażejczyk and Kunert 2011. Weather information is coded by seven digits which represent the following: 1st-thermal sensations, 2ndradiation stimuli, 3rd-physiological strain, 4th-sultriness, 5th-daily thermal contrast, 6th-rain (snow) fall (>1 mm), 7th-snow cover (>10 cm) ( Table 2). Weather characteristics from 1 to 4 refer to bio-thermal conditions during midday hours. Daily thermal contrasts are assessed by the difference between daily maximum and minimum temperatures. In case of precipitation, the best information would be the time of its occurrence during the day. Unfortunately, numerical databases provide only their daily totals and not the time of its occurrence. Thus, rainy/snowy days (i.e. days with precipitation >1 mm) were considered as suitable with limitations for outdoor tourism. The snow cover classification uses criterion of 10-cm depth which is considered as suitable for ski tourism, both in mountain and lowland resorts. For example, the code -2_2C0_011 indicates cold weather (-2) with moderate radiation stimuli (2), cold physiological strain (C), non-sultry (0) with weak daily thermal contrast (0), rain(snow) fall (1) and snow cover (1).
WSI provides evaluation of each individual weather conditions (with daily resolution) from the point of view of particular forms of recreation: sun bathing (WSI_SB), air bathing (WSI_AB), mild physical activity (e.g. walking, light plays, shopping-WSI_MR), intensive physical activity (e.g. football, biking, climbing, jogging, etc.-WSI_AR), Because of its nature, WSI is not calculated with a mathematical formula but by using a lookup table (see WSI configuration file in BioKlima©2.6 software package, https://www. igipz.pan.pl/bioklima.html). Such a table contains the WSI_ xx (where xx = SB, AB, MR, AR, ST) values corresponding to every weather class. Every weather situation is assessed using WSI as follows: 0-unfavourable, 1-favourable (with limitations), 3-very favourable. Table S1 in supplementary  materials gives an example of lookup table. For the particular periods (months, year) individual WSI values for different tourism activities are averaged and categorised as follows (Table 3).
In the present study, we have considered monthly and yearly number of days grouped from categories: moderately suitable, suitable and very suitable. We also verified if mean annual values of this group of categories are statistically  Table 1). Source: own elaboration at the map of San Jose-own map, based on the Generic Mapping Tools and ETOPO2, CC BY-SA 3.0, https://commons.wikimedia.org/w/index. php?curid=676986, maps of particular countries at the background of Google maps significant within the entire country and within different types of landscapes, namely coastal, lowland, lowland river valleys, mountain and mountain food. To discriminate among the means of WSI_SB, WSI_AB, WSI_MR, WSI_AR and WSI_ST, the Fisher's least significant difference (LSD) procedure at the 95.0% confidence level was applied. The statistical analysis were made with the use of STATGRAPHICS Centurion XVI software package.
The mean monthly values of WSI_SB, WSI_AB, WSI_MR, WSI_AR and WSI_ST were applied to distinguish periods of different suitability (unsuitable, moderately suitable, suitable and very suitable) for various forms of tourism. Such periods are compared between regions and stations.
For the general assessment of weather suitability for tourism, the average value of WSI is applied (Table 4). WSIavg is calculated with the use of following equation:

General features of weather suitability
In Poland, sun bathing can be practised from about 51 days a year in mountain resorts (HG) up to ca. 139 days in coastal locations (HEL) and 140 days in central Poland (WAR). In Ukraine, sun bath weather varies from about 112 days a year at Chornohora mountain ridge (POZ) to 149 days at coastal Weather type −3 (Very cold) Weather subtype Weather indicators set in italics Source: Błażejczyk (2007) area of Crimea peninsula (YAL). In Serbia sun baths can be practised from about 77 days a year in Pannonian lowland (NS) up to 111 days in Mt. Zlatibor. Number of days with weather suitable for air baths varies in Poland from 93 in mountains to about 175-180 in the central lowland (TOR, WAR). In Serbia, such weather is more frequent than in Poland and varies from 146 in Belgrade up to 181 days in Loznica. In Ukraine, occurrence of days of air bath weather is the greatest and change from 164 in the northern part of the country (SVI) to more than 190 days at southeastern locations (AN, YAL).
Mild recreation activity in Poland can be practised during 274 days in the mountains (HG), 308-312 days in lowland locations (WAR, WRO, TOR) and up to 320 days in coastal resorts (HEL). In Serbia, mild activity days vary from 138 in northern lowland (NS) and up to 188 days a year in mountain resort (ZLA). The most frequent are MR days in Ukrainian tourism areas, and they vary from 270 yearly at the mountain foot (KHU) to 299 at the Black Sea coast (ODE).
Regarding days with weather suitable for active recreation (AR), their greatest number occurs in Poland in the mountains (HG-almost 361 and JG-342 per year). In other resorts, it fluctuates from 324 (KRA) to 337 (HEL). Less frequent AR days are noted in Ukraine and vary from 260 at the Crimea coast (YAL) to 337 at the mountain station (POZ). The smallest number of AR days occur in Serbian resorts, from 112 yearly in Niš to 159 in Mt. Zlatibor.
Ski tourism can be practised mostly in elevated mountain locations. For example, ST suitable weather occurs at the HG station 154 days, at POZ, 115 days and at ZLA, 33 days a year. In mountain valleys (JG), such weather occurs about 31 days a year. Similar number of ST days is observed at eastern locations in Ukraine (KYI, KHU, UMA, LVI). However, at southern locations, both in Serbia and in Ukraine, there are less than 8 days a year with ST weather (Table 5).
There were analyses on how weather conditions for tourism differ inside studied areas. In this purpose, differences between annual numbers of days suitable for particular forms of outdoor tourism were compared. The most unified weather conditions occur in Serbia which is the smallest, compact country. Frequency of days suitable for various forms of tourism differ significantly only between 20 and 60% of stations. Only the weather in Zlatibor Mt. is significantly different than in other stations. The greatest spatial variation (for 60-100% of stations) is observed in the case of MR and AR weather frequency. The most differentiated weather occurs in Ukraine, the biggest country with great longitudinal and latitudinal extent. This is the cause that in every station, weather differs significantly in comparison to 60-100% of the remaining locations. The greatest weather individuality, in comparison to 82% of other stations, is observed in the mountain area of Eastern Carpathians. The biggest spatial conformity (only 53% differences with remaining stations) occurs in the case of MR weather. In Poland, the most unique weather is observed in mountain areas (HG, significantly different than in all other stations). Comparing types of tourism, the greatest individuality was found for sun bathing weather (Table 6). For a detailed information of inter-station differences in the number of days suitable for SB, AB, MR, AR and ST, see Tables S2, S3 and S4 in supplementary materials.

Seasonal cycle of weather suitability
Seasonal patterns of weather differ between regions and between locations. In general, the greatest seasonality is well seen for sun (SB) and air (AB) bathing. In winter months (DJF) in majority of the locations, there are not days suitable for sun and air baths. They sporadically appear only at stations situated in Cfa and Cfb climate zones (YAL, NS, BEL, LOZ, NIS). However, in winter months, very frequent are days with weather suitable for mild (MR) and active (AR) recreation (20-30 days per month, depending on station). In the summer season (JJA), frequency of sun and air bathing days increases. In majority of the stations in the summer months, one can observe a great decrease in days suitable for MR and AR activities. It is caused by increased air temperatures during particular summer days. Ski tourism (ST) can be practiced mostly in mountain tourism resorts from December till March or April.
The above remarks indicate only main seasonal features of weather in temperate and cold climate zones. Most of them is typical for stations located at lowland landscape (both, at river valleys and flat areas). The seasonal changes of SB and AB   (2007) weather is the best seen in most northern locations (e.g. WAR, Poland). In most southern locations (e.g. BEL, Serbia), a very strong reduction of days with weather suitable for all forms of recreation is pronounced in summer. At lowland Ukrainian stations, seasonal patterns of weather differ depending on longitudinal location. While in Kyiv (KYI, central part of the country), seasonal weather patterns are similar to those observed in Poland then in the eastern location (MAR), the winter to summer differences in SB, AB, MR and AR are the greatest. It is probably caused by exposition of that region to severe cold arctic air mass in winter and to severe hot tropical air mass in summer months (Fig. 2). Increased reduction of days with MR and AR weather is also observed at resorts located in mountain foot areas in Ukraine and Serbia (KHU, LOZ). Moreover, in most southern resorts (LOZ), a great decrease of SB and AB days is pronounced as well (Fig. 3).
Regarding mountain stations, there are similarities between Polish Tatras (HG) and Ukrainian Chornohora (POZ). In the summer season, there is an observed culmination of days with SB (12-20 days monthly), AB (17-20 days) and MR (28-31 days monthly) weather. A slight reduction of AR days (to about 22 in August) is observed only in Ukrainian Carpathians (POZ). A contrasting seasonal pattern of weather suitability occurs in Serbian mountains (ZLA). SB and AB weather are most frequent in spring (March-April) and in autumn (September-October), up to 30 days monthly. In summer, such weather is considerably reduced to 15-20 days monthly. At the same time, active recreation can be practised only   For abbreviations of stations see Table 1 Source: own elaboration For abbreviations of stations, see Table 1 Source: own elaboration  Table 1) Source: own elaboration   Table 1) Source: own elaboration  Table 1) Source: own elaboration ST conditions occur also in December. In Serbia (ZLA), one can plan an ST activity only during 20 days in January, 10 days in February, 8 days in December and 5 days in March (Fig. 4). Seasonal patterns of weather at coastal resorts can be analysed only for Poland and Ukraine, while Serbia has no access to sea coast. When comparing Polish (Hel) and Ukrainian (Odesa, Yalta) stations, one can see great differences in the annual run of weather suitable for various forms of recreation. At the Baltic coast (HEL), summer decline is not observed in mild recreation weather types (MR), and reduction of the AR weather is relatively small. There is also experienced summer peak of the frequency of SB and AB weather. At the Black Sea coast, summer decline of AR and MR weather is great (due to high temperatures), and in Yalta, active forms of recreation can be practised only 5-10 days a month. Slightly less severe summer weather is observed in Odesa. Typical for this region is also summer reduction of weather suitable for sun and air bathing. In ODE, SB and AB can be applied for about 20 days a month and in YAL, only 12-15 days monthly (Fig. 5).

Types of landscape and weather suitability
The analysis of annual cycles of the frequency of days suitable for various forms of tourism suggest significant differences between stations located in similar landscape types but in different areas. Thus, statistical significance between annual frequency of SB, AB, MR, AR and ST weathers was verified. At coastal stations, weather conditions differ significantly between Polish (HEL) and Ukrainian (ODE, YAL) resorts for almost all forms of tourism. The same is the concern in lowland, mountain foot and mountain summit stations in studied countries ( Table 7). The exception is weather suitable for ski tourism which often does not differ between mountain foot Differences statistically significant are marked in yellow (SB sun bathing, ABair bathing, MRmild physical activity, ARintensive physical activity, STski tourism) between coastal, lowland, mountain foot and mountain summit stations, 2000-2017 (for abbreviations of stations see Table 1) Source: own elaboration stations in compared areas. An opposite situation was found in the case of stations located in lowland river valleys. Very often, the annual number of SB, AB, MR, AR and ST days do not differ significantly. There are surprisingly similar weathers at stations located far away from each other (e.g. TOR and MAR or TOR and BEL). This needs more detailed research in the future. Mostly differentiated are weather conditions suitable for active forms of tourism (AR) which do not differ only between TOR and WAR but also between WAR and WRO (Table 8).

Sun bathing
The longest period of sun bathing was found for Serbian stations. In majority of tourism resorts, it lasts 12 months a year with different suitability. In all stations, very suitable weather occurs in spring (March-April). In autumn, sun bathing period covers September and October (in Niš, October and November). Sun bathing is unavailable only in three resorts: NS and ZLA in December-January (because of the increased number of cold days in Mt Zlatibor and frequent winds in the Pannonian region during the winter) and in NIS in July (due to increased number of very hot days). In Poland and in Ukraine, sun bathing period is considerably shorter than in Serbia and lasts 7-8 months (mostly from March till October), except elevated parts of mountains (HG) where moderately suitable period covers only 4 months (May-August). It is necessary to add that in the southern part of Poland dominates the moderately suitable sun bathing weather. In central and northern Poland, weather suitable for sun baths covers the period from April till September. Such suitable weather in Ukrainian tourism resorts is mostly observed in spring (April-May) and in autumn (September-October) (Fig. 6).

Air bathing
The longest AB period occurs in Serbia, and it lasts the whole year round. However, the best months for air bath are April and May as well as September and October. In winter Differences statistically significant are marked in yellow (SBsun bathing, ABair bathing, MRmild physical activity, ARintensive physical activity, STski tourism) between stations located in lowland river valleys, 2000-2017 (for abbreviations of stations see Table 1) Source: own elaboration (December-January), the conditions are only moderately suitable due to relatively great frequency of rainy and cool weather. Again, a very similar seasonal pattern of SB weather was found for Poland and Ukraine. In general, air bathing can be practised there from March till October. The shortest (from May till September) is AB period at elevated mountain resorts  Table 1) Source: own elaboration (HG) and the longest one at the Crimea coast (YAL), from February till November (Fig. 6).

Mild recreation
In Poland, mild recreation can be practised the whole year round, and only in winter (December-January), WSI_MR values indicate moderately suitable conditions (in high mountains, such period lasts from October till March). In the northern part of the region, the period of very suitable conditions occurs in April and September. Similar to Poland, seasonal pattern of mild recreation period is observed in Ukraine. Unsuitable conditions occur only in very southern-located Yalta region in July and August due to very high air temperature.
In summer (mostly in July and in some locations in JJA), weather unsuitable for mild recreation is observed in Serbian tourism resorts, except mountain station at Mt. Zlatibor. However, the period of very suitable MR weather in Serbia lasts generally from October till April and is considerably longer then in Poland and Ukraine (Fig. 7).

Active recreation
The best conditions to practice active forms of recreation occur in the Polish tourism resorts. In the period from May till August, weather is suitable for active recreation. However, very suitable weather lasts generally from September till April.
Slightly worse AR conditions occur in the Ukrainian resorts. Very suitable weather lasts from October till April (1 month shorter than in Poland). In the summer months (JJA), the conditions are only moderately suitable for active recreation (the exception is Chornohora mountains, POZ). In two stations (KHU and YAL), July and August are unsuitable for AR.
The shortest AR period was found for the Serbian tourism resorts. Very suitable conditions occur only from November till February and unsuitable weather predominates from May till August or even September (Fig. 7).

Ski tourism period
Ski tourism period can be considered only in the case of mountain resorts. In Poland (HG) and in Ukraine (POZ), ST period lasts from December till April. In Serbia (ZLA), weather conditions suitable or moderately suitable for ski tourism are 1 month shorter (December-March, Fig. 8).

General assessment of weather suitability
The results presented in the previous section show great seasonal differentiation of weather in 3 studied areas: Poland, Serbia and Ukraine. The differences are seen in the case of every form of tourism. To compare how in general weather differs between 3 studied areas, the spatially averaged WSIavg values were calculated for particular countries. In every country, there are observed two maxima (spring and autumn) and two minima (winter and summer) of weather suitability (Fig. 9). The greatest annual variability occurs in Serbian tourism resorts. From October till April, with peaks of excellent weather in March and November, conditions for tourism are much more suitable than in Poland and Ukraine. However, in summer months (June, July, August), general conditions are unsuitable for tourism, especially its MR and AR forms. In Polish and Ukrainian resorts, annual variations of WSIavg are less manifested than in Serbia. The best, very suitable conditions occur in April and September/October. However, in winter months, weather is at least moderately suitable or suitable for tourism and in summer, suitable (in Poland) and moderately suitable (in Ukraine). These seasonal variabilities of weather allow for reasonable planning of touristic visits depending of individual preferences of tourists.

Discussion
The results of present research show great seasonal and spatial differences in weather suitability for tourism in compared  regions. They are in concordance with the results of previous climate-tourism research carried out in Mediterranean and Central Europe countries. As an illustration, great summer heat stress in Greece and other Balkan countries that cause dramatically unsuitable conditions for tourism was reported by Matzarakis and Mayer (1997), Matzarakis et al. (2014), Zaninović and Matzarakis (2007), and Zaninović et al. (2006. In the cited research, thermal conditions were assessed using PET and CTIS indicators. Similar finding was also reported by Basarin et al. (2018) for Novi Sad where the number of the days above certain thresholds is increased together with the number of heat waves mostly during the last 15 years. Also, Stojićević et al. (2016) observed a growing trend of annual and seasonal values of PET for 1961-2014 in Banja Koviljaca Spa. Summer decrease in weather suitability for active tourism, both using PET and CTIS indicators, was reported for Odessa (Black Sea coast) by Katerusha and Matzarakis (2015) where period suitable for tourism lasts only from the end of April to mid-May and from the end of September to mid-October (five 10-day periods). Shevchenko et al. (2020) reported that in Kyiv, the period with comfortable bioclimatic conditions is much longer than in Odessa and continues about nine 10-day periods.
The stressful thermal conditions in the summer months in southern Europe were also confirmed when using other biometeorological indices. For instance, analysis of heat load index in Serbia shows that heat stress condition for July 2011 was reported for six different weather stations (Pecelj et al. 2013). Increased heat load values were also observed in 2000 and 2007 in Loznica, Belgrade, Niš, Vranje, Zlatibor and Novi Sad (Мilovаnоvić et al. 2017). Lastly, in the framework of the geo-ecological evaluation of the city of Loznica for the purpose of recreational tourism, the UTCI (2000-2016) shows a marked stress in 2007 in July (Pecelj et al. 2018).
Summer decline in weather suitability for mild and active recreation founded in the present research for Poland and northern Ukraine was previously reported in many papers for Polish tourism resorts (Błażejczyk and Kunert 2011;Mąkosza 2013;Mąkosza andMichalska 2010, 2011;Miszuk 2008;Radzka and Dragańska 2015). Similar summer worsening of MR and AR weather were also detected for other European countries (Błażejczyk and Matzarakis 2008) located mainly in Cfb and Dfb climate zones, e.g. for Austria (Koch et al. 1992;Matzarakis et al. 2007Matzarakis et al. , 2010, Slovenia and Croatia (Zaninović et al. 2006) and for Netherlands (Moreno et al. 2009). All authors link such summer decrease in weather suitability in southern and central Europe with advections of hot tropical and continental air (Gajić-Čapka and Zannović 1997; Koch et al. 1992;Kolendowicz et al. 2018;Nowosad et al. 2013) and with geopotential height (Tomczyk and Bednorz 2019). It makes disturbances in the regulation of heat budget, especially during physical activity (Błażejczyk and Kunert 2011;Epstein and Moran 2006;Fanger 1970;Höppe 1999;Jendritzky 1990).
In general, sun bathing is one of the most expected forms of summer recreation and tourism. However, everyone must remember that in European region, there exists the problem of high values of UVI in summer which can cause health problems, including skin cancer Global Solar UV Index 2002). In tourism, resorts should have available information about actual UVI intensity and recommendations on what to do in case of high UVI values.
Our results confirm that tourism is a climate-dependent sector of economy (Andrade et al. 2007, Hamilton 2007, Mansfeld et al. 2007, Scott et al. 2007, Vrtačnik Garbas 2007. Thus, in its future development, the consequences of recently observed climate changes must be taken into consideration (Amelung and Nicholls 2014;Nicholls and Amelung 2008;Endler et al. 2010). The period of climate records used in present research is too short for analysis of changes in trends and to formulate possible projections. However, comparison of data from regions of different climate, especially thermal, regimes make it possible to speculate that if global warming will continue current trends, it can significantly impact weather conditions which influence thermal comfort in tourism resorts. In the studied region, we should expect increase in heat stress days, as is seen in predictions made for Poland (Błażejczyk et al. 2018). Increase in heat stress days is also predicted globally by Pappenberger et al. (2015). This can lead to reduction of weather suitability for active forms of outdoor tourism and this changes must be taken into consideration by tourists and tour operators.

Conclusions
As outdoor tourism and recreational activities strongly depend on actual meteorological conditions, it is very important to evaluate how weather can influence seasonal and regional streams of tourists and recreationists. According to the results of this study, outdoor recreation should be presented through the active (MR, AR and ST) or passive (SB and AB) forms of recreation.
The present research shows great seasonal variability of weather assessed from the point of view of outdoor recreation. The passive forms of recreation (sun and air bathing) are preferred mostly in months from May till August or September. This is the season usually used for any holiday travels and stays. However, for the active forms of recreation during the summer months, weather is very oppressive, especially in the resorts located in the south of the studied region (Serbia, southern Ukraine). Thus, active forms of recreation are preferred there in autumn, winter and spring months.
Taking into account the weather and climate information that may be relevant for tourists, the results of the research suggest that the southern part of the region is a very good place for passive summer holiday stays. Otherwise, tourists preferring active forms of recreation should consider summer holiday visits in resorts located in northern Ukraine and Poland.
The present research dealt with regional differentiation of climate features for climate. However, in planning of future development of tourism resorts, local authorities should consider also local factors (e.g. green, forests and water areas, shadow, etc.) which can modify the microclimate and help in better adaptation to climate change. This should be the topic of future research.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.