Elevation Angle Characterization for LEO Satellites: First and Second Order Statistics

Abstract

The elevation angle $\theta$ is relevant for the Low Earth orbit (LEO) satellite communications since it is always changing its relative position with respect to fixed Earth stations (ES’s), and this affects the link length and received power, $P_R$. This article provides a new methodology to compute the probability density function (PDF) and cumulative distribution function (CDF) of the elevation angle, $\theta$, for diverse ES locations. This methodology requires as input parameters an ES latitude, $\varphi$, an orbit inclination value, i, and an orbit altitude, h. The elevation angle is characterized through a well known random variable, which facilitates the computation of the first and second-order statistics, and helps to determine the expected value and measures of dispersion of the angle $\theta$ for a particular ES location. The proposed methodology allows an easy and quick calculation of the elevation angle’s CDF, facilitating comparisons against CDF’s of more ES’s located at different latitudes, and longitudes, $\lambda$; as well as the comparisons of CDF’s of the elevation angle produced by different orbits. Extensive simulation results are summarized in a small table, which allows computation of the elevation angle’s CDF and PDF for multiple ES locations without requiring of simulations and statistical fitting. Finally, the proposed methodology is validated through an extensive error analysis that show the suitability of the obtained results to characterize the elevation angle.

1. Introduction

The characterization of the elevation angle, $\theta$ is relevant for low Earth orbit (LEO) satellite communications since this parameter is directly related to the varying distance between the satellite and Earth station (ES), and affects the link total attenuation. The elevation angle description through an analytical expression is a difficult problem addressed in [1,2] which has received less attention in the literature.

Nonetheless, this parameter is directly related to the link performance and channel characterization of LEO satellites. The LEO channel characterization has also received less attention than the geostationary (GEO) satellite channel and just few models such as [3] have been specifically developed and published to consider the elevation angle variations introduced at LEO. The lack of LEO channel models accounting for the always-changing elevation angle have resulted in just a few channel models available in the literature for LEO satellites and specially for small satellites as mentioned in [4].

LEO satellites are increasing their numbers and role as an enabling technology for Internet of Things (IoT), 5G, 6G, and next generation wireless networks aimed to provide global coverage with very low latency [5]. Then it is relevant to develop methodologies to analyze and compare the link performance and channel characteristics considering the variations introduced by the always-changing elevation angle.

The always-changing elevation angle condition has been a limitation to analyze the link budget and channel of LEO satellites, and common approaches to characterize it have followed segmentation in best and worst-case of the elevation angle [6], instead of analyzing the short and long term behavior of the elevation angle as an analytical function. However, some emerging problems, such as efficient power management [7], related to LEO satellites have made evident the necessity of having a way to characterize the elevation angle as an analytical function.

The observed elevation angle, $\theta$, is always-varying for LEO satellites and can be described as a function of time, $\theta \left(t\right)$, for intervals in which a satellite is visible. The calculation of $\theta$ for a given position can be deterministically obtained with mathematical procedures developed in orbital mechanics books [8] and widely implemented in software for space dynamics simulations. However, it is important to note that even with accurate predictions of $\theta$ (accounting for contacts in several days or months), the relations between one contact and another, as well as the long term behavior of the elevation angle can be hardly described without randomness. Thus, appropriate use of probability theory becomes a great tool to analyze the behavior of the elevation angle, $\theta$, in the long term.

The elevation angle, $\theta$, is usually defined for an ES as the angle (above the local horizon) at which the satellite is visible, and within the interval of 0° to 90°, $\theta \in (0^{°},{90}^{°}]$, regardless of its azimuth. The minimum value of $\theta$ at which communication is possible is often called ${\theta }_{min}$ (subject to ${\theta }_{min}\ge 0^{°}$), similarly, the maximum value of $\theta$ that can be observed from the ES is called ${\theta }_{max}$ (subject to ${\theta }_{min}\le {\theta }_{max}\le {90}^{°}$). We can define a random variable, r.v., to take possible values of $\theta$, such that $\theta \in ({\theta }_{min},{\theta }_{max}]$ as $\mathsf{\Theta }$.

The probability density function, PDF, and cumulative distribution function, CDF, of the elevation angle, $f_{\mathsf{\Theta }}\left(\theta \right)$ and $F_{\mathsf{\Theta }}\left(\theta \right)$, respectively, are useful functions to characterize the LEO channel since the elevation angle affects the received power level, $P_R$. The reasoning behind that is that $\theta$ depends on the distance from the satellite to the ES, $r_{S,E}$. Then, at greater link distances, it can be expected to have a lower received power, $P_R$, (of the transmitted signal from the satellite) than at shorter distances.

The distance between a LEO satellite and an ES, $r_{S,E}$, can vary several thousands of kilometers from a low value of $\theta$ to a high value of $\theta$. Figure 1a shows the extreme cases for the elevation angle in a LEO satellite link. Those extreme cases are not necessarily met at every contact, but represent the best and worst length-scenarios in a long period. Path lengths between the ES and satellite, $r_{S,E}$, are described in this figure using the variable $r_i,i\in \{1,2,3\}$. The figure also shows that at low values of $\theta$ the path length is larger and the received power, $P_R$ is at its minimum value. This last implication between the elevation angle, $\theta$, the received power, $P_R$, and the link length, $r_{S,E}$, assumes some constraints explained in detail in [9].

Figure 1b shows the typical-case for the elevation angle, where the value of $\theta$ is not necessarily at its minimum nor at its maximum, but it can be at any value within the range $\theta \in [{\theta }_{min},{\theta }_{max}]$. Since the elevation angle of a LEO satellite appears as always-varying from an ES’s, it is convenient to determine statistical indicators of its behavior, such as its expected value, $E[·]$, median, $MED[·]$, standard deviation, $SD[·]$; and how often does the elevation angle will be above or below a threshold, or within a region of interest, for example, using its quantiles $nQ$. In addition, since the elevation angle variations are directly related to the link length, the elevation angle for LEO satellites is also directly related to the variations of the received power $P_R$ at an ES, as described in [9].

1.1. Contributions

The elevation angle characterization is relevant for LEO satellites because it is an always-changing variable for those communication systems. The effects of the elevation angle variations are observable in the received power, link quality, and channel behavior; then, an accurate characterization of this variable is a topic of interest for planning and implementation of LEO satellite systems.

In this article, we have addressed the characterization of the elevation angle for LEO satellites, by obtaining its PDF and CDF, as well as the derivation of its first and second order statistics. First and second order statistics are relevant to evaluate the suitability of different LEO orbits and to determine which orbits are more convenient to provide coverage for a particular application or Earth station location.

In addition, we have developed an extensive analysis to validate our results, and we include supplementary materials containing elevation angle times series. The supplementary materials will facilitate reproducibility of our work, and will allow future research based on our proposed methodology and results.

1.1.1. Contribution 1

This document shows the feasibility of using a random variable to characterize the elevation angle behavior for LEO satellites with different orbit configurations. The suitability of using the proposed random variable is verified through an extensive error analysis with diverse orbit configurations.

1.1.2. Contribution 2

The PDF and CDF of the elevation angle can be obtained as proposed in [1,2]. Nonetheless, this document describes a methodology to obtain the PDF and CDF parameters of the elevation angle distribution for different orbit configurations and ES locations using a well-known random variable. The proposed random variable facilitates analytical manipulation as well as computation of the probabilities of occurrence of the elevation angle at specific values. A small table containing the resultant parameters to characterize multiple LEO orbits configurations as observed from multiple ES’s is included.

1.1.3. Contribution 3

The expected value of the elevation angle and its measures of dispersion are relevant for planning and implementation of LEO satellite systems, since the received power varies according to the elevation angle. This document describes an analytical methodology to obtain the first and second order statistics of the elevation angle for different orbit configurations and ES locations, which, to the best of our knowledge are not available in the literature and were just utilized for a specific case in [9].

1.2. Outline

The remaining of the article is organized as follows: Section 2 introduces the theoretical fundamentals required to understand the subsequent sections; Section 3 describes the developed methodology to characterize the elevation angle behavior through $f_{\mathsf{\Theta }}\left(\theta \right)$ and $F_{\mathsf{\Theta }}\left(\theta \right)$; Section 4 contains the main results; and Section 5 concludes analyzing the obtained results and opportunities for future work. Figure 2 shows the main structure of the document.

2. Materials and Methods

2.1. Elevation Angle Definition

The elevation angle, $\theta$, for an ES can be defined as the angle between the local horizon and the satellite. Several sources such as [10,11,12] contain expressions to calculate $\theta$ from geometrical relations between the satellite and the ES instantaneous positions. One of those definitions is as follows

$$\theta =arctan\left({\displaystyle \frac{cos\Delta cos{\varphi }_{\mathrm{ES}}-(r_{\mathrm{E},\mathrm{O}}/r_{\mathrm{S},\mathrm{O}})}{\sqrt{1-{cos}^2\Delta {cos}^2{\varphi }_{\mathrm{ES}}}}}\right),$$

(1)

where $r_{\mathrm{E},\mathrm{O}}$ and $r_{\mathrm{S},\mathrm{O}}$ are the distances from the center of the Earth, O, to the ES and the satellite, respectively; ${\varphi }_{\mathrm{ES}}$ is the latitude at which the ES is located (in degrees); M is the subsatellite point, which corresponds to the latitude and longitude of the satellite instantaneous position; and $\Delta$ is the difference in longitude between the ES and M (in degrees). Figure 3a shows a derivation of $\theta$ based on (1).

2.2. Relevance of the Elevation Angle for the Satellite Channel

LEO satellites differ to those at GEO in its relative position as seen from a fixed ES. Whereas GEO satellites appear as fixed, those at LEO appear as always-moving points in the sky or are absent causing link unavailability. Link availability can be seen as an ON-OFF process, where the link is ON if the satellite is visible in the sky, and it is OFF if the satellite is absent. Figure 3b shows the ON-OFF process for both a LEO and GEO satellite.

In addition to the short visibility that LEO satellites have, conditions in this lapse are always varying as a consequence of the changing position of the satellite. When the satellite is first visible from an ES, it is at the largest distance supported by the link, then, the satellite starts approaching the ES until reaching the shortest distance in that particular contact. Finally, the satellite moves away after being at its minimum contact distance. It is important to note that the minimum distance from an ES to a LEO satellite, $min\left(r_{E,S}\right)$, varies from one contact to another, as well as $max\left(r_{E,S}\right)$ does.

Differences between the Earth’s rotation rate and LEO orbiting velocities cause an always-varying link length (as seen from a fixed ES), and then, an always-varying elevation angle.

An important relation that needs to be accounted for the elevation angle, $\theta$, and the link length, $r_{S,E}$ is given by the implication that when the angle is minimum, the $r_{S,E}$ distance is maximum, i.e.,

$$\theta \to min\left(\theta \right)⟺r_{S,E}\to max\left(r_{S,E}\right),$$

(2)

similarly

$$\theta \to max\left(\theta \right)⟺r_{S,E}\to min\left(r_{S,E}\right);$$

(3)

where the maximum link length, $max\left(r_{S,E}\right)$, occurs at some value of $\theta$ such that $min\left(\theta \right)<\theta \le {90}^{°}$; and $min\left(r_{S,E}\right)$ occurs at some value of $\theta \ge 0^{°}$ lower than ${\theta }_{max}$. The implications in (2) indicate that when values of $\theta$ diminish, $r_{S,E}$ increases; similarly, in (3), when values of $\theta$ increase then $r_{S,E}$ decreases.

As a satellite moves away from an ES, $r_{S,E}$ and the free-space path loss, $L_{FS}$, increases. Longer paths occur for lower values of $\theta$, then, the atmospheric attenuation, $A_{Atm}$, increases for larger paths (assuming similar atmospheric conditions for different link paths). On the other hand, when $\theta$ takes values close to ${90}^{°}$, $r_{S,E}$ will be at its minimum, and so does $A_{Atm}$. Then, low elevation angles cause a greater $L_{FS}$, and more atmospheric attenuation, $A_{atm}$ (assuming that the atmospheric conditions are approximately the same for two distinct paths).

For LEO it is important to consider link interference, I, usually addressed as noise in link-budget calculations, and also dependent on the elevation angle as mentioned in [13]. Furthermore, low values of $\theta$ are also associated in practice with non-line-of-sight (NLOS) conditions, increasing ground interference [14], and increasing multipath fading in the land mobile satellite (LMS) channel [15,16,17].

There are several models available in the literature to characterize the received signal from satellite systems. Most of these works have been elaborated for the LMS channel, which describes the received signal at a land-moving ES; nonetheless, most of those works focus on GEO systems. Extensive reviews for the LMS channel can be found at [16,18,19].

Among the available channel models for LMS systems, a few of them focus on characterizations for LEO and non-GEO (NGEO) satellites; but, as noted by [4,20], those are few and there are still many challenges. Most of channels for NGEO systems have been developed based on previous models for GEO satellites, specially on well-known models, such as those by Loo [21] and Lutz [22].

Even though the Loo and Lutz models were not originally intended for LEO satellites, those were later used as a base for much of the non-geosynchronous (NGSO) channel models. The Loo channel assumes that the complex envelope of the received signal, $r_T$, contains a LOS, $r_D$, and a multipath component, $r_M$, which in its phasor notation can be represented as

$$r_{T}exp\left(j{\varphi }_T\right)=r_{D}exp\left(j{\varphi }_D\right)+r_{M}exp\left(j{\varphi }_M\right)$$

(4)

where ${\varphi }_T$, ${\varphi }_D$, and ${\varphi }_M$ indicate the phase of the total received signal, of the direct or LOS component, and of the multipath component, respectively.

The channel developed by Lutz [22] models a varying received signal affected by different levels of shadowing. This channel is characterized by a Markov chain with a state corresponding to light shadowing, G, and another corresponding to deep shadowing conditions B. Figure 4 illustrates the Lutz channel.

The Loo channel [21] was adapted by Corazza and Vatalaro [3] to formulate one of the best-known channel models for LEO satellites; similarly, it was modified by Abdi et al. [23], to achieve straightforward analytical expressions for the first and second order statistics of the received signal based on the Nakagami distribution. However, current models for GEO and NGEO systems have recognized the limitations of a single distribution and a single state to describe the received signal, and mutistate models with different distributions at each state have became popular.

Channel models for LEO with multiple states as in [22] have been developed containing the same kind of distribution at each state. For example, a combination of the Lutz and Loo models was implemented for NGSO by Perez-Fontan et al. [24], focusing on time series generation. This model used a three-state Markov chain, as shown in Figure 5, to indicate different levels of shadowing; the first state, $S_1$, indicates deep shadow; the second, $S_2$, indicates moderate shadowing; and the third state, $S_3$, indicates line-of-sight (LOS) conditions. Each state describes the received signal according to a Loo distribution with a different Loo triplet, $\alpha$, $\psi$, $MP$; and a different Markov chain for distinct elevation angle values. The so-called Loo triplet, which includes the mean of the direct signal amplitude, $\alpha$, the standard deviation of the direct signal amplitude, $\psi$, and the multipath power, $MP$, is the set of parameters required by the Loo distribution.

Figure 6 shows two Loo distributions to characterize the satellite channel at two elevation angles; the CDF’s were computed assuming presence of a LOS component, and a multipath power component very close to the LOS level. The Loo distributions in Figure 6 can be understood as two of the states of Figure 5 at different elevation angles; e.g., a state transition from $S_{3,{\theta }_1}$ to $S_{3,{\theta }_2}$.

In addition to the shadowing level at a receiver environment, which is mainly determined by natural and human made objects at the surroundings of the ES, the received signal level will be changing according to the elevation angle value. For LEO satellites, the elevation angle will change rapidly, and with a different rate at every contact.

As observed in Figure 6, the Loo CDF curves of the received signal rely on the elevation angle value since those depend on the LOS component of the received signal, $\alpha$. Then, in order to predict different curves we need to know the expected value and measures of dispersion of the elevation angle.

Even though we illustrate the received signal through the well-known Loo model in Figure 6, there are more LEO channel models depending on a line-of-sight component, and thus, directly depending on the elevation angle. The, elevation angle characterization proposed in this article can be a valuable tool to determine CDF curves of the received signal, based on the first and second order order statistics of the elevation angle. For example, to determine the expected CDF of the Loo model based on the elevation angle expected value, and determine how far apart other curves will be based on the measures of dispersion of the elevation angle.

2.3. Analytical Characterization of the Elevation Angle

The elevation angle was analytically characterized in [2] through its probability density function (PDF), $f_{\mathsf{\Theta }}$, which is defined as a marginal distribution from the joint PDF of $\mathsf{\Theta }$ and the maximum value of the elevation angle, ${\mathsf{\Theta }}_{max}$, as follows

$$f_{\mathsf{\Theta }}\left(\theta \right)={\int }_{\theta }^{{\theta }_M}{f}_{\mathsf{\Theta },{\mathsf{\Theta }}_{max}}(\theta ,{\theta }_{max})d{\theta }_{max}$$

(5)

where $f_{\mathsf{\Theta },{\mathsf{\Theta }}_{max}}$ is given by

$$f_{\mathsf{\Theta },{\mathsf{\Theta }}_{max}}(\theta ,{\theta }_{max})=f_{\mathsf{\Theta }|{\mathsf{\Theta }}_{max}}\left(\theta \right|{\theta }_{max})f_{{\mathsf{\Theta }}_{max}}\left({\theta }_{max}\right)$$

(6)

or equivalently

$$f_{\mathsf{\Theta },{\mathsf{\Theta }}_{max}}(\theta ,{\theta }_{max})={\displaystyle \frac{G\left(\theta \right)sin\gamma \left(\theta \right)}{\sqrt{{cos}^2\gamma \left({\theta }_{max}\right)-{cos}^2\gamma \left(\theta \right)}}}·{\displaystyle \frac{f_{{\mathsf{\Theta }}_{max}}\left({\theta }_{max}\right)}{{\int }_{{\theta }_{min}}^{{\theta }_M}{f}_{{\mathsf{\Theta }}_{max}}\left(x\right){cos}^{-1}\left(\frac{cos\gamma \left({\theta }_{min}\right)}{cos\gamma \left(x\right)}\right)dx}}$$

(7)

The integral dividing the right hand term is a constant, say $C_1$, then, (7) can be rewritten as

$$f_{\mathsf{\Theta },{\mathsf{\Theta }}_{max}}(\theta ,{\theta }_{max})={\displaystyle \frac{f_{{\mathsf{\Theta }}_{max}}\left({\theta }_{max}\right)G\left(\theta \right)sin\gamma \left(\theta \right)}{C_1\sqrt{{cos}^2\gamma \left({\theta }_{max}\right)-{cos}^2\gamma \left(\theta \right)}}}$$

(8)

and the auxiliary functions $\gamma (·)$ and $G(·)$ require an extra parameter a defined as $\alpha =r_E/r_S$, where $r_E=6378$ km is the radius of the Earth; and $r_S=r_E+h$ depends on the altitude h (from the Earth’s surface) of the circular LEO orbit in kilometers. The functions $\gamma (·)$ and $G(·)$ are defined as follows

$$\gamma \left(\theta \right)={cos}^{-1}(\alpha cos\theta )-\theta$$

(9)

$$G\left(\theta \right)=\frac{1+{\alpha }^2-2\alpha cos\gamma \left(\theta \right)}{1-\alpha cos\gamma \left(\theta \right)}$$

(10)

The term $f_{{\mathsf{\Theta }}_{max}}\left({\theta }_{max}\right)$ is defined as follows

$$\begin{array}{c}\hfill f_{{\mathsf{\Theta }}_{max}}\left({\theta }_{max}\right)=\frac{G\left({\theta }_{max}\right)}{K_2}·f_{\mathsf{\Phi }}({\varphi }_0-\gamma \left({\theta }_{max}\right)),\mathrm{for}{\theta }_{min}\le {\theta }_{max}<{\theta }_c\end{array}$$

(11)

$$f_{{\mathsf{\Theta }}_{max}}\left({\theta }_{max}\right)=\frac{G\left({\theta }_{max}\right)}{K_2}·\left[f_{\mathsf{\Phi }}\right({\varphi }_0-\gamma \left({\theta }_{max}\right)+f_{\mathsf{\Phi }}({\varphi }_0+\gamma \left({\theta }_{max}\right))],\mathrm{for}{\theta }_c<{\theta }_{max}<{\theta }_M\hfill$$

(12)

where $f_{\mathsf{\Phi }}\left(\varphi \right)$ is defined from the orbital inclination i; and ES latitude, $\varphi$, as follows

$$f_{\mathsf{\Phi }}\left(\varphi \right)=\frac{cos\varphi }{\pi \sqrt{{sin}^{2}i-{sin}^2\varphi }},\mathrm{for}\left|\varphi \right|<i$$

(13)

and ${\theta }_M=\pi /2$, and $K_2$ is given by

$$K_2=\frac{1}{2}-\frac{1}{\pi }{sin}^{-1}\left(\frac{sin(\varphi +\gamma \left({\theta }_{min}\right))}{sin\left(i\right)}\right)$$

(14)

where ${\theta }_{min}$ is the minimum elevation angle that will be considered and is in the range $0^{°}<{\theta }_{min}<{90}^{°}$.

Although results shown in [2] are highly accurate with respect to the actual elevation angle PDF, the analytical expression $f_{\mathsf{\Theta }}$ is cumbersome to obtain and evaluate. Furthermore, analytical expressions for second or higher order statistics, are a missing result in the literature to the best of our knowledge.

3. Methodology for the Elevation Angle Characterization

In this section, we describe the procedure to characterize the elevation angle behavior through a random variable in a manageable analytical expression. Also, we include a procedure to obtain useful expressions of the elevation angle, such as measures of central tendency and dispersion.

The elevation angle characterization for different LEO orbit’s configuration will help to predict which orbits are more convenient for a particular ES location. In addition, the measures of central tendency and dispersion will help to choose suitable orbits based on the expected elevation angle, and to calculate the link attenuation based on the distribution of the elevation angle, and on its expected value and variance.

3.1. Simulations Configuration

Simulations were performed for a LEO satellite with characteristics mentioned in

Table 1. Satellite characteristics for the ephemeris generation.

Dry Mass	Drag Area	Solar Radiation Pressure Area
5 kg	1 m${}^2$	1 m${}^2$

, and for circular orbit configurations listed in

Table 2. Orbit characteristics for the ephemeris generation.

	Semi-Major Axis, a	Orbit Inclinations, i
1	7378 km	20${}^{°}$, 25${}^{°}$, …, 85${}^{°}$
2	7578 km	20${}^{°}$, 25${}^{°}$, …, 85${}^{°}$
3	7778 km	20${}^{°}$, 25${}^{°}$, …, 85${}^{°}$
4	7978 km	20${}^{°}$, 25${}^{°}$, …, 85${}^{°}$
5	8178 km	20${}^{°}$, 25${}^{°}$, …, 85${}^{°}$
6	8378 km	20${}^{°}$, 25${}^{°}$, …, 85${}^{°}$

. The LEO region goes from a few hundred kilometers above the Earth’s surface up to 2000 km; but the performed simulations and methodology cover just the upper part of this region since orbit perturbations effects are much lower at those altitudes than at orbital heights closer to the surface, and facilitates operation for longer periods of time without a complex propulsion system.

Three initial values of orbital parameters not included in Table 2 are right ascension of the ascending node, $\mathsf{\Omega }$, argument of perigee, $\omega$, and true anomaly, $\nu$, which were all set to zero. From those simulations, ephemeris files with the satellite position and velocity in five-second steps were generated. Then, ES’s were placed (simulated) at 18 different and arbitrarily chosen locations listed in

Table 3. ES locations for the performed simulations specified through their latitude, $\varphi$, and longitude, $\lambda$.

	Latitude	Longitude		Latitude	Longitude
ES${}_1$	0${}^{°}$	276.7121${}^{°}$	ES${}_{10}$	45${}^{°}$	259.7121${}^{°}$
ES${}_2$	5${}^{°}$	276.7121${}^{°}$	ES${}_{11}$	50${}^{°}$	259.7121${}^{°}$
ES${}_3$	10${}^{°}$	276.7121${}^{°}$	ES${}_{12}$	55${}^{°}$	259.7121${}^{°}$
ES${}_4$	15${}^{°}$	267.7121${}^{°}$	ES${}_{13}$	60${}^{°}$	259.7121${}^{°}$
ES${}_5$	20${}^{°}$	276.7121${}^{°}$	ES${}_{14}$	65${}^{°}$	259.7121${}^{°}$
ES${}_6$	25${}^{°}$	261.7121${}^{°}$	ES${}_{15}$	70${}^{°}$	265.7121${}^{°}$
ES${}_7$	30${}^{°}$	259.7121${}^{°}$	ES${}_{16}$	75${}^{°}$	265.7121${}^{°}$
ES${}_8$	35${}^{°}$	259.7121${}^{°}$	ES${}_{17}$	80${}^{°}$	265.7121${}^{°}$
ES${}_9$	40${}^{°}$	259.7121${}^{°}$	ES${}_{18}$	−85${}^{°}$	259.7121${}^{°}$

. Finally, time series of the observed elevation angle from the ES to the satellite were calculated for each ground location using the previous generated ephemeris files and geometrical relations as in (1).

The orbital mechanical equations required to calculate the position and velocity of the satellite as a function of time, were solved using an open source program developed and maintained by the NASA, GMAT [25]. The time series for the elevation angle discussed in this paper can be generated in GMAT using the data of Table 1, Table 2 and Table 3; also, these can be generated using other software, for example STK as in [2]. The time series mentioned in this paper were generated by the authors specially for this work, and are available at [26].

3.2. Elevation Angle Time Series Analysis

Each elevation angle data set (obtained from the simulations) was individually analyzed to observe its statistical values of interest (mean, median, standard deviation, quantiles). Those time series were fitted to several distributions using the maximum likelihood estimation method. The goodness of fit of the proposed distributions was evaluated for each elevation angle time series using the Kolmogorov-Smirnov test, and it was observed a better performance of the gamma distribution, $Gamma(a,b)$, for most of the data sets. Other distributions that showed good fit based on the Kolmogorov-Smirnov test were the beta, $B(a,b)$, and Weibull, $Weib(a,b)$. Figure 7 shows the values of $\theta$ as well as the PDF and CDF for one elevation angle time series, $\theta \left(t\right)$ using the proposed distributions. Figure 8 shows the Kolmogorov-Smirnov statistic, often indicated as D, for the time series obtained for different altitudes. The Kolmogorov-Smirnov statistic is an indicator of the maximum distance between the actual time series distribution and the proposed gamma distribution. A low value of the Kolmogorov-Smirnov is often desirable to indicate a better goodness-of-fit.

The gamma PDF, $f_{\mathsf{\Theta }}\left(\theta \right)$, and CDF, $F_{\mathsf{\Theta }}\left(\theta \right)$, are defined as follows

$$f_{\mathsf{\Theta }}\left(\theta \right)={\displaystyle \frac{1}{b^a\mathsf{\Gamma }\left(a\right)}}{\theta }^{a-1}exp(-\theta /b),$$

(15)

and

$$F_{\mathsf{\Theta }}\left(\theta \right)={\displaystyle \frac{1}{\mathsf{\Gamma }\left(a\right)}}\gamma \left(a,{\displaystyle \frac{\theta }{b}}\right),$$

(16)

respectively, where $\mathsf{\Gamma }(·)$ is the gamma function, $\gamma (·)$ is the incomplete gamma function, a is the shape parameter, and b is the scale parameter. By definition both a and b are greater than zero.

The shape, a, and scale, b, parameters for the gamma distribution were obtained for each of the elevation angle time series of each ES. Those parameters were organized into several matrix arrays as shown in

Table 4. Rate parameter, a, of the gamma distribution to characterize the elevation angle distribution at ES${}_1$, for different altitudes, h, and orbit inclinations, i.

ES${}_1$	20${}^{°}$	25${}^{°}$	30${}^{°}$	35${}^{°}$	40${}^{°}$	45${}^{°}$	50${}^{°}$	55${}^{°}$	60${}^{°}$	65${}^{°}$	70${}^{°}$	75${}^{°}$	80${}^{°}$	85${}^{°}$
7378 km	1.39	1.50	1.48	1.31	1.30	1.34	1.35	1.35	1.35	1.37	1.37	1.36	1.36	1.37
7578 km	1.62	1.79	1.80	1.65	1.56	1.61	1.62	1.63	1.63	1.64	1.64	1.64	1.64	1.64
7778 km	1.57	1.77	1.82	1.73	1.55	1.61	1.63	1.64	1.64	1.65	1.65	1.66	1.66	1.66
7978 km	1.63	1.74	1.83	1.79	1.62	1.60	1.63	1.65	1.65	1.66	1.66	1.66	1.66	1.67
8178 km	1.70	1.66	1.83	1.82	1.70	1.55	1.63	1.65	1.66	1.67	1.67	1.67	1.70	1.67
8378 km	1.75	1.67	1.82	1.84	1.75	1.58	1.63	1.66	1.67	1.68	1.68	1.68	1.68	1.68

, which contains the shape parameter for the simulated orbits as seen from the ES${}_1$, and as in

Table 5. Scale parameter of the gamma distribution, $b\times {10}^{-2}$, to characterize the elevation angle distribution at ES${}_1$, for different altitudes, h, and orbit inclinations, i.

ES${}_1$	20${}^{°}$	25${}^{°}$	30${}^{°}$	35${}^{°}$	40${}^{°}$	45${}^{°}$	50${}^{°}$	55${}^{°}$	60${}^{°}$	65${}^{°}$	70${}^{°}$	75${}^{°}$	80${}^{°}$	85${}^{°}$
7378 km	7.51	8.71	9.34	8.88	8.40	8.42	8.40	8.36	8.34	8.39	8.38	8.35	8.35	8.36
7578 km	7.86	9.13	9.91	9.74	9.05	9.02	8.99	8.97	8.93	8.95	8.94	8.93	8.91	8.91
7778 km	7.45	8.55	9.41	9.56	8.83	8.73	8.70	8.67	8.64	8.66	8.64	8.64	8.63	8.64
7978 km	7.41	8.08	8.98	9.33	8.88	8.48	8.46	8.44	8.40	8.41	8.39	8.39	8.38	8.39
8178 km	7.38	7.59	8.60	9.07	8.88	8.18	8.26	8.24	8.22	8.22	8.20	8.18	8.26	8.18
8378 km	7.34	7.47	8.27	8.82	8.83	8.21	8.09	8.08	8.06	8.05	8.03	8.02	8.00	8.02

, which contains the scale parameter also for ES${}_1$.

3.3. Orbit Coverage

Figure 9a illustrates (without scale) passes of a satellite in each one of the six altitudes and for all the inclinations mentioned in Table 2. This figure also shows dashed lines indicating the ES’s latitudes mentioned in Table 3.

In Figure 9a as in Table 2 the orbit inclination parameter, i, for each altitude goes from ${20}^{°}$ to ${85}^{°}$, then, this figure can be redrawn as the grid shown in Figure 9b, where the curved lines of Figure 9a are replaced by vertical lines indicating the orbit inclination. Additionally, the ES latitude lines were inverted to start with the lowest value at the top.

From Figure 9a it can be observed that satellites with certain orbit inclination i, reach (in its orbit trajectory) at most the latitude that coincide with the inclination value. For example, a satellite with $i={20}^{°}$ will not be able to pass over and ES located at $\varphi ={80}^{°}$, since the maximum reached latitude in that orbit will be $\varphi \approx i\approx {20}^{°}$, then, that satellite will be best suited to provide coverage at ES’s located at latitudes $\varphi \le i$ or at most $\varphi \approx i$. This coincides with basic knowledge of satellite coverage; whereas polar orbits can cover almost all the globe, lower inclination orbits cover smaller portions of the Earth. Then, we can redraw Figure 9b to have a coverage grid by discarding points corresponding to cases of ${\varphi }_{ES}>i_{SAT}$, this resultant grid is shown in Figure 9c.

From the simulation results obtained for the grid points of Figure 9c, matrices ${\mathit{A}}_{\mathit{h}}$ and ${\mathit{B}}_{\mathit{h}}$ containing the shape and scale parameters of $f_{\mathsf{\Theta }}\left(\theta \right)$ and $F_{\mathsf{\Theta }}\left(\theta \right)$ can be defined as follows

$${\mathit{A}}_{\mathit{h}}={\left[\begin{array}{ccccc}{a}_{1,1}& a_{1,2}& a_{1,3}& \dots & a_{1,14}\\ a_{2,1}& a_{2,2}& a_{2,3}& \ddots & a_{2,14}\\ a_{3,1}& a_{3,2}& a_{3,3}& \ddots & ⋮\\ a_{4,1}& a_{4,2}& a_{4,3}& \ddots & ⋮\\ a_{5,1}& a_{5,2}& a_{5,3}& \ddots & ⋮\\ NaN& a_{6,2}& a_{6,3}& \ddots & ⋮\\ NaN& NaN& a_{7,3}& \ddots & ⋮\\ ⋮& ⋮& ⋮& \ddots & ⋮\\ NaN& NaN& NaN& \dots & a_{18,14}\end{array}\right]}_{18\times 14},{\mathit{B}}_{\mathit{h}}={\left[\begin{array}{ccccc}{b}_{1,1}& b_{1,2}& b_{1,3}& \dots & b_{1,14}\\ b_{2,1}& b_{2,2}& b_{2,3}& \ddots & b_{2,14}\\ b_{3,1}& b_{3,2}& b_{3,3}& \ddots & ⋮\\ b_{4,1}& b_{4,2}& b_{4,3}& \ddots & ⋮\\ b_{5,1}& b_{5,2}& b_{5,3}& \ddots & ⋮\\ NaN& b_{6,2}& b_{6,3}& \ddots & ⋮\\ NaN& NaN& b_{7,3}& \ddots & ⋮\\ ⋮& ⋮& ⋮& \ddots & ⋮\\ NaN& NaN& NaN& \dots & b_{18,14}\end{array}\right]}_{18\times 14},$$

(17)

where $NaN$’s are placed for the case ${\varphi }_{ES}>i_{SAT}$, and $h\in h_s$ corresponds to the simulation altitudes set defined in Table 2 as

$$h_s=\{1000\mathrm{km},1200\mathrm{km},\dots ,2000\mathrm{km}\}.$$

(18)

Table A1. Shape, a, and rate, b parameter for $f_{\mathsf{\Theta }}\left(\theta \right)$ and $F_{\mathsf{\Theta }}\left(\theta \right)$ for $h_k=1000$ km.

Shape parameter, a
		Orbit inclination, i
		$20{}^{°}$	$25{}^{°}$	$30{}^{°}$	$35{}^{°}$	$40{}^{°}$	$45{}^{°}$	$50{}^{°}$	$55{}^{°}$	$60{}^{°}$	$65{}^{°}$	$70{}^{°}$	$75{}^{°}$	$80{}^{°}$	$85{}^{°}$
ES latitude, $\varphi$	$0{}^{°}$	1.6063	1.5846	1.4659	1.4345	1.4584	1.4742	1.4817	1.4854	1.4864	1.4976	1.4981	1.4990	1.4954	1.5031
	$5{}^{°}$	1.7662	1.4853	1.4151	1.4548	1.4667	1.4770	1.4859	1.4898	1.4966	1.4968	1.4978	1.5015	1.5017	1.5040
	$10{}^{°}$	1.5036	1.6270	1.6291	1.4689	1.4303	1.4646	1.4756	1.4785	1.4779	1.4956	1.4867	1.4936	1.4935	1.4983
	$15{}^{°}$	1.3946	1.5570	1.6595	1.6390	1.4649	1.4340	1.4618	1.4690	1.4724	1.4938	1.4921	1.4914	1.4881	1.4920
	$20{}^{°}$	1.3749	1.4220	1.5688	1.6689	1.6382	1.4680	1.4307	1.4545	1.4573	1.4891	1.4848	1.4851	1.4828	1.4832
	$25{}^{°}$	NaN	1.3944	1.4332	1.5776	1.6720	1.6477	1.4676	1.4252	1.4474	1.4854	1.4892	1.4803	1.4776	1.4816
	$30{}^{°}$	NaN	NaN	1.3925	1.4369	1.5763	1.6771	1.6428	1.4573	1.4084	1.4733	1.4677	1.4715	1.4702	1.4730
	$35{}^{°}$	NaN	NaN	NaN	1.3925	1.4347	1.5815	1.6756	1.6354	1.4426	1.4446	1.4608	1.4641	1.4622	1.4682
	$40{}^{°}$	NaN	NaN	NaN	NaN	1.3905	1.4380	1.5786	1.6669	1.6189	1.4831	1.4273	1.4526	1.4554	1.4612
	$45{}^{°}$	NaN	NaN	NaN	NaN	NaN	1.3940	1.4356	1.5708	1.6451	1.6605	1.4663	1.4221	1.4444	1.4527
	$50{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	1.3910	1.4271	1.5550	1.6894	1.6367	1.4608	1.4121	1.4398
	$55{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.3824	1.4100	1.5894	1.6625	1.6276	1.4507	1.4099
	$60{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.3660	1.4435	1.5637	1.6466	1.6055	1.4458
	$65{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.4001	1.4173	1.5411	1.6072	1.5573
	$70{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.3734	1.3893	1.4656	1.5816
	$75{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.3180	1.3604	1.6773
	$80{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.5359	1.6471
	$85{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6787
Rate parameter, $b\times {10}^{-2}$
		Orbit inclination, i
		$20{}^{°}$	$25{}^{°}$	$30{}^{°}$	$35{}^{°}$	$40{}^{°}$	$45{}^{°}$	$50{}^{°}$	$55{}^{°}$	$60{}^{°}$	$65{}^{°}$	$70{}^{°}$	$75{}^{°}$	$80{}^{°}$	$85{}^{°}$
ES latitude, $\varphi$	$0{}^{°}$	9.3564	9.7586	9.5066	8.9569	8.8711	8.8567	8.8380	8.8193	8.8029	8.8260	8.8154	8.8086	8.7865	8.8199
	$5{}^{°}$	10.5403	10.0061	9.0549	8.9535	8.8763	8.8544	8.8469	8.8317	8.8354	8.8231	8.8158	8.8208	8.8155	8.8212
	$10{}^{°}$	7.7606	8.9977	9.8007	9.4661	8.8894	8.8741	8.8451	8.8116	8.7805	8.8244	8.7752	8.7949	8.7836	8.8062
	$15{}^{°}$	6.6219	7.7633	9.0428	9.8124	9.4342	8.8982	8.8559	8.8128	8.7869	8.8387	8.8202	8.8029	8.7816	8.7902
	$20{}^{°}$	6.6564	6.6032	7.7437	9.0479	9.7937	9.4531	8.8882	8.8237	8.7689	8.8436	8.8078	8.7884	8.7699	8.7666
	$25{}^{°}$	NaN	6.6326	6.5841	7.7402	9.0362	9.8353	9.4583	8.8648	8.8023	8.8671	8.8434	8.7874	8.7628	8.7717
	$30{}^{°}$	NaN	NaN	6.5839	6.5680	7.7144	9.0589	9.8236	9.4244	8.8083	8.8931	8.8074	8.7832	8.7553	8.7587
	$35{}^{°}$	NaN	NaN	NaN	6.5567	6.5410	7.7312	9.0548	9.7985	9.3706	8.9406	8.8534	8.7989	8.7558	8.7643
	$40{}^{°}$	NaN	NaN	NaN	NaN	6.5304	6.5513	7.7229	9.0267	9.7458	9.5196	8.8811	8.8282	8.7753	8.7717
	$45{}^{°}$	NaN	NaN	NaN	NaN	NaN	6.5434	6.5459	7.7016	8.9585	9.8790	9.4584	8.8695	8.8088	8.7893
	$50{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	6.5410	6.5295	7.6665	9.0970	9.7983	9.4411	8.8441	8.8240
	$55{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.5306	6.4975	7.7688	9.0234	9.7734	9.4087	8.8794
	$60{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.5095	6.6064	7.7225	8.9947	9.7121	9.4440
	$65{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.6420	6.5845	7.7079	8.9262	9.6598
	$70{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.6549	6.6089	7.6339	9.2649
	$75{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.7200	6.9875	8.5786
	$80{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	7.7001	7.3090
	$85{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.6598

,

Table A2. Shape, a, and rate, b parameter for $f_{\mathsf{\Theta }}\left(\theta \right)$ and $F_{\mathsf{\Theta }}\left(\theta \right)$ for $h_k=1200$ km.

Shape parameter, a
		Orbit inclination, i
		$20{}^{°}$	$25{}^{°}$	$30{}^{°}$	$35{}^{°}$	$40{}^{°}$	$45{}^{°}$	$50{}^{°}$	$55{}^{°}$	$60{}^{°}$	$65{}^{°}$	$70{}^{°}$	$75{}^{°}$	$80{}^{°}$	$85{}^{°}$
ES latitude, $\varphi$	$0{}^{°}$	1.9367	1.8476	1.7839	1.6868	1.6955	1.7417	1.7409	1.7492	1.7447	1.7597	1.7591	1.7601	1.7634	1.7602
	$5{}^{°}$	2.1177	1.8521	1.6327	1.7080	1.7098	1.7486	1.7488	1.7558	1.7608	1.7694	1.7714	1.7650	1.7716	1.7719
	$10{}^{°}$	1.6926	1.8959	1.9305	1.8002	1.6702	1.7141	1.7376	1.7431	1.7430	1.7589	1.7569	1.7578	1.7570	1.7563
	$15{}^{°}$	1.6207	1.7863	1.9154	1.9534	1.7777	1.6495	1.7212	1.7317	1.7364	1.7557	1.7562	1.7539	1.7587	1.7589
	$20{}^{°}$	1.6095	1.6490	1.8163	1.9346	1.9398	1.7980	1.6558	1.7163	1.7235	1.7514	1.7477	1.7507	1.7465	1.7494
	$25{}^{°}$	NaN	1.6232	1.6627	1.7997	1.9360	1.9544	1.7983	1.6558	1.7068	1.7490	1.7432	1.7479	1.7488	1.7452
	$30{}^{°}$	NaN	NaN	1.6252	1.6872	1.8215	1.9492	1.9548	1.7931	1.6375	1.7321	1.7337	1.7392	1.7343	1.7354
	$35{}^{°}$	NaN	NaN	NaN	1.6447	1.6618	1.8221	1.9537	1.9512	1.7807	1.6748	1.7177	1.7272	1.7325	1.7302
	$40{}^{°}$	NaN	NaN	NaN	NaN	1.6225	1.6601	1.8240	1.9460	1.9366	1.8135	1.6560	1.7160	1.7235	1.7234
	$45{}^{°}$	NaN	NaN	NaN	NaN	NaN	1.6254	1.6633	1.8164	1.9322	1.9702	1.7978	1.6528	1.7079	1.7157
	$50{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	1.6238	1.6543	1.8021	1.9626	1.9488	1.7952	1.6478	1.6998
	$55{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6163	1.6431	1.8301	1.9366	1.9394	1.7864	1.6401
	$60{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6051	1.6666	1.8028	1.9230	1.9178	1.7669
	$65{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6296	1.6402	1.7810	1.8767	1.8103
	$70{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6007	1.6094	1.6590	1.9132
	$75{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.5039	1.6502	1.9769
	$80{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.8388	1.9151
	$85{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.9243
Rate parameter, $b\times {10}^{-2}$
		Orbit inclination, i
		$20{}^{°}$	$25{}^{°}$	$30{}^{°}$	$35{}^{°}$	$40{}^{°}$	$45{}^{°}$	$50{}^{°}$	$55{}^{°}$	$60{}^{°}$	$65{}^{°}$	$70{}^{°}$	$75{}^{°}$	$80{}^{°}$	$85{}^{°}$
ES latitude, $\varphi$	$0{}^{°}$	9.8678	10.1172	10.1891	9.5261	9.2793	9.3729	9.3072	9.3066	9.2476	9.2919	9.2763	9.2667	9.2673	9.2552
	$5{}^{°}$	10.9490	10.7976	9.5396	9.4617	9.4702	9.4077	9.3339	9.3232	9.3145	9.3424	9.3382	9.2878	9.3186	9.3206
	$10{}^{°}$	7.9755	9.3073	10.1848	10.2797	9.5051	9.3438	9.3330	9.3051	9.2626	9.2973	9.2773	9.2654	9.2527	9.2481
	$15{}^{°}$	7.0189	8.0187	9.2511	10.2487	10.0220	9.3057	9.3599	9.3083	9.2744	9.3104	9.3002	9.2706	9.2890	9.2896
	$20{}^{°}$	7.1070	6.9667	8.0353	9.3656	10.2409	10.1571	9.3472	9.3401	9.2772	9.3222	9.2813	9.2749	9.2438	9.2536
	$25{}^{°}$	NaN	7.0470	6.9514	7.9165	9.3034	10.2353	10.1621	9.3559	9.2997	9.3585	9.2960	9.2845	9.2656	9.2505
	$30{}^{°}$	NaN	NaN	7.0018	6.9659	7.9953	9.3214	10.2443	10.1469	9.2910	9.3912	9.3140	9.2875	9.2399	9.2356
	$35{}^{°}$	NaN	NaN	NaN	7.0394	6.8984	8.0022	9.3392	10.2286	10.1082	9.4209	9.3499	9.2977	9.2695	9.2465
	$40{}^{°}$	NaN	NaN	NaN	NaN	6.9521	6.8918	7.9991	9.3174	10.1895	10.2072	9.3678	9.3487	9.2934	9.2639
	$45{}^{°}$	NaN	NaN	NaN	NaN	NaN	6.9573	6.9004	7.9858	9.2827	10.2757	10.1639	9.3584	9.3328	9.2961
	$50{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	6.9589	6.8851	7.9537	9.3612	10.2210	10.1526	9.3583	9.3438
	$55{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.9545	6.8728	8.0376	9.3045	10.1982	10.1379	9.3821
	$60{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.9538	6.9553	7.9965	9.2952	10.1634	10.1416
	$65{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	7.0605	6.9440	8.0014	9.2414	9.9957
	$70{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	7.0826	6.9882	7.8742	9.7894
	$75{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	7.1230	7.5669	8.9129
	$80{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	8.1780	7.6105
	$85{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.9258

,

Table A3. Shape, a, and rate, b parameter for $f_{\mathsf{\Theta }}\left(\theta \right)$ and $F_{\mathsf{\Theta }}\left(\theta \right)$ for $h_k=1400$ km.

Shape parameter, a
		Orbit inclination, i
		$20{}^{°}$	$25{}^{°}$	$30{}^{°}$	$35{}^{°}$	$40{}^{°}$	$45{}^{°}$	$50{}^{°}$	$55{}^{°}$	$60{}^{°}$	$65{}^{°}$	$70{}^{°}$	$75{}^{°}$	$80{}^{°}$	$85{}^{°}$
ES latitude, $\varphi$	$0{}^{°}$	2.0032	1.8456	1.8402	1.6716	1.7100	1.7340	1.7449	1.7511	1.7543	1.7637	1.7632	1.7651	1.7702	1.7679
	$5{}^{°}$	2.1655	1.9754	1.6522	1.6934	1.7231	1.7409	1.7490	1.7540	1.7591	1.7614	1.7649	1.7663	1.7667	1.7666
	$10{}^{°}$	1.6677	1.8992	1.9435	1.8657	1.6721	1.7127	1.7353	1.7436	1.7499	1.7599	1.7588	1.7621	1.7586	1.7636
	$15{}^{°}$	1.6088	1.7424	1.9101	1.9636	1.8729	1.6705	1.7134	1.7323	1.7393	1.7562	1.7543	1.7568	1.7675	1.7563
	$20{}^{°}$	1.6203	1.6962	1.7952	1.9295	1.9741	1.8764	1.6701	1.7079	1.7269	1.7506	1.7481	1.7527	1.7548	1.7517
	$25{}^{°}$	NaN	1.5976	1.6651	1.8075	1.9402	1.9817	1.8801	1.6663	1.7024	1.7427	1.7421	1.7498	1.7327	1.7498
	$30{}^{°}$	NaN	NaN	1.6369	1.6635	1.8108	1.9413	1.9814	1.8736	1.6573	1.7178	1.7295	1.7404	1.7402	1.7417
	$35{}^{°}$	NaN	NaN	NaN	1.6445	1.6715	1.8127	1.9434	1.9763	1.8665	1.6778	1.7050	1.7302	1.7423	1.7363
	$40{}^{°}$	NaN	NaN	NaN	NaN	1.6456	1.6718	1.8117	1.9378	1.9673	1.8869	1.6650	1.7061	1.7206	1.7275
	$45{}^{°}$	NaN	NaN	NaN	NaN	NaN	1.6449	1.6710	1.8075	1.9279	1.9885	1.8715	1.6685	1.7048	1.7163
	$50{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	1.6440	1.6645	1.7968	1.9466	1.9668	1.8700	1.6689	1.6933
	$55{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6368	1.6538	1.8118	1.9217	1.9582	1.8582	1.6591
	$60{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6265	1.6684	1.7868	1.9054	1.9326	1.8204
	$65{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6403	1.6414	1.7611	1.8477	1.8296
	$70{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6074	1.5970	1.6574	1.9844
	$75{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.5302	1.7400	2.0124
	$80{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.9051	1.9381
	$85{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.9298
Rate parameter, $b\times {10}^{-2}$
		Orbit inclination, i
		$20{}^{°}$	$25{}^{°}$	$30{}^{°}$	$35{}^{°}$	$40{}^{°}$	$45{}^{°}$	$50{}^{°}$	$55{}^{°}$	$60{}^{°}$	$65{}^{°}$	$70{}^{°}$	$75{}^{°}$	$80{}^{°}$	$85{}^{°}$
ES latitude, $\varphi$	$0{}^{°}$	9.4854	9.6587	9.9384	9.3116	9.0728	9.0317	8.9969	8.9756	8.9583	8.9661	8.9488	8.9440	8.9534	8.9475
	$5{}^{°}$	10.3295	10.6508	9.6070	9.1546	9.0625	9.0309	8.9989	8.9776	8.9675	8.9546	8.9515	8.9460	8.9355	8.9351
	$10{}^{°}$	7.6674	8.8435	9.6689	9.9236	9.2595	9.0497	9.0065	8.9747	8.9588	8.9634	8.9429	8.9396	8.9208	8.9386
	$15{}^{°}$	6.7932	7.5227	8.7334	9.6790	9.9223	9.2504	9.0394	8.9900	8.9562	8.9731	8.9424	8.9348	8.9586	8.9182
	$20{}^{°}$	6.9480	6.8614	7.6001	8.7492	9.6886	9.9382	9.2493	9.0196	8.9703	8.9873	8.9450	8.9408	8.9327	8.9139
	$25{}^{°}$	NaN	6.7756	6.7266	7.5953	8.7531	9.7130	9.9496	9.2345	9.0001	9.0167	8.9594	8.9542	8.8746	8.9231
	$30{}^{°}$	NaN	NaN	6.8307	6.6879	7.5795	8.7557	9.7114	9.9298	9.2092	9.0479	8.9779	8.9616	8.9292	8.9203
	$35{}^{°}$	NaN	NaN	NaN	6.8201	6.6873	7.5802	8.7575	9.6934	9.9105	9.2742	9.0104	8.9880	8.9710	8.9328
	$40{}^{°}$	NaN	NaN	NaN	NaN	6.8046	6.6837	7.5755	8.7420	9.6708	9.9675	9.2362	9.0237	8.9654	8.9487
	$45{}^{°}$	NaN	NaN	NaN	NaN	NaN	6.8013	6.6834	7.5681	8.7196	9.7269	9.9219	9.2539	9.0292	8.9823
	$50{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	6.8058	6.6760	7.5510	8.7727	9.6706	9.9225	9.2661	9.0351
	$55{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.8030	6.6671	7.6017	8.7233	9.6600	9.8991	9.2876
	$60{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.8055	6.7306	7.5743	8.7167	9.6255	9.8749
	$65{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.8898	6.7259	7.5850	8.6605	9.6210
	$70{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.9155	6.7657	7.6365	9.4222
	$75{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	7.1459	7.5283	8.4633
	$80{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	7.9448	7.2608
	$85{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.6159

,

Table A4. Shape, a, and rate, b parameter for $f_{\mathsf{\Theta }}\left(\theta \right)$ and $F_{\mathsf{\Theta }}\left(\theta \right)$ for $h_k=1600$ km.

Shape parameter, a
		Orbit inclination, i
		$20{}^{°}$	$25{}^{°}$	$30{}^{°}$	$35{}^{°}$	$40{}^{°}$	$45{}^{°}$	$50{}^{°}$	$55{}^{°}$	$60{}^{°}$	$65{}^{°}$	$70{}^{°}$	$75{}^{°}$	$80{}^{°}$	$85{}^{°}$
ES latitude, $\varphi$	$0{}^{°}$	2.0512	1.9235	1.8537	1.7406	1.6844	1.7333	1.7450	1.7537	1.7570	1.7618	1.7630	1.7656	1.7689	1.7680
	$5{}^{°}$	2.1942	2.0607	1.7717	1.6601	1.7164	1.7388	1.7508	1.7558	1.7613	1.7647	1.7677	1.7701	1.7681	1.7687
	$10{}^{°}$	1.7309	1.8164	1.9409	1.9105	1.7446	1.6908	1.7336	1.7464	1.7500	1.7590	1.7580	1.7620	1.7657	1.7720
	$15{}^{°}$	1.5853	1.7549	1.8968	1.9701	1.9155	1.7528	1.6897	1.7309	1.7402	1.7591	1.7563	1.7607	1.7604	1.7762
	$20{}^{°}$	1.6239	1.6532	1.7886	1.9184	1.9776	1.9304	1.7523	1.6878	1.7279	1.7522	1.7538	1.7591	1.7573	1.7647
	$25{}^{°}$	NaN	1.6489	1.6724	1.7993	1.9237	1.9910	1.9325	1.7512	1.6846	1.7395	1.7450	1.7540	1.7550	1.7582
	$30{}^{°}$	NaN	NaN	1.6547	1.6754	1.7984	1.9330	1.9935	1.9317	1.7473	1.6970	1.7291	1.7440	1.7478	1.7502
	$35{}^{°}$	NaN	NaN	NaN	1.6560	1.6715	1.8062	1.9326	1.9893	1.9235	1.7578	1.6838	1.7267	1.7376	1.7407
	$40{}^{°}$	NaN	NaN	NaN	NaN	1.6526	1.6789	1.8052	1.9302	1.9815	1.9365	1.7465	1.6836	1.7237	1.7222
	$45{}^{°}$	NaN	NaN	NaN	NaN	NaN	1.6606	1.6794	1.8026	1.9224	1.9943	1.9211	1.7483	1.6818	1.7154
	$50{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	1.6590	1.6753	1.7936	1.9316	1.9736	1.9156	1.7457	1.6728
	$55{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6550	1.6660	1.8010	1.9081	1.9610	1.9026	1.7332
	$60{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6450	1.6730	1.7752	1.8866	1.9276	1.8312
	$65{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6509	1.6425	1.7422	1.7973	1.9083
	$70{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6107	1.5687	1.7232	2.0355
	$75{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6165	1.8079	2.0373
	$80{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.9495	1.9564
	$85{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.9366
Rate parameter, $b\times {10}^{-2}$
		Orbit inclination, i
		$20{}^{°}$	$25{}^{°}$	$30{}^{°}$	$35{}^{°}$	$40{}^{°}$	$45{}^{°}$	$50{}^{°}$	$55{}^{°}$	$60{}^{°}$	$65{}^{°}$	$70{}^{°}$	$75{}^{°}$	$80{}^{°}$	$85{}^{°}$
ES latitude, $\varphi$	$0{}^{°}$	9.1342	9.5174	9.5944	9.3135	8.7818	8.7807	8.7426	8.7223	8.6959	8.6881	8.6741	8.6742	8.6742	8.6766
	$5{}^{°}$	9.7966	10.3752	9.7881	8.8842	8.8110	8.7746	8.7474	8.7192	8.7038	8.6957	8.6905	8.6867	8.6685	8.6762
	$10{}^{°}$	7.6274	8.2279	9.2022	9.6381	9.2692	8.7674	8.7567	8.7267	8.6895	8.6886	8.6647	8.6670	8.6692	8.6931
	$15{}^{°}$	6.6080	7.3062	8.2988	9.2230	9.6287	9.2905	8.7538	8.7376	8.6938	8.7173	8.6821	8.6788	8.6626	8.7096
	$20{}^{°}$	6.8179	6.5901	7.3007	8.3061	9.2175	9.6659	9.2847	8.7465	8.7247	8.7338	8.7024	8.6934	8.6663	8.6864
	$25{}^{°}$	NaN	6.7566	6.5695	7.2854	8.3000	9.2485	9.6723	9.2833	8.7361	8.7599	8.7128	8.7018	8.6785	8.6748
	$30{}^{°}$	NaN	NaN	6.7136	6.5383	7.2623	8.3150	9.2536	9.6715	9.2719	8.7749	8.7319	8.7107	8.6841	8.6771
	$35{}^{°}$	NaN	NaN	NaN	6.6823	6.5095	7.2720	8.3100	9.2416	9.6468	9.3038	8.7408	8.7284	8.6969	8.6821
	$40{}^{°}$	NaN	NaN	NaN	NaN	6.6578	6.5203	7.2661	8.3021	9.2209	9.6820	9.2726	8.7454	8.7268	8.6719
	$45{}^{°}$	NaN	NaN	NaN	NaN	NaN	6.6751	6.5234	7.2610	8.2885	9.2550	9.6393	9.2808	8.7496	8.7354
	$50{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	6.6781	6.5207	7.2526	8.3177	9.2060	9.6285	9.2832	8.7699
	$55{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.6827	6.5200	7.2868	8.2773	9.1899	9.6123	9.3061
	$60{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.6921	6.5679	7.2625	8.2650	9.1523	9.5221
	$65{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.7600	6.5619	7.2685	8.1612	9.4608
	$70{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.7813	6.5664	7.5990	9.0780
	$75{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	7.2701	7.4458	8.0980
	$80{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	7.7161	6.9900
	$85{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.3842

,

Table A5. Shape, a, and rate, b parameter for $f_{\mathsf{\Theta }}\left(\theta \right)$ and $F_{\mathsf{\Theta }}\left(\theta \right)$ for $h_k=1800$ km.

Shape parameter, a
		Orbit inclination, i
		$20{}^{°}$	$25{}^{°}$	$30{}^{°}$	$35{}^{°}$	$40{}^{°}$	$45{}^{°}$	$50{}^{°}$	$55{}^{°}$	$60{}^{°}$	$65{}^{°}$	$70{}^{°}$	$75{}^{°}$	$80{}^{°}$	$85{}^{°}$
ES latitude, $\varphi$	$0{}^{°}$	2.0864	1.9884	1.8169	1.7943	1.6412	1.7248	1.7438	1.7550	1.7609	1.7664	1.7714	1.7707	1.7703	1.7721
	$5{}^{°}$	2.2116	2.1213	1.8764	1.5897	1.7065	1.7368	1.7513	1.7599	1.7645	1.7665	1.7747	1.7742	1.7748	1.7745
	$10{}^{°}$	1.7923	1.7348	1.9270	1.9337	1.8142	1.6412	1.7269	1.7468	1.7539	1.7655	1.7619	1.7680	1.7676	1.7691
	$15{}^{°}$	1.5297	1.7356	1.8821	1.9664	1.9519	1.8180	1.6389	1.7278	1.7445	1.7576	1.7543	1.7642	1.7655	1.7646
	$20{}^{°}$	1.6239	1.6565	1.7779	1.9041	1.9797	1.9592	1.8191	1.6388	1.7249	1.7473	1.7510	1.7577	1.7601	1.7588
	$25{}^{°}$	NaN	1.6591	1.6767	1.7915	1.9167	1.9911	1.9649	1.8231	1.6361	1.7307	1.7462	1.7516	1.7553	1.7561
	$30{}^{°}$	NaN	NaN	1.6676	1.6810	1.7985	1.9204	1.9936	1.9658	1.8162	1.6417	1.7218	1.7392	1.7478	2.1573
	$35{}^{°}$	NaN	NaN	NaN	1.6709	1.6870	1.8007	1.9238	1.9959	1.9608	1.8265	1.6329	1.7194	1.7363	2.6343
	$40{}^{°}$	NaN	NaN	NaN	NaN	1.6724	1.6879	1.8011	1.9240	1.9903	1.9688	1.8155	1.6332	1.7164	1.7308
	$45{}^{°}$	NaN	NaN	NaN	NaN	NaN	1.6746	1.6876	1.8010	1.9156	1.9949	1.9541	1.8130	1.6334	1.7091
	$50{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	1.6731	1.6855	1.7926	1.9201	1.9743	1.9451	1.8060	1.6291
	$55{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6700	1.6776	1.7929	1.8943	1.9580	1.9234	1.7826
	$60{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6603	1.6760	1.7682	1.8697	1.9093	1.7970
	$65{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6585	1.6442	1.7219	1.7171	1.9728
	$70{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6097	1.5173	1.7822	2.0709
	$75{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6951	1.8594	2.0537
	$80{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.9798	1.9692
	$85{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.9458
Rate parameter, $b\times {10}^{-2}$
		Orbit inclination, i
		$20{}^{°}$	$25{}^{°}$	$30{}^{°}$	$35{}^{°}$	$40{}^{°}$	$45{}^{°}$	$50{}^{°}$	$55{}^{°}$	$60{}^{°}$	$65{}^{°}$	$70{}^{°}$	$75{}^{°}$	$80{}^{°}$	$85{}^{°}$
ES latitude, $\varphi$	$0{}^{°}$	8.8196	9.3308	9.1861	9.2502	8.5489	8.5670	8.5293	8.5090	8.4907	8.4811	8.4812	8.4660	8.4518	8.4567
	$5{}^{°}$	9.3468	10.0676	9.8322	8.6449	8.6084	8.5685	8.5380	8.5164	8.4966	8.4795	8.4907	8.4766	8.4677	8.4626
	$10{}^{°}$	7.5795	7.7597	8.7975	9.3397	9.2387	8.5089	8.5396	8.5159	8.4881	8.4885	8.4597	8.4644	8.4493	8.4529
	$15{}^{°}$	6.4075	7.0547	7.9508	8.8251	9.3577	9.2367	8.4880	8.5318	8.5000	8.4924	8.4577	8.4679	8.4561	8.4484
	$20{}^{°}$	6.7047	6.4674	7.0536	7.9492	8.8321	9.3703	9.2355	8.4831	8.5173	8.5034	8.4701	8.4654	8.4547	8.4423
	$25{}^{°}$	NaN	6.6531	6.4384	7.0389	7.9545	8.8556	9.3818	9.2454	8.4746	8.5307	8.5031	8.4748	8.4606	8.4516
	$30{}^{°}$	NaN	NaN	6.6111	6.4100	7.0327	7.9549	8.8569	9.3822	9.2265	8.4932	8.5107	8.4837	8.4704	11.6161
	$35{}^{°}$	NaN	NaN	NaN	6.5857	6.4036	7.0290	7.9578	8.8601	9.3696	9.2579	8.4709	8.5067	8.4846	17.2642
	$40{}^{°}$	NaN	NaN	NaN	NaN	6.5714	6.4000	7.0262	7.9577	8.8481	9.3920	9.2259	8.4777	8.5106	8.4929
	$45{}^{°}$	NaN	NaN	NaN	NaN	NaN	6.5756	6.3990	7.0300	7.9416	8.8626	9.3530	9.2231	8.4927	8.5222
	$50{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	6.5769	6.4042	7.0208	7.9641	8.8185	9.3373	9.2172	8.5267
	$55{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.5872	6.4044	7.0426	7.9226	8.7991	9.3054	9.2177
	$60{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.5948	6.4369	7.0282	7.9128	8.7370	9.1203
	$65{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.6509	6.4375	7.0166	7.7056	9.2740
	$70{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.6702	6.3810	7.5445	8.7644
	$75{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	7.3375	7.3430	7.7934
	$80{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	7.5042	6.7684
	$85{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.2083

and

Table A6. Shape, a, and rate, b parameter for $f_{\mathsf{\Theta }}\left(\theta \right)$ and $F_{\mathsf{\Theta }}\left(\theta \right)$ for $h_k=2000$ km.

Shape parameter, a
		Orbit inclination, i
		$20{}^{°}$	$25{}^{°}$	$30{}^{°}$	$35{}^{°}$	$40{}^{°}$	$45{}^{°}$	$50{}^{°}$	$55{}^{°}$	$60{}^{°}$	$65{}^{°}$	$70{}^{°}$	$75{}^{°}$	$80{}^{°}$	$85{}^{°}$
ES latitude, $\varphi$	$0{}^{°}$	2.1133	2.0407	1.8663	1.8277	1.6851	1.7151	1.7434	1.7543	1.7634	1.7690	1.7730	1.7757	1.7747	1.7801
	$5{}^{°}$	2.2236	2.1664	1.9608	1.6623	1.6910	1.7328	1.7511	1.7601	1.7679	1.7742	1.7770	1.7779	1.7768	1.7844
	$10{}^{°}$	1.8464	1.7634	1.9054	1.9463	1.8627	1.6881	1.7180	1.7460	1.7568	1.7711	1.7679	1.7719	1.7718	1.7758
	$15{}^{°}$	1.5715	1.7152	1.8664	1.9610	1.9725	1.8702	1.6873	1.7189	1.7443	1.7574	1.7629	1.7669	1.7658	1.7698
	$20{}^{°}$	1.6176	1.6582	1.7717	1.8939	1.9824	1.9815	1.8733	1.6862	1.7153	1.7485	1.7528	1.7588	1.7612	1.7578
	$25{}^{°}$	NaN	1.6675	1.6827	1.7860	1.9111	1.9917	1.9878	1.8782	1.6867	1.7171	1.7421	1.7518	1.7559	1.7612
	$30{}^{°}$	NaN	NaN	1.6781	1.6874	1.7972	1.9143	1.9936	1.9895	1.8745	1.6893	1.7104	1.7379	1.7453	1.7527
	$35{}^{°}$	NaN	NaN	NaN	1.6799	1.6963	1.7986	1.9170	1.9978	1.9867	1.8777	1.6819	1.7081	1.7317	1.7418
	$40{}^{°}$	NaN	NaN	NaN	NaN	1.6883	1.6955	1.7993	1.9211	1.9909	1.9881	1.8674	1.6807	1.7453	1.7237
	$45{}^{°}$	NaN	NaN	NaN	NaN	NaN	1.6861	1.6955	1.7996	1.9118	1.9915	1.9746	1.8627	1.6781	1.6952
	$50{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	1.6849	1.6942	1.7925	1.9092	1.9715	1.9623	1.8499	1.6713
	$55{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6842	1.6870	1.7886	1.8857	1.9512	1.9312	1.8113
	$60{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6740	1.6828	1.7609	1.8523	1.8820	1.8489
	$65{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6644	1.6463	1.6979	1.7474	2.0230
	$70{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.6026	1.5620	1.8330	2.0953
	$75{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.7623	1.8999	2.0647
	$80{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	2.0023	1.9790
	$85{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	1.9584
Rate parameter, $b\times {10}^{-2}$
		Orbit inclination, i
		$20{}^{°}$	$25{}^{°}$	$30{}^{°}$	$35{}^{°}$	$40{}^{°}$	$45{}^{°}$	$50{}^{°}$	$55{}^{°}$	$60{}^{°}$	$65{}^{°}$	$70{}^{°}$	$75{}^{°}$	$80{}^{°}$	$85{}^{°}$
ES latitude, $\varphi$	$0{}^{°}$	8.5432	9.1399	9.0974	9.1202	8.5964	8.3881	8.3606	8.3318	8.3159	8.3014	8.2988	8.2944	8.2826	8.2926
	$5{}^{°}$	8.9698	9.7735	9.7793	8.8345	8.4312	8.3942	8.3640	8.3383	8.3238	8.3219	8.3099	8.3004	8.2846	8.3080
	$10{}^{°}$	7.5197	7.6751	8.4453	9.0607	9.1342	8.5590	8.3611	8.3452	8.3173	8.3285	8.2929	8.2901	8.2799	8.2844
	$15{}^{°}$	6.4748	6.8508	7.6660	8.4985	9.0872	9.1384	8.5429	8.3518	8.3260	8.3161	8.2966	8.2894	8.2711	8.2777
	$20{}^{°}$	6.6023	6.3675	6.8675	7.6707	8.5195	9.1022	9.1422	8.5359	8.3333	8.3363	8.2948	8.2840	8.2721	8.2523
	$25{}^{°}$	NaN	6.5665	6.3426	6.8486	7.6805	8.5340	9.1141	9.1546	8.5354	8.3373	8.3136	8.2950	8.2791	8.2811
	$30{}^{°}$	NaN	NaN	6.5273	6.3099	6.8477	7.6785	8.5321	9.1163	9.1432	8.5477	8.3180	8.3065	8.2834	8.2823
	$35{}^{°}$	NaN	NaN	NaN	6.4935	6.3074	6.8422	7.6793	8.5404	9.1072	9.1546	8.5241	8.3179	8.2998	8.2902
	$40{}^{°}$	NaN	NaN	NaN	NaN	6.4955	6.3003	6.8407	7.6895	8.5230	9.1128	9.1241	8.5273	8.2834	8.2992
	$45{}^{°}$	NaN	NaN	NaN	NaN	NaN	6.4891	6.3005	6.8438	7.6709	8.5293	9.0800	9.1167	8.5360	8.3301
	$50{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	6.4918	6.3065	6.8397	7.6762	8.4912	9.0590	9.1003	8.5634
	$55{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.5068	6.3106	6.8529	7.6475	8.4680	9.0132	9.0711
	$60{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.5164	6.3391	6.8386	7.6242	8.3726	9.0403
	$65{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.5571	6.3402	6.8044	7.6273	9.0766
	$70{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.5688	6.4544	7.4789	8.4834
	$75{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	7.3598	7.2366	7.5378
	$80{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	7.3162	6.5865
	$85{}^{°}$	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	6.0765

in Appendix B, show the shape and scale parameters obtained from simulations with configurations shown in Table 1, Table 2 and Table 3. From those values, matrices ${\mathit{A}}_{\mathit{h}}$ and ${\mathit{B}}_{\mathit{h}}$ can be constructed, and linear interpolation can be applied to those in order to obtain shape and scale parameters for different orbit inclinations, ES’s latitudes, and altitudes that were not simulated, but are in the range of $h_s$, $i_s$, ${\varphi }_s$.

3.4. Reducing the Shape and Rate Matrices

From the grid of Figure 9c, an array as shown in Figure 9d can be created for each simulated orbit altitude. Figure 9d shows the diagonal numbers of a rectangular matrix ${\mathit{D}}_{\mathit{h}}$ defined as

$${\mathit{D}}_{\mathit{h}}={\left[\begin{array}{ccccc}{d}_0& d_1& d_2& \dots & d_{13}\\ d_{-1}& d_0& d_1& \ddots & d_{12}\\ d_{-2}& d_{-1}& d_0& \ddots & ⋮\\ d_{-3}& d_{-2}& d_{-1}& \ddots & ⋮\\ d_{-4}& d_{-3}& d_{-2}& \ddots & ⋮\\ NaN& d_{-4}& d_{-3}& \ddots & ⋮\\ NaN& NaN& d_{-4}& \ddots & ⋮\\ ⋮& ⋮& ⋮& \ddots & ⋮\\ NaN& NaN& NaN& \dots & d_{-4}\end{array}\right]}_{18\times 14},$$

(19)

for each altitude $h\in h_s$, with $NaN$’s are placed for cases where ${\varphi }_{ES}>i_{SAT}$.

3.5. PDF and CDF of the Elevation Angle

In addition to $h_s$, we can define two additional sets containing the orbit inclinations (as in Table 2) and ES’s latitudes (as in Table 3) at which simulations were performed, as 

$$i_s=\{{20}^{°},{25}^{°},\dots ,{85}^{°}\}$$

(20)

and

$${\varphi }_s=\{0^{°},5^{°},\dots ,{85}^{°}\},$$

(21)

respectively.

Depending on the input values of ${\varphi }_k$, $h_k$, and $i_k$, if those are part of the previously defined sets ($i_s$, $h_s$, ${\varphi }_s$) eight cases can arise as follows:

• Case 1. $i_k\in i_s\mathrm{and}{\varphi }_k\in {\varphi }_s\mathrm{and}{h}_k\in h_s$
• Case 2. $i_k\in i_s\mathrm{and}{\varphi }_k\in {\varphi }_s\mathrm{and}{h}_k\notin h_s$
• Case 3. $i_k\in i_s\mathrm{and}{\varphi }_k\notin {\varphi }_s\mathrm{and}{h}_k\in h_s$
• Case 4. $i_k\in i_s\mathrm{and}{\varphi }_k\notin {\varphi }_s\mathrm{and}{h}_k\notin h_s$
• Case 5. $i_k\notin i_s\mathrm{and}{\varphi }_k\in {\varphi }_s\mathrm{and}{h}_k\in h_s$
• Case 6. $i_k\notin i_s\mathrm{and}{\varphi }_k\in {\varphi }_s\mathrm{and}{h}_k\notin h_s$
• Case 7. $i_k\notin i_s\mathrm{and}{\varphi }_k\notin {\varphi }_s\mathrm{and}{h}_k\in h_s$
• Case 8. $i_k\notin i_s\mathrm{and}{\varphi }_k\notin {\varphi }_s\mathrm{and}{h}_k\notin h_s$

Case 1 occurs when all the input values ${\varphi }_k$, $i_k$ and $h_k$ are contained in ${\varphi }_s$, $i_s$, and $h_s$, respectively; and Case 8 occurs when none of ${\varphi }_k$, $i_k$ and $h_k$ parameters coincide with previously simulated values. Cases 2 to 7 occur when at least one of the input values (any of ${\varphi }_k$, $i_k$, and $h_k$) coincides with a value within its corresponding sets ${\varphi }_s$, $i_s$ and $h_s$. Those cases illustrate all possible scenarios, ranging from not-required interpolation to interpolation inside a four-point mesh.

We can represent the ES latitude, ${\varphi }_k$, orbit inclination, $i_k$, or both (for which we want to determine their PDF) as a query point between two-known points, inside a three-point mesh, or inside a four-point mesh. Then we can use an interpolation method, such as those proposed in Appendix C, to find an appropriate value for the gamma distribution parameters at that query coordinates. In addition to those proposed in Appendix C, more interpolation methods are widely available in the literature.

3.6. First and Second Order Statistics of $f_{\mathsf{\Theta }}\left(\theta \right)$

The elevation angle expected value can be calculated from (15) as

$$E\left[\mathsf{\Theta }\right]={\int }_{\theta }\theta f_{\mathsf{\Theta }}\left(\theta \right)\mathrm{d}\theta ,$$

(22)

Nonetheless, sometimes a satellite link requires to operate above a minimum value of $\theta$, ${\theta }_{min}$, then, elevation angle values below ${\theta }_{min}$ are not of interest. Recent satellite systems operating above the Ku band consider ${\theta }_{min}$ values above some tens of degrees, e.g., the Starlink constellation considered a ${\theta }_{min}={40}^{°}$ as mentioned in [27]. We can obtain the conditional expected value for values above ${\theta }_{min}$ using conditional probability as follows

$$E\left[\mathsf{\Theta }\right|\mathsf{\Theta }\ge {\theta }_{min}]={\int }_{\theta }\theta f_{\mathsf{\Theta }}\left(\theta \right|\theta \ge {\theta }_{min})\mathrm{d}\theta ,$$

(23)

which can be rewritten as in [28] as

$$E\left[\mathsf{\Theta }\right|\mathsf{\Theta }\ge {\theta }_{min}]={\displaystyle \frac{{\int }_{{\theta }_{min}}^{{\theta }_{max}}\theta f_{\mathsf{\Theta }}\left(\theta \right)\mathrm{d}\theta }{F_{\mathsf{\Theta }}\left({\theta }_{max}\right)-F_{\mathsf{\Theta }}\left({\theta }_{min}\right)}},$$

(24)

The variance of $\mathsf{\Theta }$, is obtained as follows

$$\mathrm{Var}\left[\mathsf{\Theta }\right|\mathsf{\Theta }\ge {\theta }_{min}]=E\left[{\mathsf{\Theta }}^2\right|\mathsf{\Theta }\ge {\theta }_{min}]-E{\left[\mathsf{\Theta }\right|\mathsf{\Theta }\ge {\theta }_{min}]}^2,$$

(25)

And the standard deviation of $\mathsf{\Theta }$ is then defined as

$$\begin{array}{c}\hfill \mathrm{SD}\left[\mathsf{\Theta }\right|\mathsf{\Theta }\ge {\theta }_{min}]=\mathrm{Var}{\left[\mathsf{\Theta }\right|\mathsf{\Theta }\ge {\theta }_{min}]}^{1/2}\end{array}$$

(26)

3.7. Choosing Orbits to Maximize Mean Value of $\theta$

It is possible to choose an orbit configuration to maximize the elevation angle expected value, $E\left[\mathsf{\Theta }\right|\mathsf{\Theta }\ge {\theta }_{min}]$, for given ranges of $\varphi$, h, and i. The problem can be stated as

$$\begin{array}{c}\hfill \underset{i,\varphi ,h}{max}E\left[\mathsf{\Theta }\right|\mathsf{\Theta }\ge {\theta }_{min}]\\ \hfill \mathrm{s}.\mathrm{t}.& i_1\le i\le i_2\hfill \\ \hfill & {\varphi }_1\le \varphi \le {\varphi }_2\hfill \\ \hfill & h_1\le h\le h_2\hfill \end{array}$$

(27)

where $i_1\le i_2$, ${\varphi }_1\le {\varphi }_2$, and $h_1\le h_2$. This problem can be solved by optimization methods or using computational tools to perform an iterative evaluation within a for or while cycle.

Reduction of the Elements in the Diagonals

It was found that functions $f_{\mathsf{\Theta }}\left(\theta \right)$ and $F_{\mathsf{\Theta }}\left(\theta \right)$ were very similar for some diagonals, d shown in Figure 9d. Then, from the diagonals as shown in Figure 9d, we propose a reduction first, by obtaining the expected value $E[·]$ of each diagonal, $d_l$, (note the use of subindex $l\in \{-4,-3,\dots ,13\}$ as an indicator of the diagonal number) as $E\left[d_l\right]$. Then, (19) can be rewritten as the following row vector

$$E\left[{\mathit{D}}_{\mathit{h}}\right]={[E\left[d_{-4}\right],E\left[d_{-3}\right],\dots ,E\left[d_{13}\right]]}_{1\times 18}$$

(28)

where the diagonals with just $NaN$’s elements are omitted.

From observation, we find that some adjacent elements in $E\left[{\mathit{D}}_{\mathit{h}}\right]$ are approximately equal, then, we reduced $E\left[{\mathit{D}}_{\mathit{h}}\right]$ by taking the mean value of adjacent similar values according to

$${\displaystyle \frac{1}{|L-l|}}\sum _l^{L}E\left[d_l\right],L>l$$

(29)

where l, and L are, respectively, the lower and upper diagonal numbers for the interval in which the values $E\left[d_l\right],E\left[d_{l+1}\right],\dots ,E\left[d_L\right]$ were observed approximately the same (with a maximum-distance criteria between the fitted curves of 5%). After this reduction $E\left[{\mathit{D}}_{\mathit{h}}\right]$ ends as a row vector with only four elements, and we call this vector ${\mathit{p}}_{\mathit{h}}$ for $h\in h_s$. This procedure can be repeated for each altitude in $h_s$, and after grouping all the row vectors ${\mathit{p}}_{\mathit{h}}$ for all simulated altitudes h we arrive to the matrix $\mathit{P}$ given by

$$\mathit{P}=\left[\begin{array}{c}{\mathit{p}}_{\mathbf{1000}\mathrm{km}}\\ {\mathit{p}}_{\mathbf{1200}\mathrm{km}}\\ {\mathit{p}}_{\mathbf{1400}\mathrm{km}}\\ {\mathit{p}}_{\mathbf{1600}\mathrm{km}}\\ {\mathit{p}}_{\mathbf{1800}\mathrm{km}}\\ {\mathit{p}}_{\mathbf{2000}\mathrm{km}}\end{array}\right]=\left[\begin{array}{cccc}{p}_{1,1}& p_{1,2}& p_{1,3}& p_{1,4}\\ p_{2,1}& p_{2,2}& p_{2,3}& p_{2,4}\\ p_{3,1}& p_{3,2}& p_{3,3}& p_{3,4}\\ p_{4,1}& p_{4,2}& p_{4,3}& p_{4,4}\\ p_{5,1}& p_{5,2}& p_{5,3}& p_{5,4}\\ p_{6,1}& p_{6,2}& p_{6,3}& p_{6,4}\end{array}\right]$$

(30)

Matrix $\mathit{P}$ will be called ${\mathit{P}}_{\mathit{a}}$ if it contains the shape parameters, a, for $f_{\mathsf{\Theta }}\left(\theta \right)$, and ${\mathit{P}}_{\mathit{b}}$ if it contains the scale parameter, b, for $f_{\mathsf{\Theta }}\left(\theta \right)$.

From matrices ${\mathit{P}}_{\mathit{a}}$ and ${\mathit{P}}_{\mathit{b}}$, it is possible to recover any of the points of ${\mathit{A}}_{\mathit{h}}$ and ${\mathit{B}}_{\mathit{h}}$ and then estimate the CDF for a query orbit altitude $h_k$, orbit inclination, $i_k$, and ES’s latitude, ${\varphi }_k$. However, we observed that some $F_{\mathsf{\Theta }}\left(\theta \right)$ curves for high latitudes were not fitting very well after the reduction from matrices ${\mathit{A}}_{\mathit{h}}$ and ${\mathit{B}}_{\mathit{h}}$ to ${\mathit{P}}_{\mathit{a}}$ and ${\mathit{P}}_{\mathit{b}}$, then, two weight matrices were obtained through manual tuning to improve the correspondence between ${\mathit{P}}_{\mathit{a}}$ and ${\mathit{A}}_{\mathit{h}}$, and between ${\mathit{P}}_{\mathit{b}}$ and ${\mathit{B}}_{\mathit{h}}$. The weight matrix for the shape parameter, ${\mathit{w}}_{\mathit{a}}$ and ${\mathit{w}}_{\mathit{b}}$, is defined as

$${\mathit{w}}_{\mathit{a}}={\left[\begin{array}{ccccc}{w}_{a_{1,1}}& w_{a_{1,2}}& w_{a_{1,3}}& \dots & w_{a_{1,14}}\\ w_{a_{2,1}}& w_{a_{2,2}}& w_{a_{2,3}}& \ddots & w_{a_{2,14}}\\ w_{a_{3,1}}& w_{a_{3,2}}& w_{a_{3,3}}& \ddots & ⋮\\ w_{a_{4,1}}& w_{a_{4,2}}& w_{a_{4,3}}& \ddots & ⋮\\ w_{a_{5,1}}& w_{a_{5,2}}& w_{a_{5,3}}& \ddots & ⋮\\ NaN& w_{a_{6,2}}& w_{a_{6,3}}& \ddots & ⋮\\ NaN& NaN& w_{a_{7,3}}& \ddots & ⋮\\ ⋮& ⋮& ⋮& \ddots & ⋮\\ NaN& NaN& NaN& \dots & w_{a_{18,14}}\end{array}\right]}_{18\times 14},{\mathit{w}}_{\mathit{b}}={\left[\begin{array}{ccccc}{w}_{b_{1,1}}& w_{b_{1,2}}& w_{b_{1,3}}& \dots & w_{b_{1,14}}\\ w_{b_{2,1}}& w_{b_{2,2}}& w_{b_{2,3}}& \ddots & w_{b_{2,14}}\\ w_{b_{3,1}}& w_{b_{3,2}}& w_{b_{3,3}}& \ddots & ⋮\\ w_{b_{4,1}}& w_{b_{4,2}}& w_{b_{4,3}}& \ddots & ⋮\\ w_{b_{5,1}}& w_{b_{5,2}}& w_{b_{5,3}}& \ddots & ⋮\\ NaN& w_{b_{6,2}}& w_{b_{6,3}}& \ddots & ⋮\\ NaN& NaN& w_{b_{7,3}}& \ddots & ⋮\\ ⋮& ⋮& ⋮& \ddots & ⋮\\ NaN& NaN& NaN& \dots & w_{b_{18,14}}\end{array}\right]}_{18\times 14},$$

(31)

Both ${\mathit{w}}_{\mathit{a}}$ and ${\mathit{w}}_{\mathit{b}}$ are the same for all altitudes, and their numerical values are shown in Appendix B.

4. Results

The shape, a, and rate, b, parameters for matrices ${\mathit{P}}_{\mathit{a}}$ and ${\mathit{P}}_{\mathit{b}}$, respectively, are shown in

Table 6. Shape parameter, a, for matrix ${\mathit{P}}_{\mathit{a}}$, using the diagonal notation as in (19), and as shown in Figure 9d.

	Diagonals Range
Altitude	(${\mathit{d}}_{-4},{\mathit{d}}_{-2}$)	(${\mathit{d}}_{-1},{\mathit{d}}_0$)	(${\mathit{d}}_1,{\mathit{d}}_2$)	(${\mathit{d}}_3,{\mathit{d}}_{13}$)
1000 km	1.470954	1.638425	1.450794	1.488802
1200 km	1.707323	1.935321	1.721572	1.752709
1400 km	1.712923	1.951520	1.746205	1.754760
1600 km	1.723436	1.955254	1.782190	1.757023
1800 km	1.731110	1.952362	1.796549	1.774827
2000 km	1.741450	1.954794	1.836833	1.760709

and

Table 7. Rate parameter, $b\times {10}^{-2}$, for matrix ${\mathit{P}}_{\mathit{b}}$, using the diagonal notation as in (19), and as shown in Figure 9d.

	Diagonals Range
Altitude	(${\mathit{d}}_{-4},{\mathit{d}}_{-2}$)	(${\mathit{d}}_{-1},{\mathit{d}}_0$)	(${\mathit{d}}_1,{\mathit{d}}_2$)	(${\mathit{d}}_3,{\mathit{d}}_{13}$)
1000 km	7.037052	9.459459	9.075805	8.805617
1200 km	7.391963	9.840010	9.640436	9.284945
1400 km	7.130587	9.319436	9.412412	8.952740
1600 km	6.943998	8.879456	9.242146	8.693598
1800 km	6.791048	8.510588	9.041604	8.632867
2000 km	6.675182	8.224503	8.941424	8.305444

for the orbit altitudes $h\in h_s$ and for the diagonals as in ${\mathit{D}}_{\mathit{h}}$, (19).

Figure 10a–c show the PDF obtained from the empirical data and the ones obtained from the parameters listed in Table 6 and Table 7.

Table 8. First and second order statistics of the elevation angle for the same orbit configurations and ES’s locations that in Figure 10a–d.

	Empirical Result	Proposed Method’s Result	Empirical Result	Proposed Method’s Results	Empirical Result	Proposed Method’s Results	Empirical Result	Proposed Method’s Results
ES latitude, ${\varphi }_{ES}$	30${}^{°}$		45${}^{°}$		20${}^{°}$		85${}^{°}$
Orbit inclination, $i_{SAT}$	30${}^{°}$		65${}^{°}$		75${}^{°}$		85${}^{°}$
Orbit altitude, $h_{SAT}$	1000 km
Shape parameter, a	–	1.4710	–	1.6384	–	1.4888	–	1.8387
Rate Parameter, $b\times {10}^{-2}$	–	7.0371	–	9.4595	–	8.8056	–	7.0371
$E[\mathsf{\Theta }>5{}^{°}]$	24.55${}^{°}$	23.09${}^{°}$	19.43${}^{°}$	19.74${}^{°}$	20.24${}^{°}$	19.70${}^{°}$	27.70${}^{°}$	26.66${}^{°}$
$\mathrm{SD}[\mathsf{\Theta }>5{}^{°}]$	18.40${}^{°}$	16.86${}^{°}$	14.19${}^{°}$	13.13${}^{°}$	14.94${}^{°}$	13.48${}^{°}$	18.12${}^{°}$	18.90${}^{°}$
$E[\mathsf{\Theta }>15{}^{°}]$	34.77${}^{°}$	31.13${}^{°}$	29.77${}^{°}$	27.96${}^{°}$	30.55${}^{°}$	28.24${}^{°}$	35.48${}^{°}$	33.23${}^{°}$
$\mathrm{SD}[\mathsf{\Theta }>15{}^{°}]$	17.46${}^{°}$	16.26${}^{°}$	14.24${}^{°}$	12.45${}^{°}$	14.49${}^{°}$	12.93${}^{°}$	16.32${}^{°}$	17.96${}^{°}$
$E[\mathsf{\Theta }>25{}^{°}]$	44.33${}^{°}$	39.66${}^{°}$	40.12${}^{°}$	37.03${}^{°}$	40.39${}^{°}$	37.38${}^{°}$	42.73${}^{°}$	40.91${}^{°}$
$\mathrm{SD}[\mathsf{\Theta }>25{}^{°}]$	15.93${}^{°}$	15.88${}^{°}$	13.40${}^{°}$	12.05${}^{°}$	13.50${}^{°}$	12.61${}^{°}$	14.79${}^{°}$	17.29${}^{°}$

shows the shape, a, and scale, b, parameters as well as the first and second order statistics for the $f_{\mathsf{\Theta }}\left(\theta \right)$ functions of Figure 10a–d.

A case study was developed for a satellite with an altitude of 1500 km, an orbit inclination of 43${}^{°}$, and ES latitude of 22${}^{°}$. Using the numerical results of Table 6 and Table 7, a mesh as shown in Figure 11a can be created for altitudes of 1400 km and 1600 km (shape and scale values were obtained for those altitudes and are shown in Table 6 and Table 7). Then, using an interpolation method as that described in Appendix C, parameters a and b can be obtained to characterize the elevation angle curve.

The shape and scale parameters for the case study are shown in

Table 9. Shape and scale parameters for the elevation angle PDF and CDF for ${\varphi }_{ES}={22}^{°}$ and $i_{SAT}={43}^{°}$ at $h_k=1500$ km, corresponding to the four-point mesh shown in Figure 11a.

Four-Point Mesh Setup
	1400	1600		1400	1600
$a({\varphi }_1,i_1)$	1.9515	1.9552	$b({\varphi }_1,i_1)$	9.3194 $\times {10}^{-2}$	8.8794 $\times {10}^{-2}$
$a({\varphi }_1,i_2)$	1.9515	1.9552	$b({\varphi }_1,i_2)$	9.3194 $\times {10}^{-2}$	8.8794 $\times {10}^{-2}$
$a({\varphi }_2,i_1)$	1.9515	1.9552	$b({\varphi }_2,i_1)$	9.3194 $\times {10}^{-2}$	8.8794 $\times {10}^{-2}$
$a({\varphi }_2,i_2)$	1.7462	1.7822	$b({\varphi }_2,i_2)$	9.4112 $\times {10}^{-2}$	9.2215 $\times {10}^{-2}$
Obtained parameters from interpolation
$a(\varphi ,i)$	1.8776	1.8929	$b(\varphi ,i)$	9.3529 $\times {10}^{-2}$	9.0025 $\times {10}^{-2}$
$a(\varphi ,i)$	1.8853		$b(\varphi ,i)$	9.1777 $\times {10}^{-2}$
Empirical parameters
$a(\varphi ,i)$	1.9804		$b(\varphi ,i)$	9.5987 $\times {10}^{-2}$

. First, parameters for the four-points (as in Figure 11b) are listed for each altitude (1400 km and 1600 km). Then, the resultant parameters obtained from interpolation were computed. Finally, the shape and scale parameters for the query altitude can be obtained using a weighted arithmetic mean were the weights can be obtained from the distance of the query altitude to the closest characterized orbit altitudes. In this case, the distances from the query altitude to the closest characterized altitudes are the same, then we obtained the shape and scale for the query altitude using the arithmetic mean of the shape and scale values at 1400 km and 1600 km.

Both the elevation angle CDF obtained with the proposed methodology and the empirical CDF are shown in Figure 11b. The distance magnitude between those two curves is also shown in Figure 11b and labeled as absolute error ($\left|\mathrm{Error}\right|$). The distance (or error) between the gamma CDF and the empirical CDF is very small, and its maximum value is around 0.025 (2.5%).

4.1. Results Validity

4.1.1. Effects of the Orbit Configuration in the Elevation Angle Distribution

The Kepler orbital elements were chosen to describe the orbits, those parameters include the semi-major axis, a, eccentricity, e, orbit inclination, i, right ascension of the ascending node, $\mathsf{\Omega }$, argument of perigee, $\omega$, and true anomaly, $\nu$. The values of the right ascension of the ascending node, $\mathsf{\Omega }$, argument of perigee, $\omega$, and true anomaly, $\nu$, do not affect the statistical properties of the elevation angle, then, those values are not further discussed. Figure 12a,b show the PDF and CDF elevation angle curves observed from the same location for two simulated satellites with different orbit configurations which share the same initial semi-major axis, a, eccentricity, e, and orbit inclination, i; but, differ in the initial values of the right ascension of the ascending node, $\mathsf{\Omega }$, argument of perigee, $\omega$, and of the true anomaly, $\nu$.

4.1.2. Elevation Angle CDF for Different Longitudes

ES’s located at different longitudes share approximately the same elevation angle distribution for the same LEO satellite configuration. The reasoning behind this fact is discussed in [2], and it is illustrated in Figure 13a, which shows the elevation angle PDF for three ES’s located at the same latitude, $\varphi$, but at different longitudes, $\lambda$’s. Figure 13a shows that the elevation angle PDF is the same for ES located at different longitudes, $\lambda$’s, when they share the same latitude, $\varphi$, and those are served by satellites with the same orbit configuration (semi-major axis, a, eccentricity, e, and, inclination, i). The PDF’s shown in Figure 13a were generated for a satellite with a circular orbit and an inclination of $i={60}^{°}$, and an altitude of $h=1800$ km; in addition, all ES’s were located at the same latitude $\varphi ={25}^{°}$ N.

4.1.3. Long Term Validity of $f_{\mathsf{\Theta }}\left(\theta \right)$ and $F_{\mathsf{\Theta }}\left(\theta \right)$

From a few days $f_{\mathsf{\Theta }}\left(\theta \right)$ can be observed very similar compared with the same PDF is a much greater period, e.g., one year. Figure 13b shows the functions $F_{\mathsf{\Theta }}\left(\theta \right)$ for three different simulation periods ranging from one week to one year. The CDF curve is almost identical for all the simulation periods. The orbit configuration and ES’s location were choose to be the same that for Figure 13a.

4.1.4. Elevation Angle CDF for Different Satellite Characteristics

Simulations for this paper where performed for a satellite with the mass and drag areas listed in Table 1, nonetheless, satellites with different values in those characteristics will share very similar functions $F_{\mathsf{\Theta }}\left(\theta \right)$. Figure 14 shows the empirical CDF’s for two satellites with the same orbital parameters but with different mass and drag areas. For both satellites, the orbit altitude, h, was set to 1800 km and the orbit inclination, i, to 60${}^{°}$, and both ES’s were located at 25${}^{°}$ N. As shown in Figure 14, variations in the satellite mass and drag area characteristics will have a very little impact in the elevation angle CDF’s. For lower altitudes the effect is greater, but not significant enough for the altitudes covered in this paper.

Table 10. Satellites with different mass and drag area characteristics for additional performed simulations.

	Satellite 1	Satellite 2
Dry mass	5 kg	200 kg
Drag area	1 m${}^2$	10 m${}^2$
Solar radiation pressure area	1 m${}^2$	10 m${}^2$

show the mass and areas for satellites shown in Figure 14.

4.2. Error Analysis

The error between the empirical CDF curves of the elevation angle and the ones obtained with the parameters in Table 6 and Table 7, was quantified by means of the absolute error, $ϵ$, and its CDF. Figure A1a–f in Appendix A show the mean absolute error, $\overline{ϵ}$, between the empirical CDF’s for the simulated points of Figure 9c, and for the orbit altitudes $h_k\in h_s$.

The maximum observed mean error for the individual characterizations shown in Appendix A was about 5% (${\overline{ϵ}}_{max}\approx 0.05$) for the worst case. However, the greatest error values were observed to occur at lower elevation angles (below 20${}^{°}$ were the link is often unreliable since line-of-sight is harder to find) that are often not considered by satellite communication systems due to greater attenuation.

The absolute errors $ϵ$ were also analyzed through their individual CDF’s, showing that most of the errors where below 5%. Figure 15a–d show the CDF’s of the absolute error, $ϵ$, corresponding to Figure 10a–d. In addition, $ϵ$ was analyzed for all the simulated satellite orbits and according to their elevation angle range.

The PDF, CDF, and the complementary CDF (CCDF) showing the resultant absolute errors, $ϵ$, for all the simulations are shown in Figure 16, where it can be observed that around 95% of the values of $ϵ$ are below 0.05. Figure 17 shows the same absolute errors as in Figure 16b, but now expanded and classified by elevation angle range; in this figure, it can be observed that most of the absolute errors are below 5% for all the elevation angle ranges, and that errors’ CDF distributions are approximately the same for all the elevation angle ranges. Figure 18 expands the 90th-percentile of the CDF’s in Figure 17, showing that in fact, most of the errors are below 2%; the 90th-percentile mean value is shown as a red cross for each boxplot.

5. Conclusions

This article demonstrated the feasibility of using a gamma random variable to characterize the elevation angle PDF and CDF. The proposed methodology allowed PDF and CDF calculations of the elevation angle for LEO satellites for altitudes between 1000 km and 2000 km. The characterization of elevation angle through a gamma random variable allowed an easy computation of first and second order statistics, which, to the best of our knowledge, have not been previously addressed in the literature.

The proposed methodology was validated with an error analysis against the empirical CDF’s. The results showed that errors are low, with a mean absolute error much below 5% in most of the cases, and with around 95% of the total errors being below 5%.

Furthermore, the proposed methodology allows an easy comparison between multiple orbits in order to determine the most suitable orbital parameters to provide coverage above a minimum elevation angle at a particular latitude.