Next Article in Journal
Analysis of Various Pickup Coil Designs in Nonmodule-Type GaN Power Semiconductors
Next Article in Special Issue
A Novel Error Correction Approach to Improve Standard Point Positioning of Integrated BDS/GPS
Previous Article in Journal
Autonomous Deployment of Underwater Acoustic Monitoring Devices Using an Unmanned Aerial Vehicle: The Flying Hydrophone
Previous Article in Special Issue
Ionosphere-Constrained Single-Frequency PPP with an Android Smartphone and Assessment of GNSS Observations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Analysis of the Eurasian Tectonic Plate Motion Parameters Based on GNSS Stations Positions in ITRF2014

by
Marcin Jagoda
* and
Miłosława Rutkowska
Faculty of Civil Engineering Environmental and Geodetic Sciences, Koszalin University of Technology, Śniadeckich 2, 75-453 Koszalin, Poland
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(21), 6065; https://doi.org/10.3390/s20216065
Submission received: 10 September 2020 / Revised: 22 October 2020 / Accepted: 23 October 2020 / Published: 25 October 2020
(This article belongs to the Special Issue GNSS Signals and Sensors)

Abstract

:
The article is the fourth part of our research program concerning an analysis of tectonic plates’ motion parameters that is based on an observation campaign of an array of satellite techniques: SLR, DORIS, VLBI, and now GNSS. In this paper, based on the International Terrestrial Reference Frame 2014 (ITRF2014) for observations and using the GNSS technique, the Eurasian tectonic plate motion was analyzed and the plate motion parameters Φ, Λ (the position of the rotation pole), and ω (the angular rotation speed) were adjusted. Approximately 1000 station positions and velocities globally were obtained from the GNSS campaign over a 21-year time interval and used in ITRF2014. Due to the large number of data generated using this technique, the analyses were conducted separately for each tectonic plate. These baseline data were divided into a number of parts related to the Eurasian plate, and are shown in this paper. The tectonic plate model was analyzed on the basis of approximately 130 GNSS station positions. A large number of estimated station positions allowed a detailed study to be undertaken. Stations that agree with the plate motion were selected and plate parameters were estimated with high accuracy. In addition, stations which did not agree with the tectonic plate motion were identified and removed. In the current paper, the influence of the number and location of stations on the computed values and accuracy of the tectonic plate motion parameters is discussed. Four calculation scenarios are examined. Each scenario contains 30 stations for the common solution of the European and Asiatic part of the Eurasian plate. The maximum difference between the four calculation scenarios is 0.31° for the Φ parameter and 0.24° for the Λ parameter, indicating that it is at the level of the value of the formal error. The ω parameter has the same value for all the scenarios. The final stage of the analysis is the estimation of parameters Φ, Λ, and ω based on all of the 120 stations used in the four calculation scenarios (i.e., scenario 1 + scenario 2 + scenario 3 + scenario 4). The following results are obtained: Φ = 54.81° ± 0.37°, Λ = 261.04° ± 0.48°, and ω = 0.2585°/Ma ± 0.0025°/Ma. The results of the analysis are compared with the APKIM2005 model and another solution based on the GNSS technique, and a good agreement is found.

1. Introduction

The earth, due to its dynamics, requires constant observation. Continuous changes occurring inside the earth caused by various geophysical and geological phenomena are reflected on the surface of the globe [1], such as in the form of tectonic plate movements. The theory of the movement of tectonic plates (the so-called continental drift hypothesis) was first published in 1915 by A. Wegener in his fundamental work "Die Entstehung der Kontinente und Ozeane" [2]. In the following years of the 20th century a number of scholars addressed these issues, e.g., geologists O.C. Hilgenberg [3], F.A. Vening-Meinesz [4], S.W. Carey [5], and X. Le Pichon [6], until the adoption by the scientific community in 1968 of the so-called Theory of Plate Tectonics, whose assumptions were described in [7].
During the 1980’s and 1990’s of the last century a dynamic development of satellite observation techniques (GPS, SLR, DORIS, VLBI) occurred, providing an unprecedented opportunity to conduct research on a global scale with very high accuracy. These techniques were also applied in the study of tectonic plate motion. The accuracy of determining the position of points on the earth’s surface using satellite methods has continuously increased, and this increase in accuracy has allowed the numerical values describing the movement of tectonic plates to be determined. Continuity of research in this field is a noticeable trend, e.g., GNSS technique: [8,9,10,11,12]; SLR technique: [13,14,15,16,17,18]; DORIS technique: [10,19,20]; and VLBI technique: [21,22,23,24,25,26].
Initially, geological methods based on the description of the mechanism of tectonic plate motion were used to study the movement of tectonic plates; only the satellite techniques mentioned above made it possible to precisely estimate the numerical values of this movement, as discussed in [27]. Thorough knowledge and the ability to mathematically describe the phenomena causing the movement of the tectonic plates, and the possibility of taking precise measurements using satellite techniques, allows points on the earth’s surface and their annual movements to be determined with high accuracy. This is the basis for determining the parameters that describe the movement of tectonic plates. Tectonic plates move in relation to each other on the astenospheric surface at a speed ranging from a few to several centimeters per year. The lithosphere is constantly destroyed in the subduction zones and is renewed on mid-ocean ridges. Continents usually form a section of a plate and are moved with it. The following types of boundaries between the plates exist: divergent, convergent, and horizontally sliding. Issues related to these boundaries have been previously addressed, e.g., [28,29]. Large stresses can be present on the edges of the plates. These are discharged in the form of earthquakes. In the case of tectonic plates, it is difficult to speak of rigid boundaries between plates and to indicate where one plate ends and another begins; at the plate boundaries, there is often an area in which displacements are incompatible with the movement of the whole plate.
The movement of the tectonic plate is described by the rotation vector Ω, whereas the parameters of this movement are described by the geographical position of the rotation pole—Φ and Λ—and the angular velocity of rotation—ω (or by individual components of the rotation pole) [30,31]. According to [30], the displacement of the observational station as a function of the parameters of plate motion is expressed by Equation (1):
Δ ϕ = ω Δ t cos Φ sin ( λ Λ ) Δ λ = ω Δ t ( sin Φ cos ( λ Λ ) tan ϕ cos Φ ) ,
where:
  • Δϕ, Δλ—displacement of the observational station position in latitude and longitude;
  • ϕ, λ—observational station position (GNSS, SLR, DORIS, VLBI);
  • Φ, Λ, ω—plate motion parameters (the position of the rotation pole in latitude and longitude, and the angular rotation speed, respectively).
The aim of this study is to estimate and analyze the parameters describing the motion of the Eurasian plate based on GNSS station positions, and to identify and eliminate from the calculation those stations whose motion is not consistent with the motion of the plate, hence “polluting” the results of the motion parameters and increasing the value of the formal error.

2. Materials and Methods

There are approximately 1000 GNSS stations located around the world whose positions are provided in the International Terrestrial Reference Frame 2014 (ITRF2014) [32]. From this baseline data, about 130 stations located on the Eurasian continent were analyzed. These stations are not uniformly distributed. Approximately 70% of the stations are located in the European part of the Eurasian plate. Generally, however, the number of estimated station positions and velocities is sufficiently large, and allows plate motion parameters to be estimated with high accuracy and in detail. It also allows stations that are not compatible with the Eurasian tectonic plate motion to be identified and removed, and only stations that are consistent with the solution to be selected.
The baseline data were divided into four calculation scenarios. Each scenario contained 30 randomly distributed stations for a common solution of the Eurasian plate. The stations used in the four scenarios of the calculations are tabulated in Appendix A, Appendix B, Appendix C, Appendix D, separately for each scenario, and are shown in Figure 1. Detailed information on individual stations, including their positions and velocities, is presented on the ITRF website: http://itrf.ensg.ign.fr/ITRF_solutions/2014/. The final solution of the calculation was determined from all of the 120 stations (scenario 1 + scenario 2 + scenario 3 + scenario 4). The locations of the GNSS stations that are not consistent with the Eurasian plate motion are shown in Figure 2 and Figure 3. All calculations were carried out using the authors’ own software.
The description of the method applied in this paper to determine the plate motion parameters, and the theoretical basis for this motion, has been extensively presented by the authors in their earlier papers on the analysis of the plate motion based on movements of observation stations: SLR [33,34]; DORIS [34,35]; and VLBI [34,36]. Thus these aspects were omitted herein.

3. Tectonic Plate Theory for the Eurasian Plate

The Eurasian plate, which is the subject of this analysis, is the most complex part of the earth in terms of geological structure [37]. On the western, northern, and north-eastern sides, the plate borders the North American plate. In the west and north, this boundary is a divergent one, on which the Mid-Atlantic Ridge was formed. In the east it is a convergent boundary [38]. In the north, the area of Sweden and southern Norway has been subject to continuous uplift movements for about 1.5 billion years to date [39]. The north-eastern part of the Asian continent belongs to the North American plate, whereas the area comprising the Koryak Mountains and the Kamchatka area belongs to the Eurasian plate; however, the course of the boundary itself has not yet been sufficiently studied. The Eurasian plate also borders on the Bering and Pacific plates, which are being pulled under the Eurasian plate, constituting is a convergent boundary. Two small plates are positioned between the Pacific and Eurasian plates: the Okhotsk plate, including the island of Sakhalin on which tectonic movements are ongoing, and the Amur plate. These plates are often treated as part of the Eurasian plate because their motion is similar to that of the latter plate [40]. A zone of islands with heterogeneous crust structure extends from Hokkaido in the north to Taiwan in the south. This string lies at the meeting point of the Pacific, Philippine, Okhotsk, and Amur plates. The boundary between the Philippine and Pacific plates is a convergent one, running through the island of Honshu. In Indonesia, eastern Asia joins the Pacific range. The southern boundary of the Eurasian plate with the African, Anatolian, Arabic, and Indo-Australian plates is defined by approximately 12,000 km of the Alpine-Himalayan range, the Alpine fold zone connecting through the Strait of Gibraltar with the Alpine structures of the Atlas mountains, and the Alpine ranges in Asia. This boundary includes the Atlas, the Pyrenees, the Alps, the Apennines, the Carpathians, the Rhodopes, the Caucasus, the Tien Shan, Tibet, the Himalayas as far as Burma, and the Andaman Sea. The Apennines were created in the process of subduing the Tyrrhenian and Adriatic plates. The Caucasus was created as a result of a collision between the Anatolian and European plates. The area is characterized by high seismicity and is still subject to tectonic processes. In this region, the plates overlap and are squeezed and folded, the continents collide, and the phenomenon of orogenesis takes place. Now, this phenomenon occurs at the place where the Indian and Eurasian plates meet. This collision began about 40 million years ago. Since that contact, the Indian Peninsula has moved 2000 km into Eurasia. The consequences of this collision include the creation of the Himalayas, Tibet, and Tien Shan, in addition to earthquakes in China and north-western Iran [41,42]. The processes associated with the formation of the continental structure continue to this day. The Indian Peninsula moves northwards at a speed of 3.3–4.8 cm per year, resulting in a continuous uplift of the Himalayas. The Arabian Peninsula moves to the northeast at a speed of 1.4–1.8 cm per year. The Pacific plate is subdued by the Asian continental lithosphere at a speed of 6.7–7 cm per year. These processes result in strong earthquakes, tsunamis, and volcanic eruptions [43].
For the Eurasian, African, and Arabic plates, the whole area of southern Europe and the Mediterranean area is considered a boundary zone. The lithosphere here is very cracked and forms a large number of micro plates [44]. The African plate is rubbing against the Eurasian and Arabic plates. The area is highly seismically active and the lithosphere is strongly cracked and divided into the Anatolian, Black Sea, Aegean, Macedonian, Adriatic, Corsican, Atlas, and Iberian micro plates. The north-eastern part of Corsica belongs to the Alpine structures and came into existence as a result of collisions with the micro plate of the Tyrrhenian Sea (the Apulian plate) and the African plate.

4. Results and Discussion

This section presents the results of the estimation of the Eurasian tectonic plate motion parameters based on GNSS station positions in ITRF2014 [32]. Initially, the solution was derived for two stations located on a plate. Then, single stations were added individually (sequential solution) for a maximum of 29 steps for 30 stations (in each of the four scenarios of calculations). The selection of the stations in the scenarios was based on the following criterion: first stations located in the stable region of the plate were included; then, stations located near the plate boundary but which did not disturb the solution accuracy were included. Results for the four calculation scenarios are given in Figure 4, Figure 5 and Figure 6 and are summarized in tables (Appendix A, Appendix B, Appendix C, Appendix D) for values of plate motion parameters Φ, Λ, and ω. The stability of the parameters and their errors for scenarios 1–4 is observed for about 18–20 stations. For the subsequently added stations, the estimated plate parameter values differ by less than the value of the formal error. This proves that the necessary number of stations is about 20. The final solutions based on 30 stations for the four calculation scenarios are as follows: parameter Φ for scenario 1 equals 54.87° ± 0.47°, for scenario 2 equals 54.84° ± 0.38°, for scenario 3 equals 54.56° ± 0.36°, and for scenario 4 equals 54.69° ± 0.41°. The maximum difference amounts to 0.31°, which is consistent with the level of the value of the formal error. Parameter Λ for scenario 1 equals 261.02° ± 0.50°, for scenario 2 equals 260.98° ± 0.57°, for scenario 3 equals 261.03° ± 0.46°, and for scenario 4 equals 261.22° ± 0.53°. The maximum difference between the scenarios amounts to 0.24°, which is less that the value of the formal error. Similarly, for ω the result equals 0.2585°/Ma for all of the scenarios. This indicates the high consistency of these solutions after elimination of erroneous stations. First, erroneous stations and those not compatible with the Eurasian plate motion or not located on the Eurasian plate were identified and deleted from the analysis. The analysis concerning selection of these stations was performed for calculation scenario 1 (given in Appendix A) and is presented in Table 1. The second line of this table shows the result of the final solution for calculation scenario 1. Next, “doubtful” stations were analyzed. Lines 3 to 13 of Table 1 provide the results obtained after the addition of 30 individually selected “doubtful” stations given in scenario 1 and the new solution resulting from 31 stations. Adding Bilibino (BILI) station located on the North American plate on which part of the Asian continent is situated (geodetic position is approx. B = 67°, L = 166°) to the solution changes the plate parameters by about 1.8° and 1.4° in latitude and longitude, respectively. This station should not be included in the solution because differences are greater than the values of error. Next, Magadan (MAGO) station is located on the boundary of two plates, namely, the North American plate and the small Okhotsk plates. The geodetic position of the MAGO station is equal to approx. B = 59°, L = 150°. This station changes the plate parameters by about 1.2° and 1.3° in latitude and longitude, respectively. Differences are greater than the values of error, therefore this station should not be included in the solution. Station Petropavlovsk (PETP) is located on the Kamchatka Peninsula, and its geodetic position is equal to approx. B = 52°, L= 158°. This station changes the plate parameters significantly, by about 4.5° and 6.8° in latitude and longitude, respectively. The differences are about 9 times greater than the values of error; hence, this station is eliminated from the solution. Teheran (TEHN) station is located near the boundary of the Iran plate in the geodetic position of approx. B = 35°, L = 51°. This station changes plate parameters insignificantly, by about 0.2° and 0.4° in latitude and longitude, respectively, but also increases the values of error of the determined parameters and is not be included in the solution. Next, Kunming (KUNM) station is located in South Himalaya in a very seismically active region of China on the Eurasian plate, and has a geodetic position of approx. B = 24°, L = 102°. The station changes plate parameters by about 0.8° and 1° in latitude and longitude, respectively, and is not used in the solution. Lhasa (LHAS) station is located on the boundary of the Eurasian plate and the Deccan plate in South Himalaya, and its geodetic position is approx. B = 29°, L = 91°. The station significantly increases the error of the plate parameters, by about 2.8° and 4° in latitude and longitude, respectively. Due to these values this station is eliminated from the solution. The Japanese stations of Mizusawa (MIZU), Koganei (KGNI), Usuda (USUD), Tsukuba (TSKB), Mitaka (MTKA), and Abashiri (P202) are located at B = 35–44°, L = 138–144° near the boundary of the Pacific, Okhotsk, Amur, and Philippine plates. Shifts of these stations are not consistent with the Eurasian plate motion and values of motion vectors are significantly smaller, thus, these stations are not included in the solution.
The stations with velocity directions and values that are not consistent with the Eurasian plate motion are not included in the solution (these are shown in Figure 2 and Figure 3). Inclusion of these stations results in a description of the Eurasian plate motion that is incompatible with the real movement of the whole plate rotation around poles Φ and Λ.
The final stage of the analysis was the estimation of parameters Φ, Λ, and ω based on all of the 120 stations used in the four calculation scenarios (scenario 1 + scenario 2 + scenario 3 + scenario 4). The following results were obtained: Φ = 54.81° ± 0.37°, Λ = 261.04° ± 0.48°, and ω = 0.2585°/Ma ± 0.0025°/Ma. These results were compared with the APKIM2005 model by H. Drewes [45] and with another solution based on the GNSS technique given in Larson et al. [8]. The comparison is shown in Table 2. The APKIM2005 model [45] is based on weekly solutions for SLR, DORIS, VLBI, and GNSS techniques observed for the time interval from 1993 to 2004. The database for these techniques allows estimation of the velocity vectors caused by plate motion and plate motion parameters for the following 17 plates: African, Amur, Antarctic, Arabian, Anatolian, Australian, Caribbean, Eurasian, Indian, Nazca, North American, South American, Okhotsk, Pacific, Somalia, Sunda, and Yangtze. In [8] authors analyzed the GNSS data from January 1991 to March 1996. All of the presented data were analyzed using the GIPSY/OASIS II software. On the basis of the computed network, velocities for 38 sites located on the African, Antarctic, Australian, North American, South American, Pacific, and Eurasian plates were estimated. For the Eurasian plate, the parameters were estimated on the basis of eight stations; six of these are located on the European part of the plate, and two are located on the Asian plate.
Plate motion parameters obtained by Drewes for the Eurasian plate equal: Φ = 53.4° ± 0.4°, Λ = 264.3° ± 0.5°, and ω = 0.259°/Ma ± 0.001°/Ma [45]. However, in Larson’s determination [8] they are: Φ = 56.3°, Λ = 257.2°, and ω = 0.26°/Ma. The differences between these three solutions are approximately 3 and 7° for Φ and Λ, respectively, whereas for ω parameter the difference amounts to 0.0015°/Ma.

5. Conclusions

The computations and results of the analysis given in this paper allow the following conclusions to be drawn:
  • The plate motion parameters for the final solution based on 120 (scenario 1 + scenario 2 + scenario 3 + scenario 4) GNSS station positions taken from ITRF2014 are equal to Φ = 54.81° ± 0.37° for latitude, Λ = 261.04° ± 0.48° for longitude, and ω = 0.2585°/Ma ± 0.0025°/Ma for rotation speed.
  • The convergence and stability of the solutions for the four calculation scenarios based on 30 randomly distributed stations for each scenario are obtained for about 18–20 stations, as presented in Figure 4 for latitude Φ, Figure 5 for longitude Λ, and Figure 6 for rotation speed ω.
  • The selection of suitable stations for determining the parameters allows the determination on the basis of approximately 20 stations to be made, as shown in the four calculation scenarios. Adding more stations to the calculations results in a change in the value of the determined parameters by a value that does not exceed the formal error.
  • The geological model of the Eurasian plate is the most complicated on earth. The Eurasian plate includes a subduction zone in the convergent eastern boundary. Part of the North American plate, on which Bilibino (BILI) station is situated, is located in the eastern section of the Eurasian continent. The annual shift of this station is not consistent with the Eurasian plate motion. This station cannot be included in the solution because of a change in the pole rotation of the plate parameters by approx. 2 degrees in latitude, which corresponds to 200 km, and 1.5 degrees in longitude, which corresponds to 150 km. Magadan (MAGO) station is located on the boundary of the American and Okhotsk plates. The annual shift of the station is not consistent with the Eurasian plate motion by approx. 1 degree in latitude, which equals 100 km, and 1 degree in longitude, which also equals 100 km, therefore this station cannot be used in the solution. Petropavlovsk (PETP) station is located on the Okhotsk plate on Kamchatka Peninsula. The shift of this station is also not consistent with the Eurasian plate motion, and amounts to 4.5 degrees in latitude, which corresponds to 450 km, and 6.8 degrees in longitude, which equals 680 km, hence the station cannot be used in the solution. Teheran (TEHN) station is also not used in the solution because it is located near the boundary of the Iran plate in a high seismic activity region, despite a relatively small change in the pole rotation of plate parameters, by approx. 0.2 degrees in latitude, which corresponds to 20 km, and 0.4 degrees in longitude, which equals 40 km.
  • The Eurasian plate contains areas of high seismic activity and cracked boundaries, both convergent and divergent, with other tectonic plates, such as the Pacific, Deccan, Iran, Arabian, and Anatolian plates. Often, shifts of the analyzed stations located on the Eurasian plate are not compatible with the tectonic plate motion; for example, Lhasa (LHAS) station, which is located on the boundary of the Eurasian plate and the Deccan plate in the South Tibet and Himalaya mountains, is in a very active and cracked region. This station increases the error of the estimated plate parameters significantly, by approx. 3 degrees in latitude, equaling 300 km, and 4 degrees in longitude, equaling 400 km. The station is not included in the solution.
  • The Japanese stations of Mizusawa (MIZU), Koganei (KGNI), Usuda (USUD), Tsukuba (TSKB), and Abashiri (P202) are located in a very active and cracked area on the boundary of the Philippine plate. In this region, three cracked plates come into contact: the Pacific plate with the Okhotsk and Amur plates. Shifts of these stations are not consistent with the Eurasian plate motion and are not included in the solution.
  • Wakkanai station (used in calculation scenario 4) is located on the boundary of the Amur and Okhotsk plates. The shift of this station is consistent with the Eurasian plate motion to a high degree, therefore it is included in the solution.
  • The application of the sequential calculation method allows stations whose movement is not consistent with that of the entire plate to be identified and eliminated from the solution.
  • Analysis and elimination of the selected stations, as shown in Table 1, is indispensable because it allows results to be obtained only on the basis of the stations with shifts that are consistent with the Eurasian plate.
  • The Mediterranean Sea and surrounding areas, i.e., the Anatolian and Arabic plates, are excluded from the analysis in this work. These regions will be analyzed separately because they include a significant number of micro plates, which each has a characteristic motion that differs from that of the Eurasian plate. In the Asian part, the Deccan plate is omitted because it does not belong to the Eurasian plate.

Author Contributions

Conceptualization, M.J.; methodology, M.J.; software, M.R.; validation, M.J., M.R.; formal analysis, M.J.; investigation, M.J.; resources, M.J. and M.R.; data curation, M.J. and M.R.; writing—original draft preparation, M.J.; writing—review and editing, M.J. and M.R.; visualization, M.J.; supervision, M.J.; project administration, M.J.; funding acquisition, M.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Computed plate motion parameters, and their errors, of the Eurasian plate for the GNSS network using the sequential method. Scenario 1 of the calculations.
Table A1. Computed plate motion parameters, and their errors, of the Eurasian plate for the GNSS network using the sequential method. Scenario 1 of the calculations.
No.Name of the Station Φ [°] Λ [°] ω [°/Ma]
2Ny-Alesund+Bucarest53.04 ± 2.64262.45 ± 1.790.2659 ± 0.0183
32 + Paris53.61 ± 2.35261.37 ± 1.490.2643 ± 0.0151
43 + Yerevan54.00 ± 1.90262.84 ± 1.330.2642 ± 0.0136
54 + Simeiz53.72 ± 1.69261.46 ± 1.200.2637 ± 0.0109
65 + Copenhaguen53.21 ± 1.45261.73 ± 1.100.2609 ± 0.0098
76 + Uzhgorod53.49 ± 1.28262.01 ± 1.010.2635 ± 0.0077
87 + Graz53.82 ± 1.10262.22 ± 0.940.2634 ± 0.0070
98 + Borowiec54.00 ± 0.97262.09 ± 0.860.2627 ± 0.0067
109 + Venezia54.06 ± 0.87262.44 ± 0.810.2630 ± 0.0059
1110 + Wettzell54.28 ± 0.80262.11 ± 0.78 0.2625 ± 0.0052
1211 + Lviv54.11 ± 0.77261.94 ± 0.730.2621 ± 0.0048
1312 + Sofia54.35 ± 0.74261.62 ± 0.710.2615 ± 0.0045
1413 + Novossibirsk54.21 ± 0.72261.27 ± 0.690.2598 ± 0.0040
1514 + Marseille54.40 ± 0.68261.36 ± 0.680.2597 ± 0.0040
1615 + Shanghai54.62 ± 0.67261.49 ± 0.640.2594 ± 0.0038
1716 + Pecny54.73 ± 0.66261.42 ± 0.630.2595 ± 0.0038
1817 + Jozefoslaw54.70 ± 0.63261.24 ± 0.570.2593 ± 0.0037
1918 + Domen54.70 ± 0.60261.20 ± 0.550.2590 ± 0.0035
2019 + Ulan Bator54.76 ± 0.60261.15 ± 0.530.2589 ± 0.0033
2120 + Perugia54.78 ± 0.57261.07 ± 0.530.2588 ± 0.0032
2221 + Kitab54.69 ± 0.57261.10 ± 0.530.2587 ± 0.0031
2322 + Grasse54.63 ± 0.51261.14 ± 0.530.2585 ± 0.0029
2423 + Madrid Robledo54.46 ± 0.50261.25 ± 0.530.2585 ± 0.0028
2524 + Zimmerwald54.55 ± 0.48261.30 ± 0.520.2585 ± 0.0028
2625 + Istanbul54.65 ± 0.48261.18 ± 0.520.2583 ± 0.0027
2726 + Klaipeda 54.70 ± 0.48261.09 ± 0.520.2585 ± 0.0026
2827 + Irkoutsk54.84 ± 0.47261.03 ± 0.510.2585 ± 0.0026
2928 + Herstmonceux54.90 ± 0.47261.07 ± 0.500.2585 ± 0.0026
3029 + Chumysh54.87 ± 0.47261.02 ± 0.500.2585 ± 0.0026

Appendix B

Table A2. Computed plate motion parameters, and their errors, of the Eurasian plate for the GNSS network using the sequential method. Scenario 2 of the calculations.
Table A2. Computed plate motion parameters, and their errors, of the Eurasian plate for the GNSS network using the sequential method. Scenario 2 of the calculations.
No.Name of the Station Φ [°] Λ [°] ω [°/Ma]
2Metsahovi+Matera53.64 ± 2.75262.33 ± 3.300.2644 ± 0.0259
32 + Trondheim53.18 ± 2.25261.88 ± 2.770.2621 ± 0.0176
4 3 + Potsdam53.55 ± 1.99261.65 ± 2.330.2593 ± 0.0136
5 4 + Toulouse53.74 ± 1.54261.34 ± 1.860.2577 ± 0.0123
65 + La Rochelle54.13 ± 1.31261.00 ± 1.550.2586 ± 0.0093
76 + Villafranca54.88 ± 1.00260.88 ± 1.230.2574 ± 0.0085
87 + Torino I55.10 ± 0.86260.93 ± 1.150.2601 ± 0.0080
98 + Changchun55.32 ± 0.79260.52 ± 1.100.2587 ± 0.0076
109 + Wuhan55.06 ± 0.74261.07 ± 1.020.2589 ± 0.0069
1110 + Onsala55.00 ± 0.71261.43 ± 1.000.2582 ± 0.0063
1211 + Vaasa54.95 ± 0.66261.05 ± 0.970.2585 ± 0.0059
1312 + Penc54.89 ± 0.62260.83 ± 0.960.2581 ± 0.0054
1413 + Wroclaw54.79 ± 0.62260.82 ± 0.960.2580 ± 0.0052
1514 + Nicea54.72 ± 0.60260.96 ± 0.950.2582 ± 0.0048
1615 + Chize54.80 ± 0.54260.92 ± 0.920.2587 ± 0.0045
1716 + Boras54.87 ± 0.53261.03 ± 0.880.2596 ± 0.0044
1817 + Stavanger54.86 ± 0.50261.20 ± 0.860.2590 ± 0.0042
1918 + Mendeleevo54.81 ± 0.49261.14 ± 0.860.2594 ± 0.0040
2019 + Port Vendres54.84 ± 0.44261.05 ± 0.810.2592 ± 0.0039
2120 + Cagliari54.86 ± 0.43261.05 ± 0.720.2587 ± 0.0036
2221 + Teddington54.86 ± 0.42260.95 ± 0.690.2585 ± 0.0036
2322 + Roquetes54.81 ± 0.41260.92 ± 0.660.2586 ± 0.0034
2423 + Cannes54.81 ± 0.40260.96 ± 0.640.2586 ± 0.0033
2524 + Kootwijk54.82 ± 0.40261.01 ± 0.620.2585 ± 0.0031
2625 + Terschelling54.82 ± 0.39260.97 ± 0.620.2584 ± 0.0031
2726 + Westerbork54.85 ± 0.39261.01 ± 0.610.2585 ± 0.0029
2827 + Braunschweig54.83 ± 0.38261.00 ± 0.590.2585 ± 0.0029
2928 + Poltava54.80 ± 0.38261.02 ± 0.570.2585 ± 0.0029
3029 + Tixi Seismic54.84 ± 0.38260.98 ± 0.570.2585 ± 0.0029

Appendix C

Table A3. Computed plate motion parameters, and their errors, of the Eurasian plate for the GNSS network using the sequential method. Scenario 3 of the calculations.
Table A3. Computed plate motion parameters, and their errors, of the Eurasian plate for the GNSS network using the sequential method. Scenario 3 of the calculations.
No.Name of the Station Φ [°] Λ [°] ω [°/Ma]
2Brest+Joensuu51.62 ± 0.85255.11 ± 0.980.2333 ± 0.0245
32 + Padova53.26 ± 0.77260.16 ± 0.920.2559 ± 0.0230
43 + Redu53.45 ± 0.61260.43 ± 0.740.2552 ± 0.0145
54 + Yebes54.17 ± 0.60260.87 ± 0.620.2552 ± 0.0140
65 + Wladyslawowo54.30 ± 0.59260.76 ± 0.580.2551 ± 0.0135
76 + Warnemuende54.23 ± 0.59261.01 ± 0.560.2564 ± 0.0130
87 + Sassnitz54.52 ± 0.58261.26 ± 0.550.2567 ± 0.0120
98 + Oberpfaffenhofe54.52 ± 0.58261.44 ± 0.530.2567 ± 0.0115
109 + Osan air base54.63 ± 0.57261.72 ± 0.520.2567 ± 0.0105
1110 + Badary54.77 ± 0.56261.79 ± 0.520.2568 ± 0.0100
1211 + Hsinchu54.71 ± 0.56261.61 ± 0.500.2568 ± 0.0095
1312 + Zwenigorod54.88 ± 0.55261.40 ± 0.500.2572 ± 0.0090
1413 + Daejeon55.00 ± 0.55261.45 ± 0.500.2571 ± 0.0073
1514 + Suwon-Shi54.81 ± 0.54261.28 ± 0.490.2571 ± 0.0057
1615 + Vilhelmina54.66 ± 0.54261.52 ± 0.490.2576 ± 0.0055
1716 + Shaanxi54.62 ± 0.53261.34 ± 0.480.2577 ± 0.0040
1817 + Beijng54.80 ± 0.52261.14 ± 0.480.2626 ± 0.0037
1918 + Borkum54.83 ± 0.49261.17 ± 0.480.2626 ± 0.0035
2019 + Helgoland Islan54.79 ± 0.46261.12 ± 0.480.2626 ± 0.0028
2120 + Wabern54.67 ± 0.43261.06 ± 0.480.2626 ± 0.0021
2221 + Cascais54.78 ± 0.42261.08 ± 0.470.2600 ± 0.0020
2322 + Palma De mallor54.80 ± 0.40261.16 ± 0.470.2594 ± 0.0020
2423 + Alicante54.68 ± 0.39261.07 ± 0.470.2585 ± 0.0019
2524 + Morpeth54.65 ± 0.38261.05 ± 0.470.2585 ± 0.0019
2625 + Brussels Ukkle54.59 ± 0.37260.99 ± 0.460.2585 ± 0.0018
2726 + Moscow54.56 ± 0.37261.12 ± 0.460.2585 ± 0.0018
2827 + Khabarovsk54.56 ± 0.36261.00 ± 0.460.2585 ± 0.0018
2928 + Esbjerg54.55 ± 0.36261.07 ± 0.460.2585 ± 0.0018
3029 + Krasnoyarsk54.56 ± 0.36261.03 ± 0.460.2585 ± 0.0018

Appendix D

Table A4. Computed plate motion parameters, and their errors, of the Eurasian plate for the GNSS network using the sequential method. Scenario 4 of the calculations.
Table A4. Computed plate motion parameters, and their errors, of the Eurasian plate for the GNSS network using the sequential method. Scenario 4 of the calculations.
No.Name of the Station Φ [°] Λ [°] ω [°/Ma]
2Ajaccio+Saint Jean des55.09 ± 1.86260.45 ± 2.060.2507 ± 0.0119
32 + Martsbo55.30 ± 0.75262.61 ± 1.290.2541 ± 0.0056
43 + Visby55.56 ± 0.69262.08 ± 1.250.2566 ± 0.0052
54 + Innsbruck Haf55.22 ± 0.68261.42 ± 1.190.2583 ± 0.0050
65 + Borowa Gora55.03 ± 0.67262.95 ± 1.120.2592 ± 0.0048
76 + Lamkowko55.28 ± 0.64262.20 ± 1.000.2610 ± 0.0044
87 + Riga55.16 ± 0.62261.73 ± 0.900.2622 ± 0.0040
98 + Kharkiv54.99 ± 0.62261.06 ± 0.840.2637 ± 0.0037
109 + Tashkent54.83 ± 0.61261.58 ± 0.700.2641 ± 0.0033
1110 + Mikolajev54.75 ± 0.61261.35 ± 0.670.2646 ± 0.0033
1211 + Svetloe55.00 ± 0.61261.15 ± 0.660.2654 ± 0.0033
1312 + Zelenchukskaya55.24 ± 0.61261.49 ± 0.650.2585 ± 0.0032
1413 + Redzikowo55.15 ± 0.61262.00 ± 0.650.2586 ± 0.0031
1514 + Arti55.00 ± 0.60262.11 ± 0.650.2605 ± 0.0031
1615 + Poligan Bishk54.76 ± 0.60262.05 ± 0.640.2615 ± 0.0031
1716 + Morpeth54.99 ± 0.58261.84 ± 0.640.2611 ± 0.0030
1817 + Genova54.90 ± 0.57261.40 ± 0.620.2602 ± 0.0029
1918 + Lowestoft54.78 ± 0.55261.36 ± 0.610.2586 ± 0.0029
2019 + North Shields54.73 ± 0.55261.31 ± 0.600.2594 ± 0.0028
2120 + Ganovce54.64 ± 0.51261.27 ± 0.600.2594 ± 0.0027
2221 + Vacov54.54 ± 0.49261.21 ± 0.590.2585 ± 0.0025
2322 + Bellmunt de Seg54.56 ± 0.47261.33 ± 0.580.2585 ± 0.0025
2423 + Cantabria54.71 ± 0.46261.35 ± 0.570.2586 ± 0.0024
2524 + Valencia 54.60 ± 0.44261.25 ± 0.550.2585 ± 0.0024
2625 + Thessaloniki54.52 ± 0.43261.22 ± 0.540.2585 ± 0.0024
2726 + Kiev54.72 ± 0.42261.35 ± 0.540.2585 ± 0.0024
2827 + Liverpool54.67 ± 0.42261.23 ± 0.530.2585 ± 0.0024
2928 + Norilsk54.64 ± 0.42261.29 ± 0.530.2585 ± 0.0024
3029 + Wakkanai54.69 ± 0.41261.22 ± 0.530.2585 ± 0.0024

References

  1. Rutkowska, M.; Jagoda, M. SLR technique used for description of the Earth elasticity. Artif. Satell. 2015, 50, 127–141. [Google Scholar] [CrossRef] [Green Version]
  2. Wegener, A. Die Entstehung der Kontinente und Ozeane; Friedr. Vieweg & Sohn Akt.-Ges: Braunschweig, Germany, 1915. [Google Scholar]
  3. Hilgenberg, O.C. Vom Wachsenden Ertball; Selbstverlag: Berlin, Germany, 1933. [Google Scholar]
  4. Vening-Meinesz, F.A. Major tectonic phenomena and the hypothesis of convection currents in the Earth. J. Geol. Soc. 1948, 103, 191–207. [Google Scholar] [CrossRef]
  5. Carey, S.W. The Tectonic Approach to Continental Drift; Geology Department Symposium University: Tasmania, Australia, 1956. [Google Scholar]
  6. Le Pichon, X. Sea floor spreading and continental drift. J. Geophys. Res. 1968, 73, 3661–3697. [Google Scholar] [CrossRef]
  7. Cox, A.; Hart, R.B. Plate Tectonics: How It Works; John Wiley & Sons: Hoboken, NJ, USA, 1986. [Google Scholar]
  8. Larson, K.M.; Freymueller, J.T.; Philipsen, S. Global plate velocities from the Global Positioning System. J. Geophys. Res. Solid Earth 1997, 102, 9961–9981. [Google Scholar] [CrossRef] [Green Version]
  9. Jade, S. Estimates of plate velocity and crustal deformation in the Indian subcontinent using GPS. Curr. Sci. 2004, 86, 1443–1448. [Google Scholar]
  10. Bettinelli, P.; Avouac, J.P.; Flouzat, M. Plate Motion of India and Interseismic Strain in the Nepal Himalaya from GPS and DORIS Measurements. J. Geod. 2006, 80, 567–589. [Google Scholar] [CrossRef]
  11. Bastos, L.; Bos, M.; Fernandes, R.M. Deformation and Tectonics: Contribution of GPS Measurements to Plate Tectonics-Overview and Recent Developments; Sciences of Geodesy—I: Advances and Future Directions; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  12. Elliott, J.; Freymueller, J.T.; Larsen, C.F. Active tectonics of the St. Elias orogen, Alaska, observed with GPS measurements. J. Geophys. Res. Solid Earth 2013, 118, 5625–5642. [Google Scholar] [CrossRef]
  13. Van Gelder, B.H.; Aardoom, L. SLR Network Designs in View of Reliable Detection of Plate Kinematics in the East Mediterranean; Reports of the Department of Geodesy; Delft University of Technology: Delft, The Netherlands, 1982; pp. 1–24. [Google Scholar]
  14. Christodoulidis, D.C.; Smith, D.E.; Kolenkiewicz, R.; Klosko, S.M.; Torrence, M.H.; Dunn, P.J. Observing tectonic plate motions and deformations from satellite laser ranging. J. Geophys. Res. 1985, 90, 9249–9263. [Google Scholar] [CrossRef]
  15. Smith, D.E.; Kolenkiewicz, R.; Dunn, P.J.; Robbins, J.W.; Torrens, M.H.; Klosko, S.M.; Williamson, R.G.; Pavlis, E.C.; Douglas, N.B.; Fricke, S.K. Tectonic motion and deformation from Satellite Laser Ranging to LAGEOS. J. Geoph. Res. 1990, 95, 22013–22041. [Google Scholar] [CrossRef]
  16. Sengoku, A. A plate motion study using Ajisai SLR data. Earth Planets Space 1998, 50, 611–627. [Google Scholar] [CrossRef] [Green Version]
  17. Alothman, A.O.; Schillak, S. Recent Results for the Arabian Plate Motion Using Satellite Laser Ranging Observations of Riyadh SLR Station to LAGEOS-1 and LAGEOS-2 Satellites. Arab. J. Sci. Eng. 2014, 39, 217–226. [Google Scholar] [CrossRef]
  18. Atanasova, M.; Georgiev, I.; Chapanov, Y. Global Tectonic Plate Motions from Slr Data Processing; University of Architecture Civil Engineering and Geodesy: Sofia, Bulgaria, 2018; Volume 51, pp. 109–114. [Google Scholar]
  19. Crétaux, J.F.; Soudarin, L.; Cazenave, A.; Bouillé, F. Present-day tectonic plate motions and crustal deformations from the DORIS space system. J. Geophys. Res. Solid Earth 1998, 103, 30167–30181. [Google Scholar] [CrossRef]
  20. Soudarin, L.; Crétaux, J.F. A model of present-day tectonic plate motions from 12 years of DORIS measurements. J. Geod. 2006, 80, 609–624. [Google Scholar] [CrossRef]
  21. Kroger, P.M.; Lyzenga, G.A.; Walace, K.S.; Davidson, J.M. Tectonic motion in the western United States inferred from Very Long Baseline Interferometry measurements 1980–1986. J. Geoph. Res. 1987, 92, 14151–14153. [Google Scholar] [CrossRef]
  22. Sato, K. Tectonic plate motion and deformation inferred from very long baseline interferometry. Tectonophysics 1993, 220, 69–87. [Google Scholar] [CrossRef]
  23. Haas, R.; Gueguen, E.; Scherneck, H.G.; Nothnagel, A.; Campbell, J. Crustal motion results derived from observations in the European geodetic VLBI network. Earth Planets Space 2000, 52, 759–764. [Google Scholar] [CrossRef] [Green Version]
  24. Haas, R.; Nothnagel, A.; Campbell, J.; Gueguen, E. Recent crustal movements observed with the European VLBI network: Geodetic analysis and results. J. Geodyn. 2003, 35, 391–414. [Google Scholar] [CrossRef]
  25. Krásná, H.; Ros, T.; Pavetich, P.; Böhm, J.; Nilsson, T.; Schuh, H. Investigation of crustal motion in Europe by analysing the European VLBI sessions. Acta Geod. Geophys. 2013, 48, 389–404. [Google Scholar] [CrossRef]
  26. Kurt, O. Monitoring Movements of Tectonic Plates by Analyzing VLBI Data via QGIS. In Proceedings of the Scientific Congress of the Turkish National Union of Geodesy and Geophsysics (TNUGG-SC), Izmir, Turkey, 30 May–2 June 2018; pp. 198–201. [Google Scholar]
  27. Argus, D.F.; Gordon, R.G.; Heflin, M.B.; Eanes, R.J.; Willis, P.; Peltier, W.R.; Owen, S.E. The angular velocities of the plates and the velocity of Earth’s centre from space geodesy. Geophys. J. Int. 2010, 180, 913–960. [Google Scholar] [CrossRef]
  28. Le Pichon, X.; Angelier, J.; Sibuet, J.C. Plate boundaries and extensional tectonics. Tectonophysics 1982, 81, 239–256. [Google Scholar] [CrossRef]
  29. Sen, S. Earth: The Extraordinary Planet; Avant Publishing Services Pvt. Ltd.: Noida, Uttar Pradesh, India, 2018. [Google Scholar]
  30. Drewes, H. A geodetic approach for the recovery of global kinematic plate parameters. Bull. Géodésique 1982, 56, 70–79. [Google Scholar] [CrossRef]
  31. Drewes, H. Global Plate Motion Parameters Derived from Actual Space Geodetic Observations. In Proceedings of the International Association of Geodesy Symposia 101, Global and Regional Geodynamics, Edinburgh, UK, 3–5 August 1989. [Google Scholar]
  32. Altamimi, Z.; Rebischung, P.; Métivier, L.; Collilieux, X. ITRF2014: A new release of the International Terrestrial Reference Frame modeling nonlinear station motions. J. Geophys. Res. 2016, 121, 6109–6131. [Google Scholar] [CrossRef] [Green Version]
  33. Kraszewska, K.; Jagoda, M.; Rutkowska, M. Tectonic Plate Parameters Estimated in the International Terrestrial Reference Frame ITRF2008 Based on SLR Stations. Acta Geophys. 2016, 64, 1495–1512. [Google Scholar] [CrossRef] [Green Version]
  34. Jagoda, M.; Rutkowska, M.; Suchocki, C.; Katzer, J. Determination of the tectonic plates motion parameters based on SLR, DORIS and VLBI stations positions. J. Appl. Geod. 2020, 14, 121–131. [Google Scholar] [CrossRef]
  35. Kraszewska, K.; Jagoda, M.; Rutkowska, M. Tectonic plates parameters estimated in International Terrestrial Reference Frame ITRF2008 based on DORIS stations. Acta Geophys. 2018, 66, 509–521. [Google Scholar] [CrossRef] [Green Version]
  36. Jagoda, M.; Rutkowska, M. Use of VLBI measurement technique to determination of the tectonic plates motion parameters. Metrol. Meas. Syst. 2020, 27, 151–165. [Google Scholar]
  37. Schonenberg, R.; Neugebauer, J. Einfuhrung in Die Geologie Europas; Rombach Verlag: Freiburg im Breisgau, Germany, 1994. [Google Scholar]
  38. Parfenow, L.M.; Badarch, G.; Berzin, N.A.; Hwang, D.H.; Khanchuk, A.I.; Kuzmin, M.I.; Nokleberg, W.J.; Obolenskiy, A.A.; Ogasawara, M.; Prokopiev, A.V.; et al. Introduction to Region al Geology, Tectonics and Metallogenesis on North Coast Asia; United States Geological Survey: Reston, VA, USA, 2007.
  39. Rey, P.; Burg, J.P.; Casey, M. The Scandinavian Caledonides and Their Relationship to the Varisean Belt; Geological Society of London: London, UK, 1997. [Google Scholar]
  40. Moores, E.M.; Fairbridge, R.W. (Eds.) Encyclopedia of European and Asian Regional Geology; Chapman & Hall: London, UK; New York, NY, USA; Tokyo, Japan; Melbourne, Australia; Madras, India, 1997. [Google Scholar]
  41. Yin, A.; Harrison, M. (Eds.) The Tectonic Evolution of Asia; Cambridge University Press: Cambridge, UK, 1996. [Google Scholar]
  42. Molnar, P.; Tapponier, P. Cenozoic tectonics of Asia: Effects of continental collision. Science 1975, 189, 419–426. [Google Scholar] [CrossRef] [PubMed]
  43. Auboutin, J. Some aspects of the tectonics of subduction zones. Tectonophysics 1989, 160, 1–21. [Google Scholar] [CrossRef]
  44. McCann, T. The Geology of Central Europe: Precambrian and Paleozoik; Geological Society of London: London, UK, 2008. [Google Scholar]
  45. Drewes, H. The actual plate kinematic and crustal deformation model APKIM2005 as basis for a non-rotating ITRF. In Proceedings of the International Association of Geodesy Symposia 134, Berlin/Heidelberg, Germany, 9–14 October 2009; pp. 95–101. [Google Scholar]
Figure 1. The location of the GNSS stations used in the four scenarios of the calculations.
Figure 1. The location of the GNSS stations used in the four scenarios of the calculations.
Sensors 20 06065 g001
Figure 2. The location and motion of the GNSS stations that are not consistent with the Eurasian plate motion.
Figure 2. The location and motion of the GNSS stations that are not consistent with the Eurasian plate motion.
Sensors 20 06065 g002
Figure 3. The location and motion of the GNSS stations that are not consistent with the Eurasian plate motion; area of Section 1 (Figure 2).
Figure 3. The location and motion of the GNSS stations that are not consistent with the Eurasian plate motion; area of Section 1 (Figure 2).
Sensors 20 06065 g003
Figure 4. Results of the sequential solutions of the Φ parameter for four calculation scenarios.
Figure 4. Results of the sequential solutions of the Φ parameter for four calculation scenarios.
Sensors 20 06065 g004
Figure 5. Results of the sequential solutions of the Λ parameter for four calculation scenarios.
Figure 5. Results of the sequential solutions of the Λ parameter for four calculation scenarios.
Sensors 20 06065 g005
Figure 6. Results of the sequential solutions of the ω parameter for four calculation scenarios.
Figure 6. Results of the sequential solutions of the ω parameter for four calculation scenarios.
Sensors 20 06065 g006
Table 1. Eurasian plate motion parameters obtained after adding individually selected stations (given in Appendix A) that are not consistent with the Eurasian plate motion or are not located on the Eurasian plate to calculation scenario 1.
Table 1. Eurasian plate motion parameters obtained after adding individually selected stations (given in Appendix A) that are not consistent with the Eurasian plate motion or are not located on the Eurasian plate to calculation scenario 1.
No.Name of the Station
and Geodetic Position
Φ [°] Λ [°] ω [°/Ma]
30 (scenario 1 - final solution) 54.87 ± 0.47 261.02 ± 0.50 0.2585 ± 0.0026
3130 (from scenario 1) + Bilibino (BILI),
B = 67°, L = 166°
53.06 ± 0.47 259.60 ± 0.51 0.2500 ± 0.0026
3130 (from scenario 1) + Magadan (MAGO), B = 59°, L = 150°53.66 ± 0.56 259.71 ± 0.57 0.2505 ± 0.0030
3130 (from scenario 1) + Petropavlovsk (PETP), B = 52°, L = 158°50.41 ± 1.48 254.19 ± 1.34 0.2335 ± 0.0066
3130 (from scenario 1) + Teheran (TEHN),
B = 35°, L = 51°
55.11 ± 0.83 261.40 ± 0.77 0.2551 ± 0.0038
3130 (from scenario 1) + Kunming (KUNM), B = 24°, L = 102°54.03 ± 0.65 260.04 ± 0.68 0.2502 ± 0.0036
3130 (from scenario 1) + Lhasa (LHAS),
B = 29°, L = 91°
57.63 ± 0.98 265.05 ± 1.170.2666 ± 0.0062
3130 (from scenario 1) + Mizusawa (MIZU), B = 38°, L = 141°50.20 ± 1.84 254.67 ± 1.670.2334 ± 0.0082
3130 (from scenario 1) + Koganei (KGNI),
B = 35°, L = 139°
52.23 ± 1.32 257.04 ± 1.39 0.2330 ± 0.0069
3130 (from scenario 1) + Usuda (USUD),
B = 35°, L = 138°
52.30 ± 1.27 257.10 ± 1.24 0.2336 ± 0.0062
3130 (from scenario 1) + Tsukuba (TSKB),
B = 35°, L=140°
51.39 ± 1.44256.25 ± 1.37 0.2333 ± 0.0068
3130 (from scenario 1) + Abashiri (P202),
B = 44°, L = 144°
52.42 ± 1.26 256.35 ± 1.22 0.2340 ± 0.0057
Table 2. Comparison of the Eurasian plate motion parameters estimated in this paper with solutions derived by Drewes [45] and Larson et al. [8].
Table 2. Comparison of the Eurasian plate motion parameters estimated in this paper with solutions derived by Drewes [45] and Larson et al. [8].
No.Solution Φ [°] Λ [°] ω [°/Ma]
1APKIM2005 model [45]53.4 ± 0.4264.3 ± 0.50.259 ± 0.001
2Larson et al. [8]56.3257.20.26
3Solution given in this paper (based on 120 GNSS station positions)54.81 ± 0.37261.04 ± 0.480.2585 ± 0.0025
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jagoda, M.; Rutkowska, M. An Analysis of the Eurasian Tectonic Plate Motion Parameters Based on GNSS Stations Positions in ITRF2014. Sensors 2020, 20, 6065. https://doi.org/10.3390/s20216065

AMA Style

Jagoda M, Rutkowska M. An Analysis of the Eurasian Tectonic Plate Motion Parameters Based on GNSS Stations Positions in ITRF2014. Sensors. 2020; 20(21):6065. https://doi.org/10.3390/s20216065

Chicago/Turabian Style

Jagoda, Marcin, and Miłosława Rutkowska. 2020. "An Analysis of the Eurasian Tectonic Plate Motion Parameters Based on GNSS Stations Positions in ITRF2014" Sensors 20, no. 21: 6065. https://doi.org/10.3390/s20216065

APA Style

Jagoda, M., & Rutkowska, M. (2020). An Analysis of the Eurasian Tectonic Plate Motion Parameters Based on GNSS Stations Positions in ITRF2014. Sensors, 20(21), 6065. https://doi.org/10.3390/s20216065

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop