Temporal variation of allergenic potential in urban parks during the vegetation period: a case study from Bratislava, Slovakia

Park greenery represents an oasis for urban residents; however, during the flowering period of trees that produce allergenic pollen grains, these areas threaten individuals suffering from seasonal allergic respiratory diseases. In this study, we evaluated the temporal distribution of the allergenic potential of three most important urban parks in Bratislava over the vegetation period, using a modification of the Urban Green Zone Allergenicity Index (IUGZA) and Individual-Specific Allergenic Potential Index (IISA) designed as a running index — rIUGZA and rIISA. We found that rIUGZA gives better information for park management and revitalization, since it considers the potential size of woody plants, while rIISA, considering the actual size of the vegetation, provides more relevant information for pollen-allergy sufferers. Based on rIISA, the allergenic potential was highest in May for the Grassalkovich Garden (formal baroque garden) and Janko Kráľ Park (English landscape park) and in April for the Medic Garden (repurposed baroque garden). We also found differences in the duration of the period of increased allergenic potential in these parks, ranging from 1 to 3 months. Based on the total annual sums of rIISA, we found the highest allergenic potential in the Medic Garden and lowest in the Janko Kráľ Park. This variance is caused mainly by the different density of trees and percentage of allergenic species. The biggest contributors to the allergenic potential were Platanus, Acer and Tilia. Based on the information on temporal variation of the allergenic potential during the vegetation period provided by the running indices, it is possible to improve the planning of park revitalization based on the flowering period of allergenic species and provide better information to the pollen-allergy sufferers for minimizing the allergenic effect of urban green areas on their health during a particular month.


MZ
Medická záhrada, Medic Garden SJK Sad Janka Kráľa, Janko Kráľ Park GZ Grasalkovičova záhrada, Grassalkovich Garden I UGZA Urban Green Zone Allergenicity Index I ISA Individual-Specific Allergenic Potential Index PAV Potential allergenic value of a given area ap allergenic potential of a given plant taxon tp type of pollination ppp length of pollen production period cmf calendar months of flowering I UGZA ′ Urban Green Zone Allergenicity Index calculated with cmf instead of ppp I ISA ′ Individual-Specific Allergenic Potential Index calculated with cmf instead of ppp PAV′ potential allergenic value of a given area calculated with cmf instead of ppp rI UGZA Running Urban Green Zone Allergenicity Index rI ISA Running Individual-Specific Allergenic Potential Index Responsible Editor: Lotfi Aleya

Highlights
• Running allergenicity indices monitor change of allergenic potential during vegetation season.
• rI UGZA considers potential size of vegetation, which is better for park management.
• rI ISA considers actual size of vegetation, which is better for pollen-allergy sufferers.

Introduction
Green spaces are a desirable component of the urban environment (Kothencz et al. 2017), but can also harm human health as green spaces are sources of allergenic pollen grains (Cariñanos et al. 2019).Allergic diseases are rising globally, with a vast clinical and economic impact.The World Health Organization (WHO) estimates the number of cases of allergic asthma to be 300 million and allergic rhinitis to be 400 million (Bousquet and Kaltaev 2007).The total cost of treatment of allergic diseases in the European Union stands between 55 and 151 billion EUR annually (Zuberbier et al. 2014).Allergenic airborne bioparticles, often embodied in longrange atmospheric transmission, are in urban environments subject to long-term monitoring by a network of aerobiological stations throughout Europe (Buters et al. 2018).However, a specific feature of larger green areas with a significant tree cover, such as urban parks, is the trapping of different air pollutants including pollen, so that the concentration of these particles in such an area may differ profoundly from the outdoor situation (Escobedo et al. 2011).Most of the allergenic pollen grains in this area originate from local vegetation, primarily trees and grasses.Therefore, pollen concentrations provided by a network of aerobiological stations cannot be considered an appropriate proxy for the allergenic potential of urban parks.Depending on the phenophase, the allergenic potential of urban parks may be lower or higher than their surroundings.
The allergenic potential of different green areas in the urban environment can be quantified using the Urban Green Zone Allergenicity Index (I UGZA ) (Cariñanos et al. 2014), firstly used in Spain to assess allergenic potential of different urban parks.This method has so far been used to calculate the allergenic potential of urban greenery principally in areas with Mediterranean climate in Spain (Cariñanos et al. 2016, 2017, Cariñanos and Marinangeli 2021;Velasco-Jiménez et al. 2020;Sabariego et al. 2021), Italy and Portugal (Cariñanos et al. 2019;Suanno et al. 2021) and Turkey (Kara and Aşık 2022).In temperate regions, it has only been applied less often, e.g. in Poland (Kasprzyk et al. 2019) and Germany (Jochner-Oette et al. 2018).Jochner-Oette et al. (2018) also used a method taking into account the actual measurements of the trees instead of their potential size, creating the Individual-Specific Allergenic Potential Index (I ISA ).
The problem with insufficient green areas in the city centres is common for many cities, including Bratislava (Belčáková et al. 2022).The oldest parts of these cities follow the original building plan from centuries ago, leaving no space for new green areas, while the plots in the newer parts of the city centres are highly sought after by investors and building parks here is economically less profitable (Winkler et al. 2023).Often the only type of greenery in these areas are historical parks and gardens, from which some could be kept in their original state while others were modernised.With the increase of sedentary lifestyle, the function of urban parks as areas of active recreation is getting more important (Zhang and Zhou 2018).Paths are used by runners and bigger grassy areas can serve for sports like yoga, martial arts or badminton.Even children need more physical activity in their schedule, and playgrounds in parks can fulfil this need (Floyd et al. 2011).Connected to this is the need of expressing the changes in the allergenic potential of these spaces continuously throughout the vegetation season so that the health benefits of outdoor activity are not negated by the allergenicity of the area.Individuals sensitive to respiratory allergic diseases need this information for making a decision about the most suitable place for their recreational activity, and the park management can use it to issue warnings and recommendations as well as planning to avoid peak in allergenicity that may be undetectable by I UGZA or I ISA .For example, a park with a low number of allergenic species can have a low value of I UGZA /I ISA and still be highly allergenic at a particular time of the year due to the simultaneous flowering of some of these species.For this reason, we adjusted both of these indices to create a running index offering the information about the changes of the allergenic potential of these parks through the whole vegetation season.
The aim of this study was to determine the temporal variation of allergenic potential in the three most important urban parks in the city centre of Bratislava during the vegetation season, using running indices based on I UGZA and I ISA .Additionally, we identified the woody plant species that contribute the most to the allergenic potential of these areas and proposed measures to mitigate the negative impact of these spaces on individuals with pollen allergies.

Study area
The three following public parks in the city centre of Bratislava (Fig. 1) were picked as the study locations.
Medická záhrada (MZ, Medic Garden) was founded in 1770 as a baroque garden, but nowadays almost nothing is left from the original historical disposition after a last major reconstruction in 1985 (Steinhübel 1990).It is surrounded by the highly urbanised area of the city centre, and as Fig. 1 The position of the studied public parks in the city centre of Bratislava, Slovakia (source: https:// mu-basm.gispl an.sk/ mapa/ histo ria) 1 3 such has a high number of visitors seeking its recreational function, which it was unable to satisfy in its historical shape.This led to a radical reconstruction, although some trees were left from the original disposition.It contains a more formal part with aesthetic purpose centred around a fountain, shaded lanes with trees, free grassy areas and a playground for children.In our study, it represents historical greenery repurposed into a more modern form in late twentieth century.It has an area of 3.14 ha, and the GPS coordinates are 48° 08′ 59.4″ N 17° 07′ 10.4″ E.
Sad Janka Kráľa (SJK, Janko Kráľ Park, literally Orchard) is the oldest public park in Central Europe, founded in the 1770s.It is situated on the right shore of the Danube River in the Petržalka City part.Paths were cut in the original floodplain forest and planted by alleys in the shape of an eight-rayed star.Its current disposition dates back to a reconstruction starting in 1832 in the style of an English landscape park, showcasing several newly introduced species of trees and shrubs (Steinhübel 1990).Until present, it remains one of the most important recreational green areas in Bratislava.Its area is 23 ha, with GPS coordinates 48° 08′ 05.4″ N 17° 06′ 34.6″ E.
Grasalkovičova záhrada (GZ, Grassalkovich Garden) is a baroque garden belonging to a palace that is now the seat of the president of Slovak republic, known as the Presidential Palace.The palace was built in 1760 by a Hungarian aristocrat Anton Grassalkovich (Tomaško 2004).The current disposition of the garden mostly corresponds with the historical one, with a regular parterre typical for French gardens.It is located in the urbanised area of the city centre, but the representative function it still has prevented any major changes in its disposition towards a less formal one.It is open for public but surrounded by a high solid wall that can present a psychological barrier from entering it for casual recreation.The only exception from the historical shape are wooden playground elements placed in the garden.The last reconstruction of the park was in the 1990s.The area of the garden is 3.68 ha, and the GPS coordinates are 48° 09′ 01.3″ N, 17° 06′ 28.6″ E.

Dendrological inventory
Trees and shrubs in all parks were determined at the level of genus or species, and their phenophases were observed in the vegetation period of 2021 and 2022.In dioecious species, only male individuals were counted.Measurements were taken at the same time: crown diameter and height (m) for trees and height (m) and ground coverage (m 2 ) for shrubs.Using these data, we calculated the crown volume of each specimen using the formula for an ellipsoid (broad-leaved species) or a cone (conifers).We did not consider the volume of the grassy areas since the mowing regime is the same in all studied localities and aims to prevent their flowering and emission of pollen as much as possible.The nomenclature of taxa and syntaxa are based on WFO (2023).

Calculating allergenic potential
In I UGZA (Cariñanos et al. 2014), the allergenic potential is expressed on a scale from 0 to 1, where 1 represents the highest possible potential allergenic value (PAV max ) of a given area.The value of this index is a sum of the contributions of all species in the studied area and depends from their allergenic potential (ap), type of pollination (tp), the length of the pollen production period (ppp) as per Table 1 and the potential volume of an individual of each species (V), determined by its height and crown diameter at maturity multiplied by the number of individuals of the species in the studied area (n).It is compared to a potential area with maximal allergenic value and maximal height of vegetation cover (H max ) in the whole area (S T ) of the studied locality.The resulting formula is: I ISA uses the same method of calculation, but instead of the potential volume of an individual of each species uses actual measurements of the dendroflora in the studied area and determines its real volume.
Both indices depend on the maximum height of the vegetation in the study area, which might differ in various climatic conditions.We used the maximum height of 30 m, which corresponds to the height of several broad-leaved native species in ideal conditions.The parameters for determining the PAV in both indices are the type of pollination, duration of the pollination period and allergenic potential.The increase in any of the parameters by just 1 can double the resulting PAV, so the border cases need to be considered carefully.
The allergenic potential of each species depends on the frequency and severity of pollinosis caused by pollen grains of that species and is available in databases and literature sources.However, these sources are only applicable to a particular region, e.g.Cariñanos and Marinangeli (2021) for the Mediterranean area, since the concentration and repeated exposure to different allergens can influence the severity of pollinosis.For a species to be considered allergenic, it must produce a high quantity of airborne pollen containing allergenic molecules, but also be abundant in the region (Hrubiško 1998).We used data considering pollen-related allergic aggressiveness of local woody plant species (Jurko 1990a, b, c;Hrubiško 1998) to determine the allergenic potential of the air in the studied area.Only when not available did we use data from as geographically and climatically close regions as possible (Jochner-Oette et al. 2018;Kasprzyk et al. 2019).
We consider the principal pollination period parameter the most unreliable of the ones used to calculate PAV.The difference of just a few days can double the PAV here, while most sources only give the flowering period in months, e.g.March-April.The solution of this problem could be using palynological data from local sources.However, as a provider of such data, we noticed big variation in the duration of the pollination period between individual years (Ščevková et al. 2010).For example, the pollination season for the genus Populus lasted 60 days in 2008, but only 17 days in 2006.The timing of pollen release in the atmosphere depends strongly on the weather, especially on the temperature and humidity (Bartková-Ščevková 2003), and can drastically differ between years.Other problems with using aeropalynological data is the absence of species with biotic pollination from them, and the inability to distinguish between individual species of a genus or even a family or group of families like Cupressaceae-Taxaceae.
With regard to these limitations of determining the principal pollination period and low information value of singlenumber indices for the general public (they are useful for planning of urban vegetation but not for informing people with pollinosis about the time window when it is safe or not safe for them to visit a particular public park), we propose using a Running I UGZA (rI UGZA ) and Running I ISA (rI ISA ) instead.
rI UGZA and rI ISA are calculated similarly to I UGZA and I ISA , with omission of ppp from the formula.Instead, they are calculated separately for every month, using the value 1 instead of ppp if the species is flowering during that time and 0 if it is not.The overall allergenicity index for the studied location can be calculated as a sum of rI UGZA or rI ISA values for all months of the year.Therefore, our PAV value is calculated by a slightly different formula, using cmf (calendar months of flowering): To mark the distinction, we mark it with an apostrophe.Since the usual flowering period is between 1 and 3 months, the values are comparable with the original I UGZA and I ISA , which also uses values between 1 and 3 for the ppp parameter, therefore we mark them as I UGZA ′ and I ISA ′.They also correspond quite well to the original values since the flowering time of 4-6 weeks (ppp 2) means the plant flowers over the course of 2 calendar months.Unlike ppp, cmf is easier to find in guides to local flora (for native species) and ornamental plant catalogues (for introduced ornamental species) and does not have such a big variability between years as ppp.
The calculated values of I UGZA ′ and I ISA ′ were compared with various parameters of the parks (area, number of woody plants, their growth forms, density of trees and Shannon's diversity index; Shannon and Weaver 1949).The values of rI UGZA and rI ISA for each month were compared with the parameters changing in time (the number of flowering woody plants, trees, shrubs, allergenic and wind pollinated woody plants in each month as well as the number of their species) by Spearman non-parametric test.
Figure 3 shows the calculated values of the rI UGZA and rI ISA for the three studied public parks over the vegetation period, which in Bratislava lasts from February to October.However, all values were negligible before March and after July.The difference between these values is caused by the fact that rI UGZA takes into account the potential size of the woody plants at their maturity, while rI ISA works with their actual size at present.Therefore, rI UGZA is better suited for long-term planning of park management and revitalization, Table 2 The species of woody plants and their numbers identified at the locations MZ -Medic Garden, SJK -Janko Kráľ Park and GZ -Grassalkovich Garden in Bratislava with their type of pollina-tion (tp), allergenic potential (ap), flowering in calendar months and potential allergenic value (PAV′) valid for the area of Slovakia (Central Europe) An analysis of the greatest contributors to allergenic potential in the studied areas overall and in individual months are in Tables 3 and 4 and Figs. 4 and 5. Table 3 shows the taxa that would be top contributors to the longterm allergenic potential in the studied area when reaching their maximum potential size, which is a crucial information   for future planning, represented by I UGZA ′.Table 4 shows the actual allergenic potential in the studied area at present, represented by I ISA ′.This information is most valuable for people with pollen-related allergies.If we would consider all studied parks as one area, the greatest contributors to both I UGZA ′ and I ISA ′ are the genera Platanus, Acer and Tilia.
In MZ, the greatest contributors to the overall pollenrelated allergenicity (both I UGZA ′ and I ISA ′) are the plane trees (Platanus × hispanica) flowering from April to June (Table 4).There are 15 of them, with both tp and ap = 3.
In June, another big contributor to the actual allergenicity of the area (expressed by I ISA ′) is the amphiphilous species Tilia platyphyllos with tp = 2 and ap = 3, since it is the most numerous species in the area (107 trees).These trees form an alley consisting of a globose cultivar of this species which does not reach the full size of the wild form, which is reflected in I UGZA ′ where, despite its abundance, this particular cultivar is not among the principal contributors to its value.Thanks to its bigger potential size, which needs to be taken into account for future planning, this place belongs to In SJK, 16 specimens of Platanus × hispanica are also the greatest contributor to the allergenicity in the area from April to June (Table 4).163 specimens of most abundant genus Acer with tp = 2 and ap = 3 cause another significant increase of allergenicity in April and May.These two taxa are the biggest contributors to both I UGZA ′ and I ISA ′, since most of their individuals are fully grown.
In GZ, 19 of the 34 identified taxa flower in May, causing a peak of allergenicity in this month (Table 4).The greatest contributors to both I UGZA ′ and I ISA ′ are 64 specimens of Tilia × europaea (tp = 2, ap = 3) flowering in June and July and 36 specimens of Quercus robur with tp = 3 and ap = 4, flowering in May.
Figures 4 and 5 show the principal contributors to the allergenic potential during the individual months in absolute numbers, which enable comparing the values between the Fig. 5 The principal contributors to the actual allergenic potential in the studied parks: MZ -Medic Garden, SJK -Janko Kráľ Park and GZ -Grassalkovich Garden in individual months and their absolute contribution to the I ISA ′ value studied areas.The most noticeable peaks in the total contribution to I UGZA ′ are all in the GZ, caused by the flowering of Quercus robur in May and Tilia × europaea in June-July.However, most of these trees have been planted recently and will not reach their full size for several more decades.This is reflected in the absolute contribution to I ISA ′, where the greatest peak is caused by the simultaneous flowering of 107 specimen of Tilia platyphyllos in MZ in May, although the peaks caused by the above mentioned species in GZ are also noticeable.
Table 5 shows the different parameters of the parks that can influence the values of I UGZA ′ and I ISA ′.The composition of species does not differ much between the three studied areas: from the 84 identified taxa, only 23 are unique for a single area, with most of them ( 17) in SJK due to its much bigger size.The composition of growth forms is also comparable, with the ratio of shrubs to trees between 1:4.38 in GZ and 1:5.21 in SJK.However, there is one significant difference: the density of woody plants.The values of I ISA ′ best correspond with the density of trees, allergenic woody plants and wind-pollinated woody plants, while the values of I UGZA ′ show negative correlation with Shannon's biodiversity index.
Spearman's correlation between the values of rI UGZA , rI ISA and several woody plant-related parameters of the parks in the individual months (Table 6) shows a significant positive correlation of rI ISA with the number of flowering allergenic woody plants and their species in the given month, as well as with the number of tree and shrub species and density of flowering trees.Besides these, rI UGZA also has a significant positive correlation with the number of flowering woody plants and their species in the given month, as well as the number of wind pollinated woody plants and their species.The two indices are also significantly correlated with each other.

Discussion
Our results show a difference of one magnitude order in using I UGZA , which takes into account the potential size that each species of woody plants in the study area can achieve, compared to I ISA operating with the actual measurements of each specimen.The cause for this becomes evident when looking at the potential size of some tree species in the central European region.Several of the species can achieve height of 30 m and a 25-m diameter of the crown.However, if we look at the trees in urban parks, we see they achieve this size only rarely.The dendroflora in an urban environment has to deal with several stress factors (Wong et al. 2021), causing lowered life expectancy of trees.Even the planting plans of urban parks do not usually consider the potential size of trees under optimal conditions but plant

Table 5
Chosen parameters of the studied parks in Bratislava: MZ -Medic Garden, SJK -Janko Kráľ Park and GZ -Grassalkovich Garden (bold type represents the greatest value) them much closer to fill in the space in the early years of the tree's development.As the trees mature in the time horizon of decades, such dense planting plans lead to the necessity to cut some of them to give the others more space to develop.For this reason, the difference between I UGZA and I ISA is especially noticeable in recently planted areas.In GZ, where most of the trees were planted in the last 30 years, I UGZA is 16 times higher compared to I ISA .In SJK, it is 10 times higher and in MZ only 7 times.Jochner-Oette et al. ( 2018) calculated the values of I UGZA = 0.173 and I ISA = 0.018 in a public park in Eichstätt, Germany, showing a similar difference of one magnitude order between the two values.
When comparing the three studied localities, we see that SJK has the lowest value of both I UGZA ′ (0.41) and I ISA ′ (0.043).It is an English landscape park, with the highest area of the three studied parks, but lowest density of woody plants and lowest representation of highly allergenic species.The highest value of I UGZA ′ (0.985) was found in GZ, while the MZ has the highest value of I ISA ′ (0.091).We can conclude that the baroque garden, which is a more formal type of historical greenery, has a higher allergenic potential than the looser composition of an English landscape park.This is true even if the original composition of the area has been changed and modernised, like in MZ.The value of I UGZA ′ in GZ is actually very close to 1, which would mean the same allergenic potential as an equal area with a full vegetation cover by fully grown anemophilous trees with highly allergenic pollen.When looking at the area, this is obviously not true: we can see ornamental lawns and paths lined with trees of small size, mostly young, but some even cut to maintain the shape.If all of these trees were fully grown to their potential size, their volume might even exceed the available space in the garden.In such a case of recently planted or renovated parks or highly ornamental historical gardens where trees are cut to maintain small shape, I UGZA ′ can give a potential value that may never be achieved.It is a good tool for planning, but for comparing the actual allergenicity of the urban green spaces, I ISA ′ should be considered despite the time-consuming process of taking the actual measurements of trees.For large areas, it may be possible to use approximate size categories while creating the dendrological inventory.The difference between I UGZA and I ISA also shows that using the potential size of trees in ideal conditions can skew its value into unrealistic numbers and therefore it might be better to use an average size the particular species achieves at maturity in urban conditions, considering the space limitations of the particular area.Building a database of these measurements for different climatic regions would be useful to achieve higher precision in the allergenic potential index used for planning of urban greenery.
However, a single value is not enough for people suffering from pollinosis to decide if they should avoid a certain area or not, and neither is it enough for management of the Based on a high PAV value, we can identify species that should be avoided in public spaces in general, but even taxa with a moderate PAV value, like Acer and Tilia, can cause a peak in allergenicity if their abundance is high.rI UGZA makes it possible to plan the plantings and cuttings in the parks without causing such peaks of allergenicity in a particular month.Based on its values, we identified May as the month with highest allergenic potential in all studied parks, meaning that future plantings should avoid allergenic species that flower in this month and when cutting is needed, such trees should be considered before others.These include the species Platanus × hispanica, Celtis occidentalis and Acer sp. in MZ, P. × hispanica and Acer sp. in SJK and Quercus robur in GZ.However, in the case of SJK, we see the necessary individual approach to each park.The plane trees, contributing the most to the value of rI UGZA here, are also some of the most valuable dendrological objects in the park.They are 200-250 years old and the biggest one has a circumference of 752 cm (Zahradníková et al. 2020).For this reason, future plantings should avoid trees that flower in the same months, and if cutting is necessary, it should be focused on younger individuals or on Acer sp.
The values of rI ISA give information about the actual allergenic potential of the parks based on the actual measurements of the woody plant.Thanks to them, the public can make an informed decision about visiting parks and doing activities in them in the particular time of the year.The individual parks differ not only by the values of rI ISA , but also by the period of increased allergenic potential, making some more suitable for allergy sufferers than others.In MZ, the highest value of rI ISA (not just in this park but from all of the studied locations) was rescored in April, but the period of increased allergenic potential lasts until June, caused first by the flowering of Platanus × hispanica and later Tilia platyphyllos.This makes MZ the least suitable for outdoor activities for pollen-related allergy sufferers from the studied parks.The second highest peak in rI ISA was recorded in GZ in May, caused by the flowering of Quercus robur, which is one of the main allergens in the area (Hrubiško 1998).However, this is the only month of increased allergenic potential in this park, and so it is more advisable to visit it in the other parts of vegetation season, at least until the linden trees are fully grown and the values in June and July become higher, as we see by the values of rI UGZA .In SJK, rI ISA is constantly lowest from all studied parks, with the exception of April.
The period of increased allergenic potential lasts from April to May, so sensitive individuals should avoid this park in this period, but overall it is the best choice for outdoor activities during the year.Even in April and May, the particular allergenic species (again predominantly P. × hispanica) can be avoided due to the big size of the park.
As has been mentioned before, the pollination period can differ between years and so can the beginning and end of the phenophases, with an observed shift towards earlier timing with the ongoing climatic change (Cho et al. 2017), so the running allergenicity indices may not be fully accurate each year, but still offer a good approximation of the actual allergenicity for planning a visit to a particular location.However, it is important to note that the actual levels or airborne allergenic particles also depend on the weather, with highest values during a warm, sunny weather.Increased humidity can lead to lower release of pollen from anthers and faster sedimentation of pollen grains, but also to rupture of pollen caused by osmotic shock and releasing of small respirable allergenic molecules in the air (D 'Amato et al. 2021).
Besides the pollen of woody plants, the pollen of grasses (Poaceae family) and fungal spores are the greatest contributors to the allergenicity in urban parks.Ornamental flowers usually grown in the parks are mostly entomophilic, and so their contribution to the allergenic potential is negligible.Grasses are considered one of the main aeroallergens, with tp = 3 and ap = 4 (Sabariego et al. 2021), and with flowering from May to August can considerably increase the allergenic potential of an area.However, the intensive mowing regime in public parks prevents flowering in most of the grasses that are a part of the lawns there.As such, the lawn maintenance purposedly aims to decrease the risk of pollinosis during the pollen season of grasses.Fungal spores are also abundant in urban green areas, especially those of the genus Cladosporium and Alternaria, reaching their peak in July-August (Kasprzyk et al. 2021).However, it is not possible to detect these without aerobiological sampling, which may not be available to the management of these areas.This is why we assume the contribution of fungal spores to be approximately equal in all parks and relatively low thanks to the immediate removal of dead plant matter (moved grass, woody plant cuttings and leaves, old declining trees) which is the primary source for fungal growth here.However, all of these factors would be necessary to consider in urban green areas with less intensive management.
Similarly to our results, Cariñanos et al. (2017) found a significant positive correlation between I UGZA and the number of trees/ha and a negative correlation with Shannon's biodiversity index.From the two indices we considered, we found I ISA to be better correlated with the density of trees and I UGZA with Shannon's index (negatively).Regarding the running indices, rI UGZA was better correlated with the parameters taking into account the number of flowering woody plants of particular type in the given month.However, these correlations can differ in various parks depending on the percentage of juvenile trees and site-specific growth conditions which may not allow some of the trees to reach their potential size even in maturity.The more the size of the trees in park nears their potential size, the closer these two indices will be.
Allergenicity indices have been used and tested in several regions of Europe, especially in the Mediterranean.Like with I UGZA and I ISA , the running indices are applicable for determining the allergenic potential of urban green zones in other cities.However, when comparing the potential allergenicity of green spaces in different regions, it is necessary to consider local data in determining the PAV or every plant species and adjust the values in all indices accordingly, as well as adjust the maximum potential height of the vegetation in I UGZA .Aerobiological data and pollen calendars can help with determining whether a species is anemophilic or entomophilic, local phenological data give information about the time of pollination and allergological data are needed to determine the allergenicity of the species.This is especially important because the allergenicity of some species can differ in various regions based on their abundance, pollen production and sensibilisation of the population (Hrubiško 1998).When these local values are taken into account, the running indices provide an accurate information about the temporal variation of the allergenic potential in any urban green area.

Conclusions
The allergenic potential of public green spaces varies during the vegetation period, and information about this variation is important for people suffering from pollinosis making a decision about their outdoor activities, as well as for the planning and management of green public spaces.Using a modified I UGZA and I ISA with calendar months of flowering instead of the duration of the pollination period, we created a running index, showing the variation of allergenic potential during the vegetation season in three important public parks in Bratislava.It enables to see the temporal of potential allergenicity during the year and identify the period of increased allergenic potential.The running indices are applicable and comparable for green areas in different locations, on the condition that the parameters determining the PAV of each species are checked and adjusted for the local conditions.Based on the identified woody plant species that most contribute to it, targeted measures can be implemented to reduce or avoid allergenicity within these areas.These measures may include selective plant removal, the introduction of low-allergenic species and strategic park design to minimise the dispersion of allergenic pollen, as well as avoidance of individuals sensitive to pollen of certain areas in particular months.

Fig. 4
Fig.4The principal contributors to the long-term allergenic potential in the studied parks: MZ -Medic Garden, SJK -Janko Kráľ Park and GZ -Grassalkovich Garden in individual months and their absolute contribution to the I IUGZA ′ value

Table 1 The parameters of I UGZA according to Cariñanos et al. (2014)
ISA provides more relevant information about the actual situation during the vegetation season for pollenallergy sufferers.The data show the main allergenic season in these parks lasting from March to July.For the rest of the year, the values are negligible for both rI UGZA and rI ISA .The highest value of rI UGZA (0.362) was noted in May in GZ and the highest value of rI ISA (0.029) in April in MZ.The sum of the monthly indices gives the value of I UGZA ′ and I ISA ′ for these locations.For MZ, I UGZA ′ is 0.67 and I ISA ′ 0.091.SJK has I UGZA ′ 0.41 and I ISA ′ 0.043 and GZ has I UGZA ′ 0.985 and I ISA ′ 0.061.The graph also gives information about the duration of the period of increased allergenic potential in the individual parks, defined as the period when the value of rI UGZA /rI ISA reaches over 25% of the total I UGZA ′/I ISA ′ over the whole vegetation period.Taking into account the value of rI UGZA , this period is April-May in both MZ and SJK and May-June in GZ.For rI ISA , it lasts from April to June in MZ, April-May in SJK and May in GZ.

Table 6
Spearman's correlations between rI UGZA , rI ISA , and several woody plant-related parameters of the parks in the individual months (March-July) green areas.A situation is possible where the allergenicity in a certain area is exceptionally high during 1 or 2 months, but low values during the rest of the year push the potential allergenicity index lower and show no problem with the current composition of allergenic species.For this reason, we find it more informative to calculate a running allergenicity index for each month.