3.1. Tornado Outbreak in Nebraska on 5 May 2021
The NOAA National Weather Service (NWS) Storm Prediction Center (SPC) provides daily weather information including storm reports, convective outlooks, and mesoscale analysis at
https://www.spc.noaa.gov/ (accessed on 10 October 2022) According to the SPC storm report for 5 May 2021 (
https://www.spc.noaa.gov/climo/reports/210505_rpts.html, accessed on 10 September 2022) several tornadoes occurred in Nebraska; for example, the towns of Meadow Grove and Clay Center were hit by tornadoes at 2110 UTC and 2130 UTC, respectively. The NWS SPC mesoscale discussion at 2100 UTC noted that convection continues to intensify ahead of a slowly east–southeastward advancing cold front across portions of the Central Plains. A weather watch was not issued for eastern Nebraska, but severe hail and strong winds were not ruled out. The atmospheric changes that happened in the pre-convective environment of the Clay Center tornado, not included in the NWS SPC mesoscale analysis, are the focus of our first case study. Specifically, coincident GOES-16 ABI infrared radiances with a 2 km spatial resolution are combined with the CrIS single field-of-view (SFOV ~14 km) moisture retrievals from one LEO overpass in a spatial fusion approach. Thereafter, these moisture retrievals—here in terms of relative humidity (RH), which is expressed as a percentage [%]—are transferred to a 2 km spatial and a 10 min temporal resolution at the subsequent (or earlier) ABI measurement times. This spatial plus temporal fusion is an attempt to capture the vertical detail from the hyperspectral sounder as though it was from a geostationary perspective. Could this additional information have enabled forecasters to increase the lead time for the tornado watch and possibly warning?
Our fusion investigation starts at the time of the LEO overpass. The SNPP and the NOAA-20, which is hereafter referred to as JPSS-1 (or simply J1), overpasses occurred over Nebraska at approximately 1900 UTC and 1950 UTC, respectively, on 5 May 2021. The ABI radiances are fused with the relative humidity (RH) retrievals derived from the SNPP CrIS overpass at 1900 UTC to produce one sequence of ABI/CrIS-SNPP (or A/C-SN) fused RH data. Then, the ABI is also fused with the JPSS-1 CrIS overpass at 1950 UTC to produce another independent sequence of ABI/CrIS-J1 (or A/C-J1) fused RH data. These two data sets are compared to demonstrate the strengths and weaknesses of the fusion method. Since the SNPP overpass is earlier, it contains a more subtle indication of the moisture gradients that suggest the upcoming convective conditions; this makes it more challenging to transfer the initial information forward to the tornado touchdown time. On the other hand, the JPSS-1 CrIS retrievals of the RH at 500 and 850 hPa at 1950 UTC describe the foreboding convective storms credibly, albeit at a relatively coarse spatial resolution of 14 km, shown in
Figure 2a,e and then as zoomed-in views in
Figure 2c,g. We focus on the lower and upper tropospheric RH (at 850 and 500 hPa) to detect low-level moisture movement into drier regions as well as upward moisture movement in updrafts, which are associated with severe thunderstorm development.
It is noted that the retrieval gaps occur below thick and opaque clouds, since our retrieval algorithm is IR only; nonetheless, retrievals at levels above the clouds are available. The ABI/CrIS spatial fusion results where the RH is now at a 2 km spatial resolution are shown in
Figure 2b,f, with enlargements over the area of interest in
Figure 2d,h (marking the locations of the Clay Center and Meadow Grove). The increased detail in the RH depiction from the ABI/CrIS fusion retrievals is readily apparent, as is the relatively dry mid-level over the moist lower-level air.
To investigate the changes over time, the spatial fusion from the overpass time is enhanced by several temporal fusion steps to create a time series of the RH from 1800 to 2200 UTC in 10 min intervals.
Figure 3a–c shows the GOES-16 ABI brightness temperatures (BTs) for band 15 (12.3 µm), together with the 500 hPa RH fusion results for the SNPP (overpass time is 1900 UTC) and NOAA-20 (overpass time is 1950 UTC). The changes in moisture low and aloft are seen in both fusion time series, but they are more distinct in the ABI/CrIS-J1 series. Note that the temporal fusion goes forward as well as backward in time to construct these images.
To lend confidence to the spatial and temporal fusion process, we investigated the changes in the RH between the two consecutive ABI time steps from the ABI/CrIS-SNPP fusion forward in time (starting at 1900 UTC) and compared them with those from the ABI/CrIS-J1 backward in time (starting at 1950 UTC) over the fifty minutes between their overpasses.
Figure 4 shows the comparison. Some 10 min changes evident in the SNPP going forward are also found in the JPSS-1 going backward and vice versa, indicating that that space–time fusion agrees for the two ABI/CrIS fusion streams. Going forward in time with the ABI/CrIS-SNPP fusion, and backward with the ABI/CrIS-J1 fusion, features not in the initial overpass fusion appear after 50 min, coming into agreement with the opposite fusion loop. While the fusion results forward and backward in time do not entirely match, the gradient changes are reasonably comparable. Differences in the magnitude can be attributed to the different instruments making the measurements as well as the differing initial atmospheric conditions at the times of the two LEO overpasses. This includes an overall weaker moisture signal found in the retrievals at the earlier SNPP overpass time, which were stronger in the later JPSS-1 retrievals (see
Figure 3).
Figure 5 shows the humidity information available from the ABI/CrIS-J1 fusion at 500 and 800 hPa in comparison with the 2-hour (2 h) forecast RAP (Rapid Refresh) [
16] model moisture fields at a 13 km spatial resolution at 2000, 2100, and 2200 UTC. The scale of the moisture features resolved by the RAP model at 850 and 500 hPa is too coarse to resolve the convective development at Meadow Grove or the Clay Center. The area of upper-level moistening and lower-level drying is broadly suggested but not resolved. In addition, the RAP availability at only hourly intervals limits its ability to report more rapid changes. It is of note that these fusion results are solely based on satellite observations and are independent of the model data, whereas the model forecast data (i.e., the RAP) incorporates multiple data sources. Fusion soundings enhance the coverage of soundings in a close proximity to the supercell environments; an increased sounding coverage has been shown to improve the prediction of the location where a tornado is likely to occur [
17].
To further illustrate the type of atmospheric detail fusion can provide, the ABI/CrIS fusion results at the 500 and 850 hPa pressure levels are shown superimposed on the ABI Band 15 (12.3 µm) cloud background every 30 min from 20 UTC to 22 UTC in
Figure 6. Additionally, shown are the corresponding Next Generation Weather Radar (NEXRAD,
https://www.ncei.noaa.gov/products/radar/next-generation-weather-radar, accessed on 10 September 2022) composite reflectivities from nearby radar stations. Signs of convective instability (i.e., dry mid-level air over warm moist low-level air) are already evident at 20 UTC in the ABI/CrIS fusion, confirmed by strong NEXRAD signals, especially within and east of the line of convective clouds (i.e., eastern half of the plots); this situation intensifies quickly, indicating an increasing possibility of severe thunderstorms and tornadoes. A convective boundary layer, possibly capped by a temperature inversion, is also apparent from the 850 hPa panels, where cold dry air moves in from the west towards the Clay Center, whereas the moisture content at the higher 500 hPa level stays relatively large (one example is marked by an oval in the 2100 UTC panel). The NEXRAD shows the location of increased reflectivity but does not indicate the strength of the forcing dynamics.
Next, the ABI/CrIS fusion approach is applied to every tropospheric level to confirm the presence of a convective boundary layer; the resulting fusion humidity profiles for one specific location (about 45 km west of the Clay Center) are shown in
Figure 7. Starting at 1950 UTC (i.e., time of the JPSS-1 or NOAA-20 overpass), the profiles corroborate the existence of low-level cold dry air advection over a moistening boundary layer, which is most pronounced at 2100 UTC, just 30 min before the tornado touched down in Clay Center. At 2030 UTC, the sounding shows evidence of a temperature inversion; here, in terms of the RH, this is seen as a dent at 900 hPa above the moist air near the surface. Half an hour later at 2100 UTC, this feature—and the impact of cold dry air advection together with the further moistening boundary layer—is even more pronounced. This sequence of soundings implies that moisture, heat, and instability in the PBL built up to the point when strong updrafts must have been released, resulting in severe thunderstorms, which happened at 2130 UTC in CC. The hourly RAP 2 h forecast RH profiles valid at 20 and 21 UTC are also shown (in red); however, the RAP profiles show only slight vertical changes from one hour to the next, whereas the ABI/CrIS fusion profiles successfully capture the changes in the vertical moisture distribution from one ABI measurement time to the next.
Next, when zooming in a bit further on the map,
Figure 8 shows the images and profiles for three timesteps starting at 2040 UTC, which are only 10 min apart. In
Figure 8, the rapid atmospheric changes that occurred within a 20 min period are well captured, for example, in the 500 hPa images, the high moisture content moves quickly eastward towards the CC in an area that feeds the convection indicated in the NEXRAD images. The fusion profiles (
Figure 8d), from a location approx. 60 km west of the Clay Center, describe the same environment as mentioned previously (
Figure 7), that is that the atmosphere becomes very unstable with respect to the deep convection caused by the surface air being more moist and warmer compared to the air aloft. The associated radar trends (
Figure 8a) pinpoint the activity within the cloud cover and complement the clear sky ABI/CrIS fusion images.
It may be concluded, from this 5 May 2021 case study, that the fusion of the CrIS high vertical resolution profiles with the high frequency of the ABI radiance measurements allows for a timely detection of the rapid horizontal and vertical changes in the atmosphere, that can cause severe thunderstorms. The NWS SPC convective outlook, issued that day at 1630 UTC, indicated a 10% probability of thunderstorms and less than a 2% probability for tornadoes for all of Nebraska. At 2100 UTC, a weather watch was still not issued, but a locally confined and brief severe threat area east of the front was mentioned (in the NWS SPC mesoscale discussion). The additional ABI/CrIS fusion information may have suggested an adjustment to the probabilities and enabled forecasters to consider issuing a tornado watch. These results would not have been possible with either the ABI data (which lacks the vertical resolution) or the CrIS data (which lacks the temporal and horizontal resolution) alone.
3.2. High Winds and Hail in Texas on 24 May 2022
A strong mesoscale convective system (MCS) produced strong winds and severe hail from the southeast New Mexico/Texas border into parts of west-central Texas on 24 May 2022 (
https://www.spc.noaa.gov/climo/reports/220524_rpts.html, accessed on 10 September 2022). An early projection of the severe weather outlook is shown in
Figure 9a, where the area with an enhanced risk of severe storms (highlighted in orange) was placed too far eastward. The NWS (National Weather Service) Storm Prediction Center (SPC) issues an enhanced severe weather outlook for an area when numerous severe storms are possible and likely more persistent and widespread. The focus of this second case study is whether the additional ABI/CrIS (or geostationary hyperspectral sounder) information would have suggested a better location of the enhanced outlook area. This severe weather event has been evaluated as part of the 2022 NOAA Hazardous Weather Testbed (HWT) experiment, which is a joint project of the NWS (National Weather Service) and the National Severe Storm Laboratory (
https://hwt.nssl.noaa.gov/. accessed on 10 October 2022).
Figure 9a shows the Storm Prediction Center outlook issued at 1630 UTC on 24 May 2022 along with the subsequent wind and hail reports that started coming in around 2100 UTC (there was also one tornado sighting at 2300 UTC near Garden City, TX). Most of the severe weather reports, including those in western Texas near southeastern New Mexico, appear in the ‘slight risk’ severe weather outlook area (yellow), but not in the enhanced outlook area (orange). The GOES-16 ABI images at 2300 UTC in
Figure 9b,c further illustrate the severity, location, and extent of the thunderstorms across western Texas and far southeastern New Mexico at 2100 UTC.
We investigate the potential for GEO/LEO fusion for improving the location and extent of the enhanced outlook area on 24 May 2022. The GOES-16 ABI infrared radiances with a 2 km spatial resolution are fused with the JPSS-1 (i.e., NOAA-20) CrIS SFOV (~14 km) relative humidity (RH) retrievals from the 1950 UTC LEO overpass and then transferred to a 2 km spatial and a 10 min temporal resolution. Additionally, low-level moisture gradients are studied using the information contained in the rotational water vapor lines in the IR window near 12 µm.
The spatially coarse CrIS RH retrievals and the associated 2 km ABI/CrIS fusion results at the time of the overpass, are shown in
Figure 10. Although convective clouds in the northeastern portion of the region reduce the overall retrieval yield, the remaining retrieval data points are sufficient to adequately describe the atmospheric conditions at the overpass time. The CrIS RH retrievals show a strong moisture gradient parallel to the southwestern border of Texas, and the ABI/CrIS fusion further sharpens the gradient depiction and fills in the gaps at both the 500 and 850 hPa pressure levels. Furthermore, the ABI/CrIS fusion captures the changing distribution of the moisture with drying below and moistening aloft (e.g., just east of Fort Stockton, TX).
Next, temporal fusion is applied to the spatial fusion results that were shown in
Figure 10 (right column). Going backwards in time in 10 min interval yields the temporal fusion results, i.e., the changes in 500 and 850 hPa RH, which are shown for every hour from 1600 to 2000 UTC in
Figure 11. Clearly, the tropospheric moisture distribution seen in the fusion RH products at 1600 UTC suggest that the enhanced outlook issued for severe weather at 1630 UTC could have been extended further west, even though the NEXRAD images show little activity. At 2000 UTC, the NEXRAD does confirm the strong convective activity that is coming from the west. Additionally, shown are the fusion RH profiles (in black) and the 13 km RAP 2 h forecast RH profiles (in red) for Fort Stockton at these times (
Figure 11d). Whereas the RAP profiles show little vertical change over the 4 h period, the fusion profiles describe the hourly changes in the water vapor content more realistically. The upward movement of moisture expected in developing convective storms is captured in the fusion profiles. While the JPSS-1 overpass time comes after the 1630 UTC outlook, the GEO/LEO fusion suggests retroactively that the enhanced outlook could have been extended further west. In the future, a geostationary hyperspectral sounder will be able to provide this information in a timely fashion.
Another perspective is offered by the depiction of low-level moisture through the difference of the brightness temperatures in the micro window at 786.25 cm
−1, minus that in the rotational water vapor line at 784.38 cm
−1. This micro split window is sensitive to the moisture gradients in the lowest 100 hPa of the boundary layer [
18].
Figure 12 shows the hourly progression from 1600 to 2300 UTC of low-level (i.e., below 900 hPa) moisture from southeastern New Mexico into the region below the Texas panhandle, where high winds and hail occurred, starting at 2100 UTC. Larger differences (greater than 25 K) are evident in clear skies while smaller differences (less than 15 K) appear in cloudy skies. The strong low-level moisture movement starting at 1600 UTC into the enhanced weather outlook area also suggests that the risk area could have been extended further west at that time.