the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 16 Apr 2024
Modelling of atmospheric concentrations of fungal spores: a two-year simulation over France using CHIMERE
Abstract. Fungal spore organic aerosol emissions have been recognised as a significant source of particulate matter as PM10; however, they are not widely considered in current air quality models. In this work, we have implemented the parametrisation of fungal spore organic aerosol (OA) emissions introduced by Heald and Spracklen (2009) (H&S) and further modified by Hoose et al. (2010) in the CHIMERE regional chemistry-transport model. This simple parametrisation is based on two variables which are leaf area index (LAI) and specific humidity. We have validated the geographical and temporal representativeness of this parametrisation on a large scale by using yearly polyol observations and primary biogenic organic aerosol factors from PMF analysis at 11 French measurement sites. For a group of sites in northern and eastern France, the seasonal variation of fungal spore emissions, displaying large summer and small winter values, is correctly depicted. However, the H&S parametrisation fails to capture fungal spore concentrations for a smaller group of Mediterranean sites with less data availability both in terms of absolute values as well as seasonal variability, leading to strong negative biases especially during the autumn and winter seasons occur. Two years of CHIMERE simulations with the H&S parametrisation have shown a significant contribution of fungal spore OA to PM10 mass, lower than 10 % during winter, and reaching up to 20 % during summer in high emission zones, especially over large forested areas. In terms of contribution to organic matter (OM) concentrations, the simulated fungal spore contribution in autumn is as high as 40 % and reaches at most 30 % of OM for other seasons. As a conclusion, the fungal spore OA contribution to total OM concentrations is shown to be substantial enough to be considered as a major PM10 fraction and shall then be included in state-of-the-art chemistry transport models. The H&S parametrisation shows satisfactory results over northern and eastern France, but may underestimate concentrations for Mediterranean areas that may indicate missing factors influencing emissions or a missing source of spores.
-
Notice on discussion status
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
-
Preprint
(1705 KB)
-
Supplement
(976 KB)
-
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
- Preprint
(1705 KB) - Metadata XML
-
Supplement
(976 KB) - BibTeX
- EndNote
- Final revised paper
Journal article(s) based on this preprint
Interactive discussion
Status: closed
-
RC1: 'Comment on egusphere-2024-698', Anonymous Referee #1, 23 Apr 2024
This manuscript explores fungal spore concentrations over France. The paper describes efforts: (1) to implement an existing scheme for fungal spores into the regional CHIMERE model (2) to evaluate the resulting simulation with 2 years of surface observations from 11 diverse sites in France and (3) to estimate the importance of this source to total OA and PM10. The objective and methods are straight-forward and the results are generally clearly described. I provide suggestions largely to clarify the methods and results.
- Line 43: suggest replace “shall” with “should”
- Lines 44-47: This abstract text is repetitive, and I recommend it be cut.
- Line 51: Inappropriate references for the text. The two references use AMS observations which measure particles less than ~1 um, not PM10. Pai et al. do not show an underestimate of in situ observations.
- Line 90: typo? Remove “the variability”
- Lines 113-114: The number of data points seems like excessive detail that bogs down the introduction
- Line 179-182: PMF is also commonly used for source apportionment of online measurements (such as AMS). Suggest you re-phrase.
- Table 1: Provide units or more details in caption for the column “PMF” and “Polyols”. Are these number of samples?
- Line 231: The reference Fountoukis and Nenes (2007) refers to ISORROPIA II – either correct the text to refer to this version or cite the original reference for ISORROPIA I
- Section 2.2.1: should also describe the simulation of other PM10 components (e.g. BC, dust, sea salt) since PM10 simulations are included in the manuscript.
- Table 2: This provides extraneous or duplicative information from the text. Suggest cut.
- Line 271: reference for this sentence?
- Equation 1: why are there multiple constants (c, 5, 1.5x10-2) rather than integrating these into a single constant?
- Line 291: please include units for density
- Lines 290-292: What size of fungal spore is assumed for this conversion? Are the spores treated as a size distribution or monodisperse?
- Line 298: are the fungal spores treated as soluble?
- Section 3.1 discussion of lifetime: Lifetime depends on the assumed size and solubility of the aerosol. These model assumptions should be described (see previous comments) and the lifetime diagnosed for the specific simulation. In particular, line 341 “Despite little transport” has not been shown. Are there any spores aloft in the simulation? Any transported to the eastern boundary of the domain?
- Line 445: It seems odd to end the discussion of Figure 4 without discussing the outliers Aix and Nice. This comes later in the text, so the authors may wish to explicitly tell the reader that they will discuss these outliers in Section 3.2.2, or move this discussion forward.
- Line 460: why do the authors use MFE here when using MFB in Figure 4. I suggest harmonizing for consistent statistical metrics.
- Figure 5 and lines 539-546: Please include all sites with long-term measurements discussed in the text in Figure 5.
Citation: https://doi.org/10.5194/egusphere-2024-698-RC1 -
CC1: 'Reply on RC1', Matthieu Vida, 28 May 2024
Publisher’s note: this comment is a copy of AC2 and its content was therefore removed.
Citation: https://doi.org/10.5194/egusphere-2024-698-CC1 -
AC2: 'Reply on RC1', Gilles Foret, 01 Jul 2024
Thank you so much for the detailed review work to improve the clarity and understanding of the paper. Responses to all your comments can be found below, together with the amended paper containing the track changes.
• Line 43: suggest replace “shall” with “should”
--> done
• Lines 44-47: This abstract text is repetitive, and I recommend it be cut.--> We have cut the sentence, it is now:
“The H&S parameterisation shows satisfactory results over northern and eastern France, but may underestimate concentrations for Mediterranean areas that may indicate missing factors influencing emissions or a missing source of spores.”
• Line 51: Inappropriate references for the text. The two references use AMS observations which measure particles less than ~1 um, not PM10. Pai et al. do not show an underestimate of in situ observations.--> We have replaced the sentence “Modelling of the organic matter (OM) fraction of PM10 chronically underestimates in situ observations (Ciarelli et al., 2016; Pai et al., 2020). This underestimation can be attributed to several causes such as the complexity of the organic matter composition, which is not yet fully known, incomplete emission inventories or their inherent uncertainties, and poorly parametrised atmospheric chemical transformations.” by “Modelling of the organic matter (OM) fraction of PM10 with chemical transport models can be complex due to the varied composition of organic matter, which is not yet fully known, incomplete emission inventories or their inherent uncertainties, and poorly parameterised atmospheric chemical transformations.”
• Line 90: typo? Remove “the variability”--> It was a typo. We have changed the sentence
“amaké et al. (2019a) identified the parameters responsible the variability for up to 82 % of the annual variability of polyols […]” by “Samaké et al. (2019a) identified the parameters responsible for up to 82 % of the annual variability of polyols […]”
• Lines 113-114: The number of data points seems like excessive detail that bogs down the introduction.--> We have replaced the sentence “This comparison was carried out using 1,536 hourly data points, 1,200 of which came from the German and Finnish stations, each with 600 data points.”by “This comparison was carried out using 1,536 hourly data points, that most of which came from the German (600) and Finnish (600) stations.”
• Lines 179-182: PMF is also commonly used for source apportionment of online measurements (such as AMS). Suggest you re-phrase.--> We have replaced the sentence “Positive Matrix Factorisation (PMF) is one of the most widely used techniques for identifying factors contributing to aerosol concentrations using field measurements” by “Positive Matrix Factorisation (PMF) is one of the most widely used techniques for identifying factors contributing to aerosol concentrations using online and offline measurement data”
• Table 1: Provide units or more details in caption for the column “PMF” and “Polyols”. Are these number of samples?--> We have replaced the caption “Summary of organic matter and polyols filter-based observations as well as primary biogenic factor […]” by “Summary of the number of daily filters analysed for polyols as well as OM from primary biogenic factor […]”
• Line 231: The reference Fountoukis and Nenes (2007) refers to ISORROPIA II – either correct the text to refer to this version or cite the original reference for ISORROPIA I.--> We have replaced the text “The ISORROPIA thermodynamic model […]” by “The ISORROPIA II thermodynamic model […]”
• Section 2.1: should also describe the simulation of other PM10 components (e.g. BC, dust, sea salt) since PM10 simulations are included in the manuscript.--> A description of the simulation of the chemical composition of PM10 and its evaluation is provided in a separate paper. For spores, no interaction with other species was taken into account. This is why we have chosen not to go into detail on the simulation of all the other species.
• Table 2: This provides extraneous or duplicative information from the text. Suggest cut.--> We have cut the Table 2
• Line 271: reference for this sentence?--> We have cut the sentence, it is now “This yields to large emissions especially over Southern Europe which are not confirmed by measurements.”
• Equation 1: why are there multiple constants (c, 5, 1.5x10-2) rather than integrating these into a single constant?--> We have rewritten the equation from Hoose et al. (2010). We have added the simplified version “30 867 × LAI × qv”.
• Line 291: please include units for density.--> We have modified the sentence “Fungal spore number concentrations are transformed into mass using an aerosol density of 1 […]” by “Fungal spore number concentrations are transformed into mass using an aerosol relative density of 1 kg kg-1 […]”
• Lines 290-292: What size of fungal spore is assumed for this conversion? Are the spores treated as a size distribution or monodisperse?--> We have added the sentence “Fungal spores are treated as non-soluble particles and considered as monodispersed with a diameter of 3 µm in a model size class closest to 3 µm. In the model configuration with 10 size bins, spores are included in the size bin 8 corresponding to sizes between 2.5 and 5 µm, with an average diameter of 3.5 µm.”
• Line 298: are the fungal spores treated as soluble?--> No, this is now mentioned in the sentence added in response to the previous comment.
• Section 3.1 discussion of lifetime: Lifetime depends on the assumed size and solubility of the aerosol. These model assumptions should be described (see previous comments) and the lifetime diagnosed for the specific simulation. In particular, line 341 “Despite little transport” has not been shown. Are there any spores aloft in the simulation? Any transported to the eastern boundary of the domain?--> There are indeed concentrations of spores at higher altitudes (> 20 metres) but obviously much lower than in the first 20 metres. On the other hand, the transport at the eastern boundary of the domain is negligible. As mentioned in the text, deposition is fairly rapid, given the large size of the particles. Hummel et al (2015) calculated a lifetime of about 5h for biogenic fluorescent aerosols (FBAP). If we compare the simulated fungal spore emission and concentration maps, we see a very strong similarity indicating a short lifetime.
We have added the sentence “Since only the transport and deposition of fungal spores are taken into account, without interactions with other species, the conditions are similar to those of Hummel et al. (2015), with a lifetime of the same order of magnitude (5 hours), further confirming the low transport of spores.” and changed the sentence “Despite little transport” by “Despite these conditions of transport and deposition”
• Line 445: It seems odd to end the discussion of Figure 4 without discussing the outliers Aix and Nice. This comes later in the text, so the authors may wish to explicitly tell the reader that they will discuss these outliers in Section 3.2.2, or move this discussion forward.--> We have added the sentence “The complicated results observed at stations on the Mediterranean coast (Aix-en-Provence, Marseille, Nice, Port-de-Bouc) will be discussed in the next section.”
• Line 460: why do the authors use MFE here when using MFB in Figure 4. I suggest harmonizing for consistent statistical metrics.--> Both statistical metrics are different and used in parallel. We have replaced the sentence “Finally, for comparison with the same polyol data set, daily and monthly mean fractional error (MFE) are respectively 0.79 and 0.56 at all sites” by “Finally, for comparison with the same polyol data set, daily and monthly mean fractional bias (MFB) are respectively -11 % and -11 % at all sites (Table S1, Figure S1). The root mean square error (RMSE) and the mean fractional error (MFE) was also calculated for estimating the error (Table S1, Figure S2). Daily and monthly MFE are respectively 79 % and 56 % at all sites; for the median RMSE the results are respectively 0.04 µg m-3 and 0.03 µg m-3.”
• Figure 5 and lines 539-546: Please include all sites with long-term measurements discussed in the text in Figure 5.--> We have included in Figure 5 all the timeseries for non-Mediterranean sites. The timeseries for OMpb and polyols for all sites are in the supplement.
Citation: https://doi.org/10.5194/egusphere-2024-698-AC2
-
RC2: 'Comment on egusphere-2024-698', Anonymous Referee #2, 07 May 2024
The manuscript describes and evaluates the implementation of fungal spores in a chemistry transport model (CTM) over France. The study is interesting and pioneering in the CTMs and shows that fungal spores can account for a large fraction of PM10 and OM in the area. The manuscript is well-written and is easy to follow. I have the below comments that should be addressed before it can be published in ACP.
Comments:
Line 172: Out of 699 stations in France, fungal spore tracers are observed in 13. If the full station network is not used to evaluate model results in the present study, I find it confusing to read about the 699 stations in the beginning.
Line 188: Are these 7 sites from the 13 sites in total used in this study? And why only these 7?
Lines 190-195: is there an estimate of how significant this contribution can be?
Table 1: It would be great if you could provide the altitude information.
Figure 1: Is it possible to have the map showing the land use types and denote the measurement type, ie polyol vs PMF?
Line 243: Could you comment on the 9x9 km spatial resolution representing the traffic stations.
Lines 274-276: Would be good if you can show the statistics of these 3 different parameterizations (H&S, Hummel, and Janssen) in the supplement? It would be interesting to see the impact of considering wind as used in Janssen, in relation to the hypothesis in lines 323-325. The correlations may be close but how about the magnitudes?
Line 286: So, the fungal spores are emitted in only one size (3 mm)? It would be good to discuss this assumption, rather than using a size distribution?
Line 330: Could you infer the lifetime through your simulations? It would be interesting to compare with Hummel.
Line 435: Does Chimere calculate OM or OC? If the latter, how is the OC converted to OM?
Figure 4: Please explain in the caption what the box plot and the symbols represent.
The poorest model performance seems to be on the urban sites (traffic and industrial). Is this due to model resolution? Or the meteorological model not performing well for the urban areas?
Lines 495-499: Is this because of meteorology? Winds, precipitation, etc? It would be good to write in the methodology if these are nudged simulations. There might be an issue with transport, which could be identified by evaluating the model winds (speed and direction) that can explain the large underestimation in the southern sites?
Citation: https://doi.org/10.5194/egusphere-2024-698-RC2 -
CC2: 'Reply on RC2', Matthieu Vida, 28 May 2024
Publisher’s note: this comment is a copy of AC1 and its content was therefore removed.
Citation: https://doi.org/10.5194/egusphere-2024-698-CC2 -
AC1: 'Reply on RC2', Gilles Foret, 01 Jul 2024
Thank you so much for your careful review. I am honoured that you found the study interesting and that the paper is for you well written and easy to follow. Responses to all your comments can be found below, together with the amended paper containing the track changes.
• Line 172: Out of 699 stations in France, fungal spore tracers are observed in 13. If the full station network is not used to evaluate model results in the present study, I find it confusing to read about the 699 stations in the beginning.--> Initially this paper included an assessment of PM10 and we forgot to delete this sentence. In the end, it was moved to a separate paper that will be submitted in the next few days. We have removed the sentence “A total of 699 air quality stations performed measurements in metropolitan France during the period of the study, restricted to 2013 and 2014, including fixed and mobile stations.”
• Line 188: Are these 7 sites from the 13 sites in total used in this study? And why only these 7?--> As mentioned in Table 1, the sites with OM data for the primary biogenic PMF factor are among those with polyol measurements. However, not all sites were subject to a PMF study since the complete chemical composition was not carried out. This explains why not all the sites measuring polyols have PMF data. The same applies to the number of PMF vs. polyol data, which can be explained by the fact that the chemical composition of PM10 was not carried out on all the filters.
• Lines 190-195: is there an estimate of how significant this contribution can be?--> Yes, there are estimations, and we modified the following sentence to include one. We have replaced the sentence “PMF results at all sites include a factor which can be attributed to PBOA because of the large concentrations of the two polyols in this factor, representing more than 90 % of them.” by “PMF results at all sites include a factor which can be attributed to PBOA because of the large concentrations of the two polyols in this factor, representing more than 90 % of the polyols total mass in this factor (Samaké et al., 2019a).”
• Table 1: It would be great if you could provide the altitude information.--> We have included altitude information in Table 1.
• Figure 1: Is it possible to have the map showing the land use types and denote the measurement type, ie polyol vs PMF?--> We have included the number of data for each site in Figure 1 respectively for OMpb and polyols.
• Line 243: Could you comment on the 9x9 km spatial resolution representing the traffic stations.--> The aim of this study is to look at the concentrations of primary biogenic aerosols. Since no interaction with other species is taken into account and these aerosols are emitted by vegetation, resolution at traffic sites is likely not an issue there. Furthermore, our results show good correlations and low biases at traffic sites for fungal spores. There could be an impact on resuspension, but this is not taken into account in the model for organic aerosols within the model.
• Lines 274-276: Would be good if you can show the statistics of these 3 different parameterizations (H&S, Hummel, and Janssen) in the supplement? It would be interesting to see the impact of considering wind as used in Janssen, in relation to the hypothesis in lines 323-325. The correlations may be close but how about the magnitudes?--> Thank you for this suggestion. As these results are already published in the literature and not part of our work (we rely on them), we refrain from reproducing them in our work, but they are easily accessible in the open scientific literature. Please note that the evaluation of the H&S and Hummel parametrisations are performed with the same biogenic fluorescent aerosol data (Hummel et al., 2015) and that the Janssen parametrisation has been evaluated alone in a separate paper on a spore dataset (Janssen et al., 2021). The differences in magnitudes can be highlighted with the discussion of biases.
• Line 286: So, the fungal spores are emitted in only one size (3 mm)? It would be good to discuss this assumption, rather than using a size distribution?--> We have added the sentence “Fungal spores are treated as non-soluble particles and considered as monodispersed with a diameter of 3 µm in a model size class closest to 3 µm. In the model configuration with 10 size bins, spores are included in the size bin 8 corresponding to sizes between 2.5 and 5 µm, with an average diameter of 3.5 µm.”
• Line 330: Could you infer the lifetime through your simulations? It would be interesting to compare with Hummel.--> Calculating a spore lifetime is complex. We only used the deposition and transport processes on the spores like in Hummel et al. (2015), so we assume that the spore lifetime should be similar. As mentioned in the text, deposition is fairly rapid, given the large size of the particles. Hummel et al (2015) calculated a lifetime of about 5h for biogenic fluorescent aerosols (FBAP). If we compare the simulated fungal spore emission and concentration maps, we see a very strong similarity indicating a short lifetime.
• Line 435: Does Chimere calculate OM or OC? If the latter, how is the OC converted to OM?--> Chimere simulates OM. In section 2.1.2 we have added the sentence “The organic carbon of the primary biogenic PMF was multiplied by a factor of 1.8 to obtain the organic matter concentrations of this PMF factor (OMpb) (Favez et al., 2010; Petit et al., 2015).”
• Figure 4: Please explain in the caption what the box plot and the symbols represent.--> We have added the sentence “Correlations and mean fractional bias (MFB) boxplots are obtained using the scores for each site (1 point per site) and are illustrated respectively at left and right side. The points corresponding to the sites represent the scores obtained on all the data for the site.”
• The poorest model performance seems to be on the urban sites (traffic and industrial). Is this due to model resolution? Or the meteorological model not performing well for the urban areas?--> It's not a question of the type of site, but of the location. The most complicated sites to model are those close to the Mediterranean. This may be due to the resolution of the model, but above all to the different climate and possibly sources, which can have an impact on spore concentrations. This is already largely discussed in the paper and summarised at the end of the discussions section (section 3.3).
• Lines 495-499: Is this because of meteorology? Winds, precipitation, etc? It would be good to write in the methodology if these are nudged simulations. There might be an issue with transport, which could be identified by evaluating the model winds (speed and direction) that can explain the large underestimation in the southern sites?--> The meteorological data come from the WRF model, which was run in parallel with the CHIMERE model for the simulations and indeed nudged. The resolution of WRF meteorological data are consistent with the resolution chosen for the model, i.e. 9 x 9 km² in our study. For the 9 x 9 km² simulations, the domain and the simulations are not the result of nesting; the 9 x 9 km² European domain is our only domain and our only simulation. We have identified the meteorology as a source of problems in the Mediterranean area, with night-time temperatures thought to play a major role.
We have replaced the sentence “The WRF model is used for meteorological forcing (Skamarock et al., 2008).” by “The WRF 3.7.1 model (Skamarock et al., 2008) is used for meteorological simulation coupled to chimere with no aerosol effect. The spectral nudging (Von Storch et al., 2000) is used with NCEP temperature, wind, humidity, pressure for wavelengths higher than 2 000 km whereas in the boundary layer the WRF model freely generates its own dynamic.”
Citation: https://doi.org/10.5194/egusphere-2024-698-AC1
-
CC2: 'Reply on RC2', Matthieu Vida, 28 May 2024
Interactive discussion
Status: closed
-
RC1: 'Comment on egusphere-2024-698', Anonymous Referee #1, 23 Apr 2024
This manuscript explores fungal spore concentrations over France. The paper describes efforts: (1) to implement an existing scheme for fungal spores into the regional CHIMERE model (2) to evaluate the resulting simulation with 2 years of surface observations from 11 diverse sites in France and (3) to estimate the importance of this source to total OA and PM10. The objective and methods are straight-forward and the results are generally clearly described. I provide suggestions largely to clarify the methods and results.
- Line 43: suggest replace “shall” with “should”
- Lines 44-47: This abstract text is repetitive, and I recommend it be cut.
- Line 51: Inappropriate references for the text. The two references use AMS observations which measure particles less than ~1 um, not PM10. Pai et al. do not show an underestimate of in situ observations.
- Line 90: typo? Remove “the variability”
- Lines 113-114: The number of data points seems like excessive detail that bogs down the introduction
- Line 179-182: PMF is also commonly used for source apportionment of online measurements (such as AMS). Suggest you re-phrase.
- Table 1: Provide units or more details in caption for the column “PMF” and “Polyols”. Are these number of samples?
- Line 231: The reference Fountoukis and Nenes (2007) refers to ISORROPIA II – either correct the text to refer to this version or cite the original reference for ISORROPIA I
- Section 2.2.1: should also describe the simulation of other PM10 components (e.g. BC, dust, sea salt) since PM10 simulations are included in the manuscript.
- Table 2: This provides extraneous or duplicative information from the text. Suggest cut.
- Line 271: reference for this sentence?
- Equation 1: why are there multiple constants (c, 5, 1.5x10-2) rather than integrating these into a single constant?
- Line 291: please include units for density
- Lines 290-292: What size of fungal spore is assumed for this conversion? Are the spores treated as a size distribution or monodisperse?
- Line 298: are the fungal spores treated as soluble?
- Section 3.1 discussion of lifetime: Lifetime depends on the assumed size and solubility of the aerosol. These model assumptions should be described (see previous comments) and the lifetime diagnosed for the specific simulation. In particular, line 341 “Despite little transport” has not been shown. Are there any spores aloft in the simulation? Any transported to the eastern boundary of the domain?
- Line 445: It seems odd to end the discussion of Figure 4 without discussing the outliers Aix and Nice. This comes later in the text, so the authors may wish to explicitly tell the reader that they will discuss these outliers in Section 3.2.2, or move this discussion forward.
- Line 460: why do the authors use MFE here when using MFB in Figure 4. I suggest harmonizing for consistent statistical metrics.
- Figure 5 and lines 539-546: Please include all sites with long-term measurements discussed in the text in Figure 5.
Citation: https://doi.org/10.5194/egusphere-2024-698-RC1 -
CC1: 'Reply on RC1', Matthieu Vida, 28 May 2024
Publisher’s note: this comment is a copy of AC2 and its content was therefore removed.
Citation: https://doi.org/10.5194/egusphere-2024-698-CC1 -
AC2: 'Reply on RC1', Gilles Foret, 01 Jul 2024
Thank you so much for the detailed review work to improve the clarity and understanding of the paper. Responses to all your comments can be found below, together with the amended paper containing the track changes.
• Line 43: suggest replace “shall” with “should”
--> done
• Lines 44-47: This abstract text is repetitive, and I recommend it be cut.--> We have cut the sentence, it is now:
“The H&S parameterisation shows satisfactory results over northern and eastern France, but may underestimate concentrations for Mediterranean areas that may indicate missing factors influencing emissions or a missing source of spores.”
• Line 51: Inappropriate references for the text. The two references use AMS observations which measure particles less than ~1 um, not PM10. Pai et al. do not show an underestimate of in situ observations.--> We have replaced the sentence “Modelling of the organic matter (OM) fraction of PM10 chronically underestimates in situ observations (Ciarelli et al., 2016; Pai et al., 2020). This underestimation can be attributed to several causes such as the complexity of the organic matter composition, which is not yet fully known, incomplete emission inventories or their inherent uncertainties, and poorly parametrised atmospheric chemical transformations.” by “Modelling of the organic matter (OM) fraction of PM10 with chemical transport models can be complex due to the varied composition of organic matter, which is not yet fully known, incomplete emission inventories or their inherent uncertainties, and poorly parameterised atmospheric chemical transformations.”
• Line 90: typo? Remove “the variability”--> It was a typo. We have changed the sentence
“amaké et al. (2019a) identified the parameters responsible the variability for up to 82 % of the annual variability of polyols […]” by “Samaké et al. (2019a) identified the parameters responsible for up to 82 % of the annual variability of polyols […]”
• Lines 113-114: The number of data points seems like excessive detail that bogs down the introduction.--> We have replaced the sentence “This comparison was carried out using 1,536 hourly data points, 1,200 of which came from the German and Finnish stations, each with 600 data points.”by “This comparison was carried out using 1,536 hourly data points, that most of which came from the German (600) and Finnish (600) stations.”
• Lines 179-182: PMF is also commonly used for source apportionment of online measurements (such as AMS). Suggest you re-phrase.--> We have replaced the sentence “Positive Matrix Factorisation (PMF) is one of the most widely used techniques for identifying factors contributing to aerosol concentrations using field measurements” by “Positive Matrix Factorisation (PMF) is one of the most widely used techniques for identifying factors contributing to aerosol concentrations using online and offline measurement data”
• Table 1: Provide units or more details in caption for the column “PMF” and “Polyols”. Are these number of samples?--> We have replaced the caption “Summary of organic matter and polyols filter-based observations as well as primary biogenic factor […]” by “Summary of the number of daily filters analysed for polyols as well as OM from primary biogenic factor […]”
• Line 231: The reference Fountoukis and Nenes (2007) refers to ISORROPIA II – either correct the text to refer to this version or cite the original reference for ISORROPIA I.--> We have replaced the text “The ISORROPIA thermodynamic model […]” by “The ISORROPIA II thermodynamic model […]”
• Section 2.1: should also describe the simulation of other PM10 components (e.g. BC, dust, sea salt) since PM10 simulations are included in the manuscript.--> A description of the simulation of the chemical composition of PM10 and its evaluation is provided in a separate paper. For spores, no interaction with other species was taken into account. This is why we have chosen not to go into detail on the simulation of all the other species.
• Table 2: This provides extraneous or duplicative information from the text. Suggest cut.--> We have cut the Table 2
• Line 271: reference for this sentence?--> We have cut the sentence, it is now “This yields to large emissions especially over Southern Europe which are not confirmed by measurements.”
• Equation 1: why are there multiple constants (c, 5, 1.5x10-2) rather than integrating these into a single constant?--> We have rewritten the equation from Hoose et al. (2010). We have added the simplified version “30 867 × LAI × qv”.
• Line 291: please include units for density.--> We have modified the sentence “Fungal spore number concentrations are transformed into mass using an aerosol density of 1 […]” by “Fungal spore number concentrations are transformed into mass using an aerosol relative density of 1 kg kg-1 […]”
• Lines 290-292: What size of fungal spore is assumed for this conversion? Are the spores treated as a size distribution or monodisperse?--> We have added the sentence “Fungal spores are treated as non-soluble particles and considered as monodispersed with a diameter of 3 µm in a model size class closest to 3 µm. In the model configuration with 10 size bins, spores are included in the size bin 8 corresponding to sizes between 2.5 and 5 µm, with an average diameter of 3.5 µm.”
• Line 298: are the fungal spores treated as soluble?--> No, this is now mentioned in the sentence added in response to the previous comment.
• Section 3.1 discussion of lifetime: Lifetime depends on the assumed size and solubility of the aerosol. These model assumptions should be described (see previous comments) and the lifetime diagnosed for the specific simulation. In particular, line 341 “Despite little transport” has not been shown. Are there any spores aloft in the simulation? Any transported to the eastern boundary of the domain?--> There are indeed concentrations of spores at higher altitudes (> 20 metres) but obviously much lower than in the first 20 metres. On the other hand, the transport at the eastern boundary of the domain is negligible. As mentioned in the text, deposition is fairly rapid, given the large size of the particles. Hummel et al (2015) calculated a lifetime of about 5h for biogenic fluorescent aerosols (FBAP). If we compare the simulated fungal spore emission and concentration maps, we see a very strong similarity indicating a short lifetime.
We have added the sentence “Since only the transport and deposition of fungal spores are taken into account, without interactions with other species, the conditions are similar to those of Hummel et al. (2015), with a lifetime of the same order of magnitude (5 hours), further confirming the low transport of spores.” and changed the sentence “Despite little transport” by “Despite these conditions of transport and deposition”
• Line 445: It seems odd to end the discussion of Figure 4 without discussing the outliers Aix and Nice. This comes later in the text, so the authors may wish to explicitly tell the reader that they will discuss these outliers in Section 3.2.2, or move this discussion forward.--> We have added the sentence “The complicated results observed at stations on the Mediterranean coast (Aix-en-Provence, Marseille, Nice, Port-de-Bouc) will be discussed in the next section.”
• Line 460: why do the authors use MFE here when using MFB in Figure 4. I suggest harmonizing for consistent statistical metrics.--> Both statistical metrics are different and used in parallel. We have replaced the sentence “Finally, for comparison with the same polyol data set, daily and monthly mean fractional error (MFE) are respectively 0.79 and 0.56 at all sites” by “Finally, for comparison with the same polyol data set, daily and monthly mean fractional bias (MFB) are respectively -11 % and -11 % at all sites (Table S1, Figure S1). The root mean square error (RMSE) and the mean fractional error (MFE) was also calculated for estimating the error (Table S1, Figure S2). Daily and monthly MFE are respectively 79 % and 56 % at all sites; for the median RMSE the results are respectively 0.04 µg m-3 and 0.03 µg m-3.”
• Figure 5 and lines 539-546: Please include all sites with long-term measurements discussed in the text in Figure 5.--> We have included in Figure 5 all the timeseries for non-Mediterranean sites. The timeseries for OMpb and polyols for all sites are in the supplement.
Citation: https://doi.org/10.5194/egusphere-2024-698-AC2
-
RC2: 'Comment on egusphere-2024-698', Anonymous Referee #2, 07 May 2024
The manuscript describes and evaluates the implementation of fungal spores in a chemistry transport model (CTM) over France. The study is interesting and pioneering in the CTMs and shows that fungal spores can account for a large fraction of PM10 and OM in the area. The manuscript is well-written and is easy to follow. I have the below comments that should be addressed before it can be published in ACP.
Comments:
Line 172: Out of 699 stations in France, fungal spore tracers are observed in 13. If the full station network is not used to evaluate model results in the present study, I find it confusing to read about the 699 stations in the beginning.
Line 188: Are these 7 sites from the 13 sites in total used in this study? And why only these 7?
Lines 190-195: is there an estimate of how significant this contribution can be?
Table 1: It would be great if you could provide the altitude information.
Figure 1: Is it possible to have the map showing the land use types and denote the measurement type, ie polyol vs PMF?
Line 243: Could you comment on the 9x9 km spatial resolution representing the traffic stations.
Lines 274-276: Would be good if you can show the statistics of these 3 different parameterizations (H&S, Hummel, and Janssen) in the supplement? It would be interesting to see the impact of considering wind as used in Janssen, in relation to the hypothesis in lines 323-325. The correlations may be close but how about the magnitudes?
Line 286: So, the fungal spores are emitted in only one size (3 mm)? It would be good to discuss this assumption, rather than using a size distribution?
Line 330: Could you infer the lifetime through your simulations? It would be interesting to compare with Hummel.
Line 435: Does Chimere calculate OM or OC? If the latter, how is the OC converted to OM?
Figure 4: Please explain in the caption what the box plot and the symbols represent.
The poorest model performance seems to be on the urban sites (traffic and industrial). Is this due to model resolution? Or the meteorological model not performing well for the urban areas?
Lines 495-499: Is this because of meteorology? Winds, precipitation, etc? It would be good to write in the methodology if these are nudged simulations. There might be an issue with transport, which could be identified by evaluating the model winds (speed and direction) that can explain the large underestimation in the southern sites?
Citation: https://doi.org/10.5194/egusphere-2024-698-RC2 -
CC2: 'Reply on RC2', Matthieu Vida, 28 May 2024
Publisher’s note: this comment is a copy of AC1 and its content was therefore removed.
Citation: https://doi.org/10.5194/egusphere-2024-698-CC2 -
AC1: 'Reply on RC2', Gilles Foret, 01 Jul 2024
Thank you so much for your careful review. I am honoured that you found the study interesting and that the paper is for you well written and easy to follow. Responses to all your comments can be found below, together with the amended paper containing the track changes.
• Line 172: Out of 699 stations in France, fungal spore tracers are observed in 13. If the full station network is not used to evaluate model results in the present study, I find it confusing to read about the 699 stations in the beginning.--> Initially this paper included an assessment of PM10 and we forgot to delete this sentence. In the end, it was moved to a separate paper that will be submitted in the next few days. We have removed the sentence “A total of 699 air quality stations performed measurements in metropolitan France during the period of the study, restricted to 2013 and 2014, including fixed and mobile stations.”
• Line 188: Are these 7 sites from the 13 sites in total used in this study? And why only these 7?--> As mentioned in Table 1, the sites with OM data for the primary biogenic PMF factor are among those with polyol measurements. However, not all sites were subject to a PMF study since the complete chemical composition was not carried out. This explains why not all the sites measuring polyols have PMF data. The same applies to the number of PMF vs. polyol data, which can be explained by the fact that the chemical composition of PM10 was not carried out on all the filters.
• Lines 190-195: is there an estimate of how significant this contribution can be?--> Yes, there are estimations, and we modified the following sentence to include one. We have replaced the sentence “PMF results at all sites include a factor which can be attributed to PBOA because of the large concentrations of the two polyols in this factor, representing more than 90 % of them.” by “PMF results at all sites include a factor which can be attributed to PBOA because of the large concentrations of the two polyols in this factor, representing more than 90 % of the polyols total mass in this factor (Samaké et al., 2019a).”
• Table 1: It would be great if you could provide the altitude information.--> We have included altitude information in Table 1.
• Figure 1: Is it possible to have the map showing the land use types and denote the measurement type, ie polyol vs PMF?--> We have included the number of data for each site in Figure 1 respectively for OMpb and polyols.
• Line 243: Could you comment on the 9x9 km spatial resolution representing the traffic stations.--> The aim of this study is to look at the concentrations of primary biogenic aerosols. Since no interaction with other species is taken into account and these aerosols are emitted by vegetation, resolution at traffic sites is likely not an issue there. Furthermore, our results show good correlations and low biases at traffic sites for fungal spores. There could be an impact on resuspension, but this is not taken into account in the model for organic aerosols within the model.
• Lines 274-276: Would be good if you can show the statistics of these 3 different parameterizations (H&S, Hummel, and Janssen) in the supplement? It would be interesting to see the impact of considering wind as used in Janssen, in relation to the hypothesis in lines 323-325. The correlations may be close but how about the magnitudes?--> Thank you for this suggestion. As these results are already published in the literature and not part of our work (we rely on them), we refrain from reproducing them in our work, but they are easily accessible in the open scientific literature. Please note that the evaluation of the H&S and Hummel parametrisations are performed with the same biogenic fluorescent aerosol data (Hummel et al., 2015) and that the Janssen parametrisation has been evaluated alone in a separate paper on a spore dataset (Janssen et al., 2021). The differences in magnitudes can be highlighted with the discussion of biases.
• Line 286: So, the fungal spores are emitted in only one size (3 mm)? It would be good to discuss this assumption, rather than using a size distribution?--> We have added the sentence “Fungal spores are treated as non-soluble particles and considered as monodispersed with a diameter of 3 µm in a model size class closest to 3 µm. In the model configuration with 10 size bins, spores are included in the size bin 8 corresponding to sizes between 2.5 and 5 µm, with an average diameter of 3.5 µm.”
• Line 330: Could you infer the lifetime through your simulations? It would be interesting to compare with Hummel.--> Calculating a spore lifetime is complex. We only used the deposition and transport processes on the spores like in Hummel et al. (2015), so we assume that the spore lifetime should be similar. As mentioned in the text, deposition is fairly rapid, given the large size of the particles. Hummel et al (2015) calculated a lifetime of about 5h for biogenic fluorescent aerosols (FBAP). If we compare the simulated fungal spore emission and concentration maps, we see a very strong similarity indicating a short lifetime.
• Line 435: Does Chimere calculate OM or OC? If the latter, how is the OC converted to OM?--> Chimere simulates OM. In section 2.1.2 we have added the sentence “The organic carbon of the primary biogenic PMF was multiplied by a factor of 1.8 to obtain the organic matter concentrations of this PMF factor (OMpb) (Favez et al., 2010; Petit et al., 2015).”
• Figure 4: Please explain in the caption what the box plot and the symbols represent.--> We have added the sentence “Correlations and mean fractional bias (MFB) boxplots are obtained using the scores for each site (1 point per site) and are illustrated respectively at left and right side. The points corresponding to the sites represent the scores obtained on all the data for the site.”
• The poorest model performance seems to be on the urban sites (traffic and industrial). Is this due to model resolution? Or the meteorological model not performing well for the urban areas?--> It's not a question of the type of site, but of the location. The most complicated sites to model are those close to the Mediterranean. This may be due to the resolution of the model, but above all to the different climate and possibly sources, which can have an impact on spore concentrations. This is already largely discussed in the paper and summarised at the end of the discussions section (section 3.3).
• Lines 495-499: Is this because of meteorology? Winds, precipitation, etc? It would be good to write in the methodology if these are nudged simulations. There might be an issue with transport, which could be identified by evaluating the model winds (speed and direction) that can explain the large underestimation in the southern sites?--> The meteorological data come from the WRF model, which was run in parallel with the CHIMERE model for the simulations and indeed nudged. The resolution of WRF meteorological data are consistent with the resolution chosen for the model, i.e. 9 x 9 km² in our study. For the 9 x 9 km² simulations, the domain and the simulations are not the result of nesting; the 9 x 9 km² European domain is our only domain and our only simulation. We have identified the meteorology as a source of problems in the Mediterranean area, with night-time temperatures thought to play a major role.
We have replaced the sentence “The WRF model is used for meteorological forcing (Skamarock et al., 2008).” by “The WRF 3.7.1 model (Skamarock et al., 2008) is used for meteorological simulation coupled to chimere with no aerosol effect. The spectral nudging (Von Storch et al., 2000) is used with NCEP temperature, wind, humidity, pressure for wavelengths higher than 2 000 km whereas in the boundary layer the WRF model freely generates its own dynamic.”
Citation: https://doi.org/10.5194/egusphere-2024-698-AC1
-
CC2: 'Reply on RC2', Matthieu Vida, 28 May 2024
Peer review completion
Journal article(s) based on this preprint
Viewed
| HTML | XML | Total | Supplement | BibTeX | EndNote | |
|---|---|---|---|---|---|---|
| 406 | 119 | 36 | 561 | 52 | 23 | 24 |
- HTML: 406
- PDF: 119
- XML: 36
- Total: 561
- Supplement: 52
- BibTeX: 23
- EndNote: 24
Viewed (geographical distribution)
| Country | # | Views | % |
|---|---|---|---|
| United States of America | 1 | 195 | 34 |
| France | 2 | 105 | 18 |
| China | 3 | 37 | 6 |
| Germany | 4 | 31 | 5 |
| United Kingdom | 5 | 31 | 5 |
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
- 195
Matthieu Vida
Gilles Foret
Guillaume Siour
Florian Couvidat
Olivier Favez
Arineh Cholakian
Sébastien Conil
Matthias Beekmann
Jean-Luc Jaffrezo
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
- Preprint
(1705 KB) - Metadata XML
-
Supplement
(976 KB) - BibTeX
- EndNote
- Final revised paper