Vulnerability of Plantation Carbon Stocks to Defoliation under Current and Future Climates

Plantation species globally are susceptible to a range of defoliating pests, but pest damage is rarely considered when estimating biomass C sequestered by these forests. We examined the impacts of defoliation on Eucalyptus globulus plantation C stocks under current and future climates using Mycospharella Leaf Disease (MLD) as a case study, hypothesising that biomass C sequestered in plantations would decrease with a warming and drying climate, and that impacts of defoliation would be strongly site dependent. Six E. globulus plantation sites with varying productivity were selected for the study. Current (1961–2005) and future (2030 and 2070) severity and frequency of MLD were estimated for each site using the bioclimatic niche model CLIMEX, and used as inputs to the process-based forest productivity model CABALA. CABALA was used to develop annual estimates of total living and dead biomass for current, 2030 and 2070 climate scenarios. Averaged annual biomass outputs were used to initialise the carbon accounting model FullCAM for calculation of C sequestered in living and dead biomass over a growing cycle. E. globulus plantations were predicted to sequester between 4.8 and 13.4 Mg C·ha·year over 10 years under current climatic conditions. While our estimates suggest that overall this is likely to increase slightly under future climates (up to a maximum of 17.2 Mg C·ha·year in 2030, and a shift in minimum and maximum values to 7.6 and 17.6 respectively in 2070), we predict considerable between-site variation. Our results suggest that biomass C sequestration will not necessarily be enhanced by future climatic conditions in all locations. We predict that biomass C sequestration may be OPEN ACCESS Forests 2014, 5 1225 reduced considerably by defoliation meaning that any gains in C sequestration associated with changing climate may be substantially offset by defoliation. While defoliation has a generally small impact under current climatic conditions in these plantations, the impact is likely to increase in the future, with reductions of up to 40% predicted for some sites under future climates. We conclude that the combined impacts of climate change on pest frequency and severity, and on host responses to defoliation, may reduce biomass C sequestration in E. globulus plantations in the future.


Introduction
Forests are vulnerable to a range of disturbance events, including attack from insects and pathogens, and altered environmental conditions associated with changing climate.This vulnerability is well illustrated by recent outbreaks of mountain pine beetle (Dendroctonus ponderosae) in northern America, where extensive tree mortality and reductions in forest net primary production associated with D. ponderosae activity converted the forests from a small carbon (C) sink to a large C source [1] and changed species composition [2].Similarly, in Eucalyptus delegatensis forest in Australia, the interaction of drought and pest attack (Cardiaspina spp.) reduced the C gain of the forest, converting it from a sink to a source of C in the short term [3].In both of these examples, climate perturbations were the primary trigger for the pest outbreaks.
The impact of attack by insects and pathogens on forests and the ecosystem services they provide is related to forest vigour, pest dynamics, the interactions between host and pest, and natural enemies [4][5][6].Because of this complexity, studies of the impacts of pest attack on the productivity and C stores of managed forests are rare, with studies often short term in nature [7][8][9][10][11][12][13][14], and not extending to at-harvest consequences.Where modelling has been used to examine impact, results suggest large variation in impact related to the forest management strategies imposed [15] or variation in site conditions [16].These studies point to a likely increase in impact of pests on forest production under future climates [15,16].
In addition to possible changes in the way trees respond to pest attack under changing climates, there is considerable scope for pest distribution and activity to change.There is ample evidence that changing climate, and particularly warmer mean temperatures, will affect the range of pest species, as well as patterns of mortality and lifecycle traits [17][18][19][20][21], although predicting with certainty the likely effects of climate change on specific insects or pathogens is difficult due to the complex way in which they interact with their host species, enemies, competitors and symbionts [22].Bioclimatic species distribution models provide a tool for examining potential range shifts and changes in activity of pest species associated with changing climate [23][24][25], and the outputs of such models can be used as inputs to forest productivity models to provide greater realism to modelled impact [16,26].
When planted on cleared land, trees can provide a substantial additional carbon (C) sink [27], and C sequestration in reforestation is an acceptable method of mitigation under Kyoto Article 3.3 [28].Carbon sequestered in living and dead biomass and soil is allowable under the Kyoto Protocol, and under reforestation the living and dead biomass commonly accounts for the greatest C component (e.g., ~60% for a range of eucalypt species growing in 26-year-old stands, [28]).Estimates of annual biomass C sequestration by plantation forests in Australia are around 21 Mt CO 2 -e.Published values for individual plantations in low rainfall or disturbed environments in Australia range between 1.9 and 4.7 Mg C•ha −1 •year −1 [29,30], but published values for plantations in higher rainfall zones are rare.When defining the capacity of plantations to sequester C, it is important to understand the impacts of disturbance events such as pest damage on biomass C sequestration, in order to define the risks posed by disturbance and to develop strategies for managing those risks.
One of the main temperate plantation species globally, Eucalyptus globulus Labill., is susceptible to a range of defoliating pests [31].Many of these pests are likely to exhibit a polewards shift in range under future climates, and increased activity within that range [32].The extent to which E. globulus plantation C sequestration is currently affected by pest activity and how this might change in the future, are largely unknown.
We addressed the question of the vulnerability of biomass C stocks to disturbance by examining the impacts of defoliation under current and future climates.We used Mycosphaerella leaf disease (MLD) in Australia as a case study because Mycosphaerella species are the most significant fungal pathogens causing defoliation of eucalypts in temperate plantations globally [33][34][35][36][37]. Symptoms of MLD include red discoloration and the development of necrotic lesions on leaves followed by premature leaf senescence [38,39].It can reduce tree growth [37,40], and in severe cases causes tree mortality.Projected changes in distribution of MLD, modelled using the bioclimatic niche model CLIMEX [25] suggest increased climatic suitability for MLD with changing climate across much of the E. globulus plantation estate in Australia [32].The impacts of projected changes in MLD severity, frequency and seasonality on E. globulus wood production have been previously examined by incorporating the CLIMEX outputs as inputs into the forest growth model CABALA [16,41].
In this paper, we used CLIMEX outputs to develop defoliation scenarios in the forest productivity model CABALA [42], and CABALA outputs to initialise the carbon accounting model FullCAM [43] to estimate the impacts of MLD defoliation on E. globulus plantation C sequestration under current and future climates.Because MLD can strongly affect stand biomass, we focused on its effects on biomass C while recognising that disturbance by defoliation may also affect soil C. We hypothesised that C sequestered in living and dead biomass would increase with changing climate, and that vulnerability of that C to a disturbance such as defoliation would be strongly site dependent.We focused on six sites across a rainfall gradient in southern Australia.

Methods
The framework used for the study is summarised in Figure 1.We examined the interaction of exposure of plantations to pests and climate, and sensitivity of plantations in terms of the way they responded to that exposure, to determine the impact on biomass C sequestration.
Specifically, we took published estimates of changes in MLD severity and frequency developed for the six E. globulus sites described below (Pinkard et al., 2010, [16]), and used these to develop regime files for those sites in the forest productivity model CABALA [42].The relative error between observed and predicted stand volume for the six sites was between 1% and 6%.We developed climate files as described below and used these as CABALA inputs.We then used CABALA estimates of stand biomass production as inputs into the carbon accounting model FullCAM to develop site-level estimates of living and dead biomass C [43].

Figure 1.
The framework used for the analysis of impact of defoliation on biomass carbon sequestration in eucalypt plantations.Exposure to pests was defined using previously published data on severity and frequency of pest outbreaks at the sites selected for the study [16].

Climate Data/Inputs
Daily current climate data were obtained for the study sites using a SILO data drill [44] for years 1961-1991 (Hadley Mk2) and 1975-2005 (CSIRO Mk3).Patterns of monthly-average change in temperature and rainfall were generated for 30-year periods centred on 2030 and 2070 using the Hadley Mk2 and CSIRO Mk3 climate models and used to modify the current climate base data.The models were selected to allow analysis of the impacts of a more moderate (i.e., less warming and drying-Hadley) and a more extreme (i.e., warmer and drier-CSIRO) future climate on biomass C sequestration, to provide indications of best and worst case scenarios.The A1FI emission scenario was applied to the Hadley model, with atmospheric CO 2 rising from 360 ppm under current conditions to 450 ppm for 2030 and 660 ppm for 2070 climates.The A2 emission scenario was applied to the CSIRO model, with atmospheric CO 2 rising from 360 ppm under current conditions to 450 ppm for 2030 and 550 ppm for 2070.This provided best and worst-case scenarios given our generally poor understanding of levels of photosynthetic acclimation to elevated CO 2 in E. globulus growing under a range of site conditions [45,46].

Modelling the Impact of Pest Attack on Forest Biomass
Six E. globulus sites with predicted productivity varying from final stand volume of 121 to 404 m 3 •ha −1 (Table 1) were selected to examine the potential impacts of climate change and defoliation on at-harvest carbon sequestration of standing biomass and litter.A single pest species, Mycosphaerella leaf disease (MLD), was selected.Severity and frequency of MLD epidemics under current, 2030 and 2070 climate projections were estimated using the bioclimatic niche model CLIMEX [47], as outlined in Pinkard [48] and summarised in Table 2.These CLIMEX outputs were used as inputs to the process-based forest productivity model CABALA (Figure 1) [41,42].CABALA links carbon, nitrogen, water and nitrogen cycles so that simultaneous impacts of leaf area loss on biomass allocation, transpiration and water use, and nitrogen reallocation can be assessed.It was originally calibrated for E. globulus.Modifications to CABALA to extend the light interception and photosynthetic sub-modules to deal with discontinuous crowns make the model ideally suited to assess the impacts on productivity of leaf area damage and loss.CABALA allows the user to specify where in the crown defoliation and discoloration occur [42].Defoliation associated with MLD tends to be first observed in the inner crown, with concomitant necrosis occurring in the remaining crown, and is mainly observed on juvenile foliage that is present in the crown until 3-5 years of age [49].This was mimicked as closely as possible in the input files.For low MLD severity scenarios, MLD was specified to develop in year 1.For moderate severity, MLD was specified to develop in years 1 and 3, and for severe MLD, infection was specified to develop in years 1, 2 and 3.Because of some uncertainty around the level of leaf loss that may be experienced for each of the MLD severity classes, defoliation scenarios were generated assuming either 10%, 20% or 30% leaf loss for low severity, 40%, 50% or 60% defoliation for moderate severity, and 60%, 70% or 80% defoliation for the severe class.Clearfell was set at 10 years, which is a routine rotation length for E. globulus pulpwood plantations in Australia.CABALA was used to develop annual estimates of total standing biomass, and the fraction of that biomass in leaves, branches, stem, coarse and fine roots.The litter component also was estimated.The model was run for current, 2030 and 2070 climate scenarios, with 10 different starting dates per site and climate scenario to allow for inter-annual variation in weather.The CSIRO-Mk3-A2 and Hadley-Mk2-A1FI modelled weather sequences were used.In addition, estimates were made assuming either no acclimation of photosynthesis to elevated CO 2 (sensu Ainsworth and Long [50]) or full acclimation to elevated CO 2 , with assumed atmospheric [CO 2 ] of 360 ppm for current conditions, 450 ppm for 2030 conditions and for 2070 either 550 or 600 ppm for the CSIRO and Hadley GCM's respectively.

Modelling Impact of Pest Attack on Sequestered Carbon
Averaged annual CABALA biomass outputs over a rotation length were used to initialise FullCAM (Figure 1) [51] for calculation of C sequestered in living biomass over a rotation length.FullCAM also provided estimates of C in litter, calculated from substrate type and climate [52,53].Based on published values, FullCAM applies C percentages ranging between 48.4% and 53.3% of total biomass to calculate the C percentage of biomass.

Statistical Analyses
The means and standard errors of the final living and dead biomass for the ten model runs per defoliation scenario were calculated.The percentage difference in biomass C (living and dead) between defoliated and undefoliated stands was calculated.Group regression analysis was used to determine differences between current, 2030 and 2070 climates in the relationship between mean annual rainfall and the percentage difference in biomass C between defoliated and undefoliated stands.Regression analysis was also performed to determine whether the relationship between biomass (living and dead) C sequestered under current climate and 2030 or 2070 climate differed significantly from the 1:1 line.Frequency histograms were generated to examine the direction of change in defoliation impact between current and future climates.All analyses were performed using GENSTAT 15 (VSNI International, Hemel Hempstead, UK).

Impact of Defoliation on Biomass C Sequestration under Current Climate
FullCAM estimates of C sequestered in living biomass of E. globulus varied between 4.8 and 13.4 Mg•ha −1 •year −1 under current climate (Figure 2).There was considerable between-site variation, reflecting variable growth rates (Table 1).Estimated C in litter was very low (0.06-0.17 mg•ha −1 •year −1 ) (data not presented).
At only two of the six E. globulus sites was defoliation predicted to reduce biomass C by more than 10% under current climate conditions (Table 3), and of those sites, FullCAM estimated maximum reductions in C sequestration of around 13%.The defoliation levels are those predicted for each site and climate scenario using CLIMEX, whereby L = 10%-30% leaf loss and 20% frequency of years affected; M = 40%-60% leaf loss and 60% frequency of years affected; and S = 60%-80% leaf loss and80% frequency of years affected (see Table 2).Defoliation impacts were modelled in 10% defoliation increments in each severity class and averaged.Forests 2014, 5 1233

Biomass C Responses to Changing Climate in the Absence of Defoliation
Biomass C sequestration (defined as Mg C•ha −1 •year −1 ) in E. globulus plantations was predicted to increase, up to a maximum of 17.2 mg C•ha −1 •year −1 in 2030, and a shift in minimum and maximum values to 7.6 and 17.6 respectively in 2070 (Figure 2).
Under 2030 climate conditions, this increase was greatest compared to current conditions at high productivity sites and where it was assumed that there was no acclimation of photosynthesis to elevated CO 2 (i.e., a CO 2 fertilisation effect) (Figure 2A).Biomass C under 2030 climate was significantly greater than that under current climate, for all climate and CO 2 scenarios.Both the intercept and the slope of the relationship differed significantly from the 1:1 line (p < 0.001), indicating that biomass C would increase proportionally more at higher productivity sites, and decrease slightly at lower productivity sites, in 2030 compared with current climate.
Under 2070 climate and CO 2 conditions, increases in biomass C became more pronounced, with increases compared with current biomass C predicted across the productivity gradient under conditions where no photosynthetic acclimation to elevated CO 2 was assumed (Figure 2B).Under this scenario, greater C sequestration was predicted for the Hadley GCM than the CSIRO GCM.As for 2030, both the intercept and the slope of the relationship between current and future biomass C differed significantly from the 1:1 line (p < 0.001).
Under 2070 conditions, an increase in biomass C also was predicted for lower productivity sites, when photosynthetic acclimation to elevated CO 2 was assumed (Figure 2B), but only under the CSIRO GCM, and the increases were much smaller than those predicted for the scenario assuming no acclimation of photosynthesis.

Impact of Defoliation on Biomass C under Future Climates
The impact of defoliation on biomass C sequestration of E. globulus under future climates varied considerably between sites, assumptions of CO 2 acclimation, and the climate model applied (Table 3).In 2030, when the CSIRO GCM was used and photosynthetic acclimation to elevated CO 2 was assumed, an increase in defoliation impact was predicted in 80% of cases (Figure 3A), while in 20% of cases, defoliation impact was less than under current conditions.Using the Hadley GCM, however, defoliation impact was predicted to decrease in 66% of cases, and of those cases where impact of defoliation decreased, 15% showed a reduction in impact of between 10% and 20%.
If we assumed no photosynthetic acclimation to elevated CO 2 , in the majority of cases the impact of defoliation on C sequestration of E. globulus was predicted to increase (Figure 3B).In only 15% of cases was a decrease in defoliation impact predicted.
When the 2070 CSIRO GCM was used and it was assumed that there was photosynthetic acclimation to elevated CO 2 , the impact of defoliation on C sequestration was predicted to decrease in nearly 40% of cases (Figure 3C).Using the Hadley GCM, this figure was close to 60%.If we assumed no photosynthetic acclimation, defoliation impact was predicted to increase in 70% of cases (Figure 3D).A positive value on the x-axis indicates a reduction in defoliation impact compared to current; a negative value indicates an increase in defoliation impact compared to impact predicted for current condition.

Does Abiotic Stress Change the Impact of Defoliation on Biomass C?
The question of whether abiotic stress influences the impact of defoliation on C sequestration was examined using rainfall as a case study.There was a positive relationship between mean annual rainfall and total C (Figure 4A), although this was not statistically significant.There was a significant negative relationship between mean annual rainfall and the loss of biomass C due to defoliation (R 2 = 0.39, p = 0.008).Hence defoliation had the greatest impact on biomass C of E. globulus at higher rainfall sites (Figure 4B).

Biomass C sequestration in the Absence of Defoliation
In the absence of defoliation, we predicted increases in C sequestration in 2030 and 2070 compared with current climate, but only if no acclimation of photosynthesis to rising atmospheric CO 2 was assumed.At lower productivity sites, the assumption of photosynthetic acclimation resulted in lower predicted C sequestration in 2030.Photosynthesis is inherently limited by atmospheric CO 2 [54], and increased photosynthetic and growth rates are observed when plants are grown at elevated CO 2 [50].Species commonly acclimate to growth at elevated CO 2 , as a result of limitations from resources other than CO 2 that may arise as stands grow and competition develops between individuals [50,55].The degree to which trees acclimate photosynthesis at elevated CO 2 is highly variable [50].Studies in Australia suggest that E. globulus may respond strongly to enhanced atmospheric CO 2 concentrations under conditions of adequate water and nutrient availability [45], but photosynthesis may acclimate under conditions of resource limitation as are commonly observed under field conditions [46].Our results suggest that, if photosynthetic acclimation is observed in the field, then biomass C sequestration in E. globulus plantations is unlikely to increase under the future climate scenarios examined, and may in fact decline slightly at lower productivity sites.

Defoliation Reduced Biomass C Sequestration
Damage to forest ecosystems from pests is widespread but in Australia is generally considered to be of low severity [56].Occasional outbreaks occur but are generally localised in nature.This generally low severity of impact is indicated in the simulations where under current conditions, our simulated defoliations decreased biomass C sequestration by up to 13%.
Our results give clear indications that the impact of defoliation on biomass C sequestration will increase under future climates at some sites, and that there will be considerably more variation in responses than is observed under current conditions.Our results suggest that defoliation has the potential to reduce future biomass C sequestration to current levels or lower, depending on the degree to which photosynthesis acclimates to higher atmospheric CO 2 .Hence the interaction of climate change and defoliation may pose an important risk to C sequestration in the future.
The defoliation scenarios that were included in this study were developed based around our understanding of key periods of vulnerability to defoliation by MLD.In E. globulus this is up to age 4-5 years.However eucalypts are susceptible to defoliation from a range of pests throughout their lifecycle, and severe levels of defoliation have been observed in later-age stands [7].If changing climate results in increased activity of defoliating pests during other stages of plantation development, then our predictions of the impacts of defoliation on biomass C sequestration will be underestimates.
Severe epidemics reduce biomass C via both decreases in forest net primary production and increases in tree mortality [1,3].The modelling approach used in this study focused on changes in forest productivity rather than tree mortality.CABALA cannot currently track tree mortality.E. globulus is susceptible to tree mortality at severe levels of defoliation (e.g., [10]).Hence our predictions of the impact of defoliation on C sequestration may be underestimates under conditions likely to result in tree death.
While we have demonstrated that defoliating pests such as MLD may reduce biomass C sequestration in plantations of E. globulus, other pests, particularly stem borers, are likely to have a much more severe impact on biomass C sequestration because of the extensive tree mortality they can cause.Stem borers are favoured by stress and particularly drought conditions [15,[57][58][59].Our level of understanding of the host's physiological responses to stem borer attack is insufficient for modelling the impacts of stem borers on biomass C at this stage.While leaf chewers that target both juvenile and adult foliage throughout the life of a rotation also may result in serious reductions in productivity, we currently lack the CLIMEX models to allow us to apply our approach to define impacts of these pests on C sequestration.

Influence of Site Factors on Defoliation Impact
Our study suggests that, for E. globulus, impacts of defoliation on biomass C are likely to increase with increasing stand productivity.This is consistent with observations from dry sites in Western Australia, where decreases in stand leaf area were not found to affect stand productivity [55].The result implies that defoliation may improve plant water status under some dry conditions, thereby imparting an advantage to defoliated stands.The increased impact on high productivity sites reflects the marginal gain in production per unit leaf area from these sites being markedly higher than the marginal gain per unit leaf area on lower productivity sites.

Model Assumptions and Uncertainties
Our approach of estimating biomass C sequestration assuming either no or full acclimation to elevated CO 2 is one way of dealing with uncertainty associated with responses to elevated CO 2 , and provides us with an upper and lower boundary of likely responses.In our study the degree of photosynthetic acclimation to elevated CO 2 that was assumed had a large impact on estimates, with differences of up to 25% predicted depending on the assumption.It is likely that, under field conditions, photosynthetic acclimation will be observed to some degree in older stands, although in young, rapidly growing stands with less competition for resources, acclimation may not be observed to the same extent [44].
The assumptions around site-specific climatic conditions under which plantations are growing also had a small impact on predicted biomass C sequestration, with predictions differing between the CSIRO and Hadley GCM's by up to 10%.At most sites the impact of defoliation were predicted to be less under the Hadley than the CSIRO GCM, although irrespective of the GCM used, mean annual rainfall was predicted to be a strong determinant of defoliation impact.
As we move from 2030 to 2070 predictions of biomass C sequestration, there is increasing uncertainty around how well the model predicts at the higher levels of atmospheric CO 2 likely to be experienced.Little is known about how E. globulus will respond to high levels of CO 2 , in terms of degree of acclimation of photosynthesis to elevated CO 2 , feedbacks in the plant related to nutrient and water availability [45,46]; and any direct effects of elevated CO 2 on stomata.
When determining the appropriate scale at which to predict impacts of defoliation on biomass C sequestration, it is important to consider the uncertainties associated with different scales of prediction.Our results highlight the considerable variability in responses to defoliation that might be expected due to between-site variation.This variation is unlikely to be captured if predicting at the regional scale.However there is considerable uncertainty around future changes to rainfall and temperature [60] and this uncertainty escalates when moving from the regional to the local scale.
The modelling tools available for this type of analysis also bring their own scaling problems.CLIMEX and FullCam are both commonly used at a regional or coarser scale, although they can be applied to specific locations which is the approach we used to ensure predictions were at the same scale as CABALA which predicts at a stand scale.

Conclusions
E. globulus plantations were predicted to sequester between 4.8 and 13.4 mg C•ha −1 •year −1 under current conditions.While this may increase in the future, our results suggest this is only likely to occur if photosynthesis does not acclimate to higher atmospheric CO 2 , and thus we cannot assume that biomass C stocks will be enhanced by future climates.In addition, biomass C sequestration may be considerably reduced by defoliation events, with reductions of up to 35% predicted for some sites under future climates.Reductions were generally greatest at higher rainfall sites, although at the lowest rainfall site (Lotons) large reductions also were predicted.It seems likely that, while defoliation has a generally small impact under current conditions, this impact will increase in the future.We conclude that biomass C sequestration is likely to become more vulnerable under future climates.
In this study, we took the approach of trying to develop realistic defoliation scenarios using real sites and linking predicted changes in pest species distribution to plantation productivity.This limits the conclusions that can be drawn, but highlights the importance of understanding variability in responses for improving accuracy of predictions.

Figure 2 .
Figure 2. Relationship between estimates of C sequestration under current climate and (A) 2030 and (B) 2070 climate scenarios for E. globulus stands in southern Australia, for two climate models and assuming either no or full photosynthetic acclimation to elevated CO 2 .Values are mg•ha −1 •year −1 .

Figure 3 .
Figure 3. Changes in defoliation impact for E. globulus, compared to predictions under current climate, for (A) 2030 climate and photosynthetic acclimation; (B) 2030 climate and no photosynthetic acclimation; (C) 2070 climate and photosynthetic acclimation; and (D) 2070 and no photosynthetic acclimation.The frequency of occurrence is shown for both the CSIRO and Hadley GCMs.A positive value on the x-axis indicates a reduction in defoliation impact compared to current; a negative value indicates an increase in defoliation impact compared to impact predicted for current condition.

Figure 4 .
Figure 4. Relationship between (A) total C (standing biomass + litter) and rainfall, and (B) between the difference in total C between defoliated and undefoliated trees and mean annual rainfall for E. globulus.Estimates were made for current, 2030 and 2070 climate conditions using the CSIRO-Mk3 GCM.

Table 1 .
Names and grid references of the E. globulus sites used for the prediction of pest impacts on C sequestration under current and changing climate.Long-term (1961Long-term ( -1991)  )mean annual rainfall (mm) and mean annual temperature (°C) and entire stand volume at age 10 years are also presented.

Table 2 .
[48]anticipated severity and frequency of defoliation by Mycosphaerella leaf disease at six E. globulus sites in southern Australia based on climatic suitability modelled using CLIMEX[48], and the anticipated leaf loss associated with each defoliation severity.Defoliation severity was estimated using the Hadley Mk2 climate model.This information was used to develop defoliation scenarios in CABALA.L = low, M = moderate, and S = severe defoliation.

Table 3 .
•year −1 .Difference (%) in biomass carbon (living and dead) between defoliated and undefoliated Eucalyptus globulus growing at six sites in southern Australia.Estimates are presented for a stand age of 10 years, for the two climate models (CSIRO and Hadley) for current, 2030 and 2070 conditions, and assuming either photosynthetic acclimation to elevated CO 2 , or no acclimation of photosynthesis to elevated CO 2 .