Comparing the Cell Dynamics of Tree-Ring Formation Observed in Microcores and as Predicted by the Vaganov-Shashkin Model

New insights into the intra-annual dynamics of tree-ring formation can improve our understanding of tree-growth response to environmental conditions at high-resolution time scales. Obtaining this information requires, however, a weekly monitoring of wood formation, sampling that is extremely time-intensive and scarcely feasible over vast areas. Estimating the timing of cambial and xylem differentiation by modeling thus represents an interesting alternative for obtaining this important information by other means. Temporal dynamics of cambial divisions can be extracted from the daily tree-ring growth rate computed by the Vaganov-Shashkin (VS) simulation model, assuming that cell production is tightly linked to tree-ring growth. Nonetheless, these predictions have yet to be compared with direct observations of wood development, i.e.,viamicrocoring, over a long time span. We tested the performance of the VS model by comparing the observed and predicted timing of wood formation in black spruce [Picea mariana(Mill.)]. We obtained microcores over 15 years at 5 sites along a latitudinal gradient in Quebec (Canada). The measured variables included cell size and the timing of cell production and differentiation. We calibrated the VS model using daily temperature and precipitation recorded by weather stations located on each site. The predicted and observed timing of cambial and enlarging cells were highly correlated (R-2= 0.8); nonetheless, we detected a systematic overestimation in the predicted timing of cambial cells, with predictions delayed by 1-20 days compared with observations. The growth rate of cell diameter was correlated with the predicted growth rate assigned to each cambial cell, confirming that cell diameter developmental dynamics have the potential to be inferred by the tree-ring growth curve of the VS model. Model performances decrease substantially in estimating the end of wood formation. The systematic errors suggest that the actual relationships implemented in the model are unable to explain the phenological events in autumn. The mismatch between the observed and predicted timing of wood formation in black spruce within our study area can be reduced by better adapting the VS model to wet sites, a context for which this model has been rarely used.