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.
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)
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).
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).
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).
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).
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).
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).
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.
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).
Therefore, using formula (16), the observer position of altitude above sea level h,
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.
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.
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).
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).
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).
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)
Hence, the solar noise temperature observed from actual earth station antennas is as shown in formula (25).
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).
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).
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
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).
TU Media Ground Station Information.
TU Media Ground Station Information.
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.
Solar interference time graph from TU media 9.2m antenna calculated through a model (Oct 09, 2007).
Spring quarter calculation result of 2010 (KST 11:20:00 - 11:30:00).
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.
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.
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).
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 -