ProbSevere v3 in west Texas

Severe convection was rather marginal and sparse on 3 May, 2023. However, a few storms of note demonstrated better predictions for ProbSevere v3 (PSv3) over ProbSevere v2 (PSv2). 

A storm in the Texas Panhandle, southwest of Dumas, TX, was in an environment with moderate MLCAPE (1200 J/kg) and low-to-moderate effective bulk shear (30 kt). This storm was around for over an hour before it finally produced severe hail reports (up to 1.5" in diameter). Approximately an hour prior, however, PSv3 was much greater than PSv2; 62% vs. 4% (see Figure 1). In this case, PSv2 was too reliant on the ENI flash rate, which was only 4 fl/min. While lightning is important in PSv3 models, it is less affected by storms with a dearth of lightning than version 2, especially when other parameters show increased severe potential (in this case, MESH and mid-level azshear). 

Figure 1: ProbSevere contours, MRMS MergedRef, and NWS severe thunderstorm warning for a storm near Dumas, TX.

Figure 2: Time series of PSv3 and PSv2 probabilities for the storm near Dumas, TX, which produced 1.5" hail. 

The 00Z sounding from Amarillo, TX (Figure 3) might be able to shed some light on why the MESH was so high (up to 1.6") yet lightning was quite low (~6 fl/min) from both ENI and GLM sensors. In this case, the precipitable water in the column was quite low (about 0.68"). This, combined with the low environmental RH in the mixed phase region of the troposphere probably combined to produce relatively few collisions of liquid and frozen water particles needed to produce charge separation necessary for lightning production. Ice nuclei, on the other hand, seemed to accrete supercooled droplets efficiently, producing large hail stones.

Figure 2: Observed sounding from Amarillo, TX at 00Z on 4 May 2023.

Another storm in southwest Texas, near Pecos, also produced severe hail (1" diameter).  The lightning was also low for this storm, though PSv3 had 51% while PSv2 has 37%. The MRMS azimuthal shear parameters were zero here. This was likely because the average beam height of the lowest tilt of the closest radar was 10-11 km above ground level! The MRMS products compute azimuthal shear between 0-2 km and 3-6 km AGL. The missing azimuthal shear values for this storm reduced the probabilities of severe, but not quite as much in PSv3 as in PSv2. 

Figure 3: ProbSevere contours, MRMS MergedRef, and NWS severe thunderstorm warning for a storm near Pecos, TX.

East Coast Tornadoes

The past few days have seen several strong tornadoes along the U.S. east coast. A shortwave trough with ample upper-level diffluence provided a forcing mechanism for severe storms from Florida to Virginia.

Near Juno Beach, FL, a tornado damaged power lines, homes, buildings, and cars. Maximum wind speeds were estimated at 130 mph (rated EF2). Oddly enough, this tornado was only about 20 miles north of a weaker tornado from the day before.

Figure 1: ProbSevere v3 contours, MRMS MergedReflectivity, and NWS severe weather warnings for a storm near Juno Beach, FL. Outer contours on ProbSevere objects are colored by the probability of tornado.

ProbTor v3 (PTv3) is better calibrated than its v2 counterpart. There was a distinct ramp up in the tornado probability for this storm prior to tornadogenesis, compared to PTv2 (Figure 2). Part of this ramp up was due to higher 0-1 km storm-relative helicity depicted in the HRRR (~160 J/kg), which was much higher than the RAP. Storm rotation was also slowly increasing. Interestingly, this occurred at the same time that lightning and reflectivity-based parameters were decreasing. Despite low overall probability for tornado (20-30%), the ramp up, coupled with the fact that PTv3 remains on the low end overall (max of ~60%) could perhaps have tipped off users to look more closely at this developing storm.

Figure 2: Time series of PTv3 and PTv2 for a tornadic storm on the Florida coast, along with severe reports and NWS severe weather warnings.

The next day, Virginia Beach, VA was hit with an EF3 tornado, with peak winds estimated between 140 and 150 mph. Remarkably, no injuries were reported despite damage to 100 homes. In this case, PTv3 exceeded PTv2, and even hit 60%, which is a very high value for v3. The dip in probability shortly before the tornado was likely due to a pronounced reduction in mid-level azimuthal shear, which quickly rebounded (the 1-3 km mean wind also dropped from 37 kt to 30 kt during that time). 

Figure 3: ProbSevere v3 contours, MRMS MergedReflectivity, and NWS severe weather warnings for a tornadic storm near Virigina Beach, VA. Outer contours on ProbSevere objects are colored by the probability of tornado.

Figure 4: Time series for ProbSevere v3 probabilities, along with reports and NWS severe weather warnings.

Eastern Oregon storm

An elongated trough drew in enough elevated CAPE for a few storms to pop up in eastern Oregon. Very strong lapse rates (≥ 8-9.5 C/km) mixed down strong momentum from roughly 4-5 km AGL, producing downed trees and measured gusts to 65 mph. ProbSevere v3 (PSv3) showed rapid increases in the probability of severe at around 00:20 UTC and 00:50 UTC (Figures 1 and 2). 

Figure 1: ProbSevere v3, MRMS MergedReflectivity, and NWS severe weather warnings for a storm in eastern Oregon.

Figure 2: Time series of PSv3 probabilities for the eastern Oregon storm

The increasing MESH (up to around 1"), composite reflectivity (up to 63 dBZ), and ENI lightning (up to 10 fl/min) contributed to the first rapid rise in the probability of severe. A very strong satellite growth rate at 00:46 led to the next jump, up from 49% to 67%. One can see the cooling cloud tops in the GOES-18 visible/IR sandwich imagery that preceded 3 of the 4 severe wind reports (the first wind report was recorded at 00:41 UTC). 

Figure 3: GOES-18 10.3-µm brightness temperature and 0.64-µm reflectance for the storm in eastern Oregon.

In this example, the probability of severe hail seemed to be the strongest value; however, wind was the only hazard reported. While PSv3 generally has improved discernment among severe weather hazards, this would be a good case to look into to see if we can better incorporate model predictors in PWv3 that take into account the mixing of momentum very high in the atmosphere down to the surface, which is relatively common in the western U.S. 

Elevated hailers

Some elevated storms in Wisconsin have been producing one-inch-diameter hail today. One such storm traveled just west and north of Madison, Wisconsin (Figure 1).

Figure 1: ProbSevere contours and MRMS MergedReflectivity for an elevated storm in Dane Co., Wisconsin.

The storms have been forming in a low-CAPE environment (all of the CAPE is from the most-unstable parcel), with decent effective bulk shear. In ProbSevere version 3, ProbHail has been handling this event better than version 2. While ProbWind v2 is elevated (40%), this isn't really a wind threat, but solely a hail threat. We've seen that ProbSevere v3 consistently discerns the most prominent hazard type better than v2, generally speaking.

One-inch-diameter hail was reported at 15:28 UTC, just a few minutes after a large increase in ProbHail v3 (Figure 2), which was partially due to an increase in the maximum composite reflectivity.

GLM and Lightning Cast

The GLM parallax showed up again Thursday, June 16th, over the PBZ area. This was even more evident than yesterday’s event in WI that was written about in a blog post.  Figure 1 has ProbSevere, LightningCast, GLM Flash Extent Density, and ENTLN data overlaid in a 4-panel.  This case was fairly simple to “self correct” the parallax as the GLM was clearly displaced to the north of ProbSevere (as well as the base reflectivity).  Really once you get a few cases under your belt recognizing the parallax, it’s not too challenging to keep that “self correct” in the back of your mind.  One interesting thing to note about Figure 1 is the storm just outside the PBZ CWA just south of Mount Veron, Ohio (See bottom left in Figure 1).  The ProbSevere and GLM FED are lined up perfectly and this is a great example of utilizing the lower threshold in the colormap. The bullseye shows up much nicer than the larger thresholds in the top two images.

Figure 1: GLM with ProbSevere and LightningCast

- Podium

What’s the best way to look at NUCAPS?

So far this week, I’ve been impressed with the performance of NUCAPS and for the most part Modified NUCAPS with providing thermodynamic profile information in areas that there are overpasses.  However, what would be the best setup to figure out how well NUCAPS is doing.  Let’s dive in:

Today we’re located in Binghamton, NY (BGM) office with a Enhanced risk over the region.  We’ve also got a NUCAPS sounding directly over the region which can help us with evaluating how well things like the SPC Mesoanalysis graphics or other models are doing:

But, visible satellite shows that there are some clouds that could be interfering with the retrievals:

After playing around a bit, it may be a good idea to load the Day Cloud Phase RGB, NUCAPS points, and surface observations into one pane:

This is useful for three reasons:

1. We can see that the red in the Day Cloud Phase are mainly high clouds that are pretty thin and has enough gaps to allow for good retrievals underneath the cirrus.
2. The low clouds are pretty thin except over the southern part of the CWA where some of the retrievals are yellow indicating that caution should be taken when looking at the profiles
3. The surface observations can be used to give an idea of how well the surface T/Td are in the soundings which will impact all of the convective parameter calculations (especially CAPE/CIN)

So what about the soundings?  Looking at the NUCAPS and Modified NUCAPS here:

We can see how things are handled by the soundings.

Here is the original NUCAPS sounding at this point:

Of note, the surface conditions in the NUCAPS sounding is too warm (85 vs 80 at the nearest ob) while the Td were fairly close 69 vs 70.  How did the modified NUCAPS do?

80/70 in the sounding which is a much more likely scenario based on surface observations closest to the time of the overpass/sounding.

Why is this important?  Well, all your convective parameters are based on these two conditions.  MLCAPE drops by around 400 J/kg (~3900 J/kg vs ~3500 J/kg) AND CIN increases from -19 J/kg in the unmodified to -49 J/kg in the modified NUCAPS.  This too lines up with what the visible satellite shows; a lack of boundary layer CU potentially meaning the CAP is holding strong.  

-Pym

ProbSevere V3 Adjusting During New Cell Formation

 ProbSevere V3 has been showing much improvement over version two, with this particular day featuring somewhat pusley mixed-mode storms across upstate New York. This particular group of cells had habitual new updraft development, with storms nearing or briefly becoming severe as high dBZ cores develop aloof then make their way downward. These particular cells were able to form 1.00” hail. The ProbSevere V3 time series reflected the marginally severe nature of these storms very well, with values peaking at 60 to 65%. The very encouraging sign was the peaks and valleys in the ProbSevere V3 time series that showed this group of cells peaking at over 50% severe probabilities, dipping below 50% as the new updraft takes over, and then once again peaking above 50% once the new updraft strengthened. I’d definitely recommend forecasters to take a look at the timeseries to build confidence in cell trends during warning applications.

Shown above is the group of cells, with the initial one showing a core of 60+ dBZ that begins to drift off to the northwest while a newer updraft takes over and moves more east-northeast.

- aerobeaver

ProbSevere Time Series

 I utilized the ProbSevere Time Series in warning operations today (June 16, 2022) within the Pittsburgh, PA CWA.  With only using basic radar products and no other products  (that I’d normally use during warning operations), it was extremely helpful in warning decision making.  After issuing several warnings already and monitoring a storm quickly developing in the western portions of the CWA that was getting close to severe, I issued a warning based on an uptick in the ProbSevere Time Series (See Figure 1 vs Figure 2).  The uptick in ProbHail occurred within about 5-6 minutes or roughly 2 radar scans of the PBZ WSR-88D prompted a warning.  The storm continued for several minutes and regenerated about 30 minutes later that prompted another warning downstream (See Figure 3).

Figure 1: ProbSevere Time Series at 2056z

Figure 2: ProbSevere Time Series at 2100z

Figure 3: ProbSevere Time Series Plots.

- Podium

NUCAPS sounding comparison vs observed 19z PIT sounding.

 With the threat for severe weather today, WFO Pittsburg launched a 19z special sounding to access the atmosphere. The atmosphere featured a moist boundary layer with dewpoints the lower 70s and decent lapse rates throughout much of the atmosphere of around 7 C/km. The low-level wind profile in the observed sounding indicates veering winds with height and low level winds generally from the west southwest to west. The wind profile was helpful in determining the severe weather threat. Based on this low-level wind profile, one can conclude that the tornado threat appears to be fairly low.

There was also a NUCAPS pass across the region around 1741 UTC. I selected a point just to the west of the Pittsburgh office based on the general atmosphere advection to the west based on the observed sounding. Due to this, it would seem that the point just to the west of Pittsburgh at 1741 UTC would be a good approximation of the atmosphere near the Pittsburgh area at the time of the 19 UTC sounding.

The Observed sounding indicates major fluctuations in the dewpoint temperature throughout the atmosphere but the NUCAPs sounding indicates a much smoother moisture profile. Overall, these different methods indicate rather similar PW values to the observed around 1.31 inches and NUCAPS PW of 1.27 inches. Overall, a fairly consistent moisture profile. There are some noticeable differences in the boundary layer moisture profile. The modified NUCAPS is included as the last image. It seems that the modified NUCAPS introduces a bit too much moisture with PW values of 1.43 inches. These higher PW values also lead to higher values of instability as well.

Overall, the NUCAPS sounding was helpful in diagnosing the thermodynamic profile of the atmosphere. It also arrived over an hour earlier than the observed sounding. The sounding indicates a moderately unstable atmosphere across the region capable of leading to severe weather. The moisture profile is smoother in the NUCAPs but it doesn’t seem to affect the overall moisture with PW values being the same. It would be nice to add a model wind profile to the NUCAPs sounding. This way it would be possible to access the tornado threat as well because there is not much available in the thermodynamic profile alone to make this assessment. However, I’m very impressed by the similarity between the NUCAPs and observed sounding.

Observed 19 UTC PIT sounding and NUCAPS from 1741 UTC.

19 UTC PIT Observed Sounding

Location of Observed 19z PIT sounding.

1741 UTC NUCAPS Sounding.

1741 UTC NUCAPS soundings point, brighter green dot.

1741 UTC Modified NUCAPS.

- Marty McFly

Observations from BUF on Thursday June 16th

 Optical Flow Wind Storm Top Divergence Can Aid in Warning Operations

The storm of the day produced a 2.5” hail stone near Cato, NY.

Strong storm top divergence signals an intense or intensifying thunderstorm. Matched with upper level radar scans, satellite interrogation (clean IR), and ProbSevere, optical flow wind products may be another tool to aid the warning forecaster and/or storm scale mesoanalysis,  Below are the corresponding optical flow wind storm top divergence images approaching 4 PM EDT (3 PM CDT) when the sig severe hail was reported. The thunderstorm of interest is centered just south of eastern Lake Ontario.

LightningCast as Graphical Messaging/IDSS Tool

As the DSS and graphics person for BUF, I took the opportunity to highlight the utility of LightningCast for Graphical NowCasts and DSS graphics. The fairly broad brushed nature of the lightning probability contours at the timesteps utilized for graphics I think is a positive for a few reasons: they show actionable probabilities (10+, 25+, 50+, 75+) that users can understand; the contours are akin to annotations on a radar graphic, so would be a time saver; and finally, broad brushed is a way to prompt users and partners to seek out more information, like real time zoomed in radar data, perhaps calling or sending us a chat on NWSChat, or even simply keeping an eye on the sky and listening for thunder. The parallax corrected product in the images below I feel would be more useful for graphics of the nature shown below. Ultimately, for operational use, there would probably need to be some work done to ensure they are not confusing to the user and aesthetically pleasing.

Graphical NowCasts for our websites and social media

Graphic for IDSS Event

A Case for the Sharpened GLM FED Color Scale Used this Week

Over the course of satellite product interrogation this week, the GLM FED stood out as one of the most useful products.  Below are examples from in/near the BUF CWA at 2043z and 2047z. Note that the top of the color scale was set to 128 flashes/5 min over the grid point vs. the default 256 flashes/5 min in AWIPS. The FED also paired well with the MFA and TOE on this 4-pane procedure.

The color curve used this week really popped and correlated well with frequent ENTLN detections. Another example below from shortly after 6 PM EDT (2202 and 2207z)  shows that there will be instances where perhaps an even small color bar range would be useful.

While the strongest storm at this time still showed up well on the AWIPS default color bar range, the 128 flashes/5 min top end of the range (top right) helped it pop even more and corresponded well conceptually with the ENTLN cloud flash detections. Furthermore, the thunderstorms to the north of the strongest thunderstorm showed up better from an SA perspective, if you have been focusing on the strongest storm for warning ops. The smallest color bar range on the bottom left further enhances the above described effect.

Addendum: GLM Flashes and a more smoothed FED

The excellent COD NexLab Satellite and Radar page added within the past year the GLM flash centroids to their GOES derived overlays. From an apples to apples perspective with respect to the ground based lightning detection network displays, this may be a useful product to add into AWIPS.

\

Addendum 2: Would a more smoothed FED be preferable?

Below are a few web based examples of FED displays.

STAR NESDIS Viewer

Severe Storm Tracker

A case can be made that these smoothed FED examples would be somewhat less obtrusive than the default AWIPS FED display, especially for storm scale interrogation by the warning operator or storm scale mesoanalyst.

- Hurricane84