The next stage in my earthquake analysis project is to partition the events into groups with epicentre over land or water.
Since our existing catalog contains the latitude and longitude for the epicentres, it was a relatively simple matter to pipe these into gmtselect and label the events accordingly. The resulting data when sucked into R looks like this:
land.ocean[sample(1:nrow(land.ocean), 10),]
id surface
9776 pde20071224224102830_10 water
14462 pde20050127160223910_30 water
12454 pde20060716182649380_35 water
16737 pde20111122021806520_97 land
19442 pde20050110184730180_31 land
10613 pde20070802032440340_14 water
9761 pde20071227155344630_84 water
18537 pde20070825220549600_10 land
4428 pde20110214102127320_34 water
16914 pde20110521220626390_37 land
These data were then merged into the existing catalog using the id field.
A quick sanity check at this point would be good.
That looks about right: a clear distinction between land and ocean events.
Do the magnitude and depth of earthquakes on land differ significantly from those which occur under water? Looking at the plots below suggests that there is no qualitative difference.
Certainly it seems reasonable that an event happening deep beneath the surface of the Earth should not be affected by whether or not there is a (quite) thin layer of water on top! However, for my analysis of secondary effects it is probably going to be rather significant. But more on that later.