Spin parameters of nanosatellite BLITS determined from Graz 2 kHz SLR data

The nanosatellite BLITS (Ball Lens In The Space) is the first object designed as a passive, spherical retroreflector of the Luneburg type, dedicated for Satellite Laser Ranging (SLR). The 2 kHz SLR station Graz measures spin parameters of this satellite, providing information about the rotational dynamics of the body. The measurements obtained during the period from September 26, 2009 to November 24, 2010 show a significant change of the spin configuration. The spin axis was dynamically precessing since the launch and currently is sinus-like behaving between coordinates RA 120°…150°, Dec 30°…60° (J2000 inertial reference frame). The angle between the symmetry axis and the spin axis of BLITS is not constant, but is decreasing since the launch, while its spin period is rather stable with a mean value of 5.613 s (clockwise rotation). The satellite was dynamically changing its attitude during the first three months after deployment; after this time the spin parameters are relatively stable.     

Spin axis orientation of Ajisai determined from Graz 2 kHz SLR data

The Graz 2 kHz Satellite Laser Ranging (SLR) measurements allow determination of the spin axis orientation of the geodetic satellite Ajisai. The high repetition rate of the laser makes it possible to determine the epoch time when the laser is pointing directly between two corner cube reflector (CCR) rings of the satellite. Identification of many such events during a few (up to 3) consecutive passes allows to state the satellite orientation in the celestial coordinate system. Six years of 2 kHz SLR data (October 2003–October 2009) delivered 331 orientation values which clearly show precession of the axis along a cone centered at 14^h56^m2.8^s in right ascension and 88.512° in declination (J2000.0 celestial reference frame) and with an half-aperture angle θ of 1.405°. The spin axis precesses with a period of 117 days, which is equal to the period of the right ascension of the ascending node of Ajisai’s orbit. We present a model of the axis precession which allows prediction of the satellite orientation – necessary for the envisaged laser time transfer via Ajisai mirrors.     

Spin rate and spin axis orientation of LARES spectrally determined from Satellite Laser Ranging data

Satellite Laser Ranging (SLR) is a powerful technique able to measure spin rate and spin axis orientation of the fully passive, geodetic satellites. This work presents results of the spin determination of LARES – a new satellite for testing General Relativity. 529 SLR passes measured between February 17 and June 9, 2012, were spectrally analyzed. Our results indicate that the initial spin frequency of LARES is f₀ = 86.906 mHz (RMS = 0.539 mHz). A new method for spin axis determination, developed for this analysis, gives orientation of the axis at RA = 12^h22^m48^s (RMS = 49^m), Dec = −70.4° (RMS = 5.2°) (J2000.0 celestial reference frame), and the clockwise (CW) spin direction. The half-life period of the satellite’s spin is 214.924 days and indicates fast slowing down of the spacecraft.     

10 Years of LAGEOS-1 and 15 years of LAGEOS-2 spin period determination from SLR data

Satellite Laser Ranging (SLR) stations measure distance to the satellites equipped with Corner Cube Reflectors (CCRs). These range measurements contain information about spin parameters of the spacecraft. In this paper we present results of spin period determination of two passive satellites from SLR data only: 10 years of LAGEOS-1 (10426 values), and 15 years of LAGEOS-2 (15580 values). The measurements have been made by standard 10 Hz SLR systems and the first 2 kHz SLR system from Graz (Austria). The obtained data allowed calculation of the initial spin period of the satellites: 0.61 s for LAGEOS-1 and 0.906 s for LAGEOS-2. Long time series of the spin period values show that the satellite’s slowing down rate is not constant but is oscillating with a period of 846 days for LAGEOS-1 and 578 days for LAGEOS-2. The results presented here definitely prove that the SLR is a very efficient technique able to measure spin period of the geodetic satellites.     

New results on spin determination of nanosatellite BLITS from High Repetition Rate SLR data

The nanosatellite BLITS (Ball Lens In The Space) demonstrates a successful design of the new spherical lens type satellite for Satellite Laser Ranging (SLR). The spin parameters of the satellite were calculated from more than 1000 days of SLR data collected from 6 High Repetition Rate (HRR) systems: Beijing, Changchun, Graz, Herstmonceux, Potsdam, Shanghai.     

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.     

Spectral filter for signal identification in the kHz SLR measurements of the fast spinning satellite Ajisai

The high repetition rate satellite laser ranging (SLR) measurements to the fast spinning satellites contain a frequency signal caused by the rotational motion of the corner cube reflector (CCR) array. The spectral filter, developed here, is based on the Lomb algorithm, and is tested with the simulated and the observed high repetition rate SLR data of the geodetic satellite Ajisai (spin period ∼2 s). The filter allows for the noise elimination from the SLR data, and for identification of the returns from the single CCRs of the array – even for the low return rates. Applying the spectral filter to the simulated SLR data increases the S/N ratio by a factor 40–45% for all return rates. Filtering out the noise from the observed data strengthens the frequency signal by factor of ∼25 for the low return rates, which significantly helps to determine the spin phase of the satellite. The spectral filter is applied to the Graz SLR data and the spin rates of Ajisai are determined by two different methods: the frequency analysis and the phase determination of the spinning retroreflector array.The analysis of more than 8 years of the Graz SLR measurements indicates an exponential spin rate trend: f = 0.67034 exp(−0.0148542 Y) [Hz], RMS = 0.085 mHz, where Y is the year since launch. The highly accurate spin rate information demonstrates periodic changes related to the precession of the orbital plane of Ajisai, as it determines the amount of energy received by the satellite from the Sun. The rate of deceleration of Ajisai is not constant: the half life period of the satellite’s spin oscillates around 46.7 years with an amplitude of about 5 years.     

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.     

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):     

Attitude analysis of space debris using SLR and light curve data measured with single-photon detector

Attitude is the important parameter for active debris removal and collision avoidance. This paper deduced the spin axis orientation and spin period of the rocket body, CZ-3B R/B (NORAD ID 38253), using the satellite laser ranging and light curve data measured with single-photon detector at Graz station. The epoch method and LC & SLR residuals fitting were combined to determine these values. The derived right ascension angle was around 220°, the declination angle was near 64° and the sidereal period was calculated to be 117.724 s, for 2017-07-03. The results derived from the two distinct methods were mutually validated. Rocket bodies are a major contributor to space debris and this work provides a reference for attitude determination and attitude modelling.     

Design and performances of laser retro-reflector arrays for Beidou navigation satellites and SLR observations

Beidou is the regional satellite navigation system in China, consisting of three kinds of orbiting satellites, MEO, GEO and IGSO, with the orbital altitudes of 21500–36000 km. For improving the accuracy of satellites orbit determination, calibrating microwave measuring techniques and providing better navigation service, all Beidou satellites are equipped with laser retro-reflector arrays (LRAs) to implement high precision laser ranging. The paper presents the design of LRAs for Beidou navigation satellites and the method of inclined installation of LRAs for GEO satellites to increase the effective reflective areas for the regional ground stations. By using the SLR system, the observations for Beidou satellites demonstrated a precision of centimeters. The performances of these LRAs on Beidou satellites are very excellent.     

Laser measurements to space debris from Graz SLR station

In order to test laser ranging possibilities to space debris objects, the Satellite Laser Ranging (SLR) Station Graz installed a frequency doubled Nd:YAG pulse laser with a 1 kHz repetition rate, a pulse width of 10 ns, and a pulse energy of 25 mJ at 532 nm (on loan from German Aerospace Center Stuttgart – DLR). We developed and built low-noise single-photon detection units to enable laser ranging to targets with inaccurate orbit predictions, and adapted our standard SLR software to include a few hundred space debris targets. With this configuration, we successfully tracked – within 13 early-evening sessions of each about 1.5 h – 85 passes of 43 different space debris targets, in distances between 600 km and up to more than 2500 km, with radar cross sections from >15 m² down to <0.3 m², and measured their distances with an average precision of about 0.7 m RMS.     

ETALON spin period determination from kHz SLR data

Using kHz Satellite LASER Ranging (SLR) data of the SLR station Graz only, we determined the spin periods of the two ETALON satellites – launched into high orbits of about 20,000 km – and their spin period increase during 3 years. The determined spin period values and spin period increase rates at 2004-01-01 are: T_ET1 = 63 s + 0.484 s/year, and T_ET2 = 65.5 s + 0.401 s/year.     

Contribution of satellite laser ranging to combined gravity field models

In the framework of satellite-only gravity field modeling, satellite laser ranging (SLR) data is typically exploited to recover long-wavelength features. This contribution provides a detailed discussion of the SLR component of GOCO02S, the latest release of combined models within the GOCO series. Over a period of five years (January 2006 to December 2010), observations to LAGEOS-1, LAGEOS-2, Ajisai, Stella, and Starlette were analyzed. We conducted a series of closed-loop simulations and found that estimating monthly sets of spherical harmonic coefficients beyond degree five leads to exceedingly ill-posed normal equation systems. Therefore, we adopted degree five as the spectral resolution for real data analysis. We compared our monthly coefficient estimates of degree two with SLR and Gravity Recovery and Climate Experiment (GRACE) time series provided by the Center for Space Research (CSR) at Austin, Texas. Significant deviations in C₂₀ were noted between SLR and GRACE; the agreement is better for the non-zonal coefficients. Fitting sinusoids together with a linear trend to our C₂₀ time series yielded a rate of (−1.75 ± 0.6) × 10⁻¹¹/yr; this drift is equivalent to a geoid change from pole to equator of 0.35 ± 0.12 mm/yr or an apparent Greenland mass loss of 178.5 ± 61.2 km³/yr. The mean of all monthly solutions, averaged over the five-year period, served as input for the satellite-only model GOCO02S. The contribution of SLR to the combined gravity field model is highest for C₂₀, and hence is essential for the determination of the Earth’s oblateness.     

Achievable debris orbit prediction accuracy using laser ranging data from a single station

Earlier studies have shown that an orbit prediction accuracy of 20 arc sec ground station pointing error for 1–2 day predictions was achievable for low Earth orbit (LEO) debris using two passes of debris laser ranging (DLR) data from a single station, separated by about 24 h. The accuracy was determined by comparing the predicted orbits with subsequent tracking data from the same station. This accuracy statement might be over-optimistic for other parts of orbit far away from the station. This paper presents the achievable orbit prediction accuracy using satellite laser ranging (SLR) data of Starlette and Larets under a similar data scenario as that of DLR. The SLR data is corrupted with random errors of 1 m standard deviation so that its accuracy is similar to that of DLR data. The accurate ILRS Consolidated Prediction Format orbits are used as reference to compute the orbit prediction errors. The study demonstrates that accuracy of 20 arc sec for 1–2 day predictions is achievable.     

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.     

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.     

GSFC DORIS contribution to ITRF2008

The NASA GSFC DORIS analysis center has provided weekly DORIS solutions from November 1992 to January 2009 (839 SINEX files) of station positions and Earth Orientation Parameters for inclusion in the DORIS contribution to ITRF2008. The NASA GSFC GEODYN orbit determination software was used to process the orbits and produce the normal equations. The weekly SINEX gscwd10 submissions included DORIS data from Envisat, TOPEX/Poseidon, SPOT-2, SPOT-3, SPOT-4, SPOT-5. The orbits were mostly seven days in length (except for weeks with data gaps or maneuvers). The processing used the GRACE-derived EIGEN-GL04S1 gravity model, updated modeling for time-variable gravity, the GOT4.7 ocean tide model and tuned satellite-specific macromodels for SPOT-2, SPOT-3, SPOT-4, SPOT-5 and TOPEX/Poseidon. The University College London (UCL) radiation pressure model for Envisat improves nonconservative force modeling for this satellite, reducing the median residual empirical daily along-track accelerations from 3.75 × 10⁻⁹ m/s² with the a priori macromodel to 0.99 × 10⁻⁹ m/s² with the UCL model. For the SPOT and Envisat DORIS satellite orbits from 2003 to 2008, we obtain average RMS overlaps of 0.8–0.9 cm in the radial direction, 2.1–3.4 cm cross-track, and 1.7–2.3 cm along-track. The RMS orbit differences between Envisat DORIS-only and SLR & DORIS orbits are 1.1 cm radially, 6.4 cm along-track and 3.7 cm cross-track and are characterized by systematic along-track mean offsets due to the Envisat DORIS system time bias of ±5–10 μs. We obtain a good agreement between the geometrically-determined geocenter parameters and geocenter parameters determined dynamically from analysis of the degree one terms of the geopotential. The intrinsic RMS weekly position repeatability with respect to the IDS-3 combination ranges from 2.5 to 3.0 cm in 1993–1994 to 1.5 cm in 2007–2008.     

Satellite and ground detection of very dense smoke clouds produced on the islands of the Paraná river delta that affected a large region in Central Argentina

Intense fires were produced on the Paraná river delta islands, Argentina, during most part of 2008, by a combination of an exceptionally dry period and the farmers’ use of a fire land-cleaning technique. In April 2008, those fires significantly affected the nearby regions and their inhabitants, from Rosario city to Buenos Aires mega-city. In this work we present satellite as well as ground Aerosol Optical Depth (AOD) at 550 nm data obtained during the propagation of pollution clouds to the central zone of Argentina. The highest value (1.18) was registered at Buenos Aires by atmospheric remote sensing, using the satellite instrument MODIS/Terra on April 18th 2008 at 10:35 local time (= UT − 3 h). On the same day, ground air quality detectors also measured in this city the highest Total Suspended Particle (TSP) value of the month, 2.02 mg/m³. The AOD(550) daily variation at Rosario Astronomical Observatory, which is located near the Paraná riverside, was derived by combining solar ultraviolet erythemal irradiance data (measured with a YES biometre) with model calculations. On April 25th 2008, from 12:00 to 15:30 local time, a rather high and constant AOD(550) value was registered, with a mean value of (0.90 ± 0.21). Cities located on the side of the Rosario–Buenos Aires highway (San Nicolás, Baradero and San Pedro) were also affected, showing a mean AOD(550) between the Rosario and Buenos Aires values. The particulate matter was collected with gridded samplers placed on the Paraná river islands as well as at the Rosario Observatory. They were analysed with a Scanning Electron Microscope (SEM) and mainly showed a biological origin. Even if normally large particles travel small distances from the source, organic aerosol in the range of 40–100 μm and complex asymmetric structures were registered several kilometres away from the aerosol sources on the islands. Another event of intense UV index attenuation (98.6%) occurred on September 18th 2008, due to very dense smoke clouds that extended over the Rosario area for several hours. The clouds were driven away from the fires by East–northeast and East–southeast winds. The minimum value of this index measured around noon allows to derive a maximum AOD(550)_max = (3.65 ± 0.90) at 12:45 local time. Soot clouds extended over the Paraná river, transporting Burned Biomass Debris (BBD) that deposited on Rosario. In particular, burned leaves and small branches with dimensions of 1–20 cm were collected. The mean (BBD) particles deposited on the ground from 7:00 to 19:00 local time were (0.92 ± 0.20) BBD/(m² h).     

Improved kHz-SLR tracking techniques and orbit quality analysis for LEO missions

Satellite gravity field missions such as CHAMP, GRACE and GOCE are designed as low Earth orbiting spacecraft (LEO) with orbit heights of about 250–500 km. The challenging mission objectives require a very precise knowledge of the satellite orbit position in space. For these missions precise orbit information is typically provided by GPS satellite-to-satellite tracking (SST) observations supported by satellite laser ranging (SLR).