Tidal strain observations carried out in the Sopronbánfalva Geodynamic Observatory (SGO) in Hungary are used to test thirteen different ocean tide loading models for diurnal and semi-diurnal tidal harmonic constituents O1, K1 and M2. Strain data with one minute sampling rate were corrected for temperature and barometric pressure and decimated to one hour sampling rate. Strain data, corrected in this way, were subjected to correction for ocean load. In the case of the diurnal tidal constituents O1 and K1 the measured amplitude factors of nearly 0.5 became close to the theoretical as a result of the correction, while in the case of the M2 semi diurnal wave the measured amplitude factor of almost 1 hardly changed due to correction. It was only found a negligible difference between the individual global ocean tide loading models mainly due to using different Earth models and Green functions. The effect of the diurnal (O1 and K1) and the semidiurnal (M2) ocean tide loading components is in the same order of magnitude at the SGO. The large residual vectors after the correction suggest that local effects need further investigation.
Earth tide, Extensometer, Tidal parameters, Ocean tide loading, Strain
Ocean tide surface loading causes both radial and tangential displacements of the Earth surface and changes of gravity. This latter comprises effects resulting from radial displacement in the Earth's gravity field, the internal redistribution of mass, and the direct gravitational attraction of the tidal water mass [1-5].
Various authors have been dealing with testing and comparison of recent global ocean tide models on the basis of gravimetric [6-17], displacement (GPS) [18-25], extensometric [15,16,26,27] and tilt measurements [28-30].
In this paper, the efficiency of ocean loading corrections are compared using thirteen different ocean tide loading models for the diurnal and semi-diurnal tidal harmonic constituents O1, K1 and M2 on the basis of strain data observed in the SGO in 2017.
Extensometric measurements were carried out in the SGO which is located on the Hungarian-Austrian border in the Sopron Mountains. The coordinates of the observatory are: Latitude 47°40'55'' N; longitude 16°33'32" E; the altitude is 280 m a.s.l. The yearly mean temperature in the gallery is 10.4 ℃ and the yearly and daily temperature variations are less than 0.5 ℃ and 0.05 ℃, respectively. The instrument is a 22 m long quartz-tube extensometer with capacitive transducer. The azimuth of the extensometer is 116°, and its scale factor is 2.093 ± 0.032 nm mV-1. Observation site and construction of the extensometer and its calibration are described in detail by Mentes [31,32].
Strain, temperature and barometric pressure data recorded with a sampling rate of one minute were used for correction by the T-soft program [33]. Data series were despiked and ungapped. The long-term constituent of the strain and temperature data were approximated by fitting a polynomial of 9th order to the raw data series and were subtracted from strain and temperature data, respectively. Theoretical tide was subtracted from the remaining strain data and then strain data were corrected for the temperature and barometric pressure by simple linear regression method and after the correction the theoretical tide was added back. During the correction procedure time lags between strain and temperature and barometric pressure data were taking into consideration. Then the data were low-pass filtered and decimated to one hour sampling and processed by ETERNA 3.40 Earth tide data processing program package [34] using the Wahr-Dehant Earth model [35,36] and the HW95 tidal potential catalogue [37].
For the ocean tide load prediction the SPOTL routines [38] and the ocean load provider service [39] were used. The name of the SPOTL codes are: gr.mmmmmm.www.pnn,c[e|m], The mmmmmm string denotes the Earth model. The SPOTL uses three different Earth models, which are denoted namely the Gutenberg-Bullen Model A average Earth (gbaver) and two extreme models of the Earth's crust and mantle structure [40], one oceanic (gbocen) and one continental shield (gbcont). www denotes: Who computed the function: the source code is wef for W.E. Farrell [1], who computed and tabulated all Green functions. pnn denotes: 01 coarse grid; 02 fine grid (about 20 m for the innermost range) running to much closer to the center (about 100 m). The ce suffix denotes a reference frame with the center of mass of solid Earth and the cm a reference frame with the common center of mass of the load and the Earth.
For ocean load correction of strain data ten ocean tide models have been used with SPOTL, supplemented by the local model osu.mediterranean.2011: EOT11a [41], HAMTIDE11a [42], OSU.TPXO72atlas, OSU.TPXO72, TPX070 [43], DTU10 [44], CSR4.0 [45], FES2004 [46], FES95.2.1 [47], SCHW1 [48] and three other models were chosen from the Free Ocean Tide Loading Provider created by Scherneck and Bos (http://holt.oso.chalmers.se/loading/): FES2012 [49], FES2014b [50,51], GOT00.2 [52].
Without external forces the common centre of mass of oceans and the solid Earth will remain fixed in space. Since the ocean tides cause water mass displacements, its centre of mass moves periodically and it is compensated by an opposite motion of the centre of mass of the solid Earth. Accordingly, in these three cases the ocean load was calculated relative to the fixed common mass center of the ocean and the solid Earth and the moving center of mass of the solid Earth. In both cases the calculations were carried out using elastic [1] and visco-elastic Earth model STW105 [53].
The ocean loading correction (see Figure 1) was carried out according to Neumeyer, et al. [54]. The L amplitude and λ phase of the ocean tide load vectors were determined from the above-mentioned different ocean tide loading models and were subtracted from the observed strain tidal vectors (A→m(Am,αm)) in case of the O1, K1 and M2 tidal waves to obtain the corrected tidal strain vectors (Ac−→(Ac,αc)). The remaining residual X depends on the accuracy of the instrument calibration, the local effects, such as the cavity effect, the inaccurately corrected temperature, air pressure and on the accuracy of the ocean tide model.
Results of the tidal evaluation of the measured strain data in case of the tidal constituents O1, K1 and M2 without ocean tide loading corrections are shown in Table 1. In the SGO the obtained amplitude factors for the diurnal waves (O1, K1,) are about 0.5, half of the theoretical value, while for the semidiurnal wave M2 it is about 1. Table 2, Table 3, Table 4 and Table 5 show the amplitudes (L) and phases (λ) of the ocean tide load, the amplitudes (Ac) and phases (αc) of strain corrected for the ocean tide load, the amplitudes (X) and phases (ξ) of the residues and the corrected amplitude factors (ηc) in case of different Earth and global ocean tide loading models.
Table 2 shows the O1 tidal constituents corrected for ocean tide loading with SPOTL in case of different Earth models and 10 global ocean tide models. The uncorrected amplitude factor is 0.5323 (see Table 1), while the corrected amplitude factors (average: 1.025 ± 0.001) are somewhat higher than one in the case of solid Earth model with the center of mass (ce suffix denotes the reference frame with the center of mass of solid Earth) and somewhat lower (average: 0.928 ± 0.003) in the case of the Earth model with the common center of mass of the load and the Earth (cm suffix denotes the reference frame with the common center of mass of the load and the Earth). Similar results were obtained for K1 (see Table 3) for which the uncorrected amplitude factor was 0.5283. The corrected average amplitude factors are 1.026 ± 0.001 (ce) and 0.966 ± 0.002 (cm). The measured amplitude factor for M2 is 1.0036 (see Table 1), while the corrected average values (see Table 4) are 1.059 ± 0.004 (ce) and 1.039 ± 0.007 (cm). In both cases the corrected values are higher than the measured value, but similarly to O1 and K1 the corrected values are higher in case of reference frame with the center of Earth (ce) than in the case of common center of load mass and the center of Earth.
Table 5 shows the results of correction of the O1, K1, M2 waves for ocean tide loading derived from three ocean loading models (FS2012, FS214b and GOT00.2) with different Earth models calculated by the ocean load provider service. In case of the elastic and visco-elastic Earth models almost identical corrected amplitude factors (ηc) were obtained for O1, K1 and M2 waves but these were somewhat higher when motion correction was applied in case of the diurnal waves (O1 and K1), while for the semidiurnal wave M2 the opposite values were obtained.
Comparing the residues in Table 2, Table 3, Table 4 and Table 5, it can be seen that the amplitudes (X) of the residual vectors changed only slightly compared to the measured ones (B), while the phase angles decreased significantly. Large residues remaining after the correction suggest that local effects (e.g. cavity) affect the measurement site, which requires further study. The amplitudes (L) of the ocean load vectors are in the same order of magnitude for the diurnal (O1, K1) and semidiurnal (M2) waves at the site of the SGO.
Strain measurement was used to test thirteen global ocean tide loading models. Tidal parameters corrected for ocean tide loading were calculated. All models provided virtually the same result. In the case of the diurnal tidal constituents O1 and K1 the measured amplitude factors of nearly 0.5 became close to 1 as a result of the correction, while in the case of the M2 semi diurnal wave, the measured amplitude factor of almost 1 hardly changed due to correction. It was only found a negligible difference between the individual tide loading models mainly due to the use of different Earth models, and Green functions. The effect of the diurnal (O1 and K1) and the semidiurnal (M2) ocean tide loading components is in the same order of magnitude at the SGO. The large residual vectors after the correction suggest that local effects need further investigation.
This work was funded by the Hungarian National Research Fund (OTKA) under project K 109060. Special thanks to Tibor Molnár for his careful maintenance of the instruments.