Estimation of the Love and Shida numbers: h2, l2 using SLR data for the low satellites

In this paper we present results for the global elastic parameters: Love number h₂ and Shida number l₂ derived from the analysis of Satellite Laser Ranging (SLR) data. SLR data for the two low satellites STELLA (H = 800 km) and STARLETTE (H = 810 km) observed during 2.5 years from January 3, 2005 until July 1, 2007 with 18 globally distributed ground stations were analyzed. The analysis was done separately for the two satellites. We do a sequential analysis and study the stability and convergence of the estimates as a function of length of the data set used.     

The impact on tsunami detection from space using GNSS-reflectometry when combining GPS with GLONASS and Galileo

The devastating Sumatra tsunami in 2004 demonstrated the need for a tsunami early warning system in the Indian Ocean. Such a system has been installed within the German-Indonesian Tsunami Early Warning System (GITEWS) project. Tsunamis are a global phenomenon and for global observations satellites are predestined. Within the GITEWS project a feasibility study on a future tsunami detection system from space has therefore been carried out. The Global Navigation Satellite System Reflectometry (GNSS-R) is an innovative way of using GNSS signals for remote sensing. It uses ocean reflected GNSS signals for sea surface altimetry. With a dedicated Low Earth Orbit (LEO) constellation of satellites equipped with GNSS-R receivers, densely spaced sea surface height measurements could be established to detect tsunamis. Some general considerations on the geometry between LEO and GNSS are made in this simulation study. It exemplary analyzes the detection performance of a GNSS-R constellation at 900 km altitude and 60° inclination angle when applied to the Sumatra tsunami as it occurred in 2004. GPS is assumed as signal source and the combination with GLONASS and Galileo signals is investigated. It can be demonstrated, that the combination of GPS and Galileo is advantageous for constellations with few satellites while the combination with GLONASS is preferable for constellations with many satellites. If all three GNSS are combined, the best detection performance can be expected for all scenarios considered. In this case an 18 satellite constellation will detect the Sumatra tsunami within 17 min with certainty, while it takes 53 min if only GPS is considered.     

Estimation of the elastic Earth parameters (h2, l2) using SLR data

We present results for the global elastic parameters h₂ and l₂ derived from the analysis of Satellite Laser Ranging (SLR) data. SLR data for the two satellites LAGEOS 1 and LAGEOS 2 observed during 2.5 years from January 3, 2005 until July 1, 2007 with 18 globally distributed ground stations were analysed using different approaches. The analysis was done separately for the two satellites and approaches to estimate the two elastic parameters independently and together were performed. We do a sequential analysis and study the stability of the estimates as a function of length of the data set used. The adjusted final values for h₂ equal to 0.6151 ± 0.0008 and 0.6152 ± 0.0008, and those for l₂ equal to 0.0886 ± 0.0003 and 0.0881 ± 0.0003 for LAGEOS 1 and LAGEOS 2 tracking data are compared to other independently derived estimates. These parameters and their errors achieve stability at about the 24 and 27 month time interval for h₂ and l₂, respectively.     

Modeling and initial assessment of the inter-frequency clock bias for COMPASS GEO satellites

The COMPASS system is a project established by China to develop an independent global satellite navigation system, which has five GEO (Geostationary Orbit) satellites and thirty Non-GEO satellites. An apparent inter-frequency clock bias (IFCB) for COMPASS GEO satellites is investigated using the real data. The bias also is modeled by the different models. Based on the 15 months (DOY 121, 2011–214, 2012) single-day-estimated results, the periodic variation of IFCBs of the COMPASS GEO satellite is studied using a harmonic analysis. The notable periods of 12 h and 8 h are noted. The harmonics-based models with different periods and different orders and quadratic function based model are used to describe the IFCB. The performances show that the 4-order harmonics-based model with the periods of 24, 12, 8 and 6 h is most optimal than others for describing the IFCB of COMPASS GEO satellite. Its amplitudes and phases estimated from a least square fit are used to study the features of the IFCB. The results show that the current amplitudes and phases do not present special features. Although the irregular amplitudes and phases of the model are disadvantageous for the long-term prediction of IFCB, it is obvious that the modeling IFCB can simple its service and a few of coefficients can replace the IFCB series. The performance of the model in short-term prediction IFCB is tested using the ten-day data (DOY 215-224, 2012).     

Estimation of GPS instrumental biases from small scale network

With 4 GPS receivers located in the equatorial anomaly region in southeast China, this paper proposes a grid-based algorithm to determine the GPS satellites and receivers biases, and at the same time to derive the total electron content (TEC) with time resolution of 15 min and spatial resolution of 1° by 3.5° in latitude and longitude. By assuming that the TEC is identical at any point within a given grid block and the biases do not vary within a day, the algorithm arranges unknown biases and TECs with slant path TEC from the 4 receivers’ observations into a set of equations. Then the instrumental biases and the TECs are determined by using the least squares fitting technique. The performance of the method is examined by applying it to the GPS receiver chain observations selected from 16 geomagnetically quiet days in four seasons of 2006. It is found that the fitting agrees with the data very well, with goodness of fit ranging from 0.452 TECU to 1.914 TECU. Having a mean of 0.9 ns, the standard deviations for most of the GPS satellite biases are less than 1.0 ns for the 16 days. The GPS receiver biases are more stable than that of the GPS satellites. The standard deviation in the 4 receiver bias is from 0.370 ns to 0.855 ns, with a mean of 0.5 ns. Moreover, the instrumental biases are highly correlated with those derived from CODE and JPL with independent methods. The typical precision of the derived TEC is 5 TECU by a conservative estimation. These results indicate that the proposed algorithm is valid and qualified for small scale GPS network.     

Determination of station positions and velocities from laser ranging observations to Ajisai,Starlette and Stella satellites

The positions and velocities of the four Satellite Laser Ranging (SLR) stations: Yarragadee (7090), Greenbelt (7105), Graz (7839) and Herstmonceux (7840) from 5-year (2001–2005) SLR data of low orbiting satellites (LEO): Ajisai, Starlette and Stella were determined. The orbits of these satellites were computed from the data provided by 20 SLR stations. All orbital computations were performed by means of NASA Goddard’s GEODYN-II program. The geocentric coordinates were transformed to the topocentric North–South, East–West and Vertical components in reference to ITRF2005. The influence of the number of normal points per orbital arc and the empirical acceleration coefficients on the quality of station coordinates was studied. To get standard deviation of the coordinates determination lower than 1 cm, the number of the normal points per site had to be greater than 50. The computed positions and velocities were compared to those derived from LAGEOS-1/LAGEOS-2 data. Three parameters were used for this comparison: station coordinates stability, differences from ITRF2005 positions and velocities. The stability of coordinates of LEO satellites is significantly worse (17.8 mm) than those of LAGEOS (7.6 mm), the better results are for Ajisai (15.4 mm) than for Starlette/Stella (20.4 mm). The difference in positions between the computed values and ITRF2005 were little bit worse for Starlette/Stella (6.6 mm) than for LAGEOS (4.6 mm), the results for Ajisai were five times worse (29.7 mm) probably due to center of mass correction of this satellite. The station velocities with some exceptions were on the same level (≈1 mm/year) for all satellites. The results presented in this work show that results from Starlette/Stella are better than those from Ajisai for station coordinates determination. We can applied the data from LEO satellites, especially Starlette and Stella for determination of the SLR station coordinates but with two times lower accuracy than when using LAGEOS data.     

Observation of the lunar topography by the laser altimeter LALT on board Japanese lunar explorer SELENE

The SELENE Laser Altimeter (LALT) is designed to map the Moon’s topography and will be launched in summer 2007. LALT incorporates Q-switched Cr doped Nd:YAG laser (1064 nm) with an output energy of 100 mJ and 1 Hz repetition frequency for about one year mission period. The laser pulse travels to the Moon’s surface and reflections from the surface are detected by a silicon avalanche photo-diode. The ranging distance is 50–150 km with about 5 m accuracy. Several corrections for accurate ranging data are investigated. The flight hardware has been qualified and passed all the integration tests. A principal goal of the LALT instrument is to obtain a much more detailed lunar topographic map which is superior in global coverage, measurement accuracy and number of data points to previous observations and models. The overall science objectives of LALT are (1) determination of lunar global figure, (2) internal structure and surface processes, (3) exploration of the lunar pole regions, and (4) reduction of lunar occultation data.     

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.     

Long-term observations of the trapped high-energy proton population (L < 4) by the NOAA Polar Orbiting Environmental Satellites (POES)

The Space Environment Monitor (SEM) onboard the NOAA POES satellites has been measuring the near-Earth charged particle environment since 1978, providing an extensive database that can be used for studying the long-term behavior of this population of trapped particles. POES stands for Polar Orbiting Environmental Satellite. These satellites orbit at ∼840 km altitude and at an inclination of 98°. The SEM-1 instrument was flown on the POES satellites beginning in 1978 with TIROS-N and NOAA-6 in 1979 and continuing to NOAA-14 launched in 1995 with the exception of NOAA-9 and NOAA-11 (NOAA-13 failed shortly after launch). Its replacement, SEM-2, has flown on the POES NOAA-15, -16, and -17 satellites (from 1998). Here we present the results of a study on the statistical variations of the high-energy trapped proton environment. Among the detectors in SEM, the four SEM-2 omni-directional proton detectors for energies >16 MeV, >36 MeV, >70 MeV, and >140 MeV provide the data most relevant to this study.     

Relativistic electrons high doses at International Space Station and Foton M2/M3 satellites

The paper presents observation of relativistic electrons. Data are collected by the Radiation Risk Radiometer-Dosimeters (R3D) B2/B3 modifications during the flights of Foton M2/M3 satellites in 2005 and 2007 as well as by the R3DE instrument at the European Technology Exposure Facility (EuTEF) on the Columbus External Payload Adaptor at the International Space Station (ISS) in the period February 20 – April 28, 2008. On the Foton M2/M3 satellites relativistic electrons are observed more frequently than on the ISS because of higher (62.8°) inclination of the orbit. At both Foton satellites the usual duration of the observations are a few minutes long. On the ISS the duration usually is about 1 min or less. The places of observations of high doses due to relativistic electrons are distributed mainly at latitudes above 50° geographic latitude in both hemispheres on Foton M2/M3 satellites. A very high maximum is found in the southern hemisphere at longitudinal range 0°–60°E. At the ISS the maximums are observed between 45° and 52° geographic latitude in both hemispheres mainly at longitudes equatorward from the magnetic poles. The measured absolute maximums of dose rates generated by relativistic electrons are found to be as follows: 304 μGy h⁻¹ behind 1.75 g cm⁻² shielding at Foton M2, 2314 μGy h⁻¹ behind 0.71 g cm⁻² shielding at Foton M3 and 19,195 μGy h⁻¹ (Flux is 8363 cm⁻² s⁻¹) behind les than 0.4 g cm⁻² shielding at ISS.     

GRAS radio occultation on-board of Metop

The GRAS radio occultation instrument is flying on Metop-A and belongs to the EPS (EUMETSAT Polar System). GRAS observes GPS satellites in occultation. Within this work, validation of GRAS closed-loop bending angle data against co-located ECMWF profiles extracted from model fields and occultations from the COSMIC constellation of radio occultation instruments is shown. Results confirm the high data quality and robustness, where GRAS shows lower bending angle noise against ECMWF than COSMIC and in terms of occultations per day, one GRAS ≈ two COSMIC satellites. This is partly due to the operational setup of EPS. For the investigation we focus on two observation periods where updates in the ECMWF (March 2009) and COSMIC processing (October 2009) have improved the statistics further. Bending angles biases agree to within 0.5% against ECMWF and to within 0.1% against COSMIC after the updates for altitudes between 8 and 40 km. In addition, we also analyze the impact of the Metop orbit processing on the derived GRAS bending angle data, where different GPS and Metop orbit solutions are analyzed. Results show that a batch based orbit processing would improve in particular the bending angle bias behavior at higher altitudes. Requirements for the operational processing of GRAS data are briefly outlined, options to ease the use of other positioning system satellites in the near future are discussed. A simplified analysis on the observation of several of these systems, e.g. GPS and Galileo, from one platform shows that about 16% of occultations are found within 300 km, ±3 h, thus providing similar information. A constellation of 2 GRAS like instruments would have only about 10% close-by.     

Generating precise and homogeneous orbits for Jason-1 and Jason-2

Driven by the GMES (Global Monitoring for Environment and Security) and GGOS (Global Geodetic Observing System) initiatives the user community has a strong demand for high-quality altimetry products. In order to derive such high-quality altimetry products, precise orbits for the altimetry satellites are a necessity. With the launch of the TOPEX/Poseidon mission in 1992 a still on-going time series of high-accuracy altimetry measurements of ocean topography started, continued by the altimetry missions Jason-1 in 2001 and Jason-2/OSTM in 2008. This paper contributes to the on-going orbit reprocessing carried out by several groups and presents the efforts of the Navigation Support Office at ESA/ESOC using its NAPEOS software for the generation of precise and homogeneous orbits referring to the same reference frame for the altimetry satellites Jason-1 and Jason-2. Data of all three tracking instruments on-board the satellites (beside the altimeter), i.e. GPS, DORIS, and SLR measurements, were used in a combined data analysis. About 7 years of Jason-1 data and more than 1 year of Jason-2 data were processed. Our processing strategy is close to the GDR-C standards. However, we estimated slightly different scaling factors for the solar radiation pressure model of 0.96 and 0.98 for Jason-1 and Jason-2, respectively. We used 30 s sampled GPS data and introduced 30 s satellite clocks stemming from ESOC’s reprocessing of the combined GPS/GLONASS IGS solution. We present the orbit determination results, focusing on the benefits of adding GPS data to the solution. The fully combined solution was found to give the best orbit results. We reach a post-fit RMS of the GPS phase observation residuals of 6 mm for Jason-1 and 7 mm for Jason-2. The DORIS post-fit residuals clearly benefit from using GPS data in addition, as the DORIS data editing improves. The DORIS observation RMS for the fully combined solution is with 3.5 mm and 3.4 mm, respectively, 0.3 mm better than for the DORIS-SLR solution. Our orbit solution agrees well with external solutions from other analysis centers, as CNES, LCA, and JPL. The orbit differences between our fully combined orbits and the CNES GDR-C orbits are of about 0.8 cm for Jason-1 and at 0.9 cm for Jason-2 in the radial direction. In the cross-track component we observe a clear improvement when adding GPS data to the POD process. The 3D-RMS of the orbit differences reveals a good orbit consistency at 2.7 cm and 2.9 cm for Jason-1 and Jason-2. Our resulting orbit series for both Jason satellites refer to the ITRF2005 reference frame and are provided in sp3 file format on our ftp server.     

Dayside ionospheric uplift during strong geomagnetic storms as detected by the CHAMP,SAC-C,TOPEX and Jason-1 satellites

Ionosphere response to severe geomagnetic storms that occurred in 2001–2003 was analyzed using data of global ionosphere maps (GIM), altimeter data from the Jason-1 and TOPEX satellites, and data of GPS receivers on-board CHAMP and SAC-C satellites. This allowed us to study in detail ionosphere redistribution due to geomagnetic storms, dayside ionospheric uplift and overall dayside TEC increase. It is shown that after the interplanetary magnetic field turns southward and intensifies, the crests of the equatorial ionization anomaly (EIA) travel poleward and the TEC value within the EIA area increases significantly (up to ∼50%). GPS data from the SAC-C satellite show that during the main phase of geomagnetic storms TEC values above the altitude of 715 km are 2–3 times higher than during undisturbed conditions. These effects of dayside ionospheric uplift occur owing to the “super-fountain effect” and last few hours while the enhanced interplanetary electric field impinged on the magnetopause.     

Ambiguity resolution in precise point positioning with hourly data for global single receiver

Integer ambiguity resolution (IAR) can improve precise point positioning (PPP) performance significantly. IAR for PPP became a highlight topic in global positioning system (GPS) community in recent years. More and more researchers focus on this issue. Progress has been made in the latest years. In this paper, we aim at investigating and demonstrating the performance of a global zero-differenced (ZD) PPP IAR service for GPS users by providing routine ZD uncalibrated fractional offsets (UFOs) for wide-lane and narrow-lane. Data sets from all IGS stations collected on DOY 1, 100, 200 and 300 of 2010 are used to validate and demonstrate this global service. Static experiment results show that an accuracy better than 1 cm in horizontal and 1–2 cm in vertical could be achieved in ambiguity-fixed PPP solution with only hourly data. Compared with PPP float solution, an average improvement reaches 58.2% in east, 28.3% in north and 23.8% in vertical for all tested stations. Results of kinematic experiments show that the RMS of kinematic PPP solutions can be improved from 21.6, 16.6 and 37.7 mm to 12.2, 13.3 and 34.3 mm for the fixed solutions in the east, north and vertical components, respectively. Both static and kinematic experiments show that wide-lane and narrow-lane UFO products of all satellites can be generated and provided in a routine way accompanying satellite orbit and clock products for the PPP user anywhere around the world, to obtain accurate and reliable ambiguity-fixed PPP solutions.     

22 Years of AJISAI spin period determination from standard SLR and kHz SLR data

Satellite Laser Ranging (SLR) is a powerful and efficient technique to measure spin parameters of satellites equipped with corner cube reflectors. We obtained spin period determination of the satellite AJISAI from SLR data only: 17246 pass-by-pass estimates from standard 1–15 Hz SLR data (14/Aug/1986–30/Dec/2008) and 1444 pass-by-pass estimates (9/Oct/2003–30/Dec/2008) from data of the first 2 kHz SLR system from Graz, Austria. A continuous history of the slowing down of AJISAI spin is derived from frequency analysis, and corrected for the apparent effects. The apparent corrections, elaborated here, allowed very accurate determination of AJISAI initial spin period: 1.4855 ± 0.0007 [s]. The paper identifies also non-gravitational effects as a source of the periodical changes in the rate of slowing down of the satellite.     

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.     

Total electron content measurements in ionospheric physics

With the advent of modern global networks of dual-frequency Global Positioning System (GPS), total electron content (TEC) measurements along slant paths connecting GPS receivers and satellites at 22,000 km have become the largest data set available to ionospheric scientists. The TEC can be calculated from the time and phase delay in the GPS signal using the GPS Toolkit, but an unknown bias will remain. In addition, UHF/VHF radio beacons on board low-Earth-orbiting satellites can also be used to measure the electron content. However, the TEC measurements are obtained by integrating TEC differences between slant paths, but also contain biases. It is often necessary to use data assimilative algorithms like the Ionospheric Data Assimilation Three-Dimensional (IDA3D), and to treat both GPS- and LEO-beacon TEC measurements as relative data in order to conduct ionospheric studies.     

The GLONASS-M satellite yaw-attitude model

The proper modeling of the satellites’ yaw-attitude is a prerequisite for high-precision Global Navigation Satellite System (GNSS) positioning and poses a particular challenge during periods when the satellite orbital planes are partially eclipsed. Whereas a lot of effort has been put in to examine the yaw-attitude control of GPS satellites that are in eclipsing orbits, hardly anything is known about the yaw-attitude behavior of eclipsing GLONASS-M satellites. However, systematic variations of the carrier phase observation residuals in the vicinity of the orbit’s noon and midnight points of up to ±27 cm indicate significant attitude-related modeling issues. In order to explore the GLONASS-M attitude laws during eclipse seasons, we have studied the evolution of the horizontal satellite antenna offset estimates during orbit noon and orbit midnight using a technique that we refer to as “reverse kinematic precise point positioning”. In this approach, we keep all relevant global geodetic parameters fixed and estimate the satellite clock and antenna phase center positions epoch-by-epoch using 30-second observation and clock data from a global multi-GNSS ground station network. The estimated horizontal antenna phase center offsets implicitly provide the spacecraft’s yaw-attitude. The insights gained from studying the yaw angle behavior have led to the development of the very first yaw-attitude model for eclipsing GLONASS-M satellites. The derived yaw-attitude model proves to be much better than the nominal yaw-attitude model commonly being used by today’s GLONASS-capable GNSS software packages as it reduces the observation residuals of eclipsing satellites down to the normal level of non-eclipsing satellites and thereby prevents a multitude of measurements from being incorrectly identified as outliers. It facilitates continuous satellite clock estimation during eclipse and improves in particular the results of kinematic precise point positioning of ground-based receivers.     

Spin parameters of High Earth Orbiting satellites Etalon-1 and Etalon-2 determined from kHz Satellite Laser Ranging data

The high repetition rate Satellite Laser Ranging (SLR) system developed in Graz, Austria, measures ranges to the High Earth Orbiting satellites Etalon-1 and Etalon-2 with the millimeter accuracy. The 2 kHz repetition rate of the laser and the relatively high return rates allow to use the SLR data to calculate the spin parameters of the Etalon satellites. The analysis of the 10 years (October 2003–September 2013) of the SLR data gives trends of the spin axes orientation (J2000 Inertial Reference Frame):     

利用COSMIC数据分析全球对流层顶温度和高度的变化特性

利用COSMIC卫星无线电掩星观测的温度数据分析了全球范围的对流层顶参量的变化特性.COSMIC多颗卫星连续观测的温度数据有很好的全球覆盖和很好的垂直分辨率,是分析全球范围对流层顶变化的重要探测数据.经过数据筛选,得到日平均温度的全年分布数据.采用大气温度递减率判据得到对流层顼温度和对流层顶高度的全球分布和变化.结果表明,全球范国内的对流层顶温度和对流层顶高度分布的基本结果与利用无线电探空资料和欧洲气象局再分析数据得到的结果基本一致;对流层项温度和对流层顶高度呈现明显的季节变化,赤道附近和中纬度地区的对流层顶高度变化与高纬地区的对流层顶高度变化出现反相:对流层顶温度和高度的变化呈现明显的南北半球不对称性.