Delphinium barbeyi is a common herbaceous wildflower in montane meadows at 2,900 m near the Rocky Mountain Biological Laboratory, and its flowers are important nectar resources for bumblebees and hummingbirds. During the period 1977-1999 flowering was highly variable in both timing (date of first flower ranged from 5 July to 6 August, mean=17 July) and abundance (maximum open flowers per 2×2-m plot ranged from 11.3 to 197.9, mean=82). Time and abundance of flowering are highly correlated with the previous winter's snowpack, as measured by the amount of snow remaining on the ground on 15 May (range 0-185 cm, mean=67.1). We used structural equation modeling to investigate relationships among snowpack, first date of bare ground, first date of flowering, number of inflorescences produced, and peak number of flowers, all of which are significantly correlated with each other. Snowpack depth on 15 May is a significant predictor of the first date of bare ground (R =0.872), which in turn is a significant predictor of the first date of flowering (R=0.858); snowpack depth is also significantly correlated with number of inflorescences produced (R =0.713). Both the number of inflorescences and mean date of first flowering are significant predictors of flowers produced (but with no residual effect of snowpack). Part of the effect of snowpack on flowering may be mediated through an increased probability of frost damage in years with lower snowpack - the frequency of early-season "frost events" explained a significant proportion of the variance in the number of flowers per stem. There is significant reduction of flower production in La Niña episodes. The variation in number of flowers we have observed is likely to affect the pollination, mating system, and demography of the species. Through its effect on snowpack, frost events, and their interaction, climate change may influence all of these variables.

In this paper, we test the mid-domain hypothesis as an explanation for observed patterns of flowering diversity in two sub-alpine communities of insect-pollinated plants. Observed species richness patterns showed an early-season increase in richness, a mid-season peak, and a late-season decrease. We show that a "mid-domain" null model can qualitatively match this pattern of flowering species richness, with R(2) values typically greater than 60%. We find significant or marginally significant departures from expected patterns of diversity for only 3 out of 12 year-site combinations. On the other hand, we do find a consistent pattern of departure when comparing observed versus null-model predicted flowering diversity averaged across years. Our results therefore support the hypothesis that ecological factors shape patterns of flowering phenology, but that the strength or nature of these environmental forcings may differ between years or the two habitats we studied, or may depend on species-specific characteristics of these plant communities. We conclude that mid-domain null models provide an important baseline from which to test departure of expected patterns of flowering diversity across temporal domains. Geometric constraints should be included first in the list of factors that drive seasonal patterns of flowering diversity.