🌟 Dropping code!

Add all code for producing figures for the manuscript:

Lithospheric anisotropy characterises the post-subduction setting of
northern Borneo: new results from XKS splitting analysis

LICENSE

@@ -0,0 +1,121 @@
+Creative Commons Legal Code
+
+CC0 1.0 Universal
+
+    CREATIVE COMMONS CORPORATION IS NOT A LAW FIRM AND DOES NOT PROVIDE
+    LEGAL SERVICES. DISTRIBUTION OF THIS DOCUMENT DOES NOT CREATE AN
+    ATTORNEY-CLIENT RELATIONSHIP. CREATIVE COMMONS PROVIDES THIS
+    INFORMATION ON AN "AS-IS" BASIS. CREATIVE COMMONS MAKES NO WARRANTIES
+    REGARDING THE USE OF THIS DOCUMENT OR THE INFORMATION OR WORKS
+    PROVIDED HEREUNDER, AND DISCLAIMS LIABILITY FOR DAMAGES RESULTING FROM
+    THE USE OF THIS DOCUMENT OR THE INFORMATION OR WORKS PROVIDED
+    HEREUNDER.
+
+Statement of Purpose
+
+The laws of most jurisdictions throughout the world automatically confer
+exclusive Copyright and Related Rights (defined below) upon the creator
+and subsequent owner(s) (each and all, an "owner") of an original work of
+authorship and/or a database (each, a "Work").
+
+Certain owners wish to permanently relinquish those rights to a Work for
+the purpose of contributing to a commons of creative, cultural and
+scientific works ("Commons") that the public can reliably and without fear
+of later claims of infringement build upon, modify, incorporate in other
+works, reuse and redistribute as freely as possible in any form whatsoever
+and for any purposes, including without limitation commercial purposes.
+These owners may contribute to the Commons to promote the ideal of a free
+culture and the further production of creative, cultural and scientific
+works, or to gain reputation or greater distribution for their Work in
+part through the use and efforts of others.
+
+For these and/or other purposes and motivations, and without any
+expectation of additional consideration or compensation, the person
+associating CC0 with a Work (the "Affirmer"), to the extent that he or she
+is an owner of Copyright and Related Rights in the Work, voluntarily
+elects to apply CC0 to the Work and publicly distribute the Work under its
+terms, with knowledge of his or her Copyright and Related Rights in the
+Work and the meaning and intended legal effect of CC0 on those rights.
+
+1. Copyright and Related Rights. A Work made available under CC0 may be
+protected by copyright and related or neighboring rights ("Copyright and
+Related Rights"). Copyright and Related Rights include, but are not
+limited to, the following:
+
+  i. the right to reproduce, adapt, distribute, perform, display,
+     communicate, and translate a Work;
+ ii. moral rights retained by the original author(s) and/or performer(s);
+iii. publicity and privacy rights pertaining to a person's image or
+     likeness depicted in a Work;
+ iv. rights protecting against unfair competition in regards to a Work,
+     subject to the limitations in paragraph 4(a), below;
+  v. rights protecting the extraction, dissemination, use and reuse of data
+     in a Work;
+ vi. database rights (such as those arising under Directive 96/9/EC of the
+     European Parliament and of the Council of 11 March 1996 on the legal
+     protection of databases, and under any national implementation
+     thereof, including any amended or successor version of such
+     directive); and
+vii. other similar, equivalent or corresponding rights throughout the
+     world based on applicable law or treaty, and any national
+     implementations thereof.
+
+2. Waiver. To the greatest extent permitted by, but not in contravention
+of, applicable law, Affirmer hereby overtly, fully, permanently,
+irrevocably and unconditionally waives, abandons, and surrenders all of
+Affirmer's Copyright and Related Rights and associated claims and causes
+of action, whether now known or unknown (including existing as well as
+future claims and causes of action), in the Work (i) in all territories
+worldwide, (ii) for the maximum duration provided by applicable law or
+treaty (including future time extensions), (iii) in any current or future
+medium and for any number of copies, and (iv) for any purpose whatsoever,
+including without limitation commercial, advertising or promotional
+purposes (the "Waiver"). Affirmer makes the Waiver for the benefit of each
+member of the public at large and to the detriment of Affirmer's heirs and
+successors, fully intending that such Waiver shall not be subject to
+revocation, rescission, cancellation, termination, or any other legal or
+equitable action to disrupt the quiet enjoyment of the Work by the public
+as contemplated by Affirmer's express Statement of Purpose.
+
+3. Public License Fallback. Should any part of the Waiver for any reason
+be judged legally invalid or ineffective under applicable law, then the
+Waiver shall be preserved to the maximum extent permitted taking into
+account Affirmer's express Statement of Purpose. In addition, to the
+extent the Waiver is so judged Affirmer hereby grants to each affected
+person a royalty-free, non transferable, non sublicensable, non exclusive,
+irrevocable and unconditional license to exercise Affirmer's Copyright and
+Related Rights in the Work (i) in all territories worldwide, (ii) for the
+maximum duration provided by applicable law or treaty (including future
+time extensions), (iii) in any current or future medium and for any number
+of copies, and (iv) for any purpose whatsoever, including without
+limitation commercial, advertising or promotional purposes (the
+"License"). The License shall be deemed effective as of the date CC0 was
+applied by Affirmer to the Work. Should any part of the License for any
+reason be judged legally invalid or ineffective under applicable law, such
+partial invalidity or ineffectiveness shall not invalidate the remainder
+of the License, and in such case Affirmer hereby affirms that he or she
+will not (i) exercise any of his or her remaining Copyright and Related
+Rights in the Work or (ii) assert any associated claims and causes of
+action with respect to the Work, in either case contrary to Affirmer's
+express Statement of Purpose.
+
+4. Limitations and Disclaimers.
+
+ a. No trademark or patent rights held by Affirmer are waived, abandoned,
+    surrendered, licensed or otherwise affected by this document.
+ b. Affirmer offers the Work as-is and makes no representations or
+    warranties of any kind concerning the Work, express, implied,
+    statutory or otherwise, including without limitation warranties of
+    title, merchantability, fitness for a particular purpose, non
+    infringement, or the absence of latent or other defects, accuracy, or
+    the present or absence of errors, whether or not discoverable, all to
+    the greatest extent permissible under applicable law.
+ c. Affirmer disclaims responsibility for clearing rights of other persons
+    that may apply to the Work or any use thereof, including without
+    limitation any person's Copyright and Related Rights in the Work.
+    Further, Affirmer disclaims responsibility for obtaining any necessary
+    consents, permissions or other rights required for any use of the
+    Work.
+ d. Affirmer understands and acknowledges that Creative Commons is not a
+    party to this document and has no duty or obligation with respect to
+    this CC0 or use of the Work.

README.md

@@ -0,0 +1,31 @@
+# Lithospheric anisotropy characterises the post-subduction setting of northern Borneo: new results from XKS splitting analysis
+Data and code required to reproduce the results presented in the manuscript:
+
+Lithospheric anisotropy characterises the post-subduction setting of northern Borneo: new results from XKS splitting analysis
+
+Bacon, C.A., Rawlinson, N., Pilia, S., Gilligan, A., Wehner, D., Cornwell, D.G., Tongkul, F.
+
+## Steps to reproduce figures
+1. Install the packages listed in the environment.yml file, either manually, or using (for example) conda:
+
+```
+conda install -f environment.yml
+```
+
+2. Add `splitracer_summary.mplstyle` (a matplotlib stylesheet) to your `mpl_configdir` (usually found at `~/.config/matplotlib`)
+
+3. Optional: Install Helvetica font for Matplotlib
+
+4. Run `download_grd_datafiles.gmt` script to download the DEM data used in Figure 1:
+
+```
+cd figures
+bash download_grd_datafiles.gmt
+```
+
+5. Navigate to each figure directory and run the `.gmt` or `.py` scripts.
+
+## Notes
+These figures were prepared using Linux 20.04. A limited number were produced using Affinity Designer for iPad, a licensed piece of software for graphic design. `.afdesign` files are provided.
+
+For Supplementary Figure S6, one must install `AnisotroPy`, an open-source package for the analysis of seismic anisotropy. This can be downloaded from GitHub at https://github.com/hemmelig/AnisotroPy, where one can also find instructions on how to install the package.

environment.yml

@@ -0,0 +1,13 @@
+name: epsl_manuscript_CABacon
+channels:
+  - conda-forge
+  - defaults
+dependencies:
+  - matplotlib=3.5.1
+  - numpy=1.21.5
+  - obspy=1.2.2
+  - pandas=1.4.1
+  - pygmt=0.6.0
+  - pyproj=3.3.0
+  - python=3.8.13
+  - scipy=1.8.0

figures/README

@@ -0,0 +1,21 @@
+Main text
+---------
+Figure 1 - Network and context map showing Sabah and the wider northern Borneo region.
+Figure 2 - Cartoon showing tectonic events in the last 20 Ma.
+Figure 3 - Joint analysis of splitting measurements at MetMalaysia station KKM.
+Figure 4 - Map of splitting results, overlaid on lithospheric thickness estimates from two-plane-wave tomography.
+Figure 5 - Interpretation cartoon linking observed splitting trends with tectonic history of northern Borneo.
+Figure 6 - Radial anisotropy derived from FWI tomographic model, SASSY21.
+
+Supplementary material
+----------------------
+Figure S1 - Simplified geological map of northern Borneo.
+Figure S2 - Annotated network map.
+Figure S3 - Maps showing earthquake catalogues and azimuthal coverage.
+Figure S4 - An example of a good splitting measurement.
+Figure S5 - An example of a null splitting measurement.
+Figure S6 - n-layer modelling for KKM.
+Figure S7 - Map of splitting results, overlaid on crustal thickness estimates from VDSS.
+Figure S8 - Map of splitting results, overlaid on crustal thickness estimates from receiver function analysis.
+Figure S9 - Map of splitting results, overlaid on horizontal 100-km depth slice through radial anisotropy model derived from SASSY21.
+Figure S10 - Various depth slices through radial anisotropy model derived from SASSY21.

figures/data/earthquake_catalogues/supp_file_3_2018-2020_earthquakes.csv

@@ -0,0 +1,129 @@
+2020-01-20T06:51:37.820000Z -59.229 -26.289 93.0 6.1 Mww
+2020-01-17T06:37:52.714000Z -11.637 -13.782 10.0 5.8 Mww
+2020-01-11T12:54:45.910000Z 17.949 -66.851 5.0 5.9 Mwp
+2020-01-07T11:18:43.540000Z 18.022 -66.776 9.0 5.8 Mww
+2020-01-07T08:24:25.143000Z 17.867 -66.828 8.07 6.4 mww
+2020-01-06T10:32:18.770000Z 17.867 -66.819 6.0 5.8 mwp
+2020-01-05T04:40:49.807000Z 16.450 -94.360 87.0 5.9 Mww
+2019-12-25T20:20:12.751000Z -25.223 -115.844 10.0 6.1 Mww
+2019-12-25T03:36:01.578000Z 50.608 -129.966 6.38 6.3 mww
+2019-12-24T19:03:52.514000Z 3.527 -74.050 11.0 6.0 Mww
+2019-12-24T16:43:33.782000Z -26.991 -63.375 581.63 6.0 Mww
+2019-12-23T20:56:23.555000Z 50.517 -129.757 10.0 6.0 mww
+2019-12-23T19:49:43.101000Z 50.608 -129.927 10.0 6.0 mww
+2019-12-23T15:47:56.102000Z 13.764 -92.861 10.0 5.8 Mww
+2019-12-19T12:35:32.201000Z 13.690 -91.314 31.22 5.9 Mww
+2019-12-03T08:46:36.374000Z -18.560 -70.650 32.44 6.0 Mww
+2019-11-27T07:23:42.552000Z 35.727 23.267 71.76 6.0 Mww
+2019-11-26T02:54:12.620000Z 41.515 19.521 20.0 6.4 Mww
+2019-11-20T04:27:05.120000Z 13.886 -93.213 14.16 6.3 Mww
+2019-11-05T20:52:01.296000Z -57.964 -9.270 10.0 6.3 Mww
+2019-11-04T21:53:25.320000Z -31.836 -71.375 53.0 6.1 Mww
+2019-11-02T18:08:41.079000Z -55.666 -26.234 8.0 6.1 Mww
+2019-10-23T16:08:14.497000Z -1.099 -13.837 10.0 6.3 mww
+2019-10-01T03:56:07.917000Z -56.005 -122.470 10.0 5.9 Mww
+2019-09-29T15:57:53.236000Z -35.476 -73.163 11.0 6.7 mww
+2019-09-27T15:11:29.878000Z -35.922 -102.774 10.0 5.9 Mww
+2019-09-26T16:36:18.134000Z -40.815 -72.002 129.0 6.1 Mww
+2019-09-24T03:23:40.000000Z 19.077 -67.270 10.0 6.0 mw
+2019-09-05T15:02:45.260000Z 43.713 -127.827 10.0 5.9 Mww
+2019-09-02T22:45:48.807000Z 23.777 -45.217 10.0 5.9 Mww
+2019-08-29T15:07:57.794000Z 43.544 -127.876 4.74 6.3 Mww
+2019-08-27T23:55:19.187000Z -60.215 -26.580 16.0 6.6 mww
+2019-08-14T21:35:18.158000Z 20.460 -109.361 10.0 5.9 Mww
+2019-08-05T00:40:46.037000Z 1.039 -27.875 10.0 5.9 Mww
+2019-08-02T05:50:55.239000Z -49.731 -113.833 10.0 6.0 Mww
+2019-08-01T18:28:07.272000Z -34.237 -72.308 25.0 6.8 mww
+2019-07-31T05:54:55.289000Z 13.266 -89.338 72.5 5.9 Mww
+2019-07-13T07:59:33.543000Z -49.551 -116.744 10.0 5.8 Mww
+2019-07-06T03:19:53.040000Z 35.770 -117.599 8.0 7.1 mw
+2019-07-04T17:33:49.040000Z 35.705 -117.506 10.71 6.4 mw
+2019-07-04T04:30:44.297000Z 51.221 -130.515 10.0 6.2 Mww
+2019-06-26T05:23:50.890000Z 8.449 -82.768 26.22 6.2 Mww
+2019-06-17T13:47:48.112000Z -65.607 -176.702 10.0 5.8 Mww
+2019-06-14T00:19:12.401000Z -30.056 -72.082 11.0 6.4 Mww
+2019-05-31T11:57:24.579000Z 18.732 -107.169 10.0 5.8 Mww
+2019-05-30T09:03:28.971000Z 13.146 -89.366 25.0 6.6 Mww
+2019-05-26T07:41:15.058000Z -5.813 -75.278 122.4 8.0 Mww
+2019-05-16T22:52:43.938000Z -4.588 -105.731 10.0 5.8 Mww
+2019-05-16T16:22:16.602000Z 12.629 -87.843 62.86 5.9 Mww
+2019-05-12T19:24:50.395000Z 8.623 -82.833 19.0 6.0 Mww
+2019-04-22T14:49:05.376000Z -56.263 -27.351 79.0 6.0 Mww
+2019-04-09T17:53:59.096000Z -58.606 -25.256 44.83 6.5 Mww
+2019-04-05T18:46:42.673000Z 1.451 -85.266 10.0 5.8 Mww
+2019-04-05T16:14:16.754000Z -55.921 -27.856 58.6 6.4 mww
+2019-03-31T07:04:04.800000Z -1.944 -80.809 18.0 6.2 Mww
+2019-03-23T19:21:18.010000Z 4.563 -76.223 122.0 6.1 Mww
+2019-03-15T05:03:50.060000Z -17.874 -65.907 359.0 6.3 Mww
+2019-03-10T18:04:59.100000Z 58.351 -31.950 10.0 5.9 Mww
+2019-03-01T08:50:42.620000Z -14.702 -70.135 267.0 7.0 mww
+2019-02-22T10:17:22.410000Z -2.199 -77.023 132.36 7.5 mww
+2019-02-14T19:57:05.020000Z 35.422 -36.079 10.0 6.2 Mww
+2019-02-04T17:45:48.120000Z -2.405 -79.835 70.0 5.8 Mww
+2019-02-01T16:14:13.090000Z 14.763 -92.298 67.93 6.6 Mww
+2019-01-30T15:31:33.440000Z -4.633 -105.480 10.0 5.9 Mww
+2019-01-20T01:32:51.850000Z -30.071 -71.420 54.82 6.7 mww
+2019-01-18T16:40:43.360000Z 8.401 -103.248 10.0 5.9 Mww
+2019-01-05T19:25:38.580000Z -8.149 -71.583 571.18 6.8 Mww
+2018-12-19T01:37:40.500000Z -36.138 -101.072 10.0 6.3 Mww
+2018-12-12T13:13:56.150000Z -55.676 -128.684 10.0 6.3 Mww
+2018-12-11T02:26:32.730000Z -58.598 -26.466 164.66 7.1 mww
+2018-11-30T17:29:29.330000Z 61.346 -149.955 46.7 7.0 mww
+2018-11-27T12:03:28.000000Z -1.029 -13.663 10.0 5.8 Mww
+2018-11-25T03:40:50.810000Z 13.182 -81.093 10.0 6.0 Mww
+2018-11-22T16:07:05.380000Z -54.205 -137.499 10.0 5.9 Mww
+2018-11-15T23:09:01.060000Z -56.236 -122.044 10.0 6.3 Mww
+2018-11-15T20:02:22.920000Z -56.706 -25.546 15.0 6.4 mww
+2018-11-11T14:03:59.560000Z 15.565 -49.872 10.0 6.3 Mww
+2018-11-09T01:49:40.060000Z 71.626 -11.250 10.0 6.7 Mww
+2018-11-01T22:19:51.540000Z -19.588 -69.292 102.0 6.2 Mww
+2018-11-01T19:30:20.940000Z -58.074 -25.206 29.0 5.8 Mww
+2018-10-29T23:26:09.160000Z -4.573 -105.910 10.0 5.8 Mww
+2018-10-29T20:17:22.520000Z -57.550 -66.304 10.0 5.8 Mww
+2018-10-29T06:54:21.440000Z -57.404 -66.409 10.0 6.3 Mww
+2018-10-28T22:23:54.080000Z 12.949 -90.385 24.69 6.1 Mww
+2018-10-25T22:54:52.590000Z 37.515 20.564 14.0 6.8 Mww
+2018-10-22T06:22:48.810000Z 49.314 -129.673 10.0 6.5 Mww
+2018-10-22T06:16:28.160000Z 49.347 -129.214 10.0 6.8 mww
+2018-10-22T05:39:40.330000Z 49.250 -129.478 10.96 6.6 mww
+2018-10-07T00:11:51.220000Z 20.034 -73.013 24.0 5.9 Mww
+2018-10-07T00:10:04.880000Z 7.661 -37.708 10.0 5.8 Mww
+2018-10-05T14:59:44.140000Z -5.691 -11.525 10.0 5.8 Mww
+2018-09-07T02:39:17.730000Z -28.892 -69.954 92.4 5.8 Mww
+2018-09-07T02:12:04.350000Z -2.345 -78.895 93.53 6.2 Mww
+2018-09-03T23:37:17.590000Z 1.439 -85.261 10.0 5.9 Mww
+2018-08-24T09:04:06.780000Z -11.042 -70.817 609.48 7.1 Mww
+2018-08-22T13:27:48.490000Z 10.659 -62.929 108.0 5.8 Mww
+2018-08-22T09:31:47.050000Z 43.645 -127.603 10.0 6.2 Mww
+2018-08-21T21:31:47.540000Z 10.779 -62.907 146.18 7.3 Mww
+2018-08-17T23:22:24.900000Z 8.769 -83.153 15.0 6.1 Mww
+2018-08-14T03:29:53.440000Z -58.111 -25.259 35.0 6.1 Mww
+2018-08-12T21:15:01.841000Z 69.520 -144.360 1.7 6.0 Mww
+2018-08-12T14:58:54.286000Z 69.562 -145.300 2.2 6.3 ml
+2018-08-03T18:50:53.990000Z -0.915 -21.978 10.0 5.8 Mww
+2018-07-23T10:36:00.330000Z -0.350 -19.234 10.0 5.9 Mww
+2018-07-19T13:31:53.880000Z 17.948 -97.736 49.75 5.8 Mww
+2018-06-30T03:56:49.400000Z 19.016 -105.117 15.0 5.9 Mww
+2018-06-24T18:43:43.340000Z -9.019 -108.467 10.0 5.9 Mww
+2018-06-18T07:09:41.680000Z -29.929 -111.877 10.0 5.9 Mww
+2018-05-04T22:32:54.370000Z 19.313 -154.998 2.06 6.9 Mw
+2018-05-02T06:32:48.880000Z -24.260 -111.654 10.0 6.0 Mww
+2018-05-01T17:15:37.770000Z -56.640 -142.011 10.0 5.9 Mww
+2018-04-10T10:19:34.030000Z -30.986 -71.557 76.08 6.2 Mww
+2018-04-02T23:23:17.670000Z 13.107 -88.687 50.0 5.9 Mww
+2018-04-02T13:40:34.840000Z -20.659 -63.006 559.0 6.8 Mww
+2018-03-10T21:45:36.340000Z -1.423 -15.210 10.0 5.9 Mww
+2018-02-19T06:57:01.830000Z 16.499 -97.579 40.17 5.9 Mww
+2018-02-17T00:36:52.230000Z 16.158 -97.882 4.26 5.8 Mww
+2018-02-16T23:39:39.700000Z 16.389 -97.979 25.98 7.2 mww
+2018-02-09T14:05:48.170000Z 18.934 -104.649 12.0 5.8 Mww
+2018-02-02T11:37:53.320000Z -65.811 -175.635 10.0 6.0 Mww
+2018-01-28T16:03:03.970000Z -53.062 9.684 10.0 6.6 Mww
+2018-01-23T09:31:42.940000Z 56.046 -149.073 25.0 7.9 Mww
+2018-01-21T01:06:42.890000Z -18.852 -69.456 116.0 6.3 Mww
+2018-01-19T16:17:42.570000Z 26.680 -111.107 10.0 6.3 Mww
+2018-01-16T15:03:53.300000Z 10.929 -86.347 22.0 5.8 Mww
+2018-01-14T09:18:45.540000Z -15.759 -74.710 39.0 7.1 mww
+2018-01-10T02:51:31.860000Z 17.474 -83.519 10.0 7.5 mww
+2018-01-09T00:05:54.210000Z 1.675 -90.822 10.0 5.8 Mww
+2019-08-08T11:25:31.024000Z 37.948 29.697 10.0 5.8 Mww