3.2. Aerobiological study
Regarding the total annual pollen records, no clear pattern was observed for all the sampling years (Fig. 3). In 2017, a greater annual total integral was obtained in the periphery, while in the other years the opposite occurred. A longer time series would be necessary to exclude the particularities of each year and clearly discern if a sampling point generally registers higher integrals than the other. However, the directionality of the interannual changes in these values was the same at both sampling points (i.e. they increased or decreased at the same time at the two sampling points).
The 11 main pollen types constituted, as a whole, the 89.43% of the total pollen collected during the study period. However, despite these pollen types were the same during the 3 years at both sampling points (with the exception of Casuarina), slight differences were detected in the order of abundance from one year to another. Furthermore, the order of abundance of the different pollen types changed depending on the sampling location. In general, almost all pollen types had a higher relative abundance in the city center, compared to the periphery (Fig. 4).
As previous studies highlighted, the abundance of the pollen types detected at certain city area is highly conditioned by the abundance of the emission sources in the surroundings, by ornamental trees and local meteorological conditions (Charalampopoulos et al., 2021; Fuentes Antón et al., 2020a; Pecero-Casimiro et al., 2019; Velasco-Jiménez et al., 2013; Werchan et al., 2018). This may explain the differences observed in the order of abundance of the pollen types detected (Fig. 4). Besides that, the relative abundances are also influenced by the amounts of total pollen detected in each sampling location, so the analysis of these results must be done jointly with those of the annual integrals (Fig. 5).
As can be seen in Fig. 5, the values of the annual pollen integrals showed differences from year to year, depending on the sampling station and the pollen type. Parietaria stands out, always presenting much higher values in the city center, probably due to the greater abundance of ruined buildings, which favors the proliferation of nitrophilous herbaceous plants, as is the case of Parietaria that usually grows in abandoned lots and old walls, which explains not only its higher annual integrals in the center of the city, but also the greater relative abundance (Figs. 4 and 5). On the contrary, Amaranthaceae, Casuarina and Plantago always presented higher values in the periphery. Amaranthaceae and Plantago are abundant plants in agro-forestry and crop transition zones that are more abundant in the outskirts of the city (Fig. 2), which may explain their higher integrals detected at this sampling point. On the other hand, Casuarina is a common ornamental tree on the Teatinos university campus, where the pollen trap is installed, which would have given rise to the highest integrals recorded in the periphery. Previous studies also found similar patterns in the distribution of the airborne pollen concentrations of these pollen types inside cities due to differences in land uses and in the distribution of ornamental trees (Fuentes Antón et al., 2020a; Velasco-Jiménez et al., 2013).
Regarding the remaining pollen types, they reached higher annual pollen integrals in one or another sampling stations depending on the year and did not present patterns as clear as in the case of the mentioned above, despite the fact that some of them presented important differences in terms of relative abundance. This may be explained by a similar abundance of the emission sources of these pollen types at both sampling sites, and wind dynamics during each year being the determinants of the airborne pollen concentrations detected (Fuentes Antón et al., 2020a).
In the case of Cupressaceae, the annual integrals presented scarce differences between the sampling points, being slightly higher in the city center in the years 2017 and 2019, and slightly higher in the periphery in 2018. This similarity is probably due to the fact that cypresses are common ornamental plants in the vicinity of both sampling sites, and they have been previously reported to influence the local airborne pollen concentrations (Charalampopoulos et al., 2021).
The annual integrals of Olea europaea were similar in both sampling points in 2018 and 2019, but much higher in 2017 in the periphery. Given that there are no great differences in terms of the area occupied by olive groves in both sampling sites and that the percentage of olive grove use increases with the distance in a similar way, the observed differences could be caused by punctual variations in the meteorological conditions, especially wind dynamics. It has been observed that when the northwest wind blows, coming from inland, the levels of olive pollen increase, being the peripheral area of the city the most exposed to such winds. On the contrary, winds blowing from the sea (southeast), tent to decrease pollen concentrations, being the pollen trap situated in the city center more exposed to these winds. Something similar happened with Quercus pollen type: its annual integral was higher on the periphery in 2017, similar in 2018 at both sampling points, and higher in the center in 2019. The forests that contribute to increase this pollen type are found in similar abundances around both pollen collectors (Fig. 2), so the differences detected among the years are probably also conditioned by wind dynamics.
In the case of the Pinus, the annual integrals were similar in both locations in 2017, higher in the periphery in 2018 and higher in the center in 2019. Despite their proximity, the repopulation pine forests located in the Montes de Malaga, towards the northwest of the city do not seem to have had a great influence on the quantities detected in the city center pollen trap. This may be due to the fact that the prevailing winds in Malaga are those blowing from SE and NW, and they would not contribute to the transport of pollen from these pine forests to the sampling sites (Recio et al., 2018). In this case, the main source of pollen would be small copses located within the parks and peri-urban areas, frequent in the vicinity of both sampling sites.
The annual integral of Platanus was higher in 2017 in the city center but in 2018 and 2019 it was higher in the periphery. The integral of Poaceae was practically similar in the two sampling points throughout the period studied, and the integral of U. membranacea was higher in the center in 2017, higher in the periphery in 2018 and similar in both sampling sites in 2019. Platanus is represented in the area by ornamental trees whose abundance is similar in both sampling points. Grasses, due to their wide distribution and high number of species, are present in a large number of different land uses, therefore, its pollen concentrations in the air tends to be similar at both locations. Something similar happens in the case of U. membranacea, as it is a nitrophilous ruderal plant, it does not present a special predominance in one or another sampling point, apart from the influence that atmospheric conditions and wind trajectories may exert.
Although several studies have evidenced the influence that different heights above ground level could have on the pollen concentrations detected in the atmosphere, these differences tend to be insignificant when the collectors are located more than 10 m above ground level (Rojo et al., 2019a, 2020). Therefore, we consider that the influence of the scarce difference in the heights of the pollen traps used in this study, has been of little relevance for the results obtained. On the other hand, in other localities, land uses and pollen emission sources located in the vicinity of cities have also been identified as the main variables responsible for the differences detected between nearby sampling points (Fuentes Antón et al., 2020b; Pecero-Casimiro et al., 2019; Picornell et al., 2019; Rojo et al., 2020; Velasco-Jiménez et al., 2013).
At both sampling points, about 85% of the total annual pollen is concentrated from February to June (both inclusive), the period in which the highest daily pollen concentrations of most of the main pollen types took place, with the exception of Casuarina, which is typically autumnal (Figs. 6 and 7). In general, May is the month in which the highest concentrations of the total pollen are detected at both sampling points, a typical behavior in the Mediterranean basin. Although almost all the pollen types presented a single period of maximum concentrations, some of them such as Cupressaceae or Amaranthaceae, presented two peaks of different intensity. In the case of Cupressaceae, two periods of high concentrations were detected, the first, more intense, from January to May and, the second, of lesser intensity, from September to December, due to the presence of different species, some of them with autumn flowering. Something similar occurs in the case of Amaranthaceae, whose maximum concentrations occur between April and May and, after a marked decrease, they slightly increase again at the end of summer (August-September). In any case, sporadic detection of pollen types is frequent throughout the year, due fundamentally to re-buoyancy phenomena.
As previously observed for relative abundances and annual integrals, mean daily concentrations, in general, were also higher in the periphery for the pollen types Amaranthaceae, Casuarina, and Plantago, and higher in the city center for Parietaria (Figs. 6 and 7). This last pollen type showed the greatest differences between the two sampling points, especially during the months of March to July. The rest of the pollen types presented similar daily average concentrations throughout the year at both sampling points, despite the fact that some peaks of greater intensity were sporadically detected at one or another. In general, the timing and intensity of the daily average concentrations detected were similar to those of other nearby cities located in the south of the Iberian Peninsula such as Granada, Ronda or Cordoba (De Linares et al., 2019; Picornell et al., 2019; Velasco-Jiménez et al., 2018).
If we compare the graphs of seasonal behavior of the different pollen types in both sampling points, we observe that there is a great coincidence in terms of their profiles. In fact, the correlation studies carried out between the mean daily concentrations of the different pollen types at both sampling points (center and periphery), resulted in positive and significant correlation coefficients for p ≤ 0.001 in all cases. (Table 1). However, after applying Mann-Whitney-Wilcoxon tests, significant differences (p < 0.05) were observed in 6 of the 11 pollen types studied (i.e. Amaranthaceae, Casuarina, Olea europea, Plantago, Quercus, and Parietaria).
Table 1
Results of the statistical analysis (Spearman correlation and Wilkoxon tests) between the values of the daily mean pollen concentrations of the two sampling sites during the main pollen season (MPS). Significant p-values are marked with an asterisk (p ≤ 0.05).
| | Spearman Correlation | Wilcoxon test |
| Pollen types | Spearman coefficient | P-value | P-value |
| Amaranthaceae | | 0.536 | < 0.001* | < 0.001* |
| Casuarina | | 0.536 | < 0.001* | < 0.001* |
| Cupressaceae | 0.767 | < 0.001* | 0.477 |
| Olea europaea | 0.863 | < 0.001* | 0.006* |
| Parietaria | 0.539 | < 0.001* | < 0.001* |
| Pinus | 0.704 | < 0.001* | 0.652 |
| Plantago | 0.619 | < 0.001* | < 0.001* |
| Platanus | 0.628 | < 0.001* | 0.862 |
| Poaceae | 0.678 | < 0.001* | 0.152 |
| Quercus | 0.890 | < 0.001* | < 0.001* |
| U. membranacea | 0.487 | < 0.001* | 0.060 |
Given the proximity of both sampling points, we can expect that the meteorological conditions that influence pollen concentrations in the atmosphere are similar and therefore, the concentrations oscillate in the same direction, even when the absolute values be different due to differences in emission sources or air flows. The significant differences observed according to the Mann-Whitney-Wilcoxon tests were probably due to the fact that, repeatedly, concentration peaks detected in one of the sampling points were not detected in the other, and vice versa.
Regarding the MPS of the pollen types studied (Table 2), differences between both sampling stations have been observed in the dates of start and end, in the duration of the whole MPS and the periods pre- and post-peak, as well as in the dates and intensity of the peak days. Actually, few coincidences in dates have been found, and sometimes these differences are of up to 30 days between one and another sampling sites. It is possible that differences in the beginning, end and duration of the MPS are related to the methodology used, which is based on percentages and, therefore, it is dependent on the value reached by the annual integral of the pollen type in question, but this would not be the case of the peak day or their intensity and they have also been variable.
Some pollen types present their maximum daily peaks earlier in the periphery, as occurred in the case of Amaranthaceae, Casuarina, Olea europaea and Quercus. On the contrary, Parietaria always presented its maximum peaks earlier in the center. The rest of the pollen types did not present a clear behavior regarding the date in which the peak day take place, occurring first in the center or in the periphery depending on the year. For some pollen types, i.e., Casuarina, Plantago and Quercus, the differences are less than 10 days. However, for the rest of the pollen types, the fluctuations were greater. For example, in Platanus the fluctuations oscillated between 1 and 19 days; in U. membranacea between 6–24 days; and in Parietaria, between 10–80 days (Table 2).
Table 2
Phenological parameters related to the main pollen season (MPS) for the main pollen types during the studied years in the city center and the periphery sampling sites.
| Year | Station | Start date | Final date | MPS Duration (days) | Pre-peak Duration (days) | Post-peak Duration (days) | Peak date | Peak day concentration (p.g./m3) |
| Amaranthaceae |
| 2017 | Center | 16/02 | 29/11 | 287 | 76 | 211 | 02/05 | 14 |
| Outskirts | 21/02 | 09/10 | 231 | 54 | 177 | 15/04 | 59 |
| 2018 | Center | 14/02 | 09/10 | 238 | 79 | 114 | 03/05 | 21 |
| Outskirts | 03/04 | 02/10 | 183 | 25 | 158 | 27/04 | 61 |
| 2019 | Center | 30/01 | 30/10 | 274 | 100 | 174 | 09/05 | 16 |
| Outskirts | 21/02 | 19/10 | 241 | 54 | 187 | 15/04 | 17 |
| Casuarina |
| 2017–2018 | Center | 15/09/17 | 20/12/17 | 97 | 57 | 40 | 10/11/ 17 | 22 |
| Outskirts | 26/10/17 | 16/01/18 | 83 | 11 | 72 | 10/11/ 17 | 43 |
| 2018–2019 | Center | 24/10/18 | 25/12/18 | 63 | 1 | 62 | 24/10/ 18 | 17 |
| Outskirts | 27/10/18 | 02/12/18 | 37 | 10 | 27 | 05/11/ 18 | 63 |
| Cupressaceae |
| 2017–2018 | Center | 14/10/17 | 24/04/18 | 193 | 131 | 62 | 21/02/18 | 427 |
| Outskirts | 18/10/17 | 17/04/18 | 182 | 161 | 21 | 27/03/18 | 449 |
| 2018–2019 | Center | 22/09/18 | 03/05/19 | 224 | 173 | 51 | 13/03/19 | 212 |
| Outskirts | 08/10/18 | 23/04/19 | 198 | 145 | 53 | 01/03/19 | 273 |
| Olea europaea |
| 2017 | Center | 11/04 | 31/08 | 143 | 27 | 116 | 07/05 | 282 |
| Outskirts | 29/03 | 11/06 | 75 | 29 | 46 | 26/04 | 612 |
| 2018 | Center | 03/05 | 10/07 | 69 | 32 | 37 | 03/06 | 1048 |
| Outskirts | 01/05 | 06/07 | 67 | 35 | 32 | 04/06 | 1047 |
| 2019 | Center | 05/04 | 05/07 | 92 | 43 | 49 | 17/05 | 943 |
| Outskirts | 03/04 | 30/06 | 89 | 45 | 44 | 17/05 | 893 |
| Parietaria |
| 2017 | Center | 04/02 | 01/12 | 301 | 46 | 255 | 21/03 | 94 |
| Outskirts | 22/01 | 23/11 | 306 | 45 | 261 | 07/03 | 35 |
| 2018 | Center | 07/02 | 04/10 | 240 | 92 | 148 | 09/05 | 212 |
| Outskirts | 21/01 | 17/12 | 331 | 28 | 303 | 17/02 | 32 |
| 2019 | Center | 27/01 | 12/10 | 196 | 103 | 156 | 09/05 | 207 |
| Outskirts | 08/01 | 14/10 | 280 | 34 | 246 | 10/02 | 29 |
| Pinus |
| 2017 | Center | 17/02 | 28/10 | 254 | 21 | 233 | 09/03 | 15 |
| Outskirts | 18/02 | 02/09 | 197 | 10 | 187 | 27/02 | 60 |
| 2018 | Center | 24/02 | 05/08 | 163 | 22 | 141 | 17/03 | 32 |
| Outskirts | 26/02 | 01/08 | 157 | 8 | 149 | 05/03 | 50 |
| 2019 | Center | 08/02 | 21/07 | 164 | 21 | 143 | 28/02 | 74 |
| Outskirts | 13/02 | 19/07 | 157 | 17 | 140 | 01/03 | 90 |
Table 3
Parameters related to the main pollen season for the main pollen types during the period 2017–2019 in the city center and the periphery sampling sites.
| Year | Station | Start date | Final date | MPS Duration (days) | Pre-peak Duration (days) | Post-peak Duration (days) | Peak date | Peak day concentration (g.p./m3) |
| Plantago |
| 2017 | Center | 06/03 | 10/08 | 158 | 41 | 117 | 15/04 | 21 |
| Outskirts | 13/03 | 27/07 | 137 | 34 | 103 | 15/04 | 55 |
| 2018 | Center | 21/02 | 11/09 | 203 | 54 | 149 | 15/04 | 19 |
| Outskirts | 28/03 | 30/06 | 95 | 46 | 49 | 12/05 | 41 |
| 2019 | Center | 08/03 | 24/10 | 231 | 59 | 172 | 05/05 | 17 |
| Outskirts | 13/03 | 16/09 | 188 | 57 | 131 | 08/05 | 37 |
| Platanus |
| 2017 | Center | 07/03 | 12/04 | 37 | 4 | 33 | 10/03 | 140 |
| Outskirts | 06/03 | 03/05 | 59 | 6 | 53 | 11/03 | 117 |
| 2018 | Center | 07/03 | 14/04 | 39 | 4 | 35 | 10/03 | 85 |
| Outskirts | 10/03 | 12/05 | 64 | 20 | 44 | 29/03 | 90 |
| 2019 | Center | 05/03 | 18/04 | 45 | 22 | 23 | 26/03 | 48 |
| Outskirts | 04/03 | 10/04 | 38 | 8 | 30 | 11/03 | 102 |
| Poaceae |
| 2017 | Center | 27/02 | 12/11 | 259 | 81 | 178 | 18/05 | 24 |
| Outskirts | 18/02 | 02/10 | 227 | 68 | 159 | 26/04 | 66 |
| 2018 | Center | 03/04 | 02/08 | 122 | 41 | 81 | 13/05 | 71 |
| Outskirts | 02/04 | 14/08 | 135 | 58 | 90 | 29/05 | 91 |
| 2019 | Center | 04/03 | 10/09 | 191 | 75 | 116 | 17/05 | 144 |
| Outskirts | 11/03 | 05/09 | 179 | 72 | 107 | 21/05 | 143 |
| Quercus |
| 2017 | Center | 20/03 | 26/09 | 191 | 27 | 164 | 15/04 | 79 |
| Outskirts | 13/03 | 09/08 | 150 | 34 | 116 | 15/04 | 330 |
| 2018 | Center | 31/03 | 28/06 | 90 | 29 | 61 | 28/04 | 336 |
| Outskirts | 07/04 | 06/07 | 91 | 22 | 69 | 28/04 | 231 |
| 2019 | Center | 02/04 | 20/07 | 110 | 13 | 97 | 14/04 | 1022 |
| Outskirts | 02/04 | 29/08 | 150 | 3 | 147 | 04/04 | 448 |
| Urtica membranacea |
| 2017 | Center | 04/02 | 25/04 | 81 | 24 | 57 | 27/02 | 43 |
| Outskirts | 04/02 | 18/05 | 104 | 36 | 68 | 11/03 | 33 |
| 2018 | Center | 24/01 | 26/07 | 184 | 43 | 141 | 07/03 | 29 |
| Outskirts | 16/02 | 02/06 | 107 | 47 | 60 | 03/04 | 18 |
| 2019 | Center | 30/01 | 13/06 | 135 | 37 | 98 | 07/03 | 34 |
| Outskirts | 27/01 | 29/05 | 123 | 34 | 89 | 01/03 | 28 |
The peak day values of the pollen types Urtica membranacea and Parietaria were higher during the entire study period in the city center, while Amaranthaceae, Casuarina, Pinus and Plantago obtained higher peak concentrations in the periphery, which seems to be related to the higher abundance of these species in the vicinity of pollen sampler. The remaining pollen types did not show a homogeneous pattern, so the values of their peaks seem to be more conditioned by wind dynamics than by any specific land use pattern.
Previous studies also reported differences in the start, peak and end dates of the MPS when sampling at different locations inside the same city (Charalampopoulos et al., 2021; Fuentes Antón et al., 2020b; Velasco-Jiménez et al., 2013). In accordance with our results, these differences were related to the distribution of the pollen emission sources and wind dynamics.
The study carried out showed differences in the results obtained at both sampling points, not only in timing but also in the values reached by the different pollen types and their relative percentages. The different land uses of the surrounding places in which the two samplers are located, including parks and gardens and proximity to the sea, seem to have influenced the results obtained. In fact, the influence that the vegetation close to the pollen traps has on pollen concentrations already has been highlighted by other authors such as Kasprzyk et al. (2019) or Rojo et al. (2020). However, the influence of the prevailing winds or the structure of the city itself, could explain some of the differences found since the position of the buildings constitutes barriers that can modify the local air flow and affect pollen dynamics (Damialis et al., 2005; Gonzalo-Garijo et al., 2006; Peel et al., 2013; Ciani et al., 2020).