不同海潮模型对中国大陆构造环境监测网络GNSS数据解算的影响

选取2018年中国大陆构造环境监测网络GNSS基准站及中国周边IGS站数据,基于CSR4、FES2004、GOT00、NAO99b等海潮负荷改正模型,采用单一变量的解算方式,分析不同海潮模型对我国GNSS基准站基线、三维坐标、时间序列等造成的影响。结果表明：利用不同海潮模型解算基线精度,水平方向相差不大于1 mm,垂向相差不大于3 mm;各测站NEU坐标精度均在毫米级以下,以耿马站时间序列为例,与其他模型相比,基于FES2004模型解算的各向偏差趋势波动较弱,其年序列图的精度及连续性、稳定性更佳。简言之,利用FES2004海潮模型解算GNSS站点位置,效果较好。     

Influence of the modified global ocean tide model with local tides of East and South China Seas on load gravity in China and its neighbor area

Introduction With the development of science and technology, the accuracy of gravity measurement is im-proved. The gravity observation with FG5 gravimeter has reached accuracy of μ magnitude. And the high accurate superconducting gravimeter can detect the tiny signal of 10?2 μ magnitude in frequency domain. With the high-accuracy gravity observation on Earth′s surface, the Earth′s tidal parameters can be determined precisely. And the observations can also be used to invert the struc-ture…     

LAGEOS Sensitivity to Ocean Tides

Satellite Laser Ranging (SLR) to LAGEOS has a remarkable contribution to high-precise geodesy and geodynamics through deriving and validating various global geophysical models. This paper validates ocean tide models based on the analysis of satellite altimetry data, coastal tide gauges, and hydrodynamic data, i.e., CSR3.0, TOPEX4.0, CSR4.0A, FES2004, GOT00.2, and the CSRC Schwiderski model. LAGEOS orbits and SLR observation residuals from solutions based on different ocean tide models are compared and examined. It is found that LAGEOS orbits are sensitive to tidal waves larger than 5 mm. The analysis of the aliasing periods of LAGEOS orbits and tidal waves reveals that, in particular, the tidal constituent S₂ is not well established in the recent ocean tide models. Some of the models introduce spurious peaks to empirical orbit parameters, which can be associated with S₂, S_a, and K₂ tidal constituents, and, as a consequence, can be propagated to fundamental parameters derived from LAGEOS observations.     

基于近岸海岛GPS数据反演海洋分潮负荷影响

本文探索了海潮负荷特征值的反演方法.基于中国近岸海岛GPS站(平潭与闸坡)观测数据,采用FFT方法提取了海潮负荷特征值,分析了反演特征值与全球海潮模型FES2004、NAO.99b和GOT4.7计算出的相关特征值之间的差异,评估了反演海潮分潮频率项的精度,并利用平潭站得到的海潮负荷特征值对附近的三沙站进行海潮改正以评价反演效果.结果表明:(1)频谱分析可精确提取4个半日分潮和4个全日分潮负荷的频率信息;与已知频率相比,P1分潮的反演频率误差为1.4%,其他7种分潮负荷反演频率误差均小于1%.(2)两个海岛分属不同的潮波入侵通道,反演分潮振幅和初相存在差异,但反演分潮频率几乎一致,间接证实它们属于同一潮波系统,也表明其潮差有别.(3)反演振幅与三种全球模型具有较好的一致性;其中S2、O1、P1、Q1四个分潮在水平方向互差为1~2 mm,高程方向上的互差均小于3 mm;K2、K1、M2、N2在水平方向振幅互差多数小于2 mm,个别差异高达4 mm,高程方向互差多为5~6 mm,个别超过10 mm.(4)反演得到海潮负荷改正模型相对于3种全球模型在三沙站的改正效果略佳,间接表明反演结果有效、可靠.(5)动态PPP结果中虽然存在多种误差,其时间序列仍可分离并提取海潮负荷的影响.     

Transcontinental tidal transect: European Atlantic coast-Southern Siberia-Russian Pacific coast

The paper presents results of measurements with digital tidal LaCoste-Romberg gravimeters on the European Atlantic coast-Southern Siberia-Russian Pacific coast transect in 1995–2005. The transect includes four West European (Chizé, Ménesplet, Mordelles, and Wikle), two South Siberian (Klyuchi and Talaya), and two Far Eastern (Zabakalskoe and Yuzhno-Sakhalinsk) stations. Gravimetric measurements at the Talaya station (SW Baikal rift zone) are supplemented by long-term laser extensometer observations. The position of the stations within the rectangle (45°–55°N, 0.4°–142°E) allows one to assess existing tidal strain models (WD93 and DDW99) and various ocean tide models (SCW80, CSR3, FES95, ORI96, CSR4, FES02, GOT00, NAO99, and TPX06). Data of intracontinental stations (with a small ocean effect at distances of 2000–3000 km) agree well with the DDW99 tidal strain model (with regard to the mantle viscosity). The uncertainty of digital tidal gravity measurements is 0.25%. Results of laser extensometer measurements are at the same accuracy level. Then, the Love and Shida numbers calculated at midlatitudes of the intracontinental zone of Eurasia from combined data are h = 0.6077 ± 0.0008, k = 0.3014 ± 0.0001, and l = 0.0839 ± 0.0001. The analysis of results of Pacific and Atlantic stations located at distances of 30–300 km from the ocean showed that the FES02, CSR4, GOT00, NAO99, and TPX06 ocean tide models are preferable.     

Parameters of the fluid core resonance inferred from superconducting gravimeter data

We have estimated the parameters of fluid core resonance (FCR) due to the nearly diurnal free wobble of the Earth's core based on the superconducting gravimeter (SG) data obtained at the following four observation sites; Esashi and Matsushiro in Japan, Canberra in Australia and Membach in Belgium. By fitting the tidal admittances normalized with the O1 wave at each site to a model of the damped harmonic oscillator, we obtained values of 429.66 ± 1.43 sidereal days, 9350–10,835, −4.828E−4 ± 3.4E−6, −3.0E−5 ± 4.5E−6 for the eigenperiod, the Q-value and the real and imaginary parts of the resonance strength, respectively. Our values obtained from only using the gravity data are very consistent with those inferred from the VLBI nutation data. Our study strongly indicates that the systematic difference between two estimations from the gravity and the nutation in particular for the Q-value, which has been shown in previous works, is mainly caused by the inaccurate correction for the ocean tide effects. The error in the ocean tide correction is discussed based on the comparison among four global ocean tide models; Schwiderski model (1980), NAO.99b (Matsumoto et al., 2000), CSR4.0 (Eanes and Bettadpur, 1994) and GOT99.2b (Ray, 1999).     

Loading effect of a self-consistent equilibrium ocean pole tide on the gravimetric parameters of the gravity pole tides at superconducting gravimeter stations

The gravimetric parameters of the gravity pole tide are the amplitude factor δ, which is the ratio of gravity variations induced by polar motion for a real Earth to variations computed for a rigid one, and the phase difference κ between the observed and the rigid gravity pole tide. They can be estimated from the records of superconducting gravimeters (SGs). However, they are affected by the loading effect of the ocean pole tide. Recent results from TOPEX/Poseidon (TP) altimeter confirm that the ocean pole tide has a self-consistent equilibrium response. Accordingly, we calculate the gravity loading effects as well as their influence on the gravimetric parameters of gravity pole tide at all the 26 SG stations in the world on the assumption of a self-consistent equilibrium ocean pole tide model. The gravity loading effect is evaluated between 1 January 1997 and 31 December 2006. Numerical results show that the amplitude of the gravity loading effect reaches 10⁻⁹ m s⁻², which is larger than the accuracy (10⁻¹⁰ m s⁻²) of a SG. The gravimetric factor δ is 1% larger at all SG stations. Then, the contribution of a self-consistent ocean pole tide to the pole tide gravimetric parameters cannot be ignored as it exceeds the current accuracy of the estimation of the pole tide gravity factors. For the nine stations studied in Ducarme et al. [Ducarme, B., Venedikov, A.P., Arnoso, J., et al., 2006. Global analysis of the GGP superconducting gravimeters network for the estimation of the pole tide gravimetric amplitude factor. J. Geodyn. 41, 334–344.], the mean of the modeled tidal factors δ_m = 1.1813 agrees very well with the result of a global analysis δ_CH = 1.1816 ± 0.0047 in that paper. On the other hand, the modeled phase difference κ_m varies from −0.273° to 0.351°. Comparing to the two main periods of the gravity pole tide, annual period and Chandler period, κ_m is too small to be considered. Therefore, The computed time difference κ_L induced by a self-consistent ocean pole tide produces a negligible effect on κ_m. It confirms the results of Ducarme et al., 2006, where no convincing time difference was found in the SG records.     

Earth tides observed by gravity and GPS in southeastern Alaska

We analyzed gravity data obtained in Juneau and global positioning system (GPS) data obtained from three PBO sites in southeastern Alaska (SE-AK), which are part of a US research facility called ‘EarthScope’, and we compared the obtained tidal amplitudes and phases with those estimated from the predicted tides including both effects of the body tide and ocean tide. Global tide models predict the ocean tides in this region of complex coastline and bathymetry. To improve the accuracy of prediction, we developed a regional ocean tide model in SE-AK.Our comparison results suggest: (1) by taking into account the ocean tide effect, the amplitude differences between the observation and the predicted body tide is remarkably reduced for both the gravity and displacement (e.g. for the M₂ constituent, 8.5–0.3 μGal, and 2.4–0.1 cm at the AB50 GPS site in Juneau in terms of the vector sum of three components of the north–south, east–west and up–down), even though the ocean tide loading is large in SE-AK. (2) We have confirmed the precise point positioning (PPP) method, which was used to extract the tidal signals from the original GPS time series, works well to recover the tidal signals. Although the GPS analysis results still contain noise due to the atmosphere and multipath, we may conclude that the GPS observation surely detects the tidal signals with the sub-centimeter accuracy or better for some of the tidal constituents. (3) In order to increase the accuracy of the tidal prediction in SE-AK, it is indispensable to improve the regional ocean tide model developed in this study, especially for the phase.     

Earth tides observed by gravity and GPS in southeastern Alaska

We analyzed gravity data obtained in Juneau and global positioning system (GPS) data obtained from three PBO sites in southeastern Alaska (SE-AK), which are part of a US research facility called ‘EarthScope’, and we compared the obtained tidal amplitudes and phases with those estimated from the predicted tides including both effects of the body tide and ocean tide. Global tide models predict the ocean tides in this region of complex coastline and bathymetry. To improve the accuracy of prediction, we developed a regional ocean tide model in SE-AK.Our comparison results suggest: (1) by taking into account the ocean tide effect, the amplitude differences between the observation and the predicted body tide is remarkably reduced for both the gravity and displacement (e.g. for the M₂ constituent, 8.5–0.3 μGal, and 2.4–0.1 cm at the AB50 GPS site in Juneau in terms of the vector sum of three components of the north–south, east–west and up–down), even though the ocean tide loading is large in SE-AK. (2) We have confirmed the precise point positioning (PPP) method, which was used to extract the tidal signals from the original GPS time series, works well to recover the tidal signals. Although the GPS analysis results still contain noise due to the atmosphere and multipath, we may conclude that the GPS observation surely detects the tidal signals with the sub-centimeter accuracy or better for some of the tidal constituents. (3) In order to increase the accuracy of the tidal prediction in SE-AK, it is indispensable to improve the regional ocean tide model developed in this study, especially for the phase.     

Assessment of GPS data for meteorological applications over Africa: Study of error sources and analysis of positioning accuracy

The aim of this study is to assess the availability and quality of data from the International GNSS Service (IGS) Global Positioning System (GPS) network in Africa, especially for retrieving zenith tropospheric delay (ZTD), from which precipitable water vapour (PWV) can be derived, in view of application to the African Monsoon Multidisciplinary Analysis (AMMA) project. Three major error sources for the GPS data analysis evaluating PWV in Africa are the accuracy of the satellite orbits, the correction for the radio delay induced by the ionosphere and the vertical site displacements due to ocean loading. The first part of this study examines these error sources and the validity of GPS data for meteorological applications in Africa in dedicated analyses spanning the year 2001. These analyses were performed using the IGS precise orbits. Weak degradation of baseline precision with increasing baseline lengths suggests that the average orbital error is not limiting the GPS analysis in Africa. The impact of the ionosphere has been evaluated during a maximum of solar activity in 2001. The loss of L₂ data has actually been observed. It amounts to 2% on average for 2001, with maxima of 8% during magnetic storm events. A slight decrease in formal accuracy of ZTD seems to be related to the loss of L₂ data at the end of the day. This indicates that scintillation effects are present in the GPS observations but however are not a major limitation. The impact of ocean loading is found to be significant on ZTD estimates (up to ±2 mm in equivalent PWV). The use of a proper ocean loading model eliminates this effect.The second aspect of this study concerns the IGS analysis quality for the African stations. The accuracy has been assessed through position dispersion between individual solutions and the most recent version of the IGS combined solution IGb00, and residuals from the transformation of the IGS combined solution into the International Terrestrial Reference Frame 2005. The positioning performance of the IGS analysis is consistent with an accuracy in ZTD of ±6 mm (±1 mm in PWV), as requested for meteorological applications such as planned in AMMA.     

On the comparison of tidal gravity parameters with tidal models in central Europe

Two accurately calibrated superconducting gravimeters (SGs) provide high quality tidal gravity records in three central European stations: C025 in Vienna and at Conrad observatory (A) and OSG050 in Pecný (CZ). To correct the tidal gravity factors from ocean loading effects we compared the load vectors from different ocean tides models (OTMs) computed with different software: OLFG/OLMP by the Free Ocean Tides Loading Provider (FLP), ICET and NLOADF. Even with the recent OTMs the mass conservation is critical but the methods used to correct the mass imbalance agree within 0.1 nm/s². Although the different software agrees, FLP probably provides more accurate computations as this software has been optimised. For our final computation we used the mean load vector computed by FLP for 8 OTMs (CSR4, NAO99, GOT00, TPX07, FES04, DTU10, EOT11a and HAMTIDE). The corrected tidal factors of the 3 stations agree better than 0.04% in amplitude and 0.02° in phase. Considering the weighted mean of the three stations we get for O1 δ_c = 1.1535 ± 0.0001, for K1 δ_c = 1.1352 ± 0.0003 and for M2 δ_c = 1.1621 ± 0.0003. These values confirm previous ones obtained with 16 European stations. The theoretical body tides model DDW99/NH provides the best agreement for M2 (1.1620) and MATH01/NH for O1 (1.1540) and K1 (1.1350). The largest discrepancy is for O1 (0.05%). The corrected phase α_c does not differ significantly from zero except for K1 and S2. The calibrations of the two SG's are consistent within 0.025% and agree with Strasbourg results within 0.05%.     

Determination of the long period tidal waves in the GGP superconducting gravity data

Several results about the long period (LP) tidal waves are obtained by the analysis of series of superconducting gravity data, provided by the Global Geodynamics Project (GGP). The most important result is the determination of a single group called LPMF, composed by all LP tides but representing accurately the parameters of the M_f wave. As the LP tidal generating potential is vanishing at latitudes ±35°15′52″ we cannot determine accurate tidal amplitude factors for the stations located between ±40° and ±30°. However, it is still possible to obtain tidal residual vectors and compare them with oceanic tidal loading computations. For 15 stations the NAO99 oceanic model is giving a coherent picture. For nine stations with M_f amplitude larger than 3 μgal (1 μgal = 10 nm s⁻²) a global analysis is obtained by introducing the loading effect of the ocean directly in the observation equations. The mean amplitude factor obtained for LPMF is larger than expected from the models and there is a significant phase lag, showing the imperfection of the tidal oceanic models for M_f. Other new result is the first separate estimation of the parameters of the LP tides, generated by the tidal potential of third degree, dominated by a Lunar declinational monthly wave, called here 3M_md. Due to their small amplitudes (under 1 μgal) these waves are practically hidden by the noise. Nevertheless, the quality of the data and the flexibility of the VAV analysis method [Venedikov, A.P., Arnoso, J., Vieira, R., 2001. Program VAV/2000 for tidal analysis of unequally spaced data with irregular drift and colored noise. J. Geodetic Soc. Jpn. 47 (1), 281–286; Venedikov, A.P., Arnoso, J., Vieira, R., 2003. VAV: a program for tidal data processing. Comput. Geosci. 29, 487–502.] allow getting significant results, in agreement with the theory of the Earth deformation by the tidal potential of third degree.     

Loading effects in Metsähovi from the atmosphere and the Baltic Sea

Loading by atmosphere and by the Baltic Sea cause gravity change at Metsähovi, located 15 km from the open sea. Gravity is changed by both the Newtonian attraction of the loading mass and by the crustal deformation. We have performed loading calculations using appropriate Green's function for both gravity and deformation, for both atmospheric and Baltic loading. The loading by atmosphere has been computed using a detailed surface pressure field from high resolution limited area model (HIRLAM) for north Europe up to 10° distances. Baltic Sea level is modelled using tide gauge records. Calculations show that 1 m of uniform layer of water corresponds to 31 nm s⁻² in gravity and −11 mm in height. Modelled loading is compared with observations of the superconducting gravimeter T020 for years 1994–2002. The combination of HIRLAM and a tide gauge record decreases RMS of gravity residuals by 14% compared to single admittance in air pressure corrections without sea level data. Regression of gravity residuals on the tide gauge record at Helsinki (at 30 km distance) gives a gravity effect of 26 nm s⁻² m⁻¹ for Baltic loading.The gravity station is co-located with a permanent GPS station. We have also associated the loading effects of the atmosphere and of the Baltic Sea with temporal height variations. The range of modelled vertical motion due to air pressure was 46 mm and that due to sea level 18 mm. The total range was 38 mm. The effects of the Baltic Sea and of the atmosphere partly cancel each other, since at longer periods the inverse barometer assumption is valid. Regression of the modelled height on local air pressure gives −0.37 mm hPa⁻¹, corresponding approximately to width 6° for pressure system.We have tested the models using one year of daily GPS data. Multilinear regression on local air pressure and sea level in Helsinki gives the coefficient −0.34 mm hPa⁻¹ for pressure, and −11 mm m⁻¹ for sea level. These match model values. Loading by air pressure and Baltic Sea explains nearly 40% of the variance of daily GPS height solutions.     

The lunar tides in the Antarctic mesosphere and lower thermosphere

Horizontal winds in the mesosphere and lower thermosphere over the Antarctic have been measured by a meteor radar at Rothera (67.5°S, 68.0°W) and MF radar at Davis (68.6°S, 78.0°E). Data from Rothera recorded over a 20-month interval in 2005–2006 and data from Davis recorded over the 13-year interval 1994–2006 are examined to investigate the monthly mean behaviour of the lunar semidiurnal tide. Both data sets show a clear signal of the 12.42-h lunar semidiurnal (M₂) tide. The amplitude reaches values as large as 8 m s⁻¹. The vertical wavelengths of the tide vary seasonally from 10 to 65 km. Comparisons of the phase of the tide measured over the two sites reveals that it does not purely consist of a migrating wavenumber 2 mode. This suggests that other, non-migrating, modes are likely to be present.     

Modification of fault slip models of the Mw 9.0 Tohoku Earthquake by far field GPS observations

GPS data from Crustal Movement Observation Network of China (CMONOC) are used to derive far-field co-seismic displacements induced by the M_w 9.0 Tohoku Earthquake. Significant horizontal displacements about 30 mm, 10 mm, and 20 mm were caused by this large event in northeast China, north China, and on the Korean peninsula respectively. Vectors of relatively large horizontal displacements with dominant east components pointed to the epicenter of this earthquake. The east components show an exponential decay with the longitude, which is characteristic of the decay of the co-seismic horizontal displacements associated with earthquakes of thrust rupture. The exponential fit of the east components shows that the influence of the co-seismic displacements can be detected by GPS at a distance of about 3200 km from the epicenter of the earthquake. By considering the capability of the far field displacements for constraining the inversion of the fault slip model of the earthquake, we use spherically stratified Earth models to simulate the co-seismic displacements induced by this event. Using computations and comparisons, we discuss the effects of parameters of layered Earth models on the results of dislocation modeling. Comparisons of the modeled and observed displacements show that far field GPS observations are effective for constraining the fault slip model. The far field horizontal displacements observed by GPS are used to modify the slips and seismic moments of fault slip models. The result of this work is applicable as a reference for other researchers to study seismic source rupture and crustal deformation.     

Seasonal variation of the <Emphasis Type="Italic">M</Emphasis><Subscript>2</Subscript> tide

The seasonal cycle of the main lunar tidal constituent M ₂ is studied globally by an analysis of a high-resolution ocean circulation and tide model (STORMTIDE) simulation, of 19 years of satellite altimeter data, and of multiyear tide-gauge records. The barotropic seasonal tidal variability is dominant in coastal and polar regions with relative changes of the tidal amplitude of 5–10 %. A comparison with the observations shows that the ocean circulation and tide model captures the seasonal pattern of the M ₂ tide reasonably well. There are two main processes leading to the seasonal variability in the barotropic tide: First, seasonal changes in stratification on the continental shelf affect the vertical profile of eddy viscosity and, in turn, the vertical current profile. Second, the frictional effect between sea-ice and the surface ocean layer leads to seasonally varying tidal transport. We estimate from the model simulation that the M ₂ tidal energy dissipation at the sea surface varies seasonally in the Arctic (ocean regions north of 60°N) between 2 and 34 GW, whereas in the Southern Ocean, it varies between 0.5 and 2 GW. The M ₂ internal tide is mainly affected by stratification, and the induced modified phase speed of the internal waves leads to amplitude differences in the surface tide signal of 0.005–0.0150 m. The seasonal signals of the M ₂ surface tide are large compared to the accuracy demands of satellite altimetry and gravity observations and emphasize the importance to consider seasonal tidal variability in the correction processes of satellite data.     

Global surface mass from a new combination of GRACE,modelled OBP and reprocessed GPS data

Weekly surface loading variations are estimated from a joint least squares inversion of load-induced GPS site displacements, GRACE gravimetry and simulated ocean bottom pressure (OBP) from the finite element sea-ice ocean model (FESOM).In this study, we directly use normal equations derived from reprocessed GPS observations, where station and satellite positions are estimated simultaneously. The OBP weight of the model in the inversion is based on a new error model, obtained from 2 FESOM runs forced with different atmospheric data sets.Our findings indicate that the geocenter motion derived from the inversion is smooth, with non-seasonal RMS values of 1.4, 0.9 and 1.9 mm for the X, Y and Z directions, respectively. The absolute magnitude of the seasonal geocenter motion varies annually between 2 and 4.5 mm. Important hydrological regions such as the Amazon, Australia, South-East Asia and Europe are mostly affected by the geocenter motion, with magnitudes of up to 2 cm, when expressed in equivalent water height.The chosen solar radiation pressure model, used in the GPS processing, has only a marginal effect on the joint inversion results. Using the empirical CODE model slightly increases the annual amplitude of the Z component of the geocenter by 0.8 mm. However, in case of a GPS-only inversion, notable larger differences are found for the annual amplitude and phase estimates when applying the older physical ROCK models. Regardless of the used radiation pressure model the GPS network still exhibits maximum radial expansions in the order of 3 mm (0.45 ppb in terms of scale), which are most likely caused by remaining GPS technique errors.In an additional experiment, we have used the joint inversion solution as a background loading model in the GPS normal equations. The reduced time series, compared to those without a priori loading model, show a consistent decrease in RMS. In terms of the annual height component, 151 of the 189 stations show a reduction of at least 10% in seasonal amplitude.On the ocean floor, we find a positive overall correlation (0.51) of the inversion solution with time series from globally distributed independent bottom pressure recorders.Even after removing a seasonal fit we still find a correlation of 0.45. Furthermore, the geocenter motion has a significant effect on ocean bottom pressure as neglecting it causes the correlation to drop to 0.42.     

近海潮汐效应对测站位移的负荷影响

利用TPXO6全球海潮模型和中国东海和南海潮汐资料研究了近海潮汐效应对我国沿海部分台站负荷位移的影响.结果表明,近海效应对位移的水平和垂直分量的影响均达到了mm的量级,其中在厦门站的影响最大,水平和垂直分量的振幅分别接近3 mm和15 mm,因此近海潮汐效应是沿海地区GPS动力学应用中应该考虑的负荷效应中的一个因素.另外,各台站负荷位移水平分量随时间变化的规律是不一样的,最大值对于不同的台站出现在不同的方位,并且对于垂直分量,在上海站和在其他沿海站的影响存在着明显的反相关系.本文结果可为GPS资料处理提供有益参考.     

由热带气旋MATTHEW引起的非潮汐负荷效应影响

热带气旋能在短期内造成海水和大气质量的重新分布,使得近海地表受力发生变化,进而产生非潮汐负荷形变,对现今高精度大地测量的影响已不容忽视.为了保证空间大地测量结果的精度和稳定性,热带气旋引起的地表形变必须进行有效的估计.因此本文联合NOS、GLOSS验潮站数据与海潮模型,通过获取非潮汐残余量分析了热带气旋“MATTHEW”引起的风暴潮.利用ECCO海洋环流模型、ERAin大气再分析模型、HUGO-m海洋动力学模型,分别估计了“MATTHEW”引起非潮汐海洋负荷、非潮汐大气负荷、动力学响应下非潮汐海洋负荷对地表位移的影响,结果表明热带气旋影响下的广大区域地表都不同程度受到非潮汐负荷的作用,最大位移分别达到-9.13 mm、3.31 mm、-6.11 mm,并且加入动力学响应的非潮汐海洋负荷要普遍大于IB(Inverted Barometer)响应下的结果.在对比不同位置站点所受负荷差异时,发现“大陆站”非潮汐海洋负荷形变普遍大于“岛屿站”,而“岛屿站”更易受非潮汐大气负荷的影响.     

武汉超导重力仪观测最新结果和海潮模型研究

利用武汉台站GWR_C032超导重力仪观测资料,在对原始数据进行有效预处理的基础上作调和分析,获得反映地球内部介质特征的重力潮汐参数.基于卫星测高技术和有限元方法同时考虑验潮站数据作约束条件获得的多个全球海潮模型,利用负荷理论和数值褶积积分技术计算了重力负荷,对周日和半日频段内的重力潮汐参数实施负荷改正,提出了“负荷改正有效性”概念,研究了全球海潮模型适应性.数值结果说明,海潮改正的有效性高达91%（O1,NAO99）和92%（M2,ORI96）.基于11个海潮模型对主波（O1,K1,M2和S2）的负荷改正说明平均有效性为（86%,70%,73%和84%）,振幅因子与理论模型间的差异分别从（212%,155%,116%和080%）降到（031%,039%,034%和008%）,同时还说明利用NAO99和ORI96全球海潮模型能获得比其他模型更佳的负荷改正效果.文章还利用国际地球动力学计划网络其他7个台站的超导重力仪观测研究了全球海潮模型的适定性问题,结果说明不同模型中不同潮波具有明显的区域特点,早期构制的SCW80全球海潮模型仍可作为大地测量研究中的重要参考模型.