the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Growth and actual leaf temperature modulate CO2 responsiveness of monoterpene emissions from holm oak in opposite ways
Michael Staudt
Juliane Daussy
Joseph Ingabire
Nafissa Dehimeche
Download
- Final revised paper (published on 26 Oct 2022)
- Supplement to the final revised paper
- Preprint (discussion started on 13 Jul 2022)
- Supplement to the preprint
Interactive discussion
Status: closed
-
RC1: 'Comment on bg-2022-142', Anonymous Referee #1, 02 Aug 2022
Growth and actual leaf temperature modulate CO2 -responsiveness of monoterpene emissions from Holm oak in opposite ways
Staudt et al.
General Comment:
The authors did a very thorough investigation on a specific scientific question that certainly is of relevance for the evaluation of climate change impacts on biogenic emissions and feedbacks on air chemistry. In my opinion the experiment has been well set up and carried out. The interpretation is supported by a number of ancillary measurements so that some interesting ideas about the potential underlying mechanisms could be developed. Also, the authors revealed a well-founded knowledge about the topic and the relevant literature.
On the downside, I noticed that wording and style could be improved. Many sentences are inconveniently complicated or long and selected expression are often unfamiliar or imprecise. I would recommend to check, shorten, and involve an English native to improve the text. Also some shifts between results and discussion sections and a better description of the equations used for sensitivity analysis should be considered at the appropriate places.
Specific Comments:
Abstract
It seems unclear to me, what the cool and warm growth regimes look like. Indicating only the 5-degree difference is not sufficient. Compared with the quite extensively discussed results and conclusion, the description of the outcome is relatively meager.
Introduction
L65: I assume that MTs are not synthetized but only stored in resin ducts.
L76: superfluous ‘very’ (remove)
L85: superfluous ‘before’ (remove)
Description
There is a bit of a mix between description and discussion, check (e.g. L200-203)
Could you please indicate the equation used for emission factor reduction in MEGAN here (and not in the results as a caption text)?
Define G400, A400
Results
Figure 1: It is a bit irritating that the emission factor (per unit m2) should increase with the number of leaves. I see that the latter is meant as a growth indicator, which should, however, be better illustrated (e.g. final number of leaves? Number of leaves in the end of the growth period?)
Figure 2: Better use the same design for Ci in each of the graphs (i.e. that which shows relative NPQ)
Figure 3: You probably mean key relations instead of key correlations. Actually, I have difficulties to see understand both, the explanations of how this is calculated and the reason why it has been done.
L364-366: The difference between the explanations for the two different responses to temperature are unclear. Rephrase and consider to elaborate the arguments.
L371ff: Should this really be one figure caption? Generally, I expect a short, clear and consistent description of what I see. This is violated at least since line 376. Instead, take care that the abbreviations are all clear (e.g. chloro, growth?). It could also be considered to use this figure as a basis for discussion and put into chapter 4, possibly in several stages in order to better support the reasoning in the different chapter.
Discussion
What I am missing is a discussion in how far the results can be assumed general findings or are specific for Quercus ilex? Is it likely that conifers, evergreens, broadleaves or Mediterranean plants react similar? Do you think the BVOC emission groups should then be differentiated by their degree of genetic relatedness or to site conditions typical for the species?
L513ff: With the summary here, the paragraph tends to be lengthy and repetitive. I would suggest to take the essence from this paragraph to the conclusions (and delete it here).
L550ff: Here, for the first time if I am not mistaken, the authors declare that they also run some simulations to test the sensitivity of the found mechanisms. While I am not against such exercises, this comes as a surprise and should have been mentioned and described before (and shorten it here). Also Fig. 6 is a result and only part of its description belongs into discussions.
Conclusion
L599: concentrations instead of variations; “hardly effect emissions” or “affect emissions only marginally” or similar instead of “affect little emissions”. (good example for wrong wording)
L615ff: Missing knowledge as well as stating additional references is not something, that should be put into a conclusion. Please consider to shift it towards the discussion.
Citation: https://doi.org/10.5194/bg-2022-142-RC1 -
AC1: 'Reply on RC1', Michael Staudt, 05 Aug 2022
Dear referee 1
First of all, I wish to thank you for the fast and thorough revision of our manuscript. I am aware that the manuscript is quite long and a bit complicated, hence not easy to review. I appreciate and greatly acknowledge all comments made by both referees that will help me to improve the manuscript. In the following I will just respond to the points in order to keep the discussion going (answers are in italics).
Growth and actual leaf temperature modulate CO2 -responsiveness of monoterpene emissions from Holm oak in opposite ways
Staudt et al.
General Comment:
The authors did a very thorough investigation on a specific scientific question that certainly is of relevance for the evaluation of climate change impacts on biogenic emissions and feedbacks on air chemistry. In my opinion the experiment has been well set up and carried out. The interpretation is supported by a number of ancillary measurements so that some interesting ideas about the potential underlying mechanisms could be developed. Also, the authors revealed a well-founded knowledge about the topic and the relevant literature.
Answer: I am very pleased about these very positive and encouraging comments.
On the downside, I noticed that wording and style could be improved. Many sentences are inconveniently complicated or long and selected expression are often unfamiliar or imprecise. I would recommend to check, shorten, and involve an English native to improve the text.
Answer: I apologize for my poor English. Unfortunately, I do not have a native speaker available for language proofreading. I suggest that the final manuscript version (if accepted) be reviewed by a professional language editor. Perhaps Copernicus and the associate editor can advise me on this. I will also carefully re-examine the text myself, for example, by doing multiple forward and backward translations with free Internet translators. I usually use DeepL, but perhaps there are better ones.
Also some shifts between results and discussion sections and a better description of the equations used for sensitivity analysis should be considered at the appropriate places.
Answer: The MEGAN equation will be placed in the main text.
Specific Comments:
Abstract
It seems unclear to me, what the cool and warm growth regimes look like. Indicating only the 5-degree difference is not sufficient. Compared with the quite extensively discussed results and conclusion, the description of the outcome is relatively meager.
Answer: I understand that the description of the results in the ABSTRACT may seem meagre compared to the relatively long DISCUSSION. However, in my opinion, it contains all the results presented in the RESULTS and the DISCUSSION. I also believe that an abstract should be concise and avoid speculative, very far-reaching conclusions, especially if they must include an additional contextual introduction. However, if the referee feels that a particular result of the study is missing in the abstract, I will try to include it.
Regarding the temperature-growth regimes, I will add details on the day/night temperatures that were applied.
Introduction
L65: I assume that MTs are not synthetized but only stored in resin ducts.
Answer: The terpenes in the resin ducts are synthesised in the glandular epithelium surrounding the cavity of the resin duct into which they are secreted (see e.g. https://nph.onlinelibrary.wiley.com/doi/full/10.1111/nph.15984 and references therein). Generally, the synthesis is particularly intense during resin duct development. Hence stored resin MTs are not synthesized in the photosynthetic parenchyma (photosynthetic source tissues) of leafs/needles and then transported into the resin ducts. However, MT synthesis in the resin duct epithelium may rely on photosynthates (essentially sucrose) provided by the photosynthetic parenchyma of source leaves (i.e. leaves that produce more photosynthates than they use for their own respiration and maintenance).
L76: superfluous ‘very’ (remove)
Answer: ok
L85: superfluous ‘before’ (remove)
Answer: ok
Description
There is a bit of a mix between description and discussion, check (e.g. L200-203)
Answer: Yes, I suggest removing the sentence “Previous studies showed that …”.
Could you please indicate the equation used for emission factor reduction in MEGAN here (and not in the results as a caption text)?
Answer: yes, it will be included in the section M&M
Define G400, A400
Answer: The definition of G and A is found in L130 and the meaning of the subscript number 400 is explained in L205. Briefly, a variable with the subscript 400 is the value that resulted from the first measurement of the CO2 response curve made at 400 ppm CO2, 1000 PPFD and at an assay temperature of 30°C or 35°C. Accordingly, in our study, the term "emission factor" corresponds to the temperature-normalized E400 values. I can add an explanatory example if the reviewer finds the explanation given on L205 unclear or insufficient.
Results
Figure 1: It is a bit irritating that the emission factor (per unit m2) should increase with the number of leaves. I see that the latter is meant as a growth indicator, which should, however, be better illustrated (e.g. final number of leaves? Number of leaves in the end of the growth period?)
Answer: Yes the total of leaves per plant is taken as a measure of the plant’s growth performance. To better understand this it should be noted that young potted QI plants in greenhouse culture do not exhibit a fixed period of leaf growth as it occurs in the field. Under field conditions QI tress show typically only one leaf flush lasting from late spring to early summer (onset of drought), though under certain circumstances there can be a second one in the same year either from the buds formed in the spring and/or from dormant buds. In our experiment, the QI saplings kept under none-stress conditions (well watered, no extreme temperatures) continued more or less growing in repeated cycles or even in indeterminate growth manner (central apex) until the end of the experiment. The acorns were potted at the same time and the number of leaves (plus few other morphological features) were determined at the end of the experiment in September. Also, there was no apparent difference in leaf size among the four growth treatments with only moderate differences in specific leaf weight (LMA). These facts allowed us to consider the number of leaves as a proxi for foliage growth. The plants were not immediately harvested after the experiment, because we wanted keeping them alive for eventual additional measurements. Finally these were not made, due to the lack of time, manpower and because the plants had to leave rapidly the greenhouse compartments. Regarding the results: Growth at elevated CO2 had a fertilizing effect on leaf growth, while growth under the warmer conditions had rather a negative effect. However, individual plants in each population differed appreciably in plant size and leaf mass, even though they were always maintained under the same growth conditions. This variability was positively related to the emission factor, A400 and ETR400, measured on a single leaf of each plant. This observation is (although not completely novel), in my opinion, one of the most interesting findings that deserved to be considered in the discussion (L400 ff).
Figure 2: Better use the same design for Ci in each of the graphs (i.e. that which shows relative NPQ)
Answer: Yes of course, all figure panels should show the same Ci scaling. This bug is probably due to a copy-paste error of the Excel figures that we overlooked during the final check of the manuscript and for which I apologize.
Figure 3: You probably mean key relations instead of key correlations. Actually, I have difficulties to see understand both, the explanations of how this is calculated and the reason why it has been done.
Answer: Yes key relations might be more appropriate term than key correlations. Pearson analysis is the analysis of the linear relationship between two quantitative variables. The result is given as Pearson correlation coefficient R, which ranges between -1 and +1 and provides information about the direction and strength of the linear relationship, or as determination coefficients R² (0-1), which is a measure of goodness of fit explaining the proportion of variance explained by the model (linear relationship in case of Pearson). I consider Pearson correlation analyses as a simple mean to check relationships between the key variables of interest (that are 3 in the present study: emission factor, relative emission at low CO2, relative emission at high CO2) to other variables, thus providing indications of the determinants of variations and mechanisms behind. For example at 35°C assay temperature we found that the emission response to low CO2 (E<400/E400) was positively correlated with the ETR response to low CO2 (ETR<400/ETR400, Fig. 3a) and with the leaf’s initial photosynthesis rate (A400, Fig 3b). The leaf’s initial photosynthesis depended much on the leaf’s initial stomatal opening G400 (R² between A400 & G400 = 0.924) and hence ETR<400/ETR400 also correlated with G400 (Fig.3c). However, neither A400 nor G400 correlated with ETR<400/ETR400, suggesting that the emissions response to low CO2 levels is determined by two independent factors (cf.), which could therefore together explain more than 80% of its variability (R2: 0.420 and 0.445; Figs 3a, b). I strongly prefer to present such results in scatter plots as in Figs. 1 and 3 because they show the input data in its original form and its distribution (possible presence of outliers, clusters, tendencies for non-linear relationships). However, when the Pearson analyses involve a large number of variables, as in the present study, two-dimensional scatter plots are less suitable for illustrating correlation networks in a clear and concise way. For this reason, I created the diagrams shown in Figure 5.
L364-366: The difference between the explanations for the two different responses to temperature are unclear. Rephrase and consider to elaborate the arguments.
Answer: I will revise this part.
L371ff: Should this really be one figure caption? Generally, I expect a short, clear and consistent description of what I see. This is violated at least since line 376. Instead, take care that the abbreviations are all clear (e.g. chloro, growth?). It could also be considered to use this figure as a basis for discussion and put into chapter 4, possibly in several stages in order to better support the reasoning in the different chapter.
Answer: This figure was thought to provide readers an overview and summary of the outcome of the Pearson correlation analyses. The caption is indeed very long and it possible to remove caption text from “The results can be summarized…” . I found it also difficult to find the right placement of this figure in the text. Referee 2 suggests removing this figure since it is hard to read. I may suggest moving this figure (after few corrections) to the supplement 2 near the corresponding correlation matrices (Table S3).
Discussion
What I am missing is a discussion in how far the results can be assumed general findings or are specific for Quercus ilex? Is it likely that conifers, evergreens, broadleaves or Mediterranean plants react similar? Do you think the BVOC emission groups should then be differentiated by their degree of genetic relatedness or to site conditions typical for the species?
Answer: This is indeed an important point with respect to emission modeling and inventories. As mentioned in the manuscript, our results show a strong similarity with isoprene emissions. Therefore, I might be tempted to conclude that all monoterpene emissions directly linked to their de-novo synthesis in photosynthetic tissues might behave similarly. However, I prefer to be very cautious about such generalizations. Even for isoprene emissions, considerable interspecies differences in CO2 responses have been observed. There is a recent paper by Niinemets et al. 2021 that specifically addresses this question (https://onlinelibrary.wiley.com/doi/abs/10.1111/pce.14131).
L513ff: With the summary here, the paragraph tends to be lengthy and repetitive. I would suggest to take the essence from this paragraph to the conclusions (and delete it here).
Answer: I will revise this part and see whether how it could be placed in the section CONCLUSIONS. This paragraph at the end of the discussion was thought to provide readers a brief summary of the multitude of metabolic responses that could have interacted and affected the emissions during CO2 ramping including others not mentioned before in the discussion. Generally, experimental studies on this topic are focusing on particular processes according their hypothesis (which is understandable) thus neglecting a bit the complexity that in reality exists.
L550ff: Here, for the first time if I am not mistaken, the authors declare that they also run some simulations to test the sensitivity of the found mechanisms. While I am not against such exercises, this comes as a surprise and should have been mentioned and described before (and shorten it here). Also Fig. 6 is a result and only part of its description belongs into discussions.
Answer: I fully understand that this part of the DISCUSSION is surprising after the previous part. Originally, when I wrote the first draft, I just wanted to make some general statements here about the degree of emission inhibition observed in our study and the known emission response to temperature (roughly what is written on LL 529-539). The temperature dependence of emission inhibition at high CO2 levels then led me to run simulations based on climate data collected at the flux tower of our forest station, combining different scenarios for maximum CO2 inhibition, warming and seasonality. In fact, on the annual scale, it is not easy to predict the extent to which a given emission inhibition could compensate for the emission increase due to X °C global warming, as many different factors interact in a non-linear way. I find the results of the simulations quite informative and have therefore retained this section in the paper in the hope of increasing its impact. It might be possible to move large parts of it to the end of the RESULTS as a small extra chapter.
Conclusion
L599: concentrations instead of variations; “hardly effect emissions” or “affect emissions only marginally” or similar instead of “affect little emissions”. (good example for wrong wording)
Answer: Thank you very much for the concrete examples, which will help me to improve the wording.
L615ff: Missing knowledge as well as stating additional references is not something, that should be put into a conclusion. Please consider to shift it towards the discussion.
Answer: The missing knowledge cited here is meant as a kind of outlook (what should be done next...?), which is a conclusion in a broader sense. Nevertheless, as aforementioned I will revise this part.
Thanks again for the helpful comments. I look forward to the next round of discussion and manuscript review.
With kind regards,
Michael Staudt
Citation: https://doi.org/10.5194/bg-2022-142-AC1 -
AC3: 'Reply on RC1', Michael Staudt, 31 Aug 2022
Dear Referees
This message is just a "corrigendum" to our first response to Referee 1's comments posted on August 5. When re-reading our first comments, we found nasty “typos” and a missing reference to the manuscript in our long answer to the referee 1 comment to Figure 3 (“Figure 3: You probably mean key relations instead of key correlations. Actually, I have difficulties to see understand both, the explanations of how this is calculated and the reason why it has been done.”)
In our answer, the following sentences should read correctly:
“For example at 35°C assay temperature we found that the emission response to low CO2 (E<400/E400) was positively correlated with the ETR response to low CO2 (ETR<400/ETR400, Fig. 3a) and with the leaf’s initial photosynthesis rate (A400, Fig 3b). The leaf’s initial photosynthesis depended much on the leaf’s initial stomatal opening G400 (R² between A400 & G400 = 0.924) and hence E<400/E400-1 also correlated with G400 (Fig.3c). However, neither A400 nor G400 correlated with ETR<400/ETR400-1, suggesting that the emissions response to low CO2 levels is determined by two independent factors (cf. L431 ff), which could therefore together explain more than 80% of its variability (R2: 0.420 and 0.445; Figs 3a, b).”
Due to a copy-paste error, the term "E<400/E400-1" was confused with the term "ETR<400/ETR400-1”, rendering our response meaningless. This error might be exemplary of a weakness of our manuscript that complicates its readability, namely the many abbreviations of ecophysiological variables and their derivatives. During the revision of our manuscript (if approved by the editor and reviewers), we will improve this by reducing, simplifying, and clarifying the terminology currently used in the manuscript (including the abstract). We will also add equations to better illustrate how the various calculations and simulations were performed. It is in our own interest to produce an article enjoyable to read in order to attract a broad readership of BIOGEOSCIENCES, many of whom may not be specialists of plant ecophysiology.
Michael Staudt, on behalf of the co-authors
Citation: https://doi.org/10.5194/bg-2022-142-AC3
-
AC1: 'Reply on RC1', Michael Staudt, 05 Aug 2022
-
RC2: 'Comment on bg-2022-142', Anonymous Referee #2, 02 Aug 2022
The manuscript describes a greenhouse experiment where the effects of elevated CO2 and growth temperature on holm oak leaf scale monoterpene emission rates are assessed. This is very relevant research topic already for decades, and the authors manage to scrutinize the experiment in a way that they can eventually conclude novel and interesting results.
The monoterpene emission responses to elevated CO2 and temperature were decoupled. Clear differences between cool- and warm-grown plants could be seen, the latter being more sensitive to CO2 inhibition. Contrasting this, a lower actual measurement temperature seemed to lead to larger CO2 inhibition compared to measurements at higher (35C) temperatures. This is rather surprising when the temperature difference is only 5C. The authors explain this with the leaf energy balance, similarly as has been shown for isoprene. Still, some explanations of the seemingly rather small temperature difference should be interesting for readers. In contrast, growth CO2 had no significant effect on emission CO2 sensitivity, although it promoted plant growth and the leaf’s emission factor.
The methods are well designed and elegantly used. Several different normalisation methods are used for assessing the uncertainties related to plant chemotype, growth conditions and measurement conditions. Finally, the obtained non-linear responses are used to upscale the short term impacts to annual emission dynamics using the MEGAN algorithm.
Overall, the ms represents an elegant experiment and is well compiled. It could be revised by removing some of the speculations and using the figures more directly to show the reader the main results, this would lead to significant shortening and clarifications of the main messages. Some linguistic errors and typos should be corrected, and a few other aspects could be clarified in the manuscript:
- how old were the measured leaves, were they of same age?
- what part of the canopy? how tall were the saplings?
- what was their rooting size?
- emission measurements: how many adsorbent tubes per CO2 level and leaf?
- was the humidity of incoming air controlled?
- Supplementary table 1 has remnants of non-english origin (mars)
- Figure 5 is an overview of the correlation network, but is does not really clarify the results and is almost impossible to read. I recommend removing it. However, I was missing a multivariate analysis where the combined effects of temperature and CO2 levels could have been assessed.
Citation: https://doi.org/10.5194/bg-2022-142-RC2 -
AC2: 'Reply on RC2', Michael Staudt, 05 Aug 2022
Dear referee 2
Thanks you very much for your compliments on our study and the efforts to carefully read and review our manuscript and the inspiring comments. I will soon start revising the manuscript. In the meantime, I would like to address some of the points you raised in your comments:
Referee 2 Comments:
Growth and actual leaf temperature modulate CO2 -responsiveness of monoterpene emissions from Holm oak in opposite ways
The manuscript describes a greenhouse experiment where the effects of elevated CO2 and growth temperature on holm oak leaf scale monoterpene emission rates are assessed. This is very relevant research topic already for decades, and the authors manage to scrutinize the experiment in a way that they can eventually conclude novel and interesting results.
The monoterpene emission responses to elevated CO2 and temperature were decoupled. Clear differences between cool- and warm-grown plants could be seen, the latter being more sensitive to CO2 inhibition. Contrasting this, a lower actual measurement temperature seemed to lead to larger CO2 inhibition compared to measurements at higher (35C) temperatures. This is rather surprising when the temperature difference is only 5C. The authors explain this with the leaf energy balance, similarly as has been shown for isoprene. Still, some explanations of the seemingly rather small temperature difference should be interesting for readers. In contrast, growth CO2 had no significant effect on emission CO2 sensitivity, although it promoted plant growth and the leaf’s emission factor.
The methods are well designed and elegantly used. Several different normalisation methods are used for assessing the uncertainties related to plant chemotype, growth conditions and measurement conditions. Finally, the obtained non-linear responses are used to upscale the short term impacts to annual emission dynamics using the MEGAN algorithm.
Overall, the ms represents an elegant experiment and is well compiled. It could be revised by removing some of the speculations and using the figures more directly to show the reader the main results, this would lead to significant shortening and clarifications of the main messages.
Answer: Thank you very much for these very positive words. I will consider all advices during my revision. However, it would be helpful for me if the reviewer could clarify the meaning of "the more direct use of figures..." and also which speculative conclusions she/he thinks should be removed.
Some linguistic errors and typos should be corrected, and a few other aspects could be clarified in the manuscript:
- how old were the measured leaves, were they of same age? what part of the canopy?
Answer: All leaves were mature leaves of the current year. The age of the leaves in months is not known. For practical reasons, leaves were selected that were not too small and were at the end of the shoots so that they could be accommodated in the LiCOR leaf chamber and covered the entire chamber surface.
- how tall were the saplings?
Answer: the height of the saplings ranged between 15 and 70 cm.
- what was their rooting size?
Answer: The size of roots or any measure requiring the harvest of the plants were (unfortunately) not made at the end of the experiment.
- emission measurements: how many adsorbent tubes per CO2 level and leaf?
Answer: Only one VOC sample per CO2 level so that the VOC sampling phase and consequently the duration of the entire CO2 ramping was not too long. Overall, we had very few sample losses due to errors in the GC-MS analysis. In contrast, we lost entire CO2 ramping series because the leaf was injured or even detached (the petiole of Q ilex leaves is very short).
- was the humidity of incoming air controlled?
Answer: The humidity of the incoming air was not controlled (H2O-scrubber was not used). Mean relative humidity during CO2-ramping was 43 +/- 5 and 32 +/- 6 at 30°C and 35°C respectively.
- Supplementary table 1 has remnants of non-english origin (mars)
Answer: thanks for this hint, I will correct it.
- Figure 5 is an overview of the correlation network, but is does not really clarify the results and is almost impossible to read. I recommend removing it. However, I was missing a multivariate analysis where the combined effects of temperature and CO2 levels could have been assessed.
Answer: I agree that the figure is not easy to read even though it shows only the main connections. However, besides showing the key correlations, it can help readers understanding the experimental protocol. Therefore, I suggest keeping the figure as supplementary material in supplement 2 (near the corresponding matrices of Pearson correlation). The scatter plots shown in Figs 1 and 3, which can be extended, plus more detailed descriptions of the Pearson results in the text should be sufficient. In addition, I will improve the coloration of the matrices towards a kind of heat map, thus providing a more easy readable overview of the correlation networks.
Regarding multivariate analyses I will check what I am able to do and whether the outputs are instructive.
Thank you again for your kind review. I look forward to your feedback.
With kind regards,
Michael Staudt
Citation: https://doi.org/10.5194/bg-2022-142-AC2