Abstract
This study on investigating change in the material properties of ash trees after infestation by emerald ash borer (EAB) (Agrilus planipennis) occurred at two locations in northeast Ohio in the summer of 2013. The trees at either site were divided into three groups based on % canopy lost from EAB (group I = 0 to 5%, group II = 6 to 25%, and group III = greater than 40%). A digital image correlation (DIC) system was used to evaluate and compare strain (tissue deformation) on ash branches that were (static) loaded to failure. Stress wave transmission times (Tm) of sound waves through stem wood and fine roots and root balls of the ash trees also were assessed. The DIC evaluations revealed that branches of ash trees that were in groups II and III exhibited significantly lower strain after static loading compared to that observed for trees in group I. Analysis of stress wave Tm revealed that group III trees had significantly higher Tm times compared to the other two groups. Fine root necrosis was significantly higher in group III trees and lowest in group I trees. Extracted root balls from group III trees had significantly higher percentage decay compared to that observed from trees in groups I and II. These data provide fundamental insight into the material properties of ash trees after infestation by EAB and can contribute to arboricultural guidelines for ash tree preservation and help develop safety protocols to address structural loss in trees after EAB infestation.
INTRODUCTION
Since the discovery of the emerald ash borer (EAB), Agrilus planipennis Fairmare, in North America in 2002 in Michigan and Ontario (Cappaert et al. 2005; Poland and McCullough 2006), EAB has expanded its distribution to affect a significant proportion of the ash tree populations in North America. As an urban tree resource, Kovacs et al. (2010) estimate that 37.5 million ash trees in 25 eastern states of the U.S.A. may be affected and have to be removed. Cappert et al. (2005) estimate ash tree removal associated with EAB infested trees may cost between $20 and $60 billion US dollars, while Hauer and Peterson (2016) report on municipalities in affected areas being impacted annually by $280.5 million US dollars for removal and for preservation efforts. Early detection is crucial to any program geared towards preserving ash trees, and cues that assist monitoring programs in urban settings (e.g., in Persad et al. 2015) may help guide tree preservation efforts. Management choices for treating ash trees after infestation are limited to relatively narrow windows through which they can be effective. If preservation is not implemented, then the gradual loss of ash tree canopy occurs, resulting ultimately in tree degradation and removal prior to the trees posing unacceptable risk of human injury or property damage (Hauer 2012; Van Natta and Hauer 2012). When working in ash trees that are in decline or killed by EAB, arborist safety and tree stability considerations, especially near targets in urban environments, have received considerable attention. Working in and around EAB infested trees may involve added risk; several studies (Persad et al. 2013, 2015) demonstrate that reduced structural integrity of EAB-infested ash trees occurs well before actual tree death. Arboricultural safety guidance programs will continue to benefit from a growing body of data that investigate further effects of EAB on the whole tree and with focus on stem, branch, and root collar integrity.
At present no data exist on strain mapping of EAB affected branches, stem material properties, root plate stability, and fine root density of ash trees that are infested with EAB. Recent use of digital image correlation (DIC) software suggests that strain can be mapped remotely (Sebera et al.2014; Sebera et al.2016; Dahle 2017). To resolve this in the context of ash canopy decline (as a visual cue) this study was conducted during the Biomechanics Week of the International Society of Arboriculture (ISA) in August 2013. The 2010 Tree Biomechanics Summit recommended research in understanding the mechanism and mode of tree failures as well as effectiveness of mitigation practices (Dahle et al. 2014).
The study was designed to assist the practitioner with data that could assist in decision making with respect to treat, remove, and/or to estimate the level of risk associated with canopy loss. The study was conducted on three groups of ash trees that had varying levels of canopy loss from infestation by EAB. The first aim was to add to the existing body of knowledge on EAB-infested ash tree branch loading using a DIC stereoscopic photogrammetry system to map strain on branches that were subjected to static loading to failure.
Secondly, a Fakoop acoustics meter was used to compare transmission times of sound waves through the stem wood (Wang et al. 2004a; Dahle et al. 2016) of ash trees in each of the three ash tree groups. Stress waves when generated through wood help to define the physical and mechanical properties of the wood material and may travel at speeds which are dependent on several variables, including direction of travel in the wood in relation to wood grain and fiber, moisture content, and decay incidence (Armstrong 1991; Pellerin 2002; Pellerin 1985).
The third objective focused on the root plate zone and root balls of ash trees in each of the three ash tree groups; root areas were examined for fine root necrosis (root necromass/biomass ratios or Nm/Bm). One technique in assessing fine root production involves assessing the measurements of live (biomass, Bm) and dead (necromass, Nm) amounts of fine roots in terms of dry weights in sequential soil cores (Persson and Stadenburg 2009). Additionally, decay pockets in the root balls were identified and compared for each of three groups of EAB-infested ash trees.
MATERIALS AND METHODS
Green ash (Fraxinus pennsylvanica M.) trees at two sites that were first infested with emerald ash borer (EAB) for at least 2 to 3 years prior to this study were selected for this study in northeast Ohio. Although EAB was present at both sites, ash trees were not uniformly affected and had lost varying amounts of canopy at each of the two sites. Trees at either site were grouped into three groups of 10 trees each, using canopy loss as the main criteria to distinguish the effects of EAB presence: group I – 0 to 5%; group II – 6 to 25%; and group III – > 40%; trees with canopy loss between 25 and 40% were not considered in this trial. The trees in each of the groups occurred in a randomly mixed arrangement within the tree populations at either site.
Green ash trees at the first site (Shalersville, OH: 41.2469° N Latitude, 81.2299° W Longitude) were planted from nursery grown seed stock approximately 30 years ago at site on 4.5 m (14.7 foot) centers (Davey Tree Experts Research Station, Shalersville, OH). Trees at this site ranged from 22.4 to 45 cm (8.8 to 17.7 in) in stem diameter at breast height (dbh).
Green ash trees at the second site grew on the fringe of a pasture in a naturally seeded woodlot (Streetsboro, OH; 41.2392° N Latitude, 81.3459° W Longitude). Trees were estimated from landowner records to be 12 to 15 years old and were varying distances from each other, from 5 m to 6 m (16.4 feet to 19.6 feet). These were smaller trees compared to site one and ranged from 13.5 to 17.2 cm (5.3 to 6.7 in) in stem dbh.
Three methods were used to compare the structural integrity of the three ash tree groups. 1) A stereoscopic photogrammetry digital image correlation (DIC) system was used to evaluate and compare strains on ash branches that were (static) loaded to failure; 2) a Fakoop microsecond timer was used to compare transmission times (Tm μs/m) of sound waves radially through the pith of stem wood at seven heights; 3) the root plate and root balls of ash trees in each of the three groups were examined for fine root necrosis (root necromass/biomass ratios) and percentage of the surface area of the root ball that was decayed, respectively.
DIC Analysis of Ash Branches with EAB Infestation
Strain (ε, relative deformation or change in length of a material, Dahle et al 2017) can be evaluated using a number of standard techniques, which may include the use of strain gauges (example Clairie et al. 2013) or transducers (Brudi and van Wassenaer 2002; Clarie et al. 2006; James and Kane 2008; Kane and Clouston 2008). A more in-depth review of measuring strain in trees can be found in James et al. 2018. Digital image correlation (DIC) algorithms have been used in commercial systems to optically measure deformation; Sebera et al. (2014) used the VIC-3D (Correlated solutions Ltd.) system, while the ARAMIS 3D (GOM GMBH) has been utilized to measure ε in wood (e.g., Oscarsson et al. 2014; Dahle et al. 2017). Percent ε = (L0 − L1) / L0 × 100 is calculated by the ARAMIS 3D DIC software, which tracks the spatial movement of the speckles in reference to each other (Windisch et al. 2007; Tanasic et al. 2012), where L0 is the length before and L1 length after loading. Unlike strain gauges that can only measure strain at a specific location, DIC can map strain ε across the area of interest at the sample, even in the case of irregular geometry, which is often the case in trees. The advantage of stereoscopic photography is that 3D DIC systems can compute displacements and strains on curved surfaces, compared to monoscopic systems which employ a single camera and only compute in-plane deformations (e.g., Windisch et al. 2007; Tanasic et al. 2012; Dahle et al. 2017). A stochastic paint speckling is applied to create a unique series of reference points that can be tracked during testing. This is accomplished by applying black speckling, at the desired size, on top of a white background (Dahle et al. 2017). In this trial, the ARAMIS 3D DIC system was placed on a sturdy scaffolding system which allowed for height placement of the camera to align with the painted branch union (of the branch that was subsequently loaded to failure). We utilized 17 mm lenses, and the DIC system was calibrated to 13 paired images of a 350 × 280 mm2 panel to develop a 3-dimenstional region for testing. Images were acquired at 2448 × 2050 pixels with a resolution of 72 pixels per inch. The working distance was 270 cm with a calibration deviation of less than 0.3. The facet size was set to 20 × 20 pixels with a step size of 15. Linear strain was computated with a 7 × 7 field with a validity quote of 55.5%. The intersection deviation was maintained below 0.03 which provided an accuracy of 0.03 ± 0.003%. A branch replicate of each of a subset of five trees in each ash tree group was selected for testing based on a defect-free union, no apparent loss of bark over more than 10% of the branch, and an attachment angle within 35 to 60 degrees with the stem.
The goal of the static loading was to induce failure and not measure a specific moment arm. The static loading was conducted by attaching a load line to the approximate center of gravity of the selected branch (Persad et al. 2013; Dahle and Grabosky 2012.) The load was applied from a 2 ton battery operated winch until failure occurred and the DIC cameras were focused on the speckled branch section from a side view (Figure 1). During loading and upon ultimate fracture, the strain map was captured in the DIC system at 3 frames per second, and strain maps for each static loading event were captured and analyzed. Major strain (tensile) was measured on the stage just prior to failure (0.33 seconds) at the failure location along the top portion of the branch in the axial direction (Figure 2).
Acoustic Measurements “Transmission Time (Tm)”
Stress wave propagation in standing trees is utilized in forest harvest wood management programs, and the Fakopp 2D Microsecond Timer (FAKOOP Enterprise, H-9423 Agfglva, Fenyo. Str 26 Hungary) used in our study is one standard tool in common use (Wang et al. 2004b) to determine wood quality in standing timber.
Two needles connected to the Fakoop Timer are driven into the bark of the standing tree on opposite sides of the stem and reach the sapwood of the tree; accelerometers connected to the needles pick up a start and stop signal. A hammer is used to tap one of the sensor needles (start sensor) which generates a stress wave into the tree stem in a radial direction. The opposite sensor works in tandem as a receiver sensor and picks up the stop signal displaying the radial stress wave velocity on an LCD screen (Wang et al. 2004a). The measured transmission times (Tr) values were converted to radial stress wave transmission times (Tm) (μs/m) (example, Mattheck and Bethge 1993; Wang et al. 2004a). In our study, needles were inserted on opposite sides and into the sapwood of the base of each tree and examined from 5 cm (2 in) from the soil line upwards in 15 cm (6 in) increments until 7 readings were obtained (final height near 95 cm [37 in] or close to a meter from the soil line). Five trees in each of the groups were evaluated at both sites. A forestry caliper tool was used to measure the distance (D mm) between the needles at each reading interval, which allowed for the conversion from the transmission time readings (Tr μs) to transmission times (Tm μs/m) using the formula: Tm= Tr/D as described by Wang et al. 2004b.
Root Plate Evaluation
Root Plate Necromass/Biomass (Nm/Bm) Ratios: Soil Cores and Fine Root Analysis to Determine (Nm/Bm) Ratios
Root Plate Necromass/Biomass (Nm/Bm) ratios represent the dead root mass divided by the amount of living root mass determined after examination under the dissecting microscope of fine root tissue taken from the root zones of each replicate tree. Fine roots from trees are pervasive in the soil in and around the root plate zone and are useful indicators of health, as fine root production improves in suitable site conditions and with tree health (Fontaine et al. 2003). Soil cores of fine roots in this trial were sampled and used as a representation of fine root production under the relatively uniform conditions of soil and topography inherent at the site. Soil cores were taken using a bulk density/slide hammer at the base of each ash tree replicate outwards in five concentric circles. Each progressive circle was a distance equivalent to the stem diameter at breast height further out than the preceding one. Six cores or soil/root samples were randomly taken on each circle, thus thirty cores were collected for each of ten ash tree replicates of the three groups of trees. Core samples were taken to a depth of 20 cm (7.8 in) and individual major roots that cores were extracted from were traced out with an air spade (365 cubic feet per minute) to ensure that samples were taken from correctly identified trees, as roots from adjacent trees often form interlocking root mats. Fine roots in the root cores were extracted from the soil in the laboratory and separated out using a wire mesh sieve (100 nm) and washed gently. The fine roots up to 2 mm were paper dried and were examined under the dissecting scope; each strand was cut longitudinally along the length using a scalpel to identify living and non-living tissue (< 1 mm up to 2 mm). After separating the living roots and dead fine roots, each sample was oven dried at 65° C (149° F) for 24 hours to constant weight (ash/debris excluded, Persson 2012). The dry weights (g) of total fine roots, living fine root tissue (biomass, Bm), and of non-living fine root tissue (necromass, Nm) observed for each core was recorded (Persson and Stadenburg 2010). Comparison of the Nm/Bm ratios were made for 30 soil cores for each ash tree replicate in each of the three ash tree groups.
Root Ball Decay Evaluation
An evaluation of the extent of decay in the root balls of the ash trees in each of the three groups of trees was conducted to investigate root ball decay (if evident) with the extent of canopy loss. The soil surface area around the base of the stem of all trees was cleared of any vegetation or top-laying rocks. Root balls were excavated by first cutting a circle (circumference) with radius that was three times the diameter of the stem at breast height using a diamond edged blade mounted on an industrial concrete cutting saw (48 cm blade). The blade was plunged into the soil and a circular cut was made along the marked circumference. Although at this site the root plates were inundated with small rocks and compacted soil, the full blade depth of 48 cm was achieved for all root balls. The root balls within this volume of earth were excavated using a track Excavator and 48 cm bucket (Mark J. Hoenigman, Busy Bee Services Ltd. Novelty, Ohio). The root balls were power-washed with tap water to remove encrusted soil and air dried. Pockets of decay, if evident on major basal roots (above 2 cm in diameter), were bark traced or the bark around the decay was cut to fully reveal the decay pockets. Percent surface area affected was estimated using a wax paper template placed tightly over the decay zone, marking and measuring the surface area. The overall percentage decay of the root ball was estimated from the ratio of the decay pockets to that of the entire surface area (which was obtained by similarly measuring and recording the surface area of the entire excavated root ball and major roots present; roots < 2 cm were not assessed).
STATISTICAL METHODS
All analyses were conducted in SAS (version 9.1 SAS Institute, Cary, NC). Statistical comparison of means between tree groups were made using the Scheffes test (alpha value = 0.05).
RESULTS AND DISCUSSION
DIC Analysis of Ash Branches with EAB Infestation
Strain at the point of failure after static loading branches to failure was significantly higher (or the change in length upon loading was higher) in group I trees compared to that observed for group II and group III, which were not significantly different from each other (Figure 3; P = 0.0059, F = 8.13, DF 2. 14; Scheffes test for Sq Rt Perc. SAS). These data indicate that as canopy loss progresses after infestation with EAB, the strain mapped by DIC decreased significantly with canopy loss (from 6% canopy loss and greater). A decrease in strain could be indicative of less stress required for failure in these branches. Unfortunately, we were not able to measure the applied loading in this study, thus we could not determine stress induced during the testing. Consequently we cannot compare stress and strain at the point of failure.
The incorporation of DIC into this arboricultural research trial is a novel way of using DIC to investigate strains and the effect of a biological entity on wood structural integrity. Further work in utilization of DIC in evaluating how insect and disease pests and other biotic factors may impact tree structure is encouraged. Additionally, the DIC system and strain mapping allows for further understanding of how plants may transfer loading or react during static or dynamic loading (wind, rain, snow, and ice) from one section of tree to another (i.e. from a branch to the central stem) and can help build our arborist safety programs and give insight to tree stability in urban settings.
Acoustic Transmission Time (Tm)
Acoustic transmission time (Tm μs/m) observed at site one was not statistically different between groups I and II; group III trees, however, with the most canopy loss had statistically longer Tm μs/m times compared to groups I and II (Figure 4) (P = 0.0034, F = 9.13, DF 2. 14; Scheffes test for Sq Rt Perc. SAS). Transmission times at site two (smaller trees) increased with canopy loss and was significantly different between all three groups (P = 0.0052, F = 8.53, DF 2. 14; Scheffes test for Sq Rt Perc. SAS). The mean transmission times obtained in this study for the group I trees are consistent with the range (725 to 861 μs/m) recorded by Mattheck and Bethge (1993) and Smulski (1991) (645 to 845 μs/m) for healthy standing ash trees. Lower transmission time (higher sound velocity) for group I could be also explained by higher moisture content. Persad et al. (2013) found that trees with similar early canopy loss ranges had significantly higher moisture compared to trees with advanced canopy dieback.
The overall data at both sites indicate that as canopy loss progresses, the transmission times of acoustical stress waves that are perpendicular to the wood grain increase. Pellerin et al. (1985) estimate a 50% decrease in wood structural strength with a 30% increase in transmission time. This study has an average of 33% increase in Tm between group I and group III trees, which may suggest an estimated > 50% in loss of stiffness between trees with < 5% canopy loss compared to trees with > 40% canopy loss. These data are consistent with the findings of an earlier study on static loading (Persad et al. 2013) on decrease in structural strength associated with age of EAB infestation in ash trees.
Root Plate Necromass/Biomass (Nm/Bm) Ratios and Root Ball Decay Evaluation
Root Plate Necromass/Biomass Ratios (Nm/Bm)
From an analysis of Nm/Bm ratios of fine roots observed from soil cores sampled in the three groups of EAB-infested ash trees at site two, group III trees had significantly higher amounts of dead fine roots (Figure 5) compared to groups I and II, which were not significantly different from each other (P = 0.004, F = 9.30, DF 2. 29; ANOVA Scheffe test, SAS). Fine roots die off for a number of reasons, including reduced carbohydrate supplies, drought, and stresses associated to changes in the soil profile (Reynolds 1970; Marshall and Waring 1985). Reduction in carbohydrate supply to the roots in group III trees may be associated with EAB activity and canopy loss (> 40%), as groups I and II, which had higher Nm/Bm values, were in the same woodlot and were spread randomly amongst the group III trees within the site. The loss of fine roots may also mean the inability to osmo-regulate effectively and could contribute to decay in larger roots. These data that align fine root death with canopy loss can help in tree preservation guidelines, as they may determine if sufficient fine roots are present to facilitate uptake of water and soil applied systemic insecticides and or fertilizers.
Root Ball Decay Evaluation
Percentage decay in the extracted root balls was significantly different in all groups, with group I trees having the least decay and group III having the highest percentage decay (P = 0.0002, F = 11.64, DF 2. 27; Scheffes test for Sq Rt Perc. SAS) (Figure 6). Root decay pockets may arise as root tissue is carbohydrate starved and/or is unable to osmo-regulate effectively and may result in water or moisture accumulation and an increase in pathogen activity. These findings on decay incidence with canopy loss can be used in risk evaluation and also provide additional criteria for tree preservation efforts, as root tissue that is compromised may lead to insufficient uptake of water and systemic insecticides resulting in lower efficacy on EAB larvae.
CONCLUSION
The tools used in this trial, the DIC system, stress waves transmission times, and root evaluations of ash trees infested by EAB, provided insight into variation of material properties related to canopy loss. Trees in groups II and III exhibited significantly lower strain during static loading to branch fracture compared to that observed for trees in group I. A decrease in strain could be indicative of less stress required for failure in these branches.
Analysis of stress wave transmission times revealed that at site one (larger trees), group III trees with the highest canopy loss had significantly higher Tm μs/m times compared to the other two groups. Transmission times increased significantly with advanced canopy loss at site two (smaller trees) between all groups. The use of stress wave technology in evaluating ash tree condition for preservation and other purposes may thus be more apparent in smaller sized trees. Fine root analyses of the root zone areas (site two) revealed that the Nm/Bm fine root ratios were significantly different between all groups, with group III trees having highest Nm/Bm ratios and group I trees having the lowest. In an evaluation of decay in extracted root balls (site two), group III trees had significantly higher percentage surface area that was decayed compared to groups I and II, which had statistically similar decay incidence. The loss of carbohydrate resources to root tissue with increasing canopy loss may be a factor in the root necromass and decay incidence observed. These data provide fundamental insight into material properties of ash trees after infestation by EAB and can contribute to arboricultural guidelines for ash tree preservation and help develop safety protocols to address structural loss in trees after EAB infestation.
ACKNOWLEDGEMENTS
The authors appreciate the help of the volunteers of Biomechanics week 2013. Mark J. Hoenigman assisted with the root ball extractions. We also thank Mike Hannebique, Katie Harper, Josh Gula, and John Siefer for help with various equipment and fieldwork throughout the trials.
- © 2019, International Society of Arboriculture. All rights reserved.