Search

Atmospheric ties are theoretically affected by the height differences between antennas at the same site and the meteorological conditions. However, there is often a discrepancy between the expected zenith delay differences and those estimated from geodetic analysis, potentially degrading a combined solution employing atmospheric ties. In order to investigate the possible effects on GNSS atmospheric delay, this study set up an experiment of four co-located GNSS stations of the same type, both antenna and receiver. Specific height differences for each antenna w.r.t the reference antenna are given. One antenna was equipped with a radome at the same height and type as a antenna close to the ground. In addition, a meteorological sensor was used for meteorological data recording. The results show that tropospheric ties from the analytical equation based on meteorological data from GPT3, Numerical Weather Model, and in-situ measurements, and ray-traced tropospheric ties, reduced the bias of zenith delay roughly by 72 %. However, the in-situ tropospheric ties yield the best precision in this study. These results demonstrate, that the instrument effects on GNSS zenith delays were mitigated by using the same instrument. In contrast, the radome causes unexpected bias of GNSS zenith delays in this study. Additionally, multipath effects at low-elevation observations degraded the tropospheric east gradients.

The next-generation, broadband geodetic very long baseline interferometry system, named VGOS, is developing its global network, and VGOS networks with a small size of 3--7 stations have already made broadband observations from 2017 to 2019. We made quality assessments for two kinds of observables in the 21 VGOS sessions currently available: group delay and differential total electron content ($\delta$TEC). Our study reveals that the random measurement noise of VGOS group delays is at the level of less than 2 ps (1 ps = 10$^{-12}$ s), while the contributions from systematic error sources, mainly source structure related, are at the level of 20 ps. Due to the significant improvement in measurement noise, source structure effects with relatively small magnitudes that are not overwhelming in the S/X VLBI system, for instance 10 ps, are clearly visible in VGOS observations. Another critical error source in VGOS observations is discrete delay jumps, for instance, a systematic offset of about 310 ps or integer multiples of that. The predominant causative factor is found to be related to source structure. The measurement noise level of $\delta$TEC observables is about 0.07 TECU, but the systematic effects are five times larger than that. A strong correlation between group delay and $\delta$TEC observables is discovered with a trend of 40 ps/TECU for observations with large structure effects; there is a second trend in the range 60 ps/TECU to 70 ps/TECU when the measurement noise is dominant., Comment: 25 pages, 11 figures, and accepted for publication in Journal of Geodesy

A new realization of the International Celestial Reference Frame (ICRF) is presented based on the work achieved by a working group of the International Astronomical Union (IAU) mandated for this purpose. This new realization, referred to as ICRF3, is based on nearly 40 years of data acquired by very long baseline interferometry. The ICRF3 includes positions at 8.4 GHz for 4536 sources, supplemented with positions at 24 GHz for 824 sources and at 32 GHz for 678 sources, for a total of 4588 sources. A subset of 303 sources among these, uniformly distributed on the sky, are identified as "defining sources" and as such serve to define the axes of the frame. Source positions are reported for epoch 2015.0 and must be propagated for observations at other epochs for the most accurate needs, accounting for the acceleration toward the Galactic center, which results in a dipolar proper motion field of amplitude 0.0058 milliarcsecond/yr (mas/yr). The frame shows a median positional uncertainty of about 0.1 mas in right ascension and 0.2 mas in declination, with a noise floor of 0.03 mas in the individual source coordinates. A subset of 500 sources is found to have extremely accurate positions at 8.4 GHz, in the range of 0.03 to 0.06 mas. Comparing ICRF3 with the Gaia Celestial Reference Frame 2 in the optical domain, there is no evidence for deformations larger than 0.03 mas between the two frames. Significant positional offsets between the three ICRF3 frequencies are detected for about 5% of the sources. Moreover, a notable fraction (22%) of the sources shows optical and radio positions that are significantly offset. There are indications that these positional offsets may be the manifestation of extended source structures. This third realization of the ICRF was adopted by the IAU at its 30th General Assembly in August 2018 and replaced the previous realization, ICRF2, on January 1, 2019., Comment: 28 pages, 22 figures

International audience; There has been growing interest in recent years in the use of homogeneously reprocessed ground-based GNSS, VLBI, and DORIS measurements for climate applications. Existing datasets are reviewed and the sensitivity of tropospheric estimates to the processing details is discussed. The uncertainty in the derived IWV estimates and linear trends is around 1 kg m−2 RMS and ± 0.3 kg m−2 per decade, respectively. Standardized methods for ZTD outlier detection and IWV conversion are proposed. The homogeneity of final time series is limited however by changes in the stations equipment and environment. Various homogenization algorithms have been evaluated based on a synthetic benchmark dataset. The uncertainty of trends estimated from the homogenized times series is estimated to ±0.5 kg m−2 per decade. Reprocessed GNSS IWV data are analysed along with satellites data, reanalyses and global and regional climate model simulations. A selection of global and regional reprocessed GNSS datasets and ERA-interim reanalysis are made available through the GOP-TropDB tropospheric database and online service. A new tropo SINEX format, providing new features and simplifications, was developed and it is going to be adopted by all the IAG services.

Very Long Baseline Interferometry (VLBI) plays an indispensable role in the realization of global terrestrial and celestial reference frames and in the determination of the full set of the Earth Orientation Parameters (EOP). The main goal of this research is to assess the quality of the VLBI observations based on the sensitivity and robustness criteria. Sensitivity is defined as the minimum displacement value that can be detected in coordinate unknowns. Robustness describes the deformation strength induced by the maximum undetectable errors with the internal reliability analysis. The location of a VLBI station and the total weights of the observations at the station are most important for the sensitivity analysis. Furthermore, the total observation number of a radio source and the quality of the observations are important for the sensitivity levels of the radio sources. According to the robustness analysis of station coordinates, the worst robustness values are caused by atmospheric delay effects with high temporal and spatial variability. During CONT14, it is determined that FORTLEZA, WESTFORD, and TSUKUB32 have robustness values changing between 0.8 and 1.3 mm, which are significantly worse in comparison to the other stations. The radio sources 0506-612, NRAO150, and 3C345 have worse sensitivity levels compared to other radio sources. It can be concluded that the sensitivity and robustness analysis are reliable measures to obtain high accuracy VLBI solutions.

Terrestrial reference frames (TRF), such as the ITRF2008, are primary products of geodesy. In this paper, we present TRF solutions based on Kalman filtering of very long baseline interferometry (VLBI) data, for which we estimate steady station coordinates over more than 30 years that are updated for every single VLBI session. By applying different levels of process noise, non-linear signals, such as seasonal and seismic effects, are taken into account. The corresponding stochastic model is derived site-dependent from geophysical loading deformation time series and is adapted during periods of post-seismic deformations. Our results demonstrate that the choice of stochastic process has a much smaller impact on the coordinate time series and velocities than the overall noise level. If process noise is applied, tests with and without additionally estimating seasonal signals indicate no difference between the resulting coordinate time series for periods when observational data are available. In a comparison with epoch reference frames, the Kalman filter solutions provide better short-term stability. Furthermore, we find out that the Kalman filter solutions are of similar quality when compared to a consistent least-squares solution, however, with the enhanced attribute of being easier to update as, for instance, in a post-earthquake period.

Local ties (LTs) at co-located sites are currently used to combine different single-technique solutions to determine global terrestrial reference frames (TRFs). We assess by simulations the impact of different LT standard deviations, biased LTs and a selection of LTs on the datum realization of global TRFs. The simulations are based on Global Positioning System (GPS), Satellite Laser Ranging, and Very Long Baseline Interferometry (VLBI) observations covering the time span 2008–2014. We find that LT standard deviations of 1 cm and better yield differences in the TRF-defining parameters below 1 mm. VLBI is most affected by altering the LT standard deviations, especially in the translations since VLBI is inherently not sensitive to the origin of the TRF. Altering the standard deviations of the LTs applied in ITRF2005, ITRF2008, ITRF2014 results in small differences reaching a maximum of 0.6 mm at the VLBI stations. Simulating technique-wise biased LT stations shows the largest differences in the TRF-defining parameters of more than 2 mm, if all GPS LT stations are biased by 1 cm, proving that GPS plays the major role in the connection of the three techniques. Simulating single biased LT stations by 1 cm in either the north, east, or height component indicates small differences of less than 0.8 mm in the TRF-defining parameters, the largest differences result at LT stations located on the southern hemisphere. The selection of LTs demonstrates that the southern hemisphere LTs are very important, especially for the realization of the scale.

The tropospheric horizontal gradients with high spatiotemporal resolutions provide important information to describe the azimuthally asymmetric delays and significantly increase the ability of ground-based GNSS (Global Navigation Satellite Systems) within the field of meteorological studies, like the nowcasting of severe rainfall events. The recent rapid development of multi-GNSS constellations has potential to provide such high-resolution gradients with a significant degree of accuracy. In this study, we develop a multi-GNSS process for the precise retrieval of high-resolution tropospheric gradients. The tropospheric gradients with different temporal resolutions, retrieved from both single-system and multi-GNSS solutions, are validated using independent numerical weather models (NWM) data and water vapor radiometer (WVR) observations. The benefits of multi-GNSS processing for the retrieval of tropospheric gradients, as well as for the improvement of precise positioning, are demonstrated. The multi-GNSS high-resolution gradients agree well with those derived from the NWM and WVR, especially for the fast-changing peaks, which are mostly associated with synoptic fronts. The multi-GNSS gradients behave in a much more stable manner than the single-system estimates, especially in cases of high temporal resolution, benefiting from the increased number of observed satellites and improved observation geometry. The high-resolution multi-GNSS gradients show higher correlation with the NWM and WVR gradients than the low-resolution gradients. Furthermore, the precision of station positions can also be noticeably improved by multi-GNSS fusion, and enhanced results can be achieved if the high-resolution gradient estimation is performed, instead of the commonly used daily gradient estimation in the multi-GNSS data processing.

The revitalized Russian GLONASS system provides new potential for real-time retrieval of zenith tropospheric delays (ZTD) and precipitable water vapor (PWV) in order to support time-critical meteorological applications such as nowcasting or severe weather event monitoring. In this study, we develop a method of real-time ZTD/PWV retrieval based on GLONASS and/or GPS observations. The performance of ZTD and PWV derived from GLONASS data using real-time precise point positioning (PPP) technique is carefully investigated and evaluated. The potential of combining GLONASS and GPS data for ZTD/PWV retrieving is assessed as well. The GLONASS and GPS observations of about half a year for 80 globally distributed stations from the IGS (International GNSS Service) network are processed. The results show that the real-time GLONASS ZTD series agree quite well with the GPS ZTD series in general: the RMS of ZTD differences is about 8 mm (about 1.2 mm in PWV). Furthermore, for an inter-technique validation, the real-time ZTD estimated from GLONASS-only, GPS-only, and the GPS/GLONASS combined solutions are compared with those derived from very long baseline interferometry (VLBI) at colocated GNSS/VLBI stations. The comparison shows that GLONASS can contribute to real-time meteorological applications, with almost the same accuracy as GPS. More accurate and reliable water vapor values, about 1.5---2.3 mm in PWV, can be achieved when GLONASS observations are combined with the GPS ones in the real-time PPP data processing. The comparison with radiosonde data further confirms the performance of GLONASS-derived real-time PWV and the benefit of adding GLONASS to stand-alone GPS processing.

This paper presents a summary of the activities of the Bureau of Products and Standards (BPS) to support IAG’s Global Geodetic Observing System (GGOS) in its goal to provide observations and consistent geodetic products needed to monitor, map, and understand changes in the Earth’s shape, rotation, and mass distribution. As a key activity the BPS has compiled an inventory of the standards and conventions currently adopted and used by the IAG and its components for the processing of geometric and gravimetric observations as the basis for the generation of IAG products. The outcome of the BPS inventory concerning numerical standards and the product-based review is summarized and recommendations for future improvements are provided. Finally, an overview about the ongoing and planned activities of the BPS is given.

We performed Monte Carlo simulations of very-long-baseline interferometry (VLBI) observations of Earth-orbiting satellites incorporating co-located space-geodetic instruments in order to study how well the VLBI frame and the spacecraft frame can be tied using such measurements. We simulated observations of spacecraft by VLBI observations, time-of-flight (TOF) measurements using a time-encoded signal in the spacecraft transmission, similar in concept to precise point positioning, and differential VLBI (D-VLBI) observations using angularly nearby quasar calibrators to compare their relative performance. We used the proposed European Geodetic Reference Antenna in Space (E-GRASP) mission as an initial test case for our software. We found that the standard VLBI technique is limited, in part, by the present lack of knowledge of the absolute offset of VLBI time to Coordinated Universal Time at the level of microseconds. TOF measurements are better able to overcome this problem and provide frame ties with uncertainties in translation and scale nearly a factor of three smaller than those yielded from VLBI measurements. If the absolute time offset issue can be resolved by external means, the VLBI results can be significantly improved and can come close to providing 1 mm accuracy in the frame tie parameters. D-VLBI observations with optimum performance assumptions provide roughly a factor of two higher uncertainties for the E-GRASP orbit. We additionally simulated how station and spacecraft position offsets affect the frame tie performance.

We investigate the possibilities to use data from water vapor radiometers (WVR) to calibrate the wet tropospheric delays in geodetic Very Long Baseline Interferometry (VLBI) observations. We test three methods: (1) direct calibration using WVR measurements aquired in the directions of the VLBI observations, (2) estimating zenith wet delays and gradients from the WVR data and use these to correct the VLBI data, and (3) including the WVR measurements as additional observations in the VLBI data analysis. Furthermore, in all cases we model the WVR calibration errors in the data analysis. We test the three methods using data from the continuous VLBI campaigns CONT02–CONT14. We find clear improvements when applying methods 1 and 3 for CONT05 campaign, however, the results are degraded for the other campaigns.

Very Long Baseline Interferometry (VLBI) is the only space geodetic technique capable of measuring all the Earth orientation parameters (EOP) accurately and simultaneously. Modeling the Earth's rotational motion in space within the stringent consistency goals of the Global Geodetic Observing System (GGOS) makes VLBI observations essential for constraining the rotation theories. However, the inaccuracy of early VLBI data and the outdated products could cause non-compliance with these goals. In this paper, we perform a global VLBI analysis of sessions with different processing settings to determine a new set of empirical corrections to the precession offsets and rates, and to the amplitudes of a wide set of terms included in the IAU 2006/2000A precession-nutation theory. We discuss the results in terms of consistency, systematic errors, and physics of the Earth. We find that the largest improvements w.r.t. the values from IAU 2006/2000A precession-nutation theory are associated with the longest periods (e.g., 18.6-yr nutation). A statistical analysis of the residuals shows that the provided corrections attain an error reduction at the level of 15 μas. Additionally, including a Free Core Nutation (FCN) model into a priori Celestial Pole Offsets (CPOs) provides the lowest Weighted Root Mean Square (WRMS) of residuals. We show that the CPO estimates are quite insensitive to TRF choice, but slightly sensitive to the a priori EOP and the inclusion of different VLBI sessions. Finally, the remaining residuals reveal two apparent retrograde signals with periods of nearly 2069 and 1034 days. This work was funded and realized in the framework of the project AYA2016-79775-P (AEI/FEDER, UE) and APOSTD/2026/079.

In this study, we assess the impact of two combination strategies, namely local ties (LT) and global ties (GT), on the datum realization of Global Terrestrial Reference Frames in view of the Global Geodetic Observing System requiring 1 mm-accuracy. Simulated Very Long Baseline Interferometry (VLBI) and Satellite Laser Ranging (SLR) data over a 7 year time span was used. The LT results show that the geodetic datum can be best transferred if the precision of the LT is at least 1 mm. Investigating different numbers of LT, the lack of co-located sites on the southern hemisphere is evidenced by differences of 9 mm in translation and rotation compared to the solution using all available LT. For the GT, the combination applying all Earth rotation parameters (ERP), such as pole coordinates and UT1-UTC, indicates that the rotation around the Z axis cannot be adequately transferred from VLBI to SLR within the combination. Applying exclusively the pole coordinates as GT, we show that the datum can be transferred with mm-accuracy within the combination. Furthermore, adding artificial stations in Tahiti and Nigeria to the current VLBI network results in an improvement in station positions by 13 and 12%, respectively, and in ERP by 17 and 11%, respectively. Extending to every day VLBI observations leads to 65% better ERP estimates compared to usual twice-weekly VLBI observations.

Precise positioning with the current Chinese BeiDou Navigation Satellite System is proven to be of comparable accuracy to the Global Positioning System, which is at centimeter level for the horizontal components and sub-decimeter level for the vertical component. But the BeiDou precise point positioning (PPP) shows its limitation in requiring a relatively long convergence time. In this study, we develop a numerical weather model (NWM) augmented PPP processing algorithm to improve BeiDou precise positioning. Tropospheric delay parameters, i.e., zenith delays, mapping functions, and horizontal delay gradients, derived from short-range forecasts from the Global Forecast System of the National Centers for Environmental Prediction (NCEP) are applied into BeiDou real-time PPP. Observational data from stations that are capable of tracking the BeiDou constellation from the International GNSS Service (IGS) Multi-GNSS Experiments network are processed, with the introduced NWM-augmented PPP and the standard PPP processing. The accuracy of tropospheric delays derived from NCEP is assessed against with the IGS final tropospheric delay products. The positioning results show that an improvement in convergence time up to 60.0 and 66.7% for the east and vertical components, respectively, can be achieved with the NWM-augmented PPP solution compared to the standard PPP solutions, while only slight improvement in the solution convergence can be found for the north component. A positioning accuracy of 5.7 and 5.9 cm for the east component is achieved with the standard PPP that estimates gradients and the one that estimates no gradients, respectively, in comparison to 3.5 cm of the NWM-augmented PPP, showing an improvement of 38.6 and 40.1%. Compared to the accuracy of 3.7 and 4.1 cm for the north component derived from the two standard PPP solutions, the one of the NWM-augmented PPP solution is improved to 2.0 cm, by about 45.9 and 51.2%. The positioning accuracy for the up component improves from 11.4 and 13.2 cm with the two standard PPP solutions to 8.0 cm with the NWM-augmented PPP solution, an improvement of 29.8 and 39.4%, respectively.

Precise transformation between the celestial reference frames (CRF) and terrestrial reference frames (TRF) is needed for many purposes in Earth and space sciences. According to the Global Geodetic Observing System (GGOS) recommendations, the accuracy of positions and stability of reference frames should reach 1 mm and 0.1 mm year $$^{-1}$$ , and thus, the Earth Orientation Parameters (EOP) should be estimated with similar accuracy. Different realizations of TRFs, based on the combination of solutions from four different space geodetic techniques, and CRFs, based on a single technique only (VLBI, Very Long Baseline Interferometry), might cause a slow degradation of the consistency among EOP, CRFs, and TRFs (e.g., because of differences in geometry, orientation and scale) and a misalignment of the current conventional EOP series, IERS 08 C04. We empirically assess the consistency among the conventional reference frames and EOP by analyzing the record of VLBI sessions since 1990 with varied settings to reflect the impact of changing frames or other processing strategies on the EOP estimates. Our tests show that the EOP estimates are insensitive to CRF changes, but sensitive to TRF variations and unmodeled geophysical signals at the GGOS level. The differences between the conventional IERS 08 C04 and other EOP series computed with distinct TRF settings exhibit biases and even non-negligible trends in the cases where no differential rotations should appear, e.g., a drift of about 20 $$\upmu $$ as year $$^{-1 }$$ in $$y_{\mathrm{pol }}$$ when the VLBI-only frame VTRF2008 is used. Likewise, different strategies on station position modeling originate scatters larger than 150 $$\upmu $$ as in the terrestrial pole coordinates.

This paper introduces the reader to our Kalman filter developed for geodetic VLBI (very long baseline interferometry) data analysis. The focus lies on the EOP (Earth Orientation Parameter) determination based on the Continuous VLBI Campaign 2014 (CONT14) data, but also earlier CONT campaigns are analyzed. For validation and comparison purposes we use EOP determined with the classical LSM (least squares method) estimated from the same VLBI data set as the Kalman solution with a daily resolution. To gain higher resolved EOP from LSM we run solutions which yield hourly estimates for polar motion and dUT1 = Universal Time (UT1) – Coordinated Universal Time (UTC). As an external validation data set we use a GPS (Global Positioning System) solution providing hourly polar motion results. Further, we describe our approach for determining the noise driving the Kalman filter. It has to be chosen carefully, since it can lead to a significant degradation of the results. We illustrate this issue in context with the de-correlation of polar motion and nutation. Finally, we find that the agreement with respect to GPS can be improved by up to 50% using our filter compared to the LSM approach, reaching a similar precision than the GPS solution. Especially the power of erroneous high-frequency signals can be reduced dramatically, opening up new possibilities for high-frequency EOP studies and investigations of the models involved in VLBI data analysis. We prove that the Kalman filter is more than on par with the classical least squares method and that it is a valuable alternative, especially on the advent of the VLBI2010 Global Observing System and within the GGOS frame work.

Closure quantities measured by very long baseline interferometry (VLBI) observations are independent of instrumental and propagation instabilities and antenna gain factors, but are sensitive to source structure. A new method is proposed to calculate a structure index based on the median values of closure quantities rather than the brightness distribution of a source. The results are comparable to structure indices based on imaging observations at other epochs and demonstrate the flexibility of deriving structure indices from exactly the same observations as used for geodetic analysis and without imaging analysis. A three-component model for the structure of source 3C371 is developed by model-fitting closure phases. It provides a real case of tracing how the structure effect identified by closure phases in the same observations as the delay observables affects the geodetic analysis, and investigating which geodetic parameters are corrupted to what extent by the structure effect. Using the resulting structure correction based on the three-component model of source 3C371, two solutions, with and without correcting the structure effect, are made. With corrections, the overall rms of this source is reduced by 1 ps, and the impacts of the structure effect introduced by this single source are up to 1.4 mm on station positions and up to 4.4 microarcseconds on Earth orientation parameters. This study is considered as a starting point for handling the source structure effect on geodetic VLBI from geodetic sessions themselves., 5 figures, 15 pages, accepted by Journal of Geodesy at 19 Dec., 2016

The Free Core Nutation (FCN) is a free mode of the Earth's rotation caused by the different material characteristics of the Earth's core and mantle. This causes the rotational axes of those layers to slightly diverge from each other, resulting in a wobble of the Earth's rotation axis comparable to nutations. In this paper we focus on estimating empirical FCN models using the observed nutations derived from the VLBI sessions between 1993 and 2013. Assuming a fixed value for the oscillation period, the time-variable amplitudes and phases are estimated by means of multiple sliding window analyses. The effects of using different a priori Earth Rotation Parameters (ERP) in the derivation of models are also addressed. The optimal choice of the fundamental parameters of the model, namely the window width and step-size of its shift, is searched by performing a thorough experimental analysis using real data. The former analyses lead to the derivation of a model with a temporal resolution higher than the one used in the models currently available, with a sliding window reduced to 400 days and a day-by-day shift. It is shown that this new model increases the accuracy of the modeling of the observed Earth's rotation. Besides, empirical models determined from USNO Finals as a priori ERP present a slightly lower Weighted Root Mean Square (WRMS) of residuals than IERS 08 C04 along the whole period of VLBI observations, according to our computations. The model is also validated through comparisons with other recognized models. The level of agreement among them is satisfactory. Let us remark that our estimates give rise to the lowest residuals and seem to reproduce the FCN signal in more detail. This work has been partly supported by two Spanish Projects from CGL2010-12153-E and AYA2010-22039-C02-01.

The atmospheric parameters, zenith delays and gradients, obtained by the DORIS, GPS, VLBI, and numerical weather models, ECMWF and NCEP, are compared at five DORIS co-located sites during the 15 days of the CONT14 campaign from 2014-05-06 until 2014-05-20. Further examined are two different solutions of GPS, VLBI and NCEP: for GPS, a network solution comparable to the TIGA reprocessing analysis strategy and a precise point positioning solution, for VLBI, a least squares and a Kalman filtered and smoothed solution, and for NCEP two spatial resolutions, 0.5° and 1.0°, are tested. The different positions of the antenna reference points at co-location sites affect the atmospheric parameters and have to be considered prior to the comparison. We assess and discuss these differences, tropospheric ties, by comparing ray-traced atmospheric parameters obtained at the positions of the various antenna reference points. While ray-traced ZHD and ZWD at the co-located antennas significantly differ, the ray-traced gradients show only very small differences. Weather events can introduce larger disagreement between atmospheric parameters obtained at co-location sites. The various weather model solutions in general agree very well in providing tropospheric ties. The atmospheric parameters are compared using statistical methods, such as the mean difference and standard deviations with repect to a weighted mean value. While GPS and VLBI atmospheric parameters agree very well in general, the DORIS observations are in several cases not dense enough to achieve a comparable level of agreement. The estimated zenith delays from DORIS, however, are competitive with the other space geodetic techniques. If the DORIS observation geometry is insufficient for the estimation of an atmospheric gradient, less than three satellites observed during the definition interval, the DORIS atmospheric parameters degrade and show small quasi-periodic variations that correlate with the number of observations and in particular with the number of satellites. An increase in the DORIS constellation concerning more satellites and in general more observations is very likely to significantly improve the quality of DORIS derived atmospheric parameters. For the first time we tested a 6 h sampling of the DORIS gradients. Where the observations are sufficiently dense, the increased sampling results in an improvement of the agreement of the DORIS gradients with the other solutions.

The CONT14 campaign with state-of-the-art VLBI data has observed the source 0642+449 with about one thousand observables each day during a continuous observing period of fifteen days, providing tens of thousands of closure delays---the sum of the delays around a closed loop of baselines. The closure delay is independent of the instrumental and propagation delays and provides valuable additional information about the source structure. We demonstrate the use of this new "observable" for the determination of the structure in the radio source 0642+449. This source, as one of the defining sources in the second realization of the International Celestial Reference Frame (ICRF2), is found to have two point-like components with a relative position offset of -426 microarcseconds in right ascension and -66 microarcseconds in declination. The two components are almost equally bright with a flux-density ratio of 0.92. The standard deviation of closure delays for source 0642+449 was reduced from 139 ps to 90 ps by using this two-component model. Closure delays larger than one nanosecond are found to be related to the source structure, demonstrating that structure effects for a source with this simple structure could be up to tens of nanoseconds. The method described in this paper does not rely on a priori source structure information, such as knowledge of source structure determined from direct (Fourier) imaging of the same observations or observations at other epochs. We anticipate our study to be a starting point for more effective determination of the structure effect in VLBI observations., 14 pages, 13 figures, Accepted by Astronomical Journal on 12 Jul, 2016

The real-time estimation of polar motion (PM) is needed for the navigation of Earth satellite and interplanetary spacecraft. However, it is impossible to have real-time information due to the complexity of the measurement model and data processing. Various prediction methods have been developed. However, the accuracy of PM prediction is still not satisfactory even for a few days in the future. Therefore, new techniques or a combination of the existing methods need to be investigated for improving the accuracy of the predicted PM. There is a well-introduced method called Copula, and we want to combine it with singular spectrum analysis (SSA) method for PM prediction. In this study, first, we model the predominant trend of PM time series using SSA. Then, the difference between PM time series and its SSA estimation is modeled using Copula-based analysis. Multiple sets of PM predictions which range between 1 and 365 days have been performed based on an IERS 08 C04 time series to assess the capability of our hybrid model. Our results illustrate that the proposed method can efficiently predict PM. The improvement in PM prediction accuracy up to 365 days in the future is found to be around 40% on average and up to 65 and 46% in terms of success rate for the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hbox{PM}}_{x}$$\end{document}PMx and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hbox{PM}}_{y}$$\end{document}PMy, respectively.

The understanding of forced temporal variations in celestial pole motion (CPM) could bring us significantly closer to meeting the accuracy goals pursued by the Global Geodetic Observing System (GGOS) of the International Association of Geodesy (IAG), i.e., 1 mm accuracy and 0.1 mm/year stability on global scales in terms of the Earth orientation parameters. Besides astronomical forcing, CPM excitation depends on the processes in the fluid core and the core–mantle boundary. The same processes are responsible for the variations in the geomagnetic field (GMF). Several investigations were conducted during the last decade to find a possible interconnection of GMF changes with the length of day (LOD) variations. However, less attention was paid to the interdependence of the GMF changes and the CPM variations. This study uses the celestial pole offsets (CPO) time series obtained from very long baseline interferometry (VLBI) observations and data such as spherical harmonic coefficients, geomagnetic jerk, and magnetic field dipole moment from a state-of-the-art geomagnetic field model to explore the correlation between them. In this study, we use wavelet coherence analysis to compute the correspondence between the two non-stationary time series in the time–frequency domain. Our preliminary results reveal interesting common features in the CPM and GMF variations, which show the potential to improve the understanding of the GMF’s contribution to the Earth’s rotation. Special attention is given to the corresponding signal between FCN and GMF and potential time lags between geomagnetic jerks and rotational variations. J.M.F was partially supported by Spanish Projects PID2020-119383GB-I00 (AEI/FEDER, UE) and PROMETEO/2021/030 (Generalitat Valenciana). S.B was supported by the Generalitat Valenciana SEJIGENT program (SEJIGENT/2021/001) and by the European Research Council (ERC) under the ERC2017-STG SENTIFLEX project (Grant Agreement 755617).