Abstract
Monetary valuation using urban tree appraisals can be performed with formulas, a common practice in many countries. This study compares twelve parametric type formulas: Amenity Valuation of Tree and Woodlands (Helliwell), Standard Tree Evaluation Method (STEM), French Method, Italian Method, Tedesco, Norma Granada, Trunk Replacement Formula (CTLA), Burnley Method, Danish Method, Swiss Method, and two Chilean formulas used in Municipalities of Concepción, La Pintana, and Maipú (COPIMA Method), and Peñalolén Method. Formulas were then applied to 30 trees located in Santiago, Talca and Concepción, Chile.
Researchers used eight appraisers divided into two groups, according to senior-level and junior-level experience. Statistical differences were determined using the Kruskal-Wallis test of non-parametric variance, while Fisher’s least significant difference test was used to identify homogeneous groups. The results show a wide dispersion of values that were high for “emblematic” trees and low for young or low-vigor trees.
Formula, type of appraisers, and inter-appraiser differences formed nine, two, and three groups, respectively. The lowest-appraised trees were obtained using the Danish and French Method, while the highest values were obtained with the Burnley, Helliwell, and STEM formulas. Although there were differences in tree value according to the type of appraiser, when comparing difference among appraisers, researchers found these were not due to experience level. Given the wide range of values found, the study authors cannot recommend any specific formula(s) for assessing urban trees, as results will depend on the variables of interest used in the formulas and their intended application and use.
Urban trees are defined as trees found in areas located in urban or peri-urban areas (Tyrväinen et al. 2003), in residential and commercial sidewalks, parks, greenbelts, industrial parks, and vacant lots, among other land uses (Cordell et al. 1984). Different authors acknowledge their contribution to the sustainability of cities and their role as an economic asset (Tyrväinen et al. 2003; Konijnendijk et al. 2004; Ponce-Donoso et al. 2009). While recognizing the difficulty in determining the monetary value of urban trees, there are several available methods for calculating this value (Caballer 1999; Tyrväinen 2001; Watson 2002; Price 2003; Grande-Ortiz et al. 2012). The most common appraisal method for monetary valuation of urban trees is using formulas (Watson 2001), which are commonly of two types: parametric and capitalization. Parametric, also known as multiplication, is defined as the quantification of one or more variables, including both structural and other subjective ones (e.g., aesthetics, botanical, location, or significance), which, as recommended, should be performed by experienced appraisers, as they are not always accurate (Price 2003). Capitalization formulas use more conventional econometric methods (Grande-Ortiz et al. 2012). Chueca (2001) has indicated that these formulations require that the subjectivity of the variables used be reduced in the applied formulas. Hence, some professionals recommend the use of capitalization type formulas for their simplicity, while others prefer the parametric formula, stating that it better reflects the true total economic value of the tree (Petersen and Straka 2011; Grande-Ortiz et al. 2012).
Currently, several studies are assessing the different ecosystem services provided by urban trees, such as property value increases, carbon sequestration, reduction of noise and pollutants, energy use savings, and others for different conditions and contexts (Dobbs et al. 2011; Roy et al. 2012; Haase et al. 2014; Escobedo et al. 2015). Likewise, available formulas that allow for the monetary valuation of these urban forest ecosystem services are commonly used in methods and models such as i-Tree (i-Tree 2012) and CAVAT (Neilan 2010).
However, these valuation formulas are based on different assumptions and approaches that vary by form, application, and in the final total monetary value of the tree (Randrup 2005). For example, some of these formulas are based on the cost of replacing a damaged or vandalized tree (i.e., replacement cost), and are adjusted according to commonly used factors, such as tree vitality, damage type, location, aesthetics, overall amenities, age, and even provision of environmental services, which finally deliver a comparable monetary value (Moore and Arthur 1992; CTLA 2000; Helliwell 2008; Ponce-Donoso et al. 2013; Östberg and Sjögren 2016). Some of these formulas include the Trunk Formula Method from the Council of Tree and Landscape Appraisers of the United States (CTLA 2000), Burnley Method of Australia (Moore and Arthur 1992), Amenity of Trees and Woodland of the United Kingdom or Helliwell Method (Helliwell 2008), Norma Granada of Spanish Association of Parks and Public Gardens (AEPJP 2007), and the Standard Tree Evaluation Method (STEM) of New Zealand (Flook 1996).
These various formulas will indeed have different levels of acceptance and validity even within their country of origin. The Helliwell Method, for example, is accepted and regularly used in the UK, as is the CTLA Method in the U.S. and Canada (Cullen 2005; Cullen 2007). But even within these countries their acceptance is not complete among all legal contexts and by individual tree appraisers because of the variability in the range of estimated values (Watson 2001; Watson 2002). Furthermore, since these formulas were developed for use in primarily temperate, industrialized, English-speaking countries, the factors used might not be relevant to different cultural, ecological, and socio-political contexts (i.e., tree forms and species from tropical environments, land-use definitions from emerging countries, culturally specific tree maintenance practices, translation of English-language variables).
Accordingly, several international studies have applied parametric and capitalization valuation formulas in different contexts, such as Argentina (Contato-Carol et al. 2008), Brazil (Leal et al. 2008), Chile (Ponce-Donoso et al. 2012; Ponce-Donoso et al. 2013), Spain (Grande-Ortiz et al. 2012), the United States (Watson 2002), Finland (Tyrväinen 2001), Hungary (Hegedüs et al. 2011), as well as for applied extension education programs in these English-speaking countries (Harris 2007; Sarajevs 2011). As such, these different formulas can be used for different applications internationally, including damage assessments, legal claims or investment values, and replacement and damage costs, among others (Grande-Ortiz et al. 2012; Östberg and Sjögren 2016).
However, the appraised values will be influenced by how, and the purposes for which the formulas are used. The STEM formula is commonly used in its existing form, while some formulas are derived from others—such as the CAVAT, which was derived from the Helliwell Method (Randrup 2005), and the French Method, which originates from the Swiss Method via an adaptation that accounts for maintenance of ornamental species (Contato-Carol et al. 2008). Similarly, the CTLA Method has inspired the Danish Method (Randrup 2005). In South America, different appraisal methods are applied inconsistently, and as in other regions, their acceptance by local judges has been varied (Contato-Carol et al. 2008; Ponce-Donoso et al. 2009). As a result, there is a need to quantitatively evaluate these different valuation formulas and their parameters or variables under different contexts for their use in urban forest management and legally related activities.
Moreover, studies highlight the subjectivity of the appraisers when applying these formulas. The subjectivity of which results in high variability in values, in particular with the Helliwell Method, whereas with CTLA and Burnley Methods, values have been reported to be lower (Watson 2002). Ponce-Donoso et al. (2012) found high values for the STEM formula, average values for CTLA, and low values with the Burnley Method. Contato-Carol et al. (2008) noted that the Swiss and Finnish Methods had higher values, while similar mid-range values were obtained with both the CTLA and French Methods, which are due to the inclusion of aesthetics and ornamental variables, as previously mentioned; although, in larger specimens, the CTLA Method presented higher values than did the French Method. The commonly used CTLA Method’s validity has been questioned because of its subjectivity, which leads to statistically significant differences (Cullen 2007) and mid-range monetary values when compared with other formulas (Contato-Carol et al. 2008; Ponce-Donoso et al. 2012). Further, the CTLA Method showed lower values when compared to other methods, indicating the need to test the different methods outside their country of origin because extreme differences are not always obvious (Watson 2002). Grande-Ortiz et al. (2012) does, however, indicate that the CTLA Method can be widely used because of its low degree of difficulty, thus providing a comparatively good and available method for international use.
Therefore, the objective of this study was to analyze the monetary values obtained by the application of these twelve urban tree assessment formulas. The study was done using eight different appraisers with varied backgrounds, separated into two groups. Researchers also used three different cities in central Chile to develop three scenarios. This quantitative approach can be used to develop a single formula for central Chile, since currently eight different formulas are being used in ten municipal courts in this region, resulting in broad and disparate results (Ponce-Donoso et al. 2012).
MATERIALS AND METHODS
The study analyzed public urban trees in three Chilean cities. The first city was the Municipality of Santiago, located in the Metropolitan Region, with an elevation of 599 m above sea level, with an area of 22.4 km2, and 200,800 inhabitants. The second city was Talca, located in the Maule Region at 102 m above sea level, with an area of 232 km2 and 201,800 inhabitants. Finally, the third city considered in the study was Concepción, located in the Bio-Bio Region, 12 m above sea level with an area of 221.6 km2 and 216,100 inhabitants (NCL 2013).
Following conventional formula methods, 30 trees were selected, representing a total of 16 tree species (Table 1). The design was based on eight different appraisers, consisting of foresters and agronomists, who were separated into two groups. One group was made up of professional experts with at least five years’ experience (Senior Group, SG), while the other group consisted of professionals with no experience in tree valuation (Junior Group, JG). All participants received the same instructions on the use of tree appraisal formulas, resulting in a total of 1,440 valuation appraisals.
The field work was conducted during the southern latitude summer months of December 2013 and February 2014, when the trees exhibited the best conditions for appraisal. Both biometric variables and those related to the aesthetic, condition, and location were measured and appraised. Selling prices were collected in local wholesale and retail nurseries, as well as annual maintenance costs reported by the Municipality of Talca, and supplemented with information from the Municipalities of Santiago and Concepción. The maintenance cost was calculated based on the annualized costs, including pruning, pest control, watering, and others. The price of the species in the nurseries was based on the average prices in both retail and wholesale markets.
The formulas analyzed in this study were selected with consideration to the best performance for valuation of a tree within a public-use area in a municipality, as well as by its speed of implementation and calculation, efficiency and effectiveness in data collection, and overall simplicity to appraisers. The formulas are described as follows:
The Municipalities of Concepción, La Pintana, and Maipú of Chile (refereed to hereafter as the COPIMA Method; Ponce-Donoso et al. 2009) use the following formula:
1
where A = price of species at the local market, B = aesthetic and condition value of the tree, C = situation index, and D = dimension index.
Municipality of Peñalolén of Chile Method (Ponce-Donoso et al. 2009) is as follows:
2
where A = location factor, B = tree condition as a percentage of the damage present, e = age of the species, and VA = tree value according to species and age.
Amenity Valuation of Tree and Woodlands, referred to as the Helliwell Method (Helliwell 2008), estimates the visual amenities based on a point range from 1.0 to 4.0, which accounts for seven factors:
3
The Standard Tree Evaluation Method, known as STEM in New Zealand (Flook 1996), uses a point system based on 20 attributes (3 to 27 points each), characterizing a tree’s condition, amenities, and special features of notability:
4
The French Method (Ferraris 1984) corresponds to a method that provides an index related to the maintenance and care of the tree. It is based on Swiss Method but includes an additional factor to set a monetary value in parks and private gardens:
5
where E = species and variety index, based on the reference price in the nursery, B = health and aesthetic index, U = location index, and D = dimension index.
The Italian Method (Fabbri 1989) is as follows:
6
where P = price of the same species in local nurseries; I = reflects the health and appearance of the tree; S = location index, rural or urban; and C = size index.
The Tedesco Method, from Italy (Bernatzky 1978), is as follows:
7
where Vb = basic value 1/10 of market price for tree 10 cm2 of basal area; ID = dimension index in function of DAP or circumference; IP = position index; IC = condition index, including spacing between trees, tree development, condition, and damage; IIA = environmental compatibility index, which considers variables, such as insertion into the landscape, compatibility with the soil type, and execution of the planting; IE = age index, which is related to the age of the tree that exceeds the average age of the species; and IR = reduction index due to stem damage.
The Granada Norm of the Spanish Association of Public Parks and Gardens (Norma Granada; AEPJP 2007) is a formula for a non-replaceable tree, and corresponds to:
8
where Vb = basic value of the tree, which is determined by the function ῳ * µ * (0.0059 * p2 + 0.0601 * p – 0.324). ῳ = updated coefficient corresponding to the species, fixed for each climate zone according to Köppen; µ = soil corrector coefficient; p = perimeter of the trunk 1 m above the ground. Els = intrinsic factors of the tree (roots, trunk, main structural branches, sub-branches and terminal, leaves), and Ele = tree extrinsic factors (aesthetic and functional, representativeness and rarity, situation).
The Trunk Replacement Formula from the Council of Tree and Landscape Appraisal, the CTLA Method of the United States (CTLA 2000), considers the area of the cross section of the trunk, 1.4 m over the ground level, multiplied by a value based on the cost of the regional species available in local nurseries. The value is then multiplied by corrector indices (species, condition, and location) to reduce or maintain this value:
9
where the trunk area is expressed in cm2 and the basic price is expressed per unit of cm2. The species factor relates attributes of the tree associated with tree growth, life expectancy, adaptability to environmental conditions, maintenance requirements, and other amenities. The condition is related to the characteristics of the health and vigor of the tree. The location factor corresponds to the location of the tree in the city.
The Burney Method of Australia (Moore and Arthur 1992) is as follows:
10
where a number of points related to the volume of the tree are assigned, which correspond to an inverted cone, including the base value, which is the cost per cubic meter in retail nurseries, and other shape factors, vigor, and location.
The Danish Method (Randrup 2005) is as follows:
11
where B = base, which is expressed as E + (Pn / Cn) * (Cd / Cn), where E = costs of establishing value, Pn = price of a new tree, Cn = circumference of a new tree, and Cd = circumference of the evaluated tree; H = the health index, which is expressed as the condition of (r + t + rp + rs + f) / 25, being roots (r), trunk (t), main branches (rp), secondary branches and twigs (rs), leaves and buds (f); L = index location, which is expressed as (n + a + ve + v + F) / 25, where natural ecological adaptation is (n), architecture (a), aesthetic value (ve), visibility (v) and environmental factors (F); and A = age index, which is expressed as [((b – a) * 2) / b] – 2, where a = current age and b = life expectancy.
The Swiss Method (Ferraris 1984) is as follows:
12
where Pb = base price; ID = dimension index in function of the circumference trunk; IP = location index, which varies from the center of the city to a rural area; IER = aesthetic index and sanitary condition, which is related to vegetative vigor; and IR = reduction due to damage index, which is applied as a percentage of the trunk.
For the statistical analysis, researchers used both mean and median values to better reduce the effects of outliers. The following hypotheses were used to account for sources of variation such as the specific valuation formulas used, type of appraisers used, and inter- and intra-appraiser comparisons:
Ho: αi = αj / i ≠ j; (i.e., there are no statistically significant differences between the medians of the variation sources).
Ho: αi ≠ αj / i ≠ j; (i.e., there are statistically significant differences between the medians of the variation sources).
Analyses of variance (ANOVAs) were used to determine if there were statistically significant differences between the sources of variation. The assumptions of homoscedasticity and normality were not met in all cases because of a high coefficient of variation (209.26%). Similarly, despite the transformation of data, the bias and standardized kurtosis were high (143.20 and 704.48, respectively), exceeding Kirk’s (1995) preset limit value of 2.0. In these cases, the non-parametric ANOVA Kruskal-Wallis (Conover 1999) was used, as it is less sensitive to the presence of atypical values.
Both the data obtained from the SG and JG experience groups were ranked, according to their position in the ascending order of the data; with 1 having the lowest valuation and ranking, and 2,880 the highest, while the intermediate rankings corresponded to intermediate values for each group. Statistically significant differences were found between sources of variation, so the least significant difference test (DMS; P ≤ 0.001) was applied (Conover 1999). Also, an analysis of nonparametric variance was separately conducted for each of the formulas to better observe the variability among them. Microsoft® Excel® Version 2003 and Statgraphics Centurion for Windows® (StatPoint Technologies, Inc., Warrenton, Virginia, U.S.) were used for all statistical analyses.
RESULTS
Table 2 shows the results of the average values for each of the analyzed formulas. According to the non-parametric Kruskal-Wallis analysis, the results show statistically significant differences between the median ranking in all sources of variation. In the formulas, the value was 1,235.23 (P ≤ 0.00); 4.30 was obtained between type of appraisers (P ≤ 0.04), and 14.97 between all appraisers (P ≤ 0.04).
Researchers also identified distinct homogeneous groups of formulas and appraisers as determined by the LSD test (Table 3). The formula conform nine groups overall, indicating heterogeneous and dispersed values according to the ranking, while there were two different groups for appraiser. When analyzing according to appraiser, three distinct groups were found, confirming that experience is indeed a source of differentiation (Table 3).
The median and the degree of dispersion of the ranking system for all formulas, experience of appraisers, and according to individual appraiser, are shown in Figure 1. Trends are evident, as the formulas are ranked according to the monetary value obtained. The Danish, French, and CTLA Methods are located in the lower area, while the Helliwell, Burnley and STEM Methods in the upper area.
The results of the Kruskal-Wallis test (Table 4) show results where each formula was analyzed separately, thus allowing for an evaluation of the formula by different appraisers. The Danish, Tedesco, Burney, and Helliwell Methods did show statistically significant differences between appraisers (P < 0.05).
In general, results also show very different values for each tree using the same formula and different appraisers. For example, Table 2 shows that the extreme average values that correspond to two trees in the Municipality of Talca. Tree #13, a Ginkgo biloba, scored the maximum value using the Helliwell Method (USD $40,347), locating it at the top of the ranking and corresponding to a locally emblematic specimen of more than 100-years-old, while the minimum was awarded to Tree #19, a Betula pendula, which was characterized as an old tree with reduced–intermediate vigor (USD $27) using the CTLA formula, and is located at the bottom of the value estimates. However, most importantly, researchers found disparate appraisal values according to individual appraisers and by their experience groups (Figure 1). The average value for each tree fluctuated between USD $13,823 and USD $636, with the maximum value obtained by Tree #13 and with a minimum rating for Tree #30, a Catalpa bignonioides located in Concepción. These same trends in valuation estimates have been reported by Watson (2002), Contato-Carol et al. (2008), Ponce-Donoso et al. (2009; 2012; 2013), and Grande-Ortiz et al. (2012).
In general, researchers found that the first quartile the highest average values were presented by the Helliwell, Burnley, and STEM formulas; with a value of USD $7,342 for the Helliwell Method. Second were Norma Granada (of USD $1,760), Tedesco, and COPIMA. In the third quartile the formulas were the Tedesco, the Italian, Swiss, and Peñalolén formulas, with an average of USD $574 for first one. In the bottom quartile, with the lowest values, were the CTLA, French, and Danish (USD $161) formulas.
DISCUSSION
The results show that in using average monetary values, the Danish Method delivered low values, with an average of USD $207, and a median USD $161. The Swiss Method, as reported by Contato-Carol et al. (2008), resulted in the fourth lowest average value (USD $548). The Burnley had the fourth highest average value (USD $2,688) in contrast to what was reported by Ponce-Donoso et al. (2012). However, findings were consistent with these same authors regarding the CTLA and French Methods (Contato-Carol et al. 2008; Ponce-Donoso et al. 2013). In regards to the two Chilean formulas analyzed (i.e., Peñalolén Method and COPIMA), the findings were similar to those reported by Ponce-Donoso et al. (2009; 2012; 2013). These authors found that intermediate values were obtained, and the COPIMA continued to perform better with regards to the Chilean formulas.
A wide dispersion in the average rankings were also found (Table 3; Figure 1), with differences of more than fourfold between the lowest and highest value, as exemplified by the Danish and French Methods with regards to the STEM formula, due to the fact they present statistically significant differences from the rest of the formulas in Groups A, B, and I. The French Method was followed by the CTLA, Swiss, and Peñalolén formulas, whereas mean values were displayed by the Italian Method and COPIMA, forming a group with no statistical differences in their medians. The Tedesco and Norma Granada formulas displayed slightly higher intermediate values. Finally, those with higher average values were the Burnley, Helliwell, and STEM formulas; the first two forming a group with no statistical differences in their medians (Group H; Table 3). Overall, the results showed a wide variability in the studied formulas, whose differences in the valuation were due to their structural characteristics, and just as important, appraiser variability (Grande-Ortiz et al. 2012; Ponce-Donoso et al. 2013).
When comparing the results and overall formula performance of the current study to those reported by Watson (2002), researchers discovered similar findings exhibiting high variability. While the Helliwell Method displayed the highest values in Watson’s (2002) study, the current findings show values from the Danish Method were followed by Helliwell, Tedesco, and Burnley Methods; and further, they showed no statistical significant difference when analyzed separately according to appraiser (P < 0.05; Table 4). The remaining formulas did show differences between appraisers, with the greatest in the Peñalolén, COPIMA, French, Italian, and Swiss methods (P > 0.8); a finding pointed out by Watson (2002) when comparing these same formulas. Conversely, the CTLA, Italian, French, Norma Granada, COPIMA, and Peñalolén methods showed the least variability, indicating that differences between appraisers were minimal.
Considering that all formulas were parametric, their application in these three cites delivered both low and high values (e.g., CTLA and STEM Methods) that differed with other studies (Watson 2002; Contato-Carol et al. 2008). It is not clear, however, which factor or parameter, either multiplicative or additive, is most influential and has the greatest weight in the total appraised value. Future research could analyze the sensitivity of each individual parameter separately, so as to recognize the weight they have in formula performance. Further, research is also needed on other key variables that could be used in these appraisal formulas, such as demographic characteristics (e.g., average monetary income of the city’s population), as well feasibility and viability of applying the formula.
Researchers were not able to identify with any degree of certainty which of the formulas was best at assessing each tree, as the analyzed formulas performed dissimilarly. However, when analyzing the range and median of the ranking as performance criteria, the Norma Granada, CTLA, and COPIMA formulas performed well, relative to the others (Figure 1), by achieving low average values, which are relevant for the context in which the study took place. This was also corroborated by the lower degree of dispersion [i.e., a high probability (Table 4)] and a median showing a distribution of values indicating a fair appraisal. The formulas that had the highest values were the Helliwell and Burnley; findings that are similar to what were reported in other studies (Watson 2002; Ponce-Donoso et al. 2012; Ponce-Donoso et al. 2013). Researchers also note that STEM formulas exhibited high values, although with only a moderate dispersion.
In comparing the experience level of appraisers (i.e., SG and JG), significant statistical differences were found, though the difference was minimal if the average ranking is considered (Table 3; Figure 1). This indicates further training could reduce the probability of statistical differences. When comparing all appraisers, in one group, statistically significant differences were found, thereby resulting in three different groups. Nevertheless, rankings were not ordered in function of the group’s experience (Table 3; Figure 1), thus the individuals’ experience would not be a differentiating factor, statistically speaking. As shown in Table 3, Group L incorporated three appraisers with experience and one without, Group M included two experienced and two inexperienced appraisers, and Group N had two experienced and four inexperienced appraisers (Table 3). Thus, the lower average value was obtained by two senior appraisers (S1 and S2), and the highest valuation was also obtained by one of these, while all JG appraisers are located within the ranking (Table 3).
On the other hand, the role of experts in appraisal procedures, provided a distinctive element in the valuations, where the most experienced appraisers consistently tended to get lower values in their valuation (Cullen 2005), which was contrary to the norm that experience is a requirement for the appraisal of the tree (Price 2003; Tyrväinen et al. 2003). This shows the need for a minimum level of instruction to obtain the professional skills necessary for implementing a formula, considering that these appraisers had more of an ad hoc professional training. Therefore, when considering this particular case, the Danish Method, French Method, CTLA Method, Swiss Method, Peñalolén Method, Italian Method, and COPIMA (excluding Norma Granada, Helliwell, Burnley, Tedesco, and STEM) have no statistically significant differences between the type of appraisers or among the eight appraisers.
CONCLUSIONS
Overall, researchers found statistically significant differences between the medians of the studied formulas, among types of appraisers, and among individual appraisers. This indicates that the valuation of urban trees depends preferably on the specific formula used (i.e., type and its constituent variables) as well as the appraiser’s experience; aspects that coincide with studies by Watson’s (2002), Contato-Carol et al. (2008), Grande-Ortiz et al. (2012) and Ponce-Donoso et al. (2012; 2013). The formulas that presented the best comparative performance as defined by dispersion and location of the median, were the Helliwell, Norma Granada, Tedesco, Burnley, and CTLA Methods, while the lower performance were found in the French, Swiss, Danish, and STEM Methods.
Researcher note that when the appraisers were organized into three homogeneous groups, findings show that experience did not seem to be a factor that differed statistically. Group N was the largest group, composed of 75% of the appraisers, all JG, and only two SG (Table 3). Group L included four appraisers and only one having junior experience, while Group M was composed of two each (Table 3). Again, the experience would be a determining factor when it comes to valuation of urban trees (Grande-Ortiz et al. 2012).
The results show that the use of parametric formulas is recommended when appraising urban trees in international contexts. This was observed by the ability of formulas to discriminate the value of trees, resulting in the graphic dispersion of monetary values presented in the rankings (i.e., high dollar amount for higher value tree, medium monetary value for mean tree, and lowest for those with a lower value monetary values). In this study, these features were identified in the CTLA, Helliwell, and Norma Granada Methods.
Acknowledgments
The authors thanks to Fondo Nacional de Desarrollo Científico y Tecnológico of Chile (FONDECYT) Project 1130264 for funding this project.
- © 2017, International Society of Arboriculture. All rights reserved.