GPlates: A Beginner’s Guide to Plate Tectonics Visualization

GPlates Workflow: From Data Import to Animated Reconstructions

Overview

This article shows a concise, practical workflow for using GPlates to import data, edit and build plate models, and create animated reconstructions. Steps assume GPlates 2.x on desktop (Windows/Linux/macOS) and common geospatial formats (shapefiles, GeoJSON, NetCDF, raster TIFF). Default choices are made so you can follow end-to-end without extra questions.

1. Prepare your data

  • Formats: vector (shapefile, GeoJSON, GML), raster (GeoTIFF), gridded data (NetCDF), and GPlates native (.gpml, .gpmlz) are supported.
  • CRS: ensure geographic coordinates (WGS84 / EPSG:4326). Reproject if needed.
  • Attributes: include time fields for time-dependent features (e.g., age, start_age, endage) or separate time-stamped files.
  • Topology: clean geometry (no self-intersections). Use GIS tools (QGIS, GDAL) to validate.

2. Import into GPlates

  • Open GPlates and choose File → Open Feature Collection or File → Add Raster Layer.
  • Load plate tectonic-related datasets:
    1. Plate polygons (if available) — used to assign plate IDs.
    2. Rotation model (.rot, .gpml) — essential for reconstructions.
    3. Time-dependent features (paleoshorelines, faults, hotspots) — include start/end ages in attributes or use separate files per time step.
    4. Raster/base maps (bathymetry, paleo topography) — optional for visual context.
  • Use the Data Manager (Layer Manager) to check layers and metadata.

3. Assign plate IDs and topological relationships

  • If you have plate polygons with plate IDs, ensure feature attributes match rotation model plate IDs.
  • If not, create plate polygons or assign plate IDs manually:
    • Use “Assign Plate IDs To Features” (Desktop → Plate Reconstruction Tools) to interactively assign plate IDs to vector features.
    • For many features, use attribute-based assignment via attribute table editing or batch scripts (Python/pygplates).
  • Verify topology: adjacent features on the same plate should share the same plate ID.

4. Load and verify rotation model

  • Load rotation files via File → Open Rotation Model. Common formats: .rot, .gpml, .gpmlz.
  • Inspect the rotation poles and ages in the Rotation Manager.
  • Check that plate IDs in the rotation file match your feature plate IDs. If not, edit either the rotation file or the feature attribute table to align IDs.

5. Time settings and reconstruction parameters

  • Set reconstruction time in the Time Slider (bottom of the GPlates window) or use the Reconstruction menu for batch runs.
  • Choose reconstruction options:
    • Reconstruction anchor plate (if needed).
    • Small-circle vs great-circle interpolation.
    • Topology preservation and feature splitting rules.
  • For animations, decide temporal resolution (e.g., 1 Myr, 0.5 Myr).

6. Preview reconstructions interactively

  • With time slider, scrub through ages to preview how features move.
  • Use the “Reconstruct” context menu on a layer to display the reconstructed state at the chosen time.
  • Check for artifacts: misplaced features, overlaps, or gaps due to incorrect plate IDs or rotation inconsistencies.

7. Edit and refine

  • Fix issues by:
    • Reassigning plate IDs.
    • Adjusting rotation poles or ages (edit rotation file).
    • Correcting geometry (split/merge features).
  • Re-run interactive preview to confirm fixes.

8. Batch reconstructions and export

  • For many time steps or full animations, use Reconstruction → Reconstruct All Feature Collections…:
    • Set start/end ages and step size.
    • Choose output format (GeoJSON, ESRI shapefile, GPlates feature collections, or GMT).
  • Export reconstructed features for each time step to a folder.

9. Create animated reconstructions

Option A — GPlates built-in screenshot sequence:

  • Use Reconstruction → Export Visualization → Export Image Sequence.
  • Configure map projection, resolution, background, and which layers to include.
  • Export a sequence of PNG/TIFF frames for each time step.

Option B — External animation (recommended for advanced control):

  • Export reconstructed vector frames (GeoJSON or Shapefile) or raster frames.
  • Use external tools (QGIS Time Manager, GMT, or scripting with Python + matplotlib/Cartopy) to render frames.
  • Assemble frames into video with FFmpeg:

    Code

    ffmpeg -framerate 24 -i frame_%04d.png -c:v libx264 -pix_fmt yuv420p gplatesrecon.mp4

10. Add overlays and annotations

  • Add base maps (present-day coastlines, graticules) or paleo-environment rasters.
  • Annotate plate names, age ticks, or scale bars using GPlates labeling or in post-processing (QGIS or video editor).

11. Reproducibility and sharing

  • Save session as a GPlates project (.gproj) to preserve layer links and visualization settings.
  • Share rotation models and feature collections (.gpml or .gpmlz).
  • Document the time steps, interpolation method, and any manual edits.

Quick checklist

  • Data in WGS84
  • Plate IDs assigned and matched to rotation model
  • Rotation model loaded and verified
  • Time range and step size set
  • Previewed and corrected reconstructions
  • Exported frames or image sequence
  • Assembled final video (FFmpeg or editor)
  • Saved project and exported data for sharing

Example minimal command-line workflow (pygplates + FFmpeg)

  1. Use pygplates to reconstruct features for each time step and write PNGs (script omitted for brevity).
  2. Assemble with FFmpeg:

Code

ffmpeg -framerate 15 -i recon_%03d.png -c:v libx264 -pix_fmt yuv420p gplates_animation.mp4

Further reading: consult GPlates documentation and pygplates tutorials for scripting examples and advanced rotation editing.

Comments

Leave a Reply

Your email address will not be published. Required fields are marked *