Impact of insect herbivory on plant stress volatile emissions from trees : A synthesis of quantitative measurements and recommendations for future research

Plants, and particularly trees, are the largest source of atmospheric volatile organic compounds globally. Insect herbivory alters plant volatile emission rates and the types of compounds that are emitted. These stress volatiles are a major contribution to total plant volatile emissions during active herbivore feeding, with important implications for atmospheric chemistry processes. However, emission models do not currently have a quantitative description of plant volatile emission rates during and after insect herbivore feeding. This review provides a brief background on plant volatile organic compounds, the urgency of including biotic stress emissions in models, and a summary of plant volatile emission models and steps they have taken to incorporate stress emissions into their framework. The review ends with a synthesis of volatile emissions from trees during insect herbivory. This synthesis highlights key gaps in studied systems and measurement approaches. We provide a summary of recommendations for future work to address these gaps, improve comparability between studies, and generate the data-sets we need to develop a descriptive model of these plant stress volatile emissions.Peer reviewe

Investigating evidence of enhanced aerosol formation and growth due to autumnal moth larvae feeding on mountain birch at SMEAR I in northern Finland

Laboratory studies have shown that heibivory-induced biogenic volatile organic compound (BVOC) emissions might enhance aerosol formation and growth. To increase understanding of the atmospheric relevance of this enhancement, we analyzed 25 years of data from SMEAR I (Station for Measuring Ecosystem-Atmosphere Relations) in northern Finland, where autumnal moth (Epirrita autumnata) larvae are prominent defoliators of mountain birch. We did not find a direct correlation between the autumnal moth density and aerosol processes, nor between the total number concentration and temperature, and hence the basal BVOC emissions. Instead, there is some evidence that the total particle concentration is elevated even for a few years after the infestation due to delayed defense response of mountain birch. The low total biomass of the trees concomitantly with low autumnal moth densities during most of the years of our study, may have impacted our results, hindering the enhancement of aerosol processes.Peer reviewe

Simulating ozone dry deposition at a boreal forest with a multi-layer canopy deposition model

A multi-layer ozone (O-3) dry deposition model has been implemented into SOSAA (a model to Simulate the concentrations of Organic vapours, Sulphuric Acid and Aerosols) to improve the representation of O-3 concentration and flux within and above the forest canopy in the planetary boundary layer. We aim to predict the O-3 uptake by a boreal forest canopy under varying environmental conditions and analyse the influence of different factors on total O-3 uptake by the canopy as well as the vertical distribution of deposition sinks inside the canopy. The newly implemented dry deposition model was validated by an extensive comparison of simulated and observed O-3 turbulent fluxes and concentration profiles within and above the boreal forest canopy at SMEAR II (Station to Measure Ecosystem-Atmosphere Relations II) in Hyytiala, Finland, in August 2010. In this model, the fraction of wet surface on vegetation leaves was parametrised according to the ambient relative humidity (RH). Model results showed that when RH was larger than 70% the O-3 uptake onto wet skin contributed similar to 51% to the total deposition during nighttime and similar to 19% during daytime. The overall contribution of soil uptake was estimated about 36 %. The contribution of sub-canopy deposition below 4.2m was modelled to be similar to 38% of the total O-3 deposition during daytime, which was similar to the contribution reported in previous studies. The chemical contribution to O-3 removal was evaluated directly in the model simulations. According to the simulated averaged diurnal cycle the net chemical production of O-3 compensated up to similar to 4% of dry deposition loss from about 06:00 to 15:00 LT. During nighttime, the net chemical loss of O-3 further enhanced removal by dry deposition by a maximum similar to 9 %. Thus the results indicated an overall relatively small contribution of airborne chemical processes to O-3 removal at this site.Peer reviewe

Boreal forest BVOC exchange : emissions versus in-canopy sinks

A multilayer gas dry deposition model has been developed and implemented into a one-dimensional chemical transport model SOSAA (model to Simulate the concentrations of Organic vapours, Sulphuric Acid and Aerosols) to calculate the dry deposition velocities for all the gas species included in the chemistry scheme. The new model was used to analyse in-canopy sources and sinks, including gas emissions, chemical production and loss, dry deposition, and turbulent transport of 12 featured biogenic volatile organic compounds (BVOCs) or groups of BVOCs (e.g. monoterpenes, isoprene+2-methyl-3-buten-2-ol (MBO), sesquiterpenes, and oxidation products of mono-and sesquiterpenes) in July 2010 at the boreal forest site SMEAR II (Station for Measuring Ecosystem-Atmosphere Relations). According to the significance of modelled monthly-averaged individual source and sink terms inside the canopy, the selected BVOCs were classified into five categories: 1. Most of emitted gases are transported out of the canopy (monoterpenes, isoprene + MBO). 2. Chemical reactions remove a significant portion of emitted gases (sesquiterpenes). 3. Bidirectional fluxes occur since both emission and dry deposition are crucial for the in-canopy concentration tendency (acetaldehyde, methanol, acetone, formaldehyde). 4. Gases removed by deposition inside the canopy are compensated for by the gases transported from above the canopy (acetol, pinic acid, beta-caryophyllene's oxidation product BCSOZOH). 5. The chemical production is comparable to the sink by deposition (isoprene's oxidation products ISOP34OOH and ISOP34NO3). Most of the simulated sources and sinks were located above about 0.2 h(c) (canopy height) for oxidation products and above about 0.4 h(c) for emitted species except formaldehyde. In addition, soil deposition (including deposition onto understorey vegetation) contributed 11-61% to the overall in-canopy deposition. The emission sources peaked at about 0.8-0.9 h(c), which was higher than 0.6 h(c) where the maximum of dry deposition onto overstorey vegetation was located. This study provided a method to enable the quantification of the exchange between atmosphere and biosphere for numerous BVOCs, which could be applied in large-scale models in future. With this more explicit canopy exchange modelling system, this study analysed both the temporal and spatial variations in individual in-canopy sources and sinks, as well as their combined effects on driving BVOC exchange. In this study 12 featured BVOCs or BVOC groups were analysed. Other compounds could also be investigated similarly by being classified into these five categories.Peer reviewe

Long-term total OH reactivity measurements in a boreal forest

Corrigendum: The legend in Fig. 6e has been mislabeled. The gray colorcorresponds to “Missing” and the other colors should havecorresponded to the same species as in Fig. 6f. The figure,which is also the key figure of the article, can be found belowwith the correct legend.Total hydroxyl radical (OH) reactivity measurements were conducted at the second Station for Measuring Ecosystem-Atmosphere Relations (SMEAR II), a boreal forest site located in Hyytiala, Finland, from April to July 2016. The measured values were compared with OH reactivity calculated from a combination of data from the routine trace gas measurements (station mast) as well as online and offline analysis with a gas chromatographer coupled to a mass spectrometer (GC-MS) and offline liquid chromatography. Up to 104 compounds, mostly volatile organic compounds (VOCs) and oxidized VOCs, but also inorganic compounds, were included in the analysis, even though the data availability for each compound varied with time. The monthly averaged experimental total OH reactivity was found to be higher in April and May (ca. 20 s(-1)) than in June and July (7.6 and 15.4 s(-1), respectively). The measured values varied much more in spring with high reactivity peaks in late afternoon, with values higher than in the summer, in particular when the soil was thawing. Total OH reactivity values generally followed the pattern of mixing ratios due to change of the boundary layer height. The missing reactivity fraction (defined as the difference between measured and calculated OH reactivity) was found to be high. Several reasons that can explain the missing reactivity are discussed in detail such as (1) missing measurements due to technical issues, (2) not measuring oxidation compounds of detected biogenic VOCs, and (3) missing important reactive compounds or classes of compounds with the available measurements. In order to test the second hypothesis, a one-dimensional chemical transport model (SOSAA) has been used to estimate the amount of unmeasured oxidation products and their expected contribution to the reactivity for three different short periods in April, May, and July. However, only a small fraction (<4.5 %) of the missing reactivity can be explained by modelled secondary compounds (mostly oxidized VOCs). These findings indicate that compounds measured but not included in the model as well as unmeasured primary emissions contribute the missing reactivity. In the future, non-hydrocarbon compounds from sources other than vegetation (e.g. soil) should be included in OH reactivity studies.Peer reviewe

Mitigation Impact of Different Harvest Scenarios of Finnish Forests That Account for Albedo, Aerosols, and Trade-Offs of Carbon Sequestration and Avoided Emissions

The pressure to increase forest and land carbon stocks simultaneously with increasing forest based biomass harvest for energy and materials emphasizes the need for dedicated analyses of impacts and possible trade-offs between these different mitigation options including also forest related biophysical factors, surface albedo and the formation of biogenic Secondary Organic Aerosols (SOA). We analyzed the change in global radiative forcing (RF) due to changes in these climatic agents as affected by the change in state of Finnish forests under increased or decreased harvest scenarios from a baseline. We also included avoided emissions due to wood material and energy substitution. Increasing harvests from baseline (65% of Current Annual Increment) decreased the total carbon sink (carbon in trees, soil and harvested wood products) at least for 50 years. When we coupled this change in carbon with other biosphere responses, surface albedo and aerosols, decreasing harvests from the baseline produced the largest cooling effect during 50 years. Accounting also for the avoided emissions due to increased wood use, the RF responses of the two lowest harvest scenarios were within uncertainty range. Our results show that the effects of forest management on SOA formation should be included in the analyses trying to deduce the net climate impact of forest use. The inclusion of the rarely considered SOA effects enforces the view that the lower the harvest, the more climatic cooling boreal forests provide. These results should act as a caution mark for policy makers who are emphasizing the increased utilization of forest biomass for short-living products and bioenergy as an efficient measure to mitigate climate change.Peer reviewe

A modelling study of OH, NO3 and H2SO4 in 2007– 2018 at SMEAR II, Finland : analysis of long-term trends

Major atmospheric oxidants (OH, O3 and NO3) dominate the atmospheric oxidation capacity, while H2SO4 is considered as a main driver for new particle formation. Although numerous studies have investigated the long-term trend of ozone in Europe, the trends of OH, NO3 and H2SO4 at specific sites are to a large extent unknown. The one-dimensional model SOSAA has been applied in several studies at the SMEAR II station and has been validated by measurements in several projects. Here, we applied the SOSAA model for the years 2007–2018 to simulate the atmospheric chemical components, especially the atmospheric oxidants OH and NO3, as well as H2SO4 at SMEAR II. The simulations were evaluated with observations from several shorter and longer campaigns at SMEAR II. Our results show that daily OH increased by 2.39% per year and NO3 decreased by 3.41% per year, with different trends of these oxidants during day and night. On the contrary, daytime sulfuric acid concentrations decreased by 2.78% per year, which correlated with the observed decreasing concentration of newly formed particles in the size range of 3– 25 nm with 1.4% per year at SMEAR II during the years 1997–2012. Additionally, we compared our simulated OH, NO3 and H2SO4 concentrations with proxies, which are commonly applied in case a limited number of parameters are measured and no detailed model simulations are available.Peer reviewe

Characterization of total ecosystem-scale biogenic VOC exchange at a Mediterranean oak–hornbeam forest

Recently, the number and amount of biogenically emitted volatile organic compounds (VOCs) has been discussed in great detail. Depending on the ecosystem, the published number varies between a dozen and several hundred compounds. We present ecosystem exchange fluxes from a mixed oak-hornbeam forest in the Po Valley, Italy. The fluxes were measured by a proton transfer reaction-time-of-flight (PTR-ToF) mass spectrometer and calculated using the eddy covariance (EC) method. Detectable fluxes were observed for up to 29 compounds, dominated by isoprene, which comprised over 60aEuro-% of the total upward flux (on a molar basis). The daily average of the total VOC upward flux was 10.4aEuro-nmolaEuro-m(-2)aEuro-s(-1). Methanol had the highest concentration and accounted for the largest downward flux. Methanol seemed to be deposited to dew, as the downward flux happened in the early morning, right after the calculated surface temperature came closest to the calculated dew point temperature. We estimated that up to 30aEuro-% of the upward flux of methyl vinyl ketone (MVK) and methacrolein (MACR) originated from atmospheric oxidation of isoprene. A comparison between two methods for the flux detection (manual and automated) was made. Their respective advantages and disadvantages were discussed and the differences in their results shown. Both provide comparable results.Peer reviewe

The importance of accounting for enhanced emissions of monoterpenes from new Scots pine foliage in models - A Finnish case study

Peer reviewe

Modelling the influence of biotic plant stress on atmospheric aerosol particle processes throughout a growing season

Most trees emit volatile organic compounds (VOCs) continuously throughout their life, but the rate of emission and spectrum of emitted VOCs become substantially altered when the trees experience stress. Despite this, models to predict the emissions of VOCs do not account for perturbations caused by biotic plant stress. Considering that such stresses have generally been forecast to increase in both frequency and severity in the future climate, the neglect of stress-induced plant emissions in models might be one of the key obstacles for realistic climate change predictions, since changes in VOC concentrations are known to greatly influence atmospheric aerosol processes. Thus, we constructed a model to study the impact of biotic plant stresses on new particle formation and growth throughout a full growing season. We simulated the influence on aerosol processes caused by herbivory by the European gypsy moth (Lymantria dispar) and autumnal moth (Epirrita autumnata) feeding on pedunculate oak (Quercus robur) and mountain birch (Betula pubescens var. pumila), respectively, and also fungal infections of pedunculate oak and balsam poplar (Populus balsamifera var. suaveolens) by oak powdery mildew (Erysiphe alphitoides) and poplar rust (Melampsora larici-populina), respectively. Our modelling results indicate that all the investigated plant stresses are capable of substantially perturbing both the number and size of aerosol particles in atmospherically relevant conditions, with increases in the amount of newly formed particles by up to about an order of magnitude and additional daily growth of up to almost 50 nm. We also showed that it can be more important to account for biotic plant stresses in models for local and regional predictions of new particle formation and growth during the time of infestation or infection than significant variations in, e.g. leaf area index and temperature and light conditions, which are currently the main parameters controlling predictions of VOC emissions. Our study thus demonstrates that biotic plant stress can be highly atmospherically relevant. To validate our findings, field measurements are urgently needed to quantify the role of stress emissions in atmospheric aerosol processes and for making integration of biotic plant stress emission responses into numerical models for prediction of atmospheric chemistry and physics, including climate change projection models, possible.Peer reviewe