Advanced Research WRF High Resolution Simulations of the
Inner Core Structure of Hurricanes Katrina, Rita and Wilma (2005)


The Hurricane Group of the Mesoscale and Microscale Meteorology (MMM) division of the National Center for Atmospheric Research (NCAR) produced real-time prediction and retrospective research simulations of Atlantic and eastern North Pacific tropical cyclones (TCs) using the Advanced Research core of the Weather Research and Forecasting model (WRF) from 2004 to 2007, inclusive.
(See Davis et al. 2008 for a full description of the model setup and parameters.)


Verification of the intensity and track errors of the WRF forecasts has shown that the 4 km simulations perform as well as, and occasionally superior to, the National Hurricane Center (NHC) official forecasts and other operational forecast systems. 

      
2005 intensity (kt) and position (n mi) errors for the Advanced Hurricane WRF (AHW) forecasts run with an inner moving nest of 4 km (AHW4).  Results  from other forecast techniques are defined as the NHC Official (OFCL), the Geophysical Fluid Dynamics Laboratory hurricane model (GFDL), the Florida State Superensemble (FSSE), the five day NHC statistical model (SHF5) and decay techniques (DSHP), the Met Office (UKMO), the NCEP Aviation Model (AVNO), the Navy Operational Global Atmospheric Prediction System (NGPS), the statistical climatology and persistence model (CLP5), and no change (NCHG). (Adapted from Figure 1 of Davis et al. 2008)

Verification of the inner core structure of the simulated storms, examined with a finest mesh of 1.33 km, has proven less successful.  Realistic eye, eyewall and inner spiral rainband structures are simulated on the finest mesh (1.33 km; middle image below), but the eye is generally too small and the eyewall convection too wide.  As an example,
compare the 4 km (left) and 1.33 km (middle) Hurricane Katrina (2005) reflectivity patterns at 0000 UTC 29 August (48h forecast) to the NOAA/AOML/HRD P-3 lower fuselage radar image (~1.5 km resolution) at 2317 UTC 28 August (right).

    
The scales of the images are nearly identical with both WRF images being 320 km on a side and the P-3 radar being 360 km.

An unusually large number of polygonal eyewall shapes were noted on both the 4 km and 1.33 km grids, as seen in the 1.33 km Katrina reflectivity plot above (above middle) and 4 km Hurricane Rita (2005) reflectivity (dBZ; left) and precipitable water (cm; right) plots shown below
at 0300 UTC 24 September (39h forecast).   
*** Additional times from both of these runs are shown in the "Loops" section below. ***

    

The vertices of the polygonal eyewalls (three in the case of the Rita triangular eyewall above) show up equally well in the fields of radar reflectivity and precipitable water (the column integrated water vapor) and thus a combination of both are shown in the animations below highlighting the ubiquitous nature of elliptical, triangular, square, and higher wavenumber polygons in the eyewalls of WRF simulated TCs across different resolutions and incarnations of the model.

Reflectivity & Precipitable Water Loops
~ Hurricane Rita, real time AHW 2.1, 00 UTC 23 September - 12 UTC 24 September (forecast hours 12-48), 4 km precipitable water
            * A triangular eyewall is evident throughout most of the run, becoming more distinct as Rita approaches the Texas coast  

~ Hurricane Rita, retrospective AHW 2.1, 18 UTC 21 September - 12 UTC 23 September (forecast hours 18-60), 1.33 km precipitable water
            * Triangles, squares and pentagons abound in this higher spatial resolution animation

~
Hurricane Wilma, real time AHW 2.1, 00 UTC 20 October - 00 UTC 24 October (forecast hours 01-96), 4 km Reflectivity
            * Note the triangular eyewall from 21 UTC 20 October - 04 UTC 21 October (frames 21-28) and again after 00 UTC 23 October (frames 72+)

~ Hurricane Katrina, retrospective AHW 2.1, 00 UTC 28 August - 11 UTC 29 August (forecast hours 24-59), 1.33 km precipitable water
            * Note the square eyewall from 1700-2000 UTC 28 August (frames 18-21) that transitions to a triangle at 2100 UTC

~ Hurricane Katrina, retrospective AHW 2.2, 00 UTC 28 August - 12 UTC 29 August (forecast hours 24-60), 1.33 km precipitable water
            * This newer version of the model has fewer squares, more rounded quadrilaterals, and a spectacular triangular eyewall on 29 August (frames 26-29)    




Polygonal eyewalls have been documented since the first radar images of 
TCs were captured (e.g. Lewis and Hawkins 1982), but only recently has 
their existence been explained using barotropic instability in a series of 
papers by Wayne Schubert and Jam
es Kossin (Schubert et al. 1999; 
Kossin and Schubert 2001, 2004; Kossin et al. 2002, 2004; ), although the 
idea was first proposed by M
uramatsu (1986).  Briefly, mature and rapidly 
intensifying TCs exhibit elevated rings of potential vorticity (or relative 
vorticity; see image at right) on the inner edge of the eyewall (K
ossin and 
Eastin 2001; Mallen et al. 2005).  Analogous to Rossby waves in the large 
scale circulation of the Earth, counter-propagating waves
with respect to 
 the flow
(deemed vortex Rossby waves or VRWs; Montgomery and 
Kallenbach 1997) exist on the oppositely signed radial gradients of 
potential vorticity (PV). 




               

   If these waves become phase-locked, they
grow in concert and lead to the exponential
   instability of the ring, whereby the eyewall vorticity pools into discrete areas, creating
   mesovortices (see Figure 3 of Schubert et al. 1999).  Depending on the initial conditions of
   the PV ring, the mesovortices either merge over time and relax to a monopole (Schubert et
   al. 1999; Chen and Yau 2003), or remain separate to form a quasi-steady,
   cyclonically-rotating lattice of vortices that give the appearance of elliptical (two
   mesovortices), triangular (three mesovortices) or polygonal (four or more mesovortices)
   eyewalls (Kossin and Schubert 2001).

   While most of the polygonal shaped eyewalls are short lived (< 2-3 hr) and change readily
  
between shapes and intervals of roundness in both nature and the WRF model, TCs (real
   and simulated) will sometimes "lock on" to shapes
for several hours such as the starfish in
   Hurricane Isabel (2003) (left; Kossin and Schubert 2004) and the triangle or squares
   and triangles in the WRF simulated Rita and Katrina, respectively.  Higher temporal
   resolution model output (every 10 min) was generated for these instances so the direction
   and rate of rotation of the vertices could be determined and compared to theory and
   previous observational studies of polygonal eyewalls.




High temporal Resolution (10 min) Precipitable Water Loops
~
Hurricane Katrina, retrospective AHW 2.1, 12-18 UTC 28 August (forecast hours 36-42, every 10 min), 1.33 km precipitable water

          * Squares and pentagons abound in this higher time resolution loop

~ Hurricane Katrina, retrospective AHW 2.2, 01-03 UTC 29 August (forecast hours 49-51, every 10 min), 1.33 km precipitable water
          * A long lived triangular eyewall appears to rotate clockwise in this newer version of the model


  The very active 2005 Atlantic TC season was fortunate enough to feature the Hurricane Rainband and Intensity Change
  Experiment (RAINEX) which was designed to explore the dynamic relationships between the eyewall, inner spiral rainbands
  and intensity change in tropical cyclones.  Extensive observations were collected in both Katrina and Rita and copious
  amounts of data are available to validate our WRF simulations against and to investigate whether the storms actually
  featured the realistic looking polygonal eyewalls simulated with the model.


Below is a series of reflectivity images from the NOAA P-3 lower fuselage radar in Katrina at 1753, 1931, 2233, and 2347 UTC 28 August as the plane was passing through the eye.  While the observations do show a few bumps and undulations along the inner edge of the eyewall reflectivity gradient (especially evident at 1753 UTC; far left), there are no indications of any polygonal eyewalls, in stark contrast to all of the WRF runs of Katrina above that show distinct squares or triangles throughout the period of study.

       

There are also no triangles or squares, but perhaps an ellipse at 2153 UTC (far right), in the reflectivity images below of Hurricane Rita at 1658, 1744 and 2153 UTC 23 September

     


 

  Hurricane Wilma also shows no indications of a triangular eyewall in this LOOP of lower fuselage
  radar images from 1842-2309 UTC 20 October.  This is in stark contrast to plot to the left which
  shows the 23 h forecast of precipitable water from the real time, 4 km WRF run of Wilma initialized
  at 00 UTC 20 October (the same run shown in the Wilma reflectivity loop above). 

  The first two images of the loop do, however, show some straight line edges to the eastern half of
  the eyewall and there is a hint of a pentagon- or hexagon-shaped eyewall at 2019 UTC.




Doppler wind analyses of the inner ~100 km are also available for Katrina and Rita at .5 and 3 km altitudes and can be compared to the model predicted wind speeds.  Two individual flights across the storm center have been overlaid in the figures below to give more complete coverage of the inner 100 km to compare to the WRF output.  The color scale of the wind speed has been matched as closely as possible between the HRD generated plots and the WRF output plotted with GEMPAK.