In this exercise I would like you to make a map to highlight the concept of "plate tectonic motion" for the present day Earth. It would be good to show the difference between the seismicity of mid ocean ridges and subduction zones, perhaps showing the strain rate patterns.
It may also help to plot plate motion arrows on the map to highlight divergent v. convergent zones. In order to do this "intuitively" you probably should choose the NNR (no net rotation) velocity field which tends to show major ridges as spreading in opposite directions.
You should start with the Working with point data notebook. This shows you how to grab seismic data into an array and plot it on a map. You could look at the obspy query functions in the Himalayas notebook if you want to get different data.
You will need to think about:
what colours to use and how to draw out the difference between subduction zones earthquakes and mid-ocean ridges (what is the main difference ?)
which basemap is the best choice (an image, a coastline, shaded ?)
Hint: plate motion arrows (see Global Plate Motions notebook) ?
Make sure you add comments to your code and intersperse the python with markdown explanations.
I expect to be able to run this notebook before I read your answer !!
In [ ]:
## Add your code here
In [ ]:
In [ ]:
## Plot your map here !
In [ ]: