Skip to main content
Download PDF

Insect herbivory in a mature Eucalyptus woodland canopy depends on leaf phenology but not CO2 enrichment

BMC Ecology volume 16, Article number: 47 (2016) Cite this article

Abstract

Background

Climate change factors such as elevated atmospheric carbon dioxide concentrations (e[CO2]) and altered rainfall patterns can alter leaf composition and phenology. This may subsequently impact insect herbivory. In sclerophyllous forests insects have developed strategies, such as preferentially feeding on new leaf growth, to overcome physical or foliar nitrogen constraints, and this may shift under climate change. Few studies of insect herbivory at elevated [CO2] have occurred under field conditions and none on mature evergreen trees in a naturally established forest, yet estimates for leaf area loss due to herbivory are required in order to allow accurate predictions of plant productivity in future climates. Here, we assessed herbivory in the upper canopy of mature Eucalyptus tereticornis trees at the nutrient-limited Eucalyptus free-air CO2 enrichment (EucFACE) experiment during the first 19 months of CO2 enrichment. The assessment of herbivory extended over two consecutive spring—summer periods, with a first survey during four months of the [CO2] ramp-up phase after which full [CO2] operation was maintained, followed by a second survey period from months 13 to 19.

Results

Throughout the first 2 years of EucFACE, young, expanding leaves sustained significantly greater damage from insect herbivory (between 25 and 32 % leaf area loss) compared to old or fully expanded leaves (less than 2 % leaf area loss). This preference of insect herbivores for young expanding leaves combined with discontinuous production of new foliage, which occurred in response to rainfall, resulted in monthly variations in leaf herbivory. In contrast to the significant effects of rainfall-driven leaf phenology, elevated [CO2] had no effect on leaf consumption or preference of insect herbivores for different leaf age classes.

Conclusions

In the studied nutrient-limited natural Eucalyptus woodland, herbivory contributes to a significant loss of young foliage. Leaf phenology is a significant factor that determines the level of herbivory experienced in this evergreen sclerophyllous woodland system, and may therefore also influence the population dynamics of insect herbivores. Furthermore, leaf phenology appears more strongly impacted by rainfall patterns than by e[CO2]. e[CO2] responses of herbivores on mature trees may only become apparent after extensive CO2 fumigation periods.

Background

Climate change and its drivers can have a significant impact on the physiology, abundance and distribution of insect herbivores [1–3]. Elevated CO2 concentrations (e[CO2]) often reduce the growth and survival of insect herbivores as a plant-mediated effect influenced by the decrease in leaf nitrogen concentrations [4–7] and an increase in secondary metabolites, such as phenolic compounds [8–10] generally observed at e[CO2]. Furthermore, some studies suggest modulation of plant hormone signalling and induced plant defence at e[CO2] [9, 10]. Climate change may also alter the timing and amount of precipitation, and this can potentially impact insect abundance and phenology both directly, and indirectly as a consequence of changes in plant phenology and productivity [11–13]. It has been demonstrated that e[CO2] can increase plant net primary production (NPP) [14, 15] and this could potentially benefit insect herbivores as a result of greater resource availability. However, this increase in NPP may be constrained or even reduced in nutrient-limited [16, 17] or water-limited forests [18], or also due to changing herbivory patterns. Furthermore, the measurement of NPP in field experiments may be underestimated if the impacts of herbivory are not measured, in particular the failure of new leaves to expand, due to herbivory on meristems and very young expanding leaves.

Eucalyptus (Myrtaceae) is both an ecologically and economically important tree genus in many parts of the world [19, 20]. Eucalyptus species are often characterised by sclerophyllous leaves and are often associated with low-fertility soils common in Australia [21–23]. As for all plants, the chemical and physical properties of Eucalyptus leaves change with age; young leaves typically have higher nitrogen concentration and moisture content, and reduced toughness compared to older leaves [11, 24, 25]. These factors increase palatability of young foliage to many herbivorous insects, and this can result in enhanced insect performance when feeding on young compared to older leaves [26–28]. Changes in the amount and occurrence of rainfall events may alter the relationship of insect herbivores with leaf phenology, potentially affecting diversity and abundance of insects within these forests. Furthermore, many plants, including Eucalyptus, invest heavily in secondary defence compounds [29, 30], and the production of these secondary compounds may vary throughout leaf development [29, 31, 32]. Herbivore induced plant defence, however, does not appear to occur in Eucalyptus [33] but see [34].

Increased consumption of leaves, or compensatory feeding, is often observed in herbivorous insects as a response to plants grown under e[CO2]—this is to compensate for the dilution of leaf nitrogen [5]. Leaf consumption by herbivorous insects at e[CO2] may result in an additional 17–40 % of leaf damage compared to current levels [7, 35]. Despite the potential for compensatory feeding at e[CO2], the survival of herbivorous insects may be reduced while developmental time is often increased [5, 36–38]. Therefore, e[CO2] may increase leaf damage caused by individual insects to forest trees, however, the abundances of insects in these forests may be reduced due to the negative effects of e[CO2] on insect survival and development. Overall this may result in no net change of leaf damage to trees at e[CO2].

It has previously been demonstrated that the production of insect herbivore excrements (or frass), a crude proxy for insect abundance and herbivory, increased after large rainfall events during the spring and summer at the Eucalyptus free-air CO2 enrichment (EucFACE) experiment [11]. This increase in insect activity coincided with an increase in leaf area index (LAI) at EucFACE [18], suggesting a direct link between the abundance of Eucalyptus-feeding insects and leaf phenology. However no e[CO2] effects on frass deposition or LAI changes within the first 2 years of EucFACE were found, suggesting that insect herbivory and canopy processes may not be impacted by early stages of [CO2] fumigation at EucFACE.

The measurement of frass deposition onto the woodland floor does not reveal which leaf age class experiences most damage from herbivory. Furthermore, e[CO2] may change leaf phenology and thereby resource availability for herbivores. It may also alter the preference of insect herbivores for different leaf stages if the relative palatability of young expanding versus fully expanded (mature) or old leaves changes under e[CO2]. Any change in the consumption of young expanding leaves may therefore affect the recruitment of new leaves in the forest canopy, and place stress on plants. For new leaves, LAI measurement methods may struggle to discriminate between insect removal of leaf area and reductions in NPP. This can result in an incorrect estimate and under-evaluation of NPP of forests, particularly if climate change factors alter the herbivory of new leaf production.

This study investigated the relationship between insect herbivory and leaf phenology of Eucalyptus tereticornis Sm., and the impacts of e[CO2] and rainfall patterns on these processes in a mature, evergreen canopy of this tree species forming a naturally established woodland at the EucFACE experimental site. We hypothesised that rates of insect herbivory would respond to new leaf production which again would vary across time based on rainfall. It has previously been demonstrated that rainfall is the key driver of Eucalyptus leaf phenology, including at the study site [18, 39, 40]. The aims of this study were to: (1) compare the monthly levels of insect leaf herbivory under ambient and e[CO2] conditions within a mature Eucalyptus canopy forming a woodland for two spring and summer periods at which herbivore activity was observed to be highest in the first and second year of EucFACE [11]; (2) provide estimates of leaf damage for different leaf age classes (young, mature, old) during the same two time periods which included the major new leaf production events of E. tereticornis [18], and (3) determine whether specific leaf age classes were preferred by insect herbivores and if this preference was altered under e[CO2] during the first 2 years of EucFACE.

Methods

Study site

This study was conducted at the Eucalyptus free-air CO2 enrichment (EucFACE) experiment located within a native Cumberland Plain woodland remnant [41, 42] in Richmond, NSW, Australia (33°37′S, 150°44′E). The vegetation at EucFACE has been undisturbed for at least 75 years and retains old-growth trees mixed with some re-growth. The vegetation community within the study site is characterised as Cumberland Shale Plains Woodland [41], with mature E. tereticornis as the only canopy forming tree species. The site has an open canopy, approximately 600 trees ha−1 [43], with a low density of forbs and occasional shrubs in the understorey, together with a diverse community of grasses. The site is on a loamy sand soil of the Richmond Formation [44], which is phosphorus-poor and limits tree growth at the site [45]. The average monthly temperature at the site during the time period of this study was 20 °C with an average monthly rainfall of 73 mm (Additional file 1: Fig. S1).

Six large 25 m diameter rings with a height of 28 m above ground, extending above the tree canopy, were constructed amongst the vegetation of the site. Adjacent to each ring stands a high canopy crane with a person basket that allows access to the canopy from above. Each ring also contains a central scaffold tower. Three rings were fumigated diurnally with CO2 enriched air via a proportional-integral-derivative control algorithm [46], while the three remaining rings were control rings, fumigated with ambient air. Beginning in September 2012, the target [CO2] in treatment rings was increased by 30 µmol mol−1 every month until February 2013; thereafter diurnal [CO2] targets within the treatment rings were 150 µmol mol−1 above ambient levels of ~400 µmol mol−1. Rainfall was recorded using automated tipping bucket gauges (Tipping Bucket Rain gauge TB4, Hydrological Services Pty Ltd, Liverpool, NSW, Australia) located 23.5 m above the ground on the central tower in three rings. Data from these sensors were logged every 15 min using CR3000 data loggers (Campbell Scientific, Townsville, Australia).

Leaf herbivory and leaf production measurement

From each of the six EucFACE rings three trees were randomly selected and marked in the first year, and a different set of three trees per ring was selected and marked in the second year. The upper canopy (approximately 17 m above-ground) was accessed using the canopy cranes to establish herbivory observation points on each selected tree. For each year, 14 branches per tree were tagged. We expected that the herbivory measured in the upper canopy was representative for the entire tree canopy because it had previously been demonstrated that Eucalyptus trees display a homogeneous pattern of herbivory throughout the crown [47]. Leaves on each branch were numbered sequentially from the base to the proximal end. A black permanent marker was used to mark the abaxial leaf surface near the petiole, and this has previously been demonstrated not to alter leaf formation or herbivory [48]. New leaves were marked as they emerged behind the shoot tip. Branches and leaves were initially marked in October 2012 and branches were then monitored monthly until February 2013 during the CO2 ramp-up phase (year 1). Three different trees per ring were selected, marked and observed monthly in the second monitoring period from August 2013 to March 2014 (year 2).

The surveyed periods coincided with the majority of chewing insect herbivore activity, as measured by frass deposition to the woodland floor [11] and the growth period of E. tereticornis during the austral spring and summer as measured by changes in LAI [18]. For each of the 2 years, the initial measurements of leaf area in October 2012 and August 2013 were used as a baseline to measure subsequent leaf consumption. For this purpose, approximately 100 leaves across the 14 branches per tree were marked for the monitoring throughout the consecutive months. Branches were selected for ease of access with crane and away from scaffolding to reduce risk of mechanical leaf damage or loss. Branches were then surveyed once each month for a period of four months (year 1) and seven months (year 2) and assessed for leaf damage that can be attributed to insect herbivores due to feeding marks and new leaf emergence. Leaf damage due to insect herbivores was recorded in two different ways: firstly leaves were classified into three age classes (see below) and monthly leaf damage was measured as damage within each of these age classes following the formula below; secondly cumulative leaf damage was calculated for individual leaves throughout their development during the two observation periods in year 1 and year 2 of EucFACE. Thus, cumulative leaf damage refers to the total amount of damage occurring within each monitoring period and not over the life of a particular leaf.

Based on size, colour, shape and texture, leaves were assigned to one of three age classes for each canopy survey point: young (new expanding leaves), mature (fully expanded) and old leaves [48]. Age class-specific herbivory was then calculated as average damage to leaves of each age class throughout the observation periods of each year. For each survey month, scaled digital photographs were taken of leaves that were still attached to branches. For this purpose, leaves were flattened between a scaled white board and a clear non-reflective plastic sheet [49]. Photographs from each month were then compared to photographs taken in the previous month. For each month, existing leaf area (LAe) was quantified using Adobe Photoshop CS5 (Adobe Systems Incorporated, California, USA) by manually tracing the leaf using scaled photographs. Potential leaf area (LAp) i.e. the extent of the leaf area if herbivory had not occurred, was determined by manually drawing and digitally reconstructing the leaf [50]. Skeletonising, mining, and leaf rolling were rarely observed on marked leaves, and therefore disregarded in this study. The recorded damage was exclusively due to removal of leaf area by chewing insects. Monthly increments of leaf consumption were determined with the formula [48]:

$$Lc_{(n + 1)} = \left( {1 {-} \left( {\left( {\frac{{LA_{e} }}{{LA_{p} }}} \right) {-} Lc_{n} } \right)} \right)*100$$

where Lc (n+1) is the proportion (% missing) of leaf area consumed within the observation period of one month; LA e is the actual leaf area recorded for that month; LA p is the potential leaf area if herbivory had not occurred; and Lc n is the proportion of herbivory that had occurred in the previous month.

Total leaf consumption at the end of each monitoring period in years 1 and 2 was determined by the sum of the monthly leaf consumption for each respective leaf age class per branch per tree for the monitored period. Loss of entire leaves due to herbivory was distinguished from leaf loss that may occur as a consequence of leaf senescence and abscission, or due to wind. A leaf that had completely disappeared was considered as lost due to herbivory if it had signs of herbivory in the previous month, or if an entire leaf that was undamaged in the previous month had disappeared except for its petiole. A leaf without signs of herbivory in the previous month was considered lost due to senescence and abscission, or due to wind, if it had disappeared together with its petiole. In this case no value of herbivory was assigned. This approach to assign complete loss of leaves due to herbivory is a conservative measure as it only covers known herbivory. Leaf production was determined by the average number of young new expanding leaves present per branch divided by the total number of leaves per branch and averaged per tree.

Statistical analysis

Linear mixed effects models were constructed using nlme [51] in R [v3.2.2, 52]. The fixed model contained [CO2], month and their interaction. The random model included ring with tree as a nested factor to account for repeated measures. An autocorrelation function was used in order to test for temporal autocorrelation within years, and an autoregressive moving average (ARMA) correlation structure was employed to model dependence among observations of leaf consumption and leaf production across months using a first-order autoregressive structure (AR1) [53]. The number of expanding leaves present was log + 1 transformed to normalise the model-standardised residuals. The relationships of monthly leaf damage with both the average number of young leaves per branch and with rainfall, were modelled using linear mixed effects models and R2 values were obtained using the r.squaredGLMM function in the MuMIn R package [54, 55].

Results

Approximately 3000 E. tereticornis leaves were measured for herbivory in each year across three different age classes. Leaf age classes had a significant effect on leaf consumption both in year 1 (F 2,30 = 245.654, P < 0.001; Fig. 1a) and in year 2 (F 2,30 = 286.435, P < 0.001; Fig. 1b). Young leaves incurred approximately ten times more leaf damage than either mature or old leaves (averages ranged between 25 and 32 % loss in leaf area for young leaves versus less than 2 % loss for mature or old leaves) in both years (Fig. 1). No significant CO2 treatment effect on leaf consumption was observed across leaf age classes (year 1: F 1,4 = 0.399, P = 0.562; year 2: F 1,4 = 0.042, P = 0.848; Fig. 1).

Fig. 1
figure 1

Leaf consumption during three different leaf age classes; young (expanding leaves, dark grey bars), mature (fully expanded, light grey bars), and old leaves (open bars) of mature Eucalyptus tereticornis trees grown at ambient or elevated [CO2] across two time periods, in year 1 (a) and year 2 (b). The figure inset indicates leaf age

Full size image

Significant temporal variation in the amount of leaf consumption during the monitoring periods was observed in both years (year 1: F 3,48 = 10.108, P < 0.001; year 2: F 6,92 = 30.998, P < 0.001; Fig. 2). Monthly leaf consumption peaked in December in the first year and in January in the second year. No significant difference in monthly leaf consumption was observed between CO2 treatments in either year (year 1: F 1,4 = 3.992, P = 0.116; year 2: F 1,4 = 0.028, P = 0.876; Fig. 2). Total cumulative leaf consumption observed during the monitoring periods did not differ between CO2 treatments (year 1: F 1,4 = 6.341, P = 0.066; year 2: F 1,4 = 1.681, P = 0.265; Table 1). The loss of young leaf production between the two years was nearing significance, with less young leaf production being lost in the second year (F 1,4 = , P = 0.051; Table 1).

Fig. 2
figure 2

Monthly leaf consumption experienced by all leaf age classes in year 1 (a) and year 2 (b), and the average number of young expanding leaves observed per branch in year 1 (c) and year 2 (d) on mature Eucalyptus tereticornis trees exposed to ambient (open circles or bars) or elevated (closed circles or bars) [CO2] at the EucFACE site. The figure insets indicate CO2 treatment

Full size image
Table 1 Mean percentage (±SE) of total cumulative leaf damage, young leaf production completely lost to herbivory and young leaves which remained undamaged during the expansion stage on mature E. tereticornis under ambient or elevated [CO2] at the EucFACE site over 2 years
Full size table

The average number of young leaves present per branch in the first year peaked in November (F 3,48 = 21.999, P < 0.001; Fig. 2c) while production of young leaves was highest in January and February of year 2 (F 6,92 = 68.570, P < 0.001; Fig. 2d). No significant differences in young expanding leaf production were observed between CO2 treatments in either year (year 1: F 1,4 = 7.075, P = 0.056; year 2 F 1,4 = 0.114, P = 0.753). Furthermore, e[CO2] did not alter the timing of flush production between year 1 and 2 (F 1,167 = 3.004, P = 0.085). Strong positive correlations were observed between the number of expanding leaves present and leaf consumption (R2 = 0.303, d.f. = 161, P < 0.001; Fig. 3a). Rainfall that occurred two months prior to the monitoring data points was positively correlated with leaf production (R2 = 0.317, d.f. = 161, P < 0.001; Fig. 3b). Temperature was also positively correlated with leaf production (R2 = 0.137, d.f. = 161, P < 0.001; Fig. 3c).

Fig. 3
figure 3

Linear mixed effect model regression of young leaf production per branch and monthly leaf consumption (a), rainfall from two months prior and leaf production per branch (b), and temperature and leaf production per branch (c) observed on Eucalyptus tereticornis trees at the EucFACE site under ambient (open symbols) or elevated (filled symbols) [CO2] in year 1 (circles) and year 2 (triangles)

Full size image

The number of new young leaves produced by trees and lost to herbivory did not differ as a result of CO2 treatment in either year (year 1: F 1,4 = 0.385, P = 0.569; year 2: F 1,4 = 0.148, P = 0.720; Table 1). Overall, approximately 37 % of new young leaves produced by trees escaped any form of chewing damage (Table 1). CO2 treatments did not affect this percentage of young leaves escaping any form of chewing damage in either years (year 1: F 1,4 = 0.043, P = 0.846; year 2: F 1,4 = 1.060, P = 0.361).

Discussion

Herbivory by chewing insects was measured over the first 2 years of the EucFACE experiment. In the first year, e[CO2] was gradually increased from ambient conditions to +150 µmol mol−1, while in the second year, e[CO2] was maintained at 550 µmol mol−1. The two monitoring periods included two major leaf production events, one in each year during spring and summer. Independent of [CO2], consumption of young expanding leaves was very high and decreased to very low levels once leaves were fully expanded. The drastically higher levels of herbivory on young leaves drove monthly variations in overall leaf damage by insect herbivores due to the variation in leaf production. Furthermore, e[CO2] did not affect damage due to herbivory and total leaf consumption. This was in line with our expectation that we would not detect compensatory feeding because of a previous study in which we detected that concentrations of foliar nitrogen and total phenolics at the EucFACE site were not affected by CO2 fumigation over the first 2 years of EucFACE [11]. Such a lack in plant responses may be due to the capacity of mature trees to compensate short term e[CO2] exposure by retrieving nutrient reserves [11].

We found that herbivory by leaf chewing insects increased with the production of new young leaves suggesting that insect feeding is synchronised with the emergence of new young leaves. Insects generally prefer young expanding leaves over mature and old leaves as a result of higher nutrient content and lower physical defences in young leaves, in particular in sclerophyllous Eucalyptus leaves [27, 36, 48, 56]. We have previously demonstrated that foliar nitrogen concentration was higher in flush growth than mature leaves at EucFACE, independent of [CO2] [11]. Herbivory decreased to very low levels in fully expanded leaves, suggesting that leaf flush chewers rather than senescent leaf chewers were the dominant feeding guild during the monitoring periods of this study. Frass production over the same time period [11] displayed a similar pattern of increased deposition around the periods of new leaf production [18]. This highlights the importance of new leaf production not only to insect herbivores but also for insect-mediated nutrient cycling within forests [11].

Herbivory on E. tereticornis in a native woodland observed in this study is within the range observed in other Australian systems. For example, in an Australian rainforest insect herbivore damage on young expanding leaves ranged between 10 and 30 % depending on tree species, and this fell to less than 5 % once leaves matured [48]. Similarly, Moles and Westoby [57] reported that leaf damage to expanding leaves from 51 woody dicotyledonous species in a coastal dry sclerophyll forest ranged between 0 and 51 %. Irrespective of leaf age, reported total leaf herbivore damage in Australian forests ranges between 5 and 44 % [23, 47, 58, 59]. Overall, at EucFACE we observed total herbivore consumption of between 9 and 17 % of leaf area, which is in the lower range of previous studies conducted in Australian forests. No leaf mining or leaf rolling was observed; therefore the damage measured in this study was exclusively due to leaf chewing insects. This is often the most common type of leaf damage observed in forests, when compared to other insect feeding types [59–61]. The feeding guild of leaf chewers is also the most likely impacted by e[CO2] [7]. A previous study of herbivory on E. tereticornis at the study site also identified leaf chewers as the dominant feeding guild [62].

The responses to eCO2 of two key leaf chewing insect herbivores found at the EucFACE site had previously been tested on young E. tereticornis trees in greenhouse experiments. The cup moth Doratifera quadriguttata (Lepidoptera: Limacodidae) and the leaf beetle Paropsis atomaria (Coleoptera: Chrysomelidae), both experienced negative effects of e[CO2] such as increased mortality, extended developmental times and reduced pupal weights, while both species also displayed signs of compensatory feeding [36–38]. However, we did not observe any differences in the amount of leaf damage occurring to mature E. tereticornis trees at EucFACE under CO2 enrichment. In other forests undergoing CO2 fumigation e[CO2] has been observed to either decrease [63, 64], cause no change [65], or increase [66, 67] the amount of leaf damage. These differences in herbivory responses to e[CO2] across a variety of forest types may be due to the complexity in interactions between biotic and abiotic factors which impact insect herbivores and therefore herbivory.

The different timing of the major flush production events between the 2 years at EucFACE may have affected the composition of herbivore communities within the woodland. Shifts in the timing of flush growth may have detrimental effects on herbivores that depend on flush with the potential outcome of an altered insect herbivore community structure. Although we did not directly assess insect populations our data indirectly suggests that such potential changes in chewing insect herbivore populations may have occurred. It appeared that less of the newly produced leaves were consumed in the second year, potentially due to the altered timing of flush growth. This may also indicate that there was a shift away from flush leaf chewers. In another study that assessed deposition of lerp (small covers produced by plant sap-feeding psyllids) to the woodland floor, a rapid outbreak succession by two psyllid species, Glycaspis sp. and Cardiaspina fiscella was detected on E. tereticornis after March 2014 [68]. This suggests that a rapid change in the overall canopy insect community composition must have occurred directly after the leaf phenology and herbivory surveys presented here. Our previous study on frass deposition by leaf chewing insects also demonstrated that less frass was produced in 2013/2014 (year 2) than in 2012/2013 (year 1) [11].

Our leaf herbivory study only focussed on leaf-chewing insects, and sap-feeding insects were not considered here. This focus on the leaf-chewing feeding guild may underestimate the impacts of herbivory occurring in the canopy and the true level of canopy biomass lost to insect herbivores in our study site. Significant numbers of plant sap-feeding psyllids were detected at EucFACE after March 2014 [68]. Some psyllid species can be significant leaf defoliators of Eucalyptus. For example, Cardiaspina sp. has caused area-wide defoliation on Eucalyptus moluccana but not any other Eucalyptus species in the Cumberland Plain Woodlands since 2009 [69]. Similarly, C. fiscella caused significant defoliation of the EucFACE site at the end of 2014 [68], after the end of the study presented here. This defoliation occurred because Cardiaspina can induce leaf senescence and defoliation [34]. However, throughout the surveying period of our herbivory study at EucFACE (from October 2012 until March 2014) we did not observe any high abundance of leaf defoliating psyllids [68]. After the ramp-up phase (from March 2013 to December 2014) the abundance of lerps deposited to the woodland floor was reduced for three psyllid species at e[CO2], and signs for compensatory feeding were detected in one of these three species, Glycaspis sp. [68]. These differences were recorded despite the absence of any measureable differences in N concentration in leaves [11].

The loss of young expanding leaves, as found in our study, may be more detrimental to plants than the loss of mature leaves, because the energy and resources invested in the production of new leaves have yet to be recovered [48]. This may impact forest growth and the estimation of carbon storage potential within forests and highlights the need to consider the role of herbivorous insects in ecosystem functioning [66]. Removal of leaf material by insect herbivores is often not accounted for in models of plant productivity as common techniques used to measure productivity, such as LAI [70, 71], fail to account for the loss of newly produced plant material to insect herbivory. This can complicate estimates of NPP and may underestimate true productivity of forests and utilisation of forest resources. Furthermore, plant NPP often increases at e[CO2] as a result of a carbon fertilisation effect [72, 73]. However, in the early stages of CO2 fumigation at the EucFACE site we did not observe an increase in leaf production. This may indicate that mature E. tereticornis trees within the site are limited by water and nutrient availability [18, 45], and increased photosynthesis due to e[CO2] may also result in increased respiration in the ecosystem [43] rather than the production of more leaves.

Conclusions

Predictions about damage inflicted by insect herbivores at e[CO2] are difficult to make under field conditions owing to the complex interactions between plants, insect herbivores and their antagonists. This uncertainty in insect herbivore responses hinders the ability to accurately determine model-specifications of their response to climate change [74]. Contrary to our original expectations, CO2 fumigation at the EucFACE site for the first 2 years did not affect total leaf consumption by herbivores or their leaf age preference, in part likely due to the lack of e[CO2] responses in concentrations of foliar nitrogen and total phenolics in insect frass [11]. However, it is clear that new and expanding leaves were heavily damaged while fully expanded leaves were not. Damage on young foliage is often not accounted for in estimates of forest productivity, yet this can amount to substantial underestimates of true forest productivity. Rainfall-mediated production of new leaves is an important regulator of insect herbivory in sclerophyllous forests due to the physical barriers to consumption present in mature leaves [75], and this will require further attention in climate change studies. Shifts in rainfall patterns, a potential outcome of climate change [76], can have significant effects on insect community composition, herbivory and insect frass deposition. This may have detrimental or positive outcomes for ecosystems and humans, for example by stimulating pest populations, or regulating ecosystem services that insects provide in forests and managed plantations. Understanding how rainfall may interact with e[CO2] in altering insect herbivory and herbivore abundances is important in predicting the impacts of climate change variables on insect herbivore population dynamics in forests and plantations.

Abbreviations

e[CO2]:

elevated carbon dioxide concentrations

EucFACE:

Eucalyptus free-air CO2 enrichment

LAe :

existing leaf area

LAI:

leaf area index

LAp :

potential leaf area

Lc n :

proportion of herbivory that had occurred in the previous month

NPP:

net primary production

References

  1. Cornelissen T. Climate change and its effects on terrestrial insects and herbivory patterns. Neotrop Entomol. 2011;40(2):155–63.

    Article  CAS  PubMed  Google Scholar 

  2. Tylianakis JM, Didham RK, Bascompte J, Wardle DA. Global change and species interactions in terrestrial ecosystems. Ecol Lett. 2008;11(12):1351–63.

    Article  PubMed  Google Scholar 

  3. Buckley LB, Kingsolver JG. The demographic impacts of shifts in climate means and extremes on alpine butterflies. Funct Ecol. 2012;26(4):969–77.

    Article  Google Scholar 

  4. Ainsworth EA, Rogers A. The response of photosynthesis and stomatal conductance to rising [CO2]: mechanisms and environmental interactions. Plant Cell Environ. 2007;30(3):258–70.

    Article  CAS  PubMed  Google Scholar 

  5. Robinson EA, Ryan GD, Newman JA. A meta-analytical review of the effects of elevated CO2 on plant–arthropod interactions highlights the importance of interacting environmental and biological variables. New Phytol. 2012;194(2):321–36.

    Article  CAS  PubMed  Google Scholar 

  6. Hunter MD. Effects of elevated atmospheric carbon dioxide on insect–plant interactions. Agric For Entomol. 2001;3(3):153–9.

    Article  Google Scholar 

  7. Stiling P, Cornelissen T. How does elevated carbon dioxide (CO2) affect plant–herbivore interactions? A field experiment and meta-analysis of CO2-mediated changes on plant chemistry and herbivore performance. Glob Chang Biol. 2007;13(9):1823–42.

    Article  Google Scholar 

  8. Ryan GD, Rasmussen S, Newman JA. Global atmospheric change and trophic interactions: are there any general responses? In: Baluska F, Ninkovic V, editors. Plant communication from an ecological perspective. Berlin: Springer-Verlag; 2010. p. 179–214.

    Chapter  Google Scholar 

  9. Zavala JA, Casteel CL, DeLucia EH, Berenbaum MR. Anthropogenic increase in carbon dioxide compromises plant defense against invasive insects. Proc Natl Acad Sci USA. 2008;105(13):5129–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. DeLucia EH, Nabity PD, Zavala JA, Berenbaum MR. Climate change: resetting plant-insect interactions. Plant Physiol. 2012;160(4):1677–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Gherlenda AN, Crous KY, Moore BD, Haigh AM, Johnson SN, Riegler M. Precipitation, not CO2 enrichment, drives insect herbivore frass deposition and subsequent nutrient dynamics in a mature Eucalyptus woodland. Plant Soil. 2016;399:29–39.

    Article  CAS  Google Scholar 

  12. Whitfeld TJS, Novotny V, Miller SE, Hrcek J, Klimes P, Weiblen GD. Predicting tropical insect herbivore abundance from host plant traits and phylogeny. Ecology. 2012;93(sp8):S211–22.

    Article  Google Scholar 

  13. Lee MA, Manning P, Walker CS, Power SA. Plant and arthropod community sensitivity to rainfall manipulation but not nitrogen enrichment in a successional grassland ecosystem. Oecologia. 2014;176(4):1173–85.

    Article  PubMed  Google Scholar 

  14. Norby RJ, DeLucia EH, Gielen B, Calfapietra C, Giardina CP, King JS, Ledford J, McCarthy HR, Moore DJP, Ceulemans R. Forest response to elevated CO2 is conserved across a broad range of productivity. Proc Natl Acad Sci USA. 2005;102(50):18052–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Zak DR, Pregitzer KS, Kubiske ME, Burton AJ. Forest productivity under elevated CO2 and O3: positive feedbacks to soil N cycling sustain decade-long net primary productivity enhancement by CO2. Ecol Lett. 2011;14(12):1220–6.

    Article  PubMed  Google Scholar 

  16. Norby RJ, Warren JM, Iversen CM, Medlyn BE, McMurtrie RE. CO2 enhancement of forest productivity constrained by limited nitrogen availability. Proc Natl Acad Sci USA. 2010;107(45):19368–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Johnson DW. Progressive N limitation in forests: review and implications for long-term responses to elevated CO2. Ecology. 2006;87(1):64–75.

    Article  PubMed  Google Scholar 

  18. Duursma RA, Gimeno TE, Boer MM, Crous KY, Tjoelker MG, Ellsworth DS. Canopy leaf area of a mature evergreen Eucalyptus woodland does not respond to elevated atmospheric [CO2] but tracks water availability. Glob Chang Biol. 2016;22:1666–76.

    Article  PubMed  Google Scholar 

  19. Paine TD, Steinbauer MJ, Lawson SA. Native and exotic pests of Eucalyptus: a worldwide perspective. Annu Rev Entomol. 2011;56:181–201.

    Article  CAS  PubMed  Google Scholar 

  20. FAO. Global forest resources assessment 2010. Rome: Food and Agriculture Organization of the United Nations; 2010.

    Google Scholar 

  21. Ohmart C, Edwards P. Insect herbivory on Eucalyptus. Annu Rev Entomol. 1991;36(1):637–57.

    Article  Google Scholar 

  22. Keith DA. Ocean shores to desert dunes: the native vegetation of New South Wales and the ACT. Hurstville: New South Wales Department of Environment and Conservation; 2004.

    Google Scholar 

  23. Landsberg J, Ohmart C. Levels of insect defoliation in forests: patterns and concepts. Trends Ecol Evol. 1989;4(4):96–100.

    Article  Google Scholar 

  24. Loney P, McArthur C, Sanson G, Davies N, Close D, Jordan G. How do soil nutrients affect within-plant patterns of herbivory in seedlings of Eucalyptus nitens? Oecologia. 2006;150(3):409–20.

    Article  PubMed  Google Scholar 

  25. Close DC, McArthur C, Hagerman AE, Fitzgerald H. Differential distribution of leaf chemistry in eucalypt seedlings due to variation in whole-plant nutrient availability. Phytochemistry. 2005;66(2):215–21.

    Article  CAS  PubMed  Google Scholar 

  26. Rapley LP, Allen GR, Potts BM. Genetic variation in Eucalyptus globulus in relation to susceptibility from attack by the southern eucalypt leaf beetle. Chrysophtharta agricola. Aust J Bot. 2004;52(6):747–56.

    Article  Google Scholar 

  27. Nahrung HF, Duffy MP, Lawson SA, Clarke AR. Natural enemies of Paropsis atomaria Olivier (Coleoptera: Chrysomelidae) in south-eastern Queensland eucalypt plantations. Aust J Ecol. 2008;47(3):188–94.

    Google Scholar 

  28. Boege K, Marquis RJ. Facing herbivory as you grow up: the ontogeny of resistance in plants. Trends Ecol Evol. 2005;20(8):441–8.

    Article  PubMed  Google Scholar 

  29. McKiernan AB, Hovenden MJ, Brodribb TJ, Potts BM, Davies NW, O’Reilly-Wapstra JM. Effect of limited water availability on foliar plant secondary metabolites of two Eucalyptus species. Environ Exp Bot. 2014;105:55–64.

    Article  CAS  Google Scholar 

  30. Moore BD, Wallis IR, Wood JT, Foley WJ. Foliar nutrition, site quality, and temperature influence foliar chemistry of tallowwood (Eucalyptus microcorys). Ecol Monogr. 2004;74(4):553–68.

    Article  Google Scholar 

  31. Silvestre AJD, Cavaleiro JAS, Delmond B, Filliatre C, Bourgeois G. Analysis of the variation of the essential oil composition of Eucalyptus globulus Labill. from Portugal using multivariate statistical analysis. Ind Crops Prod. 1997;6(1):27–33.

    Article  CAS  Google Scholar 

  32. Goodger JQD, Choo TYS, Woodrow IE. Ontogenetic and temporal trajectories of chemical defence in a cyanogenic eucalypt. Oecologia. 2007;153(4):799–808.

    Article  PubMed  Google Scholar 

  33. Henery ML, Wallis IR, Stone C, Foley WJ. Methyl jasmonate does not induce changes in Eucalyptus grandis leaves that alter the effect of constitutive defences on larvae of a specialist herbivore. Oecologia. 2008;156(4):847–59.

    Article  CAS  PubMed  Google Scholar 

  34. Steinbauer MJ, Burns AE, Hall A, Riegler M, Taylor GS. Nutritional enhancement of leaves by a psyllid through senescence-like processes: insect manipulation or plant defence? Oecologia. 2014;176(4):1061–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Coley PD. Possible effects of climate change on plant/herbivore interactions in moist tropical forests. Clim Change. 1998;39(2–3):455–72.

    Article  Google Scholar 

  36. Gherlenda AN, Haigh AM, Moore BD, Johnson SN, Riegler M. Responses of leaf beetle larvae to elevated [CO2] and temperature depend on Eucalyptus species. Oecologia. 2015;177(2):607–17.

    Article  PubMed  Google Scholar 

  37. Murray TJ, Ellsworth DS, Tissue DT, Riegler M. Interactive direct and plant-mediated effects of elevated atmospheric [CO2] and temperature on a eucalypt-feeding insect herbivore. Glob Chang Biol. 2013;19:1407–16.

    Article  CAS  PubMed  Google Scholar 

  38. Murray TJ, Tissue DT, Ellsworth DS, Riegler M. Interactive effects of pre-industrial, current and future [CO2] and temperature on an insect herbivore of Eucalyptus. Oecologia. 2013;171(4):1025–35.

    Article  CAS  PubMed  Google Scholar 

  39. Myers BA, Williams RJ, Fordyce I, Duff GA, Eamus D. Does irrigation affect leaf phenology in deciduous and evergreen trees of the savannas of northern Australia? Aust J Ecol. 1998;23(4):329–39.

    Article  Google Scholar 

  40. Pook EW, Gill AM, Moore PHR. Long-term variation of litter fall, canopy leaf area and flowering in a Eucalyptus maculata forest on the south coast of New South Wales. Aust J Bot. 1997;45(5):737–55.

    Article  Google Scholar 

  41. Tozer M. The native vegetation of the Cumberland Plain, western Sydney: systematic classification and field identification of communities. Cunninghamia. 2003;8(1):1–75.

    Google Scholar 

  42. Benson DH. The natural vegetation of the Penrith 1:100 000 map sheet. Cunninghamia. 1992;2(4):541–96.

    Google Scholar 

  43. Drake JE, Macdonald CA, Tjoelker MG, Crous KY, Gimeno TE, Singh BK, Reich PB, Anderson IC, Ellsworth DS. Short-term carbon cycling responses of a mature eucalypt woodland to gradual stepwise enrichment of atmospheric CO2 concentration. Glob Chang Biol. 2015;22:380–90.

    Article  PubMed  Google Scholar 

  44. Bannerman SM, Hazelton PA. Soil landscapes of the Penrith 1:100,000 sheet map and report. Sydney: Soil Conservation Service of NSW; 1990.

    Google Scholar 

  45. Crous KY, Ósvaldsson A, Ellsworth DS. Is phosphorus limiting in a mature Eucalyptus woodland? Phosphorus fertilisation stimulates stem growth. Plant Soil. 2015;391:293–305.

    Article  CAS  Google Scholar 

  46. Lewin KF, Nagy J, Nettles WR, Cooley DM, Rogers A. Comparison of gas use efficiency and treatment uniformity in a forest ecosystem exposed to elevated [CO2] using pure and prediluted free-air CO2 enrichment technology. Glob Chang Biol. 2009;15(2):388–95.

    Article  Google Scholar 

  47. Lowman MD. Herbivory in Australian forests-a comparison of dry sclerophyll and rain forest canopies. Proc Linn Soc NSW. 1995;115:77–87.

    Google Scholar 

  48. Lowman MD. Temporal and spatial variability in insect grazing of the canopies of five Australian rainforest tree species. Aust J Ecol. 1985;10(1):7–24.

    Article  Google Scholar 

  49. Paul GS, Montagnini F, Berlyn GP, Craven DJ, Van Breugel M, Hall JS. Foliar herbivory and leaf traits of five native tree species in a young plantation of Central Panama. New For. 2012;43(1):69–87.

    Article  Google Scholar 

  50. Unsicker SB, Mody K. Influence of tree species and compass bearing on insect folivory of nine common tree species in the West African savanna. J Trop Ecol. 2005;21(2):227–31.

    Article  Google Scholar 

  51. Pinheiro J, Bates D, DebRoy S, Sarkar D. nlme: Linear and nonlinear mixed effects models. R Package Version 3.1-122. 2015. https://cran.r-project.org/package=nlme. Accessed 19 Aug 2015.

  52. R Development Core Team. R. A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing. 2015. https://www.R-project.org. Accessed 3 July 2015.

  53. Bader MKF, Leuzinger S, Keel SG, Siegwolf RTW, Hagedorn F, Schleppi P, Körner C. Central European hardwood trees in a high-CO2 future: synthesis of an 8 year forest canopy CO2 enrichment project. J Ecol. 2013;101(6):1509–19.

    Article  CAS  Google Scholar 

  54. Johnson PCD. Extension of Nakagawa and Schielzeth’s R2 GLMM to random slopes models. Methods Ecol Evol. 2014;5(9):944–6.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Barton K. MuMIn: Multi-Model Inference. R package version 1.15.1. 2015. https://cran.r-project/package=MuMIn. Accessed 3 July 2015.

  56. Nahrung HF, Dunstan PK, Allen GR. Larval gregariousness and neonate establishment of the eucalypt-feeding beetle Chrysophtharta agricola (Coleoptera: Chrysomelidae: Paropsini). Oikos. 2001;94(2):358–64.

    Article  Google Scholar 

  57. Moles AT, Westoby M. Do small leaves expand faster than large leaves, and do shorter expansion times reduce herbivore damage? Oikos. 2000;90(3):517–24.

    Article  Google Scholar 

  58. Fox LR, Morrow PA. Estimates of damage by herbivorous insects on Eucalyptus trees. Aust J Ecol. 1983;8(2):139–47.

    Article  Google Scholar 

  59. Nooten SS, Hughes L. Patterns of insect herbivory on four Australian understory plant species. Aust J Ecol. 2013;52(4):309–14.

    Google Scholar 

  60. Landsberg J, Gillieson DS. Regional and local variation in insect herbivory, vegetation and soils of eucalypt associations in contrasted landscape positions along a climatic gradient. Aust J Ecol. 1995;20(2):299–315.

    Article  Google Scholar 

  61. Adams JM, Zhang Y. Is there more insect folivory in warmer temperate climates? A latitudinal comparison of insect folivory in eastern North America. J Ecol. 2009;97(5):933–40.

    Article  Google Scholar 

  62. Lopaticki G. The response of eucalypt arthropod communities to varying resource availability and habitat complexity BSc Honours thesis. Penrith: University of Western Sydney; 2010.

  63. Knepp RG, Hamilton JG, Mohan JE, Zangerl AR, Berenbaum MR, DeLucia EH. Elevated CO2 reduces leaf damage by insect herbivores in a forest community. New Phytol. 2005;167(1):207–18.

    Article  CAS  PubMed  Google Scholar 

  64. Stiling P, Cattell M, Moon DC, Rossi A, Hungate BA, Hymus G, Drake B. Elevated atmospheric CO2 lowers herbivore abundance, but increases leaf abscission rates. Glob Chang Biol. 2002;8(7):658–67.

    Article  Google Scholar 

  65. Hamilton JG, Zangerl AR, Berenbaum MR, Pippen J, Aldea M, DeLucia EH. Insect herbivory in an intact forest understory under experimental CO2 enrichment. Oecologia. 2004;138(4):566–73.

    Article  PubMed  Google Scholar 

  66. Couture JJ, Meehan TD, Kruger EL, Lindroth RL. Insect herbivory alters impact of atmospheric change on northern temperate forests. Nature Plants. 2015;1(3):150–60.

    Article  Google Scholar 

  67. Meehan TD, Couture JJ, Bennett AE, Lindroth RL. Herbivore-mediated material fluxes in a northern deciduous forest under elevated carbon dioxide and ozone concentrations. New Phytol. 2014;204(2):397–407.

    Article  CAS  PubMed  Google Scholar 

  68. Gherlenda AN, Esveld JL, Hall AAG, Duursma RA, Riegler M. Boom and bust: rapid feedback responses between insect outbreak dynamics and canopy leaf area impacted by rainfall and CO2. Glob Chang Biol. 2016. doi:10.1111/gcb.13334.

    PubMed  Google Scholar 

  69. Hall AAG, Gherlenda AN, Hasegawa S, Johnson SN, Cook JM, Riegler M. Anatomy of an outbreak: the biology and population dynamics of a Cardiaspina psyllid species in an endangered woodland ecosystem. Agric For Entomol. 2015;17:292–301.

    Article  Google Scholar 

  70. Reich PB. Key canopy traits drive forest productivity. Proc Biol Sci. 2012;279:2128–34.

    Article  PubMed  PubMed Central  Google Scholar 

  71. Walker AP, Hanson PJ, De Kauwe MG, Medlyn BE, Zaehle S, Asao S, Dietze M, Hickler T, Huntingford C, Iversen CM. Comprehensive ecosystem model-data synthesis using multiple data sets at two temperate forest free-air CO2 enrichment experiments: model performance at ambient CO2 concentration. J Geophys Res: Biogeo. 2014;119(5):937–64.

    Article  CAS  Google Scholar 

  72. Ainsworth EA, Long SP. What have we learned from 15 years of free-air CO2 enrichment (FACE)? A meta-analytic review of the responses of photosynthesis, canopy properties and plant production to rising CO2. New Phytol. 2005;165(2):351–72.

    Article  PubMed  Google Scholar 

  73. De Kauwe MG, Medlyn BE, Zaehle S, Walker AP, Dietze MC, Wang YP, Luo Y, Jain AK, El-Masri B, Hickler T. Where does the carbon go? A model-data intercomparison of vegetation carbon allocation and turnover processes at two temperate forest free-air CO2 enrichment sites. New Phytol. 2014;203(3):883–99.

    Article  PubMed  PubMed Central  Google Scholar 

  74. Boulanger Y, Gray DR, Cooke BJ, De Grandpré L. Model-specification uncertainty in future forest pest outbreak. Glob Chang Biol. 2016;22(4):1595–607.

    Article  PubMed  Google Scholar 

  75. Steinbauer MJ. Specific leaf weight as an indicator of juvenile leaf toughness in Tasmanian bluegum (Eucalyptus globulus ssp. globulus): implications for insect defoliation. Austral For. 2001;64(1):32–7.

    Article  Google Scholar 

  76. IPCC (ed.). Climate change 2013: The physical science basis. Contribution of working group I to the fifth assessment report of the Intergovernmental Panel on Climate Change. Cambridge: Cambridge University Press; 2013.

  77. Gherlenda AN, Moore BD, Haigh AM, Johnson SN, Riegler M. Data from: insect herbivory in a mature Eucalyptus woodland canopy depends on leaf phenology but not CO2 enrichment. BMC Ecol. 2016. doi:10.4225/35/57e49ec6dd3eb.

    PubMed  PubMed Central  Google Scholar 

Download references

Authors’ contributions

ANG and MR conceived and designed the experiments with input by BDM, AMH and SNJ. ANG performed the experiments. ANG analysed the data. ANG and MR primarily wrote the manuscript; other authors helped with drafting the text and interpreting the data. All authors have read and approved the final manuscript.

Acknowledgements

We thank Vinod Kumar and Steven Wohl for technical assistance at the EucFACE site.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets generated and analysed for this study are available in the HIEv repository: http://doi.org/10.4225/35/57e49ec6dd3eb [77].

Funding

EucFACE is supported by the Australian Commonwealth Government in collaboration with Western Sydney University. EucFACE was built as an initiative of the Australian Government as part of the Nation-building Economic Stimulus Package. This research was also supported by an Australian Postgraduate Award to ANG and DP1095972 of the Australian Research Council. The funding body had no role in the design of the study, collection, analysis, interpretation of the data and in writing of the manuscript.

Author information

Authors and Affiliations

  1. Hawkesbury Institute for the Environment, Western Sydney University, Locked Bag 1797, Penrith, NSW, 2751, Australia

    Andrew N. Gherlenda, Ben D. Moore, Scott N. Johnson & Markus Riegler

  2. School of Science and Health, Western Sydney University, Locked Bag 1797, Penrith, NSW, 2751, Australia

    Anthony M. Haigh

Authors
  1. Andrew N. Gherlenda
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Ben D. Moore
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Anthony M. Haigh
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Scott N. Johnson
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Markus Riegler
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding authors

Correspondence to Andrew N. Gherlenda or Markus Riegler.

Additional file

12898_2016_102_MOESM1_ESM.docx

Additional file 1: Fig. S1. Total monthly rainfall (bars) and average monthly temperature (circles) at the EucFACE site.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gherlenda, A.N., Moore, B.D., Haigh, A.M. et al. Insect herbivory in a mature Eucalyptus woodland canopy depends on leaf phenology but not CO2 enrichment. BMC Ecol 16, 47 (2016). https://doi.org/10.1186/s12898-016-0102-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12898-016-0102-z

Keywords

BMC Ecology

ISSN: 1472-6785

Contact us