Whereas temporal variability of plant phenology in response to climate change has already been well studied, the spatial variability of phenology is not well understood. Given that phenological shifts may affect biotic interactions, there is a need to investigate how the variability in environmental factors relates to the spatial variability in herbaceous species’ phenology by at the same time considering their functional traits to predict their general and species-specific responses to future climate change. In this project, we analysed phenology records of 148 herbaceous species, which were observed for a single year by the PhenObs network in 15 botanical gardens. For each species, we characterised the spatial variability in six different phenological stages across gardens. We used boosted regression trees to link these variabilities in phenology to the variability in environmental parameters (temperature, latitude and local habitat conditions) as well as species traits (seed mass, vegetative height, specific leaf area and temporal niche) hypothesised to be related to phenology variability. We found that spatial variability in the phenology of herbaceous species was mainly driven by the variability in temperature but also photoperiod was an important driving factor for some phenological stages. In addition, we found that early-flowering and less competitive species characterised by small specific leaf area and vegetative height were more variable in their phenology. Our findings contribute to the field of phenology by showing that besides temperature, photoperiod and functional traits are important to be included when spatial variability of herbaceous species is investigated
