Let's do something with Plate Tectonics and ask how fast plates move relative to one another? The answer can be found by using plate rotation vectors. Stick with me for a minute or two. This looks more complicated at first blush than it really is. For our purposes, we just need the ability to plug numbers into an equation --so we need to follow the parameter definitions and the equation.

### Want to learn more?

See "Simple Euler Poles."

The motion of a point on one tectonic plate relative to another plate can be described by the relative velocity vector *v*. The velocity *v* has magnitude and direction and is given by the cross product of the angular velocity vector ω and the plate rotation vector** r: v = ω x r **

For example, according to one of the accepted models for plate motion (NUVEL 1), the velocity of the North American Plate relative to the Pacific Plate is given by the rotation pole at: 48.7° N 78.2° W and angular velocity 7.8e-7 deg/year (note that you should read 7.8e-7 as 7.8 times ten raised to the negative 7, that is: 0.00000078 deg/year.) Therefore, a point on the Pacific plate near Parkfield California, which is at 35.9° N 120.5° W, is moving at 47.6 mm/yr relative to the rest of North America. How long will it take for this point to reach the present location of San Francisco?

How does this calculation work? Download this pdf file for the details. That file contains some useful background. The last page is the example above.

**NOTE:** Parkfield CA is the site of a National Science Foundation project called EarthScope that has drilled into the San Andreas Fault. See SAFOD Observatory for more details on the drilling project.

#### One of the first sketches to show tectonic plates

Note the three types of plate boundaries (compare to the figure on the previous page) and the definitions of lithosphere, asthenosphere, and mesosphere. Lithosphere means * the rigid * part and thus the bottom of it is defined by an isotherm (do you know why?). The base of the lithosphere is typically taken as 1300° C. Note that the plate thickens as it moves away from a divergent spreading center. Mid-ocean ridge systems are hot (they are volcanoes!) and thus ridges are relatively buoyant, which means that they have relatively higher elevation than regions around them. Ocean depth increases systematically with distance away from mid-ocean ridge systems. We'll look at this more closely in Activity 3.

#### Plate motion vectors

You can always use vector algebra to calculate linear velocity *v* from the position vector * r * and the angular velocity vector *ω,* but there's an easier way to get the magnitude of the velocity by using the solid angle between the pole of rotation and the location of interest (see below). The solid angle can be obtained using spherical trigonometry:

cos a = cos b cos c + sin b sin c cos A

where a is the solid angle of interest, b is the co-latitude of the location on Earth's surface, c is the co-latitude of the plate rotation pole and A is the surface angle between the pole and the location (that is: A is the difference between the longitude of the pole and the longitude of the location).

#### Some background on Spherical Trigonometry

To work with plate motion vectors, and to calculate the linear velocity of points on Earth's surface, we need to know the distances between various points on the globe. A useful analogy is that of linear and angular velocities associated with Earth's daily rotation. That is, the angular velocity is the same everywhere on Earth. All points rotate through 360° (2 pi radians) in 24 hours. But the linear velocity, on Earth's surface, depends on where you are relative to the rotation axis. If you're at the North Pole, then you cover only a small distance, whereas if you're at the equator, then you cover a distance equal to Earth's full circumference in 24 hours (2 pi R). As Earth rotates each day, the linear velocity of points at the Equator is much larger than points near the poles. The same type of thing happens with plate motions. Points that are close to the pole of rotation move with lower linear velocity than points that are farther from the pole. So, we need to calculate the distance between each point and the pole. These next two figures will help show how this works. Remember, for our purposes, we just need to be able to plug numbers into an equation, so we need to follow the parameter definitions and the equation.

In the diagram above, upper case letters refer to surface angles and lower case letters refer to solid angles, measured between lines that extend from the Earth's center to the surface. For a point X at, say, latitude 20° N, the angle b is 70°, because b is measured from the north pole along a line of longitude. In the calculation, it's standard to use the 'co-latitude' b and c. Note that it's easy to get b and c, based on their latitudes. But the same is not true for the solid angle a. That's why we need spherical trig. Surface angles are perhaps more familiar. They are obtained from latitude and are therefore nothing more than a larger-scale version of the angle between the first-base line and the third-base line on a baseball diamond.

Here's an example that will help to fix ideas. Do you follow? If not, please post a question on Canvas.

*Click for text description of the spherical trigonometry example image.*

What is the magnitude of the linear velocity of the Eurasian plate w.r.t State College?

$${}_{NA}{\omega}_{E}=2.2\times {10}^{-7}\text{deg/yr,62.4}\xb0\text{N,135.8}\xb0\text{E}$$State College PA: 40.8°N, 77.9°W (-77.9°)

Use spherical trig identity

$$\mathrm{Cos}\text{}a=\mathrm{cos}\text{}b\text{}\mathrm{cos}\text{}c+\mathrm{sin}\text{}b\text{}\mathrm{sin}\text{}c\text{}\mathrm{cos}\text{}A$$In our notation:

$$\mathrm{Cos}\Delta =\mathrm{cos}{\theta}_{x}\mathrm{cos}{\theta}_{p}+\mathrm{sin}{\theta}_{x}\mathrm{sin}{\theta}_{p}\mathrm{cos}\left({\varphi}_{p}-{\varphi}_{x}\right)$$**Check: Can you verify 73.3 deg. for this example? If not, make sure you're using co-latitude and that your answer is in degrees. Still having trouble? Then have a look at this.**

#### Magnitude of linear plate velocity

Once you have the angular distance between the points (Δ), you can get the linear velocity using * v = ω R sin Δ. * See the last page of this pdf file for a worked example.

### Activity 1

**NOTE:** For this assignment, you will need to record your work on a word processing document. Your work must be submitted as a Word (.doc), or PDF (.pdf) format so I can open it.

#### Directions

- Find a suitable analog for Earth's internal structure and measure the relative thicknesses of the layers. Here is an image that goes along with mine (a major league baseball), see details of my numbers on the previous page:

(Are you a fan? Here's a short set of ppt images showing how baseballs have changed over the years.) - Using the formula given on the previous page (Earth is not a perfect sphere), calculate Earth's radius at your location. Note, this is just a sort of 'warm up' for us. We could expand on this and account for differences in radius as a function of position, but we won't go that far here. I thought you might like to see how much the radius changes as a function of latitude.
- Now, we're ready for a real problem. Calculate the rate of motion between the Pacific Plate and the North American Plate at the following two locations. I expect that you will show enough of your work so that I can follow how you did the calculation. At the least, you should list all of the parameters and factors used in your calculation. Even better is to paste in your calculation from Google or elsewhere, as in this example: http://www3.geosc.psu.edu/~cjm38/540/sphericaltrigExample.html
- 37.8 N 122.5 W (the golden gate bridge)
- 34.1° N 118.3° W (Hollywood, CA)

- Check your answer using the plate motion calculator at the UNAVCO Web site.
**NOTE:**To start, you can enter just the latitude and longitude of the point of interest and hit submit. You'll get an answer, with default parameters. Hmm, does it work to just copy/paste in the numbers above? What happens if you write 122.5 W vs. -122.5? Can you include the ° symbols?For our example problem, you should set the Reference to "PA Pacific" when doing the Golden Gate case, and "North America" when doing Hollywood (what happens if you choose Pacific for the Hollywood case?). You can also try

**NNR (no net rotation). Play around with this a bit. It's useful!**Under "Model" select "All of the above" so that you can see the range of predictions. Tell me your thoughts on why there are differences in the predicted rates of motions. What happens if you use a different frame of reference?

- Go to the UNAVCO Web site and determine a plate motion rate using their GPS data for the Mission Viejo, CA station. The link will take you to the UNAVCO page for Mission Viejo, CA (SBCC). Start by looking at the Overview tab. Note that this page includes plots like the one you need to make. Now go to the Data Products tab and download the "Unfiltered Data." You'll probably want CSV format, but that's up to you. You may have to play around here a bit, but that's part of this exercise. Keep in mind that this station (SBCC) is moving in time --it's recording motion of the North American Plate! The goal is for you to plot that motion (distance vs. time) and to compare the measured rate to the predictions based on your calculations using plate rotation vectors. (Note: The folks at UNAVCO sometimes modify their web page. If you have trouble finding the data, use this file). Once you have the data, the idea it to plot things and fit a line to each of the components (North, East --note that the site has absolute values of position and also deviation; I would suggest using absolute values, but you can use either because we just want the rate of change). Then determine a rate of motion along the orientation of the plate boundary (which is roughly: 300° cw from North). How does this value compare with the prediction from the Plate Motion Calculator for this latitude?
- Save your document as either a Microsoft Word or PDF file in the following format:
L2_activity1_AccessAccountID_LastName.doc (or .pdf).

For example, student Elvis Aaron Presley's file would be named "L2_activity1_eap1_presley.pdf"—this naming convention is important, as it will help me make sure I match each submission up with the right student!

#### Submitting your work

Upload your paper to the "Lesson 2 - Activity 1" dropbox in Canvas by the due date indicated on our Course Schedule.

#### Grading criteria

See the grading rubric for specifics on how this assignment will be graded.