the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Evaluation of plume rise parameterizations in GEM-MACHv2 with analysis of image data using a deep convolutional neural network
Abstract. The study of plume rise from smokestacks and other pollutant point sources is extremely important for the estimation and modelling of the dispersion of pollutants on regional scales via atmospheric modelling platforms. However, the algorithms which have been used to represent plume rise were based on observations conducted nearly 50 years ago (the semi-empirical dimensional modelling framework of Briggs, 1984), and more recent measurement techniques are available which can be used to generate new data, against which pollutant plume rise theories may be evaluated. A key result of the theoretical formulations based on these past observations is the height reached by the plumes (the process by which they reach that height is known as plume rise). In this work, a previously developed deep convolutional neural network (Deep Plume Rise Network, DPRNet) for determining plume rise from visible RGB images was applied to images taken of a facility in the Athabasca oil sands and compared to the theoretical estimates of Briggs parameterizations as formulated in GEM-MACHv2. On average, the Briggs parameterizations tend to predict plume rise in stable and neutral conditions within 30 %, but consistently overpredict plume rise during unstable conditions by more than 100 %. Further, while Briggs parameterizations predicted diurnal variations in plume rise, no such variation was observed by the image analysis. The parameterizations could be improved reducing dimensionless constants by factors of 2 and 6 in neutral and unstable conditions, respectively. The plume height data have been shown to provide a significant resource for plume rise theory evaluation and development.
- Preprint
(1854 KB) - Metadata XML
- BibTeX
- EndNote
Status: final response (author comments only)
- RC1: 'Comment on egusphere-2025-4582', Anonymous Referee #1, 21 Nov 2025
-
CEC1: 'Comment on egusphere-2025-4582 - No compliance with the policy of the journal', Juan Antonio Añel, 07 Dec 2025
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.htmlFirst, one of the links that you provide to store part of the assets necessary to replicate your work does not correspond to an acceptable repository according to our policy. This is the case of Borealis.
Second, you have not provided the code for the deep convolutional neural network that you use in your manuscript, which is mandatory according to our policy.
Therefore, the current situation with your manuscript is irregular. Please, publish your code and data in one of the appropriate repositories according to our policy and reply to this comment with the relevant information (link and a permanent identifier for it (e.g. DOI)) as soon as possible, as we can not accept manuscripts in Discussions that do not comply with our policy, and given the no compliance of your manuscript, it should not be accepted for Discussions or review in our journal.
I must note that if you do not fix this problem, we cannot continue with the peer-review process or accept your manuscript for publication in our journal.
Juan A. Añel
Geosci. Model Dev. Executive EditorCitation: https://doi.org/10.5194/egusphere-2025-4582-CEC1 -
AC1: 'Reply on CEC1', Mark Gordon, 09 Jan 2026
We thank the editor for noting the omission of the required code and the use of a non-compliant repository. We have addressed these points as follows…
The Borealis dataset has been duplicated with an exact copy on Zenodo. To avoid confusion and to make the duplication as transparent as possible, the same DOI is used for the new Zenodo dataset: https://doi.org/10.5683/SP3/WZVZBV. Currently this DOI links to the Borealis site. In the Borealis metadata, there is a link to the Zenodo dataset: https://zenodo.org/records/18087309 and there is an explanation of the duplication on both the Borealis and Zenodo sites. [As an aside: The Borealis repository is working to guarantee long-term preservation for future publications.]
The Neural Network (DCNN) code has been added to the FRDR dataset in Version 2 which has the following DOI: https://doi.org/10.20383/103.01559.
We hope that this satisfies the requirements. Please let us know if further modifications are needed.
Citation: https://doi.org/10.5194/egusphere-2025-4582-AC1 -
CEC2: 'Reply on AC1', Juan Antonio Añel, 11 Jan 2026
Dear authors,
Many thanks for your reply. We can consider now the current version of your manuscript in compliance with the Code and Data policy of the journal. Please, if eventually the Topical Editor decides request further review of your manuscript or accept it for publication, modify the "Code and Data Accessibility" section in your manuscript and replace the link and information for Borealis with the information for the Zenodo repository.
Juan A. Añel
Geosci. Model Dev. Executive Editor
Citation: https://doi.org/10.5194/egusphere-2025-4582-CEC2
-
CEC2: 'Reply on AC1', Juan Antonio Añel, 11 Jan 2026
-
AC1: 'Reply on CEC1', Mark Gordon, 09 Jan 2026
-
RC2: 'Comment on egusphere-2025-4582', Anonymous Referee #2, 09 Feb 2026
This study uses a deep convolutional neural network framework developed by Koushafar et al. (2023) to analyse images of smoke plumes from an industrial site, comparing the plume rise estimates against theoretical predictions using Briggs parameterisation. The authors concluded that Briggs consistently overpredicts plume rise during unstable conditions with high wind speeds, and that neural network-based plume rise estimates show no diurnal variability compared to those from Briggs parameterisation.
The application of deep convolutional neural networks to plume image analysis represents a methodological advancement for estimating plume height over extended time series compared to previous studies relying on in-situ or flight measurements. Based on their results, the authors proposed revised values for the Briggs parameterization to improve plume height estimation accuracy. Overall, this is an interesting study that advances the methodology for estimating plume rise through the application of deep convolutional neural network. However, the manuscript would benefit from improvements in language clarity and revision of several figures to enhance readability and accessibility.
Major Comments:
- The paper would benefit from more precise language to improve readability and flow. The text contains repetitive wording in several sections. Additionally, the citation style is inconsistent throughout the manuscript. I have identified some of the specific citation issues in the minor comments below, but the authors should conduct a thorough review of all in-text citations.
- In the introduction, the authors reference several previous studies that estimated plume rise and compared results to Briggs parameterization (e.g., Gordon et al., 2015; Akingunola et al., 2018; Gordon et al., 2018; Fathi et al., 2025). However, the distinctions between these approaches are not clearly articulated. What are the potential reasons for the over- or under-estimation of plume heights compared to observations? Furthermore, the paper lacks substantive comparison with these previous studies. The limited comparisons provided (e.g., Lines 542-545) would be strengthened by incorporating more quantitative analysis.
- Section 2.2 would benefit from a brief summary of the DCNN-based framework to help readers better understand the methodology.
- Several figures are not colorblind-friendly, including Figures 1, 2, 9, and 10. The authors should revise the color schemes to avoid using red and green simultaneously in the same plot. The resolution in Figure 2 is very low and the annotations are not easy to read because of the colour scheme chosen.
- In Figures 9 and 10, why are the data for Manual/Visible/Unstable and Manual/Exp Fit/Unstable not shown?
- The authors can highlight the newly fitted dimensionless constant values in the conclusion, as these values could be valuable for future studies applying this parameterisation.
Minor comments:
- Line 33, be more specific about the dimensionless constants
- Line 34 and 37, repetitive about the observations made
- Line 40, missing citation for Moore
- Line 48, what kind of surface concentrations measurements?
- Line 52-55 used very similar wordings to Koushafar et al. (2023). The authors should consider paraphrasing the sentences.
- Line 56, missing comma in citation
- Line 60, check citation style
- Line 65, the authors commented that the studies referenced are relatively short in duration. Is it on the scale of hours to days? The authors can be more specific about “long term data collection” here, and mention seasonal variations in atmospheric conditions that may affect plume rise theory.
- Line 68, since the main stack is already defined here, suggest keeping the terminology throughout the manuscript
- Line 73, check citation style
- Line 73-74, subscripts for CO2, and also elsewhere throughout the paper (e.g., SO2)
- Line 77-78, “most of the emission” is a bit vague. What’s the percentage of emission from the main stack?
- Line 78, it is unclear what is ordinary operating condition
- Line 80, this paragraph is disconnected from the previous paragraph.
- Line 83, what are traditional methods?
- Line 83, check citation style
- Line 90, I find the wordings “determine plume rise” not very clear in this sentence.
- Line 91, the authors mentioned that Koushafar et al. (2023) used the same initial dataset to evaluate DPRNet. Is the image data used in this paper part of the training dataset or validating dataset?
- Line 99, to compare with the parametrisation -> to compare with the Briggs’ parametrisation
- Line 103-108, a lot of the description of the research site is not relevant to the analysis. The authors should consider revising this paragraph.
- Line 108, missing fullstops.
- Line 132, main stack already defined earlier
- Line 133, not necessary to mention 600 feet
- Line 151, move the sentence out of the bracket.
- Line 168, missing comma
- Line 223, main stack
- Line 236, Figure 1 instead of Figure 2
- Line 455, The numbers look comparable between DCNN-analysed and Briggs-predicted plume rise for the average values and the numbers are also within uncertainty. Why does it indicate Briggs underestimate plume rise on average?
- Line 458-459, From Table 3, Briggs under predicts the plume rise for neutral condition and over predicts the plume height for stable conditions. It is unclear why the authors concluded that Briggs generally under predicted plume rise.
Citation: https://doi.org/10.5194/egusphere-2025-4582-RC2 -
RC3: 'Comment on egusphere-2025-4582', Riccardo Simionato, 13 Mar 2026
This study highlights that the Briggs plume rise equations struggle with high-resolution imagery and parameterisation. The authors use a Deep Convolutional Neural Network to generate masks from visible-wavelengths images and compare the obtained geometrical parameters with theoretical calculations from Briggs equations. Their contribution represents a good advancement in the study of this kind of plumes which are usually analysed through more classical in industrial applications. Their manuscript is to be considered close to publication standard once a couple of concerns are addressed.
Specific comments:
- The authors should clearly state at the beginning that they are correcting the heights they measure for the wind direction. Moreover, they only refer to Koushafar et al. (2023) for the wind correction, without explaining the precess with any equations. Also in that paper, if used as reference for the wind correction, there is a lack of information. I would suggest to check that the equations used by the authors correspond to those published in Snee et al. (2023; 10.30909/vol.06.02.447458), a paper in which the wind correction is applied to correct height and width measurements of volcanic plumes. For example, in lines 200-201 the authors stated that only the wind direction and the camera-smokestack distance are necessary for the wind correction, but normally also field of view, inclination and camera orientation are normally required to perform a complete wind correction. Also, at lines 686-688, a bent-over plume with wind blowing towards the camera normally gives a plume height calculated in the wind-corrected plane that is smaller than that calculated in the image plane (see Snee et al., 2023, for a description of the wind correction applied to plume-shaped objects and for some examples). Also, I am not sure to understand the purpose of lines 682-686.
- I would not use the wording "real-life" referred to the position (lines 199, ). It is more a corrected projection on a plane corresponding to the wind direction than its actual position (as the wind direction can have uncertainties and the reconstruction is still in 2D). Moreover, I think Figure 3 might be misleading in this sense, as the plume is represented in 2D in birds-eye view, but this is not a real information obtained from the camera.
- Lines 463-464: If my understanding is correct, you are not comparing the same dataset (60 plumes vs only 26 plumes). I guess this might also play a role on this lower heights determined with exponential curve fitting. This needs at least to be stated, or, better, the comparison should be made between the same dataset.
- Lines 811-812 and 825-830: Dimensionless constants usually represent physical relationships. Reducing them by a factor of 2 or 6 implies that there is an important change in the dynamics. Also, this parameterisation is site-specific. All this might be at least stated in the manuscript, but has strong physical implications.
Minor comments:
- Line 79: I suggest to describe better the composition of these plumes. Are they only water vapour plumes? Do they contain other species?
- Lines 126-129: The sentence is maybe unnecessary, as it seems the authors are not correcting the plume rise measurements for the wind direction, which they are actually doing. I would consider erasing the sentence or already specify that the values are actually wind-corrected.
- Lines 183-184: I would use "plume position in the image plane to plume position in a wind-corrected plane", as this is still a planar correction (not a 3D reconstruction)
- Figure 2: I suggest adding the notations for plume rise height and distance on the images, as well as a white/black buffer around words to improve readability. In the caption, the notation of plume rise height seems to be wrong, while that of the distance is not in math format.
- Line 221: I would use "physical units" instead of "real-life distances", as it is about the conversion from pixel units to physical units (i.e., meters)
- Equation 4: Is dividing by L correct? Or is it just the notation that is misleading?
- Lines 287-288: Why data points are only available every three hours? What ERA5 product has been used? Also, ERA5 datasets normally have a citation.
- Lines 306-317: Could the authors motivate why they chose these particular formulations? Do other formulations exist for smoke plume rise and why they were not chosen?
- Line 357: Is it worth adding the variation for hour of day to the supplementary materials?
- Figure 5: I would add the camera orientation value on the figure, maybe also with a small arrow pointing towards that direction to show it effectively
- Line 395: I would change "plumes that are perpendicular to the camera orientation" for "having propagation axis that are close to the perpendicular to the camera orientation", as the plume has an extension and because also plumes twith axis hat are close to the perpendicular are used in the study.
- My general feeling is that a lot of tables could be merged (for example tables 2-5 and 6-9) to improve readability.
Technical corrections:
- Line 108: Missing point at the end of the sentence
- Lines 338-340: Maybe change the structure of the sentence for readability
- Lines 341, 343, 345: I would move "plume rise becuase of momentum" before the bullet point for clarity
- Lines 448-451: This sentence is redundant, as it is the same processing applied in section 3.2.1, consider removing it
- Lines 572-576: This sentence is unnecessary, as you could just mention that with the applicaiton of the plume rise distance the remaining images are 1514 (from the 3785 obtained in section 3.3.1 with the other filters)
- Line 647: missing comma after "average"
- Line 779: Missing "and" between "dawn" and "dusk"?
- Line 780: There is an extra point after "hour"
Citation: https://doi.org/10.5194/egusphere-2025-4582-RC3
Data sets
Smokestack Plume Images and Plume Identification Masks M. Gordon et al. https://doi.org/10.20383/103.01448
Model code and software
Code and data for the calculation of plume rise M. Gordon et al. https://doi.org/10.5683/SP3/WZVZBV
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 416 | 174 | 42 | 632 | 22 | 27 |
- HTML: 416
- PDF: 174
- XML: 42
- Total: 632
- BibTeX: 22
- EndNote: 27
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
The paper provides a comparison of the Briggs plume rise parameterisations against a deep convolutional neural network for determining plume rise from imagery (DPRNet).
I felt that there was a lack of understanding of the Briggs formulae. The underlying basis of the Briggs plume rise equations are the governing conservation equations of mass, momentum and heat. Solving these governing equations under specific conditions (buoyancy or momentum dominated, or particular meteorological conditions e.g., zero crosswind, constant buoyancy frequency) leads to the Briggs formulations. The ‘dimensionless constants’ are entrainment parameters. The Briggs formulae are commonly used for predicting plume rise but there are reasons why they may not work well, for example if the assumptions made in deriving the Briggs formulae from the underlying conservation equations do not hold true (e.g., the true atmospheric meteorological profiles may differ from that assumed). Indeed, the authors do make use of a modification for the interaction with the boundary layer top. Other authors have made direct use of the underlying conservation equations, considering the true atmospheric and release conditions (Webster and Thomson, 2002) and more elaborate models do also account for latent heat release from moisture in the rising plume / entrained atmospheric air (Fathi et al., 2025).
I found the text a bit vague in places, a little repetitive and verbose in other places and with many references to supplementary information and over-/under-predictions. It is worth the authors considering what key points / results they would like to present in the main paper and being both concise and precise. In addition, there are quite a lot of typographical errors which, with more care and attention, could have been avoided. All in all, this makes the paper quite hard to follow. Some further details are given in the detailed points below. The conclusion section was, however, well written and provided a good summary.
How are the Briggs formulae for momentum and buoyancy applied to calculate plume rise for plumes with both momentum and buoyancy? Are they just added together and, if so, is this appropriate? The calculation of the proportion of the plume rise due to momentum in section 3.2.3 assumes this simple addition but I’m not sure this is true in reality. Understanding the underlying conservation equations and the derivation of the Briggs equations (and the assumptions made) may shed some light on this.
Presumably DPRNet has been trained on an earlier dataset? What dataset was this and will the trained model be applicable to the images processed here? What about uncertainties / errors in the observations? Line of sight and plume direction are mentioned. Indeed, a sensitivity analysis of the wind direction on the determined observed plume height is conducted. Two methods for obtaining the plume rise from the observations are used, but conclusions are drawn on the performance of the Briggs formulae to these observations without consideration of uncertainties in the observations. Indeed, the authors state under- or overpredictions which are small compared to the uncertainties, say, in the observed plume rise height due to the wind direction presented in section 3.5.1.
Is Figure 3 in the horizontal plane? Is this a valid assumption? I can imagine that the plume height may not be at the same height as the camera. The caption mentions a similar transformation in the vertical plane to determine the plume rise height. A reference is given but it would be helpful to give more detail here on the calculation, what is measured (presumably SP’), what is calculated (presumably SP) and how.
Minor points:
Some minor points on the text are listed here: