Advanced
Prediction of Communication Outage Period between Satellite and Earth station Due to Sun Interference
Prediction of Communication Outage Period between Satellite and Earth station Due to Sun Interference
Journal of Astronomy and Space Sciences. 2010. Mar, 27(1): 31-42
Copyright ©2010, The Korean Space Science Society
  • Received : January 01, 2010
  • Published : March 15, 2010
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Yong-Jun Song
Solar Physics Lab., Dept. of Astronomy and Space Science, Kyung Hee University, Yongin 446-701, Korea
Kap-Sung Kim
School of Space Research, Kyung Hee University, Yongin 446-701, Korea
Ho Jin
School of Space Research, Kyung Hee University, Yongin 446-701, Korea
Byoung-Sun Lee
Satellite Control and Navigation Research Team, ETRI, Daejeon 305-700, Korea
stelle9@khu.ac.kr
Abstract
We developed a computer program to predict solar interference period. To calculate Sun's position, we used DE406 ephemerides and Earth ellipsoid model. The Sun's position error is smaller than 10arcsec. For the verification of the calculation, we used TU media ground station on Seongsu-dong, and MBSAT geostationary communication satellite. We analysis errors, due to satellite perturbation and antenna align. The time error due to antenna align has -35 to +16 seconds at 0.1 o , and -27 to +41 seconds at 0.25 o . The time errors derived by satellite perturbation has 30 to 60 seconds.
Keywords
1. Introduction
Geostationary satellites are located at an altitude of approximately 35,786km above the equator, and revolve in the same angular velocity as earth. Geostationary satellites can therefore, communicate with a ground earth station at all times. However, geostationary satellites also experience communication failure time, twice a year, closely one upon the other in spring and autumn quarters. The communication errors occur when ground station-satellite-the Sun are aligned closely, which occurs during spring and fall equinoxes. At such times, thermal noise emitted from the Sun's surface hits the rear side of the satellite and flows directly into the earth station antenna. This is called solar interference. Studies on duration calculation methods and prediction results of a solar interference phenomenon were implemented by many scientists (Vuong & Forsey 1983, Mohamadi & Lyon 1988, Lin & Yang 1989) abroad, and also by Lee et al. (1991) in Korea. To calculate the time of solar interference, information on precise position of the Sun and earth station antenna systems is necessary. Previous researches used the formula of Van Flandern (Van Flandern & Pulkkinen 1979) when calculating the Sun's position, but it has position error of about 1 arcmin. Using the precise ephemeris DE406, which published by NASA/JPL and the earth ellipsoid model, the study calculated the precise positioning of the Sun as causing error within 10 arcsec. For the verification of the calculation, we used TU media ground station located in Seongsu-dong and the MBSAT satellite operated by TU media.
2. Calculation method for coordinates of the Sun
- 2.1. DE406
In order to obtain the coordinates of celestial bodies belong to the solar system, the equation of motion of the entire solar system should be solved, and relative equations of motion of solar system bodies used in DE406 is shown in formula (1). (Standish et al. 1992)
Lager Image
Here, r i , r i and r i represent location, speed and acceleration of the celestial body i in the solarsystem barycentric coordinate. μ j = Gm j , where G is gravitational constant, m j is mass of celestial body j . Moreover, β, γ are PPN (Parameterized Post-Newtonian) parameters and γ = β = 1.
DE406 solves Formula (1) numerically; obtained the computed value by fitting, using Chebyshev polynomial (Newhall 1989), and saved the coefficients in 64 days interval.
The recursion formula of Chebyshev polynomial can be displayed as shown in formula (2).
Lager Image
Here, n refers to the number of Chebyshev coefficients of the respective celestial body saved in DE406; 12 of them are for the Sun, and 9 for earth. Using such coefficients and formula (2), we can obtain f ( t ) - the value of function from specific time t - as shown in formula (3) below (Press et al. 2002).
Lager Image
a n shown in formula (3) is the Chebyshev coefficient saved in DE406, and f(t), calculated from here, and is the position coordinate of the celestial body displayed in BCRS (Barycentric Celestial Reference System). Suppose we put geocentric-heliocentric position vectors to r in BCRS, which is calculated through DE406, to demonstrate the coordinates of the Sun viewed from a geocentric position, the vectors should be conversed to GCRS -Geocentric Celestial Reference System- using the Bias Matrix B and Precession Matrix P, which are similar to formulas (4) and (5) as follows. At this time, the position vector of the Sun viewed from geocentric is placed to r0. The Bias Matrix carries out a role to change BCRS to J2000.0 coordinate system, while the Precession Matrix calibrates celestial bodies of the current time in accordance with precession. The Bias Matrix and Precession Matrix can be obtained through Astronomical Almanac 2008 (USNO & HMNAO 2008).
Lager Image
Lager Image
Lager Image
Lager Image
Lager Image
Lager Image
where, ζ A , z A , μ A demonstrate precession angle (Lieske et al. 1977). J D TT is JD (Julian Date), calculated through TT (Terrestrial Time). The distance between center of the Earth and the Sun, r' calculated finally from GCRS is shown in formula (10).
Lager Image
Since the vector is shown as three-dimensional Cartesian coordinate, semi diameter of the actual Sun as viewed from the earth can be calculated using formula (11) (Meeus 1998).
Lager Image
where, s 0 is semidiameter of the Sun from the distance 1 AU apart, set to 959'' .63, and ∠ is actual distance between the Earth and the Sun, same to 'r''.
- 2.2 Earth Ellipsoid Model
Suppose we are to display the position Re of observers positioned at geographical coordinate ( le ; φ), altitude above the sea level h , and local sidereal time H using a geocentric coordinate system, we can illustrate as follows via an earth ellipsoid model as shown in Figure 1 (Kim 2005).
Lager Image
Lager Image
Lager Image
where, L is the latitude measured in the geocentric assuming earth is perfect round, Å is the Geoid latitude, representing observer latitude from earth ellipsoid. u, v is the coordinate of equatorial in the observer position and the coordinate to pole-ward located at the earth ellipsoid deemed the Geoids latitude Å. Using formula (12)∼(14), we can obtain formula (15) as follows.
Lager Image
Lager Image
The observer position of geodetic latitude Å and altitude above sea level h, which is illustrated in an earth ellipsoid model.
Formula (16) is obtained when cos L and sin L is calculated using the Pythagorean Theorem and formula (15).
Lager Image
Therefore, using formula (16), the observer position of altitude above sea level h,
Lager Image
Lager Image
where, r e is the length of earth radius assuming earth is a complete sphere, and the length is the same as the equatorial radius of actual earth. The position vector Re, which successes observers located at altitude above sea level h and in the center of the earth ellipsoid of which eccentricity is e , can be displayed as formula (19) as follows.
Lager Image
As shown in Figure 2 , we can illustrate formula (20) suppose we calculate the position vector ρ from observation place to the Sun, using the vector r', which is the distance between the Sun and center of the earth, and observer postion vector Re.
Lager Image
Lager Image
Observer-centric horizon coordinate system and the Earth-centric equatorial coordinate system.
As the vector ρ displayed in three-dimensional geocentric Cartesian coordination, coordinate conversion is required for the 3-axis H and 2-axis (π/2 - φ) direction. This can be illustrated in coordinate transformation matrix as shown in formula (14).
Lager Image
If the vector ρ is conversed to spherical coordinate, we can obtain the azimuth angle and the altitude of the Sun from an observation place. In addition, we calculate an off-axis angle (θ 0 ) between the position of a satellite and the current position of the Sun as formula (22), it is feasible via cosine law of spherical triangle (Lee et al. 1991).
Lager Image
where, El sun is elevation angle of the Sun, Az sun is azimuth angle of the Sun, El ant is elevation angle of the antenna pointing, and Az ant is azimuth angle of the antenna pointing.
3. Solar Interference Calculation
Antenna noise temperature is proportionate to the ratio of which apparent solar temperature and the Sun is included in the antenna beam. Therefore, to calculate antenna noise temperature, gain pattern of earth station antenna, showing the antenna beam and apparent solar temperature, should be defined.
As the gained pattern changes depending on the characteristic of each earth station antenna, calculation should be obtained through directive measurement. However, obtaining characteristics by measuring each antenna would be practically impossible, thus, the study was conducted based on the WARC-79 gain pattern (CCIR 1982), employed by ITU. The WARC-79 gain pattern has a parameter in accordance with the radio frequency transmitted to antennas and the diameter of antennas.
T sun - the off-season sun's noise temperature- can be obtained via formula (23) (Shimabukuro & Stacey 1968, Vuong & Forsey 1983).
Lager Image
where, f is frequency of solar noise and a GHz unit. The solar noise radiation will be extinct while passing through the atmosphere of earth, and the extinction efficiency has a different value according to frequency as shown in formula (24) (Johannsen & Titus 1986)
Lager Image
Hence, the solar noise temperature observed from actual earth station antennas is as shown in formula (25).
Lager Image
Moreover, as the increase of noise temperature in accordance with solar is proportionate depending upon how much solar disc comes into the range of the antenna, thus, it can be integral as formula (26) (Mohamadi & Lyon 1988).
Lager Image
where, p is the single polarization attenuation factor, its value is 0.5, and G (θ, φ) is WARC-79 gain pattern.
Formula (26) above can integrate simple solar disc, such as on φ direction, hence, it can be converted to an integration method as shown in formula (27) for θ direction (Lin & Yang 1989).
Lager Image
where, α s is 'semidiameter of the Sun', obtained from formula (7), θ 0 is 'off-axis angle' obtained from formula (15), and φ i is as follows formula
Lager Image
where, α s is 'semidiameter of the Sun', obtained from formula (7), θ 0 is 'off-axis angle' obtained from formula (15), and φ i is as follows formula
To compare the calculated antenna noise temperature with actual observed value, a barometer enabling quantitative comparison is required, thus, ∠(C/N) is defined as a barometer as follows (Vuong & Forsey 1983, Lin & Yang 1989).
Lager Image
TU Media Ground Station Information.
Lager Image
TU Media Ground Station Information.
Lager Image
Actual solar interference time observed from TU media ground station in Seongsu-dong (Provided by TU media dated Oct. 09, 2007) (TU media 2007).
In formula (29), T sys is the system noise temperature of an antenna system when there is no solar interference, whilst T ant is the antenna noise temperature calculates via formula (27). While T sys can calculate and measure an actual antenna system, general antenna systems has between roughly 250 ∼ 300K values (Lee et al. 1991).
4. Comparison and error analysis of solar interference time
A calculation program is created for solar interference based on the aforementioned formulas. The Visual C++ .Net 2008 was used for programming, and the Origin program was utilized to simulate graphics. In addition, we calculated the solar interference time of TU media ground station in Seongsu-dong displayed in Table 1 , to verify the accuracy of the program, and compared with actual observed value.
The actual C/N value and time observed from TU media ground station on Seoungsu-dong on October 9, 2007 is illustrated in Figure 3 . As of an observation result, the time of solar interference commenced at about 11 h 01 m a.m. on KST, of which reached its maximum at 11 h 02 m 30 s a.m., and over at about 11 h 04 m a.m.. ∠( C/N ) value of maximum interference time is 9.8 dB. Figure 4 is a graph calculated using a model. Where, the calculation is implemented assuming, system noise as 250 K, and efficiency of an antenna system as 65% (Mohamadi & Lyon 1988). The time in Figure 4 displays minute and second units only. As a calculation result, the time of sun interference had commenced at 11 h 01 m 02 s a.m., ended at 11 h 03 m 50 s a.m., showing the maximum interference time as 11 h 02 m 26 s a.m.. Maximum ∠( C=N ) of the maximum interference time is 9.2 dB, which is lower than observation value by 0.6 dB. Through this, we were able to confirm the consistency with a model and actual observation value.
Lager Image
Solar interference time graph from TU media 9.2m antenna calculated through a model (Oct 09, 2007).
Lager Image
Spring quarter calculation result of 2010 (KST 11:20:00 - 11:30:00).
Lager Image
Solar interference time change when the angle of an antenna has 0.1o error (KST 11:20:00 - 11:30:00).
Figure 5 is a graph demonstrates the time of solar interference expected in spring quarter of 2010. According to Figure 5 , solar interference starts between 5th and 6th of March, in the case of spring quarter in 2010, and the occurrence time is between 11 h 24 m a.m. - 11 h 28 m a.m.. Solar interference occurs the most in March 6, while maximum interference time is at 11 h 26 m 19 s a.m., and the maximum ∠( C=N ) value is 9.29 dB. Solar interference occurs throughout many days, and as days go by, the maximum occurrence time shows a tendency to become faster gradually.
An error analysis is implemented based on March 6, 2010. The date is forecasted as the day, which most solar interference would occur in the spring quarter of 2010. The maximum solar interference occurred time on this day was at 11 h 26 m 19 s a.m., and based on the MBSAT satellite information and earth station antenna illustrated in Table 1 , we calculated solar interference occurrence time according to each error, and compared the time, which maximum solar interference occurred, and compared with ∠( C=N ) value. The direction of each earth station antenna shown in Figures 6 and 7 shows the change of ∠( C=N ) value and solar interference time, which would occur in case of error of 0.1 o and 0.25 o , respectively. The El of Figures 6 and 7 signifies the altitude of the antenna, while Az represents the azimuth angle of the antenna.
Figure 6 represents when the error of direction of an antenna shows 0.1 o . The maximum time of solar interference occurs 9 seconds faster than 11 h 26 m 19 s a.m. when the altitude of an antenna is lower than the standard altitude 0.1 o . When the azimuth angle is as little as 0.1 o , the maximum solar interference occurs 35 seconds faster. When altitude is increased by 0.1 o , the maximum interference occurs 13 seconds faster, and when, azimuth angle is increased by 0.1 o , the interference occurs 16 seconds later. Through this, we are able to confirm that when the altitude of an antenna changes by 0.1 o , there is an error range with solar interference time at -35 ∼ +16 seconds. Figure 7 shows the difference when 0.25 o error is given to an antenna. The time of maximum solar interference when altitude is as little as 0.25 o , it occurs 22 seconds faster than standard time, and when azimuth angle is as big as 0.25 o , it occurs 41 seconds late. Thus, suppose we calculate the rest, when the error of bearing of an antenna is 0.25 o , we are able to confirm an error range of -27 ∼ +41 seconds.
Lager Image
Solar interference time change when the angle of an antenna has 0:25± error (KST 11:20:00 - 11:30:00).
When we examined Figure 6 and 7 , we can see that the intensity of maximum solar interference changes when antenna error occurs, and that changes in the width of ∠( C/N ) is even bigger in 0.25 o change, compare to 0.1 o angle change of an antenna, and the error of a solar interference occurrence time increase as much as the antenna pointing angle variance. Through this, we confirmed that the more the accurate measurement is carried on the antennas direction, the higher the accuracy of solar interference time predictions.
Figure 8 , illustrates the change in solar interference time when the position of geostationary satellite changes according to perturbation. Suppose we examine Figure 8 , when the position of a satellite moves to east by 0.05 o , that is when it is located at the latus rectum 144.05 o , the maximum interference time is 11 h 26 m 05 s a.m., which is faster by 14 seconds compare to when it is positioned at the latus rectum 144.0 o , the standard position. When a satellite moves to west by 0.05 o and positioned at 143.95 o east, the time of maximum solar interference is 11 h 26 m 32 s a.m., which is 13 seconds delayed than that of standard position. The change of maximum interference period according to the position of a satellite is 0.05 o , and has error range of -14 ∼ +13 seconds at the time of change. However, in a case where there is more than 0.1 o position change, an error occurs with range of one minute or more. In addition, we can also check the change of intensity of solar interference. However, unlike the error occurrence of antennas, we recognized that there is no significant change of in respective cases when a satellite's position changes.
Lager Image
Change of solar interference time according to the change of satellite position (KST 11:20:00-11:30:00).
5. Conclusion
The purpose of this study is to forecast precise prediction of solar interference time. The study calculated the precise position of the Sun, using DE406 and an earth ellipsoid model, and via a solar interference program, we predicted noise temperature and C/N decline of earth station antenna in accordance with solar interference of stationary satellite and earth station. As a result of applying a communication satellite ground station in Seongsu-dong and MBSAT operated by TU media, we confirmed that they are consistent with actual observation, and verified their accuracy through error analysis.
In-depth of consideration is needed for the intensity of solar interference, not only the affect of the Sun, but also on the specific system of each earth station antenna. Specifically, the size of antennas and a wavelength range of radio wave used in antennas affect the most in the time of solar interference, thus it is necessary to have precise information for a better and precise prediction.
An error according to the position change of a geostationary communication satellite will result a time difference within 30 seconds, if position maintenance is implemented within ±0.05 o . Therefore, we expect there will be no significant affect as practical purposes in precise position calculation of satellites. On the other hand, since the direction of antennas can be a significant error factor, precise calculation of the direction of antennas is imperative to the prediction of solar interference time.
In the study, we calculated the time of solar interference between geostationary communication satellite and earth station. In a case of general low orbit satellite, a solar interference time calculation can be conducted using the same method the study employed, if satellite altitude and azimuth angle is obtained from earth station through the orbital elements of satellites.
Moreover, even the precise position of the Moon can be calculated using DE406. It is known that there is noise of about 250K in the Moon, and the interference time according to the Moon can be calculated by taking advantage of this study.
Acknowledgements
This research was supported byWCU(World Class University) program through the National Research Foundation of Korea funded by the Ministry of Education, Science and Technology (R31-10016) and BK21 (KRF-2005-070-C00059).
References
CCIR 1982 Recommendations and Reports of CCIR VolIV-part I International Telecommunication Union Geneva 198 - 200
Johannsen K. G , Titus L 1986 IEEE Transactions on Instrumentation and Measurement 35 344 -
Kim K. S 2005 Fundamentals of Celestial Mechanics KyungHee University Press Seoul 83 - 105
Lee B. S , Lee H. J , Kim J. M , Chung S. J 1991 Proceedings of KITE Fall Conference 1991 14 (2) 9 -
Lieske J. H , Lederle T , Fricke W , Morando B 1977 A&A 58 1 -
Lin K. T , Yang L. J 1989 Comsat Technical Review 19 311 -
Meeus J 1998 Astronomical Algorithms Willmann-Bell inc Richmond 2nd ed 389 -
Mohamadi F , Lyon D. L 1988 International Journal of Satellite Communications 6 65 -
Newhall X. X 1989 CeMec 45 305 -
Press W. H , Teukolsky S. A , Vettering W. T , Flannery B. P 2002 Numerical Recipes in C++: the art of scientific computing 2nd ed Cambridge University Press NewYork 196-204, 2nd ed 152 - 166
Shimabukuro F. I , Stacey J. M 1968 ApJ 152 777 -
Standish E. M , Newhall X. X , Williams J. G , Yeomans D. K , Seidelmann P. K 1992 Explanatory Supplement to the Astronomical Almanac University Science Books Sausalito 281 -
TU media corp 2007 private communication
USNO & HMNAO 2008 Astronomical Almanac for the year 2008 The Staionery Office London B61-B68 B24 - B33
Van Flandern T. C , Pulkkinen K. F 1979 ApJs 41 391 -
Vuong X. T , Forsey R. J 1983 IEEE Transactions on Broadcasting 29 121 -