Dendrochronological Analyses and Whole-Tree Dissections Reveal Caliciopsis Canker ( Caliciopsis pinea ) Damage Associated with the Declining Growth and Climatic Stressors of Eastern White Pine ( Pinus strobus )

: Eastern white pine ( Pinus strobus ) is considered a signature species in eastern North America, particularly in New England. In recent years, however, white pine has experienced increased damage due to native pathogens that reduce the species’ growth, productivity, and economic value. One disease of concern is Caliciopsis canker, caused by the fungal pathogen Caliciopsis pinea , which is associated with excessive resin production, cankers, rough bark, bark ﬁssures / cracks, and reduced growth in white pine. Recent studies have documented the extent of Caliciopsis canker in New England and its association with soil and stocking conditions, yet few studies have focused on the biological impacts of the disease. This study used dendrochronology and whole-tree dissections to reconstruct Caliciopsis canker history in three New England white pine sites, quantify its impact on tree growth and vigor, identify pre-disposing factors, and assess potential silvicultural management options. Dendrochronology and whole-tree dissections provided a unique insight into canker damage throughout trees’ development. Canker damage was ﬁrst reported in New Hampshire in the mid-1990s, yet cankers were present as far back as 1967 and have steadily increased since the mid-1980s. Increased canker damage was signiﬁcantly associated with decreased live crown ratios and declining tree growth. Trees maintaining a 30% live crown ratio or greater generally experienced the least canker damage. Furthermore, peaks in canker occurrence were consistent across sites, indicating a regional synchronization of infection and damage. Canker damage was closely associated with climatic events such as droughts and a New England hurricane. The results suggest that Caliciopsis canker has been a ﬀ ecting white pine health over the last 40 years, and that the disease has become more prevalent in the past 20–30 years. Yet, our results suggest that if silvicultural prescriptions target low density thinnings that favor trees with higher live crown ratios ( > 30%) and low Caliciopsis symptom severity ratings, the risk of canker damage can be reduced in white pine stands.


Introduction
There is growing evidence that forest pests and pathogens are causing increased damage worldwide [1]. As globalization continues to increase economic trade, pests and pathogens can likewise heavy resin production, particularly between whorls, (C) thinning crowns with as little as 5% live crown, (D) roughened bark in the mid-to upper-crown, (E) bark fissures/cracks, which are associated with (F) necrotic tissue and internal cankers. The (F) internal canker with necrotic tissue is associated with the (E) large bark crack that has sealed over.
Reports of canker occurrence and associated white pine damage have been increasing, particularly in eastern North America [24,28,[33][34][35][36]. Caliciopsis-related damage was reported in 1997 in the northeastern U.S. [28], by 2003 in Québec, Canada [37], and in 2006 in the southeastern and midwestern U.S. [30,35]. Some of the most severe damage occurs in pole-sized white pine stands, with more than 70% of sites in New Hampshire being symptomatic [28], as well as pockets of tree mortality in the midwestern, southeastern, and northeastern U.S. [24,28,30,33,34]. To monitor and evaluate Caliciopsis canker severity, the U.S. Forest Service recently developed an assessment protocol that recommends dividing white pine tree stems into thirds and counting the number of Caliciopsis canker-associated resin streaks per third [24].
To quantify tree responses to forest pests and pathogens, studies have frequently used dendrochronology. However, the majority of this research has focused on insects (e.g., [38][39][40][41]). Furthermore, past research has either reconstructed pest/pathogen arrival and damage history (e.g., [39,41,42]), or has focused on changes in tree growth and volume (e.g., [38,40]). To date, dendrochronology and whole-tree dissections have not been used to reconstruct the history of a fungal pathogen and quantify tree response to infections on a whole-tree level.
Objectives: Given the ecological and economic importance of white pine throughout its range, there is strong interest in maintaining the health and productivity of the species. Despite some basic biological information on Caliciopsis canker (discussed above), very little is known about its impact on white pine or the factors that may predispose trees to this disease. Our general objective was to heavy resin production, particularly between whorls, (C) thinning crowns with as little as 5% live crown, (D) roughened bark in the mid-to upper-crown, (E) bark fissures/cracks, which are associated with (F) necrotic tissue and internal cankers. The (F) internal canker with necrotic tissue is associated with the (E) large bark crack that has sealed over.
Reports of canker occurrence and associated white pine damage have been increasing, particularly in eastern North America [24,28,[33][34][35][36]. Caliciopsis-related damage was reported in 1997 in the northeastern USA [28], by 2003 in Québec, Canada [37], and in 2006 in the southeastern and midwestern USA [30,35]. Some of the most severe damage occurs in pole-sized white pine stands, with more than 70% of sites in New Hampshire being symptomatic [28], as well as pockets of tree mortality in the midwestern, southeastern, and northeastern USA [24,28,30,33,34]. To monitor and evaluate Caliciopsis canker severity, the USA Forest Service recently developed an assessment protocol that recommends dividing white pine tree stems into thirds and counting the number of Caliciopsis canker-associated resin streaks per third [24].
To quantify tree responses to forest pests and pathogens, studies have frequently used dendrochronology. However, the majority of this research has focused on insects (e.g., [38][39][40][41]). Furthermore, past research has either reconstructed pest/pathogen arrival and damage history (e.g., [39,41,42]), or has focused on changes in tree growth and volume (e.g., [38,40]). To date, dendrochronology and whole-tree dissections have not been used to reconstruct the history of a fungal pathogen and quantify tree response to infections on a whole-tree level.
Objectives: Given the ecological and economic importance of white pine throughout its range, there is strong interest in maintaining the health and productivity of the species. Despite some basic biological information on Caliciopsis canker (discussed above), very little is known about its impact on white pine or the factors that may predispose trees to this disease. Our general objective was to fill these knowledge gaps, using dendrochronological and field data from three study sites in northern New England. Specific objectives included (1) reconstructing the sequence of Caliciopsis cankers over time; (2) quantifying the relationship between Caliciopsis cankers and white pine growth and vigor; (3) identifying pre-disposing factors that may increase the likelihood of Caliciopsis canker incidence and severity; and (4) evaluating the response of stands affected by Caliciopsis canker to thinnings. We hypothesize that Caliciopsis canker is reducing white pine growth and vigor, particularly in stressed and weakened trees. Stressors likely include soils, climatic variables such as drought, and overstocking. We further hypothesize that trees from thinned stands will have lower canker incidence and severity, and better growth recovery post-infection. Taken together, the information gained by addressing these objectives may guide silvicultural prescriptions aimed at mitigating damage from future infections.

Study Area
Study sites are located in the Northeastern Highlands Ecoregion in New England, which includes the Adirondacks in New York, the majority of Vermont and New Hampshire, and western Maine [43]. The region is characterized by cool temperatures, with many subregions being transitional between temperate and boreal zones. Soil types vary, but spodosols are common, often with sandy loam or loamy sand. Soils are typically nutrient-poor and range from well drained to poorly drained. Mean annual soil temperatures are frigid, ranging from 0 to 8 • C [43].

Site and Stand Selection
Three sites were selected in northern New England: N1 and N2 in central New Hampshire and M3 in southwestern Maine. The sites had a history of Caliciopsis canker damage and occurred on previous agricultural land. Since the 1930s, sites developed without disturbance, except for the harvesting of a few random trees. From 2000-2010, a portion of each site was thinned to approximately 250 trees per hectare. Thinnings occurred in December 1999/January 2000 at N1, 2010 at N2, and 2004 at M3. Our study selected a thinned and unthinned stand from each site, in close proximity to one another (<0.5 km apart).
The sample sites had annual mean temperatures ranging from 6.9 to 7.6 • C, minimum temperatures from 0.8 to 1.2 • C, and maximum temperatures from 13.0 to 13.9 • C (Table A1) [44]. The annual precipitation ranges from 111.2 to 117.9 cm (Table A1) [44]. In general, the sites occurred on moderately to poorly drained sandy soils, where white pine comprised at least 75% of the site basal area [24].
During site selection, the severity of Caliciopsis canker-associated symptoms on each white pine tree was assessed, resulting in two severity classes: (1) low severity, with zero to few external symptoms; and (2) high severity, with numerous symptoms (fully defined in [32]). The presence of C. pinea in trees was confirmed using both a culture-based approach and a molecular diagnostic assay ( [46], Broders and Wyka unpubl. data).

Field Sampling
At the center of each stand, plots similar to Forest Inventory and Analysis (FIA) Phase 2 plots were established; the plots consisted of four 7.32-m radius subplots in a tetrahedral shape [47].
Within the subplots, all trees ≥1.0 cm diameter at breast height (dbh, 1.37 m) were sampled for species, dbh, and crown class [47]. Additional data were recorded for white pines ≥12.7 cm dbh, including tree height, uncompacted live crown ratio (UNCR), and crown metrics such as width, transparency, density, and dieback. Tree-and stand-level metrics were calculated including basal area, volume, and trees per hectare.
Following initial data collection, six white pine trees per stand were destructively sampled (36 trees total) from May-August, 2015. In each stand, the three most affected trees in the high severity class and the three least affected trees in the low severity class were selected. The trees were felled and subsequently dissected into 0.3 m long bolts (see Figure 1B for example of bolts). Each bolt was inspected for Caliciopsis canker damage (as defined in Figure 1), from the base of the tree up to a 5 cm diameter. If cankers were present, a 2.54 cm thick cross-section was cut from the bottom of the bolt. For each canker, the data collected included: (a) canker width, measured from necrotic edge to necrotic edge, along the circumference of the annual ring; (b) canker height on the tree, measured as the specific height of the tree from soil to the exposed canker; and (c) the presence/absence of external symptoms associated with the canker, including resinosis, roughened bark, brown/red lesions, or bark fissures/cracks. From these field cross-sections, every third sample was returned to the Tree-Ring Laboratory at the University of Maine, totaling 318 cross-sections. This subset was selected to ensure a manageable sample size. An additional cross-section was taken from every tree at breast height; all 36 of these cross-sections were also returned to the University of Maine.

Cross-Section Preparation and Tree Ring Measurements
Cross-sections were dried for 30+ days. Following drying, the bottom of each cross-section was sanded with progressively finer sandpaper (100 through 1200 grit) using belt and random orbital sanders, until individual cells were easily visible under a stereoscope. Processed samples were crossdated by hand using the skeleton plot method, then verified by two additional researchers. On the 318 cross-sections with cankers, all cankers were numbered and dated. The year of canker formation was the earliest year with cell disruption and necrosis. The year of recovery was the year that ring width was no longer declining at the canker, according to visual estimation ( Figure 2). density, and dieback. Tree-and stand-level metrics were calculated including basal area, volume, and trees per hectare. Following initial data collection, six white pine trees per stand were destructively sampled (36 trees total) from May-August, 2015. In each stand, the three most affected trees in the high severity class and the three least affected trees in the low severity class were selected. The trees were felled and subsequently dissected into 0.3 m long bolts (see Figure 1B for example of bolts). Each bolt was inspected for Caliciopsis canker damage (as defined in Figure 1), from the base of the tree up to a 5 cm diameter. If cankers were present, a 2.54 cm thick cross-section was cut from the bottom of the bolt. For each canker, the data collected included: (a) canker width, measured from necrotic edge to necrotic edge, along the circumference of the annual ring; (b) canker height on the tree, measured as the specific height of the tree from soil to the exposed canker; and (c) the presence/absence of external symptoms associated with the canker, including resinosis, roughened bark, brown/red lesions, or bark fissures/cracks. From these field cross-sections, every third sample was returned to the Tree-Ring Laboratory at the University of Maine, totaling 318 cross-sections. This subset was selected to ensure a manageable sample size. An additional cross-section was taken from every tree at breast height; all 36 of these cross-sections were also returned to the University of Maine.

Cross-Section Preparation and Tree Ring Measurements
Cross-sections were dried for 30+ days. Following drying, the bottom of each cross-section was sanded with progressively finer sandpaper (100 through 1200 grit) using belt and random orbital sanders, until individual cells were easily visible under a stereoscope. Processed samples were crossdated by hand using the skeleton plot method, then verified by two additional researchers. On the 318 cross-sections with cankers, all cankers were numbered and dated. The year of canker formation was the earliest year with cell disruption and necrosis. The year of recovery was the year that ring width was no longer declining at the canker, according to visual estimation ( Figure 2). All samples were scanned on an Epson Expression 1640 XL flatbed scanner (Epson America, Long Beach, CA, USA) at 800-1200 dpi (dependent on cross-section size and ring width), in color, and saved as TIFF files. Images from the 318 cross-sections with cankers were imported into ImageJ All samples were scanned on an Epson Expression 1640 XL flatbed scanner (Epson America, Long Beach, CA, USA) at 800-1200 dpi (dependent on cross-section size and ring width), in color, and saved as TIFF files. Images from the 318 cross-sections with cankers were imported into ImageJ 1.x, an image analysis program [48]. In ImageJ, lines were drawn over each canker to measure canker width. Width was measured from necrotic edge to necrotic edge, along the circumference of the ring.
Images from the 36 breast height cross-sections were imported and measured using the image analysis software WinDENDRO v.2012b (Regent Instruments, Québec, CA, Canada) [49]. Annual ring widths were measured to the nearest 0.001 mm from pith to bark along four separate paths, at approximately 90 • angles from one another. The four paths were averaged for each cross-section, creating one annual ring width chronology per cross-section. Crossdating accuracy was verified using COFECHA [50], and discrepancies were corrected where necessary.

Analyses
From the 318 cross-sections, 835 cankers were originally identified. However, questionable damage that was not definitively associated with Caliciopsis canker damage was removed from the dataset prior to analysis (65 of 835 cankers removed). These removed cankers were assessed through visual observation of damage after cross-sections had been sanded and included damage that did not have necrosis or that was associated with other factors (e.g., insects, abiotic stress, etc.). The remainder (770 cankers) had Caliciopsis canker symptoms ( Figure 1) that were associated with C. pinea infections in other studies ( [9,32,51], Broders and Wyka unpub. Data). Furthermore, DNA from a subset of cankers was tested using Caliciopsis spp. primers, and despite difficulty in isolating DNA from older cankers, C. pinea was successfully isolated [46]. The dates assigned to each canker were used to reconstruct the history of Caliciopsis canker incidence in each stand. It was assumed that if cankers were randomly dispersed in time, then an individual driver or event would not be associated with increases in canker incidence (number of cankers) and severity (width of cankers). Conversely, if cankers were clustered in specific years, it would indicate a disturbance affecting the host or an event favoring pathogen dispersal potentially coinciding with canker occurrence, particularly if all sites experienced identical cluster years.
To quantify tree-level response to Caliciopsis canker damage, two dependent variables (the number of cankers, the width of cankers) and three independent variables (UNCR, annual basal area increment, annual decline score) were calculated. The dbh and height of the sample trees varied dramatically and likely influenced the number of cankers measured per cross-section, as well as the number of bolts that showed canker damage. Therefore, developing dependent variables that accounted for tree size was critical to accurately assess canker damage. The number of cankers reflects the number of infections. Alternatively, the cumulative width of cankers can provide a more reliable estimate of canker influence on individual tree health and vigor. As a result, both metrics were included as dependent variables. These metrics were summed across the sample cross-sections analyzed (approximately one-third of the cankered sections from each tree). To adjust metrics for tree size, each canker metric was divided by the volume (m 3 ) of each sample (volume calculations in [52]). These adjusted canker metrics will be referred to as the canker incidence index and canker severity index.
Crown metrics (UNCR) were measured in the field for every sample tree, according to FIA Phase 2 Protocol [47]. A clinometer was used to measure the tree's base, base of live crown (lowest live branch, not compacted), and tip (living top of tree). The whole tree length was calculated (top-base) and then the crown length (top-base of live crown). Finally, the crown length was divided by the whole tree length to calculate UNCR. Additionally, each canker's height on the tree was recorded in the field, and rounded to the nearest 0.3 m.
To test if there were differences in tree measurements among sites and treatments, analyses of variance (ANOVA) were run in R version 3.4.3 [53]. Thinning treatment (thinned vs. unthinned) and Caliciopsis symptom severity class (low vs. high) were fixed effects, while dbh, height, volume, and UNCR were each run as individual response variables. Site was part of the blocking design and was treated as a random factor.
To assess the relationship between canker damage and UNCR (a measureable external characteristic of tree health), linear mixed effects models were run for each dependent variable (number of cankers/m 3 and width of cankers/m 3 ) using the lme4 package in R [53,54]. Thinning treatment, Caliciopsis symptom severity class, UNCR, and their interactions were fixed effects, while site was treated as a random factor.
Annual basal area increments (BAI) were calculated using raw ring width measurements, for every cross-section collected at dbh, using the bai.out function from the dplR package in R [53,55]. From these BAI values, an annual decline score (D score ) was calculated for each tree to identify the years in which trees were declining in growth and the years in which they were recovering (Equation (1)). This equation is a modified version of that developed by Livingston et al. [56].
Equation (1): Decline score equation using a modified version of Livingston et al.'s [56] D score methodology, which calculates an annual decline score for every tree, using annual basal area increments.
X = year for which D score is being calculated, BAI = basal area increment for year in subscript, VAR x = the pooled variance of the means (Equation (2)) Equation (2): Pooled variance equation.
Within the D score equation, a short three-year period was specifically selected to identify periods of decline that develop quickly, such as those in response to C. pinea infection, which might not otherwise be detected with longer periods [56].
In order to identify potential risk factors that may increase the likelihood of Caliciopsis canker damage, the D score was plotted against major regional stress events, determined from historical records, and canker damage metrics. Stress events were graphed as three-year windows: year one was the event, while years two and three were the potential response period. Events include the 1957 drought (1957)(1958)(1959), the 1965 drought (1965)(1966)(1967), the 1978 drought (1978)(1979)(1980) [57][58][59]. We hypothesized that lower D score values would be associated with (A) stress events, (B) higher canker damage metrics [24,33,36,60], (C) both stressors and canker damage, or (D) none of the above.
Potential risk factors were then evaluated by comparing the D score with the year of canker formation and periods of climatic stress, focusing on when cankers occurred during the decline cycle. A decline is any period with a D score < 0 for three or more consecutive years. The tree is in recovery when the D score is ≥ 0 for two or more consecutive years. This period between the start of decline and recovery is the decline cycle. The lowest D score value in a decline cycle is considered the decline maximum [56]. If canker damage consistently occurs at the start of decline cycles, it is assumed that canker damage induces declining growth. Conversely, if canker damage occurs at or after the decline maximum, it is presumed that the tree was already declining before infection. If canker damage occurs independent of decline cycles, it is assumed that tree growth is not directly related to canker damage. Additionally, if canker damage corresponds with the maximum D score and occurs with or immediately following climatic stress, the stressor is assumed to have predisposed the tree to fungal infection and canker formation [61]. These evaluations will improve our understanding of how white pine interacts with C. pinea, which will shed light on the risk factors that may increase canker incidence and severity in the species.
The role of Caliciopsis symptom severity was considered by evaluating the relationship between D scores and canker damage in high severity versus low severity trees at each site. We hypothesized that high severity trees would have more years with negative D scores , and that there would be greater canker damage in both positive and negative D score years. These data will allow us to assess if external symptom assessments (Caliciopsis symptom severity class and UNCR) can accurately predict which trees experience greater canker damage, and if highly affected trees experience greater growth declines, all with the intention of developing a successful management strategy for Caliciopsis canker.

Results
Caliciopsis cankers first appeared in sampled sites decades earlier than initial reports indicated: 1968 at N1, 1989 at N2, and 1974 at M3 (Figure 3). Cankers occurred in low numbers until the mid-1980s (M3) to mid-1990s (N1 and N2), at which point canker incidence increased and maintained higher numbers at all three sites (  pine interacts with C. pinea, which will shed light on the risk factors that may increase canker incidence and severity in the species. The role of Caliciopsis symptom severity was considered by evaluating the relationship between Dscores and canker damage in high severity versus low severity trees at each site. We hypothesized that high severity trees would have more years with negative Dscores, and that there would be greater canker damage in both positive and negative Dscore years. These data will allow us to assess if external symptom assessments (Caliciopsis symptom severity class and UNCR) can accurately predict which trees experience greater canker damage, and if highly affected trees experience greater growth declines, all with the intention of developing a successful management strategy for Caliciopsis canker.

Results
Caliciopsis cankers first appeared in sampled sites decades earlier than initial reports indicated: 1968 at N1, 1989 at N2, and 1974 at M3 (Figure 3). Cankers occurred in low numbers until the mid-1980s (M3) to mid-1990s (N1 and N2), at which point canker incidence increased and maintained higher numbers at all three sites (  Canker damage occurred throughout the bole length of sample trees ( Figure 4A-B). The lowest confirmed canker occurred at 2.1 m, and the highest confirmed canker occurred at 21.3 m (on a 27.9 m tall tree). The majority of cankers (601 of 770, or 78%) occurred between 6-16 m on the tree ( Figure  4A). The largest cankers (>15 cm wide) occurred in the middle of the bole, between 5-12 m ( Figure  4B). The butt log (0-3.66 m) is typically the most commercially valuable log. These results Canker damage occurred throughout the bole length of sample trees ( Figure 4A,B). The lowest confirmed canker occurred at 2.1 m, and the highest confirmed canker occurred at 21.3 m (on a 27.9 m tall tree). The majority of cankers (601 of 770, or 78%) occurred between 6-16 m on the tree ( Figure 4A).
The largest cankers (>15 cm wide) occurred in the middle of the bole, between 5-12 m ( Figure 4B). The butt log (0-3.66 m) is typically the most commercially valuable log. These results demonstrate that only 2.3% (18 of 770) of cankers occurred in this log ( Figure 4A), and those that did were relatively small (mean width = 1.39 cm; Figure 4B). demonstrate that only 2.3% (18 of 770) of cankers occurred in this log ( Figure 4A), and those that did were relatively small (mean width = 1.39 cm; Figure 4B).  Table A2). However, when thinning treatment was accounted for, trees from thinned stands were significantly larger with greater UNCR values than trees from unthinned stands (Table A2A). There were no consistent trends in tree measurements between trees with high and low Caliciopsis canker symptom severity (Table A2B). Furthermore, there were no significant interactions between thinning treatment and Caliciopsis symptom severity class for each tree metric (Table A2A and A2B).
The total number and summed width of cankers varied widely between sample trees, in part due to individual tree volumes ranging from 0.19 to 1.95 m 3 . On each tree, the canker incidence index ranged from 0 to 182.6 cankers (mean = 50.8 cankers), while the canker severity index ranged from 0 to 357.8 cm (mean = 84.8 cm). It should be noted that one outlier tree was removed from the dataset during analysis, as the tree was declining/dead at time of sampling and was identified as an outlier from the rest of the sample population. The outlier had a canker incidence index of 140.4 cankers and a canker severity index of 662.2 cm.
Linear mixed effects models indicate that canker incidence and canker severity were significantly and negatively related to UNCR (p-values were < 0.001 and 0.002, respectively). As UNCR increased, the incidence and severity of cankers decreased ( Figure 5A-B). The median values were 27.7% for UNCR, 36.7 for canker incidence index, and 45.7 cm for canker severity index (data not shown). Of the 18 trees with UNCR values greater than 27.7%, 14 had canker damage metrics (incidence and severity) below the median ( Figure 5A-B), supporting the hypothesis that higher crown ratios are associated with reduced canker damage. No other fixed effects (thinning treatment, Caliciopsis symptom severity class, and all interactions) were statistically significant in the mixed effects models.  Table A2). However, when thinning treatment was accounted for, trees from thinned stands were significantly larger with greater UNCR values than trees from unthinned stands (Table A2A). There were no consistent trends in tree measurements between trees with high and low Caliciopsis canker symptom severity (Table A2B). Furthermore, there were no significant interactions between thinning treatment and Caliciopsis symptom severity class for each tree metric (Table A2A,B).
The total number and summed width of cankers varied widely between sample trees, in part due to individual tree volumes ranging from 0.19 to 1.95 m 3 . On each tree, the canker incidence index ranged from 0 to 182.6 cankers (mean = 50.8 cankers), while the canker severity index ranged from 0 to 357.8 cm (mean = 84.8 cm). It should be noted that one outlier tree was removed from the dataset during analysis, as the tree was declining/dead at time of sampling and was identified as an outlier from the rest of the sample population. The outlier had a canker incidence index of 140.4 cankers and a canker severity index of 662.2 cm.
Linear mixed effects models indicate that canker incidence and canker severity were significantly and negatively related to UNCR (p-values were < 0.001 and 0.002, respectively). As UNCR increased, the incidence and severity of cankers decreased ( Figure 5A,B). The median values were 27.7% for UNCR, 36.7 for canker incidence index, and 45.7 cm for canker severity index (data not shown). Of the 18 trees with UNCR values greater than 27.7%, 14 had canker damage metrics (incidence and severity) below the median ( Figure 5A,B), supporting the hypothesis that higher crown ratios are associated with reduced canker damage. No other fixed effects (thinning treatment, Caliciopsis symptom severity class, and all interactions) were statistically significant in the mixed effects models. Tree growth declines (Dscores < 0) were frequently associated with increases in the number of cankers that year ( Figure 6A-B), particularly since the mid-1970s. In general, the number of cankers increased when declines occurred, then decreased once the trees recovered ( Figure 6A). At N1 and N2, the number of positive and negative Dscore years was similar (23 each; Figure 6B), while at M3 there were 10-16 more positive Dscore years than negative ( Figure 6B). However, there were consistently greater numbers of cankers in negative Dscore years (Figure 6B), indicating a correlation between decline and canker occurrence. In all groups except low severity trees from M3, there were two to 21 times as many cankers in negative Dscore years compared with positive Dscore years ( Figure  6B). Furthermore, high severity trees always produced greater numbers of cankers in decline years than did low severity trees ( Figure 6A), supporting the hypothesis that external severity ratings can be used to assess which trees are experiencing greater damage. Additionally, these decline metrics were closely associated with trends in BAI values ( Figure A1), which is expected as Dscores were calculated from BAI values. Tree growth declines (D scores < 0) were frequently associated with increases in the number of cankers that year ( Figure 6A,B), particularly since the mid-1970s. In general, the number of cankers increased when declines occurred, then decreased once the trees recovered ( Figure 6A). At N1 and N2, the number of positive and negative D score years was similar (23 each; Figure 6B), while at M3 there were 10-16 more positive D score years than negative ( Figure 6B). However, there were consistently greater numbers of cankers in negative D score years ( Figure 6B), indicating a correlation between decline and canker occurrence. In all groups except low severity trees from M3, there were two to 21 times as many cankers in negative D score years compared with positive D score years ( Figure 6B). Furthermore, high severity trees always produced greater numbers of cankers in decline years than did low severity trees ( Figure 6A), supporting the hypothesis that external severity ratings can be used to assess which trees are experiencing greater damage. Additionally, these decline metrics were closely associated with trends in BAI values ( Figure A1), which is expected as D scores were calculated from BAI values.  (1957-1959, 1965-1967, 1978-1980, 1981-1983, 1995-1997, 1999-2001, 2007-2009  When evaluating Dscores and climatic events individually, there were no consistent trends: during each of the climatic events and the two-year response period, some stands experienced decline, some recovered, and others did not decline ( Figure 6A). However, the years with the greatest number of cankers always occurred during a climatic stress period ( Figure 6A), suggesting that climatic stress is associated with, or increases the likelihood of, canker occurrence.
Similar trends emerged when evaluating the size of cankers (canker width) and tree growth declines. In each group, the year with the highest total canker width occurred during a period of growth decline ( Figure 7A). Across all sites, high severity trees experienced higher total canker widths, particularly high severity trees at N1 ( Figure 7A). In negative Dscore years, all sites experienced nearly twice as much canker damage in high severity trees compared to low severity trees ( Figure  7B). In all groups, the year with the largest amount of canker damage occurred during a climatic disturbance period ( Figure 7A).  (1957-1959, 1965-1967, 1978-1980, 1981-1983, 1995-1997, 1999-2001, 2007-2009  When evaluating D scores and climatic events individually, there were no consistent trends: during each of the climatic events and the two-year response period, some stands experienced decline, some recovered, and others did not decline ( Figure 6A). However, the years with the greatest number of cankers always occurred during a climatic stress period ( Figure 6A), suggesting that climatic stress is associated with, or increases the likelihood of, canker occurrence.
Similar trends emerged when evaluating the size of cankers (canker width) and tree growth declines. In each group, the year with the highest total canker width occurred during a period of growth decline ( Figure 7A). Across all sites, high severity trees experienced higher total canker widths, particularly high severity trees at N1 ( Figure 7A). In negative D score years, all sites experienced nearly twice as much canker damage in high severity trees compared to low severity trees ( Figure 7B). In all groups, the year with the largest amount of canker damage occurred during a climatic disturbance period ( Figure 7A).  (1957-1959, 1965-1967, 1978-1980, 1981-1983, 1995-1997, 1999-2001, 2007-2009

Discussion
The analyses presented above document a detailed investigation into the impacts of one specific tree disease, Caliciopsis canker. Trees often respond to stress agents relatively quickly, typically within three years. Annual growth rings, therefore, capture this response and provide a detailed record of tree-and stand-level disturbance history that can be used to document the impact of specific stressors and diseases. Our study reconstructed the history of one fungal pathogen, quantified its impact on tree crown and growth metrics, and applied whole-stem analysis to document tree-level infection incidence and severity. As a result, this study successfully used dendrochronological techniques for quantifying and assessing tree diseases.
The results from tree-level metrics (dbh, volume, crown metrics) showed that trees from thinned stands were larger (dbh, height, volume) with larger crowns (UNCR values; Table A2A), likely due to larger trees being retained at the time of thinning and improved growing conditions after thinning. Landowner records indicated that all thinned and unthinned stands at the three sites were similar and unmanaged until the most recent thinning. Our data shows that tree ages were similar between N1 and N2 stands, but M3 thinned trees were approximately 10 years younger than unthinned trees. Adjustments were made to account for tree size, however, there were still no significant differences in canker incidence or canker severity indices between thinned and unthinned stands (p > 0.10 for  (1957-1959, 1965-1967, 1978-1980, 1981-1983, 1995-1997, 1999-2001, 2007-2009

Discussion
The analyses presented above document a detailed investigation into the impacts of one specific tree disease, Caliciopsis canker. Trees often respond to stress agents relatively quickly, typically within three years. Annual growth rings, therefore, capture this response and provide a detailed record of tree-and stand-level disturbance history that can be used to document the impact of specific stressors and diseases. Our study reconstructed the history of one fungal pathogen, quantified its impact on tree crown and growth metrics, and applied whole-stem analysis to document tree-level infection incidence and severity. As a result, this study successfully used dendrochronological techniques for quantifying and assessing tree diseases.
The results from tree-level metrics (dbh, volume, crown metrics) showed that trees from thinned stands were larger (dbh, height, volume) with larger crowns (UNCR values; Table A2A), likely due to larger trees being retained at the time of thinning and improved growing conditions after thinning. Landowner records indicated that all thinned and unthinned stands at the three sites were similar and unmanaged until the most recent thinning. Our data shows that tree ages were similar between N1 and N2 stands, but M3 thinned trees were approximately 10 years younger than unthinned trees. Adjustments were made to account for tree size, however, there were still no significant differences in canker incidence or canker severity indices between thinned and unthinned stands (p > 0.10 for both indices). As the dendrochronology results show (Figures 6 and 7), much of the canker formation occurred prior to or just after thinning, indicating treatments likely had little long-term impact on a tree's response to C. pinea infections.
Based on the results from canker incidence ( Figure 4A) and severity ( Figure 4B) along the length of the tree bole, it is evident that most cankers occurred between 6-16 m, and these also tended to be the largest cankers. These results indicate that stressed white pine at these shorter heights are more likely to suffer canker damage, which supports the recommendation to reduce white pine stand density when trees are taller than 6 m [62].
When comparing canker incidence and severity to crown metrics, both the number of cankers and the width of cankers reported similar negative relationships with UNCR values, indicating that external visual assessments of crown health can be indicative of internal canker damage. Trees with an UNCR value of 30% or greater had reduced canker damage ( Figure 5A,B). When managing for reduced risk of Caliciopsis canker, we therefore recommend selecting residual trees with a 30% UNCR or greater, and a low canker severity rating, which will likely result in a productive, healthy stand that experiences less damage from Caliciopsis canker. Similarly, D'Amato et al. [63] analyzed white spruce (Picea glauca [Moench] Voss) live crown ratios in relation to eastern spruce budworm (Choristoneura fumiferana Clemens) attack. They found that a crown ratio of 40% or greater was optimal for maintaining tree vigor and growth and increasing stand resistance to budworm attack [63], which lends support to our suggestion of maintaining at least 30% UNCR in white pine stands.
Based on reconstructions of canker history, sites often experienced elevated canker incidence in the same years (Figure 3), indicating regional synchrony in canker occurrence. This suggests that, in addition to site factors, an external force is likely influencing canker incidence and severity. Recent studies on other fungal pathogens have shown a marked increase in infection occurrence associated with changes in climate [64], which could also be driving the increase in Caliciopsis canker damage on white pine. In the present study, the year with the greatest number and width of cankers at each site was associated with a drought period ( Figure 6A, Figure 7A). Further research on the direct impact of climate change on both the pathogen and the host are needed, however, to fully test these relationships.
While these results are significant, these data do not indicate cause and effect. To address whether canker damage results in tree decline or if decline predisposes trees to Caliciopsis canker, we evaluated the relationship between D scores and canker damage metrics. Increases in both canker number and canker width typically occur one to three years after growth decline has already begun, frequently initiated by a climatic event ( Figure 6A, Figure 7A). This finding suggests that canker damage occurs at the time of, or just after, another stress is already acting on the tree. This can further be seen in the overlap between increased canker years and climatic stress events ( Figure 6A, Figure 7A). When stressors occur, increases in canker damage occur one to two years later, suggesting that C. pinea causes increased infection and damage in predisposed, stressed trees.
Collectively, these results illustrate distinct associations between tree growth declines and Caliciopsis canker damage, providing evidence that supports recent observations of declining white pine stands in eastern North America [28,33,36]. In turn, this demonstrates a strong negative relationship between tree health and Caliciopsis canker occurrence, particularly as damage increased. The data suggest that when trees are healthy and maintaining consistent growth, Caliciopsis canker is not a significant problem. This relationship between declining trees and disease has been well studied, and it is commonly accepted that stressed trees are more susceptible to disease agents [14,65]. In Maine and New Hampshire, one of the most common climate stressors is drought [15], which reduces tree resources and potentially increases white pine susceptibility to Caliciopsis canker. Furthermore, when drought stress is combined with other stressors (e.g., harvesting, windstorms, flooding), susceptibility may continue to increase. One option is to thin stands to reduce density, increase resource availability, increase tree vigor, and reduce susceptibility [15]. D'Amato et al. [66] found that thinning red pine stands resulted in increased drought resistance, particularly in younger stands. These results underscore the need to minimize stress in white pine stands in order to reduce the susceptibility to Caliciopsis canker.

Conclusions
Native pests and pathogens can cause significant changes in New England's white pine forests, particularly when they act in ways not previously documented. While C. pinea was first documented in New England in 1880, reports of Caliciopsis canker damage began increasing in the early 1990s. This study confirms Caliciopsis canker damage, and further supports previous reports that Caliciopsis cankers have been present and increasing in incidence and severity since the mid-1970s in Maine and New Hampshire. The number and size of cankers are also closely linked to decreasing tree vigor. Climatic events that stress white pine (e.g., droughts, flooding) appear to be associated with increases in canker damage. These stress events typically occur 1-3 years before increases in canker damage, supporting the hypothesis that C. pinea is frequently associated with already stressed trees.
To reduce Caliciopsis canker damage, our findings suggest that proactive silviculture practices can be effective, particularly when retaining trees with live crown ratios greater than 30%, a low canker severity rating, and reducing stand densities when trees exceed 6 m in height. Selecting such trees as crop trees will likely result in increased tree vigor and growth, decreased susceptibility to Caliciopsis canker, and in turn a decreased number and size of cankers. These actions can maintain the value of standing white pine [32], and in some instances can even increase the value by allowing currently infected trees to recover. One important consideration, however, is the timing and intensity of thinning. Care should be taken to minimize stress to the stand; harvesting should not be conducted if other stress events are co-occurring or have occurred recently. In sum, Caliciopsis canker is a potentially detrimental fungal disease both ecologically and economically; however, with monitoring and proactive management, damage can potentially be reduced and white pine can continue its role as a signature species in New England.  Acknowledgments: We graciously thank the individuals who helped fell trees, collect field data, and process samples: Bruce Allen and Kyle Lombard with New Hampshire Division of Forests and Lands; William Ostrofsky and Wayne Searles with Maine Forest Service; Colin Bridge-Koenigsberg, Garth Dixon, Jereme Frank, Savannah Haines, Joshua Savoy, and Matthew Todd with the University of Maine; and Olivia Box with Hamilton College. We also appreciate the guidance and statistical expertise of Bill Halteman, who assisted with data analysis and interpretation.