the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 12 Apr 2023
An optimisation method to improve modelling of wet deposition in atmospheric transport models: applied to FLEXPART v10.4
Abstract. Wet deposition plays a crucial role in the removal of aerosols from the atmosphere. Yet, large uncertainties remain in its implementation in atmospheric transport models, specifically in the parameterisation schemes that are often used. Recently, a new wet deposition scheme was introduced in FLEXPART. The input parameters for its wet deposition scheme can be altered by the user, and may be case-specific. In this paper, a new method is presented to optimise the wet scavenging rates in atmospheric transport models such as FLEXPART. The optimisation scheme is tested in a case study of aerosol-attached 137Cs following the Fukushima Daiichi nuclear power plant accident. From this, improved values for the wet scavenging input parameters in FLEXPART are suggested.
-
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
(4950 KB)
-
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
- Preprint
(4950 KB) - Metadata XML
- BibTeX
- EndNote
- Final revised paper
Journal article(s) based on this preprint
Interactive discussion
Status: closed
-
CEC1: 'Comment on egusphere-2023-357', Juan Antonio Añel, 06 May 2023
Dear authors,
Unfortunately, after checking your manuscript, it has come to our attention that it does not comply with our "Code and Data Policy".
https://www.geoscientific-model-development.net/policies/code_and_data_policy.htmlCurrently, the code of the model that you use, FlexPart, is archived in a repository that does not comply with our policy. The Topical Editor for your manuscript pointed it out before accepting it for Discussions. However, you replied with a repository for a few scripts, not for the model. Despite it, clearly because of an oversight on our side, your manuscript was accepted for Discussions, which is irregular.
Therefore, we have to request you publish the code for the FlexPart version that you use in your manuscript in one of the suitable repositories according to our policy. In this way, you must reply to this comment with the relevant information (link and DOI) as soon as possible, as it should be available for the Discussions stage. Also, please, include another repository with the relevant primary input/output data for your experiments.
Note that you must include in a potentially reviewed version of your manuscript the modified 'Code and Data Availability' section, with the DOIs of the code and data.
I have to highlight that If you fail to fix this problem (please, do it in a prompt manner, as soon as possible), we will have to reject your manuscript for publication in our journal.
Best regards,
Juan A. Añel
Geosci. Model Dev. Exec. EditorÂ
Citation: https://doi.org/10.5194/egusphere-2023-357-CEC1 -
AC1: 'Reply on CEC1', Stijn Van Leuven, 08 May 2023
Dear Juan A. Añel,
Thank you for bringing this to our attention.
The code for Flexpart version 10.4 can be found in the following Zenodo repository: https://doi.org/10.5281/zenodo.3542278. The Flexpart input/output data for our manuscript is archived at https://doi.org/10.5281/zenodo.7906927.
Kind regards,
Stijn Van Leuven
Citation: https://doi.org/10.5194/egusphere-2023-357-AC1
-
AC1: 'Reply on CEC1', Stijn Van Leuven, 08 May 2023
-
RC1: 'Comment on egusphere-2023-357', Sheng Fang, 07 May 2023
1, It seems that by adopting the optimization parameters xi, the pollutants scavenged by the multiple processes at every timestep would be post-processed to make the atmospheric concentration left in the atmosphere closer to the observation. Could you explain how to consider the influence of the changed atmospheric concentration at one timestep to the concentration at the following one?
2, It would be clearer to explain the method mentioned in section 4.4 and the independent verification setup in section 2.Â
3, a Sensitivity study could be done in independent verification, with few measurements adopted to identify the xi and more measurements left for verification.
4, The concentration is a 3D parameter, so how could this paper relate the observation with the delta c in the model? Is the c0 or delta c vertically integrated or are they the value on the ground level? Furthermore, The delta c contributed by the in-cloud scavenging is certainly lower than that contributed by below-cloud scavenging at the ground level, leading to the dominant role of below-cloud scavenging.
5. Why do this study choose the interior algorithm to solve Eq. 15? And how many data are actually used to solve the equation?Citation: https://doi.org/10.5194/egusphere-2023-357-RC1 -
AC2: 'Reply on RC1', Stijn Van Leuven, 22 May 2023
- It seems that by adopting the optimization parameters xi, the pollutants scavenged by the multiple processes at every timestep would be post-processed to make the atmospheric concentration left in the atmosphere closer to the observation. Could you explain how to consider the influence of the changed atmospheric concentration at one timestep to the concentration at the following one?
Flexpart calculates deposition at each timestep. However, our method doesn’t scale the deposition at each timestep independently. Instead, the deposition is scaled once for all locations and timesteps. The influence of the depletion on the concentrations from one timestep to the next is embedded in the factors λi(t) of Eq. (14). As an aside, Flexpart does not model any nonlinear effects on air concentration (such as chemical processes).
Â
- It would be clearer to explain the method mentioned in section 4.4 and the independent verification setup in section 2.
Thank you for this suggestion, we will adjust our manuscript accordingly.
Â
- a Sensitivity study could be done in independent verification, with few measurements adopted to identify the xi and more measurements left for verification.
Following this suggestion, we have performed a sensitivity test by using a random subset of the data. Using only 50% of the data provides consistent optimisation parameters that vary only by around 5% compared to usage of the whole dataset. Using the subset-optimised parameters on the remaining data results in a similar variation of the statistical scores. We will update our manuscript with this extra information.
Â
- The concentration is a 3D parameter, so how could this paper relate the observation with the delta c in the model? Is the c0 or delta c vertically integrated or are they the value on the ground level? Furthermore, The delta c contributed by the in-cloud scavenging is certainly lower than that contributed by below-cloud scavenging at the ground level, leading to the dominant role of below-cloud scavenging.
The concentrations c0 and Δc are not vertically integrated, but represent values as simulated for ground stations at certain locations (see Fig. 1 in the manuscript). These values are calculated during the simulation by the use of ‘receptors’ in Flexpart. This calculation involves the use of a parabolic kernel, which also extends in the vertical direction (to maximally 150m above the surface). More information about the parabolic kernel can be found in Section 8 of ‘Stohl et al. 2005: FLEXPART description’.
The dominant role of below-cloud scavenging, which we observe from the simulations, could indeed be explained by the fact that we consider ground level concentrations. Though, of course, the concentrations as detected by a measurement station may contain particles that have passed through a significant amount of in-cloud scavenging regimes on their way from the release point to the station.
Â
- Why do this study choose the interior algorithm to solve Eq. 15? And how many data are actually used to solve the equation?
We have no specific reason to choose the interior-point algorithm, other than that it is the default algorithm for the optimisation function ‘fmincon’ in MATLAB. Other algorithms such as a trust-region-reflective algorithm were also applied and showed nearly identical results. The equation was solved for all nonzero measurements available, totalling to 248 data-points across 20 locations covering the northern hemisphere (see Figure 1). This will be clarified in our manuscript.
Citation: https://doi.org/10.5194/egusphere-2023-357-AC2
-
AC2: 'Reply on RC1', Stijn Van Leuven, 22 May 2023
-
RC2: 'Comment on egusphere-2023-357', Nina Iren Kristiansen, 25 May 2023
Review of "An optimisation method to improve modelling of we deposition in atmospheric transport models: applied to FLEXPART v10.4"Â by Stijn Van Leuven et al.
General comments
The paper presents a novel method for calculating the individual contributions of several scavenging processes as modelled by an atmospheric transport model. This shows that below cloud scavenging is a much larger contributor than found in previous studies. The scavenging contributions are thereafter optimised based on measurements. Scavenging processes remain challenging to model and parametrise, and this work goes some way of accounting for this model uncertainty. The paper is well written, and all figures are clear and well described. I recommend the paper for publication after considering some of my questions and comments.
Main comments
- Met data
In Stohl et al (2011), FLEXPART simulations using ECMWF were compared with those using GFS met data, and quite large differences were found. “The agreement of model results (both using a priori and a posteriori emissions) with measurement data was better with GFS data than with ECMWF data. The fact that this was also found for 133Xe which is not affected by wet scavenging, shows that GFS-FLEXPART captured the general trans- port better than ECMWF-FLEXPART. Furthermore, the wet scavenging of 137Cs was much stronger with ECMWF data than with the GFS data, causing a strong underestimation of 137Cs concentrations at sites in North America and Europe“.- I believe that the location of a low-pressure system northeast of Japan was the main contributing factor in the difference between ECMWF and GFS for this case. Appreciating that repeating the current study using GFS data is a large task, I would encourage this if possible, or at the minimum include some discussion of the abovementioned findings from the Stohl paper. Your Figure 5 shows 1-month averaged ECMWF precipitation; however, I believe this will smooth out local differences like the position of a low-pressure system over 2-3 days which might play a significant role for scavenging.
- Also, the Stohl study found an underestimation of modelled air concentrations compared to measurements, while the current study shows an overestimating, albeit using different FLEXPART model versions. Some discussion of this in the current paper would be appreciated.
- Aerosol lifetime
In Grythe et al (2017) the aerosol lifetime is calculated, using the Fukushima case study and Cs-137 model and observation data. The modelled aerosol lifetime (~10 days) was found to be too short compared to measurements (~14 days), and comparison between FLEXPART v9 (~6 days) and v10 (~10 days) showed that the new wet deposition in FLEXPART perform better than the old scheme (See Fig 5 in Grythe paper). However, a too short aerosol lifetime indicate too quick removal in the model of which scavenging is a significant contribution. This contradicts the results in this current study which found the air concentrations too high compared to measurements and the scavenging needed to be increased. Some comparison and discussion of the current paper and the aerosol lifetime study would be good to include in the paper. Â
- Use of measurements
Including observation error, i.e., uncertainty in the measurements, is key when comparing model and measurements. Do you account for measurement errors? As FLEXPART does not simulate the background variability in the caesium concentrations this should be accounted for when comparing to the measurements, do you subtract a background value from the measurements?Â
- Deposition vs air concentrations
Your study finds that both air concentrations and deposition are overpredicted by the reference simulation (Fig 2 and Fig 9). I would expect deposition to be underpredicted if air concentrations are overpredicted. Does this reflect measurement uncertainty? Please comment.Â
- Case sensitivity
Are the scavenging coefficients you derive likely to be case specific to the Fukushima case, and therefore not easily generalisable? The Grythe et al 2017 paper indicate they are case dependent. The abstract and conclusion would benefit from a statement whether the derived coefficients are for Fukushima only or suggested for more ‘general’ use.Â
Specific comments
- Some further references in the first two paragraph of the introduction seem necessary.
- L28: (0.1 – 1 um) – please indicate if this is particle radius or diameter.
- Figure 6: The sum of the percentages on the reference simulation figure adds to 101%, while that for the optimised concentrations adds up to 99 %. Perhaps some rounding errors are introduced here? Please correct to 100%.
 Technical comments
L338: FLEPXART to FLEXPART
Citation: https://doi.org/10.5194/egusphere-2023-357-RC2 -
AC3: 'Reply on RC2', Stijn Van Leuven, 30 Jun 2023
We thank you for these comments, all of which have been incorporated into our manuscript. Below are our direct responses to the main comments.
- Met data
In Stohl et al (2011), FLEXPART simulations using ECMWF were compared with those using GFS met data, and quite large differences were found. “The agreement of model results (both using a priori and a posteriori emissions) with measurement data was better with GFS data than with ECMWF data. The fact that this was also found for 133Xe which is not affected by wet scavenging, shows that GFS-FLEXPART captured the general trans- port better than ECMWF-FLEXPART. Furthermore, the wet scavenging of 137Cs was much stronger with ECMWF data than with the GFS data, causing a strong underestimation of 137Cs concentrations at sites in North America and Europe“. - I believe that the location of a low-pressure system northeast of Japan was the main contributing factor in the difference between ECMWF and GFS for this case. Appreciating that repeating the current study using GFS data is a large task, I would encourage this if possible, or at the minimum include some discussion of the abovementioned findings from the Stohl paper. Your Figure 5 shows 1-month averaged ECMWF precipitation; however, I believe this will smooth out local differences like the position of a low-pressure system over 2-3 days which might play a significant role for scavenging.
- Also, the Stohl study found an underestimation of modelled air concentrations compared to measurements, while the current study shows an overestimating, albeit using different FLEXPART model versions. Some discussion of this in the current paper would be appreciated.
- We acknowledge the potential impact of the numerical weather model on our results. The optimisation is also not only sensitive to the precipitation, but also the cloud water fields, as this quantifies below and in-cloud scavenging. This information is read from ECMWF and can also differ between NWP’s. So our optimisation is also influenced by the choice of NWP in this regard. However, we noticed the substantial impact on the air concentrations when comparing Flexpart v9 vs. v10, using the same NWP data (ECMWF). With the main relevant difference between these versions being the wet deposition scheme, this is where we focused our attention.
- Using Flexpart v9, which is also used in Stohl (2011), we find an underestimation of the air concentrations as well. Since scavenging is generally decreased for FPv10 compared to v9, the former will result in greater air concentrations. A discussion on the results of the Stohl et al. 2011 study will be added to the manuscript.
- Aerosol lifetime
In Grythe et al (2017) the aerosol lifetime is calculated, using the Fukushima case study and Cs-137 model and observation data. The modelled aerosol lifetime (~10 days) was found to be too short compared to measurements (~14 days), and comparison between FLEXPART v9 (~6 days) and v10 (~10 days) showed that the new wet deposition in FLEXPART perform better than the old scheme (See Fig 5 in Grythe paper). However, a too short aerosol lifetime indicate too quick removal in the model of which scavenging is a significant contribution. This contradicts the results in this current study which found the air concentrations too high compared to measurements and the scavenging needed to be increased. Some comparison and discussion of the current paper and the aerosol lifetime study would be good to include in the paper.Â
This contradiction can likely be explained by the fact we use only one particle size bin (0.65 µm) while Grythe et al. 2017 uses multiple bins ranging from 0.4 µm to 6.2 µm. Table 3 in Grythe (2017) provides results for a release with all mass in the 0.65 µm bin, which shows a relative concentration bias of 11. This corresponds to an overestimation of the air concentrations, similar to the overestimation that we found. A discussion with a comparison to Grythe et al. 2017 will be added to the manuscript.
The motivation for our choice of particle diameter is that it corresponds to the main mode found in measurements following the Fukushima accident. We have chosen not to fine tune the particle diameter distribution further, since we observed the dry deposition to be more sensitive to the particle size compared to (below-cloud) scavenging. Our focus was on the wet scavenging due to the differences seen between Flexpart v9 and v10’s wet deposition scheme.
- Use of measurements
Including observation error, i.e., uncertainty in the measurements, is key when comparing model and measurements. Do you account for measurement errors? As FLEXPART does not simulate the background variability in the caesium concentrations this should be accounted for when comparing to the measurements, do you subtract a background value from the measurements?Â
We didn’t take into account measurement error explicitly, as we found this to be negligible compared to the measured values. The caesium background is not taken into account for similar reasons. This background is much smaller than the measured values, usually being on the order of 0.1-1 µBq/m³. If present, it should manifest as a bias at lower concentration values, which is not seen in the data (see Figure 2). The xenon background, on the other hand, is visible in the data (Figure 3) and is removed with a cut-off for the analysis. This information will be added to our manuscript.
- Deposition vs air concentrations
Your study finds that both air concentrations and deposition are overpredicted by the reference simulation (Fig 2 and Fig 9). I would expect deposition to be underpredicted if air concentrations are overpredicted. Does this reflect measurement uncertainty? Please comment.Â
There are two competing effects, especially manifesting at longer ranges. One is that an increase in deposition reduces air concentration, as mentioned. Another is that a higher air concentration increases scavenging and thus deposition. These two effects almost cancel out in the reference simulation, leading to a small overprediction of the deposition. The deposition measurements are taken from the U.S. which is some thousands of km’s from Fukushima.
- Case sensitivity
Are the scavenging coefficients you derive likely to be case specific to the Fukushima case, and therefore not easily generalisable? The Grythe et al 2017 paper indicate they are case dependent. The abstract and conclusion would benefit from a statement whether the derived coefficients are for Fukushima only or suggested for more ‘general’ use.Â
Thank you for this suggestion. We hypothesise that the scavenging coefficients should be valid more generally, under similar conditions. The NWP, for example, influences the quantification of below- and in-cloud scavenging. Thusly, the optimised coefficients are not necessarily applicable to other NWP’s. Similarly, according to Grythe et al. 2017, the scavenging coefficients may depend on a case-by-case basis. We will update our manuscript accordingly.
Citation: https://doi.org/10.5194/egusphere-2023-357-AC3 - Met data
- Met data
Interactive discussion
Status: closed
-
CEC1: 'Comment on egusphere-2023-357', Juan Antonio Añel, 06 May 2023
Dear authors,
Unfortunately, after checking your manuscript, it has come to our attention that it does not comply with our "Code and Data Policy".
https://www.geoscientific-model-development.net/policies/code_and_data_policy.htmlCurrently, the code of the model that you use, FlexPart, is archived in a repository that does not comply with our policy. The Topical Editor for your manuscript pointed it out before accepting it for Discussions. However, you replied with a repository for a few scripts, not for the model. Despite it, clearly because of an oversight on our side, your manuscript was accepted for Discussions, which is irregular.
Therefore, we have to request you publish the code for the FlexPart version that you use in your manuscript in one of the suitable repositories according to our policy. In this way, you must reply to this comment with the relevant information (link and DOI) as soon as possible, as it should be available for the Discussions stage. Also, please, include another repository with the relevant primary input/output data for your experiments.
Note that you must include in a potentially reviewed version of your manuscript the modified 'Code and Data Availability' section, with the DOIs of the code and data.
I have to highlight that If you fail to fix this problem (please, do it in a prompt manner, as soon as possible), we will have to reject your manuscript for publication in our journal.
Best regards,
Juan A. Añel
Geosci. Model Dev. Exec. EditorÂ
Citation: https://doi.org/10.5194/egusphere-2023-357-CEC1 -
AC1: 'Reply on CEC1', Stijn Van Leuven, 08 May 2023
Dear Juan A. Añel,
Thank you for bringing this to our attention.
The code for Flexpart version 10.4 can be found in the following Zenodo repository: https://doi.org/10.5281/zenodo.3542278. The Flexpart input/output data for our manuscript is archived at https://doi.org/10.5281/zenodo.7906927.
Kind regards,
Stijn Van Leuven
Citation: https://doi.org/10.5194/egusphere-2023-357-AC1
-
AC1: 'Reply on CEC1', Stijn Van Leuven, 08 May 2023
-
RC1: 'Comment on egusphere-2023-357', Sheng Fang, 07 May 2023
1, It seems that by adopting the optimization parameters xi, the pollutants scavenged by the multiple processes at every timestep would be post-processed to make the atmospheric concentration left in the atmosphere closer to the observation. Could you explain how to consider the influence of the changed atmospheric concentration at one timestep to the concentration at the following one?
2, It would be clearer to explain the method mentioned in section 4.4 and the independent verification setup in section 2.Â
3, a Sensitivity study could be done in independent verification, with few measurements adopted to identify the xi and more measurements left for verification.
4, The concentration is a 3D parameter, so how could this paper relate the observation with the delta c in the model? Is the c0 or delta c vertically integrated or are they the value on the ground level? Furthermore, The delta c contributed by the in-cloud scavenging is certainly lower than that contributed by below-cloud scavenging at the ground level, leading to the dominant role of below-cloud scavenging.
5. Why do this study choose the interior algorithm to solve Eq. 15? And how many data are actually used to solve the equation?Citation: https://doi.org/10.5194/egusphere-2023-357-RC1 -
AC2: 'Reply on RC1', Stijn Van Leuven, 22 May 2023
- It seems that by adopting the optimization parameters xi, the pollutants scavenged by the multiple processes at every timestep would be post-processed to make the atmospheric concentration left in the atmosphere closer to the observation. Could you explain how to consider the influence of the changed atmospheric concentration at one timestep to the concentration at the following one?
Flexpart calculates deposition at each timestep. However, our method doesn’t scale the deposition at each timestep independently. Instead, the deposition is scaled once for all locations and timesteps. The influence of the depletion on the concentrations from one timestep to the next is embedded in the factors λi(t) of Eq. (14). As an aside, Flexpart does not model any nonlinear effects on air concentration (such as chemical processes).
Â
- It would be clearer to explain the method mentioned in section 4.4 and the independent verification setup in section 2.
Thank you for this suggestion, we will adjust our manuscript accordingly.
Â
- a Sensitivity study could be done in independent verification, with few measurements adopted to identify the xi and more measurements left for verification.
Following this suggestion, we have performed a sensitivity test by using a random subset of the data. Using only 50% of the data provides consistent optimisation parameters that vary only by around 5% compared to usage of the whole dataset. Using the subset-optimised parameters on the remaining data results in a similar variation of the statistical scores. We will update our manuscript with this extra information.
Â
- The concentration is a 3D parameter, so how could this paper relate the observation with the delta c in the model? Is the c0 or delta c vertically integrated or are they the value on the ground level? Furthermore, The delta c contributed by the in-cloud scavenging is certainly lower than that contributed by below-cloud scavenging at the ground level, leading to the dominant role of below-cloud scavenging.
The concentrations c0 and Δc are not vertically integrated, but represent values as simulated for ground stations at certain locations (see Fig. 1 in the manuscript). These values are calculated during the simulation by the use of ‘receptors’ in Flexpart. This calculation involves the use of a parabolic kernel, which also extends in the vertical direction (to maximally 150m above the surface). More information about the parabolic kernel can be found in Section 8 of ‘Stohl et al. 2005: FLEXPART description’.
The dominant role of below-cloud scavenging, which we observe from the simulations, could indeed be explained by the fact that we consider ground level concentrations. Though, of course, the concentrations as detected by a measurement station may contain particles that have passed through a significant amount of in-cloud scavenging regimes on their way from the release point to the station.
Â
- Why do this study choose the interior algorithm to solve Eq. 15? And how many data are actually used to solve the equation?
We have no specific reason to choose the interior-point algorithm, other than that it is the default algorithm for the optimisation function ‘fmincon’ in MATLAB. Other algorithms such as a trust-region-reflective algorithm were also applied and showed nearly identical results. The equation was solved for all nonzero measurements available, totalling to 248 data-points across 20 locations covering the northern hemisphere (see Figure 1). This will be clarified in our manuscript.
Citation: https://doi.org/10.5194/egusphere-2023-357-AC2
-
AC2: 'Reply on RC1', Stijn Van Leuven, 22 May 2023
-
RC2: 'Comment on egusphere-2023-357', Nina Iren Kristiansen, 25 May 2023
Review of "An optimisation method to improve modelling of we deposition in atmospheric transport models: applied to FLEXPART v10.4"Â by Stijn Van Leuven et al.
General comments
The paper presents a novel method for calculating the individual contributions of several scavenging processes as modelled by an atmospheric transport model. This shows that below cloud scavenging is a much larger contributor than found in previous studies. The scavenging contributions are thereafter optimised based on measurements. Scavenging processes remain challenging to model and parametrise, and this work goes some way of accounting for this model uncertainty. The paper is well written, and all figures are clear and well described. I recommend the paper for publication after considering some of my questions and comments.
Main comments
- Met data
In Stohl et al (2011), FLEXPART simulations using ECMWF were compared with those using GFS met data, and quite large differences were found. “The agreement of model results (both using a priori and a posteriori emissions) with measurement data was better with GFS data than with ECMWF data. The fact that this was also found for 133Xe which is not affected by wet scavenging, shows that GFS-FLEXPART captured the general trans- port better than ECMWF-FLEXPART. Furthermore, the wet scavenging of 137Cs was much stronger with ECMWF data than with the GFS data, causing a strong underestimation of 137Cs concentrations at sites in North America and Europe“.- I believe that the location of a low-pressure system northeast of Japan was the main contributing factor in the difference between ECMWF and GFS for this case. Appreciating that repeating the current study using GFS data is a large task, I would encourage this if possible, or at the minimum include some discussion of the abovementioned findings from the Stohl paper. Your Figure 5 shows 1-month averaged ECMWF precipitation; however, I believe this will smooth out local differences like the position of a low-pressure system over 2-3 days which might play a significant role for scavenging.
- Also, the Stohl study found an underestimation of modelled air concentrations compared to measurements, while the current study shows an overestimating, albeit using different FLEXPART model versions. Some discussion of this in the current paper would be appreciated.
- Aerosol lifetime
In Grythe et al (2017) the aerosol lifetime is calculated, using the Fukushima case study and Cs-137 model and observation data. The modelled aerosol lifetime (~10 days) was found to be too short compared to measurements (~14 days), and comparison between FLEXPART v9 (~6 days) and v10 (~10 days) showed that the new wet deposition in FLEXPART perform better than the old scheme (See Fig 5 in Grythe paper). However, a too short aerosol lifetime indicate too quick removal in the model of which scavenging is a significant contribution. This contradicts the results in this current study which found the air concentrations too high compared to measurements and the scavenging needed to be increased. Some comparison and discussion of the current paper and the aerosol lifetime study would be good to include in the paper. Â
- Use of measurements
Including observation error, i.e., uncertainty in the measurements, is key when comparing model and measurements. Do you account for measurement errors? As FLEXPART does not simulate the background variability in the caesium concentrations this should be accounted for when comparing to the measurements, do you subtract a background value from the measurements?Â
- Deposition vs air concentrations
Your study finds that both air concentrations and deposition are overpredicted by the reference simulation (Fig 2 and Fig 9). I would expect deposition to be underpredicted if air concentrations are overpredicted. Does this reflect measurement uncertainty? Please comment.Â
- Case sensitivity
Are the scavenging coefficients you derive likely to be case specific to the Fukushima case, and therefore not easily generalisable? The Grythe et al 2017 paper indicate they are case dependent. The abstract and conclusion would benefit from a statement whether the derived coefficients are for Fukushima only or suggested for more ‘general’ use.Â
Specific comments
- Some further references in the first two paragraph of the introduction seem necessary.
- L28: (0.1 – 1 um) – please indicate if this is particle radius or diameter.
- Figure 6: The sum of the percentages on the reference simulation figure adds to 101%, while that for the optimised concentrations adds up to 99 %. Perhaps some rounding errors are introduced here? Please correct to 100%.
 Technical comments
L338: FLEPXART to FLEXPART
Citation: https://doi.org/10.5194/egusphere-2023-357-RC2 -
AC3: 'Reply on RC2', Stijn Van Leuven, 30 Jun 2023
We thank you for these comments, all of which have been incorporated into our manuscript. Below are our direct responses to the main comments.
- Met data
In Stohl et al (2011), FLEXPART simulations using ECMWF were compared with those using GFS met data, and quite large differences were found. “The agreement of model results (both using a priori and a posteriori emissions) with measurement data was better with GFS data than with ECMWF data. The fact that this was also found for 133Xe which is not affected by wet scavenging, shows that GFS-FLEXPART captured the general trans- port better than ECMWF-FLEXPART. Furthermore, the wet scavenging of 137Cs was much stronger with ECMWF data than with the GFS data, causing a strong underestimation of 137Cs concentrations at sites in North America and Europe“. - I believe that the location of a low-pressure system northeast of Japan was the main contributing factor in the difference between ECMWF and GFS for this case. Appreciating that repeating the current study using GFS data is a large task, I would encourage this if possible, or at the minimum include some discussion of the abovementioned findings from the Stohl paper. Your Figure 5 shows 1-month averaged ECMWF precipitation; however, I believe this will smooth out local differences like the position of a low-pressure system over 2-3 days which might play a significant role for scavenging.
- Also, the Stohl study found an underestimation of modelled air concentrations compared to measurements, while the current study shows an overestimating, albeit using different FLEXPART model versions. Some discussion of this in the current paper would be appreciated.
- We acknowledge the potential impact of the numerical weather model on our results. The optimisation is also not only sensitive to the precipitation, but also the cloud water fields, as this quantifies below and in-cloud scavenging. This information is read from ECMWF and can also differ between NWP’s. So our optimisation is also influenced by the choice of NWP in this regard. However, we noticed the substantial impact on the air concentrations when comparing Flexpart v9 vs. v10, using the same NWP data (ECMWF). With the main relevant difference between these versions being the wet deposition scheme, this is where we focused our attention.
- Using Flexpart v9, which is also used in Stohl (2011), we find an underestimation of the air concentrations as well. Since scavenging is generally decreased for FPv10 compared to v9, the former will result in greater air concentrations. A discussion on the results of the Stohl et al. 2011 study will be added to the manuscript.
- Aerosol lifetime
In Grythe et al (2017) the aerosol lifetime is calculated, using the Fukushima case study and Cs-137 model and observation data. The modelled aerosol lifetime (~10 days) was found to be too short compared to measurements (~14 days), and comparison between FLEXPART v9 (~6 days) and v10 (~10 days) showed that the new wet deposition in FLEXPART perform better than the old scheme (See Fig 5 in Grythe paper). However, a too short aerosol lifetime indicate too quick removal in the model of which scavenging is a significant contribution. This contradicts the results in this current study which found the air concentrations too high compared to measurements and the scavenging needed to be increased. Some comparison and discussion of the current paper and the aerosol lifetime study would be good to include in the paper.Â
This contradiction can likely be explained by the fact we use only one particle size bin (0.65 µm) while Grythe et al. 2017 uses multiple bins ranging from 0.4 µm to 6.2 µm. Table 3 in Grythe (2017) provides results for a release with all mass in the 0.65 µm bin, which shows a relative concentration bias of 11. This corresponds to an overestimation of the air concentrations, similar to the overestimation that we found. A discussion with a comparison to Grythe et al. 2017 will be added to the manuscript.
The motivation for our choice of particle diameter is that it corresponds to the main mode found in measurements following the Fukushima accident. We have chosen not to fine tune the particle diameter distribution further, since we observed the dry deposition to be more sensitive to the particle size compared to (below-cloud) scavenging. Our focus was on the wet scavenging due to the differences seen between Flexpart v9 and v10’s wet deposition scheme.
- Use of measurements
Including observation error, i.e., uncertainty in the measurements, is key when comparing model and measurements. Do you account for measurement errors? As FLEXPART does not simulate the background variability in the caesium concentrations this should be accounted for when comparing to the measurements, do you subtract a background value from the measurements?Â
We didn’t take into account measurement error explicitly, as we found this to be negligible compared to the measured values. The caesium background is not taken into account for similar reasons. This background is much smaller than the measured values, usually being on the order of 0.1-1 µBq/m³. If present, it should manifest as a bias at lower concentration values, which is not seen in the data (see Figure 2). The xenon background, on the other hand, is visible in the data (Figure 3) and is removed with a cut-off for the analysis. This information will be added to our manuscript.
- Deposition vs air concentrations
Your study finds that both air concentrations and deposition are overpredicted by the reference simulation (Fig 2 and Fig 9). I would expect deposition to be underpredicted if air concentrations are overpredicted. Does this reflect measurement uncertainty? Please comment.Â
There are two competing effects, especially manifesting at longer ranges. One is that an increase in deposition reduces air concentration, as mentioned. Another is that a higher air concentration increases scavenging and thus deposition. These two effects almost cancel out in the reference simulation, leading to a small overprediction of the deposition. The deposition measurements are taken from the U.S. which is some thousands of km’s from Fukushima.
- Case sensitivity
Are the scavenging coefficients you derive likely to be case specific to the Fukushima case, and therefore not easily generalisable? The Grythe et al 2017 paper indicate they are case dependent. The abstract and conclusion would benefit from a statement whether the derived coefficients are for Fukushima only or suggested for more ‘general’ use.Â
Thank you for this suggestion. We hypothesise that the scavenging coefficients should be valid more generally, under similar conditions. The NWP, for example, influences the quantification of below- and in-cloud scavenging. Thusly, the optimised coefficients are not necessarily applicable to other NWP’s. Similarly, according to Grythe et al. 2017, the scavenging coefficients may depend on a case-by-case basis. We will update our manuscript accordingly.
Citation: https://doi.org/10.5194/egusphere-2023-357-AC3 - Met data
- Met data
Peer review completion
Journal article(s) based on this preprint
Model code and software
An optimisation method to improve modelling of wet deposition in atmospheric transport models: applied to FLEXPART v10.4 Stijn Van Leuven https://doi.org/10.5281/zenodo.7789039
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 285 | 129 | 18 | 432 | 8 | 4 |
- HTML: 285
- PDF: 129
- XML: 18
- Total: 432
- BibTeX: 8
- EndNote: 4
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
Cited
Stijn Van Leuven
Pieter De Meutter
Johan Camps
Piet Termonia
Andy Delcloo
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
- Preprint
(4950 KB) - Metadata XML