Katrin Ziegler, Felix Pollinger, Susanne Böll, Heiko Paeth

• Plant phenology is well known to be affected by meteorology. Observed changes in the occurrence of phenological phases arecommonly considered some of the most obvious effects of climate change. However, current climate models lack a representationof vegetation suitable for studying future changes in phenology itself. This study presents a statistical-dynamical modelingapproach for Bavaria in southern Germany, using over 13,000 paired samples of phenological and meteorological data foranalyses and climate change scenarios provided by aPlant phenology is well known to be affected by meteorology. Observed changes in the occurrence of phenological phases arecommonly considered some of the most obvious effects of climate change. However, current climate models lack a representationof vegetation suitable for studying future changes in phenology itself. This study presents a statistical-dynamical modelingapproach for Bavaria in southern Germany, using over 13,000 paired samples of phenological and meteorological data foranalyses and climate change scenarios provided by a state-of-the-art regional climate model (RCM). Anomalies of severalmeteorological variables were used as predictors and phenological anomalies of the flowering date of the test plantForsythiasuspensaas predictand. Several cross-validated prediction models using various numbers and differently constructed predictorswere developed, compared, and evaluated via bootstrapping. As our approach needs a small set of meteorological observationsper phenological station, it allows for reliable parameter estimation and an easy transfer to other regions. The most robust andsuccessful model comprises predictors based on mean temperature, precipitation, wind velocity, and snow depth. Its averagecoefficient of determination and root mean square error (RMSE) per station are 60% and ± 8.6 days, respectively. However, theprediction error strongly differs among stations. When transferred to other indicator plants, this method achieves a comparablelevel of predictive accuracy. Its application to two climate change scenarios reveals distinct changes for various plants andregions. The flowering date is simulated to occur between 5 and 25 days earlier at the end of the twenty-first century comparedto the phenology of the reference period (1961–1990).…