Chapter 2 Questions and Solutions Question 1.

Explain what the terms centrifugal and centripetal mean with regard to a satellite in orbit around the earth. A satellite is in a circular orbit around the earth. The altitude of the satellites orbit above the surface of the earth is 1,400 km. (i) What are the centripetal and centrifugal accelerations acting on the satellite in its orbit? Give your answer in m/s2 . (ii) What is the velocity of the satellite in this orbit? Give your answer in km/s. (iii) What is the orbital period of the satellite in this orbit? Give your answer in hours, minutes, and seconds. Note: assume the average radius of the earth is 6,378.137 km and Keplers constant has the value 3.986004418 105 km3 /s2 Solution to question 1: In the case of a satellite orbiting the earth, the centrifugal force on the satellite is a force on the satellite that is directly away from the center of gravity of the earth (FOUT in Fig. 2.1) and the centripetal force is one directly towards the center of gravity of the earth (FIN in Fig. 2.1). The centrifugal force on a satellite will therefore try to fling the satellite away from the earth while the centripetal force will try to bring the satellite down towards the earth. (i) From equation (2.1) centripetal acceleration a = /r2 , where is Keplers constant. The value of r = 6,378.137 + 1,400 = 7,778.137 km, thus a = 3.986004418 105 / (7,778.137)2 = 0.0065885 km/s2 = 6.5885007 m/s2 . From equation (2.3), the centrifugal acceleration is given by a = v 2 /r, where v = the velocity of the satellite in a circular orbit.. From equation (2.5) v = (/r)1/2 = (3.986004418 105 / 7,778.137)1/2 = 7.1586494 km/s and so a = 0.0065885007 km/s2 = 6.5885007 m/s2 . NOTE: since the satellite was in stable orbit, the centrifugal acceleration must be equal to the centripetal acceleration, which we have found to be true here (but we needed only to calculate one of them). (ii) We have already found out the velocity of the satellite in orbit in part (i) (using equation (2.5)) to be 7.1586494 km/s (iii) From equation (2.6), the orbital period T = (2r3/2 )/(1/2 ) = (27,778.1373/2 )/( 3.986004418 105 )1/2 = (4,310,158.598)/(631.3481146) = 6,826.912916 s = 1 hour 53 minutes 46.92 seconds

Question 2 A satellite is in a 322 km high circular orbit. Determine: a. The orbital angular velocity in radians per second; b. The orbital period in minutes; and

2-

c. The orbital velocity in meters per second. Note: assume the average radius of the earth is 6,378.137 km and Keplers constant has the value 3.986004418 105 km3 /s2 . Solution to question 2: It is actually easier to answer the three parts of this question backwards, beginning with the orbital velocity, then calculating the period, and hence the orbital angular velocity. First we will find the total radius of the orbit r = 322 + 6,378.137 km = 6700.137 km (c) From eqn. (2.5), the orbital velocity v = (/r)1/2 = (3.986004418 105 / 6700.137)1/2 = 7.713066 km/s = 7,713.066 m/s. (b) From eqn. (2.6), T = (2r3/2 )/(1/2 ) = (26,700.1373/2 )/( 3.986004418 105 )1/2 = (3,445,921.604)/(631.3481146) = 5,458.037372 seconds = 90.9672895 minutes = 90.97 minutes. (a) The orbital period from above is 5,458.037372 seconds. One revolution of the earth covers 360o or 2 radians. Hence 2 radians are covered in 5,458.037372 seconds, giving the orbital angular velocity as 2 /5,458.037372 radians/s = 0.0011512 radians/s. An alternative calculation procedure would calculate the distance traveled in one orbit (2r = 26700.137 = 42,098.20236 km). This distance is equivalent to 2 radians and so 1 km is equivalent to 2 /42,098.20236 radians = 0.0001493 radians. From above, the orbital velocity was 7.713066 km/s = 7.713066 0.0001493 radians/s = 0.0011512 radians/s.

Question 3. The same satellite in question 2 above (322 km circular orbit) carries a 300 MHz transmitter. a. Determine the maximum frequency range over which the received signal would shift due to Doppler effects if received by a stationary observer suitably located in space. Note: the frequency can be shifted both up and down, depending on whether the satellite is moving towards or away from the observer. You need to determine the maximum possible change in frequency due to Doppler (i.e. 2f). b. If an earth station on the surface of the earth at mean sea level, 6,370 km from the center of the earth, can receive the 300 MHz transmissions down to an elevation angle of 0o , calculate the maximum Doppler shift that this station will observe. Note: Include the earths rotation and be sure you consider the maximum possible Doppler shift for a 322 km circular orbit. Solution to question 3 a. The highest Doppler shift would be observed in the plane of the satellite at the orbital height of the satellite: the satellite would be coming directly at the observer or directly

2-

away from the observer. The maximum Doppler shift would therefore be the sum of these two values. The orbital velocity was calculated in question 2 as 7,713.066 m/s. Using equation (2.44a), f / f T = VT / vp , where f is the Doppler frequency, f T is the frequency of the transmitter at rest, VT is the component of the transmitters velocity directed at the observer, and v p is the phase velocity of light. Since the observer is at orbital height, the component of the transmitters velocity towards the observer is the actual velocity of the satellite. Thus f = (7,713.066 300,000,000) / 2.9979 108 = 7,718.468928 Hz. The maximum Doppler shift therefore = 7,718.468928 2 = 15,436.93786 Hz = 15,436.94 Hz. b. It is best to draw a diagram to see what the set up looks like. Below is a view from above the orbit of the satellite (orthogonal to the orbital plane). The important element in this part of the question is the component of the satellites velocity towards the earth station.

O E E S S

Satellite in 322 km circular orbit

B The Earth

O = vector from the satellite to the earth station S = vector from the origin to the satellite E = vector from the origin to the earth station S = angular coordinate of the satellite E = angular coordinate of the earth station A and B are the points in the satellites orbit when the elevation angle at the earth station is zero, and hence (if the plane of the orbit takes the satellite directly over the earth station at zenith) the point where the Doppler shift is highest either positive or negative.

2-

Using the law of cosines, O2 = E2 + S2 2EScos(E S ) For an orbital height of h, S = E + h, giving O2 = E2 + (E + h)2 2E(E + h)cos(E S ) The component of the satellites velocity towards the earth station is d O/dt and we can obtain it by differentiating the above equation (remembering the E, h, and S may be assumed to be constant values), term by term, thus: 2O(dO /dt) = 2E(E + h )sin(E S )(dE/dt dS /dt) , giving (dO/dt) = ( E( E + h)sin( E S )(dE/dt dS /dt))/( E2 + S2 2EScos(E S ))1/2 dE/dt is the earths rotational angular velocity = 2 radians/24hours = 7.2722052 10-5 radians/second dS /dt = n, where n is the orbital angular velocity = 0.0011512 radians/s from question 2 part (a). The sign depends on whether the orbital motion and the earths rotation are in the same direction. To find the maximum Doppler shift, we need to find (E S ) when the elevation angle is 0o . This is drawn below for one of the geometries (the other being the mirror image on the other side). Local horizon at earth station (elevation angle = 0o ) E

In this figure = (E S ) cos = E/S = E/(E + h) = (6378.137)/(6378.137 +322) = 0.9519413 (E S ) = cos-1 (0.9519413) = 17.8352236o

2-

Thus dO (6378.137 )( 6378.137 + 322) sin( 17.835)( 7.272 10 5 m 1.15 10 3 ) = dt (( 6378.137 ) 2 + ( 6378.137 + 322) 2 2(6378.137)(6378.137 + 322) cos(17.835)) 1 / 2 = 14,100.03935, +16,003.6389 2052.096993

= -6.8710394 or +7.7986757 km/s The satellite may be rotating in the same direction as the angular rotation of the earth or against it. Thus, for the same direction rotation, 2 7.7986757 10 3 300 10 6 2 f = = 15,608.27719 = 15.6 kHz 2.9979 10 8 For the opposite direction rotation 2 f = 2 6.8710394 10 3 300 10 6 = 13,752.41469 = 13.75 kHz 2.9979 10 8

Question 4 What are Keplers three laws of planetary motion? Give the mathematical formulation of Keplers third law of planetary motion. What do the terms perigee and apogee mean when used to describe the orbit of a satellite orbiting the earth? A satellite in an elliptical orbit around the earth has an apogee of 39,152 km and a perigee of 500 km. What is the orbital period of this satellite? Give your answer in hours. Note: assume the average radius of the earth is 6,378.137 km and Keplers constant has the value 3.986004418 105 km3 /s2 . Solution to question 4 Keplers three laws of planetary motion are (see page 22) 1. The orbit of any smaller body about a larger body is always an ellipse, with the center of mass of the larger body as one of the two foci. 2. The orbit of the smaller body sweeps out equal areas in time (see Fig. 2.5). 3. The square of the period of revolution of the smaller body about the larger body equals a constant multiplied by the third power of the semimajor axis of the orbital ellipse.

2-

The mathematical formulation of the third law is T2 = (42 a3 )/, where T is the orbital period, a is the semimajor axis of the orbital ellipse, and is Keplers constant. The perigee of a satellite is the closest distance in the orbit to the earth; the apogee of a satellite is the furthest distance in the orbit from the earth. For the last part, draw a diagram to illustrate the geometry.

R = radius of earth = 6378.137

Perigee 500 km

Apogee 39,152 km

The semimajor axis of the ellipse = (39,152 + (2 6378.137) + 500)/2 = 26,204.137 km The orbital period is T2 = (42 a3 )/ = (42 (26,204.137)3 )/ 3.986004418 105 = 1,782,097,845.0 Therefore, T = 42,214.90075 seconds = 11 hours 43 minutes 34.9 seconds

2-

Question 5 An observation satellite is to be placed into a circular equatorial orbit so that it moves in the same direction as the earths rotation. Using a synthetic aperture radar system, the satellite will store data on surface barometric pressure, and other weather related parameters, as it flies overhead. These data will later be played back to a controlling earth station after each trip around the world. The orbit is to be designed so that the satellite is directly above the controlling earth station, which is located on the equator, once every 4 hours. The controlling earth stations antenna is unable to operate below an elevation angle of 10o to the horizontal in any direction. Taking the earths rotational period to be exactly 24 hours, find the following quantities: a. The satellites angular velocity in radians per second. b. The orbital period in hours. c. The orbital radius in kilometers. d. The orbital height in kilometers. e. The satellites linear velocity in meters per second. f. The time interval in minutes for which the controlling earth station can communicate with the satellite on each pass.

Solution to question 5 First, draw a schematic diagram to illustrate the question. The diagram is from the North Pole, above the earth, with the satellites orbit and the equator in the plane of the paper.

re

rs

re = earth radius, 6378.137 km rs = orbital radius, km S = satellite E = earth station Arrows above E and S show the direction of motion of the earth station and satellite, respectively

2-

a. Let s = satellite angular velocity and e = earths rotational angular velocity. This gives e = 2radians 2 = = 7.2722052 10 5 radians / s 24hours 86,400

If the satellite was directly over the earth station at time t = 0, their angular separation at time t will be = (s e) t. The first overhead pass occurs at = 0 and the second occurs at = 2 . We want = 2 to occur at t = 4 hours = 14,400 seconds. Thus 2 = (s e) 14,400, which gives (s e) = 2 /14,400 = 0.0004363. The possible values are (i) s = e +0.0004363 = 7.2722052 10-5 + 0.0004363 = 0.0005090221 radians/s and (ii) e 0.0004363 = 7.2722052 10-5 0.0004363 = 0.0003635779 radians/s. The first value corresponds to an orbital period of T = 2 /0.0005090221 = 12,343.64194 seconds = 3 hours 25 minutes 43.64 seconds. The second value corresponds to an orbital period of T = 2 /0.0003635779 = 17,281.58082 = 4 hours 48 minutes 1.58 seconds. The second value corresponds to the satellite orbiting in the opposite direction to the earths rotation. The question stated that the satellite rotated in the same direction as the earth and so the satellites angular rotation is 0.0005090221 radians/s. A simple way to check this answer is to note that in 4 hours the earth rotates 60o or /3 radians. A satellite going in the direction of the earths rotation must cover (2 + /3) in 4 hours to catch up with the earth station (that is, one full rotation plus the additional distance the earth has rotated). The angular velocity is therefore (2 + /3)/4 hours = (2 + /3)/14,400 = 0.000509 radians/s (assuming answer (i) above). b. The orbital period, T, is given by the number of radians in one orbit (2 ) divided by the angular velocity (0.0005090221 radians/s) = 12,343.64194 seconds = 3.43 hours (=3 hours 25 minutes 43.64). c. From equation (2.25) orbital angular velocity = (1/2 )/(a3/2 ), where is Keplers constant = 3.986004418 105 km3 /s2 and a is the semimajor axis of the orbit. The orbit is circular in this question and so the radius of the orbit is a. We know and the orbital angular velocity, thus a = (3.986004418 105 )1/3 /(0.0005090221)2/3 = 11,543.96203 km = 11,543.96 km d. The orbital height is the radius of the orbit (from the earths center) minus the earths radius = 11,543.96203 6378.137 km = 5,165.825030 = 5,165.83 km e. The linear velocity of the satellite can be found in two ways. From equation (2.5), the orbital velocity = ( /r)1/2 = ((3.986004418 105 )/( 11,543.96203))1/2 = 5.8761306 km/s = 5.876 km/s. Alternatively, orbital velocity = radius of orbit

2-

orbital angular velocity = 11,543.96203 0.0005090221 = 5.8761306 km/s. f. For this part it is best to draw a diagram illustrating the visibility geometry. o S

Local horizontal 10o E 10o

E = earth station S = satellite

C Communications are possible if o . We can find o by trigonometry. Angle AEC = 100o and, by the law of sines, a/sin(100) = re/sin , where a is the radius of the orbit. From this, sin = re sin(100)/a = (6,378.137 0.9848078)/11,543.96203 = 0.5441146 giving = 32.9641821 =32.96o . Given = 32.96o , angle ACE = 180 100 32.96 = 47.04o . The time the satellite goes from overhead to 10o elevation angle is the angle traveled/relative angular velocity of the satellite = 47.04o in radians/(0.000509 radians/s 0.000072722052) = 0.8209299/0.0004363 = 1,881.572083 seconds = 31.359347 minutes. The total time from one extreme to the other of the visible orbit is twice this value = 62.72 minutes.

2-

Question 6 What is the difference, or are the differences, between a geosynchronous satellite and a geostationary satellite orbit? What is the period of a geostationary satellite? What is the name given to this orbital period? What is the velocity of a geostationary satellite in its orbit? Give your answer in km/s. A particular shuttle mission released a TDRSS satellite into a circular low orbit, with an orbital height of 270 km. The shuttle orbit was inclined to the earths equator by approximately 28o . The TDRSS satellite needed to be placed into a geostationary transfer orbit (GTO) once released from the shuttle cargo bay, with the apogee of the GTO at geostationary altitude and the perigee at the height of the shuttles orbit. (i) What was the eccentricity of the GTO? (ii) What was the period of the GTO? (iii) What was the difference in velocity of the satellite in GTO between when it was at apogee and when it was at perigee? Note: assume the average radius of the earth is 6,378.137 km and Keplers constant has the value 3.986004418 105 km3 /s2 . Solution to question 6 A geostationary satellite orbit is one that has zero inclination to the equatorial plane, is perfectly circular (eccentricity is zero), and is at the correct orbital height to remain apparently stationary in orbit as viewed from the surface of the earth. A geosynchronous satellite orbit has most of the attributes of a geostationary orbit, but is either not exactly circular, not in the equatorial plane, or not at exactly the correct orbital height. From Table 2.1, the orbital period of a geostationary satellite is 23 hours, 56 minutes, and 4.1 seconds. The orbital period of a geostationary satellite is called a sidereal day. From Table 2.1, the velocity of a geostationary satellite is 3.0747 km/s. (i) The Geostationary Transfer Orbit (GTO) will have an apogee of 35,786.03 km (the geostationary altitude) and a perigee of 270 km (the release altitude of the TDRSS). The semimajor axis a = (2re + hp + ha)/2 = (2 6378.137 +270 + 35,786.03)/2 = 24,406.152 km From equation (2.27) and example 2.1.3, ro = re + hp and the eccentric anomaly E = 0 when the satellite is at perigee. From equation (2.27) ro = a(1 ecosE), with cosE = 1. Therefore, re + hp = a(1 e) and, rearranging the equation, e = 1 ( re + hp )/a = 1 (6378.137 + 270)/24,406.152 = 0.727604. The eccentricity of the GTO is 0.728. (ii) The orbital period T = ((4 2 a3 )/)1/2 = ((42 24,406.1523 )/ 3.986004418 105 )1/2 = 37,945.47102 seconds = 10 hours 32 minutes 25.47 seconds.

2- 10

Question 7 For a variety of reasons, typical minimum elevation angles used by earth stations operating in the commercial Fixed Services using Satellites (FSS) communications bands are as follows: C-Band 5o ; Ku-Band 10o ; and Ka-Band 20o . (i) Determine the maximum and minimum range in kilometers from an earth station to a geostationary satellite in the three bands. (ii) To what round-trip signal propagation times do these ranges correspond? You may assume the signal propagates with the velocity of light in a vacuum even when in the earths lower atmosphere.

Solution to question 7 Below is a schematic of one of the three frequency situations, in this case Ku-band (10o minimum elevation angle). Smin

Smax

10o

Smax = geostationary satellite position with maximum slant range to earth station Smin = geostationary satellite position with minimum slant range to earth station E = earth station C = center of earth re = earth radius = 6378.137 km = CE hs = height of satellite = ESmin = 35,786.03 km

2- 11

Part (i) By the law of sines, CSmax/sin100o = CE/sine giving (35,786.03 + 6378.137)/sin100o = 6378.137/sin and thus sin = (sin100o 6378.137)/42,164.167 = 0.1489710 and so = 8.5673001 Angle = 180 100 = 71.4326999 Again, by the law of sines, ES max/sin = CSmax /sin100o giving ESmax = (42,164.167 sin71.4326999)/sin100 = 40,586.12894 km For a minimum elevation angle of 5o , = 8.6671185, = 180 95 8.86671185 = 76.3328815o , and ES max = 41.126.78334 km For a minimum elevation angle of 20o , = 8.1720740, = 180 110 8.1720740 = 61.8279266o , and ES max = 39,554.56520 km Maximum ranges are 41,126.78 km (C-band), 40,586.13 km (Ku-band), and 39,554.57 km (Ka-band). The minimum range will be the same for all three frequencies of operation = 35,786.03 km, which assumes the earth station is on the equator and the satellite is directly above the earth station. Part (ii) The round trip propagation times, assuming the velocity of light is 2.997 108 m/s, are given by distance/velocity = 2 41,126.78 km/2.997 108 m/s (C-band), 2 40,586.13 km/2.997 108 m/s (Ku-band), and 2 39,554.57 km/2.997 108 m/s (Ka-band) = 274.5, 270.8, and 264 milliseconds, respectively.

Question 8 Most commercial geostationary communications satellites must maintain their orbital positions to within 0.05o of arc. If a geostationary satellite meets this condition (i.e. it has an apparent motion 0.05o of arc N-S and 0.05o of arc E-W, as measured from the center of the earth), calculate the maximum range variation to this satellite from an earth station with a mean elevation angle to the center of the satellites apparent motion of 5o . You may assume that the equatorial and polar diameters of the earth are the same.

Solution to question 8 No information was given on eccentricity in this question, so we will assume that it is 0.001.

2- 12

We can treat the station keeping variation of the satellite as if it were in a box with sides given by (in vertical extent) the eccentricity of the orbit and (in longitudinal and latitudinal extent) the inclination of the orbit. If r is the difference between the maximum orbit radius and the average radius = the difference between the minimum orbit radius and the average radius, then the total vertical movement is 2r. The total vertical motion 2r = rmax rmin = a(1 + e) a(1 e) = 2ae, where a is the average orbit radius and e is the eccentricity of the orbit. Thus 2r = 2 42,164.17 0.001 = 84.3283 km. Alternatively, using equation (10.4), r = eRav, where Rav = the average orbit radius, giving r = (0.001) 42,164.17 = 42.16417, and so the total vertical motion = 2 42.16417 = 84.3283 km. The lateral movement of the satellite N-S (longitude) and E-W (latitude) in its orbit is obtained by multiplying the orbit radius by the angular movement in radians. The angular movement is 0.05o = 0.1o total 0.0017453 radians and so the physical movement = 0.0017453 42,164.17 = 73.5903593 km. The movement box is shown below. E-W 73.5903593 km N-S 73.5903593 km Vertical 84.3283 km A

The center of the box is the nominal position of the satellite, with excursions up and down given by the eccentricity of the orbit and N-S and E-W given by the inclination of the orbit. Diagonal BC across the box is the maximum movement of the satellite (but note that there are also three other similar lengths within the box, created by joining the relevant opposing corners of the box.). Using Pythagorean geometry, we have BC2 = AC 2 + AB2 and AC 2 = 73.59035932 + 73.59035932 = 10,831.08196

2- 13

Thus BC 2 = 10,831.08196 + 84.32832 = 10,831.08196 + 7,111.262181 = 17,942.34414, which gives BC = 133.9490356 km = 133.95 km. The range variation observed by an earth station will depend on the precise coordinates of the earth station. While the elevation angle (5o ) was given, the earth station coordinates were not. At this elevation angle, the situation looks like that below.

5 elevation angle E C S

Station keeping box Local earth station horizon keeping box

C = Center of earth E = Earth station S = Satellite For an elevation angle of 5o , Angle CES = 95o . By the law of sines, we have sin ESC = 6378.137 sin 95 = 0.1506935 42164.17

Giving angle ESC = 8.6671179o The situation inside the station-keeping box is as below 84.3283 km 8.6671179o To earth Direction of view from earth station

2- 14

The maximum range variation in this case is 84.3283 (cos8.6671179o ) = 83.3653157 = 83.365 km.

Question 9 An interactive experiment is being set up between the University of York, England (approximately 359.5o E, 53.5o N) and the Technical University of Graz, Austria (approximately 15o E, 47.5o N) that will make use of a geostationary satellite. The earth stations at both universities are constrained to work only above elevation angles of 20o due to buildings, etc., near their locations. The groups at the two universities need to find a geostationary satellite that will be visible to both universities simultaneously, with both earth stations operating at, or above, an elevation angle of 20o . What is the range of subsatellite points between which the selected geostationary satellite must lie?

Solution to question 9 From equation (2.38), where El = elevation angle, cos( El) =

[1.02288235 0.30253825 cos( ) ]1 / 2

sin( )

Let cos() = X, C = cos (El), A = 1.02288235, B = 0.30253825, and so (1 X2 )1/2 = sin() The equation relating El to can be written as a quadratic in X and solved by the quadratic formula as follows: 1 X2 = C2 (A BX) (NOTE: this is the cos (El) equation written with the new parameters, and the subject of the equation changed) Solving gives C 2 B (C 4 B 2 4(C 2 A 1)) 1 / 2 X= 2 The correct solution is the positive root. With a minimum elevation angle of 20o , C = 0.9396926 giving X = 0.4721212, and thus = cos-1 (0.4721212) = 61.8279251o . The geometry of the linked pair of earth stations is as shown in the schematic below.

2- 15

York, UK latitude 53.5o N longitude 359.5o E

Graz, Austria latitude 47.5o N longitude 15o E

Range of sub-satellite points to provide a minimum of 20o elevation angle at both sites

Geostationary orbit plane

Elevation angle minimum of 20o from Graz

Elevation angle minimum of 20o from York

From equation (2.36) Cos() = cos(Le)cos(ls le) For York: cos(61.8279251) = cos(53.5) cos(ls le) giving cos(ls le) = 0.7937174. The difference between the York longitude and the sub-satellite point is cos-1 (0.7937174) = 37.4657212o . Since York is at a longitude of 359.5o E, there are two solutions for the placement of the satellite (one east of York and one west of York). Only the placement east of York will be visible to Graz and so the east most placement of the satellite as far as York is concerned is 36.9657212o E = 36.97o E. For Graz: cos(61.8279251) = cos(47.5) cos(ls le) giving cos(ls le) = 0.6988278. The difference between the Graz longitude and the sub-satellite point is cos-1 (0.6988278) = 45.6669685o . Since Graz is at a longitude of 15o E, there are two solutions for the placement of the satellite (one east of Graz and one west of Graz). Only the placement

2- 16

west of Graz will be visible to York and so the west most placement of the satellite as far as Graz is concerned is 329.3330315o E = 329.33o E The range of sub-satellite points that will provide an elevation angle of at least 20o to both York and Graz is from 329.33o E to 36.97o E.

Question 10 The state of Virginia may be represented roughly as a rectangle bounded by 39.5o N latitude on the north, 36.5o N latitude on the south, 76.0o W longitude on the east, and 86.3o W longitude on the west. If a geostationary satellite must be visible throughout Virginia at an elevation angle no lower than 20o , what is the range of longitudes within which the sub-satellite point of the satellite must lie?

Solution to question 10 This question is similar in concept to the previous question, and the same quadratic solution may be applied. The schematic below illustrates the geometry of the question. 86.3o W longitude State of Virginia 76o W longitude 39.5o N latitude

36.5o N latitude

Geostationary orbit plane Range of sub-satellite points giving a minimum of 20o elevation angles in Virginia

2- 17

Using the same quadratic equation in the solution to question 9 gives = 61.8279251o . The lowest elevation angles will be at earth stations located at the north most corners of Virginia. That is at latitude 39.5o N, with one at 76o W and the other at 86.3o W. The factor cos(Le) will be the same in both cases, cos(39.5) = 0.7716246, as will cos(ls le), which will be cos(61.8279251)/(0.7716246) = 0.6118535, which yields a separation between ls and le of 52.2763585o . Thus the satellite may be approximately 52.3o east or west of the earth station to remain at an elevation angle of 20o . However, to enable the other northern corner of Virginia to still see the satellite at an elevation angle of at least 20o , the satellite must be east of the west most earth station and west of the east most earth station. The east most northern earth station site is at 76o W and the west most earth station site is at 86.3o W. The subsatellite points must therefore be between 76o W + 52.2763585 = 128.2763585o W and 86.3o W 52.2763585 = 34.0236415o W

Question 11 A geostationary satellite system is being built which incorporates inter-satellite links (ISLs) between the satellites. This permits the transfer of information between two earth stations on the surface of the earth, which are not simultaneously visible to any single satellite in the system, by using the ISL equipment to link up the satellites. In this question, the effects of ray bending in the atmosphere may be ignored, processing delays on the satellites may initially be assumed to be zero, the earth may be assumed to be perfectly circular with a flat (i.e. not hilly) surface, and the velocity of the signals in free space (whether in the earths lower atmosphere or essentially in a vacuum) may be assumed to be the velocity of light in a vacuum. (i) What is the furthest apart two geostationary satellites may be so that they can still communicate with each other without the path between the two satellites being interrupted by the surface of the earth? Give your answer in degrees longitude between the sub-satellite points. (ii) If the longest, one-way delay permitted by the ITU between two earth stations communicating via a space system is 400 ms, what is the furthest apart two geostationary satellites may be before the transmission delay of the signal from one earth station to the other, when connected through the ISL system of the two satellites, equals 400 ms? The slant path distance between each earth station and the geostationary satellite it is communicating with may be assumed to be 40,000 km. (iii) If the satellites in part (ii) employ on-board processing, which adds an additional delay of 35 ms in each satellite, what is the maximum distance between the ISL- linked geostationary satellites now? (iv) If both of the two earth stations used in parts (ii) and (iii) must additionally now send the signals over a 2,500 km optical fiber line to the end-user on the ground, with an associated transmission delay in the fiber at each end of the link, what is the maximum distance between the ISL- linked geostationary satellites now? You may assume a refractive index of 1.5 for the optical fiber and zero processing delay in the earth station equipment and end- user equipment.

2- 18

Solution to question 11 A schematic of part (i) is shown below. SA = satellite A; SB = satellite B; C = Center of earth G = point where ISL path is tangential to earth CSA = GEO radius = 42,164.17 km CG = radius of earth = 6378.137 km SA SB earth-grazing ISL path between satellites

SA

G C

SB

(GSA)2 = (CSA)2 (GC)2 = 1,777,817,232 40,680,699.86 = 1,737136,562, Giving GS A = 41,678.97026 km By the law of sines, Sin(SACG) = (SAG Sin90o )/(CSA) = 41,678.97026/42,164.17 = 0.9884926, hence angle SACG = 81.2995147o . The difference in sub-satellite point between the two satellites with this geometry is twice angle SACG = 162.6o Part (ii) A schematic of part (ii) is shown below.

2- 19

SA = satellite A; SB = satellite B; C = Center of earth E1 = Earth station 1 and E2 = earth station 2 CSA = CSB = GEO radius = 42,164.17 km CE1 = CE2 = radius of earth = 6378.137 km SA SB ISL path between satellites

SA

SB

E1 C

E2

For a 400 ms delay and a signal velocity of 2.997 108 m/s, the total path length traversed must be velocity time = (2.997 108 m/s) (400 10-3 ) = 119,880,000 m = 119,880 km. E1 SA = E1 SB = 42,164.17 km 6378.137 = 35786.033 therefore SASB = 119,880 (35786.033 2) = 48,307.934 km. The question, however, stipulates that the slant path distances between the earth stations and their respective satellites (i.e. distances E1 SA and E1 SB) are 40,000 km. Thus, SASB = 119,880 80,000 = 39,880 km. The maximum physical separation between satellites to keep the end-to-end delay below 400 ms is 39,880 km, which translates into a separation between sub-satellite points of 56.4471538o = 56.45o .

2- 20

Part (iii) In this part, an additional total delay of 70 ms occurs in the satellites (35 ms at each satellite), which reduces the total path length for the free space part of the link to (2.997 108 m/s) ((400 70) 10-3 ) = 98,901,000 m = 98,901 km. With this reduced free space path length available, the maximum physical separation of the two satellites is 98,901 80,000 = 18,901 km, which translates into a separation between sub-satellite points of 25.9o . Part (iv) In this part, an additional transmission line exists at both ends of the link, totaling 5,000 km. The velocity of light in the fiber link is (2.997 108 m/s)/refractive index = 1.998 108 m/s, which leads to an additional delay due to the fiber links of 5,000,000/1.998 108 = 0.0250250 s = 25 10-3 s. The available free space path length is therefore reduced to (2.997 108 m/s) ((400 70 25) 10-3 ) = 91,408,500 m = 91,408.5 km. With the on board processing delay and extended fiber optic transmission line at each end, the maximum physical separation of the two satellites is 91,408.5 80,000 = 11,408.5 km, which translates into a separation between sub-satellite points of 7.78o . Summarizing, we have found the angular separation of satellites connected via an Inter Satellite Link to be as follows: (i) Maximum separation = 41,678.9702683,357.9 km 162.6o between satellites (ii) With a free space delay of 400 ms, separation = 39,880 km 56.45o between satellites (iii) With a free space delay of 400 ms and 70 ms of total on board processing delay, separation = 18,901 km 25.9o between satellites (iv) With a free space delay of 400 ms, 70 ms of total on board processing delay, and 25 ms total back haul delay at each end of the link, separation = 11,408.5 km 7.7o . A 2,500 km back haul link is a typical design requirement, thus on board processing delays need to be kept to a minimum if ISLs are to operate successfully in real time links.

2- 21