= GPlates Tutorial: Rotations == Authors - Grace E. Shephard^1^ - Kara J. Matthews^1^ - Jo Whittaker^1^ - R. Dietmar Muller^1^ ^1^ EarthByte Research Group, School of Geosciences, University of Sydney, Australia == Aim This tutorial is designed to teach the basics of implementing and manipulating rotations in GPlates. You will employ and build on the skills you acquired in the Plate Reconstructions tutorial. == Included Files For this part of tutorial you will need the data bundle found http://www.earthbyte.org/Resources/GPlates_tutorials/Rotations_Tutorial/SampleData/Rotations_Tutorial_Data.zip[here]. Unarchive the bundle to a folder named +GPlates_Rotations_Tutorial+, and remember the location so you can access the datasets later. The bundle includes the following datasets: . Rotation Models: +AusAnt_ExampleRotation.rot+, +Global_EarthByte_GPlates_Rotation_AusAnt_Example.rot+ . Coastline File: +Global_EarthByte_GPlates_Coastlines_20091014.gpml+ . Fracture Zone File: +AusAnt_FZs.gpml+ Additional datasets and files, including rotations, coastlines, isochrons, spreading ridges and COB files can be downloaded at the following site: http://www.earthbyte.org/Resources/earthbyte_gplates.html[] == Exercise A - Applying a Rotation In the Plate Reconstruction tutorial you learnt how to move features around on the globe using the 'Modify Reconstruction Pole' tool, we used the example of moving South America to fit together with Africa, like a jig-saw. It is more likely for tectonic reconstructions that you will want to change the rotation of a feature or apply a rotation to a feature back in time, rather than changing anything at the present-day. In this exercise we will learn how ot implement a rotation. We will start off with a very simple example, using a very small rotation file. Australia started to move away from Antarctica about 83 million years ago (Ma). According to Tikku and Canda (1999) Australia moved in a northward direction relative to a fixed Antarctica. You will implement this rotation in the rotation file you have been provided, +AusAnt_ExampleRotation.rot+. This file only contains rotations for Australia and Antarctica, for simplicity and so you can learn the basics of implementing rotations. - Open +AusAnt_ExampleRotation.rot+ in a text editor so you can see what it looks like (Figure 1): . Plate ID 000 = Spin axis . Plate ID 001 = Atlantic hotspots . Plate ID 701 = Africa . Plate ID 801 = Australia . Plate ID 850 = Tasmania . Plate ID 802 = Antarctica image::rotation_tut_image_1.JPG[] _________ *Figure 1:* The contents of our example rotation file. _________ You will notice that Australia moves relative to Antarctica, Antarctica moves relative to Africa, Africa moves relative to the hotspot reference frame which is fixed to the spin axis. - Open *GPlates* - *File -> Open Feature Collection*, select +AusAnt_ExampleRotation.rot+ and the EarthByte coastline file from the +GPlates_Rotation_Tutorial+ directory, and -> *Open* image::rotation_tut_image_2.JPG[] _________ *Figure 2:* Loading feature collections into GPlates. _________ - Rotate the globe so that it is centred on Australia. Now reconstruct backwards through time. You will notice that the feature data stay fixed in their present-day locations (this is because they have no relative rotations.) The only thing that changes is that features will disappear if you reconstruct before their "appear time". It is generally believed that Australia moved northwards, relative to a fixed Antarctica, between about 83Ma and the present (Tikku and Cande, 1999). We will implement this rotation. - Centre your globe so that Australia and the coastline of Antarctica nearest Australia are in view (Figure 3). image::rotation_tut_image_3.JPG[] _________ *Figure 3:* View of Australia and Antarctica. _________ - As we want to implement a rotation at 83 Ma, jump to this timem using the Time controls. - Use the 'Choose Feature' tool image:choose_feature.JPG[] to select Australia and then click 'Modify Reconstruction Pole' image:MRP.PNG[] - Drag Australia in a southward direction so that it approximately lines up with Antarctica (Figure 4). Once the feature attains the desired positions and orientation, clicking 'Apply' (right of the globe), this will open up 'Apply Reconstruction Pole Adjustment window', where you can review the details of your implemented rotation (Figure 5). image::rotation_tut_image_4.JPG[] _________ *Figure 4:* Australia has been dragged southward at 83 Ma to line up with Antarctica. _________ image::rotation_tut_image_5.JPG[] _________ *Figure 5:* Apply Reconstruction Pole Adjustment window, where you can review the details of your rotation implementation. _________ - In this window you can verify the new relative pole and details (Figure 5). Click *OK*, this will implement your rotation. You will notice that Australia is now positioned adjacent to Antarctica at 83 Ma (Figure 6). image::rotation_tut_image_6.JPG[] _________ *Figure 6:* Australia is now positioned adjacent to Antarctica at 83 Ma. It is south of it's present day position. _________ - Now you need to save your rotation file. *File -> Manage Feature Collections -> Save a copy* of the rotation file with a new name (this is so you can compare it to the old rotation file). Now load this new rotation file by clicking *Open File* and navigating to the directory where it is saved, -> *Open*. - Use the Animation slider to reconstruct from 83 Ma to the present. You will see Australia move in a northward direction relative to Antarctica! - However there is one more thing we need to do. If you jump to 600 Ma for example and animate back to 0 Ma, you will notice that Australia starts in its present day coordinates, moves southward to its 83 Ma position and then moves northwards again. This is because the location of Australia at 600 Ma is the same as present-day in our rotation file (Figure 7). We need to alter the rotation file so that there are not rotations between 600 Ma and 83 Ma. image::rotation_tut_image_7.JPG[] _________ *Figure 7:* A rotation has been added for Australia at 83 Ma. However notice that the latitude and longitude of Australia at 600 Ma is the same as present-day. _________ - Open the new rotation file in a text editor. Make the 600 Ma rotation data (lat, long, rotation angle) for Plate ID 801 the same as the 83 Ma rotation (ie duplicate the date, see Figure 8). This will result in no rotation between Australia and Antarctica until the period 83 Ma to the present. image::rotation_tut_image_8.JPG[] _________ *Figure 8:* Modified rotation file, note that the 600 Ma and 83 Ma rotations for Plate ID 801 are the same. _________ - Load your modified rotation file into GPlates and animate forward in time from say 150 Ma. You will notice that Australia stays attached to Antarctica until 83 Ma. == Exercise B - Modifying an Existing Rotation In this second exercise we will learn how to modify an existing rotation file. Keeping with the theme of Australia and Antarctica we will implement a new rotation for Australia at 83 Ma. Whittaker et al. (2007) proposed that 83 Ma Australia was located further eastwards with respect to Antarctica than previously thought (e.g. Tikku and Cande, 1999). They suggest that from about 83 Ma to 50 Ma Australia moved northwest relative to a fixed Antarctica, before then commencing northward motion between 50 Ma and the present. The timing of their proposed change in plate motion at 50 Ma coincides with the Hawaiian-Emperor bend and subduction of the Izanagi-Pacific ridge. - Eject all existing rotation files from GPlates. *File -> Manage Feature Collections* -> click the *Eject* symbol image:eject.JPG[], corresponding to all loaded rotation files. Keep the 'Manage Feature Collections' window open. - *Open File* -> select the rotation file for this exercise, +Global_EarthByte_GPlates_Rotation_AusAnt_Example.rot+ -> *Open* The rotation file we have just loaded is significantly more complicated that that of the last exercise. Reconstruct the globe back through time and you will see that all the plates move. If you open this rotation file in a text editor you can see how much longer and more detailed it is compared to the last exercise. - Use the 'Time Controls' to jump to 83 Ma. We will use the fracture zones to help us constrain the position of Australia at 83 Ma. - *File -> Manage Feature Collections -> Open File* -> select +AusAnt_FZs.gpml+ from the data bundle -> *Open* Following Whittaker et al. (2007) we will align the Antarctic fracture zone with the most westerly Australian fracture zone, whereby shifting Australia east relative to a fixed Antarctica. - Use the 'Choose Feature' button image:choose_feature.JPG[] to select Australia -> click 'Modify Reconstruction Pole' image:MRP.PNG[] -> drag Australia eastward so that the Antarctic fracture zone and the western Australian fracture zones are connected (Figure 9) -> click *Apply* (right of globe) -> you can then click *OK* in the 'Apply Reconstruction Pole Adjustment window' once you have reviewed the details of your new reconstruction and are satisfied. image::rotation_tut_image_9.JPG[] _________ *Figure 9:* Australia shifted east using the Modify Reconstruction Pole tool. _________ When you return to the globe you will notice that Australia is located further east than when you started (Figure 10). We now need to save this data. image::rotation_tut_image_10.JPG[] _________ *Figure 10:* Australia shifted further east 83 Ma. _________ - *File -> Manage Feature Collections* -> save your +Global_EarthByte_GPlates_Rotations_AusAnt_Example.rot+ file with a new name image:save_copy.JPG[] so that you preserve the old rotation file. - Now use the Time controls to watch Australia's motion from 83 Ma to present-day and you will see that there is northwest motion of Australia relative to a fixed Antarctica between 83 Ma and 50 Ma. Then Australia commences northward motion. - Open your modified rotation file and the original rotation file using a text editor and scroll down to the entries for Plate ID 801, and compare the two rotation files. You will see that they have different entries now for 83 Ma (Figure 11). image::rotation_tut_image_11.JPG[] _________ *Figure 11:* New (top) and old (bottom) rotation files showing entries for Plate ID 801. Entries for 83 Ma have changed. _________ NOTE: To better appreciate the change in motion of Austalia relative to a fixed Antarctica, you can specify Antarctica as the "anchored plate" rather than the spin axis (default). *Reconstruction -> Specify Anchored Plate ID* -> enter +802+. Now when you reconstruct the globe you can really notice that Australia moves in a northwest direction between 83 Ma and 50 Ma. === Things to consider: The cursor provides longitude and latitude locations to help with re-orienting. This is particularly useful when trying to replicate work from other literature. Check the existing rotation file for the time increments for the plates. Reconstructing at these times will avoid jumps between two time steps. For example if the existing rotation file has rotations at 10 Ma and 20 Ma, by creating a new rotation at 16 Ma will only change the rotation between 10 Ma and 16 Ma. Between 16 Ma and 20 Ma the plate may jump erratically according to the old pole of rotation, unless you change it or an older timestep. == Additional Exercise - Rotation of Borneo (Introducing how to reverse engineer rotations from images) In this final exercise we will move further afield to Borneo, Southeast Asia. The rotational history of Borneo (Kalimantan) is a contentious issue in the literature. Hall (2002) for example, impose a 45 degree counter-clockwise rotation of Borneo from 26 Ma to 10 Ma (Figure 12), whereas Lee and Lawyer (1995) prefer no rotation of this block, with only a small latitudinal change. Another extreme is Replumaz and Tapponier (2003) who favour a large clockwise rotation of Borneo (Figure 13), which appears to be incompatible with paleomagnetic data. However, the lack of large strike-slip faults around Borneo, which would be expected under significant rotation, suggest that Borneo experienced little rotation in the Miocene. Nevertheless this remains a debated topic. image::rot2_tut_image_6.JPG[] _________ *Figure 12:* Reconstructions of the Sundaland region at 35Ma, 23Ma, 15Ma and 5Ma from Hall (2002). Notice the large anti-clockwise rotation of Borneo (red arrow). _________ image::rot2_tut_image_7.JPG[] _________ *Figure 13:* Replumaz and Tapponnier (2003) reconstructions of SE Asia between 10 Ma and 15 Ma. _________ There are a couple of approaches to reconstruct a rotation as derived from an image. - One way is to best discern the latitudinal and longitudinal locations of the feature manually, for example, by selecting a distinguishing part of the tectonic shape on the image (e.g. a peninsula) that is also replicated in the coastline file. The rotation should be manually changed in *GPlates* (See Exercise A) until the correct coordinates and orientation are achieved. Remember to reconstruct the rotation at the correct time. - The second way is to load the image as a raster in *GPlates* (See the Importing Rasters tutorial), and reconstruct the rotation with the image as an underlying guide. - A third way is if you are provided with the three rotation pole values for a particular plate at a particular time period, with the conjugate plate specified. With this information you can manually enter it into the rotation file. .Methods: - Rotate the window to SE Asia and select Borneo. - Select the relevant reconstruction time. - Rotate Borneo, using the skills you have developed in the previous exercise, to the desired location and orientation as determined by method of reconstruction e.g. raster loading or manual adjustment. - Once the feature attains the desired position and orientation, clicking 'Apply' will transfer these changes to the rotation file. - To compare models you need to save the different reconstructions to different rotation files. == References Hall, R. 2002. Cenozoic geological and plate tectonic evolution of SE Asia and the SW Pacific: computer-based reconstructions, models and animations. Journal of Asian Earth Sciences, 20; 353-431 Lee, T.Y., and Lawver, L.A. 1995. Cenozoic plate reconstructions of Southeast Asia. Tectonophysics. 251; 85-138. Replumaz, A. and Tapponnier, P. 2003. Reconstruction of the deformed collision zone between India and Asia by backward motion of lithospheric blocks. Journal of Geophysical Research. 108; 2285.