skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: The Effects of Numerical Dissipation on Hurricane Rapid Intensification With Observational Heating
Abstract The computational fluid dynamics of hurricane rapid intensification (RI) is examined through idealized simulations using two codes: a community‐based, finite‐difference/split‐explicit model (WRF) and a spectral‐element/semi‐implicit model (NUMA). The focus of the analysis is on the effects of implicit numerical dissipation (IND) in the energetics of the vortex response to heating, which embodies the fundamental dynamics in the hurricane RI process. The heating considered here is derived from observations: four‐dimensional, fully nonlinear, latent heating/cooling rates calculated from airborne Doppler radar measurements collected in a hurricane undergoing RI. The results continue to show significant IND in WRF relative to NUMA with a reduction in various intensity metrics: (a) time‐integrated, mean kinetic energy values in WRF are ∼20% lower than NUMA and (b) peak, localized wind speeds in WRF are ∼12 m/s lower than NUMA. Values of the eddy diffusivity in WRF need to be reduced by ∼50% from those in NUMA to produce a similar intensity time series. Kinetic energy budgets demonstrate that the pressure contribution is the main factor in the model differences with WRF producing smaller energy input to the vortex by ∼23%, on average. The low‐order spatial discretization of the pressure gradient in WRF is implicated in the IND. In addition, the eddy transport term is found to have a largely positive impact on the vortex intensification with a mean contribution of ∼20%. Overall, these results have important implications for the research and operational forecasting communities that use WRF and WRF‐like numerical models.  more » « less
Award ID(s):
2426563
PAR ID:
10370879
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Advances in Modeling Earth Systems
Volume:
14
Issue:
8
ISSN:
1942-2466
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. This is the Masters thesis of Badrul Hasan from UMBC in the Department of Mechanical Engineering. The co-main advisors are Meilin Yu (UMBC) and Stephen Guimond (UMBC, now with Hampton University). Numerous aspects of human existence, both material and immaterial, can be disrupted by a hurricane. In this work, the computational fluid dynamics of hurricane rapid intensification (RI) are studied by running idealized simulations with two different codes: a community-based, finite-difference/split-explicit model (WRF) and a spectral-element/semi-implicit model (NUMA). Rapid intensification is what RI stands for, how a hurricane gets stronger quickly. The main goal of this study is to find out how implicit numerical dissipation (IND) affects the energy of the vortex's response to heating, which describes the fundamental dynamics of the hurricane RI process. The heating that is taken into account here is derived from data. These observations include four-dimensional, fully nonlinear, latent heating/cooling rates estimated using airborne Doppler radar readings acquired during RI in a hurricane. The results show that WRF has more IND than NUMA, with a decrease in several intensity parameters, such as (1) the time-integrated mean kinetic energy values that WRF predicts are 20% lower than those that NUMA predicts and (2) the peak, localized wind speeds that WRF predicts are 12 meters per second slower than those that NUMA predicts. To make a time series of intensity similar to NUMA, the eddy diffusivity values in WRF need to be less than those in NUMA by about 50%. Various analyses are conducted comparing WRF with NUMA. Kinetic energy budgets reveal that the pressure contribution is the primary cause in the model variations, with WRF generating an average 23% lower vortex energy input. The IND is associated with the low-order spatial discretization of the pressure gradient in WRF. In addition, the mean contribution of the eddy transport term to the vortex intensification is determined to be 20% positive. These findings have significant implications for the academic and operational forecasting communities that employ WRF and similar numerical methods. 
    more » « less
  2. Recent studies have shown that climate change and global warming considerably increase the risks of hurricane winds, floods, and storm surges in coastal communities. Turbulent processes in Hurricane Boundary Layers (HBLs) play a major role in hurricane dynamics and intensification. Most of the existing turbulence parameterizations in the current numerical weather prediction (NWP) models rely on the Planetary Boundary Layer (PBL) schemes. Previous studies (Zhang 2010; Momen et al. 2021) showed that there is a significant distinction between turbulence characteristics in HBLs and regular atmospheric boundary layers (ABLs) due to the strong rotational effects of hurricane flows. Nevertheless, such differences are not considered in the current schemes of NWPs, and they are primarily designed and tested for regular ABLs. In this talk, we aim to bridge this knowledge gap by conducting new hurricane simulations using the Weather Research and Forecasting (WRF) model as well as large-eddy simulations. We investigate the role of the PBL parameterizations and momentum roughness length in multiple hurricanes by probing the parameter space of the problem. Our simulations have shown that the most widely used WRF PBL schemes do not capture the hurricane intensification properly and underestimate their intensity. We will present that decreasing the roughness length close to the values of observational estimates and theoretical hurricane intensity models in high wind regimes (≳ 45 m s-1) led to significant improvements in the intensity forecasts of strong hurricanes. Furthermore, by decreasing the existing vertical diffusion values, on average more than 20% improvements in hurricane intensity forecasts were obtained compared to the default runs. Our results provide new insights into the role of turbulence parameterizations in hurricane dynamics and can be employed to improve the accuracy of real hurricane forecasts. The implications of these results and improvements for coastal resiliency and fluid-structure interactions will also be discussed. 
    more » « less
  3. Hurricanes are among the costliest natural disasters in the United States and regularly inflict severe damage on urban infrastructure. Accurate forecasts are therefore essential for preparedness and limiting these extreme events' economic toll. Numerical weather‑prediction (NWP) models—such as the Weather Research and Forecasting (WRF) system—are powerful forecasting tools. However, some of their physical parameterizations were neither designed for nor tested with real hurricanes. This thesis addresses that gap by evaluating two key parameterizations in WRF: (i) subgrid‑scale (SGS) turbulence schemes and (ii) surface‑roughness and urban canopy treatments. The first part of the study investigates how SGS eddy‑viscosity choices affect hurricane intensity, turbulence, and wind profiles. Large‑eddy simulations (LES) of five major hurricanes were run with a 1.5‑order, three‑dimensional turbulent‑kinetic‑energy (TKE) SGS scheme. Each storm was simulated under three eddy‑viscosity settings— default, halved, and doubled—yielding 15 cases. A parallel set of 10 cases employed an alternative nonlinear backscatter and anisotropy (NBA) SGS scheme. Two idealized LES runs and one fine-grid (~80 m) nested simulation brought the total to 33. Reducing SGS stress intensified storms by raising boundary‑layer wind speeds and lowering the altitude of peak winds, improving surface‑wind forecasts by ~9 % and minimum sea‑level pressure by ~29 % relative to the default setting. These results reveal that standard SGS models are overly dissipative because they overlook the rotational suppression of turbulence, underscoring the need for SGS schemes tailored to hurricane dynamics. The second part assesses how aerodynamic roughness length (z0) and urban‑canopy schemes shape near‑surface winds over cities. For four land‑falling hurricanes affecting Houston and New Orleans, increasing z0 in the Single‑Layer Urban Canopy Model (SLUCM) reduced modeled wind speeds and cut mean absolute error (MAE) by ~20 %, whereas decreasing z0 introduced large positive biases. Additional experiments compared three urban options—Bulk (no‑urban), SLUCM, and the multi‑layer Building Energy Model (BEM). The Bulk scheme delivered the most accurate surface‑wind forecasts in every nested domain, while SLUCM slightly outperformed BEM in the limited vertical‑profile data. These findings highlight the need to recalibrate urban schemes and surface‑drag parameters when applying WRF to hurricane‑force winds. 
    more » « less
  4. The FV3GFS is the current operational Global Forecast System (GFS) at the National Centers for Environmental Prediction (NCEP), which combines a finite-volume cubed sphere dynamical core (FV3) and GFS physics. In this study, FV3GFS is used to gain understanding of rapid intensification (RI) of tropical cyclones (TCs) in shear. The analysis demonstrates the importance of TC structure in a complex system like Hurricane Michael, which intensified to a category 5 hurricane over the Gulf of Mexico despite over 20 kt (10 m s−1) of vertical wind shear. Michael’s RI is examined using a global-nest FV3GFS ensemble with the nest at 3-km resolution. The ensemble shows a range of peak intensities from 77 to 159 kt (40–82 m s−1). Precipitation symmetry, vortex tilt, moisture, and other aspects of Michael’s evolution are compared through composites of stronger and weaker members. The 850–200-hPa vertical shear is 22 kt (11 m s−1) in the mean of both strong and weak members during the early stage. Tilt and moisture are two distinguishing factors between strong and weak members. The relationship between vortex tilt and humidification is complex, and other studies have shown both are important for sheared intensification. Here, it is shown that tilt reduction leads to upshear humidification and is thus a driving factor for intensification. A stronger initial vortex and early evolution of the vortex also appear to be the key to members that are able to resist the sheared environment. 
    more » « less
  5. Tropical cyclones are one of the deadliest natural disasters in the world that cause significant damage to the environment and infrastructure. The Hurricane Boundary Layer (HBL) plays a major role in hurricane dynamics and its intensification. Most of the existing vertical diffusion parameterizations in the current numerical weather prediction models rely on the Planetary Boundary Layer (PBL) schemes. Previous studies (Momen et al. 2021; Romdhani et al. 2022) showed that there is a significant distinction between turbulence characteristics in HBLs and regular atmospheric boundary layers (ABLs) due to the strong rotational effects of hurricane flows. Nevertheless, such differences are not considered in the current PBL schemes, and they are primarily designed and tested for regular ABLs. In this talk, we aim to bridge this knowledge gap by conducting real hurricane simulations using the Weather Research and Forecasting (WRF) model. We investigate the role of the PBL height and eddy momentum exchange coefficients in five intensifying hurricanes by probing the parameter space of the problem. Our simulations have shown that the most widely used WRF PBL schemes do not capture the hurricane intensification properly and underestimate their intensity. We will demonstrate how limiting the amount of the vertical transport of momentum greatly benefits the skill of forecasting in major hurricane simulations. We will also present how changing the height of the PBL significantly impacts the accuracy of the forecasts. By reducing the PBL height, simulated hurricanes become stronger and larger – representing the actual rapid intensification process much more accurately. Not only changes are seen in the predicted wind intensities, but also remarkable impacts are observed in storm size, the radius of maximum wind speed, hurricane track, and minimum sea level pressure. The results of this study provide insights into the role of vertical diffusion parameterizations in hurricane dynamics. Our findings can be used to improve the accuracy of real hurricane forecasts in numerical weather prediction models. 
    more » « less