Next Article in Journal
Design of an Integrated Micro-Viscometer for Monitoring Engine Oil
Next Article in Special Issue
Towards an Automated Approach for Monitoring Tree Phenology Using Vehicle Dashcams in Urban Environments
Previous Article in Journal
Kalman-Filter-Based Learning of Characteristic Profiles of Lithium-Ion Batteries
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterizing Seasonal Radial Growth Dynamics of Balsam Fir in a Cold Environment Using Continuous Dendrometric Data: A Case Study in a 12-Year Soil Warming Experiment

1
Centre d’Étude de la Forêt, Université du Québec à Montréal, Montreal, QC H3C 3P8, Canada
2
Direction de la Recherche Forestière, Ministère des Forêts, de la Faune et des Parcs du Québec, Quebec City, QC G1P 3W8, Canada
3
Science and Technology Branch, Environment Canada and Climate Change, Montreal, QC H2Y 2E7, Canada
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(14), 5155; https://doi.org/10.3390/s22145155
Submission received: 30 May 2022 / Revised: 5 July 2022 / Accepted: 7 July 2022 / Published: 9 July 2022
(This article belongs to the Special Issue Sensors and Their Application in Phenological Studies)

Abstract

:
Historical temperature records reveal that the boreal forest has been subjected to a significant lengthening of the thermal growing season since the middle of the last century, and climate models predict that this lengthening will continue in the future. Nevertheless, the potential phenological response of trees to changes in growing season length remains relatively undocumented, particularly for evergreen boreal tree species growing in cold environments. Here, we used the recently defined zero growth (ZG) concept to extract and characterize the metrics of seasonal radial growth dynamics for 12 balsam fir trees subjected to a 12-year soil warming experiment using high resolution radius dendrometer measurements. The ZG concept provides an accurate determination of growth seasonality (onset, cessation, duration, growth rates, and total growth) for these slow-growing trees characterized by significant shrinkage in tree diameter due to dehydration in the winter. Our analysis revealed that, on average, growth onset starts at day 152 ± 7 (±1 SE, 31 May–1 June) and ceases at day 244 ± 27 (31 August–1 September), for a growing season duration of about 3 months (93 ± 26 days) over a 12-year period. Growing season duration is mainly determined by growth cessation, while growth onset varies little between years. A large part (80%) of the total growth occurs in the first 50 days of the growing season. Given the dynamics of growth, early growth cessation (shorter growing season) results in a higher average seasonal growth rate, meaning that longer growing seasons are not necessarily associated with greater tree growth. Soil warming induces earlier growth cessation, but increases the mean tree growth rate by 18.1% and the total annual growth by 9.1%, on average, as compared to the control trees. Our results suggest that a higher soil temperature for warmed trees contributes to providing better growth conditions and higher growth rates in the early growing season, when the soil temperature is low and the soil water content is elevated because of snowmelt. Attaining a critical soil temperature earlier, coupled with lower soil water content, may have contributed to the earlier growth cessation and shorter growing season of warmed trees.

1. Introduction

The examination of historical daily temperature records reveals that the Canadian forest has been subjected to significant climate warming since the middle of the last century, with even greater warming in the boreal forest, typically characterized by colder temperatures [1]. In addition to the general warming trend, an average lengthening, by about 15 days, of the thermal growing season (days with daily mean temperature above 5 °C), starting 6 days earlier and ending 9 days later, coupled with increased cumulative growing degree days, has been reported for the 1948–2016 period [1]. These trends in thermal growing season length are predicted to continue throughout this century, potentially adding more than a month of favorable conditions for the yearly growth of trees and plants [2,3]. Together with the fertilization effect of rising atmospheric CO2 and N inputs on tree growth, the lengthening of the growing season could potentially contribute to a large increase in C sinks observed over the last decades in northern temperate and boreal forests [4,5].
Nevertheless, the potential phenological response of trees to changes in climate seasonality remains relatively undocumented, particularly for the cold environment of the boreal forest. Several approaches have been used to document the seasonal growth dynamics of trees, each with their strengths and weaknesses. Phenological metrics linked to vegetation activity onset, cessation, duration, and peak can be derived over a large spatial scale from satellite observations of the photosynthetic activity of vegetation [6,7]. These metrics are, however, generally restricted to the analyses of pixels with a strong detectable annual cycle of vegetation, filtering out evergreen boreal coniferous forest regions [8,9]. Characterizing vegetation dynamics from satellite information in cold environments characterized by a short growing season, short daily sunshine duration in winter, persistence of snowpack in early spring, and the dominance of coniferous tree species remains challenging [10]. Moreover, phenological metrics derived from satellite images do not inform about species-specific variations, nor seasonal tree growth dynamics and tree productivity. Indeed, the greening effect documented from the satellite-derived normalized difference vegetation index (NDVI) does not always translate into ground-observed phenological trends or tree growth rates [7,11,12]. The extraction of metrics of tree seasonality using remote sensing techniques also involves technical challenges for taking into account the quality of the sensors and other distortions in the satellite images, and large discrepancies exist between ground-observed and satellite phenological transitional dates [8]. Phenological metrics derived from satellite information are also restricted in terms of temporal and spatial resolution, which are, for example, 16 days and 250 m resolution for MODIS NDVI images [13].
Transitional phenological dates can also be assessed from ground-based observations, including networks of leaf and bloom measurements of representative cloned cultivars (e.g., [14]), networks of digital phenological cameras (e.g., [15]), and local monitoring initiatives [16]. However, the assessment of phenological metrics, such as leafing and blooming, are also generally limited to deciduous trees, and they do not inform on tree growth dynamics over the entire growing season and on resultant tree productivity, particularly for the evergreen boreal tree species.
At the tree level, cambial activity and xylem formation (xylogenesis) of boreal species can be characterized using cellular analysis of repeatedly sampled micro-cores of newly formed xylem throughout the growing season [17] and by the pinning method [18]. These techniques are, however, laborious and time consuming in addition to damaging to the trees, limiting the long-term monitoring of individual trees. Characterization of growth dynamics using these techniques is also limited in terms of temporal resolution by the sampling (or pinning) interval. In addition, shifting sample points along the stem during the growing season generates additional variability in the measurements.
Over the past decade, a growing number of studies have used electronic dendrometers to track variations in tree radius or circumference (radius hereafter) at high temporal resolution in order to document the critical dates of xylem phenology [19,20,21,22,23,24,25]. These metrics of tree growth phenology include the timing of growth onset, cessation, and maximum growth rate, along with growing season duration and cumulative seasonal growth.
However, there are also limitations to determining phenological transition dates and quantifying tree growth using point dendrometers. In dendrometer data, both the signal of irreversible tree growth and the reversible signal of stem size, due to water storage depletion and replenishment of wood tissues, are confounded. These variations correspond to three distinct previously documented phases (growth increment, contraction or shrinkage, and expansion or recovery) of the circadian cycle of stem radius variation [26,27,28] and may be extracted with the help of a function for stem cycles definition included in the Dendrometer R package [29]. In addition to the diurnal cycle that characterizes stem size variation, periods of stem shrinkage may occur over a longer period, from a few days during the growing season, due to dry conditions of air and soil [30,31], to seasonal (i.e., months) during the winter, due to the dehydration of wood tissues [22,32]. During these periods of stem shrinkage, the signal is still characterized by cycles of stem contraction and expansion, but at a value below a precedent maximum stem size. Zweifel et al. [30] thus defined the zero growth (ZG) concept, according to which no growth occurs during periods of stem shrinkage. Using this concept, the signal of stem size variation may be partitioned into a growth (GRO) phase, which begins when the previous maximum stem radius is exceeded and ends when stem shrinkage occurs, and a tree water deficit (TWD) phase, which begins when the stem size decreases below a precedent maximum. The GRO phase derived from the ZG concept is slightly different from the increment phase of the circadian cycle approach [28]. The beginning of the GRO phase is determined from a previous maximum stem size, while the increment phase begins when the stem size exceeds the morning maximum.
One advantage of the ZG approach is that the GRO signal can be cumulated over the growing season to estimate yearly cumulative growth, which is not possible with the three-cycle approach. In the three-cycle approach, the same variation in stem size may be considered to be an increment phase when it occurs several days in a row; this can thus greatly overestimate total growth, if the stem radius variation during the increment phases is cumulated over the whole growing season [24,28]. Nevertheless, it should be noted that both approaches can be reconciled with a slight change in definition. Indeed, the TWD signal derived from the ZG concept could be decomposed into phases of stem contraction, shrinkage or dehydration, and stem expansion, recovery, or rehydration, thus restricting the increment phase to the GRO phase.
Another advantage of the ZG concept is that the timing of seasonal growth onset and cessation can be respectively derived from the maximum stem radius recorded over the previous season and during the current growing season [22]. From these metrics, growing season duration, cumulative growth, and average growth rate can also be derived. This simplification is particularly relevant for slow-growing trees inhabiting cold environments characterized by a long period of winter dehydration for which spring growth onset may be confused with stem rehydration before the beginning of cambial growth [32,33]. In this context, the three-cycle approach does not seem appropriate for precisely identify the timing of growth onset and cessation [24], and deriving the timing of growth onset during the spring rehydration phase is practically impossible without resorting to the ZG concept, or without the use of a complementary technique, such as micro-coring [27,32,33]. Seasonal growth dynamics have already been estimated from raw dendrometer data by fitting a statistical model to the data [19,20,34,35,36]. As opposed to the ZG concept, this approach results in a continuous growth curve (CG approach) characterizing the seasonal growth dynamics. Even with this approach, it remains difficult to precisely determine growth onset and cessation, and the beginning of the growing season remains practically undistinguishable from the early spring rehydration period.
The aims of this study are to test the use of the ZG concept for extracting and characterizing the metrics of seasonal radial balsam fir growth dynamics using high temporal resolution dendrometer measurements and to assess interannual variability, as well as the effects of a 12-year soil warming experiment on radial growth. In addition to potentially promoting better growth throughout the growing season, the early warming of soil temperatures in spring can induce favorable conditions for xylogenesis when air temperatures are above critical temperatures [37], but delayed soil warming due to snow cover still prevents tree growth [38,39]. We thus hypothesized an earlier growth onset, a longer duration of the growing season, and a higher annual radial growth for warmed-soil trees as compared to controls.

2. Methods

2.1. Study Site

The experimental area is located at the boundary of the Lac Laflamme watershed (47°19′31.36″ N; 71°07′43.15″ W, 740–760 m.a.s.l.), which is one of the three study watersheds that are part of the Québec Forest Ecosystem Research and Monitoring Network (Réseau d’Étude et de Surveillance des Écosystèmes Forestiers; RESEF). The forest is dominated by mature balsam firs (Abies balsamea (L.) Mill.), with a small component of white spruce (Picea glauca (Moench) Voss) and paper birch (Betula papyrifera Marsh.). The site is on Precambrian charnockitic gneiss covered by sandy till, with a terrain slope of 15%. Climate norms (1981–2010) were −0.18 °C and 1460 mm for mean annual air temperature and mean annual precipitation, respectively. The forest is located in the Laurentian Mountains, Québec, Canada, and underwent defoliation by the hemlock looper (Lambdina fiscellaria, Guenée) from 2012–2014 [40]. The 2012 defoliation was concurrent with a severe drought in the area in July 2012. With only 23 mm of rain in July as compared to the 143 mm long-term average, it was the driest summer month ever observed since 1960 [41].

2.2. Soil Warming Experiment

The soil warming experiment was set up in autumn 2008 and the short-term (2–3 years) effects on soil ion fluxes and mineralizable carbon, as well as on xylogenesis and tree growth of mature balsam firs, have been documented [34,42,43,44]. The experiment consists of 12 balsam firs (6 heated and 6 control) around which 70 m of heating or non-heating cables were buried in a spiral pattern (0.9 to 2.5 m radius, approximately 30 cm between successive turns) at the interface of the organic and mineral soil (5–10 cm belowground). The soil temperature below each tree was measured continuously with a thermistor (model 107-L, Campbell Scientific, Logan, UT, USA) positioned between the cables and linked to a data logger (model CR1000, Campbell Scientific, Logan, UT, USA). The continuous measurements were averaged and recorded at 15 min intervals. The warming cables of the 6 heated trees were activated and deactivated with a data logger when the average difference of temperature between the heated and non-heated soils reached 3 °C and 5 °C, respectively, with the objective of maintaining an average difference of 4 °C. Since 2009, the soils were heated from the beginning of the spring snow melt until the first frost days in the fall.
Each of the twelve trees in the experiment was equipped with radius dendrometers (model DR1, Ecomatik, Munich, Germany), south oriented, and positioned at approximately 2 m height to ensure being above the snow cover. Stem radial variation was measured continuously, and data were averaged and recorded every 15 min using the same data logger used for the soil measurements.
The soil warming treatment was initially nested within two levels of artificial precipitation (no added NH4NO3 and three times the natural NH4NO3 concentration in ambient precipitation) applied to the canopy every week from mid-June to mid-September from 2009 to 2011 using a system of reservoirs, pumps, valves, pipes, and sprinklers. No short-term effects of this N-enriched precipitation on tree growth were detected [34], so we did not consider this treatment in our analysis.

2.3. Characterization of Tree Growth Seasonality from Dendrometer Data

The extraction of parameters defining the seasonal dynamics of tree radial growth was performed with the help of the treenetproc package in R, V 0.1.4 [22,45,46]. These parameters included growth onset and cessation (day of year; DOY), growing season length (days), total annual growth (µm), and the average daily growth rate (µm/day). Spring growth onset is defined as the DOY at which the stem radius surpasses the maximum value recorded in the past year, and growth cessation is considered as the DOY at which the maximum value is reached (Figure 1, [22]). The length of the growing season is the difference between these two values, and total annual growth is the difference between the radius measurements corresponding to the beginning and end of growth. Finally, the average daily growth rate is the total annual growth divided by the length of the growing season.
Prior to the seasonal tree growth parameter extraction, the raw 15 min dendrometer data were visually examined and validated for any erroneous data that could be manually removed. As input, the treenetproc R package requires raw dendrometer data and air temperature. The purpose of the latter is to identify frost shrinkage periods in order to avoid assigning dendrometer data as outliers [22]. Internal data logger temperature data were used as an indicator of air temperature. This ensures accurate synchronicity between the ambient temperature and dendrometer measurements. The treenetproc package was then used to temporally align the dataset, filling the dataset via linear interpolation for small gaps (fewer than 3 time steps), and correcting for erroneous shifts and jumps in the manually cleaned dendrometer data [22]. We adjusted the parameters controlling the rigidity of outlier (tol_out) and jump (tol_jump) detection so that only jumps in the data associated with dendrometer adjustment or replacement were applied.
The treenetproc R package, used to determine growth onset and cessation (grow_seas), includes a tolerance factor parameter to adjust the determination of growth onset and cessation (tol_seas). This factor allows us to take into account the asymptotic shape of the seasonal growth pattern at the beginning and end of the growing season, as well as the presence of erratic measurements during these periods that can influence the determination of the beginning and end of growth [22]. The default value is set to 5%, so that the start and end of the growing season corresponds to the DOY when 5% and 95% of the total annual growth is completed, respectively. We set the parameters to the default value for the main analysis and also set the timing to correspond to the time when 80% and 90% of the total annual growth was completed (tol_seas value of 10% and 20%). We varied the tol_seas parameter only for the determination of growth cessation because this parameter has little influence on the timing of growth initiation, which is generally characterized by a much steeper slope of stem radius variation compared to the growth cessation phase.
The 12 dendrometers were operational from March 2009 to December 2021. Since the calculation of growth onset depends on the maximum radius measurements of the previous year, seasonal parameters were characterized for each tree, starting in the second year of measurements (2010) until 2021.

2.4. Statistical Analyses

Interannual variability in parameters characterizing tree growth seasonality (growth onset and cessation, growing season length, total annual growth, and the average daily growth rate), along with the soil warming treatment effect, were assessed by means comparisons and from mixed between-subjects (control/heated) and within-subjects (years) analysis of variance of the trimmed mean [47]. The trimmed mean excludes the lowest and the largest values of the distribution and computes the arithmetic mean on the remaining values. Analysis of variance based on trimmed means has been proposed as a robust alternative to classical statistical methods, when distributional assumptions are not met [47]. The bwtrim function of the WRS2 R package V 1.1-1 was used for this purpose [47]. We performed one analysis considering the 12 growing seasons (2010–2021) of the 7 trees that survived the 2012 event and another by including the 3 growing seasons (2010–2012) of the 5 trees that did not survive. Additionally, the strength of the monotonic association among the parameters characterizing the seasonal dynamics of radial growth of individual trees (i.e., experimental units) were analyzed using Spearman’s rank-order correlation.

3. Results

3.1. Soil Temperature and Dendrometer Data

The snow cover melted earlier, spring soil warming occurred, on average, 25 days earlier (heated and control soils reach 2 °C on day 113 and 138, respectively, Figure 2), and the mean annual May to October temperature of the heated soils was, on average, 3.8 ± 0.2 °C (±SE) higher compared to that of the unheated soils over the period 2009–2021.
Based on the dendrometer data, the seasonal parameters could not be evaluated for some trees in some years. The defoliation and severe drought that occurred in the area in July 2012 had a marked impact on the growth of trees in the subsequent years, and 5 trees died in 2013 from these events (trees 102, 111, 115, 124, and 125). The stem of another tree failed during an episode of sustained high winds that occurred from 12 to 15 August 2021 (tree 122). The parameters used to define the growing season could not be evaluated for three other seasons due to malfunctioning dendrometers. For two of these seasons (tree 111 in 2011 and tree 124 in 2012), the dendrometers had to be replaced due to water infiltration and ice formation in the sensor that prevented linear movement of the potentiometer. This phenomenon was observed in one case at the end of March and in the other during the month of November, two periods characterized by alternating rain and snow. Following this observation, plastic bibs were stapled to the trunk of the tree above the sensor to protect it from rain and snow, which prevented the problem from recurring. Finally, the 2018 season could not be evaluated for one tree (122) because the potentiometer had reached its maximum stroke during the growing season, and the late readjustment of the sensor made it impossible to join the series of measurements before and after the period with missing measurements. In the end, the studied database includes 4.83 million soil temperature measurements and 4.85 million stem radius measurements at a temporal resolution of 15 min, from which we characterized the tree growth seasonality of 94 of the 144 seasons (Figures S1–S12).

3.2. Tree Growth Seasonality

Annual, treatment, and overall averages of parameters characterizing the seasonal dynamics of balsam fir radial growth are presented in Table 1. Over the 12 studied years, the growth onset of balsam fir trees occurred, on average, on day 152 ± 7 (31 May–1 June; Table 1 and Table S1). The earliest seasonal onset was observed in 2013, 10 days earlier than the average (day 142, 22 May), while the latest season, observed in 2019, started 11 days later than the average (day 163, 12 June). The average growth onset of control (DOY 151) and heated (DOY 152) trees did not differ significantly, and the soil warming treatment did not influence yearly variability (P ≥ 0.1 for treatment–year interaction; Table 2).
Growth cessation (95% of yearly growth completed) occurred, on average, on day 244 ± 27 (31 August–1 September). The earliest end date was observed in 2021, 33 days sooner than the average (day 211, 30 July), while the latest end date was observed in 2018, when the season ceased 26 days later than the average (day 270, 27 September). The growth cessation of warmed trees (day 240) occurred, on average, 1 week (7 days) earlier than for the control trees (day 248, P < 0.1; Table 1, Table 2 and Table S2), and the effect of the soil warming treatment varied among years (P < 0.1 for treatment–year interaction; Table 2), with maximum differences of 24 and 60 days for the years 2017 and 2018, respectively.
On average, the resultant growing season length was evaluated as 93 ± 26 days, with longer growing seasons for control trees (96 days) than for heated trees (88 days, P < 0.001; Table 1, Table 2 and Table S3) and a less significant interaction of the soil warming treatment with the years (P < 0.1, Table 2). On average, the total annual growth ranged from 0.66 mm in 2014 to 1.58 mm in 2017 and averaged 1.16 mm over the 12-year period (Table 1, Table 2 and Table S4). The total annual growth of the heated trees is 9.1% higher than for the control trees (P < 0.01; Table 1 and Table 2), and this difference does not vary significantly over the years (P > 0.1 for treatment–year interaction, Table 2). Trees grew, on average, at a mean daily growth rate of 13.7 µm per day, with an 18.1% higher growth rate for heated (15.0 µm d−1) compared to control (12.7 µm d−1) trees (P < 0.001, Table 1, Table 2 and Table S5). There is no significant difference in the effect of the warming treatment on tree growth rates within years (P > 0.1 for the treatment–year interaction, Table 2). As for total annual growth, the lowest overall average growth rates were observed in 2013 (8.0 µm d−1) and 2014 (5.9 µm d−1).
While, on average, 95% of the total annual growth is completed by day 244 (Table 1), 80% of total growth is completed 43 days earlier (DOY 201, Table S6), and 90% of the total growth is completed 25 days earlier (DOY 222, Table S7). Considering the lower threshold value for characterizing the end of the growing season (80% or 90% of the total annual growth instead of 95%) considerably reduces the between-trees and within-years variability of this seasonal parameter and the corresponding growing season length (Table 1 and Tables S6–S9). The effect of the warming treatment on the timing of growth cessation and growing season is more tenuous when considering a lower threshold for the characterization of growth cessation (Table 2). Differences in growth cessation, growing season length (95% of yearly growth completed), and daily growth rates between heated and control trees remain statistically, but marginally, significant (P < 0.1) when including the 3 growing seasons of the 5 trees that did not survive the combined defoliation and drought event of 2012 (Table S10). Nevertheless, no difference in the total annual growth between treatments was detected.

3.3. Correlations among the Seasonal Parameters

The analysis of the correlation among parameters characterizing radial growth seasonality revealed that the total annual growth is closely and positively related to the average daily growth rate (p = 0.88, P < 0.001), the latter being negatively related to the timing of growth cessation (95% of yearly growth completed, p = −0.59, P < 0.001) and to growth duration (p = −0.60, P < 0.001), but not to the timing of growth onset (p = −0.04, Table 3). The duration of the growing season is closely and positively related to the timing of growth cessation if a 95% threshold was considered (p = 0.93, P < 0.001), and the relationship remain statistically significant for the two other thresholds analyzed (80% and 90%). Finally, the timing of growth cessation is positively related to the timing of growth onset, with a tighter relationship when the 80% threshold was considered (P < 0.001, Table 3).

4. Discussion

4.1. Tree Growth Seasonality

This study is the first to document seasonal radial growth dynamics of trees in cold environments over a relatively long period from high resolution radius dendrometer measurements using the ZG concept. Our analysis over the 12-year period revealed that the onset of growth of balsam fir trees occurs, on average, on day 152 ± 7 (±1 SE, 31 May–1 June). The timing of growth onset and cessation of naturally growing balsam fir trees located at a nearby site (approximately 500 m away) has already been estimated from raw dendrometer data using the CG approach [20]. From this analysis, an earlier average growth onset (by 12 days) was identified for the 2004–2010 period. For the 2010 and 2011 growing season, D’Orangeville et al. [34], also using the CG approach, identified that growth onset occurred 25 and 28 days sooner, respectively, for the same trees studied here.
The earlier onset of growth derived from the CG approach may be induced by the confusion between spring growth onset and stem rehydration occurring before the beginning of cambial growth (Figure 1). The ZG concept by which growth onset was determined from the previous maximum stem radius recorded over the previous season apparently makes it possible to accurately identify the beginning of cambial growth, even during the phase of tree stem enlargement due to the rehydration of trees in early spring (Figure 1). A close look at Turcotte et al.’s study [32] revealed that the onset of wood formation in black spruce trees, as determined by cellular analysis of repetitively sampled micro-cores, coincides with the time when the radius of the tree exceeds the maximum value recorded by dendrometers over the previous season. This growth onset was May 27 in 2005, which is very close to our average growth onset estimate over the 12-year period (5 days earlier, DOY 147 vs. 152). Turcotte et al.’s study [32] is among the few documenting the large magnitude of tree shrinkage occurring during the winter period for trees growing in the cold area of the Canadian boreal forest. In most cases, trees are monitored during the frost-free period, but according to the ZG concept, trees should be monitored continuously for growth onset determination, including during the late fall and winter period of stem shrinkage due to stem dehydration, and during the early-spring stem expansion period due to stem rehydration [22].
Our dating of seasonal growth cessation shows much greater interindividual and interannual variability than growth onset estimates, as well as greater divergence from growth cessation dated using the CG approach (Table 1; [20,34]). This is consistent with earlier observations that this metric shows greater variability as compared to growth onset, which was not as well explained by the climatic variables over the 2004–2010 period for trees growing nearby [20]. The asymptotic curve of the late growing season growth makes the determination of growth cessation more complex when compared to the early growing season period, which is characterized by more abrupt changes in the dendrometer records (Figure 1 and Figures S1–S12; [20]). The CG approach is very sensitive to stem shrinkage due to dry conditions occurring at the end of the growing season, or to the onset of winter desiccation in late fall, as the asymptote of the modeled progressive growth curve will be underestimated by attempting to minimize the sum of the squared residuals at the end of the growing season [20]. Dating of growth cessation using the ZG approach is also a challenge. From this concept, identification of growth cessation relies on the maximum radius size recorded in a given year. Because this maximum value depends on a single high-resolution observation, very small erratic fluctuations in the signal occurring at the end of the growing season can lead to great variability in determining the end-of-season dates, given the asymptotic shape of the dendrometer data during this period. To overcome this problem, Knüsel et al. [22] proposed the use of a tolerance factor for characterizing growth onset and cessation by excluding the first and last observations corresponding to a certain proportion of the total annual growth. A sensitivity analysis of this factor for the application of the ZG concept revealed a much greater influence on dating of growth cessation than growth onset, for the reason mentioned above. The use of the lower threshold value reduces the variability in the timing of growth cessation by avoiding the slow growth phase of the late growing season, while capturing a large part of the annual growth (Figures S1–S12). The time it takes to reach a certain threshold (e.g., 50% of yearly growth completed) and to move from one threshold to another could also be studied [36].
Correlation analyses revealed that the length of the growing season was closely related to the timing of growth cessation (p = 0.93, P < 0.001, Table 3), but not to the timing of growth onset. This indicates that the late and slow-growing part of the season is the most determinant of the seasonal growth duration. However, growth cessation and duration are only marginally significantly related to total annual growth (p = −0.21 to −0.22, P < 0.1), indicating that growth phenology metrics (onset, cessation, duration) are not the only drivers of tree growth and productivity. Seasonal growth was also a result of the growth rate during the growing season, with the latter being strongly positively related to cumulative seasonal growth (p = 0.88, P < 0.001, Table 3). The average growth rate over the season is also negatively related to the timing of growth cessation (p = −0.59, P < 0.001). Given the dynamics of cell growth, a season that ends abruptly will result in a higher growth rate over a short period of time, whereas seasons that last longer will have a protracted period of lower average growth and a lower overall growth rate. The weak relationship observed with the start of the growing season underscores the importance of late-season weather conditions in determining growing season duration. Using the CG approach, Duchesne et al. [20] similarly observed that growing season duration for balsam fir trees growing in the vicinity of the present study from 2004 to 2010 was primarily dictated by growth cessation, while growth onset was much less variable between years. The synchronicity of growth onset over the years suggests that day length or photoperiod, a parameter that is much more constant from year to year compared to spring climatic conditions, may be an important factor determining the onset of the growing season [48,49]. In the context of climate change, this raises concerns about the adaptation of tree species to the future synchrony of climate seasonality and photoperiod.

4.2. Soil Warming Effect

We hypothesized earlier growth onset, longer duration of the growing season, and superior annual radial growth for warmed-soil trees compared to controls. Contrary to our hypothesis, our results show that growth onset of control and heated trees did not differ significantly, relativizing the potential effects of earlier snowpack melting and higher soil temperature on boreal tree seasonal growth onset. The latter is apparently more closely linked to photosynthetically active radiation and air temperature [20]. Moreover, contrary to our hypothesis, the growth of warmed trees ceased, on average, 1 week earlier and lasted 8 days less than that of the control trees. Finally, a part of our hypothesis was verified, given that (despite a shorter growing season) warmed-soil trees had an 18.1% higher average growth rate and 9.1% higher total annual growth. This suggests a positive effect of soil heating early in the growing season, but a negative effect later in the growing season. Among the factors potentially favoring a higher growth rate, it has been shown that soil warming increased the cumulative nutrient fluxes (measured using ion-exchange membranes) in the forest floor and the mineral horizon during the first 3 years of this experiment [43]. A higher soil temperature for warmed trees (Figure 2) can also contribute to providing better growth conditions and higher growth rates, particularly in the early growing season, when soil water content is elevated because of snowmelt [50]. The effect of soil warming in the late growing season may, however, exacerbate the general declining trend in soil water content throughout the growing season until September, when higher precipitation and lower plant water uptake contributed to the soil water replenishment at values similar to those noted in the spring [50]. In 2006, an anomalous early growth cessation (DOY 187, 6 July) and a low seasonal radial growth (550 µm) for trees growing nearby were related to an early peak in soil temperature [20,51]. No other climatic variable was found to be correlated with the timing of growth cessation over the 2004–2010 period. Therefore, exceeding a critical soil temperature (Figure 2) early in the season, coupled with lower soil water content, may have contributed to the earlier growth cessation and shorter growing season of warmed trees. Another likely explanation for the early growth cessation of fast-growing trees is the control of photosynthesis by sink activity [52]. According to this concept, higher seasonal productivity leads to earlier growth cessation due to carbon sink limitations. However, this concept is apparently not supported by consistent observations from the experiments [53].
As reported after a 3-year-long manipulation with soil warming at the studied site [42], a twin experiment at two other sites showed no significant effect of soil warming on the phenological phases of cell enlargement, wall thickening, and lignification, nor on the number of cells produced by black spruce trees (Picea mariana, Mill. BSP) after 3 [54], and 6-year-long manipulations [55]. These studies characterized the xylogenesis of trees from cellular analysis of weekly sample micro-cores. The temporal resolution of the method used makes it difficult to detect a difference of a few days in tree growth phenology in response to the soil warming treatment, as we document in the present study. Here, the relatively small treatment effect could be detected due to the high temporal resolution of the approach used to document tree growth phenology. Treatment effects may also be greater over a 12-year period compared to analyses conducted over a 3- or 6-year period. Consistent with these studies, our analysis nevertheless reveals a relatively small effect of the soil warming treatment compared to the interannual variability in tree growth seasonality induced by other climatic parameters, notably air temperature [42,54,55].

4.3. Effects of Drought and Insect Defoliation

A total of 5 of the 12 instrumented trees died as a result of insect defoliation and the concurrent drought that occurred in 2012. The results indicate that 4 of the 5 dead trees showed the lowest radial growth among the 12 trees in the years preceding the event (Table S4), suggesting that the least vigorous trees are primarily those that succumbed. The impact on the growth of the surviving trees is noticeable 2 to 4 years after the event, when growth was significantly reduced compared to other years (Table S4). The death of these trees is unfortunate for the monitoring of the ongoing experiment, but these trees continue to be monitored with dendrometers. This monitoring represents an opportunity to document the etiology of tree death. Diurnal cycles are observable in the data, even several years after tree death. Analysis of these cycles, coupled with information on sap flow in living trees, could also unravel the influence of ambient air moisture and sap flow on stem radial variations.

5. Conclusions

This study is the first to document seasonal radial growth dynamics of trees in cold environments over a relatively long period of 12 years using the ZG concept. The latter seems well suited to determine the metrics of growth seasonality from dendrometer data for slow-growing trees subjected to significant winter shrinkage in tree diameter due to dehydration. The results contribute to our knowledge on tree growth seasonality in the boreal forest by showing that growth occurs over a short period of about 3 months, with most of the growth being achieved in the first 50 days of the growing season. The timing of growth cessation is the stronger determinant of growing season duration, while growth onset has a smaller influence. The average growth rate of trees is higher during growing seasons that end early and is lower as seasons lengthen, meaning that longer seasons are not necessarily the most productive. Soil warming induced earlier growth cessation and shortened the growing season, but increased the mean growth rate and total annual growth of the trees. These results provide crucial information on the potential phenological response of trees to changes in thermal growing season length to better anticipate the effects of climate change on boreal forest productivity. However, as is often the case in this type of study, the results are based on the monitoring of one tree species at a single site. Additional studies are needed to improve our knowledge of the phenological response of boreal forest tree species to climate change.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/s22145155/s1, Table S1: Growth onset (DOY; 5% of yearly growth completed) for control and soil-heated balsam fir trees over a 12-yr period (2010–2021); Table S2: Growth cessation (DOY; 95% of yearly growth completed) for control and soil-heated balsam fir trees over a 12-yr period (2010–2021); Table S3: Growing season duration (days; from 5% to 95% of yearly growth completed) for control and soil-heated balsam fir trees over a 12-yr period (2010–2021); Table S4: Total annual growth (µm) for control and soil-heated balsam fir trees over a 12-yr period (2010–2021); Table S5: Average daily growth rate (µm/day) for control and soil-heated balsam fir trees over a 12-yr period (2010–2021); Table S6: Growth cessation (DOY; 80% of yearly growth completed) for control and soil-heated balsam fir trees over a 12-yr period (2010–2021); Table S7: Growth cessation (DOY; 90% of yearly growth completed) for control and soil-heated balsam fir trees over a 12-yr period (2010–2021); Table S8: Growing season duration (days; from 5% to 80% of yearly growth completed) for control and soil-heated balsam fir trees over a 12-yr period (2010–2021); Table S9: Growing season duration (days; from 5% to 90% of yearly growth completed) for control and soil-heated balsam fir trees over a 12-yr period (2010–2021); Table S10: Mixed between subjects (control/heated) and within subjects (years) analyze of variance of the trimmed means of parameters used to characterize growth phenology for 12 balsam fir trees over the 2010–2021 period, including the 5 trees that not survived from the combined drought and defoliation event occurring in 2012; Figures S1–S12: Characterization of growing season from stem radius variations of balsam fir trees recorded by a dendrometer at a temporal resolution of 15 min over a twelve-year period.

Author Contributions

Conceptualization: L.D. and D.H.; data curation: S.O. and L.D.; formal analysis, software, and visualisation: S.O. and L.D.; investigation: L.D.; methodology and original draft preparation: L.D. and S.O.; project administration and supervision: D.H., L.D. and D.K.; validation, writing, review, and editing: L.D., S.O., D.H. and D.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded jointly by the ministère des Forêts, de la Faune et des Parcs (Quebec, Canada, project number 142332064 conducted at the Direction de la recherche forestière and led by Louis Duchesne) and a MITACS Cluster grant (funding reference number IT22943, granted to Daniel Houle and Daniel Kneeshaw).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vincent, L.; Zhang, X.; Mekis, E.; Wan, H.; Bush, E. Changes in Canada’s Climate: Trends in Indices Based on Daily Temperature and Precipitation Data. Atmos. Ocean. 2018, 56, 332–349. [Google Scholar] [CrossRef] [Green Version]
  2. Christidis, N.; Stott, P.; Brown, S.; Karoly, D.; Caesar, J. Human Contribution to the Lengthening of the Growing Season during 1950–99. J. Clim. 2007, 20, 5441–5454. [Google Scholar] [CrossRef]
  3. Wang, J.; Guan, Y.; Wu, L.; Guan, X.; Cai, W.; Huang, J.; Dong, W.; Zhang, B. Changing Lengths of the Four Seasons by Global Warming. Geophys. Res. Lett. 2021, 48, e2020GL091753. [Google Scholar] [CrossRef]
  4. Pan, Y.; Birdsey, R.A.; Fang, J.; Houghton, R.; Kauppi, P.E.; Kurz, W.A.; Phillips, O.L.; Shvidenko, A.; Lewis, S.L.; Canadell, J.G.; et al. A Large and Persistent Carbon Sink in the World’s Forests. Science 2011, 333, 988–993. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Friedlingstein, P.; O’Sullivan, M.; Jones, M.W.; Andrew, R.M.; Hauck, J.; Olsen, A.; Peters, G.P.; Peters, W.; Pongratz, J.; Sitch, S.; et al. Global Carbon Budget 2020. Earth Syst. Sci. Data 2020, 12, 3269–3340. [Google Scholar] [CrossRef]
  6. Myneni, R.; Keeling, C.D.; Tucker, C.J.; Asrar, G.; Nemani, R. Increased plant growth in the northern high latitudes from 1981 to 1991. Nature 1997, 386, 698–702. [Google Scholar] [CrossRef]
  7. Correa-Díaz, A.; Romero-Sánchez, M.; Villanueva-Díaz, J. The greening effect characterized by the Normalized Difference Vegetation Index was not coupled with phenological trends and tree growth rates in eight protected mountains of central Mexico. For. Ecol. Manag. 2021, 496, 119402. [Google Scholar] [CrossRef]
  8. Zhu, W.; Tian, H.; Xu, X.; Pan, Y.; Chen, G.; Lin, W. Extension of the growing season due to delayed autumn over mid and high latitudes in North America during 1982–2006. Glob. Ecol. Biogeogr. 2011, 21, 260–271. [Google Scholar] [CrossRef]
  9. Park, H.; Jeong, S.-J.; Ho, C.-H.; Park, C.-E.; Kim, J. Slowdown of spring green-up advancements in boreal forests. Remote Sens. Environ. 2018, 217, 191–202. [Google Scholar] [CrossRef]
  10. Beck, P.S.; Atzberger, C.; Høgda, K.A.; Johansen, B.; Skidmore, A.K. Improved monitoring of vegetation dynamics at very high latitudes: A new method using MODIS NDVI. Remote Sens. Environ. 2006, 100, 321–334. [Google Scholar] [CrossRef]
  11. Berner, L.; A Beck, P.S.; Bunn, A.; Lloyd, A.H.; Goetz, S. High-latitude tree growth and satellite vegetation indices: Correlations and trends in Russia and Canada (1982–2008). J. Geophys. Res. Earth Surf. 2011, 116, G01015. [Google Scholar] [CrossRef] [Green Version]
  12. Richardson, A.D.; Hufkens, K.; Milliman, T.; Frolking, S. Intercomparison of phenological transition dates derived from the PhenoCam Dataset V1.0 and MODIS satellite remote sensing. Sci. Rep. 2018, 8, 5679. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Didan, K. MOD13Q1 MODIS/Terra Vegetation Indices 16-Day L3 Global 250m SIN Grid V006 [Data Set]. NASA EOSDIS Land Processes DAAC; 2015. Available online: https://lpdaac.usgs.gov/products/mod13q1v006/ (accessed on 1 March 2022).
  14. Schwartz, M.D.; Ahas, R.; Aasa, A. Onset of spring starting earlier across the Northern Hemisphere. Glob. Chang. Biol. 2006, 12, 343–351. [Google Scholar] [CrossRef]
  15. Richardson, A.D.; Hufkens, K.; Milliman, T.; Aubrecht, D.M.; Furze, M.E.; Seyednasrollah, B.; Krassovski, M.B.; Latimer, J.M.; Nettles, W.R.; Heiderman, R.R.; et al. Ecosystem warming extends vegetation activity but heightens vulnerability to cold temperatures. Nature 2018, 560, 368–371. [Google Scholar] [CrossRef] [PubMed]
  16. Keenan, T.; Gray, J.; Friedl, M.A.; Toomey, M.; Bohrer, G.; Hollinger, D.Y.; Munger, J.W.; O’Keefe, J.; Schmid, H.P.; Wing, I.S.; et al. Net carbon uptake has increased through warming-induced changes in temperate forest phenology. Nat. Clim. Chang. 2014, 4, 598–604. [Google Scholar] [CrossRef]
  17. Rossi, S.; Anfodillo, T.; Menardi, R. Trephor: A New Tool for Sampling Microcores from tree stems. IAWA J. 2006, 27, 89–97. [Google Scholar] [CrossRef] [Green Version]
  18. Seo, J.-W.; Eckstein, D.; Schmitt, U. The pinning method: From pinning to data preparation. Dendrochronologia 2007, 25, 79–86. [Google Scholar] [CrossRef]
  19. Miller, T.W.; Stangler, D.F.; Larysch, E.; Honer, H.; Seifert, T.; Kahle, H.-P. A methodological framework to optimize models predicting critical dates of xylem phenology based on dendrometer data. Dendrochronologia 2022, 72, 125940. [Google Scholar] [CrossRef]
  20. Duchesne, L.; Houle, D.; D’Orangeville, L. Influence of climate on seasonal patterns of stem increment of balsam fir in a boreal forest of Québec, Canada. Agric. For. Meteorol. 2012, 162–163, 108–114. [Google Scholar] [CrossRef]
  21. Van Der Maaten, E.; Pape, J.; Van Der Maaten-Theunissen, M.; Scharnweber, T.; Smiljanić, M.; Cruz-García, R.; Wilmking, M. Distinct growth phenology but similar daily stem dynamics in three co-occurring broadleaved tree species. Tree Physiol. 2018, 38, 1820–1828. [Google Scholar] [CrossRef]
  22. Knüsel, S.; Peters, R.; Haeni, M.; Wilhelm, M.; Zweifel, R. Processing and Extraction of Seasonal Tree Physiological Parameters from Stem Radius Time Series. Forests 2021, 12, 765. [Google Scholar] [CrossRef]
  23. Rossi, S.; Deslauriers, A.; Anfodillo, T.; Morin, H.; Saracino, A.; Motta, R.; Borghetti, M. Conifers in cold environments synchronize maximum growth rate of tree-ring formation with day length. New Phytol. 2006, 170, 301–310. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Deslauriers, A.; Rossi, S.; Anfodillo, T. Dendrometer and intra-annual tree growth: What kind of information can be inferred? Dendrochronologia 2007, 25, 113–124. [Google Scholar] [CrossRef] [Green Version]
  25. Cruz-García, R.; Balzano, A.; Čufar, K.; Scharnweber, T.; Smiljanić, M.; Wilmking, M. Combining Dendrometer Series and Xylogenesis Imagery—DevX, a Simple Visualization Tool to Explore Plant Secondary Growth Phenology. Front. For. Glob. Chang. 2019, 2, 1–13. [Google Scholar] [CrossRef]
  26. Downes, G.; Beadle, C.; Worledge, D. Daily stem growth patterns in irrigated Eucalyptus globulus and E. nitens in relation to climate. Trees 1999, 14, 102–111. [Google Scholar] [CrossRef]
  27. Deslauriers, A.; Morin, H.; Urbinati, C.; Carrer, M. Daily weather response of balsam fir (Abies balsamea (L.) Mill.) stem radius increment from dendrometer analysis in the boreal forests of Québec (Canada). Trees 2003, 17, 477–484. [Google Scholar] [CrossRef]
  28. Deslauriers, A.; Rossi, S.; Turcotte, A.; Morin, H.; Krause, C. A three-step procedure in SAS to analyze the time series from automatic dendrometers. Dendrochronologia 2011, 29, 151–161. [Google Scholar] [CrossRef] [Green Version]
  29. van der Maaten, E.; van der Maaten-Theunissen, M.; Smiljanić, M.; Rossi, S.; Simard, S.; Wilmking, M.; Deslauriers, A.; Fonti, P.; von Arx, G.; Bouriaud, O. dendrometeR: Analyzing the pulse of trees in R. Dendrochronologia 2016, 40, 12–16. [Google Scholar] [CrossRef]
  30. Zweifel, R.; Haeni, M.; Buchmann, N.; Eugster, W. Are trees able to grow in periods of stem shrinkage? New Phytol. 2016, 211, 839–849. [Google Scholar] [CrossRef] [Green Version]
  31. Zweifel, R.; Sterck, F.; Braun, S.; Buchmann, N.; Eugster, W.; Gessler, A.; Häni, M.; Peters, R.L.; Walthert, L.; Wilhelm, M.; et al. Why trees grow at night. New Phytol. 2021, 231, 2174–2185. [Google Scholar] [CrossRef]
  32. Turcotte, A.; Morin, H.; Krause, C.; Deslauriers, A.; Thibeault-Martel, M. The timing of spring rehydration and its relation with the onset of wood formation in black spruce. Agric. For. Meteorol. 2009, 149, 1403–1409. [Google Scholar] [CrossRef]
  33. Tardif, J.; Flannigan, M.; Bergeron, Y. An Analysis of the Daily Radial Activity of 7 Boreal Tree Species, Northwestern Quebec. Environ. Monit. Assess. 2001, 67, 141–160. [Google Scholar] [CrossRef] [PubMed]
  34. D’Orangeville, L.; Houle, D.; Côté, B.; Duchesne, L.; Morin, H. Increased soil temperature and atmospheric N deposition have no effect on the N status and growth of a mature balsam fir forest. Biogeosciences 2013, 10, 4627–4639. [Google Scholar] [CrossRef] [Green Version]
  35. van der Maaten, E.; Bouriaud, O.; van der Maaten-Theunissen, M.; Mayer, H.; Spiecker, H. Meteorological forcing of day-to-day stem radius variations of beech is highly synchronic on opposing aspects of a valley. Agric. For. Meteorol. 2013, 181, 85–93. [Google Scholar] [CrossRef]
  36. D’Orangeville, L.; Itter, M.; Kneeshaw, D.; Munger, J.W.; Richardson, A.D.; Dyer, J.M.; A Orwig, D.; Pan, Y.; Pederson, N. Peak radial growth of diffuse-porous species occurs during periods of lower water availability than for ring-porous and coniferous trees. Tree Physiol. 2021, 42, 304–316. [Google Scholar] [CrossRef] [PubMed]
  37. Rossi, S.; DesLauriers, A.; Griçar, J.; Seo, J.-W.; Rathgeber, C.B.K.; Anfodillo, T.; Morin, H.; Levanic, T.; Oven, P.; Jalkanen, R. Critical temperatures for xylogenesis in conifers of cold climates. Glob. Ecol. Biogeogr. 2008, 17, 696–707. [Google Scholar] [CrossRef]
  38. Alvarez-Uria, P.; Körner, C. Low temperature limits of root growth in deciduous and evergreen temperate tree species. Funct. Ecol. 2007, 21, 211–218. [Google Scholar] [CrossRef]
  39. Wilson, G.; Green, M.; Brown, J.; Campbell, J.; Groffman, P.; Durán, J.; Morse, J. Snowpack affects soil microclimate throughout the year. Clim. Chang. 2020, 163, 705–722. [Google Scholar] [CrossRef]
  40. Delisle, J.; Bernier-Cardou, M.; Labrecque, A. Extreme cold weather causes the collapse of a population of Lambdina fiscellaria (Lepidoptera: Geometridae) in the Laurentian Mountains of Québec, Canada. Can. Èntomol. 2019, 151, 311–328. [Google Scholar] [CrossRef]
  41. Houle, D.; Lajoie, G.; Duchesne, L. Major losses of nutrients following a severe drought in a boreal forest. Nat. Plants 2016, 2, 16187. [Google Scholar] [CrossRef]
  42. D’Orangeville, L.; Côté, B.; Houle, D.; Morin, H.; Duchesne, L. A three-year increase in soil temperature and atmospheric N deposition has minor effects on the xylogenesis of mature balsam fir. Trees 2013, 27, 1525–1536. [Google Scholar] [CrossRef]
  43. D’Orangeville, L.; Houle, D.; Côté, B.; Duchesne, L. Soil response to a 3-year increase in temperature and nitrogen deposition measured in a mature boreal forest using ion-exchange membranes. Environ. Monit. Assess. 2014, 186, 8191–8202. [Google Scholar] [CrossRef] [PubMed]
  44. D’Orangeville, L.; Côté, B.; Houle, D.; Whalen, J. Reduced mineralizable carbon in a boreal forest soil after three years of artificial warming. Can. J. Soil Sci. 2013, 93, 567–572. [Google Scholar] [CrossRef]
  45. Wickham, H.; François, R.; Henry, L.; Müller, K. Dplyr: A Grammar of Data Manipulation, R Package Version 0.8.3. 2020. Available online: https://dplyr.tidyverse.org/ (accessed on 1 March 2022).
  46. Haeni, M.; Knüsel, S.; Wilhelm, M.; Peters, R.L.; Zweifel, R. Treenetproc—Clean, Process and Visualise Dendrometer Data, R Package Version 0.1.4. 2020. Available online: https://Github.Com/Treenet/Treenetproc (accessed on 1 March 2022).
  47. Mair, P.; Wilcox, R. Robust statistical methods in R using the WRS2 package. Behav. Res. Methods 2020, 52, 464–488. [Google Scholar] [CrossRef] [PubMed]
  48. Adams, J. Durartion of Light and Growth. Ann. Bot. 1924, os-38, 509–523. [Google Scholar] [CrossRef]
  49. Kramer, P.J. Effect of Variation in Length of Day on Growth and Dormancy of Trees. Plant Physiol. 1936, 11, 127–137. [Google Scholar] [CrossRef] [Green Version]
  50. D’Orangeville, L.; Houle, D.; Duchesne, L.; Côté, B. Can the Canadian drought code predict low soil moisture anomalies in the mineral soil? An analysis of 15 years of soil moisture data from three forest ecosystems in Eastern Canada. Ecohydrology 2015, 9, 238–247. [Google Scholar] [CrossRef]
  51. Duchesne, L.; Houle, D. Modelling day-to-day stem diameter variation and annual growth of balsam fir (Abies balsamea (L.) Mill.) from daily climate. For. Ecol. Manag. 2011, 262, 863–872. [Google Scholar] [CrossRef]
  52. Zani, D.; Crowther, T.W.; Mo, L.; Renner, S.S.; Zohner, C.M. Increased growing-season productivity drives earlier autumn leaf senescence in temperate trees. Science 2020, 370, 1066–1071. [Google Scholar] [CrossRef]
  53. Norby, R.J. Comment on “Increased growing-season productivity drives earlier autumn leaf senescence in temperate trees”. Science 2021, 371, 1689–1699. [Google Scholar] [CrossRef]
  54. Lupi, C.; Morin, H.; DesLauriers, A.; Rossi, S. Xylogenesis in black spruce: Does soil temperature matter? Tree Physiol. 2012, 32, 74–82. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Dao, M.C.E.; Rossi, S.; Walsh, D.; Morin, H.; Houle, D. A 6-Year-Long Manipulation with Soil Warming and Canopy Nitrogen Additions does not Affect Xylem Phenology and Cell Production of Mature Black Spruce. Front. Plant Sci. 2015, 6, 877. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Example of characterization of a growing season using stem radius variations of a balsam fir tree recorded by a dendrometer at a temporal resolution of 15 min. The day of the year corresponding to the maximum value recorded the previous season (1), the day of the year at which stem radius surpasses the maximum value recorded in the past year (2), and the day of the year at which the maximum value is reached (3) are labelled.
Figure 1. Example of characterization of a growing season using stem radius variations of a balsam fir tree recorded by a dendrometer at a temporal resolution of 15 min. The day of the year corresponding to the maximum value recorded the previous season (1), the day of the year at which stem radius surpasses the maximum value recorded in the past year (2), and the day of the year at which the maximum value is reached (3) are labelled.
Sensors 22 05155 g001
Figure 2. Average annual pattern of daily mean soil temperature (2010–2021) of control and heated plots. The shaded area represents the standard error.
Figure 2. Average annual pattern of daily mean soil temperature (2010–2021) of control and heated plots. The shaded area represents the standard error.
Sensors 22 05155 g002
Table 1. Average metrics of the seasonal radial balsam fir growth dynamics for control and soil-heated trees over a 12-year period (2010–2021).
Table 1. Average metrics of the seasonal radial balsam fir growth dynamics for control and soil-heated trees over a 12-year period (2010–2021).
YearsMeans
201020112012201320142015201620172018201920202021
Growth onset (DOY; 5% of yearly growth completed)
Control147155147140149149155156157164157146151
Heated145160148146157148153153154162155145152
Overall146157147142152149154155156163156146152
Growth cessation (DOY; 95% of yearly growth completed)
Control249232236241273241264248294259242214248
Heated254232242242261226268223234242232206240
Overall252232239242268235266237270252238211244
Growing season duration (days; from 5% to 95% of yearly growth completed)
Control1027789102124921099213795856896
Heated10972959610477116718080766188
Overall106759299116861128311488816693
Total annual growth (µm)
Control157216721260679547685896129211071055110210651115
Heated10181031107982481813011513195714091504154710901216
Overall1295132211787416639491160157712281248129310731160
Average daily growth rate (µm/day)
Control15.421.616.67.44.38.08.515.17.911.313.915.812.7
Heated11.414.912.78.98.116.413.926.917.618.720.417.915.0
Overall13.417.914.88.05.911.610.820.211.814.416.716.513.7
Table 2. Mixed between-subjects (control/heated) and within-subjects (years) analysis of variance of the trimmed means of parameters used to characterize the growth phenology of 7 balsam fir trees over the 2010–2021 period, excluding the 5 trees that did not survive the combined drought and defoliation events occurring in 2012.
Table 2. Mixed between-subjects (control/heated) and within-subjects (years) analysis of variance of the trimmed means of parameters used to characterize the growth phenology of 7 balsam fir trees over the 2010–2021 period, excluding the 5 trees that did not survive the combined drought and defoliation events occurring in 2012.
EffectOnsetCessationDurationTotalRate
5%80%90%95%80%90%95%
Treatment0.22.05.0 *9.6 *2.85.6 *11.2 **8.4 **21.0 ***
Year85.1 ***38.1 ***14.4 ***68.2 ***3.5 *3.0 *19.0 ***3.8 *9.6 ***
Interaction0.50.81.027.2 ***1.61.02.9 *0.91.3
Significance: * < 0.1, ** < 0.01, *** < 0.001.
Table 3. Spearman’s rank correlation coefficient among the parameters used to characterize growth phenology of 7 balsam fir trees over the 2010–2021 period, excluding the 5 trees that did not survive the combined drought and defoliation event occurring in 2012.
Table 3. Spearman’s rank correlation coefficient among the parameters used to characterize growth phenology of 7 balsam fir trees over the 2010–2021 period, excluding the 5 trees that did not survive the combined drought and defoliation event occurring in 2012.
Effect CessationDurationTotalRate
80%90%95%80%90%95%
Onset5%0.51 ***0.26 *0.23 *−0.28 **−0.17−0.05−0.08−0.04
Cessation80% 0.49 ***0.41 ***0.59 ***0.27 **0.28 **0.02−0.10
90% 0.68 ***0.30 **0.87 ***0.60 ***−0.24 *−0.45 ***
95% 0.22 *0.55 ***0.93 ***−0.22 *−0.59 ***
Duration80% 0.46 ***0.35 ***0.11−0.06
90% 0.62 ***−0.22 *−0.44 ***
95% −0.21 *−0.60 ***
Total 0.88 ***
Significance: * < 0.1, ** < 0.01, *** < 0.001.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Oogathoo, S.; Duchesne, L.; Houle, D.; Kneeshaw, D. Characterizing Seasonal Radial Growth Dynamics of Balsam Fir in a Cold Environment Using Continuous Dendrometric Data: A Case Study in a 12-Year Soil Warming Experiment. Sensors 2022, 22, 5155. https://doi.org/10.3390/s22145155

AMA Style

Oogathoo S, Duchesne L, Houle D, Kneeshaw D. Characterizing Seasonal Radial Growth Dynamics of Balsam Fir in a Cold Environment Using Continuous Dendrometric Data: A Case Study in a 12-Year Soil Warming Experiment. Sensors. 2022; 22(14):5155. https://doi.org/10.3390/s22145155

Chicago/Turabian Style

Oogathoo, Shalini, Louis Duchesne, Daniel Houle, and Daniel Kneeshaw. 2022. "Characterizing Seasonal Radial Growth Dynamics of Balsam Fir in a Cold Environment Using Continuous Dendrometric Data: A Case Study in a 12-Year Soil Warming Experiment" Sensors 22, no. 14: 5155. https://doi.org/10.3390/s22145155

APA Style

Oogathoo, S., Duchesne, L., Houle, D., & Kneeshaw, D. (2022). Characterizing Seasonal Radial Growth Dynamics of Balsam Fir in a Cold Environment Using Continuous Dendrometric Data: A Case Study in a 12-Year Soil Warming Experiment. Sensors, 22(14), 5155. https://doi.org/10.3390/s22145155

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop