# Estimation of Lesser Antilles Vertical Velocity Fields Using a GNSS-PPP Software Comparison

Pierre Sakic, Benjamin Männel, Markus Bradke, Valérie Ballu, Jean-Bernard de Chabalier, Arnaud Lemarchand

► To cite this version:

Pierre Sakic, Benjamin Männel, Markus Bradke, Valérie Ballu, Jean-Bernard de Chabalier, et al. Estimation of Lesser Antilles Vertical Velocity Fields Using a GNSS-PPP Software Comparison. International Association of Geodesy Symposia, 2020, 10.1007/1345_2020_101. hal-03004711

HAL Id: hal-03004711

https://hal.science/hal-03004711

Submitted on 13 Nov 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Estimation of Lesser Antilles Vertical Velocity Fields Using a GNSS-PPP Software Comparison

Pierre Sakic, Benjamin Männel, Markus Bradke, Valérie Ballu, Jean-Bernard de Chabalier, and Arnaud Lemarchand

## Abstract

Vertical land motion in insular areas is a crucial parameter to estimate the relative sea-level variations which impact coastal populations and activities. In subduction zones, it is also a relevant proxy to estimate the locking state of the plate interface. This motion can be measured using Global Navigation Satellite Systems (GNSS), such as the Global Positioning System (GPS). However, the influence of the processing software and the geodetic products (orbits and clock offsets) used for the solution remains barely considered for geophysics studies.

In this study, we process GNSS observations of Guadeloupe and Martinique network (Lesser Antilles). It consists of 40 stations over a period of 18 years for the oldest site. We provide an updated vertical velocity field determined with two different geodetic software, namely EPOS (Gendt et al, GFZ analysis center of IGS–Annual Report. IGS 1996 Annual Report, pp 169–181, 1998) and GINS (Marty et al, GINS: the CNES/GRGS GNSS scientific software. In: 3rd International colloquium scientific and fundamental aspects of the Galileo programme, ESA proceedings WPP326, vol 31, pp 8–10, 2011) using their Precise Point Positioning modes. We used the same input models and orbit and clock offset products to maintain a maximum of consistency, and then compared the obtained results to get an estimation of the time series accuracy and the software influence on the solutions. General consistency between the solutions is noted, but significant velocity differences exist (at the mm/yr level) for some stations.

## Keywords

Lesser Antilles · Guadeloupe · Martinique · GNSS · Vertical velocity field · Subsidence · Processing comparison · Time series analysis

## 1 Introduction

Nowadays, the Global Navigation Satellite Systems (GNSS) have become an indispensable tool to monitor the Earth crust motion. Nevertheless, their use in some part of the world remains challenging. It is the case of the Lesser Antilles Subduction Zone, at the convergence of the Nord-American Plate under the Caribbean Plate. This subduction is singular on several aspects: it is one of the slowest in the world (~2cm/yr), the lack of emerged lands prevents the determination of a complete deformation profile like it

P. Sakic (✉) · B. Männel · M. Bradke
GFZ German Research Centre for Geosciences, Helmholtz-Zentrum
Potsdam, Potsdam, Germany
e-mail: pierre.sakic@gfz-potsdam.de

V. Ballu
Littoral Environnement et Sociétés, CNRS & University of La
Rochelle, La Rochelle, France

J.-B. de Chabalier · A. Lemarchand
Institut de Physique du Globe de Paris, Université de Paris, Paris,
France
can be done in other areas. Moreover, the islands of the volcanic arc are located too far from the trench, prohibiting most often the detection of significant velocity gradients with respect to the stable plate. Because of these reasons, the locking state deduced from GNSS observations, and thus megathrust risk is uncertain. Symithe et al. (2015) estimated with GNSS observations all along the volcanic arc a low coupling rate but did not exclude a megathrust possibility either. Since the estimation of a horizontal deformation rate is a difficult task for this area, vertical motion observations can then become a proxy to help the assessment of a potential strain accumulation (e.g. Mouslopoulou et al. 2016). Moreover, island areas are also threatened by the sea level rise, and extra subsidence can be an aggravating factor (Ballu et al. 2011). For these two reasons, measuring vertical motion in the Lesser Antilles is crucial. Some vertical movement assessments in this area were performed in the past. Paleo-geodesy based on coral reef growth tends to show a subsidence trend in Martinique and Les Saintes Islands (south of Guadeloupe archipelago) (Weil-Accardo et al. 2016; Leclerc and Feuillet 2019). This subsidence is corroborated by GNSS observations for a few stations within the vertical velocity ULR6 solution (Santamaría-Gómez et al. 2017). However, an uplift with decreasing rate for la Désirade Island (West of the Guadeloupe Archipelago) was measured (Léticée et al. 2019). Even though GNSS exploitation in this area is challenging, the islands are generally well instrumented, especially the two islands of Guadeloupe and Martinique.

Therefore, this study has a double objective, technical and scientific: on one hand, we estimate a denser vertical velocity field using a maximum of GNSS data in the area. On the other hand, we compare and quantify the differences between the results (coordinates and velocities inferred from the time series) obtained with two GNSS processing software using homogeneous inputs.

## 2 Data

We considered the GNSS observations provided by three different station networks, namely the IPGP, the IGN/SONEL, and the ORPHEON networks. Maps of these networks on the two islands are represented in Fig. 1. All three have been deployed for different purposes. The IPGP (Institut de Physique du Globe de Paris) network, maintained by the two local volcanological and seismological observatories has been deployed since the early 2000s (the first stations, HOUE and SOUF, were installed in 2000) for geodynamic purposes. Some stations are installed in the vicinity of the volcanic domes to monitor the volcanic activity and the others are deployed to observe potential subduction induced deformations. This network gathers nowadays 28

stations, 20 in Guadeloupe and 8 in Martinique. Due to the remote and/or extreme conditions for some sites, the network is heterogeneous regarding the time series completeness and the continuity in time of the equipment used. The IGN/SONEL (Institut national de l'information géographique et forestière/Système d'observation du niveau des eaux littorales) maintains 4 stations in the area (2 in Guadeloupe and 2 in Martinique) mainly for reference frame geodesy, with an application to sea level monitoring for the two stations PPTG and FFT2 (Wöppelmann et al. 2011; Santamaría-Gómez et al. 2017). These stations are installed in the vicinity of tide gauges to monitor their vertical motion. The first station of this network is ABMF and was installed in 2008. Finally, the ORPHEON network which consists in 8 stations (5 in Guadeloupe and 3 in Martinique) was installed in 2013/2014 for RTK surveying purposes, but the notable continuity of the time series without any hardware change make them suitable candidates for geophysics use. These stations are provided within the framework of RENAG (Réseau national GNSS permanent, e.g. Walpersdorf et al. 2018; Rabin et al. 2018). We processed the data from May 2000 (deployment of HOUE and SOUF) to end of August 2018 for the IPGP stations and end of November 2018 for the other networks. The timeline of the used observations is represented in Fig. 2.

## 3 Processing

The observation set described below was processed using two different GNSS processing software, namely EPOS (Gendt et al. 1998; Uhlemann et al. 2015) and GINS (Marty et al. 2011; Loyer et al. 2012), but using a similar strategy, equivalent models and identical product inputs. The underlying idea is to quantify the intrinsic differences due to different software. We used a Precise Point Positioning approach with float ambiguity resolution. A PPP processing is the most suitable one for this area because the reduced number of IGS reference stations prevents efficient differential processing. Moreover, those reference stations can be affected by the same tectonic processes as the geophysics stations. Thus, IGS stations in the area (namely ABMF and LMMF) are used here as “regular” stations. We considered only GPS observations since most of the stations recorded only this constellation during most of the period considered. The orbits and clock offset products have been generated beforehand by the GFZ Analysis Center in preparation of the IGS Repro3 campaign (Männel et al. 2020). These products are consistent with ITRF2014 (Altamimi et al. 2016). Regarding the models used, we kept a consistency between the two software configurations. The same antenna eccentricities are used for both processings based on station site logs.
<figure><img src="image_2.png" /><figcaption>**Fig. 1** Maps of the stations in Guadeloupe and Martinique Islands used in this study. The colors represent the three different networks: blue for IPGP, orange for IGN/SONEL, green for ORPHEON/RENAG. (a) Guadeloupe Archipelago. (b) Martinique Island</figcaption></figure>

<figure><img src="image_3.png" /></figure>
<figure><img src="image_4.png" /><figcaption>**Fig. 2** Timeline of the daily observations processed from the three different networks available</figcaption></figure>

Once the two daily coordinate sets are obtained, we select the intersection of both to get equivalent time series with the same daily coordinates. Indeed, some daily data were not properly computed by one or the other software. Stations STG0, SBL0 and PDB0 are completely excluded because of a lack of observations. For each time series, the corresponding velocities are determined using the trend estimation software HECTOR (Bos et al. 2013). We model the time series as combinations of a long term linear trend and an annual+semi-annual periodic signal, along with white and power-law noise. The term *trend* designates hereafter the linear component. The discontinuities introduced in the trend estimation are based on the material change site logs (an antenna change is systematically considered as a discontinuity) and on a supplementary visual detection (Sakic et al. 2019). The same discontinuities are applied for both equivalent EPOS and GINS solutions.

## 4 Coordinate and Velocity Differences

To quantify the impact of the processings, we compute the differences between the two coordinate sets for the three topocentric components and the planimetric distance (Euclidean norm of East and North components). These differences are represented as a histogram in Fig. 3 and the statistical indicators are given in Table 1. We remark that the mean difference for the three components is not centered on zero but is shifted by some millimeters. We also remark that the Up difference doesn't respect a Normal distribution, which reflects the fact that the Up component remains the hardest component to estimate with GNSS technique, mainly because this geometrical parameter is highly correlated to the clock offsets and tropospheric delay parameters. For the planimetric distance, the mean difference is 12.33 mm
<figure><img src="image_5.png" /><figcaption>**Fig. 3** Differences of the three topocentric coordinates and the planimetric distance for the common daily points of the two estimated solutions</figcaption></figure>

<table><caption>**Table 1** Mean, median and standard deviation in millimeters for the three topocentric coordinates and the planimetric distance for the common daily points of the two estimated solutions</caption><thead><tr><td>Component</td><td>Mean</td><td>Std. dev.</td><td>Median</td></tr></thead><tbody><tr><td>East</td><td>-2.89</td><td>12.69</td><td>-2.92</td></tr><tr><td>North</td><td>1.96</td><td>4.17</td><td>2.08</td></tr><tr><td>Up</td><td>4.20</td><td>19.70</td><td>0.15</td></tr><tr><td>Plani. distance</td><td>12.37</td><td>9.91</td><td>9.66</td></tr></tbody></table>

but respects a Gamma distribution. Regarding the standard deviations $\sigma$, $\sigma_{North}$ is three times smaller than the $\sigma_{East}$, which can be explained by the better resolution of this component due to the general North-South trajectory of the satellites.

We compute also the differences between the two sets of estimated velocities obtained at the end of the processing, represented in Fig. 4. For the absolute differences (Fig. 4a), we observe notable differences at the mm/yr level for the planimetric components, and for some stations a difference over $\pm 1$ mm/yr for the vertical component. Regarding the relative differences (Fig. 4b), we note differences on the order of $\pm 10\%$ for the planimetric component, but for the vertical component, these differences can vary by more than a factor of two (stations over 100%). Four stations have a negative relative variation, which means that they have an opposite velocity trend. Since the stations are sorted from the most complete set of data to the sparsest one, we observe no significative relation with the length of the time series.

# 5 Vertical Velocity Results

The vertical velocity values obtained for the two processings are given in Table 2 and represented in Fig. 5 for the Guadeloupe Archipelago, Fig. 6 for the specific area of the Soufrière Volcano, and Fig. 7 for the Martinique Island.

We note general subsidence for both islands. This tendency is consistent between the two solutions. For the Guadeloupe Archipelago, the subsidence is visible for most of the stations, nevertheless, a more complex behavior around the Soufrière area can be remarked, which might be related to local volcanic deformation but also the frequent hardware change due to the harsh conditions in the area (humidity, corrosion and frequent thunderstorms). On the Marie-Galante Island, South-West of the main island, the MGL0 station has a positive trend. Moreover, the four stations TDB0, ABD0 GOSI and DSD0, have opposite trend depending on the solution. For GOSI and DSD0, the velocities estimated in both cases are very close to zero, which make this opposite trend not significant. For GOSI, the different estimated velocities are still close to each other (with overlapping formal sigmas) but a clear velocity tendency for this station is also non-significant. The case of TDB0 is remarkable, since the time series is long and almost complete but the difference between the two solutions is important (4.6 mm/yr). A detailed view of the raw time series and the estimated tendencies are shown in Fig. 8. We observe that the different scatter for both time series lead to a completely different estimation of the trend (using the strategy we selected). Station DHS0 has a lot of corrupted raw data, which lead to a reduced amount of usable observation and an overestimated vertical velocity of almost $-2$ cm/yr for one solution. A similar statement can be made for FFE0 station on the main island, where several gaps in the time series along with several hardware changes might explain the positive trend estimated.

The mean velocity rate measured for all stations on the archipelago, with the volcano area excluded, is $-1.60 \pm 1.54$ mm/yr ($1\sigma$) using EPOS solution, and $-2.17 \pm 1.23$ mm/yr ($1\sigma$) using GINS solution.

For the Martinique Island, the two solutions are also consistent and general subsidence is observed except for SAM0 station. The mean velocity rate measured is
<figure><img src="image_6.png" /><figcaption>**Fig. 4** Velocity difference for the three topocentric coordinates of the two estimated solutions. The stations are sorted from the most complete to the sparsest one. (a) Absolute velocity differences. (b) Relative velocity differences</figcaption></figure>
<table><caption>**Table 2** Vertical velocities estimated for the Guadeloupe and Martinique network, for both EPOS and GINS solutions</caption><thead><tr><th rowspan="2">Station</th><th rowspan="2">Purpose</th><th rowspan="2">Lat.</th><th rowspan="2">Long.</th><th rowspan="2">Start</th><th rowspan="2">End</th><th rowspan="2">Total days</th><th rowspan="2">Used days</th><th rowspan="2">Ratio (%)</th><th rowspan="2">Discont. (#)</th><th colspan="2">EPOS solution (mm/yr)</th><th colspan="2">GINs solution (mm/yr)</th></tr><tr><th>V<sub>Up</sub></th><th>σV<sub>Up</sub></th><th>V<sub>Up</sub></th><th>σV<sub>Up</sub></th></tr></thead><tbody><tr><td>ABD0</td><td>T</td><td>16.47</td><td>298.51</td><td>2012-10-11</td><td>2018-08-30</td><td>2149</td><td>1843</td><td>85.76</td><td>1</td><td>-0.81</td><td>0.99</td><td>0.38</td><td>0.70</td></tr><tr><td>ABER</td><td>S</td><td>16.47</td><td>298.49</td><td>2013-01-30</td><td>2018-11-25</td><td>2125</td><td>1792</td><td>84.33</td><td>0</td><td>-1.27</td><td>0.39</td><td>-0.67</td><td>0.51</td></tr><tr><td>ABMF</td><td>R</td><td>16.26</td><td>298.47</td><td>2008-07-15</td><td>2018-11-25</td><td>3785</td><td>2926</td><td>77.31</td><td>2</td><td>-0.84</td><td>0.40</td><td>-0.39</td><td>0.42</td></tr><tr><td>ADE0</td><td>T</td><td>16.30</td><td>298.91</td><td>2003-01-31</td><td>2015-02-24</td><td>4407</td><td>2504</td><td>56.82</td><td>1</td><td>-0.96</td><td>0.62</td><td>-0.37</td><td>0.58</td></tr><tr><td>AJB0</td><td>V</td><td>14.81</td><td>298.88</td><td>2010-06-29</td><td>2016-11-14</td><td>2330</td><td>1525</td><td>65.45</td><td>1</td><td>-0.89</td><td>1.33</td><td>-0.35</td><td>1.36</td></tr><tr><td>AMC0</td><td>V</td><td>16.05</td><td>298.33</td><td>2015-05-12</td><td>2018-08-30</td><td>1206</td><td>1086</td><td>90.05</td><td>0</td><td>-0.58</td><td>0.84</td><td>-1.52</td><td>0.97</td></tr><tr><td>BIM0</td><td>T</td><td>14.52</td><td>298.93</td><td>2012-08-29</td><td>2016-12-30</td><td>1584</td><td>1399</td><td>88.32</td><td>0</td><td>-0.17</td><td>0.59</td><td>1.12</td><td>0.71</td></tr><tr><td>BOUL</td><td>S</td><td>16.13</td><td>298.23</td><td>2013-01-31</td><td>2018-11-25</td><td>2124</td><td>1680</td><td>79.10</td><td>0</td><td>-3.40</td><td>0.93</td><td>-3.46</td><td>0.95</td></tr><tr><td>CBE0</td><td>T</td><td>16.07</td><td>298.39</td><td>2012-09-30</td><td>2018-08-30</td><td>2160</td><td>1740</td><td>80.56</td><td>0</td><td>-2.17</td><td>1.34</td><td>-0.93</td><td>1.45</td></tr><tr><td>DESI</td><td>S</td><td>16.30</td><td>298.93</td><td>2013-02-04</td><td>2018-11-25</td><td>2120</td><td>1657</td><td>78.16</td><td>0</td><td>-0.41</td><td>0.42</td><td>-0.36</td><td>0.51</td></tr><tr><td>DHS0</td><td>T</td><td>16.29</td><td>298.23</td><td>2009-04-16</td><td>2018-08-30</td><td>3423</td><td>336</td><td>9.82</td><td>2</td><td>-18.20</td><td>2.41</td><td>-8.29</td><td>3.80</td></tr><tr><td>DSD0</td><td>T</td><td>16.31</td><td>298.93</td><td>2011-08-11</td><td>2018-08-30</td><td>2576</td><td>2253</td><td>87.46</td><td>2</td><td>-0.95</td><td>0.83</td><td>0.04</td><td>0.73</td></tr><tr><td>FFE0</td><td>T</td><td>16.22</td><td>298.49</td><td>2003-01-26</td><td>2017-03-13</td><td>5160</td><td>1711</td><td>33.16</td><td>3</td><td>0.56</td><td>1.28</td><td>1.44</td><td>1.42</td></tr><tr><td>FFT2</td><td>R/T</td><td>14.60</td><td>298.94</td><td>2016-01-27</td><td>2017-12-12</td><td>685</td><td>642</td><td>93.72</td><td>0</td><td>-2.84</td><td>1.28</td><td>-3.32</td><td>1.56</td></tr><tr><td>FNA0</td><td>T</td><td>15.88</td><td>298.42</td><td>2004-11-26</td><td>2012-07-25</td><td>2798</td><td>1486</td><td>53.11</td><td>2</td><td>-2.13</td><td>1.11</td><td>-0.78</td><td>0.99</td></tr><tr><td>FNG0</td><td>T</td><td>16.06</td><td>298.31</td><td>2014-12-17</td><td>2018-08-30</td><td>1352</td><td>1173</td><td>86.76</td><td>1</td><td>0.09</td><td>1.91</td><td>-1.29</td><td>1.44</td></tr><tr><td>FSDC</td><td>T</td><td>14.73</td><td>298.85</td><td>2005-02-04</td><td>2016-12-30</td><td>4347</td><td>1831</td><td>42.12</td><td>3</td><td>-5.05</td><td>0.92</td><td>-3.70</td><td>1.12</td></tr><tr><td>GOSI</td><td>S</td><td>16.21</td><td>298.52</td><td>2013-01-30</td><td>2018-11-25</td><td>2125</td><td>1796</td><td>84.52</td><td>0</td><td>-0.16</td><td>0.37</td><td>0.03</td><td>0.44</td></tr><tr><td>HOUE</td><td>T</td><td>15.98</td><td>298.30</td><td>2000-05-14</td><td>2018-08-30</td><td>6682</td><td>4543</td><td>67.99</td><td>3</td><td>-0.98</td><td>0.26</td><td>-1.08</td><td>0.30</td></tr><tr><td>ILAM</td><td>T</td><td>14.77</td><td>299.12</td><td>2012-12-19</td><td>2016-12-30</td><td>1472</td><td>1308</td><td>88.86</td><td>0</td><td>-0.55</td><td>0.66</td><td>-0.09</td><td>0.66</td></tr><tr><td>LAM0</td><td>V</td><td>14.81</td><td>298.84</td><td>2005-02-02</td><td>2016-12-30</td><td>4349</td><td>2045</td><td>47.02</td><td>6</td><td>-4.60</td><td>1.02</td><td>-3.27</td><td>0.80</td></tr><tr><td>LMMF</td><td>R</td><td>14.59</td><td>299.00</td><td>2008-11-09</td><td>2018-11-25</td><td>3668</td><td>3247</td><td>88.52</td><td>2</td><td>-2.40</td><td>0.46</td><td>-2.41</td><td>0.41</td></tr><tr><td>LORI</td><td>S</td><td>14.82</td><td>298.95</td><td>2013-02-06</td><td>2018-11-25</td><td>2118</td><td>1797</td><td>84.84</td><td>0</td><td>-1.02</td><td>0.35</td><td>-1.17</td><td>0.39</td></tr><tr><td>MAD0</td><td>V</td><td>16.01</td><td>298.36</td><td>2016-11-22</td><td>2018-08-30</td><td>646</td><td>610</td><td>94.43</td><td>0</td><td>0.92</td><td>1.78</td><td>1.47</td><td>2.52</td></tr><tr><td>MAGA</td><td>S</td><td>15.89</td><td>298.69</td><td>2013-02-02</td><td>2018-11-25</td><td>2122</td><td>1099</td><td>51.79</td><td>0</td><td>-0.93</td><td>0.53</td><td>-1.14</td><td>0.50</td></tr><tr><td>MARI</td><td>T</td><td>14.47</td><td>299.14</td><td>2013-02-06</td><td>2018-11-25</td><td>2118</td><td>1568</td><td>74.03</td><td>0</td><td>-2.43</td><td>0.40</td><td>-2.73</td><td>0.45</td></tr><tr><td>MGL0</td><td>T</td><td>15.95</td><td>298.72</td><td>2013-08-20</td><td>2017-03-13</td><td>1301</td><td>1209</td><td>92.93</td><td>0</td><td>0.88</td><td>1.41</td><td>1.98</td><td>0.94</td></tr><tr><td>MLM0</td><td>V</td><td>14.78</td><td>298.82</td><td>2011-10-06</td><td>2016-12-30</td><td>1912</td><td>1776</td><td>92.89</td><td>0</td><td>-0.85</td><td>0.93</td><td>-0.44</td><td>1.24</td></tr><tr><td>MPOM</td><td>T</td><td>14.44</td><td>299.14</td><td>2012-11-28</td><td>2016-12-30</td><td>1493</td><td>1170</td><td>78.37</td><td>0</td><td>-1.91</td><td>0.76</td><td>-0.86</td><td>0.70</td></tr><tr><td>PAR1</td><td>V</td><td>16.03</td><td>298.31</td><td>2014-11-21</td><td>2016-12-30</td><td>770</td><td>748</td><td>97.14</td><td>4</td><td>0.65</td><td>4.85</td><td>1.01</td><td>3.81</td></tr><tr><td>PPTG</td><td>R/T</td><td>16.22</td><td>298.47</td><td>2016-01-27</td><td>2018-03-09</td><td>772</td><td>718</td><td>93.01</td><td>0</td><td>-2.79</td><td>1.12</td><td>-2.83</td><td>1.06</td></tr><tr><td>PSA1</td><td>V</td><td>16.04</td><td>298.33</td><td>2011-06-16</td><td>2018-08-30</td><td>2632</td><td>1970</td><td>74.85</td><td>9</td><td>-0.80</td><td>1.14</td><td>2.07</td><td>1.13</td></tr><tr><td>SAM0</td><td>V</td><td>14.84</td><td>298.84</td><td>2013-10-30</td><td>2016-12-21</td><td>1148</td><td>1055</td><td>91.90</td><td>0</td><td>1.28</td><td>1.17</td><td>2.31</td><td>1.22</td></tr><tr><td>SOUF</td><td>V</td><td>16.04</td><td>298.34</td><td>2000-05-13</td><td>2018-08-30</td><td>6683</td><td>4427</td><td>66.24</td><td>3</td><td>-1.23</td><td>0.31</td><td>-1.59</td><td>0.34</td></tr><tr><td>TAR1</td><td>V</td><td>16.04</td><td>298.33</td><td>2016-09-09</td><td>2018-08-30</td><td>720</td><td>659</td><td>91.53</td><td>10</td><td>-2.95</td><td>2.31</td><td>-1.80</td><td>2.97</td></tr><tr><td>TDB0</td><td>T</td><td>15.85</td><td>298.36</td><td>2012-10-18</td><td>2018-08-30</td><td>2142</td><td>1762</td><td>82.26</td><td>2</td><td>3.10</td><td>3.34</td><td>-1.51</td><td>2.64</td></tr><tr><td>TRIL</td><td>S</td><td>14.54</td><td>298.97</td><td>2013-02-05</td><td>2018-11-25</td><td>2119</td><td>1732</td><td>81.74</td><td>0</td><td>-1.42</td><td>0.36</td><td>-1.17</td><td>0.37</td></tr></tbody></table>

The *used days* column refers to the common number of days correctly processed in both solutions, and used for the velocity estimation. The main purpose of each station in mentioned in the second column: we distinguish stations located in the vicinity of the volcano domes, and installed for volcanic deformation monitoring (V), stations for tectonic deformation monitoring (T), stations for reference frame definition and orbit determination (R), and stations for RTK surveying (S)
<figure><img src="image_7.png" /><figcaption>**Fig. 5** Vertical velocity field obtained for the two solutions processed for stations located in Guadeloupe. A green arrow indicates an observed uplift and a red arrow an observed subsidence. The red rectangle</figcaption></figure>

indicates the area around the Soufrière Volcano, presented in detail in Fig. 6. The volcano summit is represented with a brown triangle. Dashed arrows have been shortened to stay in the frame. (a) EPOS solution. (b) GINS solution

<figure><img src="image_8.png" /><figcaption>**Fig. 6** Vertical velocity field obtained for the two solutions processed for stations located in the vicinity of the Soufrière Volcano (Basse-Terre, southern main Island of Guadeloupe). Dashed arrow has been shortened to stay in the frame. (a) EPOS solution. (b) GINS solution</figcaption></figure>

−1.80 ± 1.36 mm/yr (1σ) using EPOS solution, and
−1.68 ± 1.23 mm/yr (1σ) using GINS solution.

# 6 Comparison with Existing Solutions

To validate the consistency of our results, we compare
the vertical velocities we determined with existing
solutions for stations LMMF and ABMF (belonging to
IGS network). We consider the velocities provided by

SONEL for ULR6a (Santamaría-Gómez et al. 2017), NGL14 (Blewitt et al. 2018), JPL14 (Heflin et al. 2019) and ITRF14 (Altamimi et al. 2016) solutions. Values are given in Table 3. For LMMF, the estimated velocities are consistent with each other ($\sigma_V$ = 0.44 mm), while for ABMF we remark more important differences between each solution ($\sigma_V$ = 0.79 mm), just like we have also differences between the values of this work's solutions. Nevertheless, the negative trend remains significant.
<figure><img src="image_9.png" /><figcaption>(a)</figcaption></figure>

**Fig. 7** Vertical velocity field obtained for the two solutions processed for stations located in Martinique. A green arrow indicates an observed uplift and a red arrow an observed subsidence. The volcano summit is

<figure><img src="image_10.png" /><figcaption>(b)</figcaption></figure>

represented with a brown triangle. Dashed arrow has been shortened to stay in the frame. (a) EPOS solution. (b) GINS solution

# 7 Discussion

Using two different solutions but based on the same geodetic products and homogeneous models, we obtain significant disparities in terms of coordinates difference repeatability, especially on the East and Up components with a standard deviation at the centimeter level. Regarding the estimated vertical velocities using the same set of points and the same discontinuities, the differences are also significant. This result tends to motivate investigation on velocity combination strategies between different processing centers, as suggested and tested by Ballu et al. (2019) for instance, where a joint least square modeling is developed to combine equivalent time series from different Analysis Centers. A combination based on a maximum likelihood estimation would be also an relevant method.

Nevertheless, for the studied area of the Guadeloupe and Martinique Islands, a negative velocity trend on the Up component is observed for most of the stations, which might suggest generalized subsidence of the area. This tendency is clear for the Martinique Island, but more complex trends for the Guadeloupe Archipelago can be observed, especially in the area around the Soufrière Volcano. This result can also be nuanced, since some stations have a positive trend, which might be due to local effects. A positive trend can also be

due to an important number of discontinuities over the time series period, like the stations PAR1 (furthermore located inside the volcano area) and FFE0 (outside the volcano area). On another hand, a large number of discontinuities for the same station seem to lead also to an overestimated negative trend, like for instance the station LAM0, with a velocity estimated over -3 mm/yr for six discontinuities referenced. This statement reveals the necessity to maintain networks with a minimum of hardware discontinuities, i.e. by reducing the number of antenna changes. MGL0 station, located on the Marie-Galante Island, presents a singular behavior. It is the only station clearly uplifting, with a quasi-complete time-series of 3.5 years, without any visible discontinuity, while the other station on Marie-Galante (MAGA) presents a subsiding trend. Unfortunately, since this station belong to the commercial ORPHEON network, we have only a few metadata that prevent us to explain clearly this behavior.

We corroborate the paleo-geodesy studies carried out in the region. The coral reef records in Martinique (Weil-Accardo et al. 2016) and Les Saintes (Leclerc and Feuillet 2019) indicate also a subsidence but with a smaller order of magnitude of few tenths of a millimeter per year, which can be explained by the difference in the observation time spans (only a few years for GNSS, ca. one century for the coral records). According to those studies, long term subsidence can have multiple origins: volcanic activity, crustal faulting,
<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /><figcaption>**Fig. 8** Up component time series for the station TDB0, where an opposite velocity trend is visible. Blue dots represent the raw component determined by both software, and green dots the corresponding estimated trend. Red vertical bars represent the discontinuities considered.</figcaption></figure>

(a) EPOS. (b) GINS

<table><caption>**Table 3** Comparison of vertical velocities for LMMF and ABMF of this study with existing solutions</caption><thead><tr><td>mm/yr</td><td>LMMF</td><td>ABMF</td><td>Solution end</td></tr></thead><tbody><tr><td>EPOS</td><td>-2.40 ± 0.46</td><td>-0.84 ± 0.40</td><td>Nov. 2018</td></tr><tr><td>GINS</td><td>-2.41 ± 0.41</td><td>-0.39 ± 0.42</td><td>Nov. 2018</td></tr><tr><td>ULR6a</td><td>-3.55 ± 0.48</td><td>N/A</td><td>Dec. 2014</td></tr><tr><td>NGL14</td><td>-2.70 ± 1.33</td><td>-2.37 ± 1.15</td><td>Apr. 2019</td></tr><tr><td>JPL14</td><td>-2.49 ± 0.54</td><td>-1.74 ± 0.88</td><td>Dec. 2019</td></tr><tr><td>ITRF14</td><td>-2.54 ± 0.21</td><td>-0.92 ± 0.37</td><td>Dec. 2014</td></tr></tbody></table>

subduction of the Tiburon ridge for the Saintes Islands (Leclerc and Feuillet 2019), and a potential deep interseismic loading for Martinique (Weil-Accardo et al. 2016).

We used only one software for velocity estimation since we mainly focussed on the GNSS processing itself, but some other velocity estimator software are available (e.g. Blewitt et al. 2016; Santamaría-Gómez 2019).

The impact of the velocity estimation software on solutions have been analysed for instance by Mazzotti et al. (2020).

## 8 Conclusion

This work brings a comparison of the coordinate time series obtained for the same dataset with two different software but using consistent parameters. New homogeneously calculated vertical velocity fields are made available for geophysical modeling, with unprecedented density for the two Martinique and Guadeloupe Islands. A general subsidence trend is observed for both islands.

**Acknowledgements** We thank the different scientific institutes for the continuously operating GNSS acquisitions, the instrumental mainte-
nance and for providing publicly and freely the observations. Data of
the IPGP network, operated on site by the Guadeloupe and Martinique
seismic and volcanologic observatories, are available on volobsis.ipgp.
fr. Data of the IGN and the SONEL networks are available on rgp.
ign.fr and sonel.org respectively. ORPHEON data were provided to the
authors for a scientific use in the framework of the agreement between
the GEODATA company (orpheon-network.fr) and the RESIF-RENAG
network maintained by the CNRS-INSU (renag.resif.fr). We acknowl-
edge the CNES (Centre national d'études spatiales) for its geodetic
processing software GINS. We thank A. Walpersdorf & an anonymous
reviewer along with the editors J. Freymueller & L. Sanchez, for their
constructive comments which improved the content of this paper.

References

Altamimi Z, Rebischung P, Métivier L, Collilieux X (2016) ITRF2014: a new release of the international terrestrial reference frame modeling nonlinear station motions. J Geophys Res Solid Earth 121(8):6109–6131. https://doi.org/10.1002/2016JB013098. http://doi.wiley.com/10.1002/2016JB013098

Ballu V, Bouin MN, Siméoni P, Crawford WC, Calmant S, Boré JM, Kanas T, Pelletier B (2011) Comparing the role of absolute sea-level rise and vertical tectonic motions in coastal flooding, Torres Islands (Vanuatu). Proc Nat Acad Sci 108(32):13019–13022. https://doi.org/10.1073/pnas.1102842108. http://www.pnas.org/lookup/doi/10.1073/pnas.1102842108

Ballu V, Gravelle M, Wöppelmann G, de Viron O, Rebischung P, Becker M, Sakic P (2019) Vertical land motion in the Southwest and Central Pacific from available GNSS solutions and implications for relative sea levels. Geophys J Int 218(3):1537–1551. https://doi.org/10.1093/gji/ggz247. https://academic.oup.com/gji/article/218/3/1537/5499028

Blewitt G, Kreemer C, Hammond WC, Gazeaux J (2016) MIDAS robust trend estimator for accurate GPS station velocities without step detection. J Geophys Res Solid Earth 121(3):2054–2068. https://doi.org/10.1002/2015JB012552. http://doi.wiley.com/10.1002/2015JB012552

Blewitt G, Hammond W, Kreemer C (2018) Harnessing the GPS data explosion for interdisciplinary science. Eos 99. https://doi.org/10.1029/2018EO104623. https://eos.org/project-updates/harnessing-the-gps-data-explosion-for-interdisciplinary-science

Bos MS, Fernandes RMS, Williams SDP, Bastos L (2013) Fast error analysis of continuous GNSS observations with missing data. J Geod 87(4):351–360. https://doi.org/10.1007/s00190-012-0605-0. http://link.springer.com/10.1007/s00190-012-0605-0

Gendt G, Dick G, Söhne W (1998) GFZ analysis center of IGS–Annual Report. IGS 1996 Annual Report, pp 169–181

Heflin M, Moore A, Murphy D, Desai S, Bertiger W, Haines B, Kuang D, Sibthorpe A, Sibois A, Ries P, Hemberger D, Dietrich A (2019) Introduction to JPL's GPS time series GPS time series. https://sideshow.jpl.nasa.gov/post/tables/GNSS_Time_Series.pdf

Leclerc F, Feuillet N (2019) Quaternary coral reef complexes as powerful markers of longterm subsidence related to deep processes at subduction zones: Insights from Les Saintes (Guadeloupe, French West Indies). Geosphere 15(4):983–1007. https://doi.org/10.1130/GES02069.1

Léticée JL, Cornée JJ, Münch P, Fietzke J, Philippon M, Lebrun JF, De Min L, Randrianasolo A (2019) Decreasing uplift rates and

Pleistocene marine terraces settlement in the central lesser Antilles
fore-arc (La Désirade Island, 16 deg. N). Quat Int 508(March):43–
59. https://doi.org/10.1016/j.quaint.2018.10.030. https://doi.org/10.
1016/j.quaint.2018.10.030

Loyer S, Perosanz F, Mercier F, Capdeville H, Marty JC (2012) Zero-difference GPS ambiguity resolution at CNES-CLS IGS Analysis Center. J Geod 86(11):991-1003. https://doi.org/10.1007/s00190-012-0559-2. http://link.springer.com/10.1007/s00190-012-0559-2

Männel B, Brandt A, Bradke M, Sakic P, Brack A, Nischan T (2020) Status of IGS reprocessing activities at GFZ. In: International Association of Geodesy Symposia. Springer, Berlin. https://doi.org/10.1007/1345_2020_98

Marty JC, Loyer S, Perosanz F, Mercier F, Bracher G, Legresy B, Portier L, Capdeville H, Fund F, Lemoine JM, Biancale R (2011) GINS: the CNES/GRGS GNSS scientific software. In: 3rd international colloquium scientific and fundamental aspects of the Galileo programme, ESA proceedings WPP326, vol 31, pp 8–10

Mazzotti S, Déprez A, Henrion E, Masson C, Masson F, Menut JL, Métois M, Nocquet JM, Rolland L, Sakic P, Socquet A, Santamaria-Gómez A, Valty P, Vergnolle M, Vernant P (2020) Comparative analysis of synthetic GNSS time series - bias and precision of velocity estimations. Research report, RESIF. https://hal.archives-ouvertes.fr/hal-02460380

Mouslopoulou V, Oncken O, Hainzl S, Nicol A (2016) Uplift rate tran-
sients at subduction margins due to earthquake clustering. Tectonics
35(10):2370–2384. https://doi.org/10.1002/2016TC004248. http://
doi.wiley.com/10.1002/2016TC004248

Rabin M, Sue C, Walpersdorf A, Sakic P, Albaric J, Fores B (2018) Present-day deformations of the Jura Arc inferred by GPS surveying and earthquake focal mechanisms. Tectonics 37(10):3782–3804. https://doi.org/10.1029/2018TC005047. http://doi.wiley.com/10.1029/2018TC005047

Sakic P, Mansur G, Kitpracha C, Ballu V (2019) The GeodeZYX toolbox: a versatile Python 3 toolbox for geodetic-oriented purposes. https://doi.org/10.5880/GFZ.1.1.2019.002. http://dataservices.gfzpotsdam.de/panmetaworks/showshort.php?id=escidoc:4754924

Santamaría-Gómez A (2019) SARI: interactive GNSS position time series analysis software. GPS Solutions 23(2):1–6. https://doi.org/10.1007/s10291-019-0846-y. http://dx.doi.org/10.1007/s10291-019-0846-y

Santamaría-Gómez A, Gravelle M, Dangendorf S, Marcos M, Spada G, Wöppelmann G (2017) Uncertainty of the 20th century sea-level rise due to vertical land motion errors. Earth Planet Sci Lett 473:24–32. https://doi.org/10.1016/j.epsl.2017.05.038. http://dx.doi.org/10.1016/j.epsl.2017.05.038

Symithe S, Calais E, de Chabalier JB, Robertson R, Higgins M (2015) Current block motions and strain accumulation on active faults in the Caribbean. J Geophys Res Solid Earth 120(5):3748–3774. https://doi.org/10.1002/2014JB011779. http://doi.wiley.com/10.1002/2014JB011779

Uhlemann M, Gendt G, Ramatschi M, Deng Z (2015) GFZ global multi-GNSS network and data processing results. In: International association of geodesy symposia, vol 12. Springer, New York, pp 673–679. https://doi.org/10.1007/1345_2015_120

Walpersdorf A, Pinget L, Vernant P, Sue C, Deprez A (2018) Does long-term GPS in the Western Alps finally confirm earthquake mechanisms? Tectonics 37(10):3721–3737. https://doi.org/10.1029/2018TC005054. http://doi.wiley.com/10.1029/2018TC005054

Weil-Accardo J, Feuillet N, Jacques E, Deschamps P, Beauducel F, Cabioch G, Tapponnier P, Saurel JM, Galetzka J (2016) Two hun-
dred thirty years of relative sea level changes due to climate and
megathrust tectonics recorded in coral microatolls of Martinique
(French West Indies). J Geophys Res Solid Earth 121(4):2873-
2903. https://doi.org/10.1002/2015JB012406. http://doi.wiley.com/
10.1002/2015JB012406

Wöppelmann G, Testut L, Créach R (2011) La montée du niveau des
océans par marégraphie et géodésie spatiale: contributions françaises
à une problématique mondiale. Annales Hydrographiques 6ème
Série 8(777):11–14

**Open Access** This chapter is licensed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence and indicate if changes were made.

The images or other third party material in this chapter are included in the chapter's Creative Commons licence, unless indicated otherwise in a
credit line to the material. If material is not included in the chapter's Creative Commons licence and your intended use is not permitted by statutory
regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.