Growth rates of Eucalyptus and other Australian native tree species derived from seven decades of growth monitoring

•••

ORIGINAL PAPER

Growth rates of Eucalyptus and other Australian native tree species derived from seven decades of growth monitoring

Michael R.Ngugi1•David Doley2•Mark Cant3•Daniel B.Botkin4

©Crown copyright 2015

There is widespread interest in estimating and forecasting individual tree and forest growth rates for restoration and carbon sequestration objectives.Outside intensively managed forests,pastattempts have been limited by the lack of accurate long-term monitoring in multi-age mixed native forests to provide estimates of both expected mean diameter increments and the statistical variation in those estimates.A dataset fromEucalyptus-dominated native forests in subtropical Queensland,Australia offers an opportunity to provide accurate estimates of tree and forest growth rates.Over86,400 treesfrom 155 native specieswere identif i ed and remeasured between 1936 and 2011 in 641 permanentsample plotsacrossa 500–2000 mm mean annual rainfall gradient.Individual tree diameter at breast height (DBH)increments observed for all species ranged mainly from 0.01 to 0.5 cm yr-1(94%of values),with consistent differences between rainfall zones(mean of 500–2000 mm yr-1),and varying differences between species (155)and stem diameter class(10–100 cm).For some species,diameter increment increased progressively with rainfall(e.g.Eucalyptus siderophloia,Eucalyptus propinqua, andLophostemon confertus),but in others(e.g.Corymbia citriodorasubsp.variegata,Corymbia intermedia,andEucalyptus biturbinata)the greatest diameter increments were recorded between 1200 and 1600 mm yr-1.Where there were suff i cientdata,mostspeciesexhibited a quadratic relationship between DBH increment and DBH class,but two species(Callitris glaucophyllaandEucalyptus crebra) native to the 500–800 mm annual rainfall zone showed linear increases in DBH increment with increasing DBH. Continued monitoring of these plots would add to their already great value.

Diameter increments⋅Ecological restoration⋅Forest management⋅Eucalypt forests⋅Natural forests

Introduction

Accurate data on growth rates of individual species are critical for:informing ecological restoration;estimating changes in tree biomass and timber production potential; facilitating species selection for seed source augmentation and ecosystem rehabilitation;and also for determining commercial cutting cycles,recovery time frames and silvicultural prescriptions(Korning and Balslev 1994;Silva et al.2002;Bowman et al.2013;Stephenson et al.2014). Furthermore,in the last decade,tree planting for ecological restoration has been used as a greenhouse gas emissions offset(carbon farming)option by individuals,governments,and corporations in efforts to mitigate climatechange(Bekessy and Wintle 2008;United Nations 2008; Booth and Williams 2012).Additionally,acquisition,protection and management of native regrowth forests complemented with the enrichment planting of native species has also been used as an offset option intended to mitigate biodiversity losses associated with mining and infrastructure development(Commonwealth of Australia 2012).

There are at least 27 forestry-based carbon offset providers in the Australian market(Low Carbon Australia and the University of Queensland 2011)and similar initiatives can be found worldwide with over 200 forestry-based carbon offset providers in the marketplace(The Ecosystem Marketplace 2011).Most carbon offset plantings in Australia use propagules of locally occurring native plant species with the aim of delivering multiple benef i ts that include carbon storage and restoration of biodiversity values of degraded ecosystems(Atyeo and Thackway 2009; Booth and Williams 2012;Booth et al.2012).

Currently,estimation of the benef i ts and risks of all these activities is hampered by a dearth of knowledge on rates of growth of native tree species.There is a vital place for models of forest dynamics(Botkin 1993;Ngugi et al. 2000;Ngugi and Botkin 2011;Ngugi et al.2011)in the prediction of tree and forest growth in different environments and under different management regimes.Forest growth models depend on models of individual tree development,and these in turn depend on a reliable knowledge of the rates of stem diameter growth and the associated changes in height,stem volume and tree biomass(Stephenson et al.2014).Therefore,stem diameter increment data are fundamental to sustainable forest management and carbon farming initiatives.

Although we present analyses for Australian forests of native tree species,these f i ndings have broad application because Australian native tree species are cultivated and studied globally.For example,there are over 20 million hectaresofEucalyptusplantingsin many countriesincluding Brazil,India,China,Uruguay,Chile,South Africa,Portugal, Spain,Vietnam,Sudan and Thailand(GFI Forestry Consulting 2009;Ferraz Filho et al.2014;Wu et al.2014).The successes of these initiatives continue to be limited by a dearth ofinformation on growth rates and growth patterns of numerous Australian native tree species that have potential but are of limited natural occurrence and cultivation(Ngugi et al.2003,2004;Dale and Dieters 2007).

In this study,we investigated mean annual diameter at breast height(DBH)increments of tree species in native forests receiving about 500–2000 mm mean annual rainfall in subtropical Queensland,Australia.We used diameter measurements made over periods up to 75 years (1936–2011)on more than 86,400 trees in 641 permanent forest plots.Our specif i c objective was to determine absolute estimates and the pattern of variation in annual stem diameter increment of selected species across vegetation types,diameter classes and among rainfall zones.

Materials and methods

Study area and location of plots

Our study area covered about 2.6 million hectares of State owned uneven-aged mixed speciesnative forestsin southern Queensland,between latitudes 21°and 29°S and longitudes 146°and 154°E.The forests ranged in elevation from 14 m.a.s.l on Fraser Island to 900 m.a.s.l at Blackdown Tableland National Park in central Queensland(Fig.1).The regionalclimate issubtropicalto semi-arid,characterised by hot humid summers(December–February)and cool,dry winters(June–August).Mean annual rainfall within the study area varies from 500 to 2000 mm.These subtropical forests in Queensland range from tall moist forests along the eastern coast(e.g.Mapleton and Ringtail State Forests,and Bauple National Park)to dry woodlands of central Queensland(Pluto and Umbercollie State Forests)(Tables 1,4). Spatialdata forrainfallzones500–600,600–800,800–1000, 1000–1200,1200–1600 and 1600–2000 mm were provided by Australian Bureau of Meteorology(BOM)based on data collected between 1961 and 1990(BOM 2012).Rainfall zones are identif i ed here by the upper value.

History of management of native forest permanent plots

In the past 100 years these forests have been used for multiple purposes,including timber production and conservation of biological diversity(Queensland Government 1959)They are now held under various tenures,including Conservation Parks(3506 ha),Forest Reserves (120,711 ha),National Parks(536,387 ha)and State Forests(1,949,200 ha).The management goals for State forests in Queensland,as in all States and Territories in Australia,include perpetual protection of watersheds and biodiversity,and the supply of forest products(Queensland Government 1959).

Silvicultural treatments such as thinning of non-commercial trees and selective logging occurred during the history of these forests.Logging practices have included selective harvesting on 20–30-year cutting cycles for sawn timber,poles,girders,railway sleepers and fencing.Silvicultural guidelines determined the stem form,spacing and animal habitat quality of trees to be retained.In the lower rainfall zones(500–1200 mm),commercial harvesting usually used the single tree selection method,and in the wetter zones(1600 mm and above)group selection wasmore commonly practiced(Florence 1996).Prescribed burning and cattle grazing have been part of the f i re risk management strategy in the eucalypt forests,but f i re has generally been purposely excluded from the f i re-sensitiveCallitrisforest communities.Since 1995,these guidelines have been consistent with the principles of the Montreal Process(Montreal Process 2009).

Fig.1 Average annual rainfall isohyets and distribution and location of permanent forest plots(white circles)in State forests and National parks (black outlines)in southern and central Queensland,Australia.A legend for forest estate names and plot counts is provided in Table 4

Broad vegetation classif i cation of native forests

The Queensland Herbarium is the agency now primarily responsible for mapping and classif i cation of vegetation in the study area(Neldner et al.2012).Their classif ication system is hierarchical,ranging from broad classes (～1:5 million scale)down to regional ecosystems(RE) suitable for mapping at about 1:100,000 scale.These RE are combinations of vegetation types which are consistently associated with particular combinations of geology,landform,soil and climate(Neldner et al.2012). We used a classif i cation of broad vegetation groups (BVG)with nominal scale of 1:2 million,which divided the sub-tropical forests in this study into 16 broad vegetation groups(BVGs),ranging from tall moist forests along the eastern coast to dry woodlands of centralQueensland(Table 1).Each plot was individually assigned to a BVG based on the dominant tree species, basal area stock of dominant species,height and foliage projected cover of the forests recorded at the site (DERM 2012).Consequently,this vegetation classif i cation provided a detailed methodology of aggregating plots on moisture availability,basal area,species combination and growth conditions.In other studies,stand basal area attribute,a component of this classif i cation, has been used in growth modelling as a proxy for intertree competition(see Weiskittel et al.2011;Prior and Bowman 2014).

Table 1 Brief broad vegetation group(BVG)description of forests represented by permanent plots in the study area,a count of number of plots located in each rainfall zone,and mean diameter increment and 95%conf i dence interval(in italics)for the each BVG and within each rainfall zone

Native forests permanent sample plot data

The tree growth data used in this study were collected from 641 long-term permanent inventory plots located in 118 forest estates(Fig.1;Table 4).The earliest plots were established in 1936 at Braemar State Forest and the last batch in the late 1990s.All permanent plots established prior to 1993 were subjectively selected from a systematic grid network of temporary forest inventory plots,to cover the range of stocking,species composition and site quality of major commercial stands within forest types(Beetson 1992).Plots established after 1993 were based on forest stratif i cation data,to cover under-represented species mixtures and geographic distributions(Cant and Mannes 1995).There were 609 rectangular plots,varying in size from 0.2 to 0.5 ha,within which every tree with DBH＞10 cm that was suitable for commercial wood production was individually numbered,tagged and measured,complying with an established measurement manual (Cant and Mannes 1995).In a separate study on Fraser Island(Applegate 1982),32 circular plots were established, within which trees with DBH＞10 cm were measured in an inner circle of radius of 17.8 m and trees with DBH＞30 cm measured in an outer circle with a radius of 39.9 m.

The measurement protocol was amended in 1993 to include sampling of all tree species irrespective of their potential for timber production.This amendment was aimed at promoting recovery of species diversity in forests used for timber production.In addition,trees＜10 cm DBH but＞3 m tall and or DBH＞2 cm were sampled using small circular plots(≤5 m radius)located at the centre of each quarter of the rectangular plot.All plots irrespective of establishment year were re-measured every 2–10 years after their establishment,with additional measurements carried out after a logging,wildf i re or silvicultural treatment.The total number of times a plot was monitored ranged from two to 18,with an average of f i ve times per plot.

Data analyses

Data analysis was undertaken using R programming language(R Development Core Team 2011).Annual diameter increment was calculated as the change in individual tree diameter between two consecutive measures as a function of time period in years(time period≥1 year).Filtering was done to remove gross typing and recording errors,and errors associated with shifts in height at which diameter was measured due to recorded development of deformities at DBH.Trees with a defect at 1.3 m(DBH measure point) were excluded during diameter increment calculations.The growth rates were calculated for each 10 cm diameter class and expressed as mean,and 95%conf i dence interval(CI) about the mean.Relationships between tree diameter class (D)and diameter increment(Di)for individual species within a rainfall zone were explored graphically.Preliminary model f i tting used different model forms(exponential, logarithmic,power function,polynomial and linear)and aggregations of data from several rainfall zones in an attempt to f i t species-specif i c generalised models.Fitted models were compared using the coeff i cient of determination(r2).Models that best f i tted the data and had r2＞0.5 were retained.

Results

Growth rates of broad vegetation groups and species representation

The ranges of mean diameter increment for forest types among rainfallzones(0.18–0.27 cm yr-1)was less than that among BVGs(0.11–0.29 cm yr-1)(Table 1).The highest mean increment among BVGs(0.29±0.18 cm yr-1)was observed in BVG 8b,moist open-forests to tall open-forests dominated byEucalyptus pilularison coastal sands,subcoastalsandstonesand basaltranges.ForallBVGs,the95% conf i dence interval was about 57%of the mean increment (data not shown).

Mean diameter growth rates and 95%conf i dence intervals of 23 of the most abundant tree species are presented in Table 2.The most abundant data were available forC.glaucophylla(84,517 increment values),C.citriodorasubsp.variegata(18,505 increment values),andE. crebra(16,409 increment values).Mean stem density for trees＞10 cm DBH ranged from 287 to 399 stems ha-1but without any consistent relationship to rainfall.Mean stand basal areas for each rainfall zone,calculated across the entire sampling period,ranged from approximately 7.5 m2ha-1(600 mm rainfall zone)to 29 m2ha-1(1600 mm rainfall zone)(Fig.2).The relatively low stand basal area(15 m2ha-1)in the 1200 mm rainfall zone ref l ects the extent of commercial utilisation in this forest zone.

Species growth rates across diameter size classes within rainfall zones

Comparisons of growth rates across diameter classes for species with the highest number of annual increment determinations within each rainfall zone are presented in Fig.3.Within the 600 mm zone,growth rates ofE.populnea,E.crebraandC.glaucophyllaremained relatively constant(between 0.15 and 0.3 cm yr-1)across the observed diameter range whileE.melanophloiashowed amaximum growth rate(0.34 cm yr-1)in the 30–40 cm diameter class(Fig.3a).In the 800 mm rainfall zone,the growth rates of all frequent species were relatively consistent in stems up to 50 cm DBH but were more variable(0.04–0.49 cm yr-1in larger stems(Fig.3b). Growth rates for the frequent species within the 1000 mm zone ranged from 0.13 to 0.41 cm yr-1(Fig.3c).Within species,there was no consistent pattern of variation in diameter increment with DBH butE.moluccanashowed generally higher increments thanE.siderophloiaandE. propinqua.

Table 2 Mean diameter growth rates(cm yr-1)of selected tree species with＞200 increment data values(n)across observed rainfall zones

Fig.2 Mean stand basal area for rainfall zones calculated from all measurement occasions on all plots.Rainfall zones(range of mean annual total):600 mm(500–600 mm);800 mm(600–800 mm); 1000 mm,(800–1000 mm);1200 mm,(1000–1200 mm);1600 mm, (1200–1600 mm);2000 mm,(1600–2000 mm)

Diameter growth rates for trees in the 1200 mm zone ranged between 0.2 and 0.52 cm yr-1for trees in the 10–70 cm DBH range,with generally consistent and higher rates inE.microcorysand inE.sphaerocarpa,a species conf i ned to Blackdown Tableland National Park(Fig.3d). Highest growth rates in other species were generally observed in the 30–50 cm diameter range.There was far greater variation in growth rate among species in the 1600 mm rainfall zone(Fig.3e).E.grandisandE.pilularisshowed maximum growth rates between about 0.8 and 0.9 cm yr-1)in contrast to about 0.4 cm yr-1in otherfrequent species in the zone.Highest growth rates were observed in the 30–70 cm DBH range.Great variability in diameter growth rates was also observed in the 2000 mm zone,whereE.microcorys,E.racemosaandE.pilularisshowed growth rates between 0.4 and 0.7 cm yr-1within the 30–70 cm DBH range(Fig.3f).

Fig.3 Comparative species diameter increments within each rainfall zone as def i ned in Fig.2

Diameter growth rates across rainfall zones

Over 45 species in our dataset occurred across at least three rainfall zones(Tables 2,5).The growth rates among species and rainfall zones ranged from 0.01 to 0.96 cm yr-1, with an overall mean of 0.25±0.13 cm yr-1.A sample comparison of diameter growth rates of 10 species that occurred across 800 mm to 2000 mm rainfall zones is presented in Fig.4.Relatively constant growth rates across rainfall zones were observed forCorymbia trachyphloiasubs.trachyphloia(0.12–0.16 cm yr-1)but for the other nine species,two patterns were generally observed.Eucalyptus siderophloia,E.propinqua,E.tereticornis,E. acmenoidesandLophostemon confertusdemonstrated increased diameter growth rates(up to 0.27–0.35 cm yr-1), with increasing rainfall(Fig.4a).Corymbia citriodorasubsp.variegata,E.biturbinata,Corymbia intermediaandAngophora leiocarpa,demonstrated peak growth rates (0.21–0.31 cm yr-1)in the 1200 or 1600 rainfall zones (Fig.4b).

Table 3 Regressions of diameter increment(Di,cm yr-1)as a function of tree diameter(D,cm)for species and rainfall zones in which the regression coeff i cient of determination(r2)was＞0.5

Growth rate variation within diameter classes and across rainfall zones

Figure 5 illustrates the different patterns of variation in growth rate with diameter class for 10 common species that were distributed across several rainfall zones.There was little difference in diameter growth rates ofAllocasuarina luehmannii,Callitris glaucophylla,Corymbia intermediaorEucalyptus crebrawithin across rainfall zones and within diameter classes(Fig.5).Eucalyptus acmenoides,E. fi brosasubsp.fi brosa,andE.siderophloiashowed consistently lower growth rates in the lowest observed rainfall zones whileCorymbia citriodorasubsp.variegataand,E. fi brosasubsp.fi brosashowed higher growth rates in the highest rainfall zone of their occurrence.The variation in growth rates with diameter class and rainfall zone inLophostemon confertuswas inconsistent(Fig.5).

Regression analyses of diameter increment as a function of tree diameter are presented in Table 3 for the species illustrated in Fig.5.Although most species were observed across several rainfall zones,only in a few rainfall zones were the coef fi cients of determination values(r2)＞0.5. Species with the highest r2wereEucalyptus fi brosasubsp.fi brosa(r2=0.91),Corymbia citriodora subsp.variegata(r2=0.87)andCorymbia intermedia(r2=0.85)growing in the 1200 mm rainfall zone.

Discussion

Although the desirability of increasing forest carbon stocks in addressing climate change is well established(Bekessy and Wintle 2008;United Nations 2008;Booth and Williams 2012)and there are many estimates and methods for estimating carbon stocks(e.g.IPCC 2006),the rates at which these stocks accumulate in native forests are by comparison known imprecisely(Stephenson et al.2014). Reliable empirical estimates of growth rates and patterns for 74 native tree species growing in natural forests in subtropical Australia were obtained from an extensive individual tree increment dataset(211,189 records). Growth rates and patterns varied markedly among broad vegetation groups and species in the observed rainfall zones and diameter size classes and the patterns were diffi cult to generalise.In the low rainfall zones(＜800 mm) growth rates among species and across diameter classes were generally similar.Much greater f l uctuations in growth rates were observed in the＞1000 mm zones.Individual growth rates among all the species ranged from 0.01 to 0.96 cm yr-1.

Mean annual diameter growth rates of species within each rainfall zone

The generally similar growth rates among frequently occurring species in the low rainfall areas in contrast to the pronounced differences observed in higher rainfall areas suggests that potential growth differences among species are small under low rainfall conditions as compared tohigher rainfall conditions.The highest growth rates for most species in all rainfall zones were observed in trees with DBH between 20 and 60 cm.In these native forests, trees in the＜20 cm DBH class are often in the understorey and their growth is often suppressed by larger trees. Moreover,even though small trees may be numerous in a stand,their contribution to stand basal area and carbon stock is small compared with that of larger trees(Stephenson et al.2014).Compared to other studies,the growth rates observed in our study were considerably lower than those of 1–5 cm yr-1reported from 4 year-old plantation trial of eucalypts in the 1000 mm rainfall zone in southeast Queensland(Ngugi et al.2004).However,unlike the natural growth conditions in our study,growth conditions under experimental trials or commercial plantations are modif i ed through intensive ground preparation procedures involving deep ripping,fertilization and weed control (Ngugi et al.2004)and the absence of competition for light,water and nutrients from established(old)trees.

Fig.4 Variation in diameter growth rate of individual species with diameter class across 800–2000 mm rainfall zones with a showing species that demonstrated highest growth at highest rainfall and b species that demonstrated a peak and a decrease in diameter growth rates

Fig.5 Diameter increments across diameter size classes for ten species with distributions across several rainfall zones.Equations for the regression lines are presented in Table 3

The observed mean diameter increment forCorymbia citriodorasubsp.variegata(0.25±0.10 cm yr-1, n=18,493)is consistent with results obtained for the species in southeast Queensland and northern New South Wales that ranged from 0.20±0.03 in the lower rainfall areas,up to 0.29±0.02 in the higher rainfall areas (Lewis et al.2010).The mean individual tree growth rate range for all trees used in this study was 0.01–0.5 cm yr-1for 94%of the data,with an overall mean of 0.25 cm yr-1.This range is consistent with that reported forEucalyptusforests in temperate Australia(Bowman et al.2014)but the mean value is higher than that for woodland(0.20 cm yr-1)and for open forests (0.16 cm yr-1)reported in north Australia savanna in a two-year study(Prior et al.2004).Our mean is also within the range of 0.2–0.37 cm yr-1reported for subtropical forest trees in Puerto Rico(Brandeis 2009)and close to the 0.27 cm yr-1reported for tropical lowland tree species in the Amazonian Ecuador where the mean annual precipitation was 3244 mm(Korning and Balslev 1994).

Tree plantings for carbon capture often target species with the highest growth prospect within a rainfall area and geographic location.However,because some species have a wide geographical distribution,it is often diff i cult to assess relative environmental requirements and growth potential(Booth et al.2012).For example,many species observed in our dataset occur across several rainfall zones (between 600 and 2000 mm).Because there is limited knowledge on where such species are likely to exhibit their highest growth rates,the selection of fast growing species for carbon accumulation for a particular local area is diffi cult.Using a sample of ten species that occurred between the 800 and 2000 mm rainfall zones,we attempted to identify the zone in which each of the species recorded the highest growth rates.For example,althoughCorymbia intermediais widespread,the growth increment recorded in the 1200 mm zone was approximately four times higher than that recorded in the 800 mm zone and about 55% higher than that recorded in the 1600 and 2000 mm rainfall zones.Contrary to the f l uctuating growth increment patterns observed for the other eight species,there was a consistent three-fold increase in growth rates ofE.siderophloiaandE.acmenoideswith increase in rainfall from 800 to 2000 mm.

Growth rates of individual species across diameter classes among rainfall zones

Growth rates ofC.glaucophyllawere relatively independent of tree diameter for trees greater than 10 cm DBH, with growth rates of trees＜10 cm DBH being about half that of the larger trees.A similar pattern was observed forEucalyptus crebragrowing between the 800 and 1000 mm zones.However,in several species,includingC.citriodorasubsp.variegata,Lophostemon confertus,E.acmenoides,E.siderophloia,Corymbia intermediaandE.f i brosasubsp.fi brosa,there was an increase in growth rate with increase in diameter up to a DBH of 50 cm,generally followed by a decrease in larger trees.The large differences in growth rates across diameter classes suggest that the use of an average increment value for these species should be adopted cautiously and only when more detailed descriptions of increment are not available.

Natural multi-age forests in Australia are characterised by a large number of species and varying growing conditions.This study did not document all possible tree species and growing conditions but instead used a dataset that was initially aimed at collecting data for the purpose of estimating timber production in State Forests in Queensland. As a result,commercial timber species are well represented in the dataset relative to those perceived to have limited timber value.Increment estimates for some non-commercial species are now being obtained as a result of the changes in 1993 to the plot measurement protocol.Some of these non-commercial species indicate growth rates much higher than the over-storey eucalypt species(e.g. 0.96 cm yr-1forAcacia melanoxylonin the 2000 mm zone,n=11;and 0.76 cm yr-1forAngophora f l oribundain the 1600 mm zone,n=6).The accuracy of estimates for understorey and short-lived species is currently limited by small sample sizes but will be enhanced as subsequent re-measurements occur.Estimates of growth rates for large trees(＞80 cm)should be used very cautiously because of measurement errors resulting from bole deformities associated with aging of predominantly habitat trees and bias related to their limited representation in the dataset.

Conclusions

Stem diameter increments for native species growing naturally in uneven-aged mixed species forests in subtropical Queensland typically range from 0.01 to 0.50 cm yr-1(94%of observed values).While tree growth rate is inf l uenced by arangeofgenetic,edaphic and climaticfactors,thisstudy related increment to growth conditions(geology, landform,soil and rainfall),stem diameter(tree size)and rainfall,which are the most commonly recorded factors in a region where temperature variation is relatively small.In most species,increment is related to both rainfall and stem diameter class,but the relationships vary substantially between species and sometimes between rainfall zones.In some species,particularly from lowerrainfallenvironments, increment does not vary consistently between either rainfall or stem diameter classes.For species that occurred across several rainfall zones,the highest growth rates were not necessarily observed in the highest rainfall zone.Within a rainfall zone,the highest growth rates within species were often observed in trees within the DBH range of 20–60 cm, following a unimodal curve within that range or with the increment either increasing linearly with stem diameter (often in lower rainfall zones and especially inCallitris glaucophyllaandEucalyptus crebra).

The dataset used for this study provided a consistent record of individual tree growth that avoids the errors and costs of repeatedde novoassessments.However,many of the plots have not been remeasured in the last decade.With continued monitoring,the network of plots could provide a robust scientif i c basis for directing sustainable forest management and assist in enhancing our understanding on carbon storage,impacts of rising atmospheric CO2concentration on tree growth,forest dynamics and improve accuracy of forest growth prediction tools.

AcknowledgmentsThe data used for this study are stored in a database maintained by the Queensland Herbarium,Department of Science,Information Technology and Innovation(DSITI)and cojointly held by the Forest Research,Department of Agriculture, Fisheries and Forestry(DAFF).The legacy and contribution of past Queensland Government Forestry Departments and staff in data collection,collation and maintenance for over seven decades is greatly acknowledged.Assistance provided by Jiaorong Li and Rosemary Niehus in GIS,Queensland Parks and Wildlife Service fi eld staff,Jian Wang,David Moore,Rosemary Niehus and Sue Philips for 2011 f i eld measurements is gratefully acknowledged.Peer review and suggestions provided by Dr Don Butler,Dr Arnon Accad and Dr John Neldner for earlier version of this manuscript are acknowledged with thanks.

Appendix

See Tables 4 and 5.

Table 4 Assigned forest identif i cation numbers(FID),gazetted area and number of native forest permanent plots located in each forest,shown in Fig.1

Table 4 continued

Table 5 Mean diameter growth rates(cm yr-1)of species across rainfall zones

Table 5 continued

Table 5 continued

Applegate GB(1982)Biomass of blackbutt(Eucalyptus pilularisSm.)forests on Fraser Island.Masters Thesis,University of New England,Armidale

Atyeo C,Thackway R(2009)Mapping and monitoring revegetation activities in Australia:towards national core attributes.Aust J Environ Manag 16:140–148

Beetson T(1992)Enhancing the native forest detailed yield plot system.Queensland Forest Service,Brisbane,p 200

Bekessy SA,Wintle BA(2008)Using carbon investment to grow the biodiversity bank.Conserv Biol 22:510–513

BOM(2012)Annual average rainfall in Queensland based on a standard 30 year climatology record(1961–1990).Australian Bureau of Meteorology,Canberra

Booth TH,Williams KJ(2012)Developing biodiverse plantings suitable for changing climate conditions 1:underpinning scientif i c methods.Ecol Manag Restor 13:267–273

Booth TH,Williams KJ,Belbin L(2012)Developing biodiverse plantings suitable for changing climate conditions 2:using the atlas of living Australia.Ecol Manag Restor 13:274–281

Botkin DB(1993)Forest dynamics:an ecological model.Oxford University Press Inc.,New York,p 309

Bowman DMJS,Brienen RJW,Gloor E,Phillips OL,Prior LD(2013) Detecting trends in tree growth:not so simple.Trends Plant Sci 18:11–17

Bowman DMJS,Williamson GJ,Keenan RJ,Prior LD(2014)A warmer world will reduce tree growth in evergreen broadleaf forests:evidence from Australian temperate and subtropical eucalypt forests.Glob Ecol Biogeogr 23:925–934

Brandeis TJ(2009)Diameter growth of subtropical trees in Puerto Rico.Forest Service,Southern Research Station,Asheville,p 39 Cant M,Mannes D(1995)Native forest permanent plot establishment.Brisbane,Department of Primary Industries,Brisbane, p 69

Commonwealth of Australia(2012)Environment Protection and Biodiversity Conservation Act 1999:environmental offset policy.Department of Sustainability,Environment,Water, Population and Communities,Canberra,p 32

Dale G,Dieters M(2007)Economic returns from environmental problems:breeding salt-and drought-tolerant eucalypts for salinity abatement and commercial forestry.Ecol Eng 31:175–182

DERM(2012)Broad vegetation groups of Queensland.Brisbane, Australia.Queensland government,Queensland Herbarium. http://www.ehp.qld.gov.au/ecosystems/biodiversity/regional-eco systems/bvg.html.Accessed 5 May 2012

Ferraz Filho AC,Scolforo JRS,Mola-Yudego B(2014)The coppicewith-standards silvicultural system as applied toEucalyptusplantations:a review.J For Res 25:237–248

Florence RG(1996)Ecology and silviculture of eucalypts forests. CSIRO,Collingwood,p 413

GFI Forestry Consulting(2009)GlobalEucalyptusmap.A cartographic information resource depictingEucalyptuscultivated forests worldwide.GIT Forestry Consulting SL,Lugo

IPCC(2006)International panel on climate change guidelines for national greenhouse gas inventories.IGES,Hayama

Korning J,Balslev H(1994)Growth rates and mortality patterns of tropical lowland tree species and the relation to forest structure in Amazonian Ecuador.J Trop Ecol 10:151–166

Lewis T,Osborne D,Hogg B,Swift S,Taylor D,Macgregor-Skinner J(2010)Tree growth relationships and silvicultural tools to assist stand management in private native spotted gum dominant forests in Queensland and northern New South Wales,technical report of forest and wood products Australia

Low Carbon Australia and the University of Queensland(2011) Carbon offsets guide Australia.http://carbonoffsetguide.com.au. Accessed July 2013

Montreal Process(2009)Criteria and indicators for the conservation and sustainable management of temperate and boreal forests, p 29.http://www.rinya.maff.go.jp/mpci/2009p_4.pdf.Accessed 6 Nov 2012

Neldner VJ,Wilson BA,Thompson EJ,Dillewaard HA(2012) Methodology for survey and mapping of regional ecosystems and vegetation communities in Queensland.Version 3.2 updated August 2012,Queensland Herbarium,Queensland Department of Science,Information Technology,Innovation and the Arts, Brisbane.http://www.ehp.qld.gov.au/plants/herbarium/publica tions/pdf/herbarium_mapping_methodology.pdf.Accessed 22 Nov 2013

Ngugi MR,Botkin DB(2011)Validation of a multispecies forest dynamics model using 50-year growth fromEucalyptusforests in eastern Australia.Ecol Model 222:3261–3270

Ngugi MR,Mason EG,Whyte AGD(2000)New growth models forCupressus lusitanicaandPinus patulain Kenya.J Trop For Sci 12:524–541

Ngugi MR,Hunt MA,Doley D,Ryan P,Dart P(2003)Dry matter production and allocation inEucalyptus cloezianaandEucalyptus argophloiaseedlings in response to soil water def i cits.New For 26:187–200

Ngugi MR,Hunt MA,Doley D,Ryan P,Dart P(2004)Selection of species and provenances for low-rainfall areas:physiological responses ofEucalyptus cloezianaandEucalyptus argophloiato seasonal conditions in subtropical Queensland.For Ecol Manag 193:141–156

Ngugi MR,Johnson RW,McDonald WJF(2011)Restoration of ecosystems for biodiversity and carbon sequestration:simulating growth dynamics of brigalow vegetation communities in Australia.Ecol Model 222:785–794

Prior LD,Bowman DMJS(2014)Across a macro-ecological gradient forest competition is strongest at the most productive sites.Front Plant Sci 5:1–12

Prior LD,Eamus D,Bowman DMJS(2004)Tree growth rates in north Australian savanna habitats:seasonal patterns and correlation with leaf attributes.Aust J Bot 52:303–314

Queensland Government(1959)Forestry Act.Queensland Government,L Queensland,p 256

R Development Core Team(2011)R:a language and environment for statistical computing.R Foundation for Statistical Computing, Vienna

Silva RP,Santos J,Tribuzy ES,Chambers JQ,Nakamura S,Higuchi N(2002)Diameter increment and growth patterns for individual tree growing in Central Amazon,Brazil.For Ecol Manage 166:295–301

Stephenson NL,Das AJ,Condit R,Russo SE,Baker JP,Beckman NG,Coomes DA,Lines ER,Morris WK,Ruger N,Alvarez E, Blundo C,Bunyavejchewin S,Chuyong G,Davis SJ,Duque A, Ewango CN,Flores O,Franklin JF,Grau HR,Hao Z,Harmon ME,Hubbell SP,Kenfack D,Lin Y,Makana JR,Malizia A, Malizia LR,Pabst RJ,Pongpattananurak N,Su S-H,Sun I-F,Tan S,Thomas D,van Mantgem PJ,Wang X,Wiser SK,Zavala MA (2014)Rate of tree carbon accumulation increases continuously with tree size.Nature 507:90–93

The Ecosystem Marketplace(2011)Carbon catalog.http://www. carboncatalog.org/about/.Accessed 10 Mar 2014

United Nations(2008)The United Nations collaborative programme on reducing emissions from deforestation and forest degradation in developing countries.http://www.un-redd.org/Home/tabid/ 565/Default.aspx.Accessed 5 May 2014

Weiskittel AR,Hann DW,Kershaw JA Jr,Vanclay JK(2011)Forest growth and yield modelling.Wiley,Chichester

Wu PF,Ma XQ,Tigabu M,Huang Y,Zhou LL,Cai L,Hou XL,Oden PC(2014)Comparative growth matter accumulation and photosynthetic rate of seven species of Eucalypt in response to phosphorus supply.J For Res 25:377–383

21 May 2014/Accepted:15 December 2014/Published online:16 July 2015

The online version is available at http://www.springerlink.com

Corresponding editor:Hu Yanbo

✉Michael R.Ngugi

michael.ngugi@dsiti.qld.gov.au

1Ecological Sciences,Queensland Herbarium,Department of Science,Information Technology and Innovation, Queensland Government,Mt Coot-tha Road,Toowong, QLD 4066,Australia

2Centre for Mined Land Rehabilitation,Sustainable Minerals Institute,The University of Queensland,Brisbane, QLD 4072,Australia

3Technical Services,Pest Management,Department of National Parks,Recreation,Sport and Racing,Queensland Government,173 Hume Street,Toowoomba,QLD 4350, Australia

4Department of Biology,University of Miami,Coral Gables, FL 33124-0421,USA