From:  Climate-driven changes in pollen dynamics: increased loads and earlier, longer exposure

 Temporal trends in pollination of the Cupressaceae family.

Reference (publication date)Region (number of sites)Study periodMethodologyTrend identification techniquePIn (grains per cubic meter of air per decade and significance)PSS (days per decade and significance)PSD (days per decade and significance)
[20] (2010)NW Italy (1)1981–2007Hirst-typeLinear regression+ (*)– (3.3)+ (6.7)
[26] (2003)Switzerland (1)1979–1999Hirst-typeLinear regression+ (*)– (12.9; *)No trend
[27] (2025)Texas, USA (1)2009–2023Hirst-typeLeast square regression+ (20.5; *) for Juniperus
[36] (2016)Spain (12)1994–2013Hirst-typeLinear regression+ at 7 sites out of 12
[37] (2021)Basel, Switzerland (1)1969–2018Hirst-typeLinear regression+ (~2,050; ***)– (~5.5; **)+ (~4.4; **)
[47] (2021)Bavaria (8)1989–2018Hirst-typeLinear regression+
[51] (2024)NE Italy (20)2001–2022Hirst-typeTheil-Sen estimator+ (916)– (8.8)+ (8)
[56] (2012)Europe (97)1977–2009 (partially)Hirst-typeLinear regression+ (***)
[57] (2020)N & Central Italy (9)2000–2016Hirst-typeNon-parametric Mann-Kendall testNo trendNo trend
[62] (2025)N Italy (1)1997–2023Hirst-typeLinear regression+ ns– ns+ ns
[73] (2002)W of the NL (1)1969–2000Hirst-typeLinear regression– (6.25, **) for Juniperus
2000s–2090sModellingUnspecified– (0.8)
[74] (2025)Iberian Peninsula (12)1994–2023Hirst-typeLinear regression+ (1,745 on average, but up to 8,535 in Granada)Slightly + (1.6)
[81] (2007)Thessaloniki (1)1987–2005Hirst-typeLinear regression+ (~1,435)No trendNo trend
[82] (2024)Rome (1)1997–2016Hirst-typeLinear regressionNo trend
[86] (2020)Andalusia, Spain (3)1994–2017Hirst-typeLinear regression– in Málaga, + in Granada (***), ns elsewhere– in Granada (***) and Málaga (*), ns in Córdoba
[88] (2021)Milan, Italy1995–2017Hirst-typeUnspecified+ (2,100, especially after 2001)
[89] (2024)Madrid region (10)1994–2022Hirst-typeLinear regression+ (especially after 2012)
[90] (2024)N Italy (1)1989–2018Hirst-typeLinear regression + Mann-Kendall test+ (***)+ (**)
[91] (2016)S Spain (1)1996–2010Hirst-typeLinear regression + local regression method (LOESS)+ (5,670; ***)
[92] (2021)Oklahoma, USA (1)1996–2020Hirst-typeLinear regression+ (~2,200; ***)– (ns)+ (***)
[93] (2023)Andalusia (1)1998–2020Hirst-typeLinear regression + Mann-Kendall test + Sen’s slope estimation+ (19; *)– (1.3; ns)
[94] (2018)S Spain (1)1994–2016Hirst-typeZ coefficient estimation + Mann-Kendall test+ (~2,000; **)– (ns)+ (~17; ***)
[95] (2022)Madrid (1)1979–2018Hirst-typeLinear regression+ (~3,000 gr/m3 per 1°C annual mean increases in temperature; *)– (~7.75)+ (~3.25; ns)
[96] (2025)SW Spain (1)1993–2022Hirst-typeLinear regression + Mann-Kendall test + Theil-Sen estimatorNo trendNo trendNo trend
[97] (2016)SW Spain (1)1993–2013Hirst-typeLinear regression+ (ns)
[98] (2007)Central Italy (5)2006–2036 vs. 1990–2005Phenological model    – (4.2)– (0.8)
2037–2068 vs. 2006–2036Phenological model– (4.2)– (0.1)
2069–2100 vs. 2037–2068Phenological model– (4.7)– (0.1)

The grey backgrounds correspond to studies that have made projections about the future. +: trend towards increased pollen concentration (PIn), earlier onset of pollination (PSS), or longer pollen season (PSD); –: trend towards decreased pollen concentration (PIn), later onset of pollination (PSS), or shorter pollen season (PSD); *: p < 0.05; **: p < 0.01; ***: p < 0.001. PIn: pollen integral; PSS: pollen season start date; PSD: pollen season duration; NW: northwest; NE: northeast; N: north; ns: non-significant; W: west; NL: Netherlands; S: south; LOESS: locally estimated scatterplot smoothing; SW: southwest.