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.     

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.     

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.     

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.     

Optical response of nanosatellite BLITS measured by the Graz 2 kHz SLR system

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 optical response of BLITS has been measured by the Graz 2 kHz SLR station and compared with the response of the classical retroreflector arrays (RRA) of the Low Earth Orbiting satellites such as ERS-2 and Stella. This work demonstrates that the optical response of BLITS is flat and featureless, comparable with the signature of a point-source or a flat target, and suggests that this innovative design will deliver a higher normal point (NP) accuracy (2.55 mm) than any other SLR target currently in orbit. The high reflectivity of the glassy BLITS (about 60% of the return rate from the multi-reflector Stella) is found to be decreasing by about 30% per year, probably due to the solar irradiation. Detailed analysis of the reflective half-shell demonstrates that a high return rate of SLR measurements can be achieved regardless of the incident angle of the laser beam, thus making the spherical lens a perfect successor of the classical RRA panels mounted on active satellites such as CHAMP, GOCE and GRACE.     

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.     

Quality Assessment of the IDS Contribution to ITRF2008

Doppler Orbitography Radiopositionning Integrated by Satellite (DORIS) is one of the four fundamental techniques contributing to the ITRF. The optimal coverage over the globe of the DORIS observing sites and sites co-located with GPS, allow a strong embedding of DORIS within the ITRF network. DORIS contributes to the access to ITRF through precise orbit determination of altimetric satellites with onboard DORIS receivers. The DORIS contribution to the ITRF2008 is enhanced by the fact that the solutions of seven analysis centers were included in the submitted combined time series of weekly station positions and daily polar motion. We evaluate the quality of the DORIS combined solution in terms of its agreement with the other techniques (VLBI, SLR, GPS) contributing to the ITRF2008 combination. We show in particular that the precisions of the current IDS products range between 1.5 to 2.6 mm for station positions (at the epochs of minimum variances); better than 1 mm/yr in velocities and between 170 and 260 micro-arc-seconds for polar motion, a significant improvement by a factor of three to five, compared to past data used in the ITRF2005 combination. This improvement is certainly due to improved analysis strategies employed by the seven IDS analysis centers that contributed to the combined weekly submitted solutions of station positions and polar motion. A spectral analysis of DORIS station height time series indicates that annual and semi-annual signals are dominant. However, TOPEX draconitic period of about 118 days is still detected in about 20% of the station position power spectra. DORIS height annual signals correlate well with GPS annual signal estimated at some co-located stations, which show that DORIS technique is able to detect loading signals.     

Impact of loading effects on determination of the International Terrestrial Reference Frame

The International Terrestrial Reference Frame (ITRF), as a realization of the International Terrestrial Reference System (ITRS), is represented by a set of station positions and linear velocities. They are intended to be used as regularized coordinates to which some corrections should be added to access instantaneous coordinates. The latest ITRS realization is the ITRF2005, which has integrated time series of station positions to form long-term solutions for the four space geodetic techniques. Currently, a purely linear model is used to parameterize station displacements in the estimation process, plus occasional discontinuities in case of earthquakes or equipment changes. However the input data have been derived without applying surface loading models and so surface loading effects are supposed to be embedded in the coordinates as measured quantities. We evaluate the effect of applying a posteriori loading corrections, which include the effect of atmospheric, non-tidal ocean, and continental water loading, to time series of positions estimated from Satellite Laser Ranging (SLR), Very Long Baseline Interferometry (VLBI), and Global Positioning System (GPS) data. We notice that they reduce about 50% or more of the annual signals in the translation and scale parameter time series of the SLR and VLBI techniques, except in SLR Z translation. In general, the estimated secular frame definition is negligibly affected and estimated positions and velocities are not significantly modified for stations that have accumulated a large number of observations. A multi-technique combination of such derived frames allows concluding that, for some cases, loading model corrections might degrade co-located station coordinates almost as much as they benefit them. However, most significant improvement of the estimated secular coordinates is observed for stations with less than 100 estimated positions as demonstrated with a multi-technique combination.     

Impact of ITRS 2014 realizations on altimeter satellite precise orbit determination

This paper evaluates orbit accuracy and systematic error for altimeter satellite precise orbit determination on TOPEX, Jason-1, Jason-2 and Jason-3 by comparing the use of four SLR/DORIS station complements from the International Terrestrial Reference System (ITRS) 2014 realizations with those based on ITRF2008. The new Terrestrial Reference Frame 2014 (TRF2014) station complements include ITRS realizations from the Institut National de l’Information Géographique et Forestière (IGN) ITRF2014, the Jet Propulsion Laboratory (JPL) JTRF2014, the Deutsche Geodätisches Forschungsinstitut (DGFI) DTRF2014, and the DORIS extension to ITRF2014 for Precise Orbit Determination, DPOD2014. The largest source of error stems from ITRF2008 station position extrapolation past the 2009 solution end time. The TRF2014 SLR/DORIS complement impact on the ITRF2008 orbit is only 1–2 mm RMS radial difference between 1992–2009, and increases after 2009, up to 5 mm RMS radial difference in 2016. Residual analysis shows that station position extrapolation error past the solution span becomes evident even after two years, and will contribute to about 3–4 mm radial orbit error after seven years. Crossover data show the DTRF2014 orbits are the most accurate for the TOPEX and Jason-2 test periods, and the JTRF2014 orbits for the Jason-1 period. However for the 2016 Jason-3 test period only the DPOD2014-based orbits show a strong and statistically significant margin of improvement. The positive results with DTRF2014 suggest the new approach to correct station positions or normal equations for non-tidal loading before combination is beneficial. We did not find any compelling POD advantage in using non-linear over linear station velocity models in our SLR & DORIS orbit tests on the Jason satellites. The JTRF2014 proof-of-concept ITRS realization demonstrates the need for improved SLR+DORIS orbit centering when compared to the Ries (2013) CM annual model. Orbit centering error is seen as an annual radial signal of 0.4 mm amplitude with the CM model. The unmodeled CM signals show roughly a 1.8 mm peak-to-peak annual variation in the orbit radial component. We find the TRF network stability pertinent to POD can be defined only by examination of the orbit-specific tracking network time series. Drift stability between the ITRF2008 and the other TRF2014-based orbits is very high, the relative mean radial drift error over water is no larger than 0.04 mm/year over 1993–2015. Analyses also show TRF induced orbit error meets current altimeter rate accuracy goals for global and regional sea level estimation.     

Refining DORIS atmospheric drag estimation in preparation of ITRF2008

In preparation of ITRF2008, all geodetic technique services (VLBI, SLR, GPS and DORIS) are generating new solutions based on combination of individual analysis centers solutions. These data reprocessing are based on a selection of models, parameterization and estimation strategy unique to each analysis center and to each technique. While a good agreement can be found for models between groups, thanks to the existence of the IERS conventions, a great diversity still exist for parameter estimation, allowing possible future improvements in this direction. The goal of this study is to focus on the atmospheric drag estimation used to generate the new DORIS/IGN ignwd08 time series prepared for ITRF2008. We develop here a method to inter-compare different processing strategies. In a first step, by analyzing single-satellite solutions for a few weeks of data but for a large number of possible analysis strategies, we demonstrate that estimating drag coefficient more frequently (typically every 1–2 h instead of previously every 4–8 h) for the lowest DORIS satellites (SPOTs and Envisat) provides better geodetic results for station coordinates and polar motion. This new processing strategy also solved earlier problem found when processing DORIS data during intense geomagnetic events, such as geomagnetic storms. Differences between drag estimation strategies can mostly be found during these few specific periods of extreme geomagnetic activity (few days per year). In such a case, when drag coefficient is only estimated every 6 h or less often for single-satellite solution, a significant degradation in station coordinate accuracy can be observed (120 mm vs. 20 mm) and significant biases arose in polar motion estimation (5 mas vs. 0.3 mas). In a second step, we reprocessed a full year of DORIS data (2003) in a standard multi-satellite mode. We were able to provide statistics on a more reliable data set and to strengthen these conclusions. Our proposed DORIS analysis is easy to implement in all software packages and is now already used by several analysis centers of the International DORIS Service (IDS) when submitting reprocessed solutions for ITRF2008.     

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.     

IDS contribution to ITRF2008

For the first time, the International DORIS Service (IDS) has produced a technique level combination based on the contributions of seven analysis centers (ACs), including the European Space Operations Center (ESOC), Geodetic Observatory Pecny (GOP), Geoscience Australia (GAU), the NASA Goddard Space Flight Center (GSFC), the Institut Géographique National (IGN), the Institute of Astronomy, Russian Academy of Sciences (INASAN, named as INA), and CNES/CLS (named as LCA). The ACs used five different software packages to process the DORIS data from 1992 to 2008, including NAPEOS (ESA), Bernese (GOP), GEODYN (GAU, GSC), GIPSY/OASIS (INA), and GINS (LCA). The data from seven DORIS satellites, TOPEX/Poseidon, SPOT-2, SPOT-3, SPOT-4, SPOT-5, Envisat and Jason-1 were processed and all the analysis centers produced weekly SINEX files in either variance–covariance or normal equation format. The processing by the analysis centers used the latest GRACE-derived gravity models, forward modelling of atmospheric gravity, updates to the radiation pressure modelling to improve the DORIS geocenter solutions, denser parameterization of empirically determined drag coefficients to improve station and EOP solutions, especially near the solar maximum in 2001–2002, updated troposphere mapping functions, and an ITRF2005-derived station set for orbit determination, DPOD2005. The CATREF software was used to process the weekly AC solutions, and produce three iterations of an IDS global weekly combination. Between the development of the initial solution IDS-1, and the final solution, IDS-3, the ACs improved their analysis strategies and submitted updated solutions to eliminate troposphere-derived biases in the solution scale, to reduce drag-related degradations in station positioning, and to refine the estimation strategy to improve the combination geocenter solution. An analysis of the frequency content of the individual AC geocenter and scale solutions was used as the basis to define the scale and geocenter of the IDS-3 combination. The final IDS-3 combination has an internal position consistency (WRMS) that is 15 to 20 mm before 2002 and 8 to 10 mm after 2002, when 4 or 5 satellites contribute to the weekly solutions. The final IDS-3 combination includes solutions for 130 DORIS stations on 67 different sites of which 35 have occupations over 16 years (1993.0–2009.0). The EOPs from the IDS-3 combination were compared with the IERS 05 C04 time series and the RMS agreement was 0.24 mas and 0.35 mas for the X and Y components of polar motion. The comparison to ITRF2005 in station position shows an agreement of 6 to 8 mm RMS in horizontal and 10.3 mm in height. The RMS comparison to ITRF2005 in station velocity is at 1.8 mm/year on the East component, to 1.2 mm/year in North component and 1.6 mm/year in height.     

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

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.     

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.     

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.     

Some comparisons of orbit determination and parameter estimation based on LAGEOS 1 and LAGEOS 2 data

The modelling accuracy of the LAGEOS 1 orbit was continously improved since its launch in 1976. In spite of these experiences the modelling accuracy of LAGEOS 2 is still about 20 per cent worse. Considering e. g. only the influence of different gravity field models it has been shown that the orbital fits for arc lengths of one month is generally about 25 mm for LAGEOS 1 and more than 30 mm for LAGEOS 2. This is mainly due to the fact that LAGEOS 2 has not yet been used for the determination of most gravity field models. The influences of different model parameters on the estimation of station coordinates, Earth rotation parameters, the geogravitational coefficient, the radiation pressure coefficient, and the empirical acceleration has been studied. The differences and peculiarities of both satellites are discussed. Although the analysis of LAGEOS 2 data still does not reach the high level of LAGEOS 1 combination solutions using both satellites allow new insights, higher accuracies, and a higher time resolution for the parameters and phenomena investigated.     

The effect of geocenter motion on Jason-2 orbits and the mean sea level

We compute a series of Jason-2 GPS and SLR/DORIS-based orbits using ITRF2005 and the std0905 standards ( Lemoine et al., 2010). Our GPS and SLR/DORIS orbit data sets span a period of 2 years from cycle 3 (July 2008) to cycle 74 (July 2010). We extract the Jason-2 orbit frame translational parameters per cycle by the means of a Helmert transformation between a set of reference orbits and a set of test orbits. We compare the annual terms of these time-series to the annual terms of two different geocenter motion models where biases and trends have been removed. Subsequently, we include the annual terms of the modeled geocenter motion as a degree-1 loading displacement correction to the GPS and SLR/DORIS tracking network of the POD process. Although the annual geocenter motion correction would reflect a stationary signal in time, under ideal conditions, the whole geocenter motion is a non-stationary process that includes secular trends. Our results suggest that our GSFC Jason-2 GPS-based orbits are closely tied to the center of mass (CM) of the Earth consistent with our current force modeling, whereas GSFC’s SLR/DORIS-based orbits are tied to the origin of ITRF2005, which is the center of figure (CF) for sub-secular scales. We quantify the GPS and SLR/DORIS orbit centering and how this impacts the orbit radial error over the globe, which is assimilated into mean sea level (MSL) error, from the omission of the annual term of the geocenter correction. We find that for the SLR/DORIS std0905 orbits, currently used by the oceanographic community, only the negligence of the annual term of the geocenter motion correction results in a – 4.67 ± 3.40 mm error in the Z-component of the orbit frame which creates 1.06 ± 2.66 mm of systematic error in the MSL estimates, mainly due to the uneven distribution of the oceans between the North and South hemisphere.     

DPOD2005: An extension of ITRF2005 for Precise Orbit Determination

For Precise Orbit Determination of altimetry missions, we have computed a data set of DORIS station coordinates defined for specific time intervals called DPOD2005. This terrestrial reference set is an extension of ITRF2005. However, it includes all new DORIS stations and is more reliable, as we disregard stations with large velocity formal errors as they could contaminate POD computations in the near future. About 1/4 of the station coordinates need to be defined as they do not appear in the original ITRF2005 realization. These results were verified with available DORIS and GPS results, as the integrity of DPOD2005 is almost as critical as its accuracy. Besides station coordinates and velocities, we also provide additional information such as periods for which DORIS data should be disregarded for specific DORIS stations, and epochs of coordinate and velocity discontinuities (related to either geophysical events, equipment problem or human intervention). The DPOD model was tested for orbit determination for TOPEX/Poseidon (T/P), Jason-1 and Jason-2. Test results show DPOD2005 offers improvement over the original ITRF2005, improvement that rapidly and significantly increases after 2005. Improvement is also significant for the early T/P cycles indicating improved station velocities in the DPOD2005 model and a more complete station set. Following 2005 the radial accuracy and centering of the ITRF2005-original orbits rapidly degrades due to station loss.     

Time-varying solar radiation pressure on Ajisai in comparison with LAGEOS satellites

This study aims to investigate solar radiation pressure acting on the spherical geodetic satellites, Ajisai, LAGEOS-1, and LAGEOS-2. The solar radiation pressure coefficients ($C_R$) are derived in two independent ways: (a) through precise orbit determination (POD) using satellite laser ranging (SLR) data, and (b) through modeling using the optical properties of the satellite surface material. The average $C_R$ value of Ajisai (1.039), as calculated from the time series of $C_R$ POD estimates every 15?days, is consistently smaller than those of LAGEOS-1 (1.140) and LAGEOS-2 (1.103). This difference can be explained by the fact that the surface of Ajisai is mostly covered by mirrors. The Ajisai $C_R$ values estimated by POD show remarkable semi-annual variation, which disagrees with the results of a previous study (Sengoku et al., 1995) predicting that the $C_R$ of Ajisai varies almost annually. We attribute this semi-annual variation to two physical reasons: the non-spherical additional cross-sectional area due to the “attached fitting ring” and the low reflectivity of the surface material in the polar regions. Furthermore, the solar radiation pressure acting on Ajisai varies annually in a direction perpendicular to the sun-satellite vector. Finally, the two independent $C_R$ values of Ajisai agree more when we assume a total solar irradiance of 1361?W/m², whereas the value 1367?W/m² has been commonly used in POD.