16 April 2013 Effect of structural recession on the performance of a satellite laser retroreflector
Author Affiliations +
Satellite laser retroreflector (LRR) array is an essential reflecting element for satellite laser ranging. Because of the requirements of manufacturing and assembling, the actual LRR reflecting face is placed underneath the top of assembly. We use the ray tracing method to build up the active retroreflecting region model for LRR with a structural recession effect. Combined with the principle of far-field diffraction intensity and velocity aberration, we give the mathematical model for expressing the structural recession effect on the received optical signal intensity (ROSI) by a ground station. As an important factor for determining the performance of LRR, the distribution of ROSI is simulated for an HY-2A LRR array. The simulation shows that the value of ROSI decreases along with the increment of structural recession size. The bigger optical incident angle should bring about further reduction of ROSI. Only considering the structural recession effect on ROSI, we obtain the regularity of difference of range correction. Before and after considering the structural recession effect, the maximum difference of range correction reaches 9.36 mm for the whole incidence condition. It demonstrates that structural recession cannot be ignored for deciding the performance of LRR.



Satellite laser retroreflector (LRR) is a solid tetrahedron composed of three mutually perpendicular reflecting surfaces and one front face, which is the indispensable component in satellite laser ranging to reflect the laser beams directly back to the ground station.1 For the purpose of preventing it from being damaged during handling, LRR is usually manufactured with mounting tabs to hold it conveniently into the protective housing.23.4 However, due to the protection manners, practical solid retroreflector designs require mounting structures that result in their reflecting surfaces being recessed with respect to their front window faces. With the structural recession effect, the active retroreflecting region (ARR) of LRR is diminished for the oblique incidence of laser beams. Considering the performance of LRR, including the received optical signal intensity (ROSI), is subject to the ARR, the analysis of structural recession effect on ROSI becomes more important.

In the past, the effect of structural recession on ARR has been expounded qualitatively.5 Evidently, investigating the relationship between structural recession size and the LRR performance quantitatively is not enough. In this paper, we apply a ray tracing method to calculate the ARR. The influence of structural recession on the mathematical model of ROSI is then described. Combined with the parameters of an HY-2A satellite LRR, the direct and indirect effects of structural recession on ROSI and range correction of LRR array are studied by computer simulation.


Active Retroreflecting Region Model

The ARR is an aperture on the front face of LRR, which is established by restricting the positions of a ray on the LRR surface.

The schematic of the LRR mounting assembly is illustrated in Fig. 1. The LRR is held by a ring assembly within a cylindrical shield. The upper and lower mounting rings fit over tabs to encircle the LRR. A retainer ring or encloser (RoE) on the top of upper mounting rings presses the LRR assembly against a structural frame. Taking into consideration that the reflecting faces of LRR lie on the bottom of the upper mounting rings, the LRR can be treated as structural recession. Figure 1 shows that the structural recession parameters include the distance d2 from the top of the tab to the front face of LRR, and the distance d1 from the front face of the LRR to the RoE.

Fig. 1

Schematic of laser retroreflector (LRR) mounting assembly; RoE=retainer ring or encloser.


The reflecting face coordinate system (RCS) (x,y,z) used to describe the geometry of the LRR is shown in Fig. 2. For a better indication of each direction of incident ray, we define the front face coordinate system (FCS) (x,y,z), where the z axis is the symmetry of the LRR, and y axis is the projection of the y-axis upon the front face.

Fig. 2

Definition of coordinate system for the reflecting face coordinate system (RCS) and front face coordinate system (FCS).


The coordinate system transformation relationship between RCS and FCS can be expressed by Ref. 6


where M is the transformation matrix, and L is the vertex length of the LRR in Fig. 1.

If the incident ray with the incident angle φ and the azimuth angle θ is incident on the LRR, the direction cosine R0 of the incident ray in the FCS is written by



Then, the incident ray vector can be rewritten as R=MR0 in the RCS. The incident ray undergoes two refractions at the front face and three reflections at the reflecting face. The position and direction at each face can be traced by a ray tracing approach. There are six different orders of reflections determined by where the ray strikes the LRR.

We consider the case in which the incident ray is refracted by the front face and reflected successively by the yz, xz, and xy planes, then refracted again by the front face. According to the vector form of the refraction, if we ignore the effect of the LRR manufacture error on the direction of ray, we can obtain the direction cosine of the refracted and reflected ray in the RCS as


where n is the refractive index of the LRR material.

Because the active reflecting region is represented in the FCS, we assume that the coordinate of the incidence point is (x1,y1,0). For ease of calculation, we get another coordinate from using Eq. (1) with (x1,y1,z1). Therefore, we can calculate the positions (0,y2,z2), (x3,0,z3), (x4,y4,0, and (x5,y5,z5) on the LRR surface by using the direction cosines, incident point coordinates of the ray, and equations for the LRR surface. To give a better restrictive description of the leaving point (x5,y5,z5), we convert its coordinate from the RCS to FCS, and it can be expressed as (x5,y50).

The positions on the front face should be distributed within the circular front face, except for the incident and outgoing rays that passed through them cannot intersect with the RoE. The positions on the reflecting face must be located within the active reflection region, which is the same distribution form for the three reflecting faces. In terms of a schematic of the LRR, Fig. 3 illustrates the region boundary condition (RBC) for the positions on the front face and xy reflecting face.

Fig. 3

Region boundary condition for the positions on laser retroreflector (LRR) surface: (a) front face; (b) x-y reflecting face. The parameter r is the radius of the front face, and Δr and ΔL are the width and height, respectively, of tabs in Fig. 1.


Based on the RBC for the positions with the above-mentioned reflections order, we can express the first ARR as


where Bf, Bxy, Byz, and Bxz indicate the regions surrounded by the front face, xy plane, yz-plane, and xz-plane RBC, respectively.

We can calculate the other five ARRs in the same manner for the remaining orders of reflections, which can be denoted as Qj (j=2,3,,6).


Effect of Structural Recession on the Performance of a Satellite Laser Retroreflector

In a satellite laser ranging system, LRR is used to retroreflect a laser beam transmitted by a ground station, and support the ground station to implement the satellite laser ranging.7 The performance of satellite LRR can be regarded as the ROSI zby the ground station. A schematic of the system can be described in Fig. 4.

Fig. 4

Working principle schematics of a satellite laser ranging system with a laser retroreflector (LRR).


Due to the relative velocity between the satellite and ground station, the reflected laser beam shifts with the slightest angle aberration, which can be given approximately as8



Here, Re is the Earth radius, g=9.8m/s2 is the gravitational acceleration at the surface, h is the satellite height above sea level, and c is the velocity of light. The plane coordinate of the velocity aberration position is written as (hΔφ,0).

Given that the distance between the LRR and ground station is far greater than the ARR, the ROSI is treated as far-field diffraction intensity (FFDI) at the velocity aberration position. According to Kirchhoff’s diffraction equation, we can solve the ROSI approximately by the following expression:


where the amplitude Ut of the field on the LRR front face is assumed to be homogeneous; λ is the wavelength of the incident laser beam; ρ is the distance between points on the LRR front face and velocity aberration position; W is the phase of the reflected beam from the LRR, which has a relationship with the dihedral angle offset and flatness error.9

The ROSI is a function of the optical incidence condition and LRR structural recession size in the case that other parameters of LRR are known. To expatiate the effect of structural recession, we introduce a new factor η, which represents the decrease in the ROSI due to d1 and d2. On the basis of Eq. (6), η can be described as



The attenuation coefficient η is determined by the optical incidence condition (φ,θ) and structural recession sizes, which falls in the range of 0–1.


Simulation for the Effect of Structural Recession on the ROSI

The HY-2A satellite is a marine remote sensing satellite developed by China, which is equipped with the LRR array as a cooperative target for satellite laser ranging.10 The HY-2A satellite LRR array consists of nine LRRs, which are equally mounted on a regular 48-deg pyramid frame with one nadir-looking LRR in the center. The structure of the HY-2A satellite LRR array is shown in Fig. 5.

Fig. 5

Diagram of HY-2A satellite laser retroreflector (LRR) array.


Each LRR is made from fused quartz glass, whose reflecting faces are silver-coated. Except for structural recession sizes, other parameters of the LRR are indicated in Table 1, which are the initial conditions for calculating the ROSI.

Table 1

Parameters of the HY-2A satellite laser retroreflector.

Radius of the front face (mm)16.5
Vertex length (mm)23.3
Tabs height (mm)5
Tabs width (mm)2
Dihedral angle offset (arcsec)1.6”
Surface flatnessλ/10 (λ=632.8  nm)
Index of refraction1.461

According to the parameters of the LRR and ARR model, the distributions of the ARR with different LRR structural recession sizes and an optical incidence condition (30 deg, 0 deg) are illustrated in Fig. 6. The d/L ratio is applied to generalize the calculation results for all values of L.

Fig. 6

Distributions of active retroreflecting regions (ARRs) with different laser retroreflector (LRR) structural recession sizes for φ=30deg, θ=0deg: (a) d1/L=0 d2/L=0; (b) d1/L=0, d2/L=0.15; (c) d1/L=0.15, d2/L=0; (d) d1/L=0.15, d2/L=0.15.


Figure 6 shows that ARRs are composed of six sections that correspond with six orders of reflections, which are indicated as six colors. The shape of all ARRs approaches an ellipse. With the increment of structural recession sizes, the dimension of ARRs diminishes gradually. The reduction effect of the structural recession parameter d1 on the ARRs is evidently much stronger than structural recession parameter d2.

Considering that the ARRs vary with the appearance of structural recession, the attenuation coefficient η associated with the ROSI will present a corresponding change. The simulation in Fig. 7 shows the distribution regularities of η with different structural recession sizes for three optical incidence conditions (10 deg, 0 deg), (20 deg, 0 deg), and (30 deg, 0 deg).

Fig. 7

Distribution of η with different structural recession sizes for three optical incidence conditions (10 deg, 0 deg), (20 deg, 0 deg), and (30 deg, 0 deg): (a) d1/L; (b) d2/L.


We observe that η falls off rapidly due to the augment of structural recession size for the above three incidence conditions. For any given structural recession size, the value of η falls more quickly with an increasing incident angle. By comparing two groups of η distribution curves in Fig. 7(a) and 7(b), we have obtained that structural recession size d1 has a stronger effect on the fall-off extent of η when compared to structural recession size d2 at the same incident angle.

The comprehensive effect of structural recession parameters d1 and d2 on the attenuation coefficient η with an incidence condition (20 deg, 0 deg) is simulated in Fig. 8. The distribution of η can be described approximately as slanted plane. The slope in the d1 direction is obviously bigger than the d2 direction, which validates again, that the influence degree of d1 on η is greater than d2. For the purpose of better assembling the HY-2A LRR, the structural recession sizes of d1 and d2 are chosen as 2.5 and 2.8 mm, respectively. We get the relationship curve between incident angle and attenuation coefficient η in Fig. 9. The attenuation coefficient η decreases rapidly with an increasing incident angle. The value of η is almost equal to zero as a result of the restriction of structural recession effect when incident angle is 40 deg. The range correction is the distance of the centroid of ROSI from the center of the satellite LRR array, which can be given by Ref. 11


where Φ and Θ are incident angle and azimuth angle, respectively, in the LRR array coordinate. The relationship between them and the incidence condition (φm,θm) of individual LRRs can be converted by the structural matrix of LRR array; num is the number of LRRs in LRR array; Sm is the distance of the apparent reflection point for the m’th LRR from the plane through the center of satellite LRR array perpendicular to the incident ray. The definition of Sm is indicated in Fig. 10.

Fig. 8

Distribution of η due to the comprehensive effect of structural recession parameter d1 and d2 for the incidence condition (20 deg, 0 deg).


Fig. 9

Relationship curve between the incident angle and η for the HY-2A laser retroreflector (LRR).


Fig. 10

Diagram of definition for the location of reflection point for m’th laser retroreflector (LRR).


The symbol ΔSm represents the distance between the incidence point and reflection point for m’th LRR.12 According to the definition of the distance Sm, it can be calculated by using a general model:13


where bm is the distance of the LRR front faces from the center of LRR array.

Considering the variety of the ROSI by the structural recession effect, the range correction of the LRR array will present a few undulations. In the general case, if we ignore the structural recession effect on Sm, we may use the following expression to indicate the differences of range correction as



In terms of Eq. (10) and the HY-2A LRR array parameters, Fig. 11 shows the contour distribution of the difference of range correction. In Fig. 11, the numbers around the circumference and radial direction represent the values of Θ and Φ.

Fig. 11

Contour distribution of the difference of range corrections for incidence conditions.


The contour distribution of the differences of range correction is irregular with the change of Θ and Φ. The differences of range correction are confined between -1.81 and 9.36 mm. As for satellite laser ranging with millimeter accuracy, the differences cannot be neglected. It demonstrates the importance of analyzing the structural recession effect on the ROSI.



We use the ray tracing method and FFDI theory of an LRR to investigate the structural recession effect on the performance of an LRR array. The simulations on the ARR and ROSI demonstrate that structural recession parameters have strong effects on them, which result in reducing the ARR and ROSI with the increment of structural recession size. The extent of reduction depends on the optical incidence condition. A greater incident angle should bring about further reduction of the ROSI. On the basis of the analysis results, we obtain the difference of range correction when the structural recession effect on the ROSI is considered for the HY-2A satellite LRR array. The results show that the structural recession effect cannot be neglected when calculating the range correction for an LRR array.


1. D. CurrieS. D. AgnellocG. D. Monache, “A lunar laser ranging retroreflector array for the 21st century,” Acta Astronaut. 68(7–8), 667–680 (2011).AASTCF0094-5765 http://dx.doi.org/10.1016/j.actaastro.2010.09.001 Google Scholar

2. S. Dell’Agnelloet al., “Creation of the new industry-standard space test of laser retroreflectors for the GNSS and LAGEOS,” Adv. Space Res. 47(5), 822–842 (2011).ASRSDW0273-1177 http://dx.doi.org/10.1016/j.asr.2010.10.022 Google Scholar

3. A. PaolozziI. CiufoliniC. Vendittozzi, “Engineering and scientific aspects of LARES satellite,” Acta Astronaut. 69(3–4), 127–134 (2011).AASTCF0094-5765 http://dx.doi.org/10.1016/j.actaastro.2011.03.005 Google Scholar

4. S. Dell’Agnelloet al., “SCF-Test of a Galileo-IOV retroreflector and thermal-optical simulation of a novel GNSS retroreflector array on a critical orbit,” in Proc. IEEE 1st AESS European Conf. on Satellite Telecommunications, pp. 1–6, IEEE, Rome (2012). Google Scholar

5. D. A. Arnold, “Method of calculating retroreflector-array transfer functions,” SAO Special Report #382 (1979). Google Scholar

6. S. LiB. TangH. Zhou, “Calculation on diffraction aperture of cube corner retroreflector,” Chin. Opt. Lett. 6(11), 833–836 (2008).COLHBT1671-7694 Google Scholar

7. R. NairK. W. Goossen, “Effect of face separation in corner-cube reflectors,” Opt. Eng. 48(12), 123003 (2009).OPEGAR0091-3286 http://dx.doi.org/10.1117/1.3275461 Google Scholar

8. J. J. Degnan, “Millimeter accuracy satellites for two color ranging,” in Proc. Int. Workshop Las. Ranging Instrument., Vol. 7, pp. 36–51, National Aeronautics and Space Administration, Annapolis, Maryland (1992). Google Scholar

9. H. ZhouS. LiG. Zheng, “The phase mode of reflected beam from laser retroreflector with dihedral angle and flatness error,” Opt. Rev. 18(1), 1–6 (2011).1340-6000 http://dx.doi.org/10.1007/s10043-011-0001-7 Google Scholar

10. B. WuM. LinZ. Zhang, “Global SLR tracking support for HY-2 satellite precise orbit determination,” in Proc. Int. Workshop Las. Ranging Instrument., Vol. 48, pp. 382–384, Verlag des Bundesamtes für Kartographie und Geodasie, Bad Koetzting, Germany (2011). Google Scholar

11. N. T. Mironovet al., “ETALON-1, -2 center of mass correction and array reflectivity,” in Proc. Int. Workshop Las. Ranging Instrument., Vol. 6, pp. 9–32, National Aeronautics and Space Administration, Annapolis, Maryland (1992). Google Scholar

12. M. W. Fitzmauriceet al., “Prelaunch testing of the laser geodynamic satellite (LAGEOS),” NASA Tech. Paper 1062 (1977). Google Scholar

13. T. OtsuboJ. AmagaiH. Kunimori, “The center of mass correction of the geodetic satellite AJISAI for single-photon laser ranging,” IEEE Trans. Geosci. Rem. Sens. 37(4), 2011–2018 (1999).IGRSD20196-2892 http://dx.doi.org/10.1109/36.774712 Google Scholar



Hui Zhou received his BE in mechanical engineering in 2001, and his PhD in laser remote sensing in 2007 from Wuhan University, China. His current research focuses on the signal process for a satellite laser altimeter and design for a satellite LRR array.

Song Li received her BE in 1986 and ME in 1989 in optical engineering from Wuhan Technical University of Surveying and Mapping, China, and her PhD in surveying engineering in 2002 from Wuhan University, China. Her current research focuses on photoelectric detection, laser interferometry, and laser atmospheric detection.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Hui Zhou, Hui Zhou, Song Li, Song Li, } "Effect of structural recession on the performance of a satellite laser retroreflector," Optical Engineering 52(4), 043001 (16 April 2013). https://doi.org/10.1117/1.OE.52.4.043001 . Submission:


A smart broadband satellite system with elastic capacity
Proceedings of SPIE (April 02 2010)
Sequential orbit determination for GPS 35 and GPS 36 using...
Proceedings of SPIE (September 07 1998)
Probabilistic microcell prediction model
Proceedings of SPIE (June 24 2002)

Back to Top