A numerical investigation on the eccentricity growth of GNSS disposal orbits

We present the results of an extensive numerical exploration performed on the eccentricity growth in MEO associated with two possible end-of-life disposal strategies for GNSS satellites. The study calls attention to the existence of values of initial inclination, longitude of ascending node, and argument of perigee that are more advantageous in terms of long-term stability of the orbit. The important role of the initial epoch and a corresponding periodicity are also shown. The present investigation is influential in view of recent analytical and numerical developments on the chaotic nature of the region due to lunisolar perturbations, but also for the upcoming Galileo and BeiDou constellations.     

An absolute calibration site for radar altimeters in the continental domain: Lake Issykkul in Central Asia

Altimetry missions such as Topex/Poseidon, Jason-1, GFO and ENVISAT have been widely used in the continental domain over lakes, rivers and wetland although they were mostly dedicated to oceanic studies. Knowledge of the instrumental biases is a key issue. Numerous sites have been dedicated to calibration purposes, either in the oceanic domain (Harvest offshore platform in California, Corsica, Bass Strait in Australia) or over lakes (Lake Erie in United States). A new site (Lake Issykkul in Kirghizstan) is proposed for calibration in the continental domain. This lake is covered by past (T/P) and current radar altimetry satellites (Jason-1, T/P, GFO, and ENVISAT). Several in situ water levels and local meteorological variables are available at the site. Located in a mountainous area, it offers an opportunity for calibration far away from all other existing sites and very different environment contexts. Two GPS campaigns have been conducted on the lake in 2004 and in 2005. They consisted of cruises with stations installed onboard a boat following the satellite ground tracks, and onshore settings. This enabled estimating a bias for each altimeter and each tracking algorithm available. Biases obtained for Envisat, GFO, T/P and Jason-1 using the default ocean tracker (respectively, 48.1 ± 6.6, 7.5 ± 4.0, 0 ± 4.3 and 7.0 ± 5.5 cm) agree with biases published at the other calibration sites. For Jason-1, there is a significant disagreement with results obtained in the ocean field (7 cm instead of 13 cm) but is coherent with bias obtained on the Lake Erie site. Erroneous estimates of the sea state bias correction from non-oceanic-like waveforms is discussed as a possible explanation. Errors in the ionospheric, wet and dry tropospheric corrections for the continental domain are also highlighted and quantified.     

Monitoring volcanoes using seismic noise correlations

In this article, we summarize some recent results of measurements of temporal changes of active volcanoes using seismic noise cross-correlations. We first present a novel approach to estimate volcano interior temporal seismic velocity changes. The proposed method allows to measure very small velocity changes (≈ 0.1%) with a time resolution as small as one day. The application of that method to Piton de la Fournaise Volcano (La Réunion Island) shows velocity decreases preceding eruptions. Moreover, velocity changes from noise cross-correlations over 10 years allow to detect transient velocity changes that could be due to long-lasting intrusions of magma without eruptive activity or to pressure buildup associated to the replenishing of the magma reservoir. We also present preliminary results of noise cross-correlation waveform perturbation associated with the occurrence of dike injection and volcanic eruption. We show that such an analysis leads us to locate the areas of dike injection and eruptive fissures at Piton de la Fournaise Volcano. These recent results suggest that monitoring volcanoes using seismic noise correlations should improve our ability to forecast eruptions, their intensity and thus potential environmental impact.     

An Analytical Solution of the Lagrange Equations Valid also for Very Low Eccentricities: Influence of a Central Potential

This paper presents an analytic solution of the equations of motion of an artificial satellite, obtained using non singular elements for eccentricity. The satellite is under the influence of the gravity field of a central body, expanded in spherical harmonics up to an arbitrary degree and order. We discuss in details the solution we give for the components of the eccentricity vector. For each element, we have divided the Lagrange equations into two parts: the first part is integrated exactly, and the second part is integrated with a perturbation method. The complete solution is the sum of the so-called “main” solution and of the so-called “complementary” solution. To test the accuracy of our method, we compare it to numerical integration and to the method developed in Kaula (Theory of Satellite Geodesy, Blaisdell publ. Co., New York. 1966), expressed in classical orbital elements. For eccentricities which are not very small, the two analytical methods are almost equivalent. For low eccentricities, our method is much more accurate.     

Late and middle Pleistocene ungulates dietary diversity in Western Europe indicate variations of Neanderthal paleoenvironments through time and space

Mesowear and microwear on enamel from 763 teeth of middle and late Pleistocene ungulates were analysed to infer the potential of dental wear analysis of faunal remains as a paleoenvironmental and paleoclimatic proxy in relation to climatic changes and diversity of vegetation available in the environment. Fossil localities including levels belonging to two glacial and two interglacial stages were selected in Germany, France, and Spain. At a temporal scale, results indicate that the dietary diversity in ungulates is higher during interglacial phases (MIS 5 and 3) than during pleniglacial phases (MIS 8 and 4). Dietary diversity is concluded to be related to climate-driven vegetation changes which during interglacials lead to increased variety of potential food items available to ungulates. At the geographical scale, during interglacials, changes in diet composition are evident along geographical gradients. The corresponding dietary gradients are proposed to be related to climate and vegetation gradients reflecting more arid climates in the Mediterranean area compared to North-Western Europe. Species consistently represented at all localities investigated are Cervus elaphus (Cervidae, Artiodactyla) and Equus ferus (Equidae, Perissodactyla). C. elaphus populations are found to consistently have less abrasive diets than E. ferus populations but dietary traits of both species varied largely, revealing a significant plasticity in the feeding adaptation of both species. Those traits are concluded to be related to differences in vegetation structure at each locality and complement the evidence that ungulates have broader dietary habits than what is usually assumed.     

The mass-to-light ratio of rich star clusters

We point out a strong time evolution of the mass-to-light conversion factor, η, commonly used to estimate masses of unresolved star clusters from observed cluster spectrophotometric measures. We present a series of gas-dynamical models, coupled with the Cambridge stellar evolution tracks, to compute line-of-sight velocity dispersions and half-light radii weighted by the luminosity. We explore a range of initial conditions, varying in turn the cluster mass and/or density, and the stellar population’s initial mass function. We find that η, and hence the estimated cluster mass, may increase by as much as a factor of three over time-scales of 50 million yr. We apply these results to an hypothetic cluster mass distribution function (d.f.), and show that the d.f. shape may be strongly affected at the low-mass end by this effect. Fitting truncated isothermal (Michie–King) models to the projected light profile leads to over-estimates of the concentration parameter, c, of δ c≈0.3 compared to the same functional fit applied to the projected mass density.     

Relative performance of future altimeter systems and tide gauges in constraining a model of North Sea high-frequency barotropic dynamics

We evaluate in this paper the ability of several altimeter systems, considered separately as well as together with tide gauges, to control the time evolution of a barotropic model of the North Sea shelf. This evaluation is performed in the framework of the particular model errors due to uncertainties in bathymetry. An Ensemble Kalman Filter (EnKF) data assimilation approach is adopted, and observing-systems simulation experiments (OSSEs) are carried out using ensemble spread statistics. The skill criterion for the comparison of observing networks is, therefore, not based on the misfit between two simulations, as done in classic twin experiments, but on the reduction of ensemble variance occurring as a consequence of the assimilation. Future altimeter systems, such as the Wide Swath Ocean Altimeter (WSOA) and satellite constellations, are considered in this work. A single WSOA exhibits, for instance, similar performance as two-nadir satellites in terms of sea-level correction, and is better than three satellites in terms of model velocity control. Generally speaking, the temporal resolution of observations is shown to be of major importance for controlling the model error in these experiments. This result is clearly related to the focus adopted in this study on the specific high-frequency response of the ocean to meteorological forcing. Altimeter systems lack adequate temporal sampling for properly correcting the major part of model error in this context, whereas tide gauges, which provide a much finer time resolution, lead to better global statistical performance. When looking into further detail, tide gauges and altimetry are demonstrated to exhibit an interesting complementary character over the whole shelf, as tide gauge networks make it possible to properly control model error in a ∼100-km coastal band, while high-resolution altimeter systems are more efficient farther from the coast.