the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Implementation of the MOSAIC Aerosol Module (v1.0) in the Canadian Air Quality Model GEM-MACH (v3.1)
Abstract. The Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) aerosol thermodynamics and sectional framework has been implemented into the Canadian operational air quality model GEM-MACH. The original aerosol sub-model in GEM-MACH is based on the Canadian Aerosol Module (CAM), which uses a single-moment (mass) sectional scheme, and inorganic thermodynamics derived from the equilibrium ISORROPIA model without base metal cations. MOSAIC features non-equilibrium inorganic thermodynamics and a two-moment (mass and number) sectional scheme. For evaluation we conduct four one-year simulations with the same emissions and meteorology over the North America domain. A reference run (REF) with the Zhang et al. (2001) aerosol dry deposition scheme and a sensitivity run (EMR) with updated parameters from Emerson et al. (2020) is conducted for each aerosol model option. The results are compared to station observations and surface monthly-mean model-observation synthesis data. MOSAIC exhibits a shift in the accumulation mode mass and number distribution compared to CAM that results in more aerosol dry deposition in the REF run and a surface PM2.5 sulfate low bias of about 15 % relative to CAM. This bias is essentially removed in the MOSAIC EMR run resulting in a better fit to aggregated urban and rural stations compared to CAM over the North America domain. Comparison with the AERONET volume size distribution inversion product shows that MOSAIC gives a much higher level of agreement in terms of location of the accumulation mode peak diameter and separation of the accumulation and coarse modes. PM2.5 nitrate and ammonium for the MOSAIC EMR run shows overall better agreement with observation station data compared to both REF and EMR CAM runs at rural stations. At urban stations MOSAIC has a high bias for nitrate relative to CAM and observations during summer but it is reduced in the EMR run compared to the REF run. The high bias in ammonium seen with CAM for both REF and EMR runs relative to aggregated rural and urban station observations is reduced with MOSAIC by about 25 % between April and November.
- Preprint
(27122 KB) - Metadata XML
-
Supplement
(36574 KB) - BibTeX
- EndNote
Status: final response (author comments only)
-
RC1: 'Comment on egusphere-2024-2958', Peter Colarco, 17 Jan 2025
Two sectional aerosol schemes are tested in the GEM-MACH modeling framework. A relatively simpler scheme based on the Canadian Aerosol Module (CAM) is compared to a more complex scheme based on Model for Simulating Aerosol Interactions and Chemistry (MOSAIC). The CAM and MOSAIC configurations are both run for a year using the same size bin structure, emissions, and meteorology. Two simulations are performed for each configuration, a baseline using the default dry deposition scheme (called REF) and a sensitivity experiment that implements a newer dry deposition scheme (called EMR). Model results are compared to surface observations of component PM2.5 concentrations from networks over Canada and the USA, as well as retrieved columnar volume size distributions from AERONET sub photometers. Performance is evaluated in terms of sulfate, nitrate, ammonium, total particular mass, and aerosol water content.
My observation is that performance is more alike than dissimilar across the suite of simulations performed. Particularly for sulfate the results of the different configurations are very similar in terms of the seasonal cycles and spatial distributions. There is more diversity between the MOSAIC and CAM nitrate comparisons (although still highly correlated). Ammonium is somewhat intermediate. There is a notable improvement in the simulated nitrate surface concentration for the MOSAIC scheme with the updated (EMR) dry deposition scheme, which seems to give overall the best performance for total PM of nitrate (Figure 9). THE MOSAIC scheme additionally better agrees with the AERONET retrieved particle size distribution in terms of the placement of the fine mode and simulation of the observed gap between the fine and coarse modes. CAM by contrast has a larger fine mode that “smears” into the coarse mode with no gap. (Neither configuration represents the coarse mode especially well.) In part the differences are attributed to numerical diffusion (CAM does not conserve particle number explicitly and so numerical diffusion is more prevalent) and the CAM result is slightly improved by calling thermodynamics routines for each bin rather than in a bulk approach (at added computational expense).
It seems that the MOSAIC scheme with the EMR dry deposition produces overall the best set of results, albeit with a computational cost 3x the CAM run. Some improvements to the CAM configuration seem to be possible. A difference in the thermodynamics (no treatment of cations in CAM while they are treated in MOSAIC) are important to the nitrate simulation. There does not seem to be sufficient information given to judge the “best” approach. For example, if the goal is 5 day forecast of PM and it takes 1 hour to produce with CAM and 3 hours to produce with MOSAIC that seems insignificant. If it is 1 day versus 3 days then it matters. I can’t judge this given what was in the paper.
I suggest accepting the paper subject to minor revisions. The paper is thorough and well written with only a few errors noted. I do suggest some slight alteration of the presentation in terms of providing some quantitative information in tables.
Specific points:
Line 70 - what is the implication here of chemical and aerosol tracers not being transported by convection? This seems like it would be important to determining the vertical profile. Can you add anything about why this is not a first order issue to be addressed?
Line 213 - I am confused about what is being discussed here, noting particularly the claim to be using GEOS-5 driving data. Later on line 264 you say you are getting dynamical information from GDPS. Is something in error here?
Captions to figures 2b and 2c need to be corrected to refer to figure 2a (not 3a)
The numerical values in the labeling on Figure 5 are impossible to read even blown up on my screen. I suggest that you tabulate these either in the main text of the supplement. This would also give a more clear presentation of “best” values. Same comment applies to figures 6-9.
Line 625 - I find the appeal to non-sphericity to be unconvincing (although it could even be right). Why not consider that you just don’t have a great simulation of coarse aerosol amounts? You don’t seem to be evaluating it particularly here and so maybe things like dust and sea salt that gets into your domain from remote sources isn't done well (or local dust size distribution).
Line 711 - type - “presented in Sections 4.1 and 4.2”
Figure 14 - State explicitly here that the difference is EMR-REF
Citation: https://doi.org/10.5194/egusphere-2024-2958-RC1 -
AC1: 'Reply to RC1', Kirill Semeniuk, 18 Feb 2025
"It seems that the MOSAIC scheme with the EMR dry deposition produces overall the best set of results, albeit with a computational cost 3x the CAM run. Some improvements to the CAM configuration seem to be possible. A difference in the thermodynamics (no treatment of cations in CAM while they are treated in MOSAIC) are important to the nitrate simulation. There does not seem to be sufficient information given to judge the “best” approach. For example, if the goal is 5 day forecast of PM and it takes 1 hour to produce with CAM and 3 hours to produce with MOSAIC that seems insignificant. If it is 1 day versus 3 days then it matters. I can’t judge this given what was in the paper."
We have revised the Introduction and Conclusion sections to clarify the point raised by the reviewer. The purpose of MOSAIC in GEM-MACH is to improve the scenario simulation capability as opposed to operational forecasts. There are numerous missing process representations in GEM-MACH such as organic thermodynamics and heterogeneous chemistry. It is not possible to include these important processes in the operational model due to computational cost and time penalty. However, GEM-MACH is used for policy scenario simulations in a research mode. Here the time penalty is not critical, and the computational cost of better process representation is justified.
"Line 70 - what is the implication here of chemical and aerosol tracers not being transported by convection? This seems like it would be important to determining the vertical profile. Can you add anything about why this is not a first order issue to be addressed?"
This is indeed an important limitation of GEM-MACH but it is mitigated to a great degree by resolved transport by synoptic and mesoscale systems (e.g. Polvani and Esler, 2007; Lyons et al., 1995). However, ventilation of the atmospheric boundary layer by shallow convection is not properly captured and is a model bias (Polavarapu et al., 2016). We have added more discussion of this aspect in the text and included pertinent references.
"Line 213 - I am confused about what is being discussed here, noting particularly the claim to be using GEOS-5 driving data. Later on line 264 you say you are getting dynamical information from GDPS. Is something in error here?"
The boundary conditions were produced by a MOZART-4 simulation which is a CTM distinct from GEM-MACH. MOZART-4 is driven by GEOS-5 data, GEM-MACH is driven by GDPS and other ECCC produced inputs based on the operational forecast model GEM. We have revised the text to make this distinction clearer.
"Captions to figures 2b and 2c need to be corrected to refer to figure 2a (not 3a)"
This has been corrected.
"The numerical values in the labeling on Figure 5 are impossible to read even blown up on my screen. I suggest that you tabulate these either in the main text of the supplement. This would also give a more clear presentation of “best” values. Same comment applies to figures 6-9."
We have increased the font size in the figures. We have also added section S3 in the supplement which tabulates the similarity metrics in Figure 5.
"Line 625 - I find the appeal to non-sphericity to be unconvincing (although it could even be right). Why not consider that you just don’t have a great simulation of coarse aerosol amounts? You don’t seem to be evaluating it particularly here and so maybe things like dust and sea salt that gets into your domain from remote sources isn't done well (or local dust size distribution)."
We have changed the text to de-emphasize this explanation and point out the underestimation aspect.
"Line 711 - type - “presented in Sections 4.1 and 4.2”
This has been revised.
"Figure 14 - State explicitly here that the difference is EMR-REF"
This correction has been added.
Citation: https://doi.org/10.5194/egusphere-2024-2958-AC1 -
AC2: 'Additional information on RC1', Kirill Semeniuk, 15 Mar 2025
As part of the reply to RC2 we discovered an error in the processing of AERONET VSD. We had Level 1.5 data plotted in Figure 10 instead of Level 2.0. The figure has been revised and the coarse mode peak is now at about 6 um instead of 8 um. The discussion of the disagreement between the model and the observationally inferred VSD has been changed accordingly.
Citation: https://doi.org/10.5194/egusphere-2024-2958-AC2
-
AC1: 'Reply to RC1', Kirill Semeniuk, 18 Feb 2025
-
RC2: 'Comment on egusphere-2024-2958', Anonymous Referee #2, 24 Jan 2025
Implementation of the MOSAIC Aerosol Module (v1.0) in the Canadian Air Quality Model GEM-MACH (v3.1) Semeniuk et al.
Overview:
The paper compares two aerosol modules as part of the GEM-CMACH model. The default CAM module and a new module: the MOSAIC module. Further, the authors compare two configurations of the dry deposition scheme for both aerosol modules. The resulting four model runs for North-America for 2016 are evaluated for nitrates, sulphates and ammonium using a model-observation composite product and speciated PM observations. Further, the AERONET column integrated volume size distribution (VSD) product is used to compare the size distribution of model simulations. Overall, the two aerosol module version do not differ hugely for most aspects. Exceptions are the size distributions over the 12 bins and nitrates over oceans. Overall, MOSAIC offers improvements over CAM, especially for nitrate and ammonium, but comes at a higher computational cost.
General remarks
The paper is very detailed and provides a lot of useful information. However, it lacks a clear message and a prioritisation of the importance of the findings for the reader, who is not directly working with the particular version of the GEM-MACH model.
Overall, the paper is too long, and I suggest omitting certain evaluation aspects for the sake of clarity.
On the other hand, the discussion of the reasons for the differences in the results and their attribution to the differences in the modelling approach can be strengthened. I strongly recommend to add a more concise overview (e.g. table) of the key differences and communalities in the aerosol modelling (order of the scheme, hydration, thermodynamic system solution regimes, role of cations etc )
The comparison of the two dry deposition configurations (i.e. four model runs have be compared rather than two) does not add a lot of useful information because the the science of the dry deposition scheme is not sufficiently discussed. I suggest to present only the runs with the EMR update of the Zhang scheme because it seems to lead to better results.
The comparison of the seasonal maps of surface nitrate, ammonium and sulphate with a model-observation composite product (Donkellaar et al., 2019) should not be called “observations” because it is based on GEOS-Chem, surface observations and uncertain satellite-based PM fields. The usefulness of that composite product for the evaluation should be properly discussed before it used. The observation included in the product are more or less the same observations used by the authors themselves in the next section, which is to some extent a duplication. However, it is good to compare the simulated fields as map to show the differences also in areas without in-situ observations (e.g. over oceans)
The comparison with the size resolved AERONET VSD product offers interesting insight but also the reliability of that data set needs to be further addressed. I also suggest including in the paper a more basic evaluation with AERONET AOD or AE observations. The same holds for a verification with standard PM2.5 observations, which should be reported on.
It would be very welcome if the authors could identify more clearly in the conclusions, which aspects of the MOSAIC module seems to be most important for achieving the reported improvements in model performance.
Citation: https://doi.org/10.5194/egusphere-2024-2958-RC2 -
AC3: 'Reply to RC2', Kirill Semeniuk, 15 Mar 2025
"The paper is very detailed and provides a lot of useful information. However, it lacks a clear message and a prioritisation of the importance of the findings for the reader, who is not directly working with the particular version of the GEM-MACH model. Overall, the paper is too long, and I suggest omitting certain evaluation aspects for the sake of clarity."
We have attempted to address the clarity aspect raised by the reviewer. The additional content required to answer the reviewer's comments as discussed below has been put in the Supplement to reduce the size of the paper. However, we believe that pruning the existing text and analysis for brevity undermines our article. In particular, the comparison between the REF and EMR runs establishes the need for a new dry deposition scheme for aerosols.
"On the other hand, the discussion of the reasons for the differences in the results and their attribution to the differences in the modelling approach can be strengthened. I strongly recommend to add a more concise overview (e.g. table) of the key differences and communalities in the aerosol modelling (order of the scheme, hydration, thermodynamic system solution regimes, role of cations etc )"
We Table 1 in Section 3.1 which describes the differences between CAM and MOSAIC.
"The comparison of the two dry deposition configurations (i.e. four model runs have be compared rather than two) does not add a lot of useful information because the the science of the dry deposition scheme is not sufficiently discussed. I suggest to present only the runs with the EMR update of the Zhang scheme because it seems to lead to better results."
We have added additional discussion and references of the improvements of the Emerson et al. (2020) scheme over the Zhang et al. (2001) scheme used in the reference GEM-MACH. The Zhang scheme has a high bias in the scavenging of the accumulation and Aitken modes. In fact, it is an outlier in having the scavenging minimum occurring above 1 µm. A detailed discussion of the science of dry deposition is beyond the scope of our paper. The Emerson paper makes sufficient justifications, based on extensive new observational data, for the updated parameters it provides. In this it is supported by Pleim et al. (2022) (doi:10.1029/2022MS003050). Our results also support the Emerson et al. (2020) scheme since the Zhang et al. (2001) scheme is negatively impacting the more accurate MOSAIC size distribution while working much better with the pathological CAM size distribution in GEM-MACH. This undermines the performance of MOSAIC compared to CAM. We consider an evaluation of the impact of the Zhang scheme to be relevant information that serves as justification for using the Emerson scheme. The inclusion of the comparison with both REF and EMR runs serves to justify the use of a new dry deposition scheme and we consider this to be a significant result of our study.
"The comparison of the seasonal maps of surface nitrate, ammonium and sulphate with a model-observation composite product (Donkellaar et al., 2019) should not be called “observations” because it is based on GEOS-Chem, surface observations and uncertain satellite-based PM fields. The usefulness of that composite product for the evaluation should be properly discussed before it used. The observation included in the product are more or less the same observations used by the authors themselves in the next section, which is to some extent a duplication. However, it is good to compare the simulated fields as map to show the differences also in areas without in-situ observations (e.g. over oceans)"
We have replaced the term “observation” with “observation product” to refer to the Donkelaar et al. synthesis product in Section 4.1 of the text. But the current text already makes it clear that this is not from direct observations. We believe that justifying the use of this product is unnecessary content. There are no direct observational maps of surface aerosol constituents and objective analysis products have value and are used routinely. The station data is point-wise analysis which we do not consider a duplication.
"The comparison with the size resolved AERONET VSD product offers interesting insight but also the reliability of that data set needs to be further addressed. I also suggest including in the paper a more basic evaluation with AERONET AOD or AE observations. The same holds for a verification with standard PM2.5 observations, which should be reported on."
We have included an evaluation for measurement uncertainty in new Section S3 of the Supplement for surface network stations where possible and refer to these values in the text. We have chosen not to revise Figures 5-9 to include these measurement uncertainties. Networks do not supply measurement error data (precision and bias) in the files they provide, or it is missing. AERONET has uncertainty analysis data which can be downloaded on a station-by-station basis which we have done for the 52 stations used in our paper. For some reason this uncertainty data is not included in the regular aggregated station data files. New Section S4 in the Supplement includes the VSD size distribution uncertainty and provides seasonal average plots with the corresponding error bars in new Figure S15.
We have added content to compare AERONET AOD, single scattering albedo and AE fields with GEM-MACH output in the article (new Section 4.4) and new Section S5 in the Supplement. This analysis is limited by our offline aerosol optical properties code which is based on Bohren and Huffman (1983) and does not consider aerosol mixing state and chemical composition. However, it serves to underscore the difference in the size distribution between CAM and MOSAIC which leads to better agreement overall with observations for MOSAIC. As with the VSD data, this is a limited sampling that is more qualitative in nature.
In the course of carrying out this additional analysis we found a mistake in our original VSD processing where we used Level 1.5 instead of 2.0 products. The error data is only available for Level 2.0 and we have updated Figure 10 accordingly. The coarse mode peak is now at about 6 µm instead of 8 µm which agrees better with the model.
"It would be very welcome if the authors could identify more clearly in the conclusions, which aspects of the MOSAIC module seems to be most important for achieving the reported improvements in model performance."
We have added a short paragraph at the beginning of the Conclusions section to summarize the discussion currently present in the rest of the section.
Citation: https://doi.org/10.5194/egusphere-2024-2958-AC3
-
AC3: 'Reply to RC2', Kirill Semeniuk, 15 Mar 2025
Viewed
HTML | XML | Total | Supplement | BibTeX | EndNote | |
---|---|---|---|---|---|---|
230 | 56 | 30 | 316 | 89 | 12 | 12 |
- HTML: 230
- PDF: 56
- XML: 30
- Total: 316
- Supplement: 89
- BibTeX: 12
- EndNote: 12
Viewed (geographical distribution)
Country | # | Views | % |
---|---|---|---|
United States of America | 1 | 120 | 37 |
Canada | 2 | 41 | 12 |
France | 3 | 29 | 9 |
China | 4 | 15 | 4 |
United Kingdom | 5 | 13 | 4 |
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1
- 120