Mechanism of error propagation from the subdaily Universal Time model into the celestial pole offsets estimated by VLBI

Within the analysis of space geodetic observations, errors of the applied subdaily Earth rotation model can induce systematic effects in different estimated parameters. In this paper, we focus on the impact of the subdaily Universal Time (UT1) model on the celestial pole offsets (CPO) estimated from very long baseline interferometry (VLBI) observations. We provide a mechanism that describes the error propagation from the subdaily UT1 into the daily CPO.In typical 24-h VLBI sessions the observed quasars are well distributed over the sky. But the observations, if looked at from the Earth-fixed frame, are not homogeneously distributed. The amount of observations performed in different terrestrial directions shows an irregularity which can be roughly compared to the case where the observations are collected in only one Earth-fixed direction. This peculiarity leads to artefacts in VLBI solutions, producing a correlation between the subdaily variations in UT1 and the position of the celestial pole. As a result errors in diurnal terms of the subdaily UT1 model are partly compensated by the estimated CPO. We compute for each 24-h VLBI session from 1990 until 2011 the theoretical response of the CPO to an error in the subdaily UT1 by setting up a least-squares adjustment model and using as input the coordinates of the observed quasars and observation epochs. Then real observed response of the estimated CPO derived from the VLBI session solutions is compared to the predicted one. A very good agreement between the CPO values estimated from VLBI and the predicted values was achieved. The presented model of error propagation from the subdaily UT1 into the daily CPO allows to predict and explain the behaviour of CPO estimates of VLBI solutions computed with different subdaily Earth rotation models, what can be helpful for testing the accuracy of different subdaily tidal models.     

A priori zenith wet delays in the analysis of VLBI CONT sessions

Tropospheric delay is one of the major sources of error in VLBI (Very Long Baseline Interferometry) analysis. The principal component of this error can be accurately computed through reliable surface pressure data —hydrostatic delay— yet there is also a small but volatile component —wet delay— which is difficult to be modelled a priori. In VLBI analysis, troposphere delay is typically modelled in the theoretical delays using Zenith Hydrostatic Delays (ZHD) and a dry mapping function. Zenith Wet Delay (ZWD) is not modelled but estimated in the analysis process. This work studies inter alia the impact of including external GNSS estimates to model a priori ZWD in VLBI analysis, as well as other models of a priori ZWD.In a first stage, two different sources of GNSS troposphere products are compared to VLBI troposphere estimates in a period of 5 years. The solution with the best agreement to VLBI results is injected in the VLBI analysis as a priori ZWD value and is compared to other options to model a priori ZWD. The dataset used for this empirical analysis consists of the six CONT campaigns.It has been found that modelling a priori ZWD has no significant impact either on baseline length and coordinates repeatabilities. Nevertheless, modelling a priori ZWD can change the magnitude of the estimated coordinates a few millimeters in the up component with respect to the non-modelling approach. In addition, the influence of a priori ZWD on Earth Orientation Parameters (EOP) and troposphere estimates —Zenith Total Delays (ZTD) and gradients—has also been analysed, resulting in a small but significant impact on both geodetic products.     

Fractal analysis of short and intermediate GRB and the fireball model

In the present work the possibility of the fractal analysis application for GRB temporal profiles was studied. We have analysed the 4B revised BATSE catalog: temporal profiles of GRB with t₉₀ < 3 s (287 short and 100 intermediate) were studied on TTE data, a sample of 278 intermediate GRB with t₉₀  3 s were studied on DISCSC data. An analysis of the background fractal dimension distributions obtained using TTE and DISCSC data (143 and 110 background regions, respectively), indicates that for both datasets background fractal dimensions D_bgr = 1.5 that the fractal dimension distributions obtained by using these data can be processed simultaneously. The change of the fractal index D_bgr for Poisson statistics – dominated sets with different coefficients of error in counting (up to 10) was studied and D_bgr = 1.5. The ranges of fractal dimension (0.80  D  2.25 for short and 0.85  D  2.01 for intermediate GRB) are shifted over range for theoretical fractal curve (1 < D < 2) due to the finite detector time resolution. There are four subgroups in fractal dimension distribution for short GRB (D = 1.05 ± 0.03, D = 1.31 ± 0.05, D = 1.51 ± 0.04, D = 1.90 ± 0.03) and six subgroups for intermediate one (D = 1.05 ± 0.09, D = 1.24 ± 0.08, D = 1.44 ± 0.07, D = 1.51 ± 0.08, D = 1.64 ± 0.07, D = 1.91 ± 0.1). Time profiles with fractal dimension smaller then background can be obtained by using models with many short chaotic processes in sources, for example, fireball model with shock waves. The range of fractal dimensions for the modelled temporal profiles is 1.213  D  1.400, which can correspond to subgroups of short and intermediate GRB with D = 1.31 and D = 1.24; moreover, the fractal dimension of a simulated indented event and GRB990208 are equal within the error limits for some model parameters and it is possible to obtain smooth temporal profiles with D = D_bgr.     

A linear stability analysis of large-Prandtl-number thermocapillary liquid bridges

A linear stability analysis is applied to determine the onset of oscillatory thermocapillary convection in cylindrical liquid bridges of large Prandtl numbers (4  Pr  50). We focus on the relationships between the critical Reynolds number Re_c, the azimuthal wave number m, the aspect ratio Γ and the Prandtl number Pr. A detailed Re_c–Pr stability diagram is given for liquid bridges with various Γ. In the region of Pr > 1, which has been less studied previously and where Re_c has been usually believed to decrease with the increase of Pr, we found Re_c exhibits an early increase for liquid bridges with Γ around one. From the computed surface temperature gradient, it is concluded that the boundary layers developed at both solid ends of liquid bridges strengthen the stability of basic axisymmetric thermocapillary convection at large Prandtl number, and that the stability property of the basic flow is determined by the “effective” part of liquid bridge.     

A VLBI baseline post-adjustment approach for station velocity estimation in Eurasian continent

Baseline lengths and their time-derivatives among 58 geodetic VLBI stations were fitted by using 4439 observing sessions from the International VLBI Service for Geodesy and Astrometry (IVS). First, the velocities of eight stations in Eurasian continent were set as unknown quantities. Then, two standard global solutions from 3523 IVS sessions and 1110 sessions from database code XA, respectively, were applied prior to all-station coordinates and the non-estimated station velocities. Finally, from the relations among the coordinates, velocities, baseline length and its time-derivative, two types of baseline post-adjustment (BPA) were used to estimate the velocities of the eight stations. We discuss the data processing details, including the effect of different prior values for the stations and the optimal solution.     

Remote sensing of the ring current using energetic neutral atoms

Energetic neutral atom (ENA) images of the storm-time ring current obtained from the ISEE-1 spacecraft provide information for a “zero-order” global model of the energetic ion distribution. With the assumption of isotropic pressure and magnetostatic, non-convective pressure balance, the global system of electrical currents driven by the ion pressure can be calculated using Euler potentials for the divergenceless current density. Radial pressure gradients drive azimuthal currents, and azimuthal pressure gradients drive radial currents. The radial currents cause current lines in the inner magnetosphere to close in the ionosphere, forming a “partial” ring current. The intensities and locations of these field-aligned currents driven into and out of the ionosphere resemble those of the observed Region 2 current system, but not all observed properties of the Region 2 system are reproduced by the “zero-order” model.