SciELO - Scientific Electronic Library Online

 
vol.26 número3Machine Translation for Low-Resource English-Mizo Pair Encountering Tonal WordsFuzzy Distribution Sets índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Computación y Sistemas

versión On-line ISSN 2007-9737versión impresa ISSN 1405-5546

Comp. y Sist. vol.26 no.3 Ciudad de México jul./sep. 2022  Epub 02-Dic-2022

https://doi.org/10.13053/cys-26-3-4359 

Articles

Tropical Cyclone Simulations with WRF Using High Performance Computing

Jorge Clouthier-López1 

Ricardo Barrón-Fernández1  * 

David Alberto Salas-de-León2 

11 Instituto Politécnico Nacional, Centro de Investigación en Computación, Mexico. clouthier@gmail.com.

22 Universidad Nacional Autónoma de México, Instituto de Ciencias del Mar y Limnología, Mexico. dsalas@unam.mx.


Abstract:

Tropical cyclone (TC) Bud occurred over the Eastern Pacific Ocean from 9 June to 16 June 2018 reaching, a category-4 hurricane status (H4) due to its strong sustained winds and gusts. In this study, we simulate this TC to reproduce its track path, direction, and strength to determine the best model physics configurations that weather agencies could use to forecast TC tracks over the Eastern Pacific Ocean, adjacent to the Mexican coastline. To achieve this goal, the sensitivity for the impact of different microphysics and cumulus parameterization schemes is carried out through high-performance computer simulations, with the WRF model using the cluster CÓDICE B2 at Centro de Investigación en Computación (CIC) of the Instituto Politécnico Nacional (IPN). The realism of the TC for the different schemes is assessed by comparing the simulations and the best track data taken from the National Hurricane Center (NHC-NOAA). The NCEP-GFS forecast data is used as initial and boundary conditions. The evolution of wind and minimum pressure at sea level for the different physics combination runs are also compared with the best track data. We found that the track paths and intensities improve when Sea Surface Temperature (SST) is allowed to evolve with the modeled atmosphere via computer simulations.

Keywords: Numerical weather prediction (NWP); weather mesoscale computer modeling; Tropical cyclone

1 Introduction

Tropical cyclones (TCs) are weather-rotating systems that form associated with thunderstorms. These systems are characterized by a low-pressure center and extreme wind velocities [1, 2]. They go through four stages, from lower to higher intensity: tropical disturbance, tropical depression, tropical storm, and hurricane. TCs originate over tropical or subtropical ocean regions with sea surface temperatures (SSTs) around 27 °C.

These weather systems can cause deaths and considerable damage over coastal areas [3, 4, 5]. When TCs occur near the coast, like TC Bud, they can indirectly impact over distant inland areas. When they occur adjacent to the continent, convective storms develop over distant inland areas. According to [6], in the tropical Americas, a single TC can cause intense precipitation over inland arid regions.

Therefore, it is of fundamental importance to track down TCs well ahead to minimize societal damage and economic impacts. In general, numerical weather prediction (NWP) models, that run over supercomputers, are used to forecast the paths and strengths of TCs. The latter are highly interrelated with other atmospheric phenomena and strongly dependent on the variations of the SST, which depends on ocean feedbacks between the atmosphere and the ocean mixing layer.

NWP models use initial and boundary conditions to solve atmospheric dynamic and thermodynamic systems of partial differential equations with the assumptions of some semi-empirical physical approximations, known as parameterizations or parameterization schemes. The latter are linked with sub-grid scale processes that cannot be resolved explicitly by the governing equations implemented in the computer models.

The simulation of TCs with NWP models is very sensitive to the parameterizations used during the model runs. Several studies on the impact of parameterization schemes, when simulating mesoscale weather events, can be found in the literature.

It has been established that not only the physics options [7] but also the initial and boundary conditions in the initialization of NWP models [8, 9, 10, 11] play a decisive role in weather simulations. NWP models are currently used for both research and operational forecasting by different research and weather forecasting agencies around the world.

A sensitivity study is the most adequate to quantify the effect of different physics options; and, in turn, to predict the track path of a TC and its strength over a determined geographical region. Previous studies have made an effort to identify the optimum parameterization schemes and customize NWP models for the simulation and forecast of TCs over specific regions of the world.

For example, [12] documented a comprehensive assessment of the performance of parameterization schemes for TC computer simulations.

[13] performed a sensitivity study of the physics parameterizations when simulating the TC Jal over the Indian Ocean.

[14] performed sensitivity studies to evaluate cumulus parameterization schemes.

[15] conducted several sensitivity experiments for five TCs over the Bay of Bengal and found that the combination of the Kain-Fritsch for convection, the Yonsei University planetary boundary layer (PBL), the LIN for microphysics, and the NOAH for land surface schemes had the best performance to reproduce both the tracks and intensities.

Additionally, data assimilation has also been applied to provide data observations during the running of computer simulations to nudge the development of TC track paths, being considered a successful approach for improving TC forecasts [16, 17, 18]. However, this technique does not let the physical processes to evolve naturally. Besides, data assimilation techniques are very expensive, and they are not suitable to be implemented for forecasting purposes due to the high demand of computer resources.

In contrast to previous studies, we allowed the SST to evolve with the modeled atmosphere, for each model configuration, via the evolution of the computer simulations. According to the ocean temperature feedback between the modeled atmosphere and the ocean mixing layer, the purpose is to get reliable performance when the SST changes are taken into account in the development of TC Bud.

In this article, two ensembles of simulations with various parameterization schemes were carried out to define its members, from 00:00 UTC 10 June to 00:00 UTC 16 June 2018. A sensitivity study for a set of cumulus and microphysics parameterizations schemes was conducted to evaluate the track path and intensity of TC Bud. Three microphysics schemes, three cumulus schemes, and a moisture-advection-based trigger scheme were assessed. In one ensemble, the SST was fixed during the simulation period, while in the other one it was allowed to evolve with the modeled atmosphere via the computer simulations. The latter ensemble showed the best results.

The Weather Research and Forecasting (WRF-ARW) model was used with a new hybrid terrain-following sigma-pressure vertical coordinate. A parent domain and two child nested domains were established to define the grids over which the governing equations are solved. In the innermost domain, the assessment was performed with a grid spacing of 2 km.

This article is organized as follows. Section 2 describes the case study. Section 3 describes the model setup and methodology. Section 4 presents the results. Finally, in section 5, the conclusion is presented.

2 Case Study

TC Bud indirectly impacted the Mexican continental territory; when it occurred over the Eastern Pacific Ocean, heavy precipitation and flooding were present over inland Mexico. In Fig. 1a, the satellite spatial distribution of precipitation is observed. Fig. 1a was obtained and produced with the Giovanni online data system, developed and maintained by the NASA GES DISC.

Fig. 1 (a) Spatial distribution of observed averaged precipitation by satellite, at 02:00 UTC, on 12 June, produced with the Giovanni online data system, developed and maintained by the NASA GES DISC. (b) TC Bud’s best track (in black) and (in red) the track from the Global Forecasting System (GFS) 

This TC occurred from 18:00 UTC 9 June to 06:00 UTC 16 June 2018, adjacent to the Mexican coastline, over the Eastern Pacific Ocean. In Fig. 1b, the best track (in black) is shown. Also, in Fig. 1b, the track obtained from the Global Forecast System (GFS) (in red) is also shown. The precursor of the TC was an easterly wave that emerged off the west coast of Africa on 29 May 2018 then entered the eastern North Pacific Ocean on 6 June. Later, a low pressure formed on 9 June, and it started to rotate.

Consequently, a tropical depression, 528 km south of Acapulco Guerrero, formed at 18:00 UTC. It became a tropical storm on 10 June at 00:00 UTC, reaching hurricane status on 10 June at 18:00 UTC. The maximum wind speeds reached 61.7 m/s at 00:00 UTC on 12 June near Manzanillo. Later, Bud became a tropical storm on 13 June at noon UTC, landfalling over Baja California Sur on 15 June at 02:00 UTC, with wind speeds of 20.5 m/s and a central pressure of 999 hPa. Subsequently, it crossed Baja California Sur, weakening due to the interaction with land.

The TC converted into a convection-free post-tropical cyclone on 15 June at noon UTC, and the isobaric contours opened by 00:00 UTC on 16. Then Bud dissipated by 06:00 UTC on 16 June. In Fig. 2, a representation of part of the evolution of the TC Bud is shown. In this figure, the maximum wind speed (MWS) and the sea level pressure (SLP), obtained from one computer simulation, are shown.

Fig. 2 Evolution of TC Bud from 06:00 UTC on June 10 to 06:00 UTC on June 12. The maximum wind speed (MWS) is at the top, and the sea level pressure (SLP) is at the bottom, according to a WRF-ARW simulation 

3 Modeling Setup and Methodology

The (WRF-ARW) model, version 4.1.3 [19], was used to simulate TC Bud. The model has a fully compressible and non-hydrostatic dynamic core. It allows to simulate and forecast mesoscale convective systems with many parameterization schemes. This mesoscale numerical model was developed for research and operational forecasting. The WRF-ARW model was implemented in the CÓDICE B2 cluster at Centro de Investigación en Computación (CIC) of the Instituto Politécnico Nacional (IPN), using 16 nodes. Each node has 28 CPU cores.

In each of the two ensembles, the outer domain was set with a grid spacing of 18 km, while the innermost domain with a grid spacing of 2 km, using a ratio of 1 to 3 and two-way nesting.

In order to assess the different cumulus and microphysics parameterizations at 2 km grid spacing (over the innermost domain, d03, shown in Fig. 3), the initial and boundary conditions (IBCs) were downscaled. In Fig. 3, the nested domains d01, d02, and d03 with 18, 6, and 2 km grid spacing, respectively, are shown.

Fig. 3 Domain setup for the study. The assessment was carried over the innermost domain, d03, with a grid spacing of 2 km 

The IBCs were obtained from the National Center for Environmental Prediction (NCEP) FNL operational global analysis and forecast dataset, which is available with a grid spacing of 0.25° at 6 h intervals.

In this dataset, the files that contain the IBCs are made with the same model used in the Global Forecast System (GFS). To obtain more inputs, the data was interpolated at intervals of 3 h and used to update the outermost domain during each simulation, every 3 h. The innermost domain, covering the MTC region and the adjacent Pacific Ocean, is comprehended by a horizontal grid of 886 × 886 points. The model was configured with 51 vertical levels, using the new hybrid sigma-pressure vertical coordinate that follows the terrain and gradually makes the transition to constant pressure surfaces, reducing the numerical noise at higher levels [20, 21].

A 3 s time-step for the innermost domain was selected via the evolution of the simulations. A time-step ratio of three was chosen for the middle and outermost domains. In comparison to greater time-steps, a short time-step ensures a better simulation performance in terms of satisfying the Courant-Friedrich-Lewy condition [22]. The model needed to be warmed up (spin up) 17 h before analysis since simulation results are sensitive to the spin-up time [23].

Three different microphysics and cumulus schemes, with and without convection triggering, were varied one at a time–while the other schemes were left unchanged in each ensemble. The microphysics schemes used were the Purdue Lin scheme [24], the WRF single-moment 6-class scheme [25], and the Thompson scheme [26]. The cumulus schemes used were the Kain-Fritsch scheme [27], the Betts-Miller-Janjic scheme [28], and the simplified Arakawa-Schubert scheme [29]. The unchanged schemes were the Yonsei University scheme [30] for the planetary boundary layer, the MM5 similarity scheme [31, 32, 33, 34, 35] for the surface layer, the unified Noah land surface model [36] for the land surface, the Dudhia shortwave scheme [37] for the shortwave radiation, and the RRTM Longwave scheme [38] for the longwave radiation. A moisture advection-based trigger scheme was also considered.

In the first ensemble, the SST was provided, from the Real-time Global Sea Surface Temperature Analysis product (RTG-SST), for the period simulated over the ocean surface. In the second ensemble, just the IBCs from NCEP FNL were used, and the names for each ensemble member will end in “-GFS.” Each ensemble member of the first ensemble, which performed better than the second, is listed in Table 1. During the numerical simulations, the SST was modified according to the response of the surface winds and changes in radiative fluxes using the scheme proposed by [39]. The change in SST on a daily scale is not considerable; however, subtle changes could indirectly impact the simulated atmosphere. During TCs, surface winds cause ocean mixing over the first few meters, changing the SSTs. According to [40], negative feedback from wind-driven ocean mixing can be reflected in numerical simulations.

Table 1 Ensemble members of the ensemble with RTG-SST forcing 

Simulation Simulation acronym Microphysics parameterization Cumulus parameterization
1 LIN-KF-RTG Purdue Lin scheme Kain-Fritsch scheme
2 LIN-BMJ-RTG Purdue Lin scheme Betts-Miller-Janjic scheme
3 LIN-SAS-RTG Purdue Lin scheme Simplified Arakawa-Schubert Scheme
4 LIN-KF-TR-RTG Purdue Lin scheme Kain-Fritsch scheme Moisture advection-based trigger scheme
5 THOM-KF-RTG Thompson scheme Kain-Fritsch scheme
6 THOM-BMJ-RTG Thompson scheme Betts-Miller-Janjic scheme
7 THOM-SAS-RTG Thompson scheme Simplified Arakawa-Schubert scheme
8 THOM-KF-TR-RTG Thompson scheme Kain-Fritsch scheme Moisture advection-based trigger scheme
9 WSM6-KF-RTG WRF Single-moment 6-class scheme Kain-Fritsch scheme
10 WSM6-BMJ-RTG WRF Single-moment 6-class scheme Betts-Miller-Janjic scheme
11 WSM6-SAS-RTG WRF Single-moment 6-class scheme Simplified Arakawa-Schubert scheme
12 WSM6-KF-TR-RTG WRF Single-moment 6-class scheme Kain-Fritsch scheme Moisture advection-based trigger scheme

4 Results

Both cumulus and microphysics schemes play an important role when forecasting TCs. Cumulus schemes are associated with redistribution of moisture and heat, in the vertical direction, but they do not depend on the latent heat generated during condensation and deposition of water vapor. Besides, latent heat is important in the vertical development of convection inside TCs. During lifting, latent heat is released due to condensation or deposition. On the other hand, microphysics schemes depend on the moisture distribution in the troposphere, where heat transfer occurs between the ocean surface and the uppermost layers of the atmosphere.

The simulated tracks with and without SST forcing were compared with the best track data. The simulated paths closer to the best track path are those obtained from the ensemble members with SST forcing, which are discussed in this section. However, in Table 2, the mean track error for the ensemble members without SST forcing is shown. The tracking error was calculated as the distance between a point in the simulated track path and a point in the best track path along the corresponding great circle for all the ensemble members. Then, the average track error for each track was obtained.

Table 2 Mean track error (in km): ensembles without SST forcing 

LIN-KF-GFS 110.062
LIN-BMJ-GFS 83.624
LIN-SAS-GFS 44.6
LIN-KF-TR-GFS 44.551
THOM-KF-GFS 64.651
THOM-BMJ-GFS 95.147
THOM-KF-TR-GFS 35.337
THOM-SAS-GFS 46.612
WSM6-KF-GFS 79.583
WSM6-BMJ-GFS 76.741
WSM6-SAS-GFS 56.172
WSM6-KF-TR-GFS 34.471

Fig. 4a provides the simulated track paths for the WSM6-KF-RTG (with SST forcing) and WSM6-KF-GFS (without SST forcing) of TC

Fig. 4 (a) Track paths with longitude (deg) and latitude (deg) coordinates for WSM6-KF-RTG (with SST forcing) and WSM6-KF-GFS (without SST forcing) ensemble members. (b) Sea Level Pressure (SLP) for the THOM-KF-GFS, WSM6-KF-GFS, and LIN-KF-GFS ensemble members (without SST forcing). Best track observations are also included; as well as the track from the Global Forecasting System (GFS) in (a), shown in red The simulated track paths, the SLP along the tracks, and the MWS for the different cumulus and microphysics schemes, for the ensemble with SST forcing are presented from Fig. 5 to Fig. 7. Although the computer simulations can reproduce the best track path, they are not capable enough to predict the MWS 

Bud. It is observed that the ensemble member with SST forcing is closer to the best track path than the one without SST forcing. Besides, in Fig. 4b, the SLP for each of the THOM-KF-GFS, WSM6-KF-GFS, and LIN-KF-GFS ensemble members, along with the best track SLP (in black), is shown. The SST was not allowed to evolve with the modeled atmosphere with the Kain-Fritsch cumulus parameterization in the latter ensemble members.

It is noticed that the mentioned ensemble members overestimated the SLP of the TC. The same occurred for all other ensemble members of the ensemble without SST forcing. However, the simulated SLP improves when SST is allowed to evolve with the modeled atmosphere for most ensemble members with SST forcing.

This can be appreciated in Fig. 5b when the SLP improves compared to Fig. 4b.

Fig. 5 (a) Track paths with longitude (deg) and latitude (deg) coordinates for THOM-KF-RTG, WSM6-KF-RTG, and LIN-KF-RTG (with SST forcing) ensemble members. (b) Sea Level Pressure (SLP) for each of the ensemble members in (a). (c) Wind Speed, the maximum values of TC Bud (MWS), for each of the ensemble members in (a). Best track observations are also included; as well as the track from the Global Forecasting System (GFS) in (a), shown in red 

The simulated track paths, the SLP along the tracks, and the MWS for the different cumulus and microphysics schemes, for the ensemble with SST forcing are presented from Fig. 5 to Fig. 7. Although the computer simulations can reproduce the best track path, they are not capable enough to predict the MWS. The mean track errors for the THOM-KF-RTG and the WSM6-KF-RTG ensemble members, shown in Table 3, when the Kain-Fritsch cumulus scheme is used, are lower than the mean track error for the LIN-KF-RTG ensemble member. This is appreciated in the track paths shown in Fig. 5a. However, in the first hours of simulation, the track from THOM-KF-RTG considerably departs from the best track, but it presents the lowest mean track error along with the whole TC event.

Fig. 6 As in Fig. 5, but for THOM-KF-TR-RTG, WSM6-KF-TR-RTG, and LIN-KF-TR-RTG ensemble members, see Table 1 for acronyms 

Fig. 7 As in Fig. 5, but for THOM-BMJ-RTG, WSM6-BMJ-RTG, and LIN-BMJ-RTG ensemble members, see Table 1 for acronyms 

Table 3 Mean track error (in km) taking into account SST forcing 

THOM-KF-RTG 47.002
WSM6-KF-RTG 47.031
LIN-KF-RTG 70.041

In Table 4, the mean track error for the THOM-KF-TR-RTG, WSM6-KF-TR-RTG, and LIN-KF-TR-RTG ensemble members are shown. These ensemble members consider a moisture-advection-based trigger scheme for the Kain-Fritsch cumulus scheme. The latter considerably improved the mean track error for the three microphysics schemes.

Table 4 Mean track error (km) taking into account SST forcing 

THOM-KF-TR-RTG 32.207
WSM6-KF-TR-RTG 26.314
LIN-KF-TR-RTG 39.551

In these three ensemble members, with three different microphysics schemes, the THOM-KF-TR-RTG was the one that performed best, appreciated in Fig. 6a. However, the LIN-KF-TR-RTG is the ensemble member that best describes pressure reduction; see Fig. 6b. The THOM-KF-TR-RTG and LIN-KF-TR-RTG ensemble members qualitatively describe wind speed increase, although they do not reach the maximum value of the best track data for this variable.

Microphysics schemes that depend on cumulus schemes at coarse grid spacing domains are important in regional weather models for providing atmospheric heat and moisture tendencies [41].

Vertical flux of cloud, precipitation, and sedimentation processes of hydrometeors are also included in the microphysics schemes. The microphysics schemes used, as aforementioned, were the Purdue Lin, the WRF single moment 6-class, and the Thompson scheme.

The first includes cloud water, cloud ice, non-precipitable water, rain snow, and graupel [24]. In the second scheme, the ice crystal number concentration is expressed as a function of the ice amount [25]. Finally, the Thompson scheme is a 6-class microphysics scheme with graupel, ice, and rain number. It includes a generalized gamma distribution for each hydrometeor species with snow parameterization, depending on the ice water content and the temperature.

When the Betts-Miller-Janjic cumulus scheme was used in the computer simulations, the mean track errors, when compared to the best track data, were worse than the previous results. In Table 5, the mean track error for each ensemble member with this cumulus scheme is shown.

Table 5 Mean track error (km) taking into account SST forcing 

THOM-BMJ-RTG 85.028
WSM6-BMJ-RTG 68.875
LIN-BMJ-RTG 78.262

The simulated MWSs was underestimated by all the different combinations that defined each ensemble member. The SLP did not exactly reproduce the pressure evolution from the best track information, but this variable described the observed minimum for each ensemble members with SST forcing. The mean SLP for the ensemble members with the Simplified Arakawa-Schubert cumulus scheme is presented in Fig. 8b.

Fig. 8 As in Fig. 5, but for THOM-SAS-RTG, WSM6-SAS-RTG, and LIN-SAS-RTG ensemble members, see Table 1 for acronyms 

The evolution of MWS for each ensemble member is shown in Fig. 8c. The mean track errors for the THOM-SAS-RTG and the LIN-SAS-RTG ensemble members did not improve, showing greater values than those without SST forcing, presented in Table 6.

Table 6 Mean track error (km) taking into account SST forcing 

THOM-SAS-RTG 57.524
WSM6-SAS-RTG 46.091
LIN-SAS-RTG 62.864

5 Conclusion

The simulation of TC Bud using the WRF-ARW model, using the cluster CÓDICE B2, is documented by studying the sensitivity of different cumulus and microphysics schemes. This event was selected because TC Bud reached category-4 hurricane status, spending its lifespan over the Eastern Pacific Ocean and dissipating over the Gulf of California without traversing Mexico. Based on our experimental results, we summarize the following findings:

  1. Overall, the computer simulations have performed well in predicting the track path of the TC.

  2. Nonetheless, in any case, regardless of the cumulus and the microphysics schemes, the simulations underestimated the strength (expressed by either the SLP or the MWS) of the cyclone.

  3. The best ensemble members were those in which the moisture-advection-based trigger scheme, with the Kain-Fritsch cumulus scheme, was used, regardless of the microphysics schemes. This was reflected in the obtained best mean track error values.

  4. The track paths calculated from the THOM-SAS-RTG and the LIN-SAS-RTG ensemble members did not show an improvement compared to those that did not include SST forcing.

Acknowledgments

This work has been funded in part by Instituto Politécnico Nacional under grant SIP-2021083. We thank the Centro de Investigación en Computación of the Instituto Politécnico Nacional, at Mexico City, for providing us with high-performance-computing resources, through the CÓDICE B2 supercomputer. The analysis and visualizations in Fig. 1a were produced with the Giovanni online data system, developed and maintained by the NASA GES DISC.

References

1. Emanuel, K. (2010). Tropical cyclone activity downscaled from NOAA-CIRES reanalysis, 1908–1958. Journal of Advances in Modeling Earth Systems (JAMES), Vol. 2, No. 1, pp. 1–12. DOI: 10.3894/james.2010.2.1. [ Links ]

2. Hamill, T. M., Whitaker, J. S., Fiorino, M., Benjamin, S. G. (2011). Global ensemble predictions of 2009’s tropical cyclones initialized with an ensemble Kalman filter. Monthly Weather Review, Vol. 139, pp. 668–688. DOI: 10.1175/2010mwr3456.1. [ Links ]

3. Emanuel, K. (2005). Increasing destructiveness of tropical cyclones over the past 30 years. Nature, Vol. 436, pp. 686–688. DOI: 10.1038/nature03906. [ Links ]

4. Bouwer, L. M. (2011). Have disaster losses increased due to anthropogenic climate change? Bulletin of the American Meteorological Society, Vol. 92, No. 1, pp. 39–46. DOI: 10.1175/2010bams3092.1. [ Links ]

5. Mendelsohn, R., Emanuel, K., Chonabayashi, S., Bakkensen, L. (2012). The impact of climate change on global tropical cyclone damage. Nature Climate Change, Vol. 2, No. 3, pp. 205–209. DOI: 10.1038/nclimate1357. [ Links ]

6. Domínguez, C., Magaña, V. (2018). The role of tropical cyclones in precipitation over the tropical and subtropical North America. Frontiers in Earth Science, Vol. 6. DOI: 10.3389/feart.2018.00019. [ Links ]

7. Rao, D. V. B., Srinivas, D. (2014). Multi-Physics ensemble prediction of tropical cyclone movement over Bay of Bengal. Natural Hazards, Vol. 70, No. 1, pp. 883–902. DOI: 10.1007/s11069-013-0852-2. [ Links ]

8. Bongirwar, V., Rakesh, V., Kishtawal, C. M., Joshi, P. C. (2011). Impact of satellite observed microwave SST on the simulation of tropical cyclones. Natural Hazards, Vol. 58, No. 3, pp. 929–944. DOI: 10.1007/s11069-010-9699-y. [ Links ]

9. Liu, J., Yang, S., Ma, L., Bao, X., Wang, D., Xu, D. (2013). An initialization scheme for tropical cyclone numerical prediction by enhancing humidity in deep-convection region. Journal of Applied Meteorology and Climatology, Vol. 52, No. 10, pp. 2260–2277. DOI: 10.1175/jamc-d-12-0310.1. [ Links ]

10. Xue, M., Schleif, J., Kong, F., Thomas, K. W., Wang, Y., Zhu, K. (2013). Track and intensity forecasting of Hurricanes: impact of convection-permitting resolution and global ensemble Kalman filter analysis on 2010 Atlantic season forecasts. Weather and Forecasting, Vol. 28, No. 6, pp. 1366–1384. DOI: 10.1175/waf-d-12-00063.1. [ Links ]

11. Anisetty, S. K. A. V. P., Huang, C. Y., Chen, S. Y. (2014). Impact of FORMOSAT-3/COSMIC radio occultation data on the prediction of super cyclone Gonu (2007): a case study. Natural Hazards, Vol. 70, No. 2, pp: 1209–1230. DOI: 10.1007/s11069-013-0870-0. [ Links ]

12. Osuri, K. K., Mohanty, U. C., Routray, A., Kulkarni, M. A., Mohapatra, M. (2012). Customization of WRF-ARW model with physical parameterization schemes for the simulation of tropical cyclones over North Indian Ocean. Natural Hazards, Vol. 63, No. 3, pp. 1337–1359. DOI: 10.1007/s11069-011-9862-0. [ Links ]

13. Chandrasekar, R., Balaji, C. (2012). Sensitivity of tropical cyclone Jal simulations to physics parameterizations. Journal of Earth System Science, Vol. 121, No. 4, pp. 923–946. DOI: 10.1007/s12040-012-0212-8. [ Links ]

14. Fadnavis, S., Deshpande, M., Ghude, S. D., Raj, P. E. (2014). Simulation of severe thunder storm event: a case study over Pune, India. Nat Hazards, Vol. 72, No. 2, pp. 927–943. DOI: 10.1007/s11069-014-1047-1. [ Links ]

15. Srinivas, C. V., Rao, D. V. B., Yesubabu, V., Baskaran, R., Venkatraman, B. (2013). Tropical cyclone predictions over the Bay of Bengal using the high-resolution Advanced Research Weather Research and Forecasting (ARW) model. Quarterly Journal of the Royal Meteorological Society, Vol. 139, No. 676, pp. 1810–1825. DOI: 10.1002/qj.2064. [ Links ]

16. Munsell, E. B., Zhang, F. (2014). Prediction and uncertainty of Hurricane Sandy (2012) explored through a real-time cloud-permitting ensemble analysis and forecast system assimilating airborne Doppler radar observations. Journal of Advances in Modeling Earth Systems (JAMES), Vol. 6, No. 1, pp. 38–58. DOI: 10.1002/2013ms000297. [ Links ]

17. Srivastava, K., Bhardwaj, R. (2014). Analysis and very short range forecast of cyclone ‘‘AILA’’ with radar data assimilation with rapid intermittent cycle using ARPS 3DVAR and cloud analysis techniques. Meteorology and Atmospheric Physics, Vol. 124, No. 1, pp. 97–111. DOI: 10.1007/s00703-014-0307-7 [ Links ]

18. Subramani, D., Chandrasekar, R., Ramanujam, K. S., Balaji, C. (2014). A new ensemble based data assimilation algorithm to improve track prediction of tropical cyclones. Natural Hazards, Vol. 71, No. 1, pp. 659–682. DOI: 10.1007/s11069-013-0942-1. [ Links ]

19. Skamarock, W. C., Klemp, J. B., Dudhia, J., Gill, D. O., Liu, Z., Berner, J., Wang, W., Powers, J. G., Duda, M.G., Barker, D. M., et al. (2019). A Description of the advanced research WRF version 4. NCAR Technical Note NCAR/TN-556+STR. DOI: 10.5065/1dfh-6p97. [ Links ]

20. Powers, J. G., Klemp, J. B., Skamarock, W. C., Davis, C. A., Dudhia, J., Gill, D. O., Coen, J. L., Gochis, D. J., Ahmadov, R., Peckham, S. E., et al. (2017). The weather research and forecasting model: Overview, system efforts, and future directions. Bulletin of the American Meteorological Society, Vol. 98, No. 8, pp. 1717–1737. DOI: 10.1175/BAMS-D-15-00308.1. [ Links ]

21. Beck, J., Brown, J., Dudhia, J., Gill, D., Hertneky, T., Klemp, J., Wang, W., Williams, C., Hu, M., James, E., et al. (2020). An Evaluation of a hybrid, terrain-following vertical coordinate in the WRF-Based RAP and HRRR models. Weather and Forecasting, Vol. 35, No. 3, pp. 1081–1096. DOI: 10.1175/WAF-D-19-0146.1. [ Links ]

22. Gnedin, N. Y., Semenov, V. A., Kravtsov, A. V. (2018). Enforcing the Courant–Friedrichs– Lewy condition in explicitly conservative local time stepping schemes. Journal of Computational Physics, Vol. 359, pp. 93–105. DOI: 10.1016/j.jcp.2018.01.008. [ Links ]

23. Bonekamp, P. N. J., Collier, E., Immerzeel, W. W. (2018). The impact of spatial resolution, land use, and spinup time on resolving spatial precipitation patterns in the Himalayas. Journal of Hydrometeorology, Vol. 19, No. 10, pp. 1565–1581. DOI: 10.1175/JHM-D-17-0212.1. [ Links ]

24. Chen, S. -H., Sun, W. -Y. (2002). A one-dimensional time dependent cloud model. Journal of the Meteorological Society of Japan, Vol. 80, No. 1, pp. 99–118. DOI: 10.2151/jmsj.80.99 [ Links ]

25. Hong, S. -Y., Lim, J. -O. J. (2006). The WRF single-moment 6-class microphysics scheme (WSM6). Journal of the Korean Meteorological Society, Vol. 42, No. 2, pp. 129–151. [ Links ]

26. Thompson, G., Field, P. R., Rasmussen, R. M., Hall, W. D. (2008). Explicit forecasts of winter precipitation using an improved bulk microphysics scheme. Part II: Implementation of a New Snow Parameterization. Monthly Weather Review, Vol. 136, No. 12, pp. 5095–5115. DOI: 10.1175/2008MWR2387.1 [ Links ]

27. Kain, J. S. (2004). The Kain-Fritsch convective parameterization: An update. Journal of Applied Meteorology and Climatology, Vol. 43, No. 1, pp. 170–181. DOI: 10.1175/1520-0450(2004)043<0170:TKCPAU>2.0.CO;2. [ Links ]

28. Janjic, Z. I. (1994). The step-mountain eta coordinate model: Further developments of the convection, viscous sublayer, and turbulence closure schemes. Monthly Weather Review, Vol. 122, No. 5, pp. 927–945. DOI: 10.1175/1520-0493(1994)122<0927:TSMECM>2.0.CO;2. [ Links ]

29. Han, J., Pan, H. -L. (2011). Revision of convection and vertical diffusion schemes in the NCEP global forecast system. Weather Forecasting, Vol. 26, No. 4, pp. 520–533. DOI: 10.1175/WAF-D-10-05038.1. [ Links ]

30. Hong, S. -Y., Noh, Y., Dudhia, J. (2006). A new vertical diffusion package with an explicit treatment of entrainment processes. Monthly Weather Review, Vol. 134, No. 9, pp. 2318–2341. DOI: 10.1175/MWR3199.1. [ Links ]

31. Paulson, C. A. (1970). The mathematical representation of wind speed and temperature profiles in the unstable atmospheric surface layer. Journal of Applied Meteorology and Climatology, Vol. 9, No. 6, pp. 857–861. DOI: 10.1175/1520-0450(1970)009<0857:TMROWS>2.0.CO;2. [ Links ]

32. Dyer, A. J., Hicks, B. B. (1970). Flux-gradient relationships in the constant flux layer. Quarterly Journal of the Royal Meteorological Society, Vol. 96, No. 410, pp. 715–721. DOI: 10.1002/qj.49709641012. [ Links ]

33. Webb, E. K. (1970). Profile relationships: The log-linear range, and extension to strong stability. Quarterly Journal of the Royal Meteorological Society, Vol. 96, No. 407, pp. 67–90. DOI: 10.1002/qj.49709640708. [ Links ]

34. Beljaars, A. C. M. (1994). The parameterization of surface fluxes in large-scale models under free convection. Quarterly Journal of the Royal Meteorological Society, Vol. 121, No. 522, pp. 255–270. DOI: 10.1002/qj.49712152203. [ Links ]

35. Zhang, D., Anthes, R. A. (1982). A high-resolution model of the planetary boundary layer sensitivity tests and comparisons with SESAME-79 data. Journal of Applied Meteorology (1962–1982), Vol. 21, No. 11, pp. 1594–1609. DOI: 10.1175/1520-0450(1982)021<1594:AHRMOT>2.0.CO;2. [ Links ]

36. Tewari, M., Chen, F., Wang, W., Dudhia, J., LeMone, M. A., Mitchell, K., Ek, M., Gayno, G., Wegiel, J., Cuenca, R. H. (2004). Implementation and verification of the unified NOAH land surface model in the WRF model. 20th Conference on weather analysis and forecasting/16th conference on numerical weather prediction, pp. 11–15. [ Links ]

37. Dudhia, J. (1989). Numerical study of convection observed during the Winter Monsoon Experiment using a mesoscale two-dimensional model. Journal of the Atmospheric Sciences, Vol. 46, No. 20, pp. 3077–3107. DOI: 10.1175/1520-0469(1989)046<3077:NSOCOD>2.0.CO;2. [ Links ]

38. Mlawer, E. J., Taubman, S. J., Brown, P. D., Iacono, M. J., Clough, S. A. (1997). Radiative transfer for inhomogeneous atmospheres: RRTM, a validated correlated-k model for the longwave. Journal of Geophysical Research, Vol. 102, No. D14, pp. 16663–16682. DOI: 10.1029/97JD00237. [ Links ]

39. Zeng, X., Beljaars, A. (2005). A prognostic scheme of sea surface skin temperature for modeling and data assimilation. Geophysical Research Letters, Vol. 32, No. 14. DOI: 10.1029/2005GL023030. [ Links ]

40. Emanuel, K., DesAutels, C., Holloway, C., Korty, R. (2004). Environmental control of tropical cyclone intensity. Journal of the Atmospheric Sciences, Vol. 61, No. 7, pp. 843–858. DOI: 10.1175/1520-0469(2004)061<0843:ECOTCI>2.0.CO;2. [ Links ]

41. Nasrollahi, N., AghaKouchak, A., Li, J., Gao, X., Hsu, K., Sorooshian, S. (2012). Assessing the impacts of different WRF precipitation physics in hurricane simulations. Weather and Forecast, Vol. 27, No. 4, pp. 1003–1016. DOI: 10.1175/waf-d-10-05000.1. [ Links ]

Received: November 04, 2021; Accepted: March 01, 2022

* Corresponding author: Ricardo Barrón-Fernández, e-mail: barron2131@gmail.com

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License