The Corrected Allan Variance: Stability Analysis of Frequency Measurements with Missing Data

Atomic clocks are essential elements in a variety of applications, such as global navigation satellite systems. Consequently, monitoring their performances is fundamental. The Allan variance is the key statistical tool for the performance characterization of atomic clocks. This paper proves that the Allan variance computed from frequency measurements with missing data is affected by a bias, which can make it dramatically different from the expected behavior in the full data case. Furthermore, it shows how to eliminate (or largely reduce) this bias by correcting the Allan variance. The corrected Allan variance is obtained for some of the most common atomic clock noise components, and it is validated through numerical simulations

Time-frequency analysis of the Galileo satellite clocks: looking for the J2 relativistic effect and other periodic variations

When observed from the ground, the frequency of the atomic clocks flying on the satellites of a Global Navigation Satellite System is referred to as apparent frequency, because it is observed through the on-board signal generation chain, the propagation path, the relativistic effects, the measurement system, and the clock estimation algorithm. As a consequence, the apparent clock frequency is affected by periodic variations of different origins such as, for example, the periodic component of the J2 relativistic effect, due to the oblateness of the earth, and the clock estimation errors induced by the orbital estimation errors. We present a detailed characterization of the periodic variations affecting the apparent frequency of the Galileo clocks, obtained by applying time-frequency analysis and other signal processing techniques on space clock data provided by the European Space Agency. In particular, we analyze one year of data from three Galileo Passive Hydrogen Masers, flying on two different orbital planes. Time-frequency analysis reveals how the spectral components of the apparent frequency change with time. For example, it confirms that the amplitude of the periodic signal due to the orbital estimation errors depends on the angle between the sun and the orbital plane. Moreover, it allows to find a more precise estimate of the amplitude of the J2 effect, in agreement with the prediction of the general theory of relativity, and it shows that such amplitude suddenly decreases when the corresponding relativistic correction is applied to the data, thus validating the analytical formula used for the correction

Time–frequency analysis of the Galileo satellite clocks: looking for the J2 relativistic effect and other periodic variations

When observed from the ground, the frequency of the atomic clocks flying on the satellites of a Global Navigation Satellite System is referred to as apparent frequency, because it is observed through the on-board signal generation chain, the propagation path, the relativistic effects, the measurement system, and the clock estimation algorithm. As a consequence, the apparent clock frequency is affected by periodic variations of different origins such as, for example, the periodic component of the J2 relativistic effect, due to the oblateness of the earth, and the clock estimation errors induced by the orbital estimation errors. We present a detailed characterization of the periodic variations affecting the apparent frequency of the Galileo clocks, obtained by applying time–frequency analysis and other signal processing techniques on space clock data provided by the European Space Agency. In particular, we analyze one year of data from three Galileo Passive Hydrogen Masers, flying on two different orbital planes. Time–frequency analysis reveals how the spectral components of the apparent frequency change with time. For example, it confirms that the amplitude of the periodic signal due to the orbital estimation errors depends on the angle between the sun and the orbital plane. Moreover, it allows to find a more precise estimate of the amplitude of the J2 effect, in agreement with the prediction of the general theory of relativity, and it shows that such amplitude suddenly decreases when the corresponding relativistic correction is applied to the data, thus validating the analytical formula used for the correction

Robustness tests for an optical time scale

Optical clocks have reached such an impressive accuracy and stability that the future redefinition of the second will be probably based on an optical transition. Consequently, building time scales based on optical clocks has become a key problem. Unfortunately, optical clocks are still laboratory prototypes and are not yet capable of long times of autonomous operation. It is hence critical to understand the impact of this limited optical clock availability on the generated time scale. In this work, after describing a simple and effective optical time scale algorithm, based on the steering of a flywheel oscillator towards the optical clock, we investigate in detail the impact of the limited availability of the optical clock on the performances of the steering algorithm and of the generated time scale through numerical simulations. In particular, we simulate a time scale generated by a hydrogen maser (with a flicker floor of 5.5 × 10−16) steered towards an optical clock, by considering six different scenarios for the availability of the latter, spanning from the ideal one, i.e. continuous operation of the optical clock, to the worst one, i.e. non-uniformly distributed frequency measurements with long unavailability periods. The results prove that the steering algorithm is robust and effective despite its very simple implementation, and it is capable of very good performances in all the considered scenarios, provided that the hydrogen maser behaves nominally. Specifically, they show that a time scale with an accuracy of a few hundreds of picoseconds can be easily realized in the ideal scenario, whereas in a more realistic scenario, with one measurement per week only, the time accuracy is nonetheless of a few nanoseconds, competing with the best time scales currently realized worldwide. The performances degradation due to a non-nominal maser behaviour is also discussed

Stress transmission in granular matter

The transmission of forces through a disordered granular system is studied by means of a geometrical-topological approach that reduces the granular packing into a set of layers. This layered structure constitutes the skeleton through which the force chains set up. Given the granular packing, and the region where the force is applied, such a skeleton is uniquely defined. Within this framework, we write an equation for the transmission of the vertical forces that can be solved recursively layer by layer. We find that a special class of analytical solutions for this equation are L\'evi-stable distributions. We discuss the link between criticality and fragility and we show how the disordered packing naturally induces the formation of force-chains and arches. We point out that critical regimes, with power law distributions, are associated with the roughness of the topological layers. Whereas, fragility is associated with local changes in the force network induced by local granular rearrangements or by changes in the applied force. The results are compared with recent experimental observations in particulate matter and with computer simulations.Comment: 14 pages, Latex, 5 EPS figure

The ac stark shift and space-borne rubidium atomic clocks

open7sìDue to its small size, low weight, and low power consumption, the Rb atomic frequency standard (RAFS) is routinely the first choice for atomic timekeeping in space. Consequently, though the device has very good frequency stability (rivaling passive hydrogen masers), there is interest in uncovering the fundamental processes limiting its long-term performance, with the goal of improving the device for future space systems and missions. The ac Stark shift (i. e., light shift) is one of the more likely processes limiting the RAFS' long-term timekeeping ability, yet its manifestation in the RAFS remains poorly understood. In part, this comes from the fact that light-shift induced frequency fluctuations must be quantified in terms of the RAFS' light-shift coefficient and the output variations in the RAFS' rf-discharge lamp, which is a nonlinear inductively-couple plasma (ICP). Here, we analyze the light-shift effect for a family of 10 on-orbit Block-IIR GPS RAFS, examining decade-long records of their on-orbit frequency and rf-discharge lamp fluctuations. We find that the ICP's light intensity variations can take several forms: deterministic aging, jumps, ramps, and non-stationary noise, each of which affects the RAFS' frequency via the light shift. Correlating these light intensity changes with RAFS frequency changes, we estimate the light-shift coefficient, K-LS, for the family of RAFS: K-LS = -(1.9 +/- 0.3) x 10(-12) /%. The 16% family-wide variation in K-LS indicates that while each RAFS may have its own individual K-LS, the variance of K-LS among similarly designed RAFS can be relatively small. Combining K-LS with our estimate of the ICP light intensity's non-stationary noise, we find evidence that random-walk frequency noise in high-quality space-borne RAFS is strongly influenced by the RAFS' rf-discharge lamp via the light shift effect. Published by AIP Publishing.openFormichella, V.; Camparo, J.; Sesia, I.; Signorile, G.; Galleani, L.; Huang, M.; Tavella, P.Formichella, V.; Camparo, J.; Sesia, Ilaria; Signorile, Giovanna; Galleani, L.; Huang, M.; Tavella, Patrizi

Demonstrator of Time Services based on European GNSS signals: the H2020 DEMETRA Project

During 2015-2016, a European Consortium of 15 partners from 8 different countries, developed the DEMETRA (DEMonstrator of EGNSS services based on Time Reference Architecture), a project funded by the European Union in the frame of the Horizon 2020 program. This project aims at developing and experimenting time dissemination services dedicated to specific users like traffic control, energy distribution, finance, telecommunication, and scientific institutions. Nine services have been developed. These services provide time dissemination with accuracy levels from millisecond to the sub-ns, and also additional services like certification, calibration, or integrity. Five of these services are based on the European GNSS. After a development phase (see PTTI 2016 presentation) the full DEMETRA system has been working during six months for demonstration. The paper will report about the experimentation results, showing performances and limits of the proposed time dissemination services, aiming to foster the exploitation of the European GNSS for timing applications

Geological and Hydrogeological Characterization of Springs in a DSGSD Context (Rodoretto Valley – NW Italian Alps)

As continuous groundwater monitoring in the upper sector of Rodoretto Valley (Germanasca Valley, Italian Western Alps) is hampered by logistical problem of data collection during winter and spring months, the only tools currently available to derive hydrogeological information are non-continuous and non-long-term dataset of spring discharge (Q), temperature (T) and electrical conductivity (EC). In order to quantity aquifer groundwater reserve, available Q dataset of a small mountain spring (Spring 1 CB) was investigated by applying the analytical solutions developed by Boussinesq (J Math Pure Appl 10:5–78, 1904) and Maillet (Essais dı’hydraulique souterraine et fluviale, vol 1. Herman et Cie, Paris, 1905); T and EC datasets were also used to provide qualitative information about the nature of the aquifer that supplies the spring. The outcomes of the elaborations highlighted the limits of applicability of these methods in the presence of a non-continuous Q dataset: both Boussinesq (J Math Pure Appl 10:5–78, 1904) and Maillet (Essais dı’hydraulique souterraine et fluviale, vol 1. Herman et Cie, Paris, 1905) estimated that discharge values as a function of recession time were found to be consistently lower than the available discharge ones and the estimated groundwater volumes stored over time above the spring level turned out to be underestimated. Continuous (hourly value) and long-term Q, EC and T values are, therefore, needful to correctly quantify and to make a proper management of groundwater resources in mountain areas

Measuring Multi-Joint Stiffness during Single Movements: Numerical Validation of a Novel Time-Frequency Approach

This study presents and validates a Time-Frequency technique for measuring 2-dimensional multijoint arm stiffness throughout a single planar movement as well as during static posture. It is proposed as an alternative to current regressive methods which require numerous repetitions to obtain average stiffness on a small segment of the hand trajectory. The method is based on the analysis of the reassigned spectrogram of the arm's response to impulsive perturbations and can estimate arm stiffness on a trial-by-trial basis. Analytic and empirical methods are first derived and tested through modal analysis on synthetic data. The technique's accuracy and robustness are assessed by modeling the estimation of stiffness time profiles changing at different rates and affected by different noise levels. Our method obtains results comparable with two well-known regressive techniques. We also test how the technique can identify the viscoelastic component of non-linear and higher than second order systems with a non-parametrical approach. The technique proposed here is very impervious to noise and can be used easily for both postural and movement tasks. Estimations of stiffness profiles are possible with only one perturbation, making our method a useful tool for estimating limb stiffness during motor learning and adaptation tasks, and for understanding the modulation of stiffness in individuals with neurodegenerative diseases