Background: In tundra ecosystems, the adjustment of phenological events, such as bud burst, to snowmelt timing is crucial to the climatic adaptation of plants. Natural small-scale variations in microclimate potentially enable plant populations to persist in a changing climate. Aims: To assess how plant phenology responds to natural differences in snowmelt timing. Methods: We observed the timing of eight vegetative and reproductive phenophases in seven dwarf-shrub species in relation to differences in snowmelt timing on a small spatial scale in an alpine environment in subarctic Finland. Results: Some species and phenophases showed accelerated development with later snowmelt, thus providing full or partial compensation for the shorter snow-free period. Full compensation resulted in synchronous occurrence of phenophases across the snowmelt gradient. In other species, there was no acceleration of development. The timing of phenophases varied between two consecutive years and two opposing mountain slope aspects. Conclusions: The results have shown three distinct patterns in the timing of phenophases in relation to snowmelt and suggest alternative strategies for adaptation to snowmelt timing. These strategies potentially apply to other species and tundra ecosystems and provide a framework, enabling one to compare and generalise phenological responses to snowmelt timing under different future climate scenarios.

Age-related effects on whole-tree hydraulics are one of the key challenges to better predicting the production and growth of old-growth forests. Previous models have described the optimal state of stomatal behaviour, and field studies have implied on age/size-induced trends in tree ecophysiology related to hydraulics. On these bases, we built a Bayesian hierarchical model to link sap flow density and drivers of transpiration directly. The model included parameters with physiological meanings and accounted for variations in leaf-sapwood area ratio and the time lag between sap flow and transpiration. The model well-simulated the daily pattern of sap flow density and the variation between tree age groups. The results of parameterization show that (1) the usually higher stomatal conductance in young than old trees during midsummer was mainly because the sap flow of young trees were more activated at low to medium light intensity, and (2) leaf-sapwood area ratio linearly decreased while time lag linearly increased with increasing tree height. Uncertainty partitioning and cross-validation, respectively, indicated a reliable and fairly robust parameter estimation. The model performance may be further improved by higher data quality and more process-based expressions of the internal dynamics of trees.

Waterlogging causes hypoxic or anoxic conditions in soils, which lead to decreases in root and stomatal hydraulic conductances. Although these effects have been observed in a variety of plant species, they have not been quantified continuously over a range of water table depths (WTD) or soil water contents (SWC). To provide a quantitative theoretical framework for tackling this issue, we hypothesized similar mathematical descriptions of waterlogging and drought effects on whole-tree hydraulics and constructed a hierarchical model by connecting optimal stomata and soil-to-leaf hydraulic conductance models. In the model, the soil-to-root conductance is non-monotonic with WTD to reflect both the limitations by water under low SWC and by hypoxic effects associated with inhibited oxygen diffusion under high SWC. The model was parameterized using priors from literature and data collected over four growing seasons from Scots pine (Pinus sylvestris L.) trees grown in a drained peatland in Finland. Two reference models (RMs) were compared with the new model, RM1 with no belowground hydraulics and RM2 with no waterlogging effects. The new model was more accurate than the reference models in predicting transpiration rate (fitted slope of measured against modeled transpiration rate = 0.991 vs 0.979 (RM1) and 0.984 (RM2), R2 = 0.801 vs 0.665 (RM1) and 0.776 (RM2)). Particularly, RM2’s overestimation of transpiration rate under shallow water table conditions (fitted slope = 0.908, R2 = 0.697) was considerably reduced by the new model (fitted slope = 0.956, R2 = 0.711). The limits and potential improvements of the model are discussed.