Search

The present contribution investigates the possibility to obtain thermospheric neutral density estimates using satellite laser ranging (SLR) observations of low Earth orbiters (LEOs). This approach is based on the analysis of the satellite atmospheric drag, driven by the fact that the drag force is the largest non-gravitational perturbation acting on LEOs. Due to the uncertainty of current thermospheric models, it is the main error source in the LEO orbit determination process. Moreover, the drag is physically related to the thermospheric density distribution, the interaction of the satellite surface with the surrounding thermosphere and thermospheric winds. For this investigation, a spherical satellite called "Atmospheric Neutral Density Experiment-Pollux" (ANDE-P) developed by the Naval Research Laboratory (USA) is adopted as a case study. The satellite flew at the very low altitude of about 350 km. The most important perturbing acceleration at this altitude, the atmospheric drag, is easier to model for a spherical satellite like ANDE-P than for a satellite of complex geometry. A precise orbit determination of ANDE-P was performed with the DGFI Orbit and Geodetic parameter estimation Software (DOGS) over a period of 49 days (16 August 2009 until 3 October 2009) using SLR observations to this satellite and different thermospheric models. In total, we tested four thermospheric models, namely CIRA86, NRLMSISE00, JB2008 and DTM2013. Correspondingly, scale factors of these reference models are estimated with a 6-h resolution. The results confirm that the estimation of force model parameters from SLR measurements of the ANDE-P satellite is sensitive to differences in the density distributions provided by different models. As a consequence, information on the discrepancies between the various models and the true density can be derived from SLR measurements. Moreover, it is found that SLR observations to LEO satellites at very low altitudes are capable to estimate corrections to (scale factors of) the integrated thermospheric density if all other perturbing accelerations are modelled with sufficient accuracy. We derived time series of estimated scale factors of the thermospheric densities provided by the models and obtained the following mean values during the processed period of time at the ANDE-P altitude: 0.65±0.26 for CIRA86, 0.65±0.25 for NRLMSISE00, 0.79±0.24 for DTM2013, and 0.89±0.27 for JB2008. This suggests that all models overestimate the true thermospheric density along the ANDE-P trajectory during the processed period to a certain extent. The thermospheric densities need to be scaled downwards to fit ANDE-P SLR observations with JB2008 requiring the least amount of scaling. [ABSTRACT FROM AUTHOR]

In this paper, we consistently estimate geodetic parameters such as weekly 3-D station coordinates, Earth orientation parameters (EOP) including daily x/y-pole coordinates and the excess length of day ΔLOD, and selected weekly Earth’s gravitational field (Stokes) coefficients up to degree and order 6 from Satellite Laser Ranging measurements to up to 11 geodetic satellites. The SLR constellation consists of LAGEOS-1/2, Etalon-1/2, Stella, Starlette, Ajisai, Larets, LARES, BLITS and WESTPAC, and its observations cover a time span of 38 years ranging from February 16, 1979, to April 30, 2017. If multiple satellites with various altitudes and orbit inclinations are combined, correlations between estimated parameters are significantly reduced. This allows us (i) to investigate the ability of satellite constellations to reduce existing correlations and (ii) to estimate reliable parameters with higher precision compared to the standard 4-satellite constellation (LAGEOS-1/2, Etalon-1/2) which is currently used by the International Laser Ranging Service for the determination of the Terrestrial Reference Frame (TRF) and EOP products. In particular, the Stokes coefficients, EOP and TRF datum parameters (three translations, three rotations, one scale factor), which are highly correlated with satellite-specific orbit parameters, are improved. From our investigations, we found for an 11-satellite solution compared to the above-mentioned 4-satellite solution a decrease in the scatter of the TRF datum parameters of up to 37%, the transformation residuals are decreased by up to 22%, the scatter of the EOP is decreased by up to 22%, and their mean values are decreased by up to 84% w.r.t. the reference solutions. The largest improvement is obtained for the Stokes coefficients which significantly benefit from a combination of multiple satellites (inclinations and orbit altitudes). In total, single coefficients are improved by up to 93% and the overall improvement is up to 74%. Moreover, it could be clearly identified that Ajisai significantly disturbs the TRF solution due to an erroneous center-of-mass correction. We further quantify the impact of specific satellites on the determination of different geodetic parameters and finally evaluate the potential of the existing SLR-tracked spherical satellite constellation to support the goals of GGOS. [ABSTRACT FROM AUTHOR]

The Celestial Reference System (CRS) is currently realized only by Very Long Baseline Interferometry (VLBI) because it is the space geodetic technique that enables observations in that frame. In contrast, the Terrestrial Reference System (TRS) is realized by means of the combination of four space geodetic techniques: Global Navigation Satellite System (GNSS), VLBI, Satellite Laser Ranging (SLR), and Doppler Orbitography and Radiopositioning Integrated by Satellite. The Earth orientation parameters (EOP) are the link between the two types of systems, CRS and TRS. The EOP series of the International Earth Rotation and Reference Systems Service were combined of specifically selected series from various analysis centers. Other EOP series were generated by a simultaneous estimation together with the TRF while the CRF was fixed. Those computation approaches entail inherent inconsistencies between TRF, EOP, and CRF, also because the input data sets are different. A combined normal equation (NEQ) system, which consists of all the parameters, i.e., TRF, EOP, and CRF, would overcome such an inconsistency. In this paper, we simultaneously estimate TRF, EOP, and CRF from an inter-technique combined NEQ using the latest GNSS, VLBI, and SLR data (2005-2015). The results show that the selection of local ties is most critical to the TRF. The combination of pole coordinates is beneficial for the CRF, whereas the combination of ΔUT1 results in clear rotations of the estimated CRF. However, the standard deviations of the EOP and the CRF improve by the inter-technique combination which indicates the benefits of a common estimation of all parameters. It became evident that the common determination of TRF, EOP, and CRF systematically influences future ICRF computations at the level of several μas. Moreover, the CRF is influenced by up to 50μas if the station coordinates and EOP are dominated by the satellite techniques. [ABSTRACT FROM AUTHOR]

The long wavelength part of the Earth's gravity field can be determined, with varying accuracy, from satellite laser ranging (SLR). In this study, we investigate the combination of up to ten geodetic SLR satellites using iterative variance component estimation. SLR observations to different satellites are combined in order to identify the impact of each satellite on the estimated Stokes coefficients. The combination of satellite-specific weekly or monthly arcs allows to reduce parameter correlations of the single-satellite solutions and leads to alternative estimates of the second-degree Stokes coefficients. This alternative time series might be helpful for assessing the uncertainty in the impact of the low-degree Stokes coefficients on geophysical investigations. In order to validate the obtained time series of second-degree Stokes coefficients, a comparison with the SLR RL05 time series of the Center of Space Research (CSR) is done. This investigation shows that all time series are comparable to the CSR time series. The precision of the weekly/monthly $$C_{21}$$ and $$S_{21}$$ coefficients is analyzed by comparing mass-related equatorial excitation functions $$\chi ^{\text {mass}}_{1,2}$$ with geophysical model results and reduced geodetic excitation functions. In case of $$\chi ^{\text {mass}}_{1}$$ , the annual amplitude and phase of the DGFI solution agrees better with three of four geophysical model combinations than other time series. In case of $$\chi ^{\text {mass}}_{2}$$ , all time series agree very well to each other. The impact of $$C_{20}$$ on the ice mass trend estimates for Antarctica are compared based on CSR GRACE RL05 solutions, in which different monthly $$C_{20}$$ time series are used for replacing. We found differences in the long-term Antarctic ice loss of $$12.3$$ Gt/year between the GRACE solutions induced by the different $$C_{20}$$ SLR time series of CSR and DGFI, which is about 13 % of the total ice loss of Antarctica. This result shows that Antarctic ice mass loss quantifications must be carefully interpreted. [ABSTRACT FROM AUTHOR]

The location of the Earth's principal axes of inertia is a foundation for all the theories and solutions of its rotation, and thus has a broad effect on many fields, including astronomy, geodesy, and satellite-based positioning and navigation systems. That location is determined by the second-degree Stokes coefficients of the geopotential. Accurate solutions for those coefficients were limited to the stationary case for many years, but the situation improved with the accomplishment of Gravity Recovery and Climate Experiment (GRACE), and nowadays several solutions for the time-varying geopotential have been derived based on gravity and satellite laser ranging data, with time resolutions reaching one month or one week. Although those solutions are already accurate enough to compute the evolution of the Earth's axes of inertia along more than a decade, such an analysis has never been performed. In this paper, we present the first analysis of this problem, taking advantage of previous analytical derivations to simplify the computations and the estimation of the uncertainty of solutions. The results are rather striking, since the axes of inertia do not move around some mean position fixed to a given terrestrial reference frame in this period, but drift away from their initial location in a slow but clear and not negligible manner. [ABSTRACT FROM AUTHOR]