diff --git a/Fig3_PyGMT_backgrounds.ipynb b/Fig3_PyGMT_backgrounds.ipynb
new file mode 100644
index 0000000..abcc0a2
--- /dev/null
+++ b/Fig3_PyGMT_backgrounds.ipynb
@@ -0,0 +1,76 @@
+{
+ "cells": [
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "d57cb626-f006-4a23-950d-c673fd04b840",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import pygmt\n",
+    "\n",
+    "region = [-25, -13, 63.2, 66.7]\n",
+    "\n",
+    "fig = pygmt.Figure()\n",
+    "with fig.subplot(\n",
+    "    nrows=2,\n",
+    "    ncols=2,\n",
+    "    subsize=(\"12c\", \"8.5c\"),\n",
+    "    margins=(\"0.3c\", \"0.2c\"),\n",
+    "    frame=\"WSen\",\n",
+    "    sharex=\"b\",\n",
+    "    sharey=\"l\",\n",
+    "    autolabel=\"(a)+o0.15c/0.3c+gwhite@30+v\",\n",
+    "):\n",
+    "    # Top left\n",
+    "    fig.basemap(region=region, projection=\"M?\", panel=0)\n",
+    "    fig.coast(shorelines=True, land=\"gray\", water=\"lightblue\", resolution=\"high\")\n",
+    "    fig.basemap(map_scale=\"n0.86/0.1+c+w100k+f+l\")\n",
+    "\n",
+    "    # Bottom left\n",
+    "    fig.tilemap(region=region, projection=\"M?\", zoom=7, panel=1)\n",
+    "\n",
+    "    # Top right\n",
+    "    grd_relief = pygmt.datasets.load_earth_relief(resolution=\"01m\", region=region)\n",
+    "    fig.basemap(region=region, projection=\"M?\", panel=2)\n",
+    "    fig.grdimage(grid=grd_relief, cmap=\"SCM/oleron\", shading=True)\n",
+    "    fig.colorbar(frame=[\"xa500f250+lElevation\", \"y+lm\"], position=\"JMR\", shading=True)\n",
+    "\n",
+    "    # Bottom right\n",
+    "    fig.basemap(region=region, projection=\"M?\", perspective=(-150, 25), panel=3)\n",
+    "    fig.grdview(\n",
+    "        grid=grd_relief,\n",
+    "        cmap=\"SCM/oleron\",\n",
+    "        surftype=\"s\",\n",
+    "        shading=True,\n",
+    "        zsize=\"1.5c\",\n",
+    "        plane=\"+ggray\",\n",
+    "        perspective=True,\n",
+    "    )\n",
+    "fig.show()\n",
+    "# fig.savefig(\"Fig3_PyGMT_backgrounds.png\")"
+   ]
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "Python 3 (ipykernel)",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.12.3"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}