the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Exploring ship track spreading rates with a physics-informed Langevin particle parameterization
Abstract. The rate at which aerosols spread from a point source injection, such as from a ship or other stationary pollution source, is critical for accurately representing subgrid plume spreading in a climate model. Such climate model results will guide future decisions regarding the feasibility and application of large-scale intentional marine cloud brightening (MCB). Prior modeling studies have shown that the rate at which ship plumes spread may be strongly dependent on meteorological conditions, such as precipitating versus non-precipitating boundary layers and shear. In this study, we apply a Lagrangian particle model (PM-ABL v1.0), governed by a Langevin stochastic differential equation, to create a simplified framework for predicting the rate of spreading from a ship-injected aerosol plume in sheared, precipitating, and non-precipitating boundary layers. The velocity and position of each stochastic particle is predicted with the acceleration of each particle being driven by the turbulent kinetic energy, dissipation rate, momentum variance, and mean wind. These inputs to the stochastic particle-velocity equation are derived from high-fidelity large-eddy simulations (LES) equipped with a prognostic aerosol-cloud microphysics scheme (UWSAM) to simulate an aerosol injection from a ship into a cloud-topped marine boundary layer. The resulting spreading rate from the reduced-order stochastic model is then compared to the spreading rate in the LES. The stochastic particle-velocity representation is shown to reasonably reproduce spreading rates in sheared, precipitating, and non-precipitating cases using domain-averaged turbulent statistics from the LES.
-
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
(16392 KB)
-
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
- Preprint
(16392 KB) - Metadata XML
- BibTeX
- EndNote
- Final revised paper
Journal article(s) based on this preprint
Interactive discussion
Status: closed
-
RC1: 'Comment on egusphere-2024-235', Anonymous Referee #1, 05 May 2024
- AC1: 'Reply on RC1', Lucas McMichael, 23 Jun 2024
-
RC2: 'Comment on egusphere-2024-235', Anonymous Referee #2, 27 May 2024
This article presents a novel approach to modeling the spreading of aerosol plumes in a marine boundary layer capped with stratocumulus clouds. The model is inspired from the stochastic processes based approaches used in turbulent flows for representing Lagrangian trajectories. The results presented in this work will be of interest to the marine cloud brightening and ship-track community (broadly speaking - ACI community) and should be considered for publication after the comments listed below are carefully addressed.
Comments:
1. Can you explain the physical meaning behind the terms in the GLM/SLM model?
2. Eq 16, why would the boundary layer flow allow isotropic coeff for the drift term? Is it because you are considering only the horizontal dimensions?3. In eq 18, are TL and Co constants?
4. What are the overbar terms in Eq 23 and 14? Shouldn't those be zero?
5. Â Eq 30, explain how TL and \sigma are obtained?6. Lines 280-290: Are the surface fluxes identical for all the grid and domain sizes? I am wondering about the differences in LWP/precip. and if they are related to the differences in the surface fluxes? So, is Eq. 32 required?
7. Line 315: the unit of shear is wrong.8. Section 3.1, the forcing used for obtaining the different shear rates are not clear? Note that other readers should be able to reproduce these results. More details would be useful.
9. Â Lines 380-395: I think the figures numbers are wrong. Please check.10. Lines 402-405: The reasoning is very vague. You need to show that there is some difference in the energy cascade. Can you?Â
11. Same for line 406.12. Line 408: what do you mean by efficiency of dissipation?
13. Line 439: why is sgs diffusivity the relevant parameter? Can't you find an effective diffusivity from the resolved and SGS stresses? That would be the appropriate parameter.
14. Line 448: which part has the superlinear growth? Can you show it in the figure? I dont find any superlinear growth. Can you clarify? And what about numerical diffusion? How do you account for that?
15. How do your growth rates compare against other studies?
16. Fig 10, what about the diffusive model with diffusivity changing with time? And can you plot the diffusive model for each case and not just the control case (maybe in Fig 12). This is required to check if the particle model is actually better.
17. Why is the plume growth rate lower in the night time? Was it seen in other studies as well?18. Sec 4.2, how is TL for the particle model determined? Is it time varying? I see that Co is from a least sq. fit. Why not make it time varying as well?
19. Lines 474-479: Please show these.20. The authors have articulated in the article that the particle based estimates are better than the Gaussian-diffusion based estimates. But the Ornstein-Uhlenbeck process used here for obtaining the velocity field is a Gaussian diffusive process (i.e., random walk based on a Gaussian distribution). Then why would the two give different results? My impression was that it's a question of having the right effective diffusivity. A discussion related to this would be useful.
21. References: I think the reference to Prabhakaran et al 2024 (ACP) is more relevant than 2023 (JAS).
22. Lines 45-50: Do the other recent LES studies agree with the conclusions of Chun et al 2023?
23. Lines 78-90: parts of this paragraph doesn't seem like introductory text. This could go into abstract or in sec. 4.2. Sec 4.2 would be preferred.
Citation: https://doi.org/10.5194/egusphere-2024-235-RC2 - AC2: 'Reply on RC2', Lucas McMichael, 23 Jun 2024
Interactive discussion
Status: closed
-
RC1: 'Comment on egusphere-2024-235', Anonymous Referee #1, 05 May 2024
- AC1: 'Reply on RC1', Lucas McMichael, 23 Jun 2024
-
RC2: 'Comment on egusphere-2024-235', Anonymous Referee #2, 27 May 2024
This article presents a novel approach to modeling the spreading of aerosol plumes in a marine boundary layer capped with stratocumulus clouds. The model is inspired from the stochastic processes based approaches used in turbulent flows for representing Lagrangian trajectories. The results presented in this work will be of interest to the marine cloud brightening and ship-track community (broadly speaking - ACI community) and should be considered for publication after the comments listed below are carefully addressed.
Comments:
1. Can you explain the physical meaning behind the terms in the GLM/SLM model?
2. Eq 16, why would the boundary layer flow allow isotropic coeff for the drift term? Is it because you are considering only the horizontal dimensions?3. In eq 18, are TL and Co constants?
4. What are the overbar terms in Eq 23 and 14? Shouldn't those be zero?
5. Â Eq 30, explain how TL and \sigma are obtained?6. Lines 280-290: Are the surface fluxes identical for all the grid and domain sizes? I am wondering about the differences in LWP/precip. and if they are related to the differences in the surface fluxes? So, is Eq. 32 required?
7. Line 315: the unit of shear is wrong.8. Section 3.1, the forcing used for obtaining the different shear rates are not clear? Note that other readers should be able to reproduce these results. More details would be useful.
9. Â Lines 380-395: I think the figures numbers are wrong. Please check.10. Lines 402-405: The reasoning is very vague. You need to show that there is some difference in the energy cascade. Can you?Â
11. Same for line 406.12. Line 408: what do you mean by efficiency of dissipation?
13. Line 439: why is sgs diffusivity the relevant parameter? Can't you find an effective diffusivity from the resolved and SGS stresses? That would be the appropriate parameter.
14. Line 448: which part has the superlinear growth? Can you show it in the figure? I dont find any superlinear growth. Can you clarify? And what about numerical diffusion? How do you account for that?
15. How do your growth rates compare against other studies?
16. Fig 10, what about the diffusive model with diffusivity changing with time? And can you plot the diffusive model for each case and not just the control case (maybe in Fig 12). This is required to check if the particle model is actually better.
17. Why is the plume growth rate lower in the night time? Was it seen in other studies as well?18. Sec 4.2, how is TL for the particle model determined? Is it time varying? I see that Co is from a least sq. fit. Why not make it time varying as well?
19. Lines 474-479: Please show these.20. The authors have articulated in the article that the particle based estimates are better than the Gaussian-diffusion based estimates. But the Ornstein-Uhlenbeck process used here for obtaining the velocity field is a Gaussian diffusive process (i.e., random walk based on a Gaussian distribution). Then why would the two give different results? My impression was that it's a question of having the right effective diffusivity. A discussion related to this would be useful.
21. References: I think the reference to Prabhakaran et al 2024 (ACP) is more relevant than 2023 (JAS).
22. Lines 45-50: Do the other recent LES studies agree with the conclusions of Chun et al 2023?
23. Lines 78-90: parts of this paragraph doesn't seem like introductory text. This could go into abstract or in sec. 4.2. Sec 4.2 would be preferred.
Citation: https://doi.org/10.5194/egusphere-2024-235-RC2 - AC2: 'Reply on RC2', Lucas McMichael, 23 Jun 2024
Peer review completion
Journal article(s) based on this preprint
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
446 | 142 | 33 | 621 | 17 | 22 |
- HTML: 446
- PDF: 142
- XML: 33
- Total: 621
- BibTeX: 17
- EndNote: 22
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1
Michael J. Schmidt
Robert Wood
Peter N. Blossey
Lekha Patel
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
- Preprint
(16392 KB) - Metadata XML