Genetic variation in wood properties of mid-rotation age Eucalyptus globoidea Ebenezer A. Iyiola, Clemens M. Altaner* and Luis A. Apiolaza School of Forestry, University of Canterbury, Christchurch 8140, New Zealand *Corresponding author: clemens.altaner@canterbury.ac.nz (Received for publication 13 July 2021; accepted in revised form 26 March 2022) Abstract Background: Eucalyptus globoidea Blakely produces ground-durable (Class 2) and stiff wood and has the potential to be grown in New Zealand to supply high-value environmentally-friendly timber for use as posts in the agricultural sector and stiff veneers for the LVL industry. The New Zealand Dryland Forests Initiative (NZDFI) has established a breeding programme for this species. The objective of this study was to identify trees with superior wood properties for commercial propagation enabling the establishment of a domestic plantation resource of ground-durable timber. Methods: The genetic variation in wood properties at mid-rotation age (8-year-old) of 141 E. globoidea families was assessed for the following traits: heartwood diameter, diameter under bark at ~0.5 m height, combined sapwood diameter, heartwood collapse, sapwood collapse, standing tree acoustic velocity and extractive content in the heartwood. Families were ranked and genotypes with large heartwood diameter, high extractive content and stiffness as well as low collapse were identified. Results: Heartwood diameter (h2 = 0.51) and extractive content (h2 = 1.16) showed good heritability, which in combination with high variation are promising traits for a breeding programme. The high heritability for extractive content indicated a closer relatedness within the population than the assumption of unrelated families of half-siblings. The unfavourable correlation between the heartwood diameter and extractive content (genetic correlation (rg) = −0.45) indicated that a compromise is required for simultaneous improvement of both traits. Heritability estimates for heartwood collapse (h2 = 0.30) and acoustic velocity (h2 = 0.36) were moderate. Conclusions: Genetic selection for wood quality traits of E. globoidea is practically feasible. Superior individuals with above average performance for multiple traits were present in the breeding populations, however, this was dependent on the intended end use of the timber. New Zealand Journal of Forestry Science Iyiola et al. New Zealand Journal of Forestry Science (2022) 52:13 https://doi.org/10.33494/nzjfs522022x172x E-ISSN: 1179-5395 published on-line: 27/04/2022 © The Author(s). 2022 Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Research Article Open Access it has high stiffness and strength, and its heartwood is rated Class 2 for in-ground durability, i.e. lasting 15 to 25 years in service (AS5604 2005). The New Zealand Dryland Forests Initiative (NZDFI) is focusing on establishing a sustainable plantation resource of naturally durable eucalyptus timbers (Millen et al. 2018). Five eucalyptus species have been selected for their breeding programme, including E. globoidea, for the establishment of a domestic naturally ground-durable hardwood industry. Such a sustainable domestic resource is an environmentally friendly high- Introduction Eucalyptus globoidea Blakely, commonly known as white stringybark, belongs to the subgenus Eucalyptus (Brooker & Kleinig 1983), and is native to south-eastern Australia’s gentle undulating hill country near the coast and mountain slopes, extending to escarpments adjacent to tablelands but not inland of the ranges (Boland et al 2006). It is a species that displays several appealing properties such as early heartwood formation, drought, frost and pest-tolerance, and grows well in the New Zealand environment (Millen et al. 2018). Moreover, Keywords: Acoustic velocity; collapse; extractives; genetic gain; heartwood; natural durability; tree breeding; white stringybark http://creativecommons.org/licenses/by/4.0/), Iyiola et al. New Zealand Journal of Forestry Science (2022) 52:13 Page 2 value alternative not only to copper chromium arsenic (CCA) treated radiata pine, which is difficult to dispose of (Love 2007; Townsend & Solo-Gabriele 2006) but also naturally durable timber imports, which often originate from endangered or unsustainably harvested tropical sources (UNEP 2012). Envisaged uses for E. globoidea timber are primarily posts for the agricultural sectors (Millen et al. 2018) but it can also be used for producing stiff veneers for laminated veneer lumber (LVL) (Guo & Altaner 2018) and sawn timber (Jones et al. 2010). For the latter, stiffness and drying without collapse are relevant wood properties (Bootle 2005; Poynton 1979; Somerville & Gatenby 1996). Collapse is wood shrinkage caused by the buckling of the cell walls and flattening of the cell lumen during drying due to negative pressure (Chafe et al. 1992). It is different from normal shrinkage, which is caused by removal of bound water from the cell wall below the fibre saturation point (FSP). Naturally durable timber has less pronounced and more variable durability than treated wood. The durability varies among trees due to a combination of genetic and environmental factors (Zobel & Jett 1995). Genetic variation can be exploited in a breeding programme, while environmental variation can be addressed by appropriate siting, silviculture and timber processing. There is also variability within trees, since sapwood is never regarded as naturally durable, and typically durability of heartwood decreases towards the pith (AS5604 2005). As a consequence, the heartwood of younger trees tends to have lower durability than heartwood from older trees (Amusant et al. 2004). Therefore, to produce ground-durable posts from young trees grown in short rotation plantations, it is necessary to deploy genetically improved tree stocks. Young trees producing durable timber have been reported in breeding trials of E. cladocalyx (Bush et al. 2011). Assessing wood durability is time consuming and expensive. On the other hand, breeding programmes benefit from large sample sizes (Raymond 2002). Consequently, there has been little focus on natural wood durability in tree breeding programmes. Natural durability assessments in breeding trials, including both sampling heartwood from the trees and assessing its durability, need to be rapid and cost effective. Felling trees and sawing wood samples is destructive and too laborious for routine use in a breeding programme. A rapid means to obtain a heartwood sample without damaging the tree is taking increment cores (Estopa et al. 2017; Jones et al. 2008). Direct assessments of natural durability by mass loss of wood samples when exposed to fungi (AWPC 2007; CEN 1989; EN-350-1 1994) are impractical in a breeding context due to time and cost (Bush & Walker 2011; Harju & Venäläinen 2006; Pâques et al. 2013). Extractive content in the heartwood are responsible for the natural durability of wood (Bush et al. 2011; Hawley et al. 1924) and have been reported to be under genetic control (Bush et al. 2011; Li et al. 2018). Near- infrared (NIR) spectroscopy has been proven effective in predicting the extractive content in wood (Bush et al. 2011; Li et al. 2018; Stackpole et al. 2011). NIR spectra can be acquired within seconds from a solid wood surface and are an inexpensive method for assessing the chemical components of wood (Estopa et al. 2017; Greaves et al. 1996; Tsuchikawa & Schwanninger 2013). NIR calibrations for the extractive content were reported to have a root mean square error (RMSE) of 0.91% to 1.16% for E. bosistoana F.Muell. (Li 2018), 0.62% for Sweitenia marcrophylla King (da Silva et al. 2013), and 0.66% for E. globoidea (Kuwabara 2017). This study examined diameter under bark, heartwood diameter, combined sapwood diameter, extractive content (EC), sapwood collapse, heartwood collapse (all at 0.5 m stem height) and acoustic velocity (a proxy for wood stiffness) of E. globoidea. Breeding values of these traits were calculated for 141 half-sib families of 8-year- old E. globoidea. The phenotypic and genetic correlations between the wood traits were estimated. Materials and Methods The Atkinson trial established in November 2011 at Wairarapa, New Zealand (latitude 41° 35' 21" S and longitude 175° 24' 4" E), with 141 E. globoidea half- sib families was used in this study. The families in this trial were grown from seed collected from across the natural range of the species in Australia and from three NZ plantation sites with a known seedlot. The unknown demarcation of the seed collected in Australia means that provenance is undeclared. The trial contained 8640 trees and was laid out in an incomplete block design with each family replicated 40 to 80 times. There were 36 trees per block, planted at a spacing of 2.4 m between rows and 1.8 m within rows. Trees in the trial were assessed for diameter at breast height (DBH) and form in 2015 at age 3.4 years, pruned in August 2018 and thinned in April 2019 to remove trees with poorer stem form. All living trees not marked for removal and with a DBH larger than 30 mm were sampled with a purpose-built corer. A bark-to-bark 14 mm diameter core including the pith was extracted ~0.5 m above the ground from 2160 trees in March 2019. Heartwood quantity and quality Stem diameter under bark and heartwood diameter at 0.5 m height were measured on the green cores after highlighting heartwood with an aqueous 0.1% solution of methyl orange. Combined sapwood diameter was calculated as the difference of the two measures. Heartwood quality was assessed by predicting ethanol soluble extractive content (%) of dry cores from NIR spectra (Tensor 37, Bruker, Germany) obtained with a fibre optic probe at wavelengths from 9000 to 4000 cm-1 at 4 cm-1 intervals as described earlier (Li & Altaner 2019). Tangential collapse After the NIR measurements, the cores were equilibrated to a stable moisture content at 60% relative humidity and 25°C. Ten cores were randomly selected, and the widest tangential diameter was measured with a Vernier calliper. These 10 values were averaged and used as reference (D1) for the collapse assessment. Collapse was assessed separately for heartwood and sapwood (Equation 1). The narrowest tangential diameter of each core was measured in the sapwood (D2), as well as the two narrowest tangential diameters for the heartwood. Two diameters were measured in the heartwood as collapse was more prominent in this region and heartwood is the main target product for durable timber production. The latter were subsequently averaged. Tangential collapse (%) = ((D1–D2)/D1) × 100 (1) Standing tree acoustic velocity Ten trees were randomly sampled per family for those families with more than ten living trees. For families containing fewer than 11 standing trees, all trees in the family were sampled. Acoustic velocity was measured on 1,147 trees with the standing tree time-of-flight TreeTap tool (University of Canterbury, Christchurch, New Zealand) (Toulmin & Raymond 2007). A total of eight measurements were averaged per tree, which were acquired by placing the probes at one stem location; probes were spaced 1 m apart centered about breast height. Statistical analyses The data were analysed using R statistical software (R Core Team 2019). A linear mixed model was used for the univariate analyses (heritability) using Equation 2. yijk = μ + fi + bj + eijk (2) where yijk is a phenotypic observation of a single trait, µ is the overall intercept, fi is the random effect of the ith family, bj the random effect of the j th block and eijk the random residual for the kth individual of the ith family in the jth block. The residuals were assumed to be identically and independently normally distributed. The expected value for the phenotype is µ while the variances for each of the random terms are: V(f) = σ2f , V(b) = σ 2 b and V(e) = σ2e. The notation can be expanded to a bivariate scenario using vectors that have two phenotypes for traits 1 and 2 for each individual. The variances are then Kronecker matrix products V(f) = F0 ⊗ Zf , V(b) = B0 ⊗ Zb , and V(e) = R0 ⊗ I where Zf , Zb and I are incidence matrices, connecting the observations to model terms and F0, B0 and R0 are illustrated below. (3) where σ2f1, σ 2 b1, and σ 2 e1 are the variances for family, block and residual; σf12, σb12, and σe12 represent the covariances for family, block and residual between traits 1 and 2. The model was fitted with the Asreml-R package (Butler et al. 2009) to estimate variance components and breeding values. The phenotypic and additive genetic variations were estimated to compute the narrow sense half-sib heritability (h2) of each trait using Equation 4. h2 = Var(A)/Var(Y) = 4σ2f /(σ 2 f +σ 2 b + σ 2 e) (4) The heritability estimated in this study assumed that families were true half-siblings with a relationship coefficient within families of one quarter. The coefficient genetic of variation (CGV) for each trait was determined using the equation below. CGV = (√4σ2f / population mean) × 100 (5) The coefficient of phenotypic variation (CPV) for each trait was determined using the equation below. CPV = (standard deviation/population mean) × 100 (6) The genetic correlation is a measure of the strength of the genetic association between the performance in one trait and performance in another trait (Bourdon 2000). rg(m,n) = σ 2 fm / √ σ 2 fm σ 2 fn (7) is the covariance between traits m and n; σ2fm is the genetic variance for trait m; σ2fn is the genetic variance for trait n and rg(m,n) is the genetic correlation. The Pearson correlation coefficient was used to express the phenotypic correlation (rp) between the traits. Results and Discussion Collapse, acoustic velocity, and heartwood properties The descriptive statistics, heritability estimates, coefficients of phenotypic variation (CPVs) and coefficients of genetic variation (CGVs) for the assessed wood traits are presented in Table 1. High coefficients of variation were observed for extractive content, combined sapwood diameter, heartwood diameter, diameter under bark and sapwood collapse. Combined sapwood diameter exhibited larger variation than diameter under bark and heartwood diameter. The 8-year-old E. globoidea trees had produced more heartwood (mean heartwood diameter 91.0 mm) and had a lower coefficient of phenotypic variation (CPV = 28.9%) than has been reported for E. bosistoana (31.7 mm to 42.1 mm with a CPV of 51% to 61%) at 7 years of age (Li et al. 2018). In 9-year-old E. globulus Labill, heartwood Iyiola et al. New Zealand Journal of Forestry Science (2022) 52:13 Page 3 TABLE 1: Description of the study sites diameter ranged from 61 to 98 mm with a CPV of 6.8% to 45.1% (Miranda et al. 2014). The mean heartwood diameter reported for 8-year-old E. cladocalyx F.Muell ranged from 58 to 82 mm with a CPV of 29% based on samples collected from three different locations (Bush et al. 2011). The mean combined sapwood diameter was 50.9 mm with a CPV of 31.8% (Table 1). The average value reported in this study was lower than that reported for 7-year-old E. bosistoana (60.4 mm and 70.3 mm) (Li et al. 2018). This confirmed the general observation that E. globoidea has a narrow sapwood band (Bootle 2005). A narrow sapwood band is beneficial for post processing operations as less material needs to be removed. The CPV for E. globoidea (31.8%) was slightly higher than for E. bosistoana (26%) (Li et al. 2018). The sapwood diameter reported for 9-year-old E. globulus ranged from 16 to 29 mm (Miranda et al. 2014) and was 50 to 74 mm for 8-year-old E. cladocalyx (Bush et al. 2011). The mean predicted extractive content was 9.4% with a CPV of 47.9%. This was comparable to 7-year-old E. bosistoana for which 7.7% and 9.6% extractive content was observed at two different sites, with a CPV of 46% (Li et al. 2018). However, the extractive content of E. globoidea was lower than the average 12% observed for similarly-aged E. cladocalyx (Bush et al. 2011). Tangential collapse in the heartwood had a mean of 18.2%, similar to the value of 18.6% reported for the collapse prone E. nitens (H.Deane & Maiden) Maiden (Kube 2005) and E. dunnii Maiden (14.9%) (Arnold et al. 2004). However, care must be taken when comparing the mean values reported for collapse from different studies, as the measurement techniques differed. This study measured the maximum tangential dimensional change, a combination of normal shrinkage and collapse at the worst part of the sample, which should be related to checking. The coefficient of variation in collapse reported in this study was high (37.91%) and corresponded to the observed value (37%) for 6.5-year-old E. dunnii (Arnold et al. 2004). Collapse in the sapwood was lower (3.2%) but more variable (CPV 168.75%) compared to collapse Iyiola et al. New Zealand Journal of Forestry Science (2022) 52:13 Page 4 in heartwood (18.2%; CPV 37.91%). It is thought that collapse is caused by negative pressure generated over gas–liquid surfaces with high curvature when water evaporates from the wood. Extractives can reduce pore size/permeability of cell walls and therefore increase the likelihood of collapse in heartwood (Chafe et al. 1992). The mean standing tree acoustic velocity ranged from 2.13 to 4.27 km/s with a mean of 2.96 km/s and a CPV of 10.47% (Table 1). It is not straightforward to compare standing tree acoustic velocity measurements with other studies, as commonly used acoustic tools appear to provide different absolute values (Dungey et al. 2012). The standing tree acoustic velocity (TreeTap) of 2.96 km/s was similar to the 2.5 km/s (ranging from 2.2 to 2.8 km/s) reported for the species at age 25 years using the IML hammer (Wiesloch, Germany) (Jones et al. 2010). This is comparable with 2.7 km/s in 10-year- old Pinus radiata D.Don in Australia measured with TreeTap (Toulmin & Raymond 2007). It should also be acknowledged that age has a significant effect on the acoustic velocity of young trees, as the microfibril angle (MFA) decreases with cambial age (Lachenbruch et al. 2011) and that the mean acoustic velocity is generally lower in softwoods than hardwoods due to their higher MFA (Lindström et al. 2002). Comparable CPVs for acoustic velocity were reported in the literature; 7.87% for 8-year-old E. fastigata H. Deane & Maiden (Suontama et al. 2018), 8.1 to 10.2% for 20-year-old Pseudotsuga menziesii (Mirb) Franco (Klápště et al. 2019) and 9.05% for P. sylvestris L. (Hong et al. 2015). An even lower CPV of 2.2% was reported for 10-year-old P. radiata (Toulmin & Raymond 2007). Heritability and coefficient of genetic variation of studied wood properties This study used a coefficient of relatedness (rc) of half- siblings (rc = 0.25), i.e. the families were assumed to be perfect half-siblings and unrelated to each other (Table 1). Heartwood diameter had a heritability estimate of 0.51 and a CGV of 20.7%. Heritability reports for TABLE 1: Descriptive statistics, heritability (h2) with the 95% confidence interval in brackets for E. globoidea wood properties at age 8 years; coefficient of phenotypic variation (CPV) and coefficient of genetic variation (CGV); rc: coefficient of relatedness. Trait Mean Standard deviation Min Max CPV (%) CGV (%) h2 (rc = 0.25) Diameter under bark (mm) 141.9 33.7 60 255 23.75 19.59 0.67 (0.49, 0.85) Heartwood diameter (mm) 91.0 26.3 0 190 28.90 20.70 0.51 (0.36, 0.66) Combined sapwood diameter (mm) 50.9 16.2 6 150 31.83 25.44 0.63 (0.45, 0.80) Predicted extractive content (%) 9.4 4.5 −4.4 31.7 47.87 51.58 1.16 (0.90, 1.39) Acoustic velocity (km/s) 2.96 0.31 2.13 4.27 10.47 6.24 0.36 (0.18, 0.54) Heartwood collapse (%) 18.2 6.9 −2.9 45.5 37.91 19.16 0.30 (0.17, 0.40) Sapwood collapse (%) 3.2 5.4 −14.1 25.1 168.75 56.52 0.12 (0.03, 0.21) heartwood diameter ranged from 0.39 to 0.61 for E. grandis W.Hill and E. cladoclayx (Bush et al. 2011; Santos et al. 2004), comparable with the result obtained in the current study. The observed heritability and variation in E. globoidea are suitable for increasing heartwood production through selective breeding; however, the coefficient of variation may decline with tree age (Bush et al. 2011). The heritability estimates for combined sapwood diameter and diameter under bark were 0.63 and 0.67 with CGVs of 25.44% and 19.59%, respectively. Similar heritability estimates (0.67 to 0.82) were reported for sapwood diameter in 7-year-old E. bosistoana (Li et al. 2018) and Larix kaempferi (Lamb.) Carrière (Pâques 2001). The heritability of diameter under bark was comparable to the range (0.29 to 0.72) reported for diameter at breast height (DBH) in previous studies (Kube 2005; Whiteman 1992) and the CGV (19.59%) for diameter under bark was at the lower end of the range (15.1% to 87.1%) reported for growth of Eucalyptus hybrid clones (Wu et al. 2017). The heritability estimate for extractive content was 1.16 and had a CGV of 51.58%. The relatedness coefficient (0.25) used in this study pushed the heritability estimate above 1, potentially indicating deviations from half- siblings and inbreeding effects (Eldridge et al. 1994; Elliott & Byrne 2003). In general, heritability estimates reported for extractive content ranged from 0.19 to 0.56 for E. bosistoana and E. globulus (Li et al. 2018; Poke et al. 2006). The estimates of heritability in this species might be too high because of the unknown relatedness and the fact that the assessments were restricted to one site (White et al. 2007). The high CGV (51.58%) combined with strong genetic control indicated potential to improving natural durability of E. globoidea through selective breeding. Collapse was moderately heritable in the heartwood (h2 = 0.30) but was not very variable (CGV of 19.16%), while it had low heritability in the sapwood (h2 = 0.12) where it was more variable (CGV = 56.52%). The narrow sapwood band in E. globoidea is unlikely to be converted into solid wood products. The heritability, although low, along with the coefficient of variation estimates indicated potential for improving tangential collapse in E. globoidea heartwood through selection. The heritability of collapse in the heartwood was similar to most values reported for other eucalypts; E. nitens (0.23 to 0.61) (Hamilton et al. 2004; Kube 2005) and E. grandis (0.23 to 0.31) (Bandara 2006). The location where E. delegatensis F.Muell. ex R.T.Baker seed was collected was reported to influence the severity of the collapse, indicating some degree of genetic control (King et al. 1993). The heritability for standing tree acoustic velocity was 0.36 with a CGV of 6.24% for E. globoidea, similar to values reported in the literature: 0.26 to 0.75 for 11- to 13-year-old P. menziesii (Douglas fir) (Dungey et al. 2012; Klápště et al. 2019) and 0.35 for an E. fastigata progeny trial (Suontama et al. 2018). The coefficient of genetic variation reported in this study was low. Even lower values of genetic variability have been reported in the literature, for example 4.8% to 5.13% for P. menziesii and E. nitens, respectively (Blackburn et al. 2014; Klápště et al. 2019). The low genetic variability of acoustic velocity limits the potential for genetic improvement of stiffness. Phenotypic and genetic correlations between traits of 8-year-old E. globoidea The phenotypic and genetic correlations between the wood traits are presented in Tables 2 and 3. There were positive phenotypic correlations between diameter under bark and sapwood diameter (rp = 0.65) as well as heartwood diameter (rp = 0.88) (Table 2). Iyiola et al. New Zealand Journal of Forestry Science (2022) 52:13 Page 5 TABLE 2: Phenotypic correlations between traits for 8-year-old E. globoidea (95% confidence interval in brackets). Trait Diameter under bark Sapwood col- lapse Heartwood diameter Combined sap- wood diameter Extractive con- tent Acoustic veloc- ity Heartwood collapse −0.03 (−0.07, 0.02) 0.44 (0.41, 0.48) 0.03 (−0.03, 0.06) −0.09 (−0.13, −0.05) 0.11 (0.06, 0.15) −0.01 (−0.07, 0.05) Diameter under bark −0.10 (−0.14, −0.06) 0.88 (0.87, 0.89) 0.65 (0.62, 0.67) −0.31 (−0.35, −0.27) 0.15 (0.10, 0.20) Sapwood collapse −0.08 (−0.13, −0.05) −0.07 (−0.11, −0.02) −0.01 (−0.06, 0.03) −0.01 (−0.07, 0.05) Heartwood diameter 0.20 (0.17, 0.25) −0.14 (−0.18, −0.10) 0.14 (0.09, 0.20) Combined sapwood diameter −0.42 (−0.46, −0.39) 0.08 (0.02, 0.13) Extractive content −0.07 (−0.13, −0.02) The strong positive correlations indicated that larger trees have wider heartwood as well as sapwood. Similar phenotypic correlations were reported between DBH and heartwood diameter for 7-year-old E. bosistoana (rp = 0.59) (Li et al. 2018), 22.5-year-old E. tereticornis Sm. (rp = 0.79) (Kumar & Dhillon 2014), Acacia melanoxylon R.Br. (rp = 0.88) (Knapic et al. 2006) and 30- to 37-year-old plantation-grown Pinus radiata D.Don (rp = 0.71) (Wilkes 1991). The genetic correlation between heartwood diameter and diameter under bark was stronger than the phenotypic correlation (rg = 0.91, CI95% 0.87, 0.95, Table 3). Strong positive genetic correlations between these traits were also reported for 9-year-old E. globulus (rg = 0.99) (Miranda et al. 2014), 7-year-old E. bosistoana (rg = 0.89 to 0.98) (Li et al. 2018), E. cladocalyx (rg = 0.44) (Bush et al. 2011), L. kaempferi (rg = 0.87 to 0.92) (Pâques 2001) and 35-year-old Juglans nigra L. (rg = 0.98) (Woeste 2002). No genetic correlation (rg = 0.022) was reported between the two traits for 25-year-old P. sylvestris (Fries & Ericsson 1998). The growth traits were favourably correlated with heartwood diameter being the most relevant measure when aiming to produce naturally durable timber. A negative phenotypic correlation was observed between heartwood diameter and extractive content at both the phenotypic (rp = −0.14, Table 2) and genetic (rg = −0.45, Table 3) level. A strong negative genetic correlation (rg = −0.86) was also reported between these traits for 7-year-old E. bosistoana (Li et al. 2018). However, a positive genetic correlation between heartwood diameter and extractive content (rg = 0.32) was found for L. eurolepis (Pâques & Charpentier 2015). With the negative correlation observed for E. globoidea, trees with more heartwood tend to have lower amounts of extractives in the heartwood, resulting in lower decay resistance (Li et al. 2020). Selection for high durability will compromise heartwood quantity; however, as the correlation applies at the population level, individual trees that excel in both traits (i.e., so-called “correlation breakers”) may exist. The relationship between collapse in the sapwood and heartwood was positive at the phenotypic level (rp = 0.44, Table 2), and the genetic level (rg = 0.64, Table 3), suggesting that collapse is related to wood anatomy and amplified by heartwood extractives. No previously published results on the correlations between the two traits in other species were found for comparison. Heartwood collapse was positively correlated with extractive content at both the phenotypic and genetic level (rp = 0.11; rg = 0.23). A positive phenotypic correlation was reported between the two traits for E. regnans F.Muell. (Chafe 1987). The positive correlation between the traits is in line with current understanding of the causes of drying collapse, as heartwood extractives reduce the permeability of the cell walls and consequently increase the negative pressure during evaporation of water. However, as these correlations were not strong other factors, e.g. density, contribute to collapse. There were only weak correlations between diameter under bark and acoustic velocity at both the phenotypic and genotypic level (rp = 0.15; rg = 0.18, Table 2 and Table 3). Weak phenotypic relationships between these traits have also been reported for 8- to 25-year-old P. radiata (0.04 to 0.18) (Chauhan & Walker 2006; Toulmin & Raymond 2007) and 25-year-old E. dunnii (0.14) (Joe et al. 2004). In the current study, acoustic velocity was weakly correlated with heartwood diameter at the phenotypic level (rp = 0.14, CI95% 0.09, 0.20) but independent at the genetic level (rg = 0.10, Table 3). Therefore, our results indicate that the two traits need to be improved independently. Furthermore, the weak unfavourable correlations between acoustic velocity and extractive content at the genetic level (rg = −0.15) and phenotypic level (rp = −0.07) (Table 2) imply the need for compromise in selecting for the two traits simultaneously. Iyiola et al. New Zealand Journal of Forestry Science (2022) 52:13 Page 6 TABLE 3: Genetic correlations between the traits for 8-year-old E. globoidea (95% confidence interval in brackets). Trait Diameter under bark Sapwood collapse Heartwood diameter Combined sapwood diameter Extractive content Acoustic velocity Heartwood collapse −0.05 (−0.30, 0.22) 0.64 (0.36, 0.90) −0.09 (−0.36, 0.19) −0.02 (−0.25, 0.30) 0.23 (−0.01, 0.45) 0.05 (−0.39, 0.48) Diameter under bark 0.18 (−0.15, 0.50) 0.91 (0.87, 0.95) 0.82 (0.75, 0.89) −0.67 (−0.81, −0.54) 0.18 (−0.09, 0.45) Sapwood collapse 0.16 (−0.17, 0.51) 0.16 (−0.17, 0.52) −0.05 (−0.36, 0.29) 0.05 (−0.55, 0.62) Heartwood diameter 0.53 (0.35, 0.71) −0.45 (−0.62, −0.28) 0.10 (−0.19, 0.43) Combined sapwood diameter −0.79 (−0.89, −0.69) 0.24 (−0.01, 0.54) Extractive content −0.15 (−0.41, 0.12) Selection of superior genetics Breeding programmes typically aim to improve multiple traits simultaneously, either by selecting the best or culling the poorest performing families. In the absence of economic weights for the traits, superior genotypes can be selected by choosing families with above average breeding values in several traits. Approximately 15% of the 141 E. globoidea families showed both above average heartwood diameter and extractive content values (Figure 1). However, in a situation where heartwood collapse as a third trait is considered, only four families met the criteria (highlighted in red in the top right quadrant). If the interest is in using wood for laminated veneer lumber (LVL), stiffness, growth (stem diameter) and collapse are important traits that need to be considered (Figure 2). A total of 18 families highlighted in red in the top right corner of the quadrant met these criteria. Heritability estimates and breeding values are dependent on the relatedness of the individuals and families. Neither between family nor within family relatedness for the E. globoidea trees was known. The presented absolute values are likely overpredictions as no provenance effect was considered and families were not true half-siblings but a mixture of full-siblings and selfed individuals. Nevertheless, the ranking of the families by their breeding values is independent of their relatedness and therefore better performing genotypes can be selected by the industry. Conclusions The traits of E. globoidea investigated in this study were under varying degrees of genetic control and showed different degrees of variability. Extractive content, an indicator of natural durability, and heartwood diameter, are the two key traits determining the value of the species for the intended use as ground-durable timber, and showed the highest coefficient of phenotypic and genetic variation (CPV = 47.87% and 28.90%). Combined with the narrow sense heritability for heartwood diameter and extractive content of 0.51 and 1.16, respectively, significant genetic gain should be possible for these traits. Iyiola et al. New Zealand Journal of Forestry Science (2022) 52:13 Page 7 FIGURE 1: Relationship between family breeding values of heartwood diameter and EC for 141 E. globoidea families at age ~8 years. Families performing above average are located in the top right quadrant. Families with superior (red) and inferior (blue) average performance for heartwood collapse are highlighted. In accordance with the literature and the general understanding of the physical causes of collapse, heartwood was more prone to collapse than sapwood. Collapse in the heartwood was under moderate genetic control and this study showed that maximal tangential collapse can be included at reasonable cost into a breeding programme if core samples are available. Selection of superior genetics among the families would be practically feasible and could be used to increase utilisation of E. globoidea for solid timber products. Acoustic velocity was under moderate genetic control (0.36) with low variability, implying improvement of this trait through breeding would be challenging. Competing interests The authors declare no competing interests. Authors’ contributions EI carried out experimental work, data analysis and drafted the manuscript. CA developed the measurement methods and conceived the study. LA designed the trials and data analysis. All authors revised the manuscript. Acknowledgements The NZDFI facilitated this study by establishing and maintaining the breeding programme. We appreciate M Sharma, G Hendriks, M Holzenkämpfer, F Guo, L Nguyen, B Disney, A Millen, and R Sheppard for their assistance in sample preparation. Abbreviations CGV: coefficient of genetic variation CPV: coefficient of phenotypic variation DBH: diameter at breast height EC: extractive content LVL: laminated veneer lumber MFA: microfibril angle NIR: Near infrared spectroscopy NZDFI: New Zealand Dryland Forests Initiative Funding This project was funded by the MBIE Specialty Wood Products Partnership (FFRX1501 contract) and a University of Canterbury PhD Scholarship awarded to Ebenezer Iyiola. Iyiola et al. New Zealand Journal of Forestry Science (2022) 52:13 Page 8 FIGURE 2: Relationship between breeding values of diameter under bark and acoustic velocity for 141 E. globoidea families at age ~8 years with the families superior for tangential collapse highlighted in red. References Amusant, N., Beauchene, J., Fournier, M., Janin, G., & Thevenon, M.-F. (2004). Decay resistance in Dicorynia guianensis Amsh.: analysis of inter-tree and intra-tree variability and relations with wood colour. Annals of Forest Science, 61(4), 373-380. https://doi.org/10.1051/forest:2004030 Arnold, R.J., Johnson, I., & Owen, J. (2004). Genetic variation in growth, stem straightness and wood properties in Eucalyptus dunnii trials in northern New South Wales. Forest Genetics, 11(1), 1-12. AS5604. (2005). Timber - Natural durability ratings, Australian Standard, 29 p. AWPC. (2007). Protocols for assessment of wood preservatives. Clayton, Australia, The Australian wood preservation committee. 30 p. Bandara, K.M.A. (2006). Genetic improvement of solid wood product value of subtropical eucalypts: a case study of Eucalyptus grandis and E. dunnii. PhD, Australia National University, Australia. Blackburn, D., Hamilton, M., Williams, D., Harwood, C., & Potts, B. (2014). Acoustic wave velocity as a selection trait in Eucalyptus nitens. Forests, 5(4), 744-762. https://doi.org/10.3390/f5040744 Bootle, K. (2005). Wood in Australia - Types, properties and uses (2nd Edition ed.). McGraw-Hill Australia Pty Ltd. Bourdon, R.M. (2000). Understanding animal breeding (Vol. 2). NJ, USA: Prentice Hall. Brooker, M.I.H., & Kleinig, D.A. (1983). Field guide to eucalypts. Volume 1. South-eastern Australia. Melbourne, Australia: Inkata Press Pty Ltd. Bush, D., McCarthy, K., & Meder, R. (2011). Genetic variation of natural durability traits in Eucalyptus cladocalyx (sugar gum). Annals of Forest Science, 68(6), 1057-1066. https://doi.org/10.1007/ s13595-011-0121-z Bush, D., & Walker, J. (2011). Selecting and breeding eucalypts for natural durability. Developing a eucalypt resource: learning from Australia and elsewhere. Blenheim, NZ: Wood Technology Research Centre, pp. 125-136. Butler, D., Cullis, B., Gilmour, A., & Gogel, B. (2009). ASReml R package version 3.0. Hemel Hempstead, UK: VSN International Ltd. CEN. (1989). EN 252. Field test method for determining the relative protective effectiveness of a wood preservative in ground contact. Brussels, Belgium: European Committee for Standardisation. Chafe, S. (1987). Collapse, volumetric shrinkage, specific gravity and extractives in Eucalyptus and other species. Wood Science and Technology, 21(1), 27- 41. Chafe, S., Barnacle J.E, Hunter A.J, Llic R.L, Northway, R., & Roza A.N. (1992). Collapse: an introduction. Melbourne, Australia: CSIRO Division of Forest Products. Chauhan, S.S., & Walker, J.C. (2006). Variations in acoustic velocity and density with age, and their interrelationships in radiata pine. Forest Ecology and Management, 229(1-3), 388-394. https://doi. org/10.1016/j.foreco.2006.04.019 da Silva, A.R., Pastore, T.C.M., Braga, J.W.B., Davrieux, F., Okino, E.Y.A., Coradin, V.T.R., Camargos, J.A.A., & Do Prado, A.G.S. (2013). Assessment of total phenols and extractives of mahogany wood by Near Infrared Spectroscopy (NIRS). Holzforschung, 67(1), 1-8. https://doi.org/10.1515/hf-2011-0207 Dungey, H., Low, C., Lee, J., Miller, M., Fleet, K., & Yanchuk, A. (2012). Developing breeding and deployment options for Douglas-fir in New Zealand: breeding for future forest conditions. Silvae Genetica, 61(1- 6), 104-115. https://doi.org/10.1515/sg-2012- 0013 Eldridge, K., Davidson, J., Harwood, C., & Wyk, G.v. (1994). Eucalypt domestication and breeding. Oxford, UK: Clarendon Press. Elliott, C., & Byrne, M. (2003). Genetic diversity within and between natural populations of Eucalyptus occidentalis (Myrtaceae). Silvae Genetica, 52(3-4), 169-173. EN-350-1. (1994). Durability of wood and wood-based products - Natural durability of solid wood. Part 1: Guide to the principles of testing and classification of the natural durability of wood. Brussels, Belgium: European Committee for Standardization. Estopa, R.A., Milagres, F.R., Oliveira, R.A., & Hein, P.R.G. (2017). NIR spectroscopic models for phenotyping wood traits in breeding programs of Eucalyptus benthamii. Cerne, 23(3), 367-375. https://doi.org/ 10.1590/01047760201723032319 Fries, A., & Ericsson, T. (1998). Genetic parameters in diallel-crossed Scots pine favor heartwood formation breeding objectives. Canadian Journal of Forest Research, 28(6), 937-941. https://doi. org/10.1139/x98-061 Greaves, B.L., Schimleck, L., Borralho, N.M., Michell, A.J., & Raymond, C.A. (1996). Genetic control and repeatability of near infrared reflectance from Eucalyptus nitens woodmeal. Appita Journal, 49(6), 423-426. Guo, F., & Altaner, C. (2018). Properties of rotary peeled veneer and laminated veneer lumber (LVL) from New Zealand grown Eucalyptus globoidea. New Zealand Journal of Forestry Science, 48:3. https:// doi.org/10.1186/s40490-018-0109-7 Hamilton, M., Potts, B., Harwood, C., Apiolaza, L., & Gore, P. (2004). Comparison of non-destructive Iyiola et al. New Zealand Journal of Forestry Science (2022) 52:13 Page 9 https://doi.org/10.1051/forest:2004030 https://doi.org/10.3390/f5040744 https://doi.org/10.1007/s13595-011-0121-z https://doi.org/10.1007/s13595-011-0121-z https://doi.org/10.1016/j.foreco.2006.04.019 https://doi.org/10.1016/j.foreco.2006.04.019 https://doi.org/10.1515/hf-2011-0207 https://doi.org/10.1515/sg-2012-0013 https://doi.org/10.1515/sg-2012-0013 https://doi.org/10.1590/01047760201723032319 https://doi.org/10.1590/01047760201723032319 https://doi.org/10.1139/x98-061 https://doi.org/10.1139/x98-061 https://doi.org/10.1186/s40490-018-0109-7 https://doi.org/10.1186/s40490-018-0109-7 Iyiola et al. New Zealand Journal of Forestry Science (2022) 52:13 Page 10 assessment techniques for shrinkage and collapse in Eucalyptus nitens. Paper presented at “Eucalyptus in a Changing World”. Proceedings of IUFRO Conference, 2 p. Harju, A.M., & Venäläinen, M. (2006). Measuring the decay resistance of Scots pine heartwood indirectly by the Folin-Ciocalteu assay. Canadian Journal of Forest Research, 36(7), 1797-1804. https://doi. org/10.1139/x06-074 Hawley, L., Fleck, L., & Richards, C.A. (1924). The relation between durability and chemical composition in wood. Industrial & Engineering Chemistry, 16(7), 699-700. https://doi.org/10.1021/ie50175a015 Hong, Z., Fries, A., Lundqvist, S.-O., Andersson Gull, B., & Wu, H.X. (2015). Measuring stiffness using acoustic tool for Scots pine breeding selection. Scandinavian Journal of Forest Research, 30(4), 363-372. https:// doi.org/10.1080/02827581.2014.1001783 Joe, W., Dickson, R., Raymond, C., Ilic, J., & Matheson, A. (2004). Prediction of Eucalyptus dunnii and Pinus radiata timber stiffness using acoustics. A report for the RIRDC/Land & Water Australia/FWPRDC/ MDBS Joint Venture Agroforestry Program. Australia. https://www.agrifutures.com.au/wp- content/uploads/publications/04-013.pdf Jones, P., Schimleck, L.R., Daniels, R.F., Clark, A., & Purnell, R.C. (2008). Comparison of Pinus taeda L. whole-tree wood property calibrations using diffuse reflectance near infrared spectra obtained using a variety of sampling options. Wood Science and Technology, 42(5), 385-400. https://doi. org/10.1007/s00226-008-0178-5 Jones, T., McConnochie, R.M., Shelbourne, T., & Low, C.B. (2010). Sawing and grade recovery of 25-year-old Eucalyptus fastigata, E. globoidea, E. muelleriana and E. pilularis. New Zealand Journal of Forestry Science, 40, 19-31. King, J.N., Burdon, R., & Young, G. (1993). Provenance variation in New Zealand-grown Eucalyptus delegatensis. 2: Internal checking and other wood properties. New Zealand Journal of Forestry Science, 23(3), 314-323. Klápště, J., Stovold, T., & Dungey, H. (2019). Evaluation of wood stiffness in Douglas-fir progeny test FR280_2 and FR280_3. Retrieved 25 April 2022 from https://www.nzffa.org.nz/system/assets/5213/ SWP-T075.pdf Knapic, S., Tavares, F., & Pereira, H. (2006). Heartwood and sapwood variation in Acacia melanoxylon R. Br. trees in Portugal. Forestry, 79(4), 371-380. https:// doi.org/10.1093/forestry/cpl010 Kube, P.D. (2005). Genetic improvement of the wood properties of Eucalyptus nitens: Breeding to improve solid wood and pulp properties. Doctor of Philosophy. (PhD), University of Tasmania, Australia. Kumar, A., & Dhillon, G. (2014). Variation of sapwood and heartwood content in half-sib progenies of Eucalyptus tereticornis Sm. Indian Journal of Natural Products and Resources, 5(4), 338-344. Kuwabara, S. (2017). Prediction of extractive content of eucalyptus globoidea heartwood using near infrared spectroscopy. BSc. Hons, University of Canterbury, New Zealand. Lachenbruch, B., Moore, J.R., & Evans, R. (2011). Radial variation in wood structure and function in woody plants, and hypotheses for its occurrence. In: Meinzer F.C., Lachenbruch, B. & Dawson, T.E. (Eds.), Size-and age-related changes in tree structure and function (pp. 121-164): Springer. https://doi. org/10.1007/978-94-007-1242-3_5 Li, Y. (2018). Use of near infrared spectroscopy to predict wood traits in Eucalyptus species. PhD, University of Canterbury, New Zealand. Li, Y., & Altaner, C. (2019). Calibration of near infrared spectroscopy (NIRS) data of three Eucalyptus species with extractive contents determined by ASE extraction for rapid identification of species and high extractive contents. Holzforschung, 73(6), 537-545. https://doi.org/10.1515/hf-2018-0166 Li, Y., Apiolaza, L.A., & Altaner, C. (2018). Genetic variation in heartwood properties and growth traits of Eucalyptus bosistoana. European Journal of Forest Research, 137(4), 565-572. https://doi. org/10.1007/s10342-018-1125-0 Li, Y., Sharma, M., Altaner, C., & Cookson, L.J. (2020). An approach to quantify natural durability of Eucalyptus bosistoana by near infrared spectroscopy for genetic selection. Industrial Crops and Products, 154, 112676. https://doi. org/10.1016/j.indcrop.2020.112676 Lindström, H., Harris, P., & Nakada, R. (2002). Methods for measuring stiffness of young trees. Holz als Roh-und Werkstoff, 60(3), 165-174. https://doi. org/10.1007/s00107-002-0292-2 Love, S. (2007). Extended producer responsibility of treated timber waste. Paper presented to Sustainable Building Conference (SB07), Auckland, 14-16 November. Millen, P., Altaner, C., & Palmer, H. (2018). Naturally durable timber posts performing well. New Zealand Tree Grower 39(1), 24-26. Millen, P., van Ballekom, S., Altaner, C., Apiolaza, L., Mason, E., McConnochie, R., Morgenroth, J., & Murray, T. (2018). Durable eucalypt forests-a multi-regional opportunity for investment in New Zealand drylands. New Zealand Journal of Forestry, 63, 11-23. Miranda, I., Gominho, J., Araújo, C., & Pereira, H. (2014). Family effects in heartwood content of Eucalyptus globulus L. European Journal of Forest Research, https://doi.org/10.1139/x06-074 https://doi.org/10.1139/x06-074 https://doi.org/10.1021/ie50175a015 https://doi.org/10.1080/02827581.2014.1001783 https://doi.org/10.1080/02827581.2014.1001783 https://www.agrifutures.com.au/wp-content/uploads/publications/04-013.pdf https://www.agrifutures.com.au/wp-content/uploads/publications/04-013.pdf https://doi.org/10.1007/s00226-008-0178-5 https://doi.org/10.1007/s00226-008-0178-5 https://www.nzffa.org.nz/system/assets/5213/SWP-T075.pdf https://www.nzffa.org.nz/system/assets/5213/SWP-T075.pdf https://doi.org/10.1093/forestry/cpl010 https://doi.org/10.1093/forestry/cpl010 https://doi.org/10.1007/978-94-007-1242-3_5 https://doi.org/10.1007/978-94-007-1242-3_5 https://doi.org/10.1515/hf-2018-0166 https://doi.org/10.1007/s10342-018-1125-0 https://doi.org/10.1007/s10342-018-1125-0 https://doi.org/10.1016/j.indcrop.2020.112676 https://doi.org/10.1016/j.indcrop.2020.112676 https://doi.org/10.1007/s00107-002-0292-2 https://doi.org/10.1007/s00107-002-0292-2 Iyiola et al. New Zealand Journal of Forestry Science (2022) 52:13 Page 11 133(1), 81-87. https://doi.org/10.1007/s10342- 013-0741-y Pâques, L. (2001). Genetic control of heartwood content in larch. Silvae Genetica, 50(2), 69-74. Pâques, L., & Charpentier, J. (2015). Perspectives for genetic improvement in heartwood size and extractive content in relation to natural durability and aesthetics in interspecific hybrid larch (Larix × eurolepis). European Journal of Forest Research, 134(5), 857-868. https://doi.org/10.1007/ s10342-015-0895-x Pâques, L., García-Casas, M.d.C., & Charpentier, J. (2013). Distribution of heartwood extractives in hybrid larches and in their related European and Japanese larch parents: relationship with wood colour parameters. European Journal of Forest Research, 132(1), 61-69. https://doi.org/10.1007/s10342- 012-0654-1 Poke, F.S., Potts, B.M., Vaillancourt, R.E., & Raymond, C.A. (2006). Genetic parameters for lignin, extractives and decay in Eucalyptus globulus. Annals of Forest Science, 63(8), 813-821. https://doi.org/10.1051/ forest:2006080 R Core Team. (2019). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R- project.org/ Raymond, C.A. (2002). Genetics of Eucalyptus wood properties. Annals of Forest Science, 59(5-6), 525- 531. https://doi.org/10.1051/forest:2002037 Santos, P., Geraldi, I.O., & Garcia, J.N. (2004). Estimates of genetic parameters of wood traits for sawn timber production in Eucalyptus grandis. Genetics and Molecular Biology, 27(4), 567-573. https://doi. org/10.1590/S1415-47572004000400017 Stackpole, D.J., Vaillancourt, R.E., Alves, A., Rodrigues, J., & Potts, B.M. (2011). Genetic variation in the chemical components of Eucalyptus globulus wood. G3: Genes, Genomes, Genetics, 1(2), 151-159. https://doi.org/10.1534/g3.111.000372 Suontama, M., Miller, M., Fleet, K., Stovold, T., & Dungey, H. (2018). Genetic analysis of Eucalyptus fastigata progeny trials and implications to selection. Report for Speciality Wood Products partnership. 10. Retrieved 25 April 2022 from https://www.nzffa. org.nz/system/assets/5189/SWP-T050.pdf Toulmin, M.J., & Raymond, C. (2007). Developing a sampling strategy for measuring acoustic velocity in standing Pinus radiata using the TreeTap time of flight tool. New Zealand Journal of Forestry Science, 37(1), 96-111. Townsend, T.G., & Solo-Gabriele, H. (2006). Environmental impacts of treated wood: CRC Press. https://doi. org/10.1201/9781420006216 Tsuchikawa, S., & Schwanninger, M. (2013). A review of recent near-infrared research for wood and paper (Part 2). Applied Spectroscopy Reviews, 48(7), 560- 587. https://doi.org/10.1080/05704928.2011.62 1079 UNEP. (2012). Green Carbon, Black Trade: Illegal logging, tax fraud and laundering in the world’s tropical forests. United Nations Environment Programme. https://wedocs.unep.org/20.500.11822/8030 White, T.L., Adams, W.T., & Neale, D.B. (2007). Forest genetics. CABI. https://doi. org/10.1079/9781845932855.0000 Whiteman, P. (1992). Genetic parameters and selection strategies for Eucalyptus nitens (Dean and Maiden) in Victoria. Silvae Genetica, 41(2), 77-81. Wilkes, J. (1991). Heartwood development and its relationship to growth in Pinus radiata. Wood Science and Technology, 25(2), 85-90. https://doi. org/10.1007/BF00226808 Woeste, K.E. (2002). Heartwood production in a 35-year- old black walnut progeny test. Canadian Journal of Forest Research, 32(1), 177-181. https://doi. org/10.1139/x01-177 Wu, S., Zhu, Y., Xu, J., Lu, Z., Chen, G., Song, P., & Guo, W. (2017). Genetic variation and genetic gain for energy production, growth traits and wood properties in Eucalyptus hybrid clones in China. Australian Forestry, 80(1), 57-65. https://doi.org/ 10.1080/00049158.2016.1275948 Zobel, B., & Jett, J.B. (1995). Genetics of wood production. Springer Science & Business Media. https://doi. org/10.1007/978-3-642-79514-5 https://doi.org/10.1007/s10342-013-0741-y https://doi.org/10.1007/s10342-013-0741-y https://doi.org/10.1007/s10342-015-0895-x https://doi.org/10.1007/s10342-015-0895-x https://doi.org/10.1007/s10342-012-0654-1 https://doi.org/10.1007/s10342-012-0654-1 https://doi.org/10.1051/forest:2006080 https://doi.org/10.1051/forest:2006080 https://www.R-project.org/ https://www.R-project.org/ https://doi.org/10.1051/forest:2002037 https://doi.org/10.1590/S1415-47572004000400017 https://doi.org/10.1590/S1415-47572004000400017 https://doi.org/10.1534/g3.111.000372 https://www.nzffa.org.nz/system/assets/5189/SWP-T050.pdf https://www.nzffa.org.nz/system/assets/5189/SWP-T050.pdf https://doi.org/10.1201/9781420006216 https://doi.org/10.1201/9781420006216 https://doi.org/10.1080/05704928.2011.621079 https://doi.org/10.1080/05704928.2011.621079 https://wedocs.unep.org/20.500.11822/8030 https://doi.org/10.1079/9781845932855.0000 https://doi.org/10.1079/9781845932855.0000 https://doi.org/10.1007/BF00226808 https://doi.org/10.1007/BF00226808 https://doi.org/10.1139/x01-177 https://doi.org/10.1139/x01-177 https://doi.org/10.1080/00049158.2016.1275948 https://doi.org/10.1080/00049158.2016.1275948 https://doi.org/10.1007/978-3-642-79514-5 https://doi.org/10.1007/978-3-642-79514-5