The Case of the 100-year-old Ash Cloud

Lost in all the commotion caused by Hurricane Sandy, a curious event occurred on the other side of the country on 30 October 2012. A cloud of ash obscured the skies of Kodiak Island, Alaska, diverting flights in the region and forcing the people of Kodiak to stay inside or wear masks. Alaska has quite a few volcanoes, so this may not be a big thing to them except, this was no ordinary volcanic eruption: it was the leftovers of a volcanic eruption from 100 years ago!

The volcano that came to be known as Novarupta erupted on 6 June 1912. It was one of the largest volcanic eruptions of recorded history. It was 10 times more powerful than Mt. St. Helens with 100 times more ash. The explosion was heard more than 1100 km (700 miles) away in Juneau. The force of the eruption caused nearby Mt. Katmai to collapse on itself (10 km away). It formed the Valley of Ten Thousand Smokes and, most importantly for us, covered the surrounding land with 150 m (500 ft) of ash.

This pile of ash – still there today – can be lifted by a stiff breeze (or, more appropriately, “strong breeze” or higher on the Beaufort wind scale), and blown pretty high off the ground (4000 ft according to the news report). This isn’t the first time this has happened. MODIS observed the same thing back in 2003.

So, what did VIIRS see? Here’s the “true color” image, the RGB composite of channels M-03 (0.488 µm, blue), M-04 (0.555 µm, green) and M-05 (0.672 µm, red):

VIIRS "true color" RGB composite of channels M-03, M-04 and M-05, taken 22:23 UTC 30 October 2012

VIIRS "true color" RGB composite of channels M-03, M-04 and M-05, taken 22:23 UTC 30 October 2012

Be sure (as with all the images) to click on the image, then on the link below the banner to see it at full resolution. (The link contains the dimensions of the full size image.)

The ash cloud (blowing right over the center of Kodiak Island) is not as obvious in this image as it was in the MODIS image in the link above, although it is visible. To be fair, the plume was much more optically thick in 2003, and there were fewer clouds and less snow to confuse it with.

Here is the false color (“pseudo-true color” or “natural color”) image, the RGB composite of channels M-05 (0.672 µm, blue), M-07 (0.865 µm, green) and M-10 (1.61 µm, red):

VIIRS false color RGB composite of channels M05, M-07 and M-10, taken 22:23 UTC 30 October 2012

VIIRS false color RGB composite of channels M05, M-07 and M-10, taken 22:23 UTC 30 October 2012

Hmmm. Once again, the ash plume is visible but not particularly noticeable. Is there a way to highlight the ash plume to make it easier to see?

EUMETSAT (the European Organisation for the Exploitation of Meteorological Satellites) has defined an RGB composite for detecting dust. Their product, which was developed primarily to detect dust storms over the Saharan desert, uses channels that are present (or similar to ones that are present) on VIIRS. This means we can apply the dust product for VIIRS as the difference between M-16 and M-15 (red), the difference between M-15 and M-14 (green) and M-15 by itself (blue), all in units of brightness temperature. If you do that, and use the same color scaling they use, you get this image:

The EUMETSAT Dust RGB composite applied to VIIRS for 22:23 UTC 30 October 2012

The EUMETSAT Dust RGB composite applied to VIIRS for 22:23 UTC 30 October 2012

The arrow points to the source region of the ash plume. In this RGB composite, dust shows up as hot pink (magenta), but it’s barely visible here. The reason is that this dust product is primarily useful where there is a large temperature contrast between the dust plume and the background surface, which we don’t have here.

A more common way to detect volcanic ash is to use the “split-window difference”. The “split-window difference” is the difference in brightness temperature between a 10.7-11.0 µm channel and a 12.0 µm channel. This difference is useful because volcanic ash has a difference of opposite sign to most everything else. Here’s what the split window difference (M-15 – M-16) looks like for this case:

VIIRS "Split-window difference" image from 22:23 UTC 30 October 2012

VIIRS "Split-window difference" image from 22:23 UTC 30 October 2012

This image has been scaled so that the colors range from -1 K (black) to +7 K (white). The ash plume stands out a bit more here by being much darker than the background. The only problem is, it isn’t perfect. Large amounts of water vapor, optically thick clouds, desert surfaces and boundary layer temperature inversions can all produce a negative difference (just like volcanic ash does).

These problems can be overcome to a certain extent by combining the “split-window difference” with a Principal Component Image (PCI) analysis technique. (This technique is too complicated to describe here but, if you have access to AMS journals, check out these journal papers.) Now, the ash plume is the only thing that’s black:

VIIRS PCI analysis image from 22:23 UTC 30 October 2012

VIIRS PCI split window analysis image from 22:23 UTC 30 October 2012. Image courtesy Don Hillger. Upside-down text courtesy McIDAS-X.

Notice the smaller plume identified by the orange arrow. This plume is not easy to identify in any of the previous images. The PCI technique works well. But, we’re not going to stop there.

Remember the dust plumes off the Cape Verde islands? They produced a strong signal in the difference between M-12 (3.7 µm) and M-15 (10.7 µm) due to solar reflection. Does a 100-year-old ash plume produce a similarly strong signal? See for yourself:

VIIRS channel difference image between M-12 and M-15 from 22:23 UTC 30 October 2012

VIIRS channel difference image between M-12 and M-15 from 22:23 UTC 30 October 2012

It does produce a signal, but it’s not as bright as the surrounding clouds. The color scale here ranges from -2 K (black) to +90 K (white).

M-06 (0.746 µm) is highly sensitive to anything that reflects solar radiation in the atmosphere or on the surface, which we learned from Hurricane Isaac. Here’s what the M-06 image looks like:

VIIRS channel M-06 image, taken 22:23 UTC 30 October 2012

VIIRS channel M-06 image, taken 22:23 UTC 30 October 2012

“Big deal,” you say. “None of those are better than the PCI analysis.” That may be true, but watch what happens when we combine M-06, the M-12 – M-15 image and the split-window difference image in a single RGB composite:

VIIRS RGB composite of M06 (blue), M12 - M15 (green) and M15 - M16 (red), taken 22:23 UTC 30 October 2012

VIIRS RGB composite of M06 (blue), M12 - M15 (green) and M15 - M16 (red), taken 22:23 UTC 30 October 2012

In this composite, blue values represent the M-06 reflectance scaled from 0 to 1.6, green values represent the brightness temperature difference between M-12 and M-15 scaled from -2 K to +90 K, and red values represent the brightness temperature difference between M-15 and M-16 scaled from -1 K to +7 K.

From a theoretical perspective, this RGB composite does exactly what you want: make the thing you’re trying to detect the only thing that is a certain color. For example, the ash plumes are the only things in this image that are green. From a practical perspective, however, this RGB composite doesn’t work so well. It only works because the ash plume is over water (otherwise M-06 wouldn’t be very useful). It only works during the day, where M-06 is available and the difference between M-12 and M-15 is significant (no solar component to M-12 at night).

Plus, the rainbow of colors is difficult to make sense of: green ash; clouds ranging from light blue to purple to orange (a function of optical thickness, particle size, and phase); bright purple snow; dark purple vegetation; maroon water. It’s not exactly pleasing to the eye. In contrast, the PCI analysis technique that uses the split-window difference works day and night, over ocean and over land. And it isn’t confusing to look at. Maybe we should have stopped when we got to the PCI technique. But then, we wouldn’t have learned anything new.

Daniel, Emilia and Fabio, oh my!

It’s been a while since we last looked at some tropical cyclones with VIIRS. If you don’t keep up to date on tropical activity, you might not know there that have been a few. Granted, since Debby dumped a bunch of rain on Florida three weeks ago, the Atlantic basin has been pretty quiet. The East Pacific basin, however, has had one storm after another. The national media has largely ignored them since they have posed no threat to any landmasses. See this article from the L.A. times. Boring! Unless you can capture video of Jim Cantore struggling to stand upright, it isn’t a hurricane, right?

Wrong! First of all, eastern Pacific hurricanes affect some major shipping lanes. Second, and this is true of all hurricanes: they transport energy and moisture and help moderate the temperature imbalance between the tropics and mid-latitudes. They are important components of global energy transport.

In this post, we are going to compare the view of hurricanes provided by VIIRS against the view provided by GOES (specifically GOES-15). On 9 July 2012, there were two storms in the East Pacific: Daniel and Emilia.

Here is the GOES-15 view of Daniel followed by the VIIRS view of Daniel in their respective visible channels:

GOES-15 visible image (channel 1) of Hurricane Daniel, taken 22:45 UTC 9 July 2012

GOES-15 visible image (channel 1) of Hurricane Daniel, taken 22:45 UTC 9 July 2012. Image courtesy John Knaff.

VIIRS visible image (channel I-01) of Hurricane Daniel, taken 22:29 UTC 9 July 2012

VIIRS visible image (channel I-01) of Hurricane Daniel, taken 22:29 UTC 9 July 2012

Both images have the same latitude and longitude lines printed on them for reference and they both use the same color scales. If you zoom in, you’ll notice that the VIIRS image, with ~375 m resolution at nadir shows a bit more detail than the 1 km (1000 m) resolution GOES image. The additional detail provided by VIIRS really stands out in the infrared (IR) window channels, where GOES has 4 km resolution and VIIRS still has ~375 m resolution:

GOES-15 IR image (channel 4) of Hurricane Daniel, taken 22:30 UTC 9 July 2012

GOES-15 IR image (channel 4) of Hurricane Daniel, taken 22:30 UTC 9 July 2012

VIIRS IR image (channel I-05) of Hurricane Daniel, taken 22:29 UTC 9 July 2012

VIIRS IR image (channel I-05) of Hurricane Daniel, taken 22:29 UTC 9 July 2012

Now, it is worth noting that the high resolution IR image of VIIRS shown above comes from channel I-05, which is centered at 11.45 µm. The GOES image was produced from Imager channel 4, which is centered at 10.7 µm, so the two channels don’t exactly have the same spectral properties. VIIRS has a 10.7 µm IR channel as one of its moderate resolution bands (M-15). Here’s what that image looks like:

VIIRS IR image (channel M-15) of Hurricane Daniel, taken 22:29 UTC 9 July 2012

VIIRS IR image (channel M-15) of Hurricane Daniel, taken 22:29 UTC 9 July 2012

There isn’t a big difference between the two VIIRS channels, although you can see a bit more detail in the higher resolution (I-05) image.

On the previous orbit, VIIRS caught images of Hurricane Emilia, which was also in the view of GOES-15. Here’s how the images compare:

GOES-15 visible image (channel 1) of Hurricane Emilia, taken 21:00 UTC 9 July 2012

GOES-15 visible image (channel 1) of Hurricane Emilia, taken 21:00 UTC 9 July 2012. Image courtesy John Knaff.

VIIRS visible image (channel I-01) of Hurricane Emilia, taken 20:48 UTC 9 July 2012

VIIRS visible image (channel I-01) of Hurricane Emilia, taken 20:48 UTC 9 July 2012

GOES-15 IR image (channel 4) of Hurricane Emilia, taken 20:48 UTC 9 July 2012

GOES-15 IR image (channel 4) of Hurricane Emilia, taken 20:48 UTC 9 July 2012

VIIRS IR image (channel I-05) of Hurricane Emilia, taken 20:48 UTC 9 July 2012

VIIRS IR image (channel I-05) of Hurricane Emilia, taken 20:48 UTC 9 July 2012

In addition to the resolution differences, there is also a time difference of ~15 minutes between the VIIRS images and the GOES images. If you were to overlap these images, you would see that Emilia rotated a bit during that time. Emilia was not willing to hold the same pose for that long when having her picture taken. Once again, the M-15 image from VIIRS looks pretty similar to the I-05 image, so there’s no pressing need to show it.

Finally, let’s compare GOES-15 with VIIRS on Hurricane Fabio, which formed about a week after Daniel and Emilia were hurricanes.

GOES visible image (channel 1) of Hurricane Fabio, taken 20:30 UTC 15 July 2012

GOES-15 visible image (channel 1) of Hurricane Fabio, taken 20:30 UTC 15 July 2012. Image courtesy John Knaff.

VIIRS visible image (channel I-01) of Hurricane Fabio, taken 20:36 UTC 15 July 2012

VIIRS visible image (channel I-01) of Hurricane Fabio, taken 20:36 UTC 15 July 2012

GOES-15 IR image (channel 4) of Hurricane Fabio, taken 20:30 UTC 15 July 2012

GOES-15 IR image (channel 4) of Hurricane Fabio, taken 20:30 UTC 15 July 2012

VIIRS IR image (channel I-05) of Hurricane Fabio, taken 20:36 UTC 15 July 2012

VIIRS IR image (channel I-05) of Hurricane Fabio, taken 20:36 UTC 15 July 2012

The GOES and VIIRS images of Fabio were taken only 6 minutes apart, so there is less movement to impede the comparison.

In all three hurricanes, you can see a lot more structure to the VIIRS images in the both the visible and IR channels. It’s as if GOES represents a standard definition TV camera, and VIIRS represents a hi-def TV camera. All those wrinkles GOES is smoothing over are showing up in VIIRS. Daniel, Emilia and Fabio are going to need more makeup. (Or, they would if they weren’t already dead.)

Cape Verde Waves and Plumes

Cape Verde is an island nation off the west coast of Africa, located in the North Atlantic. The islands are a popular initiation point for tropical storms. The original capital of the 10-island archipelago was sacked twice by Sir Francis Drake, the same one who, in his later years, would fail to sack the villages along Lake Maracaibo in Venezuela due to Catatumbo lightning. That guy really got around, and I mean that literally: he circumnavigated the globe between 1577 and 1580, sacking nearly every village and boat he came across. But, this isn’t about Francis Drake – it’s about the Cape Verde islands and the amazing view of them captured by VIIRS.

False color RGB composite of VIIRS channels I-1, I-2 and I-3 taken 14:41 UTC 6 June 2012

False color RGB composite of VIIRS channels I-1, I-2 and I-3 taken 14:41 UTC 5 June 2012

Can you see the 10 major islands? One of them (Santa Luzia) is almost obscured by clouds. If you click on the image, you’ll see each of the major islands identified. Go ahead and click on it. It will help for later.

The image above was made from the RGB composite of VIIRS high-resolution imagery channels I-01, I-02 and I-03. While it technically is a false color image (uses reflectance at 0.64 µm [blue],  0.865 µm [green] and 1.61 µm [red]), it looks realistic in many situations, so that we refer to it as “pseudo-true color”. Snow and ice show up as an unrealistic blue, however, which is the main difference between it and a “true color” image. You might also notice a few more differences between the “pseudo-true color” image above and the “true color” image below.

True color RGB composite of VIIRS channels M-3, M-4 and M-5 taken 14:41 UTC 6 June 2012

True color RGB composite of VIIRS channels M-3, M-4 and M-5 taken 14:41 UTC 5 June 2012

The true color image uses moderate resolution channels M-3 (0.48 µm, blue), M-4 (0.55 µm, green) and M-5 (0.67 µm, red), which actually observe radiation in the blue, green and red portions of the visible spectrum. Apart from differences in resolution, the vegetation on the islands shows up a bit better in the “pseudo-true color” image. The islands just look brown in the true color image.

What is particularly interesting about these images are the visible effect that the islands have on the local atmosphere. Downwind (southwest, or to the lower left) of Sal, Boa Vista, and Maio, you can see singular cloud streets, much like the flow of water around a rock. In the photograph in that link, you can see how the water dips downward on both sides of the center line downstream of the rock, and upward in the middle (along the center line). The islands are acting like rocks in the atmosphere, causing upward motion behind them, and this lift was enough to form cloud streets. On either side of these cloud streets there is downward motion and, as a result, clear skies.

Downwind of São Nicolau, São Vicente and Santo Antão, the cloud streets highlight von Kármán vortices and vortex shedding, which you can see in more-controlled lab conditions here and here.

Many of the islands appear to be producing their own aerosol plumes (i.e. dust), and if you zoom in on the area between Boa Vista and Santiago, you can see gravity waves present in some of the plumes (highlighted by the arrows in the image below).

False color RGB composite of VIIRS channels I-1, I-2 and I-3 taken 14:41 UTC 5 June 2012

False color RGB composite of VIIRS channels I-1, I-2 and I-3 taken 14:41 UTC 5 June 2012

A common way to detect dust is the “split-window difference”: the difference in brightness temperature between the 11 µm channel and the 12 µm channel. On VIIRS, this means subtracting M-16 from M-15 which, when you do that, gives you this image:

Split-window difference from VIIRS (M15 minus M16) from 14:41 UTC 5 June 2012

Split-window difference from VIIRS (M15 minus M16) from 14:41 UTC 5 June 2012

The color scale goes from -0.16 K (black) to +4.0 K (white). For some reason, the dust or aerosol plumes don’t produce a strong signal here. It may be that the dust is too low in the atmosphere and the lack of temperature contrast with the surface prevents a strong signal. Maybe water vapor absorption effects in M16 are washing out the signal. Or, there could be some other explanation waiting to be discovered.

The plumes are highly reflective in the 3.7 µm channel (M-12), as are the clouds, which show up as warm spots in the image below (not as warm as the islands, however):

Moderate resolution 3.7 µm image (M-12) from VIIRS, taken 14:14 UTC 5 June 2012

Moderate resolution 3.7 µm image (M-12) from VIIRS, taken 14:41 UTC 5 June 2012

Here, just to throw you off, the color scale has been reversed so that dark colors mean higher values. The scale ranges from 295 K (white) to 330 K (black). When you take the difference of this image and the 10.6 µm brightness temperature (M-15), the clouds and aerosol plumes really show up, along with the gravity waves and vortices:

Brightness temperature difference between VIIRS channels M-12 and M-15 from 14:14 UTC 5 June 2012

Brightness temperature difference between VIIRS channels M-12 and M-15 from 14:41 UTC 5 June 2012

In this case, the M-12 brightness temperatures are always greater than the M-15 brightness temperatures (due to the combination of Earth’s emission and solar reflection in M-12 as opposed to just surface emission in M-15), so the scale varies from +5 K (black) to +30 K (white). Higher (brighter) values on this scale show off where the most solar reflection occurs at 3.7 µm – the liquid clouds and aerosol plumes.

There are much more sophisticated ways of identifying dust and aerosol plumes. To find out more, check out this article written by one of our resident experts, Steve Miller, who is currently working on applying dust detection algorithms to VIIRS.

If you are more interested in the von Kármán vortices, NASA has put together a great page that you can visit here. If you take the original image in this post, zoom out and rotate it a little bit, you can get a sense of just how far the vortices extend from their parent islands:

False color RGB composite of VIIRS channels I-1, I-2 and I-3 taken 14:41 UTC 5 June 2012

False color RGB composite of VIIRS channels I-1, I-2 and I-3 taken 14:41 UTC 5 June 2012. This image has been rotated from the previous images to highlight the length of the vortex streets.

Coincidentally, this image has been cropped to a size that makes it suitable for use as a desktop wallpaper, should you happen to have a 16:9-ratio monitor and a desire to stare at this image all day. (You have to click on the image, then click on the “1920 x 1080” link below the header to get the full resolution image.)