Effects of defoliation by the pine processionary moth Thaumetopoea pityocampa on biomass growth of young stands of Pinus pinaster in northern Portugal
iForest - Biogeosciences and Forestry, Volume 3, Issue 6, Pages 159-162 (2010)
doi: https://doi.org/10.3832/ifor0553-003
Published: Nov 15, 2010 - Copyright © 2010 SISEF
Research Articles
Abstract
Biomass growth models for 13-year-old maritime pine tree stands (Pinus pinaster Ait.) in the north-eastern Portugal were developed and used to analyse the effects of the defoliation by the pine processionary moth, Thaumetopoea pityocampa (Den. & Schiff.) on biomass increment. For the adjustment of the models, 30 individual pine trees were destructively sampled and non-linear models were tested, using the diameter at 10 centimetre height (d0.10), the total height (h), both variables (d0.10+h) and d0.102h as preditors of biomass growth. The results showed that the best predictor was d0.10+h. Application of models to analyse tree biomass after attack by the pine processionary moth showed that the decrease of biomass increment was proportional to the severity of the insect attack, with average values of losses in biomass increment ranging from 37% to 73%, depending on defoliation intensity.
Keywords
Prediction model, Thaumetopoea pityocampa, Biomass increment, Pinus pinaster
Introduction
The most recent National Forest Inventory ([1]), shows that forest represents 37.7% of the land use in Portugal, with P. pinaster covering 25.0% of that area. In northern Portugal, where this study was carried out, forest represents 30.7% of the land use and P. pinaster stands cover 40.2% of this area. Young stands, which represent 30% of total pine stands, are extremely important nationally, not only because they represent the future of pine stands, but also because they have a high ecological and economical importance. According to Oliveira ([20]), shrub biomass should be removed 3 times before the 10-th years after plantation and the first thinning should be applied before the 15-th year, removing between 20% and 40% of the aboveground biomass. Based on results reported by Lopes ([14]), young pine stands can produce an average of 15.8 ton ha-1 year-1 of arboreal biomass and 6.32 ton ha-1 year-1 of shrub biomass. If 20% of the main stand is removed after 15 years, a yield of 4.42 ton ha-1 year-1 is expected, representing high potential for biomass production. This biomass can be sold at up to 30€ per ton ([19]), a significant source of income for forest owners. Pinus biomass from adult stands is mainly used for furniture production, pulp production and biomass for energy among other possible applications.
The pine processionary Moth, Thaumetopoea pityocampa (Den. & Schiff. - Lep., Thaumetopoeidae) is one of the most destructive insects of Pinus and Cedrus in the Middle East, North Africa and many southern European countries, including Portugal. The urticant hairs of the late instar larvae provokes serious reactions in humans and other mammals ([12], [21]) but it is also responsible for significant economic damage due to severe defoliation ([3], [6], [10]). Defoliation removes both photosynthetic material and sites where chemicals such as growth hormones are produced, affecting many vital functions ([5]). It is well known that in adult trees defoliation, though repeated over consecutive years, seldom causes death ([22], [9]), but increases susceptibility to sanitary problems such as pine weevils and bark beetle attack ([10], [17]). In spite of the capability of defoliated trees to refoliate and survive, the effects of defoliation are very significant ([7]) with losses in volume, radial growth ([5], [10], [13]) and biomass production ([17]). [10] reported a significant decrease (average 21% over four years) on annual diameter increment of Pinus brutia in Turkey while Carus ([5]) identified growth reductions on radial height and volume on P. brutia after an outbreak of T. pityocampa. Cadahia & Insua ([4]) identified a decrease in wood volume increment of 14-33% in young P. radiata as well as losses in tree volume, while Bouchon & Toth ([2]) reported that T. pityocampa attack was responsible for an about 45% volume decrease over a 50 year period.
The effect of insect action on trees can be studied using predictive models ([8]) which allow estimates of dendrometric variables (e.g., height, volume, and biomass) and can be used for inventory techniques in production studies ([11], [23]). However, the use of allometric models to estimate the impact of insects on forest dynamics is virtually unknown.
The aims of this paper are: (1) to present growth biomass models for different above-ground components (leaves, stem and total biomass) of young pine trees; and use them (2) to quantify the effects of defoliation by the pine processionary moth on P. pinaster biomass increment.
Material and methods
Study site
The study was conducted in the Natural Park of Montesinho situated in Trás-os-Montes, a mountainous region of north-east Portugal. The Park (area 74 800 ha) is situated in the “Terra Fria Transmontana” climate zone, characterized by hot and dry summers, cold winters (annual mean temperature around 11 ºC) and precipitation falling mainly during the autumn (annual mean precipitation about 900 mm). The area comprises pure stands of several pine species; P. pinaster Ait. and P. nigra Arn. are the dominant species and P. sylvestris L. and P. strobus L. the secondary species.
The experimental plots were set up in 13-year-old P. pinaster plantation with a density of 650 trees per hectare.
Biomass models adjustment
Adjustment of the biomass models was performed by destructively sampling 30 pine trees not suffering from defoliation by the pine processionary moth. Trees were cut near the soil surface and the total tree height (h) and live crown were measured. Stem diameters were measured 10 cm from the base of the tree (d0.10), diameter at breast height (dbh) and then at intervals of 10%, beginning at 10% of the total height (h) and ending at 90% of the height. Bark thickness was measured at each point. The total stem was measured and weighed.
Each branch of the crown was separated from the stem and the fresh mass of each component measured using scales with one kilogramme precision. Leaves, logs and female cones were separated for each branch to measure the relative contribution of each component to total branch weight. The total volume (v) was determined using the diameters that had been measured across the stem.
Samples of each tree component were collected in order to obtain dry density later in the lab.
Pearson’s correlation was calculated between total biomass (including stems, leaves, branches and female cones biomass) and partial biomass (leaves biomass and stems biomass) and measurements of trees size (d0.10, dbh, h and v), in order to identify the most appropriate predictor variables. Adjustment of total and partial biomass prediction models was carried out by cross validation. Sampled trees were were randomly assigned to two groups: some 24 trees were selected for the adjustment phase and 6 trees were selected for the validation phase. Based on the methodology used by Mikšys et al. ([18]) tree biomass components and tree parameters were evaluated and equations for tree biomass evaluation were derived. Several non-linear regression models were tested, using d0.10, h, both variables (d0.10+h) and d0.102 h as independents variables, applying the following equation (eqn. 1):
During adjustment ([18]), the goodness-of-fit of the model was assessed based on the coefficient of determination (R2), as follows (eqn. 2):
where SSR is the sum of squares of the residuals and SST is the total sum of squares. Validation of selected models from the adjustment phase was carried out using the average deviation (AD - eqn. 3) and absolute average deviation (AAD - eqn. 4):
where Bobs is the observed total or partial biomass, Best is the estimated total or partial biomass and n is the total number of trees.
After selection of the best fitting model based on the above validation process, the model was readjusted using the total dataset (the 30 sampled trees).
Effect of T. pityocampa on biomass increment
The effect of defoliation by the pine processionary moth on pine biomass increment was evaluated at the experimental plot on 83 individually selected pine trees with different degrees of defoliation. The degree of defoliation was visually assessed in March of 2003, 2004 and 2005 using five defoliation classes: class 0 - no defoliation (0%); class 1 - light defoliation (1-25%); class 2 - moderate defoliation (26-50%); class 3 - heavy defoliation (51-75%); and class 4 - very heavy defoliation with almost no foliage remaining (76-100%). Moreover, dendrometric measurements (d0.10, dbh and h) were also carried out in February of 2004, 2005 and 2006.
Using the selected model, total biomass was calculated for those years, biomass increment was estimated for the growing years, as well as the percentage of decrease in biomass increment for undefoliated and defoliated trees. After testing data for normality and variance homogeneity, ANOVA was applied to determine the variance accounted for by the defoliation class effect, the growth years under study and the variance within each defoliation class. The Tukey-Kramer mean separation test was applied in order to determine the biomass increment differences on the basis of defoliation class.
Results
From the biomass models adjustment
The Pearson correlation matrix showed a significant relationship between total biomass, leaf biomass, stem biomass and dendrometric variables with r values varying between 0.64 and 0.94 (Tab. 1). There was a linear relationship between total and partial biomass and the analysed dendrometric variables, mainly d0.10, dbh and h.
Tab. 1 - Pearson’s r correlation matrix between total and partial (leaves and stem) dry biomass and the available dendrometric variables.
- | d0.10 | dbh | h | v | Bstem | Bleaves | Btotal |
---|---|---|---|---|---|---|---|
d0.10 | 1 | - | - | - | - | - | - |
dbh | 0.772 | 1 | - | - | - | - | - |
h | 0.723 | 0.959 | 1 | - | - | - | - |
v | 0.822 | 0.944 | 0.906 | 1 | - | - | - |
Bstem | 0.937 | 0.831 | 0.767 | 0.854 | 1 | - | - |
Bleaves | 0.854 | 0.744 | 0.641 | 0.786 | 0.849 | 1 | - |
Btotal | 0.916 | 0.819 | 0.736 | 0.878 | 0.940 | 0.941 | 1 |
The dendrometric variable with the strongest relationship with the total and partial biomass was d0.10 (Tab. 1). The variable “h” had lower correlation values, although, combined with diameter, model fitting was enhanced. A local model including a single variable (the diameter) as predictor has a limited range of applications. The inclusion of “h” may extend its applicability at a regional scale, allowing a wider range of tree forms to be covered.
Simultaneous use of d0.10 and h as predictors in the models resulted in stronger correlations with leaves, stem and total biomass, giving R2 values ranging from 0.73 to 0.91 (Tab. 2). The most difficult variable to model was leaf biomass while the best one was stem biomass.
Tab. 2 - Coefficient of determination (R2) for different combinations of diameter at 10 cm height (d0.10) and tree height (h) for total and partial biomass for the 16 tested database.
Parameters | Total biomass | Leaves biomass | Stem biomass | ||||||
---|---|---|---|---|---|---|---|---|---|
d0.10 | d0.10 + h | d0.102 h | d0.10 | d0.10 + h | d0.102 h | d0.10 | d0.10 + h | d0.102 h | |
Maximum | 0.877 | 0.898 | 0.898 | 0.839 | 0.841 | 0.83 | 0.941 | 0.95 | 0.942 |
Average | 0.825 | 0.835 | 0.823 | 0.728 | 0.733 | 0.697 | 0.886 | 0.907 | 0.902 |
Minimum | 0.69 | 0.736 | 0.736 | 0.555 | 0.585 | 0.583 | 0.756 | 0.805 | 0.803 |
Std. dev | 0.049 | 0.042 | 0.041 | 0.074 | 0.07 | 0.065 | 0.044 | 0.037 | 0.037 |
The adjusted model tendency, measured by the average deviation (AD) showed an overestimation of true biomass values for total and leaves biomass and an underestimation for stem biomass (Tab. 3). Absolute average deviation (AAD) values showed that simultaneous use of d0.10 and h provides the best biomass predictions.
Tab. 3 - The Average deviation (AD) and the Absolute Average deviation (AAD) from the validation results across the 16 sampled groups.
Parameter | Total biomass | Leaves biomass | Stem biomass | ||||||
---|---|---|---|---|---|---|---|---|---|
d0.10 | d0.10 +h | d0.102 h | d0.10 | d0.10 +h | d0.102 h | d0.10 | d0.10 +h | d0.102 h | |
AD | 0.109 | 0.284 | 0.433 | -0.014 | 0.025 | 0.180 | -0.105 | -0.029 | 0.050 |
AAD | 1.994 | 1.930 | 1.995 | 1.113 | 1.112 | 1.152 | 0.680 | 0.641 | 0.662 |
The final models for partial and total biomass, adjusted to the total dataset, are reported below (eqn. 5, eqn. 6, eqn. 7):
where Btotal, Bleaves and Bstem refer to the total, leaves and stem biomass model, respectively, and RMSE is the root mean square error.
The effect of Thaumetopoea pityocampa on biomass increment
ANOVA results for the biomass increment of pine trees attacked by the pine processionary moth using different classes of defoliation (Tab. 4) showed that the interaction between growing years, defoliation class and year itself were not statistically significant. However, a significant difference was found for defoliation class (P<0.001), revealing that the increment of pine tree biomass was affected by the both presence of the insect and by the intensity of the attack.
Tab. 4 - ANOVA results for biomass increment of Pinus pinaster Ait. after defoliation by the pine processionary moth.
Source | Sum of squares | DF | Mean square | F-value | P-value |
---|---|---|---|---|---|
Year | 10.66 | 1 | 10.66 | 0.043 | 0.84 |
Defoliation class | 15882.37 | 4 | 3970.59 | 16.07 | 0.00 |
Year · Defoliation class | 1535.46 | 4 | 383.86 | 1.55 | 0.19 |
Residual (trees/defoliation class) | 38532.39 | 156 | 247.00 | - | - |
The average biomass increment was maximum for undefoliated trees (56.6 kg per tree) and minimum for heavily and very heavily defoliated trees (16.1 and 15.1 kg per tree respectively - Tab. 5). The results of Tukey-Kramer mean separation test indicated that biomass increment from trees suffering heavy and very heavy defoliation (classes 3 and 4 respectively) did significantly differ from those with moderate defoliation (class 2), light defoliation (class 1) and undefoliated trees (class 0). There was no significant difference between undefoliated trees and those with light defoliation. The percent decrease in biomass increment during the growing years that were studied was around 70% for heavy and very heavy defoliated trees, 50% for moderated defoliated trees and 37% for light defoliated trees.
Tab. 5 - Biomass increment (mean and SD) of Pinus pinaster and results of Tukey-Kramer mean separation test based on classes of defoliation. (*): Means followed by the same letter in the same column are not significantly different (significance level: α < 0.05).
Defoliation class |
N | Mean of biomass increment (Kg/tree) |
Percent decrease in biomass increment |
---|---|---|---|
4 | 22 | 15.14 ± 2.76 a* | 73.24 |
3 | 52 | 16.14 ± 2.84 a | 71.47 |
2 | 44 | 28.93 ± 5.28 b | 48.86 |
1 | 44 | 35.54 ± 10.50 bc | 37.17 |
0 | 4 | 56.57 ± 2.06 c | - |
Discussion
The first aim of this study was to develop biomass growth models for different above-ground components (leaves, stem and total biomass) of pine trees. Results showed that d0.10 and h were the best predictor variables when used in conjunction, in contrast to the d2h variable tested by Mikšys et al. ([18]). The validation process did confirm this option. The most difficult variable to model was leaf biomass (R2=0.74), since young pine trees are shown to have very heterogeneous crowns. Indeed, the crown of oldest trees tend to become more homogeneous. Stem biomass tended to be much more homogeneous (R2=0.84), being possible to obtain more precise models for stem and total biomass estimation. Since total biomass combines the effect of both leaves and stem in biomass calculation, the evaluation of the effect of the pine processionary moth on biomass growth was assessed using the total biomass model.
The second goal was to evaluate potential biomass increment after defoliation by the pine processionary moth. Results indicated that the degree of defoliation was a decisive factor in tree biomass increment. Losses of about 49% in biomass increment were observed in moderately defoliated trees (class 2) while losses of about 71-73% were registered in heavily attacked trees (classes 3 and 4; class 4 comprised completed defoliated trees). Moderately defoliated trees values were in agreement with those determined by Markalas ([16]), but in completely defoliated trees, the impact was higher. Our results also show that the consequences of T. pityocampa activity can be detected immediately after the attack, with biomass losses reported in the same year as the infestation. Carus ([5]) also found a sharp decline in host pine growth during and directly after a pine processionary outbreak. However [13] reported a decrease in ring growth in the year following defoliation. This may be due to tree physiology, since the total biomass responds in a different way than tree ring growth. However, we realize that these conclusions are the result of a simplistic approach of the problem. The ecosystem and its dynamics are much more complex than we had assumed. While the present study was exclusively focused on the arboreal stratum, next stages should also analyse the impact on tree defoliation on shrubs biomass dynamics. We need to understand how defoliation of trees, increasing light received by the understorey, can lead to an increment of shrubs biomass. Furthermore, we must consider that the presence of defoliating insects on the tree canopy can increases nutrients input on top soil layers, changing its composition. For example, Lovett et al. ([15]) has concluded that insect defoliation represents a major perturbation to the internal N cycle of the forest, but this perturbation primarily causes a redistribution of N within the ecosystem rather than a large loss of N. Therefore, among the topics that deserve further research are the impacts of the pine processionary moth on the biomass dynamics of the entire ecosystem.
Conclusions
In conclusion, this work clearly shows that allometric models can be used to estimate the impact of insects on forest dynamics. Furthermore, results showed that the negative effects of insect attack on the biomass growth are visible in the same year at the occurrence of defoliation with a reduction of the biomass increment that is directly proportional to the intensity of the attack.
So far, there was only a notion that attack by the pine processionary moth had an important impact on the Pinus forests biomass growth. However, our results provide a tool to economically quantify these impacts. This information is important for forest owners and managers, due to the high economic importance of Pinus forest and the potential effect of these impacts on the Portuguese economy. Results from this study indicate that after heavy defoliation, losses can represent around 100 € per hectare, that means 12.6 million € for the entire country.
However, this problem cannot be analysed only from an economic point of view since the ecological importance of these attacks is also relevant. The conservation of pine forests requires appropriate management techniques to counterbalance the negative effects of T. pityocampa. Forestry personnel should carefully plan all the new pine plantations and the ecological range of the species should be adhered to avoid additional tree stress, thereby preventing insect attacks.
Further studies should be carried out in order to better understand these phenomena, taking into account uncertainties such as the effect of climate change, which would expose forest stands to even greater stress.
Acknowledgements
The authors would like to thank CITAB and the projects PTDC/AGR-CFL/68186/ 2006 and PTDC/AGR-CFL/69733/2006 for the financial support for this study.
References
Gscholar
Gscholar
Gscholar
Gscholar
Gscholar
Gscholar
Gscholar
Gscholar
Gscholar
Gscholar
Gscholar
Gscholar
Gscholar
Authors’ Info
Authors’ Affiliation
S Chacim
D Lopes
University of Trás-os-Montes and Alto Douro, Forestry and Landscape Architecture Department, 5001-801 Vila Real (Portugal)
D Lopes
Center for the Research and Technology of Agro-Environmental and Biological Sciences, University of Trás-os-Montes and Alto Douro, 5001-801 Vila Real (Portugal)
Corresponding author
Paper Info
Citation
Arnaldo PS, Chacim S, Lopes D (2010). Effects of defoliation by the pine processionary moth Thaumetopoea pityocampa on biomass growth of young stands of Pinus pinaster in northern Portugal. iForest 3: 159-162. - doi: 10.3832/ifor0553-003
Paper history
Received: Jun 07, 2010
Accepted: Oct 23, 2010
First online: Nov 15, 2010
Publication Date: Nov 15, 2010
Publication Time: 0.77 months
Copyright Information
© SISEF - The Italian Society of Silviculture and Forest Ecology 2010
Open Access
This article is distributed under the terms of the Creative Commons Attribution-Non Commercial 4.0 International (https://creativecommons.org/licenses/by-nc/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.
Web Metrics
Breakdown by View Type
Article Usage
Total Article Views: 53152
(from publication date up to now)
Breakdown by View Type
HTML Page Views: 44991
Abstract Page Views: 2934
PDF Downloads: 3889
Citation/Reference Downloads: 26
XML Downloads: 1312
Web Metrics
Days since publication: 5073
Overall contacts: 53152
Avg. contacts per week: 73.34
Article Citations
Article citations are based on data periodically collected from the Clarivate Web of Science web site
(last update: Feb 2023)
Total number of cites (since 2010): 12
Average cites per year: 0.86
Publication Metrics
by Dimensions ©
Articles citing this article
List of the papers citing this article based on CrossRef Cited-by.
Related Contents
iForest Similar Articles
Research Articles
Abundance and impact of egg parasitoids on the pine processionary moth (Thaumetopoea pityocampa) in Bulgaria
vol. 14, pp. 456-464 (online: 02 October 2021)
Research Articles
Plant phenotype affects oviposition behaviour of pine processionary moth and egg survival at the southern edge of its range
vol. 11, pp. 572-576 (online: 01 September 2018)
Research Articles
Performances of an expanding insect under elevated CO2 and snow cover in the Alps
vol. 1, pp. 126-131 (online: 27 August 2008)
Research Articles
Role of serotiny on Pinus pinaster Aiton germination and its relation to mother plant age and fire severity
vol. 12, pp. 491-497 (online: 02 November 2019)
Research Articles
Influence of tree density on climate-growth relationships in a Pinus pinaster Ait. forest in the northern mountains of Sardinia (Italy)
vol. 8, pp. 456-463 (online: 19 October 2014)
Research Articles
Long-term implications of traditional forest regulation methods applied to Maritime pine (Pinus pinaster Ait.) forests in central Spain: a century of management plans
vol. 3, pp. 33-38 (online: 02 March 2010)
Short Communications
Upscaling the estimation of surface-fire rate of spread in maritime pine (Pinus pinaster Ait.) forest
vol. 7, pp. 123-125 (online: 13 January 2014)
Research Articles
Age and thinning effects on elemental composition of Pinus pinaster and Pinus radiata needles
vol. 14, pp. 299-306 (online: 08 June 2021)
Research Articles
Impact of wildfire on the edaphic microarthropod community in a Pinus pinaster forest in central Italy
vol. 8, pp. 874-883 (online: 05 May 2015)
Research Articles
Simplified methods to inventory the current annual increment of forest standing volume
vol. 5, pp. 276-282 (online: 17 December 2012)
iForest Database Search
Search By Author
Search By Keyword
Google Scholar Search
Citing Articles
Search By Author
Search By Keywords
PubMed Search
Search By Author
Search By Keyword