{ "cells": [ { "cell_type": "markdown", "id": "81dd17c9-2794-4759-b7e9-eaaa99bc5e1b", "metadata": {}, "source": [ "# GPS velocities for the Alpine region\n", "\n", "This is a compilation of 3D GPS velocities for the Alps by\n", "[Sánchez et al. (2018)](https://doi.org/10.1594/PANGAEA.886889).\n", "The horizontal velocities are reference to the Eurasian frame. \n", "All velocity components and even the position have error estimates,\n", "which is very useful and rare to find in a lot of datasets.\n", "\n", "License: CC-BY\n", "\n", "Here, we download the data from 3 separate files (coordinates, vertical velocity, horizontal velocities) and make sure they are aligned and represent the same stations. There are some mistakes in the station names of horizontal velocity file that we fix manually (verified by the coordinates)." ] }, { "cell_type": "code", "execution_count": 1, "id": "28905da1-9159-477f-a0e0-04f80de6bf89", "metadata": {}, "outputs": [], "source": [ "import os\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "import pandas as pd\n", "import verde as vd\n", "import pooch\n", "import pyproj\n", "import pygmt" ] }, { "cell_type": "markdown", "id": "827d30b5-158e-4fe2-9e16-1de8e7ff894d", "metadata": {}, "source": [ "## Download the data\n", "\n", "Use Pooch to download and cache the data files." ] }, { "cell_type": "code", "execution_count": 2, "id": "41d23f77-ad1c-4c8c-b1a5-a7667e461e02", "metadata": {}, "outputs": [], "source": [ "fname_position = pooch.retrieve(\n", " url=\"https://store.pangaea.de/Publications/Sanchez-etal_2018/ALPS2017_NEH.CRD\",\n", " known_hash=\"sha256:24b88a0e5ab6ea93c67424ef52542d8b8a8254a150284e1a54afddbfd93e4399\",\n", ")\n", "fname_velocity = pooch.retrieve(\n", " url=\"https://store.pangaea.de/Publications/Sanchez-etal_2018/ALPS2017_NEH.VEL\",\n", " known_hash=\"sha256:0f2eff87a39260e2b3218897763dbfecdf0f464bf877bef460eff34a70e00aa7\",\n", ")\n", "fname_velocity_eurasia = pooch.retrieve(\n", " url=\"https://store.pangaea.de/Publications/Sanchez-etal_2018/ALPS2017_REP.VEL\",\n", " known_hash=\"sha256:578677246230e893c828205391d262da4af39bb24a8ca66ff5a95a88c71fe509\",\n", ")" ] }, { "cell_type": "markdown", "id": "10bb9087-c3d1-43e8-a11a-2891767df9e1", "metadata": {}, "source": [ "## Load the data\n", "\n", "These data are in a strange format and getting pandas to read it would be more work than parsing it by hand. So that's what we're going to do.\n", "\n", "First, the horizontal velocities since there are less points and we will only want the vertical and positions of these stations." ] }, { "cell_type": "code", "execution_count": 3, "id": "18e2f731-8ce1-4449-b8d2-3ade1696a66a", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
| \n", " | station_id | \n", "velocity_east_mmyr | \n", "velocity_north_mmyr | \n", "velocity_east_error_mmyr | \n", "velocity_north_error_mmyr | \n", "
|---|---|---|---|---|---|
| ACOM | \n", "ACOM | \n", "0.2 | \n", "1.2 | \n", "0.1 | \n", "0.1 | \n", "
| AFAL | \n", "AFAL | \n", "-0.7 | \n", "0.9 | \n", "0.1 | \n", "0.2 | \n", "
| AGDE | \n", "AGDE | \n", "-0.2 | \n", "-0.2 | \n", "0.1 | \n", "0.3 | \n", "
| AGNE | \n", "AGNE | \n", "0.0 | \n", "-0.2 | \n", "0.2 | \n", "0.6 | \n", "
| AIGL | \n", "AIGL | \n", "0.0 | \n", "0.1 | \n", "0.1 | \n", "0.5 | \n", "
| ... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
| WLBH | \n", "WLBH | \n", "0.0 | \n", "-0.2 | \n", "0.1 | \n", "0.2 | \n", "
| WTZR | \n", "WTZR | \n", "0.1 | \n", "0.2 | \n", "0.1 | \n", "0.1 | \n", "
| ZADA | \n", "ZADA | \n", "0.2 | \n", "3.1 | \n", "0.2 | \n", "0.4 | \n", "
| ZIMM | \n", "ZIMM | \n", "-0.1 | \n", "0.4 | \n", "0.1 | \n", "0.1 | \n", "
| ZOUF | \n", "ZOUF | \n", "0.1 | \n", "1.0 | \n", "0.1 | \n", "0.1 | \n", "
186 rows × 5 columns
\n", "| \n", " | station_id | \n", "latitude | \n", "longitude | \n", "height_m | \n", "latitude_error_m | \n", "longitude_error_m | \n", "height_error_m | \n", "
|---|---|---|---|---|---|---|---|
| ACOM | \n", "ACOM | \n", "46.547935 | \n", "13.514900 | \n", "1774.682 | \n", "0.0009 | \n", "0.0005 | \n", "0.001 | \n", "
| AFAL | \n", "AFAL | \n", "46.527144 | \n", "12.174517 | \n", "2284.085 | \n", "0.0009 | \n", "0.0009 | \n", "0.001 | \n", "
| AGDE | \n", "AGDE | \n", "43.296383 | \n", "3.466427 | \n", "65.785 | \n", "0.0018 | \n", "0.0009 | \n", "0.002 | \n", "
| AGNE | \n", "AGNE | \n", "45.467942 | \n", "7.139620 | \n", "2354.600 | \n", "0.0036 | \n", "0.0009 | \n", "0.004 | \n", "
| AIGL | \n", "AIGL | \n", "44.121398 | \n", "3.581261 | \n", "1618.764 | \n", "0.0009 | \n", "0.0009 | \n", "0.002 | \n", "
| ... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
| WLBH | \n", "WLBH | \n", "48.415171 | \n", "7.351299 | \n", "819.069 | \n", "0.0009 | \n", "0.0005 | \n", "0.001 | \n", "
| WTZR | \n", "WTZR | \n", "49.144199 | \n", "12.878911 | \n", "666.025 | \n", "0.0005 | \n", "0.0005 | \n", "0.001 | \n", "
| ZADA | \n", "ZADA | \n", "44.113177 | \n", "15.227590 | \n", "64.307 | \n", "0.0036 | \n", "0.0018 | \n", "0.004 | \n", "
| ZIMM | \n", "ZIMM | \n", "46.877098 | \n", "7.465278 | \n", "956.341 | \n", "0.0009 | \n", "0.0005 | \n", "0.001 | \n", "
| ZOUF | \n", "ZOUF | \n", "46.557221 | \n", "12.973553 | \n", "1946.508 | \n", "0.0009 | \n", "0.0005 | \n", "0.001 | \n", "
186 rows × 7 columns
\n", "| \n", " | station_id | \n", "velocity_up_mmyr | \n", "velocity_up_error_mmyr | \n", "
|---|---|---|---|
| ACOM | \n", "ACOM | \n", "1.1 | \n", "0.1 | \n", "
| AFAL | \n", "AFAL | \n", "1.3 | \n", "0.2 | \n", "
| AGDE | \n", "AGDE | \n", "0.1 | \n", "0.3 | \n", "
| AGNE | \n", "AGNE | \n", "1.5 | \n", "0.5 | \n", "
| AIGL | \n", "AIGL | \n", "0.7 | \n", "0.5 | \n", "
| ... | \n", "... | \n", "... | \n", "... | \n", "
| WLBH | \n", "WLBH | \n", "-2.8 | \n", "0.2 | \n", "
| WTZR | \n", "WTZR | \n", "-0.1 | \n", "0.1 | \n", "
| ZADA | \n", "ZADA | \n", "-0.3 | \n", "0.4 | \n", "
| ZIMM | \n", "ZIMM | \n", "1.0 | \n", "0.1 | \n", "
| ZOUF | \n", "ZOUF | \n", "1.3 | \n", "0.1 | \n", "
186 rows × 3 columns
\n", "| \n", " | station_id | \n", "velocity_east_mmyr | \n", "velocity_north_mmyr | \n", "velocity_east_error_mmyr | \n", "velocity_north_error_mmyr | \n", "velocity_up_mmyr | \n", "velocity_up_error_mmyr | \n", "latitude | \n", "longitude | \n", "height_m | \n", "latitude_error_m | \n", "longitude_error_m | \n", "height_error_m | \n", "
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", "ACOM | \n", "0.2 | \n", "1.2 | \n", "0.1 | \n", "0.1 | \n", "1.1 | \n", "0.1 | \n", "46.547935 | \n", "13.514900 | \n", "1774.682 | \n", "0.0009 | \n", "0.0005 | \n", "0.001 | \n", "
| 1 | \n", "AFAL | \n", "-0.7 | \n", "0.9 | \n", "0.1 | \n", "0.2 | \n", "1.3 | \n", "0.2 | \n", "46.527144 | \n", "12.174517 | \n", "2284.085 | \n", "0.0009 | \n", "0.0009 | \n", "0.001 | \n", "
| 2 | \n", "AGDE | \n", "-0.2 | \n", "-0.2 | \n", "0.1 | \n", "0.3 | \n", "0.1 | \n", "0.3 | \n", "43.296383 | \n", "3.466427 | \n", "65.785 | \n", "0.0018 | \n", "0.0009 | \n", "0.002 | \n", "
| 3 | \n", "AGNE | \n", "0.0 | \n", "-0.2 | \n", "0.2 | \n", "0.6 | \n", "1.5 | \n", "0.5 | \n", "45.467942 | \n", "7.139620 | \n", "2354.600 | \n", "0.0036 | \n", "0.0009 | \n", "0.004 | \n", "
| 4 | \n", "AIGL | \n", "0.0 | \n", "0.1 | \n", "0.1 | \n", "0.5 | \n", "0.7 | \n", "0.5 | \n", "44.121398 | \n", "3.581261 | \n", "1618.764 | \n", "0.0009 | \n", "0.0009 | \n", "0.002 | \n", "
| ... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
| 181 | \n", "WLBH | \n", "0.0 | \n", "-0.2 | \n", "0.1 | \n", "0.2 | \n", "-2.8 | \n", "0.2 | \n", "48.415171 | \n", "7.351299 | \n", "819.069 | \n", "0.0009 | \n", "0.0005 | \n", "0.001 | \n", "
| 182 | \n", "WTZR | \n", "0.1 | \n", "0.2 | \n", "0.1 | \n", "0.1 | \n", "-0.1 | \n", "0.1 | \n", "49.144199 | \n", "12.878911 | \n", "666.025 | \n", "0.0005 | \n", "0.0005 | \n", "0.001 | \n", "
| 183 | \n", "ZADA | \n", "0.2 | \n", "3.1 | \n", "0.2 | \n", "0.4 | \n", "-0.3 | \n", "0.4 | \n", "44.113177 | \n", "15.227590 | \n", "64.307 | \n", "0.0036 | \n", "0.0018 | \n", "0.004 | \n", "
| 184 | \n", "ZIMM | \n", "-0.1 | \n", "0.4 | \n", "0.1 | \n", "0.1 | \n", "1.0 | \n", "0.1 | \n", "46.877098 | \n", "7.465278 | \n", "956.341 | \n", "0.0009 | \n", "0.0005 | \n", "0.001 | \n", "
| 185 | \n", "ZOUF | \n", "0.1 | \n", "1.0 | \n", "0.1 | \n", "0.1 | \n", "1.3 | \n", "0.1 | \n", "46.557221 | \n", "12.973553 | \n", "1946.508 | \n", "0.0009 | \n", "0.0005 | \n", "0.001 | \n", "
186 rows × 13 columns
\n", "| \n", " | station_id | \n", "longitude | \n", "latitude | \n", "height_m | \n", "velocity_east_mmyr | \n", "velocity_north_mmyr | \n", "velocity_up_mmyr | \n", "longitude_error_m | \n", "latitude_error_m | \n", "height_error_m | \n", "velocity_east_error_mmyr | \n", "velocity_north_error_mmyr | \n", "velocity_up_error_mmyr | \n", "
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", "ACOM | \n", "13.5149004 | \n", "46.5479352 | \n", "1774.682 | \n", "0.2 | \n", "1.2 | \n", "1.1 | \n", "0.0005 | \n", "0.0009 | \n", "0.001 | \n", "0.1 | \n", "0.1 | \n", "0.1 | \n", "
| 1 | \n", "AFAL | \n", "12.1745170 | \n", "46.5271440 | \n", "2284.085 | \n", "-0.7 | \n", "0.9 | \n", "1.3 | \n", "0.0009 | \n", "0.0009 | \n", "0.001 | \n", "0.1 | \n", "0.2 | \n", "0.2 | \n", "
| 2 | \n", "AGDE | \n", "3.4664272 | \n", "43.2963831 | \n", "65.785 | \n", "-0.2 | \n", "-0.2 | \n", "0.1 | \n", "0.0009 | \n", "0.0018 | \n", "0.002 | \n", "0.1 | \n", "0.3 | \n", "0.3 | \n", "
| 3 | \n", "AGNE | \n", "7.1396203 | \n", "45.4679419 | \n", "2354.600 | \n", "0.0 | \n", "-0.2 | \n", "1.5 | \n", "0.0009 | \n", "0.0036 | \n", "0.004 | \n", "0.2 | \n", "0.6 | \n", "0.5 | \n", "
| 4 | \n", "AIGL | \n", "3.5812609 | \n", "44.1213985 | \n", "1618.764 | \n", "0.0 | \n", "0.1 | \n", "0.7 | \n", "0.0009 | \n", "0.0009 | \n", "0.002 | \n", "0.1 | \n", "0.5 | \n", "0.5 | \n", "
| ... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
| 181 | \n", "WLBH | \n", "7.3512992 | \n", "48.4151715 | \n", "819.069 | \n", "0.0 | \n", "-0.2 | \n", "-2.8 | \n", "0.0005 | \n", "0.0009 | \n", "0.001 | \n", "0.1 | \n", "0.2 | \n", "0.2 | \n", "
| 182 | \n", "WTZR | \n", "12.8789111 | \n", "49.1441992 | \n", "666.025 | \n", "0.1 | \n", "0.2 | \n", "-0.1 | \n", "0.0005 | \n", "0.0005 | \n", "0.001 | \n", "0.1 | \n", "0.1 | \n", "0.1 | \n", "
| 183 | \n", "ZADA | \n", "15.2275901 | \n", "44.1131766 | \n", "64.307 | \n", "0.2 | \n", "3.1 | \n", "-0.3 | \n", "0.0018 | \n", "0.0036 | \n", "0.004 | \n", "0.2 | \n", "0.4 | \n", "0.4 | \n", "
| 184 | \n", "ZIMM | \n", "7.4652781 | \n", "46.8770983 | \n", "956.341 | \n", "-0.1 | \n", "0.4 | \n", "1.0 | \n", "0.0005 | \n", "0.0009 | \n", "0.001 | \n", "0.1 | \n", "0.1 | \n", "0.1 | \n", "
| 185 | \n", "ZOUF | \n", "12.9735532 | \n", "46.5572211 | \n", "1946.508 | \n", "0.1 | \n", "1.0 | \n", "1.3 | \n", "0.0005 | \n", "0.0009 | \n", "0.001 | \n", "0.1 | \n", "0.1 | \n", "0.1 | \n", "
186 rows × 13 columns
\n", "| \n", " | station_id | \n", "longitude | \n", "latitude | \n", "height_m | \n", "velocity_east_mmyr | \n", "velocity_north_mmyr | \n", "velocity_up_mmyr | \n", "longitude_error_m | \n", "latitude_error_m | \n", "height_error_m | \n", "velocity_east_error_mmyr | \n", "velocity_north_error_mmyr | \n", "velocity_up_error_mmyr | \n", "
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", "ACOM | \n", "13.514900 | \n", "46.547935 | \n", "1774.682 | \n", "0.2 | \n", "1.2 | \n", "1.1 | \n", "0.0005 | \n", "0.0009 | \n", "0.001 | \n", "0.1 | \n", "0.1 | \n", "0.1 | \n", "
| 1 | \n", "AFAL | \n", "12.174517 | \n", "46.527144 | \n", "2284.085 | \n", "-0.7 | \n", "0.9 | \n", "1.3 | \n", "0.0009 | \n", "0.0009 | \n", "0.001 | \n", "0.1 | \n", "0.2 | \n", "0.2 | \n", "
| 2 | \n", "AGDE | \n", "3.466427 | \n", "43.296383 | \n", "65.785 | \n", "-0.2 | \n", "-0.2 | \n", "0.1 | \n", "0.0009 | \n", "0.0018 | \n", "0.002 | \n", "0.1 | \n", "0.3 | \n", "0.3 | \n", "
| 3 | \n", "AGNE | \n", "7.139620 | \n", "45.467942 | \n", "2354.600 | \n", "0.0 | \n", "-0.2 | \n", "1.5 | \n", "0.0009 | \n", "0.0036 | \n", "0.004 | \n", "0.2 | \n", "0.6 | \n", "0.5 | \n", "
| 4 | \n", "AIGL | \n", "3.581261 | \n", "44.121398 | \n", "1618.764 | \n", "0.0 | \n", "0.1 | \n", "0.7 | \n", "0.0009 | \n", "0.0009 | \n", "0.002 | \n", "0.1 | \n", "0.5 | \n", "0.5 | \n", "
| ... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
| 181 | \n", "WLBH | \n", "7.351299 | \n", "48.415171 | \n", "819.069 | \n", "0.0 | \n", "-0.2 | \n", "-2.8 | \n", "0.0005 | \n", "0.0009 | \n", "0.001 | \n", "0.1 | \n", "0.2 | \n", "0.2 | \n", "
| 182 | \n", "WTZR | \n", "12.878911 | \n", "49.144199 | \n", "666.025 | \n", "0.1 | \n", "0.2 | \n", "-0.1 | \n", "0.0005 | \n", "0.0005 | \n", "0.001 | \n", "0.1 | \n", "0.1 | \n", "0.1 | \n", "
| 183 | \n", "ZADA | \n", "15.227590 | \n", "44.113177 | \n", "64.307 | \n", "0.2 | \n", "3.1 | \n", "-0.3 | \n", "0.0018 | \n", "0.0036 | \n", "0.004 | \n", "0.2 | \n", "0.4 | \n", "0.4 | \n", "
| 184 | \n", "ZIMM | \n", "7.465278 | \n", "46.877098 | \n", "956.341 | \n", "-0.1 | \n", "0.4 | \n", "1.0 | \n", "0.0005 | \n", "0.0009 | \n", "0.001 | \n", "0.1 | \n", "0.1 | \n", "0.1 | \n", "
| 185 | \n", "ZOUF | \n", "12.973553 | \n", "46.557221 | \n", "1946.508 | \n", "0.1 | \n", "1.0 | \n", "1.3 | \n", "0.0005 | \n", "0.0009 | \n", "0.001 | \n", "0.1 | \n", "0.1 | \n", "0.1 | \n", "
186 rows × 13 columns
\n", "