Process-based modelling of biogenic monoterpene emissions : sensitivity to temperature and light

Process-based modelling of biogenic monoterpene emissions: sensitivity to temperature and light G. Schurgers, A. Arneth, R. Holzinger, and A. Goldstein Lund University, Department of Physical Geography and Ecosystems Analysis, Sölvegatan 12, 223 62 Lund, Sweden Utrecht University, Institute for Marine and Atmospheric Research, Utrecht, The Netherlands University of California at Berkeley, Department of Environmental Science, Policy and Management, Berkeley, CA, USA Received: 25 September 2008 – Accepted: 29 October 2008 – Published: 6 January 2009 Correspondence to: G. Schurgers (guy.schurgers@nateko.lu.se) Published by Copernicus Publications on behalf of the European Geosciences Union.

based) method, which causes a difference in upscaling from instantaneous to daily emissions, requires roughly a doubling of emission capacities to bridge the gap to production capacities.The sensitivities to changes in temperature and light were tested for the new methods, the temperature sensitivity was slightly higher than that of the short-term temperature dependent algorithm.
Applied on a global scale, the first algorithm resulted in annual total emissions of 29.6 Tg C a −1 , the second algorithm resulted in 31.8Tg C a −1 when applying the correction factor 2 between emission capacities and production capacities.However, the exact magnitude of such a correction is spatially varying and hard to determine as a global average.

Introduction
Biogenic emissions of monoterpenes influence atmospheric composition and air quality, especially on a regional scale.Monoterpene oxidation in the atmosphere contributes to production of ozone (O 3 ) in the presence of nitrogen oxides (NO x ) (Jenkin and Clemitshaw, 2000).Monoterpenes also react directly with O 3 , forming low volatility oxidation products that are important sources for secondary organic aerosol (SOA) formation and growth (Hoffmann et al., 1997;Aumont et al., 2000;Chung and Seinfeld, 2002;Tsigaridis and Kanakidou, 2003;Simpson et al., 2007).SOA yield from monoterpene ozonolysis is considered relatively large, although knowledge on many of the processes involved is still scarce (Tsigaridis and Kanakidou, 2003).Since the annual global SOA production from terrestrial biogenic volatile organics might exceed SOA production from anthropogenic VOC by more than a factor of ten, and could be of same order of magnitude as the production of sulphate particles (Tsigaridis and Kanakidou, 2003), the role of monoterpenes for radiative transfer and cloud properties is probably significant.However, at the same time their regional and global emission patterns are not very well known, and effects of changing climate, atmospheric CO 2 concentration or human land cover and land use change are uncertain.The incorporation of process-understanding related to their cellular production in global vegetation models can help to investigate these effects, as these models are applicable to a wider range of environmental conditions, including global change related questions.
Monoterpenes are produced along the chloroplastic DXP pathway, in a reaction chain that is, except for the final steps, similar to the formation of isoprene (Lichtenthaler et al., 1997).This metabolic pathway is closely linked to photosynthesis through one of the chief precursors, glyceraldehyde-3-phosphate, originating from the chloroplastic Calvin cycle, and the requirement of energy for the reduction of the precursor carbohydrates (Lichtenthaler et al., 1997).Unlike isoprene, monoterpenes and other less volatile compounds can be stored in leaves, either as nonspecific storage (Niinemets and Reichstein, 2002)  such as glandular trichomes (e.g.Gershenzon et al., 1989;Turner et al., 2000), resin canals, or resin ducts (e.g.Franceschi et al., 2005).This specific storage of monoterpenes within a leaf is hence built up during leaf development (Gershenzon et al., 2000;McConkey et al., 2000;Turner et al., 2000).The release of stored monoterpenes is mainly driven by changes in monoterpene vapour pressure, which is primarily determined by temperature (Dement et al., 1975;Tingey et al., 1980).This temperature-driven release from storage has led to the development of an algorithm for emission of monoterpenes (Tingey et al., 1980;Guenther et al., 1993), which has been successfully applied to interpret measurements on leaf or canopy scale (e.g.Ruuskanen et al., 2005;Holzinger et al., 2006), and is generally used for estimates of global monoterpene emissions (e.g.Guenther et al., 1995;Naik et al., 2004;Lathi ère et al., 2006).
Although monoterpene emissions of many species have been shown to depend primarily on temperature on a relatively short time scale of hours to days, the seasonal variation in monoterpene emissions cannot be explained by temperature response alone (Yokouchi et al., 1984;Staudt et al., 2000;Holzinger et al., 2006).Long-term (∼annual) changes in emissions were so far represented by seasonally varying emission capacities on a local scale (Staudt et al., 2000), although it is not clear whether the observed seasonal variation is related to the dynamics of the monoterpene storage or to the rate of production.On a global scale such changes are ignored, and the temperature-dependent algorithm was used for annual emission estimates so far (e.g.Naik et al., 2004;Lathi ère et al., 2005).What is more, over recent years an increasing number of studies have identified monoterpene emissions, particularly in broadleaf species, to respond to temperature and light in a pattern similar to that found for isoprene, e.g. for Quercus ilex (Staudt and Seufert, 1995;Bertin et al., 1997;Ciccioli et al., 1997;Staudt and Bertin, 1998), Fagus sylvatica (Schuh et al., 1997;Dindorf et al., 2006), Helianthus annuus (Schuh et al., 1997), several mediterranean species (Owen et al., 2002), Apeiba tibourbou (Kuhn et al., 2004), Hevea brasiliensis (Wang et al., 2007) and other tropical plant species or land cover types (Greenberg et al., 2003; Otter et al., 2003).In these species, emission takes place directly after production, without intermediate storage within the leaf, in a pattern similar to that observed for isoprene.The observed dependences reflect those of monoterpene synthesis, which is closely linked to photosynthesis.These findings suggest that modelling of monoterpene emissions for some regions will have to be revised, which will likely affect global emission estimates as well.
A limited number of studies have attempted to express monoterpene production explicitly, linking it to processes of carbon assimilation in the chloroplast (Niinemets et al., 2002;B äck et al., 2005;Grote et al., 2006), and hence being dependent on both temperature and light.Storage of monoterpenes can then be included as an additional feature to account for the observed short-term temperature dependence of monoterpene emissions (Niinemets and Reichstein, 2002;B äck et al., 2005).The release from storage can modulate emissions over periods of days to months: Mihaliak et al. (1991) showed that monoterpenes in intact plants of Mentha×piperita are stored in a stable pool for several weeks, and Gershenzon et al. (1993) found for several monoterpenestoring species no significant amount of labelled monoterpenes to be released for 8 to 12 days after a pulse of 14 CO 2 .These long-term (∼annual) changes in emissions originating from changes in the specific storage (e.g.glands or resin ducts) have not been included in modelling studies so far.
Our chief objective here is to investigate the effects of an explicit representation of chloroplastic and leaf processes on seasonal to annual monoterpene emission patterns.We develop a model for light-and temperature-dependent monoterpene emissions by combining a process-based description of monoterpene production and a temperature-dependent residence in specific storage organs within the plant.The model is implemented in the dynamic global vegetation model (DGVM) framework LPJ-GUESS (Smith et al., 2001;Sitch et al., 2003) to investigate the sensitivity of emissions to temperature and light and the use of monoterpene storage as a measure to distinguish between production and emission of monoterpenes.The goal is to create a tool that builds on process understanding and that can be used to investigate interactions of Introduction

Conclusions References
Tables Figures

Back Close
Full climate change, vegetation dynamics, vegetation productivity and trace gas emissions over periods from years to millennia within a consistent modelling framework.In this study, we concentrate on model parameterization and evaluation using observations of monoterpene emissions from a Ponderosa pine plantation.Model sensitivities for this site and implications for application on global scale will be discussed.

Short-term monoterpene emission
In those plant species that display a light-independent, temperature-driven monoterpene emission pattern, these emissions usually originate from non-specific (e.g., dissolved in the cytosol) or specific (e.g., glands, resin ducts) storage pools within leaves.
The storage pools act as a continuous source of monoterpenes, with emissions driven by changes in monoterpene vapour pressures (Dement et al., 1975;Tingey et al., 1980), hence the clear temperature dependence.Typically, an exponential algorithm as presented by Tingey et al. (1980) and Guenther et al. (1993) is used to simulate these emissions: In this equation, M is the monoterpene emission (µg g −1 h −1 ), M s is the emission rate under standard conditions (referred to as emission capacity), β is a constant (0.09 K −1 ), T is leaf temperature (K), and T s is the standard temperature (303 K).Simulations with this algorithm were performed for a broad range of species, specifically for many conifers, e.g. for Pinus elliottii (Tingey et al., 1980), P. ponderosa (Holzinger et al., 2006), and P. sylvestris (Ruuskanen et al., 2005).The temperature-dependent algorithm in Eq. ( 1) is useful for modelling the short-term emission response to temperature, as it reflects the changes in vapour pressure due to temperature, but changes in vapour Introduction

Conclusions References
Tables Figures

Back Close
Full pressure from changes in the concentrations in the storage pool of monoterpenes are not covered by the algorithm.

Monoterpene production, storage and emission
The algorithm presented above reflects the short-term dependence of monoterpene emissions from temperature.In order to simulate both the long-term changes and the short-term changes, we split the simulation in two parts: the production of monoterpenes, following a process-based approach based on the energy requirements of monoterpene synthesis, and the emission of monoterpenes, following an approach equivalent to Eq. ( 1).Between production and emission, monoterpenes can be stored for periods of different length.
Monoterpene production is simulated following Niinemets et al. (2002), who calculate the production of monoterpenes in two Quercus species based on the chloroplastic electron transport rate required to drive terpene synthesis: In these equations, J is the photosynthetic electron flux (mol m −2 h −1 ), is the fraction of this flux that is available for monoterpene production, and α converts the electron flux into monoterpenes (g mol −1 ).The fraction depends on temperature and on a species-specific electron fraction s , which forms a similar scalar to emissions as the emission capacities or standard emission rates (M s ) that are usually reported for species do.f T is a temperature factor accounting for the higher temperature optimum of terpene production observed, as was done for isoprene (Arneth et al., 2007).s , the fraction of the electron flux under standard conditions, can be derived directly from the emission capacity by calculation of photosynthesis and hence J at standard conditions.
Apart from the standard temperature of 30 1000 µmol m −2 s −1 PAR (as is standard for isoprene), even though this is not formally determined for monoterpenes that are emitted temperature-dependently. Produced monoterpenes resulting from the process-based method of Eqs. ( 2) and (3) can be stored for shorter or longer periods in a specific storage pool with size m (in g m −2 ground area).This specific storage of monoterpenes within a leaf takes place in storage organs such as glands or resin ducts.We ignore the dynamics of nonspecific storage, as it is of minor interest to this study due to its short timescale of up to several hours.Specific storage is represented with a single storage pool, and we assume that the release from specific storage depends on temperature in a similar way as the release from non-specific storage.The size of the storage pool is determined by changes in production M prod and release M emis : The concentration of monoterpenes within needles has been shown to affect the emission of monoterpenes from a number of conifer species, e.g.Pinus ponderosa (Lerdau et al., 1994), Pseudotsuga menziesii (Lerdau et al., 1995), Picea mariana and Pinus banksiana (Lerdau et al., 1997).Therefore, the release from the storage pool is simulated depending on the pool size (which is related to monoterpene concentration) with an average residence time τ.
The average residence time τ (in days) is determined under standard temperature T s and is adjusted for other temperatures with a Q 10 -relationship: The short-term temperature response of the monoterpene emission of Eqs. ( 5) and ( 6), i.e. the response with negligible changes in m, is adopted from the short-term Introduction

Conclusions References
Tables Figures

Back Close
Full The temperature dependence in their equation 8 results in a value for Q 10 between 1.8 and 2.0 for temperatures between 0 and 30 • C. For our modelling exercises we use a constant Q 10 of 1.9.This is somewhat lower than the temperature response of Eq. ( 1), which results in a Q 10 -value of 2.5 (with β=0.09).The value for τ s was varied in a set of sensitivity tests and will be discussed below.
The presented algorithm thus calculates monoterpene production according to the availability of temperature and light, closely linked to photosynthesis.The produced monoterpenes are then emitted depending on the temperature and on the amount (or concentration) in the leaves.The storage of monoterpenes and the temperaturedependent emission from storage modulates the seasonal cycle of the emissions, but does not affect the annually produced total.

Implementation in a DGVM framework and experiment setup
For comparative analysis, both the short-term temperature-dependent monoterpene emission algorithm from Eq. ( 1) and the process-based production and emission from Eqs. (2) to (6) were implemented within the dynamic global vegetation model framework LPJ-GUESS (Smith et al., 2001;Sitch et al., 2003).LPJ-GUESS simulates vegetation distribution as well as the cycles of carbon and water within the vegetation and the soil.The model calculates photosynthesis, and hence the total electron flow J required for Eq. ( 2) on a daily timestep.LPJ-GUESS can be applied as a gapmodel (Smith et al., 2001), where several age cohorts of one species or PFT, which compete for light and water, can occur in one gridcell.In this way, canopy successional dynamics are represented in a realistic manner, and modelling of vegetation dynamics on tree species level is possible (Hickler et al., 2004;Arneth et al., 2008b).
For both methods, the extrapolation from leaf-level to canopy-level within the DGVM is done linearly with the fraction of the radiation absorbed, similar as is done in LPJ-GUESS for the calculation of gross primary productivity (GPP).The storage pool m (Eq.4) is implemented as a single pool that reflects long-term changes.Introduction

Conclusions References
Tables Figures

Back Close
Full LPJ-GUESS with the two algorithms for monoterpene emission incorporated was evaluated against measurements for a Ponderosa pine (Pinus ponderosa L.) plantation at Blodgett forest, California (38.90 • N, 120.63 • W, elevation 1315 m, Holzinger et al., 2006).Monoterpene emissions were measured between June 2003 and April 2004 using proton-transfer-reaction mass-spectrometry in combination with the eddy covariance method (see Holzinger et al., 2006, for a detailed description of the measurements).Simulations were performed by applying LPJ-GUESS in gap-model mode, averaging 100 repeated calculations for a patch, which is necessary to account for the stochastic nature that is characteristic for some of the processes that underlie vegetation dynamics.To reproduce the plantation's uniform age, seedlings were established in the simulation year representing 1990, and the density was reduced in the simulation year representing 2000 to represent a thinning.The model was spun up with the monthly climate reanalysis produced by the Climatic Research Unit of the University of East Anglia (CRU, New et al., 2000) for the period 1990-2003, corrected with the anomaly between site climate and CRU climate.

Adjustments for global scale modelling
LPJ-GUESS can be applied on a global scale as a DGVM (Sitch et al., 2003).Compared to the gap-mode, it is based on simplified vegetation dynamics of the growth and development of an average individual.The vegetation is represented by ten plant functional types.Within such a global framework, a standard fraction of the electron flux used for monoterpene production s (Eq. 3) needs to be assigned as an average value to each PFT.Similar as in Arneth et al. (2007) for isoprene, s is set in a way that the emission rate obtained at standard conditions (303 K and 1000 µmol m 2 s 1 PAR) equals the emission capacity M s .These emission capacities are applied as average value to each PFT (e.g.Naik et al., 2004;Lathi ère et al., 2006), based on recommendations given in Guenther et al. (1995).
Values for the standard emissions for monoterpenes are highly uncertain, probably even more so than for isoprene: Leaf measurements, that are the basis for recommended values of M s , were mostly analysed using the temperature dependence in Eq. ( 1) which implies that emissions take place both day and night, independent of light conditions (and hence independent of the occurrence or absence of monoterpene production in the chloroplast) or of storage pool size (which may vary seasonally).For emissions from storage, the assumption of a constant emission factor is applicable for a short period only.Due to changes in the concentrations of monoterpenes in the storage organs, partial pressures will differ throughout the year, and will thereby influence volatilization and emission, and hence the measured M s .Several studies on monoterpenes report changes in emission capacities for different dates or seasons (e.g.Staudt et al., 1997;Komenda and Koppmann, 2002;Pressley et al., 2004;Hakola et al., 2006), which might be partially assigned to this storage effect.The emission rate at standard temperature, as reflected in the value of M s , is thus not equivalent to the actual production under standard conditions, which is, in all plants, a light-dependent process independent of direct release or storage after production.Therefore, in plants where such Introduction

Conclusions References
Tables Figures

Back Close
Full a storage pool exists, to maintain this pool over a period of one day or longer, production during daylight hours must be sufficient to support release from storage over 24 h.Or in other words: a production-derived value for M s must exceed an emission-derived M s notably.Of interest in this context is that emission capacities for monoterpeneemitting broadleaf species that do not store monoterpenes, taken from studies that applied a temperature and light dependent algorithm, range from 2 to 70 µg g −1 h −1 (Table 1), and are in general on the upper end of reported monoterpene emissions when compared to emissions that are released solely temperature-driven from storage (Kesselmeier and Staudt, 1999).These observations may indeed indicate a larger rate of production of monoterpenes taking place also in species where this production rate cannot be observed directly, because of the storage acting as a buffer between production and emission.
In the model, leaf production of monoterpenes is similar for all plants, independent of presence or absence of storage.For the application of storage (Eq.4) the produced monoterpenes can be transferred into storage, depending on the plant functional type (PFT).To do so, the group of PFTs was separated as in Table 2: all broadleaved trees are considered to be non-storing, while the conifers and herbs are considered to be storing monoterpenes.Such a simplification is unavoidable in DGVMs and is based on current observations: an isoprene-like release of monoterpenes was mostly found in broadleaved species, whereas conifers tend to have monoterpene storage (Kesselmeier and Staudt, 1999).Large amount of stored monoterpenes are also observed in many herbacious species.
Two global simulations were performed using the vegetation dynamics of LPJ-GUESS in DVGM mode, applying both the short-term temperature-dependent algorithm, and the production and storage algorithm.The CRU climate data for the years 1901-2000 were used to force the model with a spinup of 300 years, using the atmospheric CO 2 concentration for 1901 (290 ppm) and a detrended series of data for 1901-1950.This was followed by 100 years of simulation representing the 20th century, using the CRU data and CO 2 concentrations from ice cores and from observations.Introduction

Conclusions References
Tables Figures

Back Close
Full Screen / Esc Printer-friendly Version
This assumption also underlies all global scale estimates to date (see overview in Arneth et al., 2008a).The second simulation calculates production of monoterpenes from electron transport (Eq.2).Monoterpenes are emitted either directly, or from a storage pool (conifers and herbs, Eq. 5).For this simulation, emission capacities M s (and thus the values of s ) are adjusted as described above with a factor 2 to reflect the difference between the production (taking place only during sunlight hours) and the need to refill the storage (with emissions taking place the entire day).

Simulated monoterpene emissions
A prerequisite for reliable simulation of BVOC fluxes with vegetation models is the reproduction of important growth characteristics.For the Blodgett forest site, simulated LAI (one-sided) was 2.6 for 2003 and 2.8 for 2004 (not shown), which compares well to the observed variation of 2 to 3 for the study period (Holzinger et al., 2006).Simulated gross primary productivity (GPP, Fig. 1b) also agreed well with observations (Misson et al., 2006), except for a short period in late July 2003, for which the simulated daily GPP was approximately reduced by 25% compared to the values derived from eddy flux data, likely due to an overestimation of drought stress in the model during that period.Measurements from other years indicate that Blodgett forest experiences drought stress during summer, but the extent is less than in other Ponderosa pine forests with comparable climatic circumstances (Panek, 2004;Misson et al., 2004).Introduction

Conclusions References
Tables Figures

Back Close
Full Simulations were performed with two different algorithms for monoterpene emissions: (1) using the temperature-dependent algorithm (Eq. 1 and 2) using the monoterpene production algorithm (Eq.2) combined with release of emissions from storage (Eq.5) with various time constants τ s : 0 days (equivalent to no storage, release of emissions directly after production), and storage time constants up to 160 days (Table 3).The emission capacities M s that gave the best fit with each of the two algorithms (Eq. 1 and Eqs.2-6 with various storage settings) are of the same order of magnitude (Table 3).Holzinger et al. (2006) report an emission capacity (using Eq. 1) and based on all-sided leaf area) of 1 µmol m −2 leaf h −1 (or 0.20 µg C g −1 h −1 ), about half of the optimised value for Eq. ( 1) in this study.Differences between the two studies could be caused by a difference in extrapolation from leaf level to canopy level, as well as by the applied leaf temperature correction in this study.The simulated seasonality in emissions (Fig. 1c, for clarity only a selection of the simulations summarized in Table 3 is displayed) was very similar in all simulations, irrespective of the presence or absence of storage.This is due to the fact that temperature is always a main contributor to the variability since temperature and radiation normally correlate well, and warm days also have large rates of electron flux and hence monoterpene production.Only at high values for τ s (≥40 d) the seasonal differences were considerably reduced (not shown).Without storage, the photosynthesis-dependent simulated emissions show a strong day-to-day variability, as was the case with GPP (Fig. 1b), due to the dependence of terpene production on photosynthetic processes.The observed monoterpene emission peaks due to rain events (Fig. 1a and c) as were observed before at the same site (Schade et al., 1999) were not captured by any of the simulation experiments.Simulated changes in seasonality were caused only by changes in weather conditions and changes in the size of the storage pool, we did not include an explicit change of seasonality of monoterpene production as is often suggested for isoprene production in relation to changes in isoprene synthase activity (Wiberley et al., 2005).
Simulations performed with the photosynthesis-dependent algorithm (Eq.2) com-Introduction

Conclusions References
Tables Figures

Back Close
Full bined with release from storage of monoterpenes (Eq.5) showed an interesting feature: The best agreement between the fitted parameterization and observations, determined from the average mean error (AME) and the root mean square error (RMSE) between the two, was obtained both with no storage at all and with high residence times in storage (τ s =80 d, Table 3).However, the ratio between simulated summer and winter emissions decreases with increasing residence times due to the delay caused by the storage.The observed summer to winter emissions ratio of 5.3 was reproduced in the simulations only at values for τ s ≥40 d.Although observations for important parts of the simulation period were absent, the reasonably good fits with both low (τ s =0 d) and high (τ s =80 d) residence times, combined with the need for long storage to fit the ratio between summer and winter emissions, merits the assumption that part of the produced monoterpenes might be emitted directly, whereas another part is stored for considerable times.Such a mixture of long-term storage and direct emission (or emission from short-term storage, which is ignored in this study) was suggested by Staudt et al. (1997Staudt et al. ( , 2000) ) for Pinus pinea, and was proposed in a more general manner by Kesselmeier and Staudt (1999).A simulation which took this assumption into account (half of the produced monoterpenes was simulated to be emitted directly, the other half was stored with a standard residence time τ s of 80 days) resulted in lower values for both AME and RMSE, with values close to those obtained with the short-term temperature-dependent algorithm, and in a ratio of summer and winter emissions close to the observed value (5.0 based on the days for which observations were available).Increasing the standardized residence time τ s caused the concentration of monoterpenes in the leaves to increase to values up to 300 µg C g −1 at a τ s of 160 days (not shown), and the maximum concentration to be delayed until later in the year compared to the simulations with smaller τ s .Storage also caused the day-to-day variability of emissions to decrease (Fig. 1c), acting as a buffer between production and emission, as is the case for non-specific storage (Niinemets et al., 2004).Observed concentrations of β-pinene in a Ponderosa pine forest in Oregon, US, ranged between 2.8×10 3 and 5.1×10 3 and µg C g −1 (Lerdau et al., 1994) for September and June, respectively, with emission Introduction

Conclusions References
Tables Figures

Back Close
Full Screen / Esc Printer-friendly Version Interactive Discussion rates of 0.2 and 1.1 µg C g −1 h −1 , which indicates a similar order of magnitude for the residence times as obtained in our simulations.
The simulated seasonality of monoterpene concentrations in the Ponderosa pine plantation for the applied split of the emissions in storage (half) and direct emissions (half) is shown in Fig. 1d.The peak in simulated leaf monoterpene concentrations was reached in autumn.For the range of time coefficients applied (Table 3), the peak in concentrations shifted from summer (τ s =2.5 d) to late autumn (τ s =160 d), and the concentrations increased with increasing τ s (not shown).Measurements of the seasonal cycle of monoterpene concentrations in other species show a wide variety of patterns: A pattern similar to the simulated one, with high concentrations in summer and autumn, was observed for Black spruce (Picea mariana) in Canada (Lerdau et al., 1997), but not so for Jack pine (Pinus banksiana) in Canada, where concentrations peaked in spring and autum (Lerdau et al., 1997).Measurements of terpene concentrations in several Mediterranean species indicated low concentrations in summer and high in winter due to higher emissions at high temperatures (Llusia et al., 2006).
Our results did not account for changes in leaf mass over the year, which would affect the maximum storage pool size and could account for some of the variation observed in the timing of peak values.However, there are likely other factors playing an important role in the timing of storage concentrations that are not considered in our vegetation model.For instance, B äck et al. (2005) were able to reproduce large spring emissions of monoterpenes in boreal Pinus sylvestris by incorporating photorespiration as a carbon source, although the link between terpenoid production and photorespiration is controversial (Hewitt et al., 1990;Pe ñuelas and Llusi à, 2002).It is also plausible that terpene synthesis could take place during winter, since conifers are able to assimilate, albeit at low rates, during warm winter periods (e.g.Suni et al., 2003).Introduction

Conclusions References
Tables Figures

Back Close
Full The sensitivities of the leaf emissions calculated with the new algorithm were tested by varying the temperature and radiation input data to the model, but keeping the other settings as in the experiments described above.Temperature changes between −5 K and +5 K and radiation changes between −100% and +100% were added over the whole year to the forcing of the models for all days in spinup and simulation period.The responses shown are thus reflecting those of the canopy, including effects on canopy properties (e.g.LAI).The long-term sensitivity to temperature, reflecting the sensitivity of monoterpene production, that resulted from this was compared with the short-term (instantaneous) sensitivities of vapourization-related emissions from storage, both for the "classical" algorithm in Eq. ( 1) and for the short-term release implemented in the new algorithm (Eqs.5 and 6).GPP varied roughly 40% over the prescribed 10 K range (Fig. 2c) with changes in LAI of 30% (not shown).GPP showed a gradual increase up to its maximum at 3 K above ambient values, and a gradual decay from 3 K onwards, reflecting increasingly enhanced evaporation and stomatal closure as temperature increases in response to a soil water deficit (not shown).Due to the higher temperature optimum of terpene production combined with the relatively small changes in electron flux (and GPP) caused by the temperature difference, the long-term temperature sensitivity of monoterpene emissions (reflecting the sensitivity of production) is dominated by an exponential increase at low temperature, which levels off slightly at higher temperature due to the decreasing electron flux related to the decay in GPP (Fig. 2a).Relative changes in monoterpene production were much higher than for GPP for the same range of temperature changes.Due to the temperature optimum for GPP, relatively lower rates of GPP can coincide with either low or high rates of monoterpene emissions, depending on the temperature (Fig. 2b).
Next to the response of monoterpene production, Fig. 2a illustrates the short-term response of emissions from storage, both from Eq. ( 1) and from Eqs. ( 5) and ( 6).Both Introduction

Conclusions References
Tables Figures

Back Close
Full short-term sensitivities show an exponential increase with temperature; the difference in steepness of the two short-term curves is the result of the difference in Q 10 values that result from the two methods.The short-term sensitivities were slightly smaller than the long-term response, although the reaction to more extreme changes differs: because of the exponential nature of the short-term functions, it tends to cause larger peaks in emissions for short periods with high temperatures, whereas the long-term sensitivity showed a levelling off at high temperatures.This difference between production and emission at high temperatures was suggested as well by Kesselmeier and Staudt (1999).
Emission response to prescribed changes in radiation are shown in Fig. 3. GPP shows a logistic increase with increasing light levels from 25% of the current level onwards (Fig. 3c), below that radiation level photosynthesis over the year is too low to sustain the vegetation.LAI varied between 0 and 5.9 for the prescribed range in radiation levels (not shown).The logistic increase in GPP reflects the expected saturation for light, that is dominating larger parts of the year with increasing light levels.
For monoterpene production, the relation is more linear, and the relative changes in monoterpene production were larger than those in GPP (Fig. 3a).This was caused by a small heating of the leaf by enhanced radiation levels, causing the temperature dependence discussed above to play a minor indirect role as well.

Implications for global simulations of monoterpene emissions
For the global simulations, the determination of the standard fraction of the electron flux s from the usually reported standard emission capacities M s emphasises some critical uncertainties in understanding the processes that determine seasonal and annual emission patterns.The difference between an emission-based and a production-based value of M s , as discussed in Sect.2.4, was estimated to be approximately a factor of two, which reflects the difference between daylight hours and 24 h as well as the additional light limitation of the production and thus emissions.Incidentally, measured Introduction

Conclusions References
Tables Figures

Back Close
Full Screen / Esc Printer-friendly Version Interactive Discussion emission capacities of broadleaved trees where emissions are light-dependent, and the standardized rates hence represent the monoterpene production, also tend to be substantial (Table 1), supporting the view that the leaf production of monoterpenes during daylight hours is larger than seen when emissions are measured from storage pool release.M s was therefore doubled compared to the simulations using Eq. ( 1).
Annual global total terrestrial emissions were 29.6 Tg C a −1 for the simulation that assumed monoterpenes to be uniformly emitted from storage (Eq.1), and 31.8Tg C a −1 for the simulation that was based on production and storage, and that accounted for the frequently observed emissions without storage in broadleaved vegetation (Eq.2).The spatial distribution of the emissions is surprisingly similar in the two cases (Fig. 4).The production and storage algorithm resulted in larger rates in temperate forest regions in the eastern US, southern Brazil and China.In these areas, the applied correction factor of two is apparently too large compared to the actual reduction by the light dependence.
In dry regions in subtropical Africa, Northern India and Australia, where temperatures are high but photosynthesis rates are relatively low, the temperature-dependent algorithm resulted in larger rates.
Our estimates of global annual monoterpene emissions are at the low end of the published global totals.Naik et al. (2004), using the temperature dependence (Guenther et al., 1995) algorithm, reported 33 Tg C a −1 , which is comparable to our estimate with the same algorithm.These two experiments are comparable in their experimental design as well: both use potential natural vegetation cover, with similar tree PFTs in both models while Naik et al. (2004) simulated two additional shrub PFTs.However, these estimates are a factor of four lower than the highest published estimates (Guenther et al. (1995) report 127 Tg C a −1 , Lathi ère et al. ( 2006) report 117 Tg C a −1 ).This emphasises the large uncertainty in global BVOC emission calculations that can be introduced through use of different basal rates, vegetation cover and phenology, climatology, temporal resolution, and the use of different algorithms (Arneth et al., 2008a).
For the global simulations, storage was applied for the coniferous and herbaceous plant functional types ( applying a standard residence time τ s of 80 days.In the absence of long-term changes in the storage pool size, the parameterization of the storage equations (Eqs.5 and 6) affects only the the seasonality of the emissions, but not the annual totals.However, the seasonality of emissions is an important feature of monoterpene emission simulation when it comes to linking these to atmospheric chemistry.Application of monoterpene storage in the model with the parameters as derived in the local simulations caused significant residence times (averaged for all PFTs) mainly at high latitudes (Fig. 5).Large latitudinal differences between simulations with and without storage occurred in spring and autumn at high latitudes.During spring, when environmental conditions allow the onset of photosynthesis and monoterpene production, the storage pool is being built up, thereby moving part of the production into this storage pool and reducing the emitted amount.Moving from pole to equator, the difference between the simulations with and without storage is diminishing due to higher temperatures and the relatively larger contribution of directly emitting PFTs (Table 2).

Conclusions
We present here an analysis of monoterpene emissions that seeks to investigate the effects of two important processes separately, namely the production in the chloroplast and the ensuing emissions that may or may not be from a storage pool.The analysis aims to provide a basis for better understanding of observed seasonal patterns as well as to take into the account the increasing evidence of a direct, production-driven emission pattern in broadleaved vegetation.
The short-term sensitivities to temperature changes for both algorithms were comparable, but also the short-term sensitivity (on volatilization) and the long-term sensitivity (on production) were shown to be remarkably similar, at least as long as small changes in temperature are considered.We did not focus here on how the different monoterpene emission algorithms would be affected in simulations that take into account future climate change.It would seem that algorithms that include solely a response to increas-Introduction

Conclusions References
Tables Figures

Back Close
Full ing temperature would be more sensitive under future warming scenarios compared to those that also include a light-limitation, but the overall effects of climate change on other important processes like changes in leaf area index or vegetation cover would also need to be considered.What is more, it is uncertain whether the response of monoterpene production to increasing atmospheric CO 2 concentration follows a similar inhibitory pattern as is shown for isoprene in an increasing number of plants (Constable et al., 1999;Staudt et al., 2001;Baraldi et al., 2004), although the similarity in the chloroplastic pathways would suggest a similar response.
It is a general problem of BVOC emission modelling that parameterizations of algorithms that seek to represent observed constraints on emissions can only be based on a very limited number of studies and that true process understanding is often lacking (Guenther et al., 2006;Arneth et al., 2008a).Accounting in a global model for entire plant functional types to have either similar storage residence time or release monoterpenes directly is an inevitable necessity, but it cannot do justice to the natural variation.While most coniferous species that have been studied to date release monoterpenes mostly from storage, there are nonetheless species that emit part of their monoterpenes light-dependently (e.g.Pinus pinea, Staudt et al., 1997Staudt et al., , 2000)).At the same time, some broadleaf monoterpene emitters may also include storage organs (e.g.emissions from Eucalyptus spp.have been shown to depend primarily on temperatures, He et al., 2000).New DGVM model developments that -at least on continental scale -are capable of representing actual tree species distribution, rather than PFTs, can be used to assess the uncertainties associated with these globally applied simplified assumptions (e.g.Arneth et al., 2008b).Additionally, the distinction between monoterpene-storing and non-monoterpene-storing plant functional types has the potential to be extended to monoterpene types that are stored or non-stored.Several studies (Staudt et al., 1997(Staudt et al., , 2000) ) have shown a clear distinction between monoterpenes that are stored and thus have mainly temperature-dependent emission (e.g.limonene, α-pinene), and monoterpenes with emissions that react more directly in response to diurnal patterns of light or to shading, and that do not exhibit long-term Introduction

Conclusions References
Tables Figures

Back Close
Full storage (e.g.trans-β-ocimene, 1,8-cineole) in Pinus pinea.Eventually, such a model setup could also provide the basis for describing emissions that can occur in response to physical damage (e.g. by wind, rain or herbivores, Banchio et al., 2005;Pichersky et al., 2006).Such an analysis would present an important step forward on regional scale when seasonal emission rates are used in atmospheric chemistry simulations.
in cellular liquid or as specific storage in storage organs, Screen / Esc Printer-friendly Version Interactive Discussion response is B äck et al. (2005) for vapourization of α-pinene from the liquid phase.
The spinup was followed by a simulation with observed daily climate data (temperature, precipitation, radiation) at this site for the period from June 2003 until April 2004.The annual atmospheric CO 2 concentration was prescribed following global observations for the spinup and simulation periods.A set of simulations was performed to study the applicability of temperaturedependent (Eq. 1) and photosynthesis-and storage-dependent (Eqs. 2 and 5) algorithms to reproduce the observed emissions at Blodgett forest.The parameterization of the release from storage (τ s in Eq. 6) was varied to determine the value of best fit to the data.The emission capacities (and thus the standard fraction of the electron flux s , Eq. 3) were determined from the measurements such that the simulated emissions reproduced the annual average measured emissions on the days of the measurements.The (one-sided) specific leaf area for Ponderosa pine was prescribed to 7.8 m 2 kg −1 C followingMisson et al. (2005), and the thickness of the model's soil layers was increased to prevent an overestimation of water stress during the growing Introduction Screen / EscPrinter-friendly Version Interactive Discussion to changes in temperature and light

Fig. 2 .Fig. 3 .
Fig. 2. Sensitivity of GPP and monoterpene emissions to temperature changes: (a) Longterm change (in blue) of monoterpene emissions following Eq.(2) (reflecting the sensitivity of monoterpene production), compared to the short-term sensitivities (reflecting those of instantenous emissions from storage) of the "classical" temperature-dependent release from Guenther et al. (1993) (Eq. 1, in red) and of the storage release implemented in the photosynthesisdependent algorithm (Eqs.5 and 6, in orange); (c) Change in GPP with temperature; (b) Resulting relation between GPP and monoterpene emissions.

Fig. 5 .
Fig. 5. Annual cycle of the average residence time (in days) of the monoterpenes in the storage pool, shown are zonal means for the period 1981-2000.All PFTs (including the non-storing PFTs) are weighted according to their leaf area index in order to calculate latitudinal averages.

Table 2 )
, with half of the produced monoterpenes being stored, Introduction

Table 2 .
Naik et al. (2004)e of long-term monoterpene storage organs (+ indicates monoterpene storing PFT, − indicates non-storing), and emission capacities for the plant functional types used for the global simulation as adopted fromNaik et al. (2004).

Table 3 .
Results from simulations for Blodgett forest: scaled emission capacity M s , average mean error (AME), root mean square error (RMSE), ratio between summer (JJA) and winter (DJF) emissions for all days (for days with observation available in brackets, n=38 for summer, n=18 for winter).