Wednesday, December 29, 2004

Banda Aceh:

Today's Sumatra -Andaman Islands event was one of the great earthquakes - we will see only one or two like this in our lifetimes. The combination of the location and orientation made for a horrendous result - several tens of thousands dead, with a rising death toll guaranteed. (This post took a while to finish - so it contains material from well after the event...)

First, the tectonic setting (click on the image for larger version):

The India and Australia plates are moving northward, and are shoving under (subducting) the Burma microplate. Here's a schematic of the setting, with the India/Australia plate on the left, moving to the right and under the Burma microplate on the right:

The movement occurs in sudden bursts -- earthquakes. Like banging a tray full of water, if the earthquake moves the seafloor, it sets off waves -- tsunamis.

Here is one of the early simulations done of the tsunami by the Japan National Institute of Advanced Industrial Science and Technology:

And here is another simulation of the Sumatra-Andaman tsunami done by NOAA. Notice the reflections off and refractions around India, Diego Garcia, the 90 East Ridge and the Maldives! These reflections and refractions can focus energy and cause the tsunami to be very high on some shores and not so high on others - so the threat is not uniform along all shorelines. Note that Australia's West coast is much closer to the initial movement than Kenya, but it got almost no resulting wave.

The Geological Survey and Caltech's Seismo Lab produced a first estimate of the amount of the slip along the fault of the great Sumatra-Andaman Islands earthquake:

The red star is the location of the hypocentre - where the earthquake began. The black line is the surface location of the subduction fault (the trench), colours show the amount of displacement in cm of the hanging (upper) wall of the fault relative to the foot (lower) wall. Note that most maps show the hypocentre as if this was where all the action was, but that a lot of the motion of the seafloor was probably several hundred kilometers to the Northeast (and in fact, later modelling shows that there was significant motion even up into the Andaman and Nicobar Islands). This difference between the location of the hypocenter and the main displacement values is typical - and misleading.

Now it becomes obvious why the North-West of Aceh got hit so hard - the largest slip is 200km farther up the peninsula than the hypocentre.

Note that even in this preliminary model, an area of about 1200 km by 100km has moved several meters. This resulted in several meters of seafloor uplift (and displaced water, creating the wave). The total energy released was estimated at 3.57X10^29 dyn.cm in 200 seconds, with a peak slip of about 20 meters. (Chen Ji, Caltech).

There were at least fourteen aftershocks larger than magnitude 6.0, and one larger than 7.0 (all of these respectable 'quakes on their own).

This was enough to cause a 22cm wave in San Diego harbor, 2.6 meters in Manzanillo, Mexico due to focusing effects, and several days of tidal height disturbances as the sloshing continues to reverberate around the ocean basins.

Movie of the Andreanof tsunami movie
- moves much more slowly in the shallow waters north of the Aleutians
- reflects off/refracts through the Hawaiian/Emperor chain - the BC coast gets hit twice: once directly, the second time as a reflection off the Hawaiian seamounts.

What it looked like on the ground in Thailand (Aftonbladet, 8.0 Mb .mov file)This is a different angle from a clip that received a lot of air time. During the second surge, off to the lower left while looking out to sea, one can see a couple behind a palm tree trying to get out from behind a pile of debris blocking the surge and into the bar area. As with the shot from the other angle, I'm not sure they made it before the debris collapsed onto them, and pinned them under tons of rushing water. This is a horrendous example of how innocent looking water can quickly turn deadly. This begins as an interesting video, and quickly becomes quite frightening.

Links:

There is a lot of talk now about an Indian Ocean Tsunami warning system, and support from the U.S. Government for a global system that will include the Atlantic too. Perhaps a good idea, but several things should be noted. The frequencies of tsunamis in these two other oceans are very much less than in the Pacific - so there is a very high probability that the systems will never be used. On the other hand, the process of preparation of such a system reduces societies' vulnerability to many different types of hazards at the same time, so while we may not use them for saving people from tsunamis, we may end up using them to save people from other types of disasters. The last thing is a comment regarding tsunamis - even if there had been a fully operating system in place (which includes not only an earthquake detection system and a wave detection system, but more critically an information dissemination and public evacuation plan), the system would have done nothing to save Aceh or possibly even Thailand, since the wave would have reached them before any possible plan could have been activated.

Here is an amazingly prescient paper published in the September 2004 issue of Geoscience Australia magazine about an Indian Ocean tsunami warning system: AusGeo News 75 p.4 (part 1, with the referenced article, is a 3.9 Mb PDF file)

The best advice: if you see the sea withdraw, you can be guaranteed that it will be coming back with a vengeance. Run for the hills.

And finally, a word about feral goats in the Andamans, who, I would bet, are still there.

Wednesday, December 22, 2004

Oh no, a physics post! It even has... formulas!

Our day-to-day experience with weights and speeds is limited to a small range. We don't expect the behaviour of objects to change very much, even when the object is far outside this limited range of weight or speed.

We know very well that a small piece of lead (between 2 and 30 grams), travelling very fast (between 300 and 1,000 meters per second), can do a great deal of damage - i.e. a bullet.

However, when we think of a large piece of polystyrene foam (about 30 cm across), then all bets are off. After all, it's hard to think of being able to hurl a large piece of foam fast enough to do damage to anything, right?

I went to a lecture the other day by Doug Osheroff that demonstrated this very clearly, using some high-school level physics.

Doug has a Nobel Prize. He's a low-temperature physicist who works at Stanford, and who recently served on NASA's Columbia Accident Investigation Board. So I paid attention to what he was saying.

He outlined the basic calculations behind the conclusion that a 0.75 kg chunk of foam from the Shuttle's external tank (ET) had irreparably damaged the leading edge of Columbia's left wing.

Video of the launch taken from various angles shows the foam breaking off and striking the wing - this allows an estimation of the source of the foam (the -Y bipod ramp, and therefore its size & weight, +/- 30cm radius, 0.75 kg), and also the amount of time it took to travel the distance from the ramp to the leading edge (19 meters in about 0.16 seconds).

The key here is to realize that since the foam fell off the external tank, it was decelerating relative to the air rushing by at Mach 2.5, and accelerating relative to the shuttle (and its onrushing wing).

So:
d0 = ½at02


Where d0 is the distance travelled, a is the acceleration of the foam, and t0 is the time between the foam breaking off and the wing strike.

Substituting the values:

19 m = ½ a (0.16 sec)2


gives

a = 1,480 m/s2


or about 150 g's.

How fast was that foam travelling when it got to the wing?

v0 = a t0


substituting values:

v0 = 1,480 m/s2 · 0.16 sec


gives

v0 = 238 m/s


Converting to something we understand better: 855 km/h or 530 mph. Pretty fast.

But it's still hard to comprehend a chunk of foam doing such tremendous damage.

In terms of kinetic energy,

Ek = ½ m v02


plugging in our values:

Ek = ½ 0.75 kg· (238 m/s)2


gives

Ek = 21,240 N·m (or J)


Say we had a bullet with equivalent kinetic energy - how fast would it have to be going? Let's choose the a fairly standard .22 bullet of about 40 grains weight, or about .0026 kg

Then:

21,240 J = ½ 0.0026 kg · v12


gives

v1 = 4,040 m/s


or about 13,260 fps, which is about five times faster than a bullet coming out of a rifle - or about twenty-five times the energy a regular 0.22 bullet would pack.

OK, obviously were talking something serious here.

And what was it crashing against? Aren't those wings tough? After all, they can take flying through plasma, right?

Umm, well, no. In fact, the leading edge material is really fragile. The RCC panels on the shuttle have to be covered with protective padding when the shuttle is being serviced to avoid damage from dropped tools or even dropped bolts.

The plasma-resistant wing never stood a chance against a chunk of foam.