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

 Temporal trends in pollination of the Oleaceae 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)
[18] (2014)Europe (13)1990–2009Hirst-typeLinear regression+ in 7 sites of 10 (ns), + in Madrid (**)– or + (ns)– in 7 sites of 10 (ns, except Derby **)
[51] (2024)NE Italy (20)2001–2022Hirst-typeTheil-Sen estimator+ (646)– (13.5)+ (13.3)
Olea
[20] (2010)NW Italy (1)1981–2007Hirst-typeLinear regression+ (*)– (17.0; *)+ (6.7)
[36] (2016)Spain (12)1994–2013Hirst-typeLinear regression+ at 6 sites out of 9 (~12,900 in Córdoba)
[52] (2021)France (2)1995–2019Hirst-typeVariable+ (~350) at Aix-en-Provence; no trend at Toulon
[56] (2012)Europe (97)1977–2009 (partially)Hirst-typeLinear regression+ (ns)
[74] (2025)Iberian Peninsula (12)1994–2023Hirst-typeLinear regression+ (3,778 on average, but up to 11,940 in Jaén)Slightly + (1.6)
[82] (2024)Rome (1)1997–2016Hirst-typeReduced major axis linear regression+ (ns)
[83] (2002)Central Italy (1)1982–2001Hirst-typeLinear regression– (4.0; *)– (5.0; *)
[89] (2024)Madrid region (10)1994–2022Hirst-typeLinear regression+ (at least locally)
[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)+ (~190,000; ***)
[93] (2023)Andalusia (1)1998–2020Hirst-typeLinear regression + Mann-Kendall test + Sen’s slope estimation+ (3,830; ns)– (0.8; ns)
[94] (2018)S Spain (1)1994–2016Hirst-typeZ coefficient estimation + Mann-Kendall test+ (~14,500; **)– (ns)+ (~13.0; *)
[95] (2022)Madrid (1)1979–2018Hirst-typeLinear regressionNo trend– (~1.75)
[96] (2025)SW Spain (1)1993–2022Hirst-typeLinear regression + Mann-Kendall test + Theil-Sen estimator+ (1,668 to 1,808; **)– (1.6; ns)– (1.9; **)
[99] (2014)S Spain (1)1982–2011Hirst-typeLinear regression + seasonal-trend decomposition procedure based on LOESS + ARIMA model+ (~7,000; ***)– (~3.6; ns)+ (~11.9; ***)
[100] (2011)NW Italy (1)1981–2010Hirst-typeLinear regression+ (***)No trendNo trend
[102] (2013)Spain (4) & Italy (6)1992–2011Hirst-typeUnspecified– (especially 2004–2011, and in Italy more than in Spain)
2046–2065Thermal time model
[103] (2005)Andalusia (5)1982–2001Hirst-type– (4; *), but + in Cordoba (4)
2000s–2099Growing degree-days– (0.3 to 2.4)
[104] (2024)Barcelona (1)2000–2019Hirst-type+ (1.2; ns)
2024–2100Forecast model of the PSS + climate scenarios– (depending on the climate scenario)
[105] (2015)Spain (4) & Italy (6)2081–2100 vs. 1993–2011Pheno-meteorological models– (~0.9 to 1.1)
[106] (2012)Calabria (3)21st centuryPhenological model– (1.0 to 3.4)
[107] (2013)Central and S Italy (16)2050Model– (2 to 19, mean 8, with 1°C of warming)
[108] (2000)W Mediterranean1990s–2030sPhenological model– (0.6 to 1.6)
Fraxinus
[24] (2005)Austria in valley area (1)1980–2001Hirst-typeLinear regression+ (~950)– (3.2)
[25] (2018)Brussels (1)1982–2015Hirst-typeLocal regression method (LOESS)+
[26] (2003)Switzerland (1)1979–1999Hirst-typeLinear regressionNo trend– (14.3; *)No trend
[27] (2025)Texas, USA (1)2009–2023Hirst-typeLeast square regression+ (ns)
[28] (2021)Benelux (5)1981–2020Hirst-typeLinear regression+ in all sites; * in 4 of 5– in 4 sites of 5; * in only 1+ (* at just 1 site out of 5); – (* at another site)
[36] (2016)Spain (12)1994–2013Hirst-typeLinear regression+ at 7 sites out of 8; – in Madrid (ns)
[37] (2021)Basel (1)1969–2018Hirst-typeLinear regression+ (~588; **)– (~3.7; **)– (~3; **)
[38] (2021)Switzerland (14)1990–2020Hirst-typeLinear regression + LOESS+ (ns)No trend– (*)
[40] (2018)Brussels (1)1982–2015Hirst-typeLinear regression + Mann-Kendall test+ (*)– (ns)
[45] (2011)S Poland (1)1991–2008Hirst-typeLinear regression+ (*)No trend
[47] (2021)Bavaria (8)1989–2018Hirst-typeLinear regression+
[56] (2012)Europe (97)1977–2009 (partially)Hirst-typeLinear regression+ (***)
[73] (2002)W of the NL (1)1969–2000Hirst-typeLinear regression– (1, ns)
2000s–2090sModellingUnspecified– (1.9)
[82] (2024)Rome (1)1997–2016Hirst-typeReduced major axis linear regression+ (ns)
[86] (2020)Andalusia (3)1994–2017Hirst-typeLinear regression+ (especially in Granada; ns)– (* in Córdoba, ns elsewhere)
[89] (2024)Madrid region (10)1994–2022Hirst-typeLinear regression+ (at least locally)
[94] (2018)S Spain (1)1994–2016Hirst-typeZ coefficient estimation + Mann-Kendall test+ (~30; **)– (ns)+ (~20; *)
[109] (2019)NW Spain (1)1997–2016Hirst-typeMann-Kendall test + Theil-Sen estimator+ (~100; **)
[110] (2024)Slovakia (2)2002–2022Hirst-typeLinear regression + Mann-Kendall test + Theil-Sen estimator+ (ns)– (ns)
[111] (2002)Switzerland1951–2000Hirst-type + phenology– (~6)

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; ns: non-significant; NE: northeast; NW: northwest; N: north; S: south; LOESS: locally estimated scatterplot smoothing; SW: southwest; ARIMA: autoregressive integrated moving average; W: west; NL: Netherlands.