Next Article in Journal
Developing a System for Integrated Environmental Information in Urban Areas: An Estimation of the Impact of Thermal Stress on Health
Previous Article in Journal
Lightning Potential Forecast Evaluation and Its Correlation with Thermodynamic Indices
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Proceeding Paper

An Ensemble Forecasting Alternative Based on Stochastic Parameter Perturbation (SPP) on Potential Vorticity Anomalies through the Identification of Weather Features as Coherent Objects †

by
Platon Patlakas
1,*,
Diego S. Carrió
2 and
Emmanouil Flaounas
3
1
Department of Physics, National and Kapodistrian University of Athens, 15784 Athens, Greece
2
Department of Physics, Universitat de les Illes Balears, 07122 Palma, Spain
3
Hellenic Center for Marine Research, 19013 Anavyssos, Greece
*
Author to whom correspondence should be addressed.
Presented at the 16th International Conference on Meteorology, Climatology and Atmospheric Physics—COMECAP 2023, Athens, Greece, 25–29 September 2023.
Published: 29 August 2023

Abstract

:
The Pacific Ocean witnesses frequent cyclonic activity. The destructive impact of these storms, including strong winds, heavy rain, and storm surge, causes flooding, landslides, and extensive damage. Understanding cyclone genesis and evolution is crucial for accurate forecasts and minimizing harm. Towards this direction, an alternative ensemble forecasting approach based on a stochastic parameter perturbation (SPP) scheme, applied in potential vorticity (PV) anomalies, was developed. Testing it on Typhoon Usagi demonstrated its effectiveness in introducing uncertainties to storm tracks and cyclone development. These findings highlight the potential of stochastic methods in regional forecasting systems.

1. Introduction

Super Typhoon Usagi was a very intense cyclone (equivalent to a category 4 hurricane on the Saffir–Simpson hurricane wind scale) developed in the Western Pacific Ocean on 16 September 2013. Its development was influenced by a combination of favorable atmospheric conditions, including warm sea surface temperatures and low vertical wind shear [1]. The cyclone exhibited a well-defined eye at its center, surrounded by concentric bands of intense thunderstorms. Analysis of satellite imagery and meteorological data revealed that Cyclone Usagi underwent rapid intensification, with the minimum sea level pressure reaching 910 hPa [2]. This powerful tropical cyclone exhibited characteristics of a mature system, featuring sustained wind speeds of up 205 km/h (according to the Japan Meteorological Agency—JMA).
Cyclone Usagi had significant impacts on several countries in its path, primarily affecting the coastal regions of the Philippines, Taiwan, and southern China. It triggered extensive flooding, landslides, and infrastructure damage, leading to 39 deaths and significant economic losses estimated at approximately USD 4.32 billion [3,4]. These impacts underscore the importance of preparedness measures, early warning systems, and resilient infrastructure in vulnerable coastal regions to mitigate the devastating effects of tropical cyclones.
In this way, ensemble forecasts are employed in order to quantify the uncertainties in cyclone paths, dynamics and impacts. The present study is an effort to propose an alternative way of producing model ensembles based on stochastic parameter perturbations (SPP) on potential vorticity anomalies through the identification of weather features as coherent objects.

2. Model Set-Up

The numerical simulations presented in this study are performed using the Advanced Weather Research and Forecasting Model (WRF-ARW, version 4.2.2, [5]). The domain is set up with horizontal grid resolution of 4 km and a hybrid 61 terrain-following η levels up to 50 hPa. Initial and boundary conditions were obtained from hourly ERA5 reanalysis [6] at 0.25° grid spacing.
The specific physical parameterization schemes common to all the simulations performed in this study are summarized in the following table (Table 1):

3. A Feature-Based Stochastic Scheme (FBS)

The future-based stochastic system (FBS) aims to perturb stochastically the grid points that dynamically describe a cyclone system. This is carried out via a four-step procedure:
Step one—PV budget calculation: A module has been developed that calculates the non-conserved PV components of the total atmospheric PV at the beginning of every model time step. This module is described in detail in [13] where it has been used to analyse the processes that contribute to the intensification of Mediterranean cyclones.
Step two—Identifying and tracking objects: A new module has been developed and implemented into WRF to identify coherent 3D objects. Each object is composed of neighboring grid points of PVdiab or PVmo that exceed the absolute value of 0.75 PVUs. However, we retrain only objects that include at least one grid point of more than 2 PVUs. From the perspective of PV invertibility, these two absolute value thresholds are deemed adequate for retaining objects which describe meso-scale systems in terms of size, and have a significant impact on the atmospheric state. Since this study focuses on cyclones, we included an additional criterion that demands from objects to be composed of grid points with negative pressure perturbation (P’, expressed in pressure anomalies from model level averages) (Figure 1).
Step three—tracking objects in time: Once identified each object is separately labeled according to the time step it was identified. If there are overlapping objects, then the oldest time label is assigned to the object.
Step four—assigning a perturbation coefficient: Finally, every object is assigned to a coefficient ct that changes in time according to the following equation:
c t = 1 d t τ c t 1 + 0.5 χ 1 1 d t τ 2
If ct < −1, then ct = −2 − ct, if ct ≥ 1, then ct = 2 − ct. χ is a random number that ranges from −1 to 1, t is the time step, dt is the model time step and τ is a constant in units of time. The choice of τ is arbitrary but nevertheless it is crucial for the frequency of changes of the perturbation coefficient. As an example, Figure 2 shows examples of the time evolution of the perturbation coefficient ct at every model time step for τ = 12 h.

4. Results

In terms of the trajectory of Usagi, we can observe that the control simulation shows a similar track compared with the ones obtained from the JMA, specially during the initiation and mature stage of Usagi (Figure 3). Although the control trajectory starts to diverge in the dissipation phase of Usagi, we can certainly conclude that the control simulation performs with accuracy enough the typhoon trajectory.
To assess the model’s sensitivity in cyclone forecasting, six simulations were conducted using the Stochastically Perturbed Parametrization Tendencies (SPPT) scheme within the WRF model, while maintaining the same model configuration. The ensemble cyclone tracks exhibited a close resemblance to the reference track (Figure 4), indicating a comparable spread. Similarly, the development of the cyclonic system, as reflected in the Mean Sea Level Pressure (MSLP) values at its center, varied around those of the control simulation.
Before implementing the perturbation coefficient on the physical tendencies of the objects, a set of experiments was performed by multiplying them with constant values. These values ranged from 0 to 2, with an increment of 0.25. This allowed for a deeper understanding of the impacts this procedure had on the system’s evolution. The findings revealed that coefficients smaller than 1 (where 1 represents the control simulation) had a more pronounced effect compared to larger coefficients (Figure 5).
For the purposes of this study, the FSB scheme utilized random coefficients within the range of −0.1 to +0.4. The application of this scheme yielded a satisfactory spread, albeit narrower than that observed with the traditional SPPT methodology. The cyclone exhibited sensitivity throughout all stages, with the minimum mean sea level pressure (MSLP) value consistently higher than the control case in most instances (Figure 6).

5. Conclusions

In this study, we present preliminary results on a new SPPT scheme where perturbations are uniquely applied to areas characterized by high PV, produced by diabatic processes. Therefore, we only perturb the grid points which are expected to have a strong impact on the component of the atmospheric state that is sensitive to inherent model uncertainties. Our results show comparable spread on the tracks and MSLP evolution to the one produced by the original SPPT method. In contrast to the original SPPT method, we consider the fact that this approach is based on PV theory and, therefore, that perturbations have a physical basis, to be an advantage

Author Contributions

Conceptualization, E.F.; Data curation and Investigation, P.P., D.S.C. and E.F.; Methodology, P.P., D.S.C. and E.F.; Resources, and Software, P.P. and D.S.C.; Supervision, E.F.; Visualization, D.S.C. and E.F.; Writing—original draft, P.P., D.S.C. and E.F.; Writing—review and editing, P.P. and E.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The ERA-5 Reanalysis are openly available via Copernicus (www.copernicus.eu, accessed on 23 May 2023).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, S.; Tao, D.; Zhao, K.; Minamide, M.; Zhang, F. Dynamics and predictability of the rapid intensification of Super Typhoon Usagi (2013). J. Geophys. Res. Atmos. 2018, 123, 7462–7481. [Google Scholar] [CrossRef]
  2. NCICS.ORG. Available online: https://ncics.org/ibtracs/index.php?name=v04r00-2013259N17132 (accessed on 24 May 2023).
  3. Final Report: Effects of Typhoon ODETTE (USAGI). NDRRMC. 2014. Available online: https://reliefweb.int/report/philippines/ndrrmc-update-sitrep-no-4-re-effects-typhoon-odette-usagi (accessed on 23 May 2023).
  4. China Meteorological Agency. Member Report: China (PDF). ESCAP/WMO Typhoon Committee: 8th Integrated Workshop/2nd TRCG Forum. ESCAP/WMO Typhoon Committee. p. 16. Archived (PDF) from the Original on 3 December 2013. Available online: https://www.typhooncommittee.org/8IWS_2TRCG/docs/Members%20Report/2013MemberReportChina.pdf (accessed on 23 May 2023).
  5. Skamarock, W.C.; Klemp, J.B.; Dudhia, J.; Gill, D.O.; Barker, D.M.; Duda, M.; Huang, X.Y.; Wang, W.; Powers, J.G. A Description of the Advanced Research WRF Version3, 603NCAR; Tech. Note; Mesoscale and Microscale Meteorology Division; National Center for Atmospheric Research: Boulder, CO, USA, 2008. [Google Scholar]
  6. Hersbach, H.; Bell, B.; Berrisford, P.; Hirahara, S.; Horányi, A.; Muñoz-Sabater, J.; Nicolas, J.; Peubey, C.; Radu, R.; Schepers, D.; et al. The ERA5 global reanalysis. Q. J. R. Meteorol. Soc. 2020, 146, 1999–2049. [Google Scholar] [CrossRef]
  7. Thompson, G.; Field, P.R.; Rasmussen, R.M.; Hall, W.D. Explicit forecasts of winter precipitation using an improved bulk microphysics scheme. Part II: Implementation of a new snow parameterization. Mon. Weather Rev. 2008, 136, 5095–5115. [Google Scholar] [CrossRef]
  8. Thompson, G.; Rasmussen, R.M.; Manning, K. Explicit forecasts of winter precipitation using an improved bulk microphysics scheme. Part I: Description and sensitivity analysis. Mon. Weather Rev. 2004, 132, 519–542. [Google Scholar] [CrossRef]
  9. Kain, J.S. The Kain–Fritsch convective parameterization: An update. J. Appl. Meteorol. 2004, 43, 170–181. [Google Scholar] [CrossRef]
  10. Iacono, M.J.; Delamere, J.S.; Mlawer, E.J.; Shephard, M.W.; Clough, S.A.; Collins, W.D. Radiative forcing by long-lived greenhouse gases: Calculations with the AER radiative transfer models. J. Geophys. Res. Atmos. 2008, 113, D13. [Google Scholar] [CrossRef]
  11. Hong, S.Y.; Noh, Y.; Dudhia, J. A new vertical diffusion package with an explicit treatment of entrainment processes. Mon. Weather Rev. 2006, 134, 2318–2341. [Google Scholar] [CrossRef]
  12. Hong, S.Y. A new stable boundary-layer mixing scheme and its impact on the simulated East Asian summer monsoon. Q. J. R. Meteorol. Soc. 2010, 136, 1481–1496. [Google Scholar] [CrossRef]
  13. Flaounas, E.; Gray, S.L.; Teubler, F. A process-based anatomy of Mediterranean cyclones: From baroclinic lows to tropical-like systems. Weather Clim. Dynam. 2021, 2, 255–279. [Google Scholar] [CrossRef]
Figure 1. An example for Identifying and tracking objects—Cyclone Usagi.
Figure 1. An example for Identifying and tracking objects—Cyclone Usagi.
Environsciproc 26 00110 g001
Figure 2. Five-day time evolution of the perturbation coefficient for characteristic length (τ in Equation (1)) equal to 12 h.
Figure 2. Five-day time evolution of the perturbation coefficient for characteristic length (τ in Equation (1)) equal to 12 h.
Environsciproc 26 00110 g002
Figure 3. Usagi’s trajectories depicted by the CNTRL simulation and the JMA from 16 UTC 9 September to 21 UTC 23 September 2013. CNTRL and JMA data are depicted every 3 h.
Figure 3. Usagi’s trajectories depicted by the CNTRL simulation and the JMA from 16 UTC 9 September to 21 UTC 23 September 2013. CNTRL and JMA data are depicted every 3 h.
Environsciproc 26 00110 g003
Figure 4. SPPT ensemble member cyclone tracks and minimum MSLP evolution.
Figure 4. SPPT ensemble member cyclone tracks and minimum MSLP evolution.
Environsciproc 26 00110 g004
Figure 5. Cyclone tracks and minimum MSLP evolution during the experiments where the identified objects are multiplied with values ranging from 0 to 2 (0% to 200%).
Figure 5. Cyclone tracks and minimum MSLP evolution during the experiments where the identified objects are multiplied with values ranging from 0 to 2 (0% to 200%).
Environsciproc 26 00110 g005
Figure 6. Cyclone tracks and minimum MSLP evolution for all ensemble members produced through the proposed FSB approach.
Figure 6. Cyclone tracks and minimum MSLP evolution for all ensemble members produced through the proposed FSB approach.
Environsciproc 26 00110 g006
Table 1. WRF parameterizations used for the study.
Table 1. WRF parameterizations used for the study.
Microphysicsthe single-moment Thompson microphysics scheme [7,8]
Cumulus ParameterizationKain–Fritsch [9]
Short- and long-wave radiation physicsRRTMG scheme [10]
Planet boundary layernon-local K Yonsei University scheme [11,12]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Patlakas, P.; Carrió, D.S.; Flaounas, E. An Ensemble Forecasting Alternative Based on Stochastic Parameter Perturbation (SPP) on Potential Vorticity Anomalies through the Identification of Weather Features as Coherent Objects. Environ. Sci. Proc. 2023, 26, 110. https://0-doi-org.brum.beds.ac.uk/10.3390/environsciproc2023026110

AMA Style

Patlakas P, Carrió DS, Flaounas E. An Ensemble Forecasting Alternative Based on Stochastic Parameter Perturbation (SPP) on Potential Vorticity Anomalies through the Identification of Weather Features as Coherent Objects. Environmental Sciences Proceedings. 2023; 26(1):110. https://0-doi-org.brum.beds.ac.uk/10.3390/environsciproc2023026110

Chicago/Turabian Style

Patlakas, Platon, Diego S. Carrió, and Emmanouil Flaounas. 2023. "An Ensemble Forecasting Alternative Based on Stochastic Parameter Perturbation (SPP) on Potential Vorticity Anomalies through the Identification of Weather Features as Coherent Objects" Environmental Sciences Proceedings 26, no. 1: 110. https://0-doi-org.brum.beds.ac.uk/10.3390/environsciproc2023026110

Article Metrics

Back to TopTop