the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Applying pySTEPS optical flow algorithms to improve convection nowcasting over the Maritime Continent
Joseph Albert Smith
Cathryn Birch
John Marsham
Simon Peatman
Massimo Bollasina
George Pankiewicz
Abstract. The Maritime Continent (MC) regularly experiences powerful convective storms that produce intense rainfall, flooding and landslides, which numerical weather prediction models struggle to forecast. Nowcasting uses observations to make more accurate predictions of convective activity over short timescales (~0–6 hours). Optical flow algorithms are effective nowcasting methods as they are able to accurately track clouds across observed image series and predict forward trajectories. Optical flow is generally applied to weather radar observations, however, the radar coverage network over the MC is not complete and the signal cannot penetrate the high mountainous regions. In this research, we apply optical flow algorithms from the pySTEPS nowcasting library to satellite imagery to generate both deterministic and probabilistic nowcasts over the MC. The deterministic algorithm shows skill up to 4 hours on spatial scales of 10 km and coarser, and outperforms a persistence forecast for all lead times. Lowest skill is observed over the mountainous regions during the early afternoon and highest skill is seen during the night over the sea. A key feature of the probabilistic algorithm is its attempt to reduce uncertainty in the lifetime of small scale convection. Composite analysis of 3-hour lead time nowcasts, initialised in the morning and afternoon, show it produces reliable ensembles but with an under-dispersive distribution, and produced area under the curve scores (i.e. ratio of hit rate to false alarm rate across all probability thresholds) of 0.80 and 0.71 over the sea and land, respectively. When directly comparing the two approaches, the probabilistic nowcast shows greater skill at ≤ 60 km spatial scales, whereas the deterministic nowcast shows greater skill at larger spatial scales ~200 km. Overall, the results show promise for the use of pySTEPS and satellite retrievals as an operational nowcasting tool over the MC.
- Preprint
(2670 KB) - Metadata XML
- BibTeX
- EndNote
Joseph Albert Smith et al.
Status: open (until 30 Oct 2023)
-
RC1: 'Comment on egusphere-2023-1404', Anonymous Referee #1, 02 Aug 2023
reply
The paper presents the results of the application and validation of the pySTEPS software package on the nowcasting of convective clouds. Papers of that type are quite useful for the scientific community because they provide the necessary insights regarding the efficiency of the package and highlight innate deficiencies and issues to be addressed.
The manuscript is well structured, does not present linguistic problems and summarizes the existing related literature at an acceptable level.
Overall, I recommend publication after minor revision. My specific and technical comments follow.
Specific comments
Title: to my opinion, the paper does not improve over an existing methodology. It presents the results of the application of the software package on the convection nowcasting over an area. Also it is important to highlight the fact that satellite data are used and not radar data. Therefore the title could be changed to something similar to “Applying pySTEPS optical flow algorithms to the convection nowcasting over the Maritime Continent using satellite data”.
Keywords: The area of application should be added. Also “satellite data” and/or “Himawari satellite”
Table 1. The NWC SAF EUMETSAT Facility could be mentioned.
Line 323: Please define clearly the term “persistence nowcast”.
Line 355: Is the term “persistence forecast” interchanged randomly in the text with the term “persistence nowcast”? If no, please define the term “persistence forecast”. In general some confusion with the terms “LK nowcast”, “persistence nowcast” and “persistence forecast” exists. Please clarify and use with consistency throughout the text.
Lines 396-397: Could you please clarify how the extension of the evaluation period could smooth out the noise over sea?
Lines 398-400: Could the better scores during overnight and morning hours compared to the rest of the times (especially over land) be attributed to the limited convective activity that minimizes the erroneous propagations too?
Lines 509-510: You state at some point that “…the number of hits continues to exceed the number of misses….”. However neither POD nor POFD directly suggests that. Overall the ROC curve is a measure the forecast’s efficiency to discriminate between the events that actually happened and those that did not. Loosely explained, one could think that it shows how good is a forecast in “finding” the upcoming events without over-forecasting them. Please elaborate on or rephrase your conclusion.
A general comment: Are Himawari BTs offered in a reprojected grid of specific horizontal resolution? If not, it would be beneficial to offer some information on how the satellite data are reprojected, since the skill analysis if provided in specific resolution scales.
Technical comments:
Line 34: I would prefer “…The “Early Warnings For All” initiative…” instead of “…The Early Warnings For All Initiative…”
Line 24: ” …show to…” could be removed.
Line 138: I would prefer “…ranging from 0.47 μm to 13.3 μm…” instead of “…ranging from 0.47 μm–13.3 μm…”.
Line 141: I think that sentence “Convective clouds can be clearly identified…” is more accurate than “Clouds can be clearly identified…”.
Line 151-152: Although not necessary, another map presenting the inner and outer domain could be useful.
Figure 4. Could the caption be changed to “Composite FSSs against lead time for 3,457 LK nowcasts (solid line) and…”?
Line 437: I think that “…time be will be…” should change to “time will be..”.
Line 440-441: Could it be “…Less growth results to higher skill…” rather than “…Less growth results higher skill…”?
Line 504: “…For example, At T+1…” should read “…For example, at T+1…”.
Figure 8: I would suggest to change the horizontal axis legend (and all related text) of the ROC graph to “Probability Of False Detection (POFD)” to avoid confusion with the False alarm Ratio (FAR).
Figure 8: I would strongly suggest to use the same scale for the two axes of the ROC and the reliability graphs (that will result in two “square” graphs).
Line 560: I would use the word “presents” instead of “predicts”.
Citation: https://doi.org/10.5194/egusphere-2023-1404-RC1
Joseph Albert Smith et al.
Joseph Albert Smith et al.
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
160 | 69 | 9 | 238 | 5 | 4 |
- HTML: 160
- PDF: 69
- XML: 9
- Total: 238
- BibTeX: 5
- EndNote: 4
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1