Satellite orbits in Levi-Civita space

In this paper we consider satellite orbits in central force field with quadratic drag using two formalisms. The first using polar coordinates in which the satellite angular momentum plays a dominant role. The second is in Levi-Civita coordinates in which the energy plays a central role. We then merge these two formalisms by introducing polar coordinates in Levi-Civita space and derive a new equation for satellite orbits which unifies these two paradigms. In this equation energy and angular momentum appear on equal footing and thus characterize the orbit by its two invariants. Using this formalism we show that equatorial orbits around oblate spheroids can be expressed analytically in terms of Elliptic functions. In the second part of the paper we derive in Levi-Civita coordinates a linearized equation for the relative motion of two spacecrafts whose trajectories are in the same plane. We carry out also a numerical verification of these equations.     

Coastal coverage of ESA’ Sentinel 2 mission

Climatologic and anthropogenic pressures in coastal areas affect the coastal zone at different scales. With the development of new missions in open-access, satellites now represent an attractive solution for a broad public to capture local-scale coastal impacts at large scales. Here, the capability of the Sentinel 2 constellation to cover coastal areas and measure coastal processes –physical and biological. We show that Sentinel 2 enables high-frequency measurements across the globe. Cloud coverage at higher latitudes is overcome by decrease revisit time-intervals. Only around the equator, the longest revisit intervals and high cloud cover probability limits coastal measurements there. Sentinel 2 based methods are capable of estimating Digital Elevation Models for mid- to high-latitude coastal zones and sporadic spots for lower latitudes where 2 orbit swaths overlap. For the majority of the world’s coastal bathymetries can be obtained with the Sentinel 2 imagery surpassing the depth of closure (beyond this offshore limit sediment transport is limited). Only in sheltered areas, wave-based bathymetry inversion is limited but at these areas inversion through colouring (light penetration) prevails. This works shows that Sentinel 2 enables coastal monitoring as never before, large spatial scale with revisits of a few days at most of the world.     

在轨卫星异常报警和故障诊断方法研究

及时发现卫星出现的异常是卫星在轨管理的一项重要工作，它关系到卫星在轨使用寿命。能够及时发现并判断卫星异常的方法有很多，本文对这些方法进行了研究，并将其归纳为2类，一类是基于遥测参数超限报警方法，另一类是故障诊断方法。分别对这2类方法从功能、效果、软件实现的难度进行比较，结果表明，在不同的管理和实现难度要求下，这些方法完成卫星管理的效果不同。另外，本文还提出了一种遥测参数相对判断算法，简化了报警门限设置，能够及时发现卫星任何参数突跳，是一种有效的卫星异常报警方法。     

Coastal sea level changes in Europe from GPS,tide gauge,satellite altimetry and GRACE, 1993–2011

Sea level changes are threatening the human living environments, particularly along the European Coasts with highly dense population. In this paper, coastal sea level changes in western and southern Europe are investigated for the period 1993–2011 using Global Positioning System (GPS), Tide Gauge (TG), Satellite Altimetry (SA), Gravity Recovery and Climate Experiment (GRACE) and geophysical models. The mean secular trend is 2.26 ± 0.52 mm/y from satellite altimetry, 2.43 ± 0.61 mm/y from TG+GPS and 1.99 ± 0.67 mm/y from GRACE mass plus steric components, which have a remarkably good agreement. For the seasonal variations, annual amplitudes of satellite altimetry and TG+GPS results are almost similar, while GRACE Mass+Steric results are a little smaller. The annual phases agree remarkably well for three independent techniques. The annual cycle is mainly driven by the steric contributions, while the annual phases of non-steric (mass component) sea level changes are almost a half year later than the steric sea level changes.     

Semi-automatic determination of the Azores Current axis using satellite altimetry: Application to the study of the current variability during 1995–2006

Satellite altimetry has been widely used to study the variability of the ocean currents such as the Azores Current (AzC) in the North Atlantic. Most analyses are performed over the region that encloses the current, thus being somehow affected by other oceanographic signals, e.g., eddies. In this study, a new approach for extracting the axis of a zonal current solely based on satellite altimetry is presented. This is a semi-automatic procedure that searches for the maximum values of the gradient of absolute dynamic topography (ADT), using the geostrophic velocity as auxiliary information. The advantage of this approach is to allow the analyses to be performed over a buffer centered on the current axis instead of using a wider region. It is here applied to the AzC for the period June 1995–October 2006.     

Analysis of the DORIS contributions to ITRF2008

In its function as an ITRS Combination Centre, DGFI is in charge with the computation of an ITRF2008 solution. The computation methodology of DGFI is based on the combination of datum-free normal equations (weekly or session data sets, respectively) of station positions and Earth orientation parameters (EOP) from the geodetic space techniques DORIS, GPS, SLR and VLBI. In this paper we focus on the DORIS part within the ITRF2008 computations. We present results obtained from the analysis of the DORIS time series for station positions, network translation and scale parameters, as well as for the terrestrial pole coordinates. The submissions to ITRF2008 benefit from improved analysis strategies of the seven contributing IDS analysis centres and from a combination of the weekly solutions of station positions and polar motion. The results show an improvement by a factor of two compared to past DORIS data submitted to ITRF2005, which has been evaluated by investigating the repeatabilities of position time series. The DORIS position time series were analysed w.r.t. discontinuities and other non-linear effects such as seasonal variations. About 40 discontinuities have been identified which have been compared with the results of an earlier study. Within the inter-technique combination we focus on the DORIS contribution to the integration of the different space geodetic observations and on a comparison of the geodetic local ties with the space geodetic solutions. Results are given for the 41 co-location sites between DORIS and GPS.     

Unsteady oscillation of a deformable airfoil section in incompressible flow

The unsteady aerodynamic loads (pressure increment and generalized forces) acting on a deformable thin airfoil section are determined for its arbitrary motion in the incompressible flow with transverse gusts. In this case, the exact Küssner solution in series is used for harmonic oscillation of an airfoil, in which the Theodorsen function is approximated by a sum of fractional functions with poles. In the time domain, these functions are replaced by the unknown functions that satisfy the one-type ordinary first order differential equations. The equations obtained are combined with the differential equations of the airfoil motion in the generalized coordinates and are used for analyzing the aeroelastic system stability and calculating its dynamic response to gust loads.     

Comparison of fragments created by low- and hyper-velocity impacts

This paper summarizes two new satellite impact experiments. The objective of the experiments was to investigate the outcome of low- and hyper-velocity impacts on two identical target satellites. The first experiment was performed at a low-velocity of 1.5 km/s using a 40-g aluminum alloy sphere. The second experiment was performed at a hyper-velocity of 4.4 km/s using a 4-g aluminum alloy sphere. The target satellites were 15 cm × 15 cm × 15 cm in size and 800 g in mass. The ratios of impact energy to target mass for the two experiments were approximately the same. The target satellites were completely fragmented in both experiments, although there were some differences in the characteristics of the fragments. The projectile of the low-velocity impact experiment was partially fragmented while the projectile of the hyper-velocity impact experiment was completely fragmented beyond recognition. To date, approximately 1500 fragments from each impact experiment have been collected for detailed analysis. Each piece has been weighed, measured, and analyzed based on the analytic method used in the NASA Standard Breakup Model (2000 revision). These fragments account for about 95% of the target mass for both impact experiments. Preliminary analysis results will be presented in this paper.     

Spin parameters of Low Earth Orbiting satellites Larets and Stella determined from Satellite Laser Ranging data

Satellite Laser Ranging (SLR) measurements contain information about the spin parameters of the fully passive, geodetic satellites. In this paper we spectrally analyze the SLR data of 5 geodetic satellites placed on the Low Earth Orbits: GFZ-1, WESTPAC, Larets, Starlette, Stella, and successfully retrieve the frequency signal from Larets and Stella only. The obtained signals indicate an exponential increase of the spin period of Larets: T = 0.860499·exp(0.0197066·D) [s], and Stella: T = 13.5582·exp(0.00431232·D) [s], where D is in days since launch. The initial spin periods calculated from the first month of the SLR observations are: Larets: T_initial = 0.8239 s, Stella: T_initial = 13.2048 s. Analysis of the apparent effects indicates the counter-clockwise spin direction of the satellites. The twice more heavy Stella lost its rotational energy more than four times slower than Larets. Fitting the spin model to the observed spin trends allows determination of the spin axis orientation evolution for Larets and Stella before their rotational period becomes equal to the orbital period.     

Assessing and minimizing collisions in satellite mega-constellations

We aim to provide satellite operators and researchers with an efficient means for evaluating and mitigating collision risk during the design process of mega-constellations. We first introduce a novel algorithm for conjunction prediction that relies on large-scale numerical simulations and uses a sequence of filters to greatly reduce its computational expense. We then use this brute-force algorithm to establish baselines of endogenous (intra-constellation), or self-induced, conjunction events for the FCC-reported designs of the OneWeb LEO and SpaceX Starlink mega-constellations. We demonstrate how these deterministic results can be used to validate more computationally efficient, stochastic techniques for close-encounter prediction by adopting a new probabilistic approach from Solar-System dynamics as a simple test case. Finally, we show how our methodology can be applied during the design phase of large constellations by investigating Minimum Space Occupancy (MiSO) orbits, a generalization of classical frozen orbits that holistically account for the perturbed-Keplerian dynamics of the Earth-satellite-Moon-Sun system. The results indicate that the adoption of MiSO orbital configurations of the proposed mega-constellations can significantly reduce the risk of endogenous collisions with nearly indistinguishable adjustments to the nominal orbital elements of the constellation satellites.