Performance of Eucalyptus tereticornis provenances in subtropical climate

Eucalyptus tereticornis has a wide natural geographic distribution from Papua New Guinea to South Australia (6o-38o S), therefore, it is expected the existence of important differences among provenances. Although differences in performance among provenances have been reported in several countries, the existing information is insufficient to allow the choice of the most appropriate provenance as a source of seeds for commercial plantations for humid subtropical climates. In order to evaluate the performance of 20 provenances of E. tereticornis in Uruguay and to generate information that contributes to the identification of the best provenances for humid subtropical climates, survival, individual tree volume, stem straightness and forking were assessed at 9-years-old in 4 sites in the center, north and northeast of Uruguay. Significant differences among provenances for all traits and significant provenance by site interaction for survival, individual tree volume and forking were found. No provenance was best in all sites. The Buckenbowra provenance had simultaneously high volume and survival. The absence of a clear regionalization, coupled with a significant provenance-by site interaction, suggests that local evaluation of different provenances is essential to identify appropriate seed sources, both for tree breeding and for use in commercial plantations.


Introduction
Eucalyptus tereticornis Sm. has the largest natural distribution range among eucalypts species, from Papua New Guinea to South Australia (6-38ºS), from sea level to 1000 m in Australia and 1800 m in Papua New Guinea (1) . It is found in coastal and continental areas, with variable climate, from tropical (North) to cold temperate (South), with an average maximum temperature of 24-36 ºC and minimums of 1-19 ºC, with up to 30 frosts per year. The annual rainfall ranges from 500 to 3000 mm, with a variable regime, from monsoonal in Papua New Guinea to summer rains in Queensland, an evenly distributed rains in New South Wales or winter rains in Victoria (2) .
At the end of the 19 th century, E. tereticornis was introduced in different countries (2) . Its cultivation has spread widely, mainly to areas with annual rainfall between 800 and 1500 mm. However, due to its great plasticity, it has been planted in areas of low rainfall (400 mm in India, 550 mm in Israel, and 580 mm in Zimbabwe) and in areas of high precipitation (2180 mm in Colombia and 3500 mm in Papua New Guinea) (3) . In India, one of the world´s largest growers of Eucalyptus, with an area of 3.94 million hectares of plantations, E. tereticornis and an interspecific hybrid of this species are the predominant species in the southern plains (4) (5) . The species has been widely used in other countries, as in Zaire, Pakistan (in agroforestry systems) and Colombia (in the Atlantic coast, mainly in reforestation programs with small landowners) (2)(3)(4) (6) . In Uruguay, the species was introduced in 1853 and due to its adaptation to a wide variety of soils, it was one of the most cultivated species until the middle of 20 th century, mainly for shade and shelter for livestock, but also for productive purposes (poles, wood for carpentry and floors) and energy (firewood and charcoal) (7) .
Given the great variability of environments in which this species occurs naturally, the existence of important differences among provenances is quite expected. Although differences in performance among provenances have been reported in several countries (8)(9)(10) (11) , the existing information is insufficient to allow the choice of the most appropriate provenance as a source of seeds for commercial plantations. For this reason, the aim of the present study was to evaluate the performance of different provenances of E. tereticornis in the subtropical climate of Uruguay to generate information that contributes to the identification of the best provenances for humid subtropical climates.

Genetic material and field trials
The data analyzed in this work comes from the breeding population of E. tereticornis of INIA (National Institute of Agricultural Research), Uruguay, which comprises 100 open-pollinated families, of 20 provenances that cover a large part of the natural distribution of this species, introduced from CSIRO (Commonwealth Scientific and Industrial Research Organization), Australia ( Table 1). The evaluation of productive performance was carried out in 4 provenance trials, established in representative sites of the country, in 2008 and 2009. The sites have a similar climate and vary in soil characteristics ( Table 2). The trials were established using a randomized complete block design, with single-tree plots, and 19 to 21 replications per site. Soil preparation consisted in tillage made in bands, and trees were planted at a spacing of 4×2.5 m (1000 trees/ha) in all sites.

Traits evaluated
Survival, diameter at 1.3 m (DBH), total tree height (H), stem form and incidence of forking were assessed at 9-year-old. All surviving trees were measured for DBH, with diametric tape, and H, with hypsometer (Vertex IV HS 102). The stem straightness and forking were both assessed visually, with grading scales of 1 (straight stem) to 3 (twisted stem), and 0 (without forks) to 2 (more than one fork), respectively. Individual tree volume with bark (VOL) was calculated for each tree, using the general equation: VOL= 3.1415/4(DBH 2 )*H*F, where F is the form factor, with a value of 0.4. Since there is no adjusted volume equation for this species, an F of 0.4 was assumed to correct taper.

Data analysis
Analysis of variance of each trait was carried out jointly for all sites, using the following mixed model: Yijk = μ + SITEi + REPj(i) + PROVk + SITE.PROVik + εijk where Yijk is the observed value in the j-th replication of the i-th site, pertaining to k-th provenance; μ is the overall mean value; SITEi is the effect of the i-th site; REPj(i) is the effect of the j-th replication within the i-th site; PROVk is the effect of the k-th provenance; SITE.PROVik is the interaction effect between the i-th site and the k-th provenance; εijk is the random error. Site and provenance effects were considered as fixed factors and the effect of replications as a random factor. The significance of the fixed effects was determined by analysis of variance (ANOVA). Due to the heterogeneity variances among treatments the contrasts of means were carried out using the Tukey HSD test with 95% confidence. The statistical analyses were carried out using the R3.4.1 software (12) . For the survival analysis, a binomial distribution was considered, assigning 0 for a dead tree and 1 for a living tree.
The study of the adaptability of the provenances in the different sites was analyzed according to Finlay and Wilkinson (13) . For each provenance present in at least three sites, a linear regression of individual volume on the mean of all provenances for each site was computed. The provenances mean and the regression lines illustrate different types of responses to the range of environments. The population mean has a regression coefficient of 1, so coefficients approximating to 1 indicate average stability. When this is associated with high mean yield, provenances have good general adaptability; when associated with low mean yield, provenances are poorly adapted to all environments. In contrast, regression coefficients above 1 describe provenances with greater adaptability to high-yield sites, and below 1 to provenances more adapted to low-yield sites.

Results
For all traits, a significant effect of site (P < 0.001) and provenance (P < 0.001) was found. The provenance by site interaction was significant for all traits, except for stem straightness (P = 0.1025), so survival, individual tree volume and forking were evaluated in each site separately.

Survival
Differences among provenances for survival were pronounced (P < 0.001) in all sites. In the Rivera trial, survival ranged between 52 to 79%. The provenance that had the greatest survival was Burdekin River, significantly overcoming the Maryland, S of Loch Sport, Yurammie SF, W of Warwick, Cardwell, and SW of Mount Garnet provenances (Figure 1a). In El Carmen, survival ranged between 51 to 92%, being Kalpower the provenance with higher survival, statistically overcoming Maryland, Crediton SF, Yurammie SF, Cardwell, SW of Mount Garnet, and W of Warwick provenances (Figure 1b). The higher differences between provenances were recorded in the Carlos Reyles trial, with survival values that varied between 19 to 95%. Bunya SF exhibited the greatest survival, although it was only significantly higher than Cardwell and W of Warwick provenances, which had the lowest survival (Figure 1c).
Large differences between provenances were also recorded in Vergara, where survival ranged between 28 to 70%. The provenance that had higher survival was Buckenbowra SF, differing only from the Selection Flat SF559, Mill Stream Archer Creek, Bunya SF, Kalpower, Marlborough, Calliope River, and Cardwell provenances (Figure 1d).

Individual volume
Significant differences (P < 0.001) among provenances for growth rate were found in all sites. In Rivera, the most productive site, the average VOL ranged from 0.34 to 0.11 m 3 . Provenance Selection Flat SF559 was the fastest growing, overcoming the rest of the provenances, except for Buckenbowra SF, Bunya SF, Yurammie SF, and Calliope River, with which it did not differ significantly (Figure 2a). The mean VOL in El Carmen ranged from 0.20 to 0.065 m 3 . Buckenbowra SF was the provenance with higher growth in this trial, significantly surpassing the other provenances, except for Yurammie SF and Bunya SF, with which there were no significant differences (Figure 2b). In Carlos Reyles, where VOL ranged from 0.23 to 0.04 m 3 , the Bunya SF provenance had the greatest growth, although it only significantly exceeded the provenances of MT Garnet, Raglan Creek, Marlborough, Burdekin River, Mill Stream Archer Creek, Cardwell, and W of Warwick (Figure 2c). Vergara was the site with larger differences in growth rate, with a VOL that ranged from 0.36 to 0.076 m 3 . Yurammie SF and Buckenbowra SF were the provenances with higher growth ( Figure  2d).
The regression analysis shows that the provenances of Buckenbowra SF, Raglan Creek, Burde-kin River, MT Garnet, Cardwell, Yurammie SF, Crediton SF and Kalpower have a regression coefficient close to 1, thus indicating a stable behavior across sites ( Figure 3 and Table 3). In contrast, some provenances have regression coefficients greater than 1 (Selection Flat SF559, Mill Stream Archer Creek and Calliope River) or less than 1 (SW of Mount Garnet, S of Loch Sport, Maryland and W of Warwick), being therefore responsible for the provenance by site interaction. Within the group of stable provenances, Buckenbowra SF and Yurammie SF have a mean volume greater than the whole provenances average, so they have good general adaptability, and Cardwell and Burdekin River have lower volume than the average, showing poor general adaptability (poor performance in all sites).

Forking
Forking scores showed marked provenance differences (P < 0.001) in all sites. The mean incidence of forking ranged from 0.4 to 1.3 in the Rivera trial. The N of Woolgoolga, W of Warwick, Yurammie SF, and Kalpower provenances showed the lowest number of forks, although they only significantly differed from the S of Loch Sport and Loch Sport provenances, which presented the highest incidence of forking ( Figure 4a). El Carmen was the site with the greatest differences among provenances, with an incidence of forking ranging from 0.25 to 1.35. The Raglan Creek, Marlborough, Cardwell and Burdekin River provenances showed the best performance in this trait (low incidence of forking), although they only overcome the MT Garnet, Bunya SF, Maryland, Buckenbowra SF and S of Loch Sport provenances (Figure 4b). In Carlos Reyles, where the mean incidence of forking ranged from 0.16 to 0.87, all provenances showed low incidence of forking, except Crediton SF, Bunya SF and Maryland provenances (Figure 4c). In the Vergara trial, the mean incidence of forking ranged from 0.14 to 0.97. The Bunya SF provenance had the lowest number of forks, although it only significantly differed from the Maryland provenance (Figure 4d).

Stem straightness
As mentioned above, there was no provenance by site interaction for stem straightness. This trait showed little differences among provenances, except for Loch Sport and Yurammie SF, which clearly had very bad form (twisted stems) ( Figure 5).

Discussion
It should be noted that the provenances evaluated in this study are represented by few mother trees, and that some, such as SW of Mount Garnet or N of Woolgoolga, come from a single mother tree, which could be a limitation when characterizing their productive performance. Even so, both in the adaptive and growth traits (survival and VOL, respectively), and in the traits related to stem form (straightness and incidence of forking) there were statistically significant differences among provenances. In addition, all traits evaluated, except stem straightness, presented significant provenance by site interaction, indicating that the relative performance of provenances varies even with minor differences in the climate and/or soil where trees were cultivated. These results confirm that the great variety of environments in which E. tereticornis developed has generated a significant variation in the adaptability and performance of different provenances.
Given that the annual average rainfall is not a limitation for the development of this species in the 4 sites evaluated (1159-1558 mm) and that the temperature is very similar between the sites, the differences in productive traits across the sites could be explained by the characteristics of the soils. In sites with deep soils, light textures and good drainage, which allow a good root development, the survival and VOL were greater than in sites with more superficial soils, heavier textures and worst drainage. In contrast, the incidence of forking was lower at sites with shallower and heavier textured soils. These results demonstrate the importance of considering the soil as a key factor for the performance of this species, mainly when the rainfall is not a limiting factor.
Although no clear regionalization was observed, in general the provenances with highest growth rate were those from latitudes between 25 and 37ºS (South Queensland and New South Wales). This region has also been reported as the most productive in Zimbabwe and Argentina (2)(11) (14) . On the other hand, in the present work it was observed that all subtropical provenances, of latitudes below 19ºN (North of Queensland), had poor performance in Uruguay. However, it has been reported that the best provenances for Bangladesh, Brazil, Nigeria, Mozambique and from South to East of India are found in that region (2)(10)(15)(16)(17) (18) .
No provenance was the best in all sites to all the traits evaluated. But in terms of individual tree volume, the Buckenbowra SF, Yurammie SF and Bunya SF showed good general adaptability. Buckenbowra SF also presented adequate survival in all sites, so it can be considered as a source of seeds of high performance in Uruguay. The provenances differed markedly in the incidence of forking, highlighting Cardwell, Kalpower, Burdekin River, and W of Warwick, which had low values in all sites. In general, the stem straightness was adequate in all sites evaluated, with small differences among provenances, except for Loch Sport and Yurammie SF, which presented very poor form in Uruguay.

Conclusions
The results of this study showed the importance of correctly choosing the source of seeds for commercial plantations of E. tereticornis. In the different trials analyzed, the best provenance had a tree growth (individual volume) 3 to 5 times larger than the worst provenance. However, the absence of a clear regionalization and the existence of a significant provenance by site interaction demonstrate the importance of carrying out the local evaluation of different provenances to be able to choose the best source of seeds. These results also suggest that it is necessary to evaluate more provenances of E. tereticornis, mainly from those areas identified as of better aptitude, like South Queensland and New South Wales, to identify sources of seeds that simultaneously present high productivity and good stem form for the cultivation of this species in humid subtropical climates such as in Uruguay.