This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

\receiveddate

XX November,2023 \reviseddateXX Month, 2023 \accepteddateXX Month, 2023 \publisheddateXX Month, 2023 \currentdate2 November,2022 \doiinfoOJVT.2023.1234567

\authornote

This work was supported in part by the High Throughput and Secure Networks (HTSN) Challenge Program of the National Research Council of Canada (NRC), and in part by MDA.

Free-Space Optical (FSO) Satellite Networks Performance Analysis: Transmission Power, Latency, and Outage Probability

JINTAO LIANG (Member    IEEE)1    AIZAZ U. CHAUDHRY (Senior Member    IEEE)1    EYLEM ERDOGAN (Senior Member    IEEE)2    HALIM YANIKOMEROGLU (Fellow    IEEE)1    GUNES KARABULUT KURT (Senior Member    IEEE)3    PENG HU (Senior Member    IEEE)4    KHALED AHMED5    STEPHANE MARTEL5 Department of Systems and Computer Engineering, Carleton University, Ottawa, ON K1S 5B6, Canada Department of Electrical and Electronics Engineering, Istanbul Medeniyet University, Istanbul, 34700, Turkey Poly-Grames Research Center, Department of Electrical Engineering, Polytechnique Montreal, Montreal, QC H3T 1J4, Canada National Research Council of Canada (NRC), Ottawa, ON K1A 0R6, Canada Technology Strategy of Satellite Systems, MDA, Sainte-Anne-de-Bellevue QC H9X 3Y5, Canada
Abstract

In free-space optical satellite networks (FSOSNs), satellites can have different laser inter-satellite link (LISL) ranges for connectivity. As the LISL range increases, the number of satellites from among all the satellites in the constellation that will be needed on the shortest path between a source and a destination ground station decrease, and thereby the number of the LISLs on the shortest path decreases. Greater LISL ranges can reduce network latency of the path but can also result in an increase in transmission power for satellites on the path. Consequently, this tradeoff between satellite transmission power and network latency should be investigated, and in this work we examine it in FSOSNs drawing on the Starlink Phase 1 Version 3 (i.e., the latest version of Starlink’s Phase 1) and Kuiper Shell 2 (i.e., Kuiper’s biggest shell) constellations for different LISL ranges and different inter-continental connections. We use appropriate system models for calculating the average satellite transmission power (i.e., the average of the transmission power of all satellites on the shortest path) and network latency (i.e., the end-to-end latency of the shortest path). The results show that the mean network latency (i.e., the mean of network latency over all time slots) decreases and mean average satellite transmission power (i.e., the mean of average satellite transmission power over all time slots) increases with an increase in LISL range. For the Toronto–Sydney inter-continental connection in an FSOSN with Starlink’s Phase 1 Version 3 constellation, when the LISL range is approximately 2,900 km, the mean network latency and mean average satellite transmission power intersect are approximately 135 ms and 380 mW, respectively. For an FSOSN with the Kuiper Shell 2 constellation in this inter-continental connection, this LISL range is around 3,800 km, and the two parameters are approximately 120 ms and 700 mW, respectively. For the Toronto–Istanbul and Toronto–London inter-continental connections, the LISL ranges at the intersection are different and vary from 2,600 km to 3,400 km. Furthermore, we analyze outage probability performance of optical uplink/downlink due to atmosphere attenuation and turbulence.

{IEEEkeywords}

Free-space optical satellite networks, network latency, optical inter-satellite link, optical uplink/downlink, satellite transmission power, tradeoff, outage probability.

1 Introduction

Free-space optical wireless communications are being developed rapidly and have a promising future. In particular, free-space optical satellite networks (FSOSNs) [1] will be realized through laser inter-satellite links (LISLs) [2], also known as optical inter-satellite links, in upcoming low Earth orbit (LEO) satellite constellations, such as SpaceX’s Starlink [3],[4], Amazon’s Kuiper [5], and Telesat’s Lightspeed [6]. These new FSOSNs will have various advantages compared to satellite networks using radio frequency based inter-satellite links (ISLs), including larger link bandwidth, license-free spectrum, higher link data rate, and lower terminal power required to establish ISLs. On the other hand, FSOSNs have some downsides, such as attenuation of the optical signal during transmission through Earth’s atmosphere [7] and a limitation on the amount of energy their satellite solar panels can obtain [8]. As the mission lifetime of a satellite is mainly determined by its battery life, maximizing a satellite’s battery life is critical for minimizing the replacement cost of the satellites. Since battery life is determined by the number of charging and recharging cycles, if the satellite consumes less power, the battery recharge time will decrease, and its lifetime can increase. The analysis of the power consumption by the communication system of the satellite with respect to the total power budget of the satellite is out of the scope of this paper. However, satellite transmission power can be an important factor that can affect battery life and, by extension, the lifetime of the satellite [9]. Therefore, in this work, we focus on the analysis of satellite transmission power for optical links.

For LISLs in FSOSNs, the optical beam from a transmitting satellite can be received without attenuation. For uplink (i.e., the link between a source ground station and an ingress satellite in the FSOSN) and downlink (i.e., the link between an egress satellite in the FSOSN and a destination ground station), we consider optical links instead of radio frequency links, since optical links provide larger link bandwidth, higher link data rate, better security, smaller aperture size, and lower terminal power consumption compared to radio frequency links [1],[10]. For optical uplink and downlink, the optical beams must propagate through Earth’s atmosphere to reach the receiver, and multiple attenuation affect the beams and thereby attenuate the transmitted optical signal. Thus, satellite transmission power for optical uplink/downlink must consider the losses caused by the atmosphere, such as those due to Mie scattering and geometrical scattering [11]. In addition, it’s essential to consider atmospheric turbulence, as it has the potential to induce signal fading and consequently lead to a deterioration in the quality of the communication, including outage probability.

Latency is an important parameter to gauge the performance of data communications over satellite networks in various applications, such as the high-frequency trading of stocks between financial stock exchanges around the world [12], information exchanges for military purposes [13], etc. In such communications networks, there are four components of network latency, namely propagation delay, transmission delay, queuing delay, and processing delay. Propagation delay only depends on the propagation distance between the source and destination over the satellite network, while the rest are affected by the link data rate, the switching and routing speed, and the traffic status [14]. The processing delay, transmission delay, and queuing delay constitute the node delay.

In FSOSNs, satellites have an LISL range for connectivity, which is the range within which a satellite can successfully establish an LISL with another satellite. For longer LISL ranges, optical links between satellites on the shortest path over an FSOSN have longer propagation distances, and higher satellite transmission power is required to maintain reliable links between satellites. On the other hand, fewer satellites are needed to establish the shortest path for an inter-continental connection between two ground stations over an FSOSN, resulting in less node delay and thereby less network latency of the shortest path. In this way, different LISL ranges lead to a tradeoff between satellite transmission power and network latency, which is the focus of this paper.

To study this tradeoff, we employ the Starlink Phase 1 Version 3 [4] and Kuiper Shell 2 [5] constellations, and we consider three different inter-continental connections, specifically Toronto to Sydney, Toronto to Istanbul, and Toronto to London, and different LISL ranges for satellites. The results show that as the LISL ranges increases, the mean of network latency (i.e., the mean of network latency of the shortest paths at all time slots for an inter-continental connection) decreases. At the same time, the mean of average satellite transmission power (i.e., the mean of average satellite transmission power of the shortest paths at all time slots for an inter-continental connection) increases, where the average satellite transmission power is the mean of the transmission power of all satellites on a shortest path at a time slot.

In addition to the observation through figures, the results further show that for the Toronto–Sydney inter-continental connection over the Starlink Phase 1 Version 3 constellation, the mean of network latency and mean average satellite transmission power intersect when the LISL range is approximately 2,900 km and are approximately 135 ms and 380 mW, respectively. For the Kuiper Shell 2 constellation in this inter-continental connection, this LISL range is around 3,800 km, and these two parameters are approximately 120 ms and 700 mW, respectively. For the Toronto–Istanbul inter-continental connection, the two parameters intersect when LISL ranges are about 2,600 km for Starlink Phase 1 Version 3 constellation and 2,900 km for Kuiper Shell 2 constellation. For the Toronto to London inter-continental connection, these parameters intersect when LISL ranges are 3,400 km and 3,000 km for the Starlink Phase 1 Version 3 and Kuiper Shell 2, respectively. Based on the results, we find that for LISL and laser uplink/downlink, the network latency for the shortest path for an inter-continental connection between two ground stations and the average satellite transmission power for the shortest path have an inverse relationship. For different LISL ranges for the two FSOSNs over different inter-continental connections, we find that there is no clear relationship between the LISL range when network latency and satellite transmission power reach balance and the length of the inter-continental connection within the same FSOSN. Furthermore, we find that the performance of the outage probability is inevitably influenced by cloud conditions.

The rest of the paper is organized as follows. Related work and motivation are discussed in Section II. Section III presents the system model, including models for laser link transmission power, models for turbulence induced fading channel, and models for link and node latency. Section IV introduces the two constellations and provides methodology to calculate LISL and ground station ranges, orbital period, outage probability, network latency, and average satellite transmission power. Section V discusses the results for the tradeoff between average satellite transmission power and network latency, the effect of atmosphere turbulence on outage probability, and mentions some insights based on these results. Conclusions and directions for future work are provided in Section VI.

2 Related Work and Motivation

Several studies have discussed the relationship between satellite power consumption and battery life and have analyzed latency for satellite networks, such as [8], [9], [15, 16, 17, 18, 19, 20, 21, 22, 23, 24] and the references therein.

In [8], the authors discussed satellite battery lifetime and indicated that typical satellite missions in a LEO system are projected to last for 3 to 5 years, which can be roughly regarded as a range of 16,500 to 27,500 partial cycles. In order to maintain such a large number of cycles, battery lifetime and approaches to reduce the satellite power consumption should be considered.

In [9], the authors conducted a study on power management in Low Earth Orbit (LEO) satellite constellations, emphasizing the critical role of the electrical power system in ensuring the Quality of Service (QoS) for LEO satellite communications. They developed a model for satellite battery characteristics based on various communication traffic scenarios and proposed several approaches to extend the battery lifetime.

In [15], the authors highlighted the constraints of limited satellite power and stressed the need for resource allocation mechanisms to take satellite power consumption during transmission into account, particularly for downlink communication. Due to the increasing importance of on-board transmission power optimization, they introduced an innovative carrier and power allocation method for satellite systems to reduce satellite power consumption.

In [16], the authors introduced a communication methodology designed to efficiently manage the communication resources of LEO satellites in order to extend the battery lifetime. They proposed a system design aimed at extending the battery lifetime on board the satellites. This system can enable the deployment of large-scale LEO satellite constellations while satisfying the communication requirements of current satellite communications.

In [17], the authors focused on the power consumption required for satellite terminal communications as well as the battery Depth of Discharge (DOD) during the communication period. They also discussed power control systems, and the results of their study showed that with the reduction of the transmission power, less DOD is required per cycle and therefore longer battery life can be achieved and, by extension, longer satellite lifetimes.

In [18], the authors simulated Starlink Phase 1 Version 1 constellation and established their study based on the fundamental characteristics of the Starlink deployment. They conducted simulations of routing designs within this satellite network. They analyzed the constellation parameters and offered several preliminary insights into the dynamics of the constellation’s evolving topology. They also discussed the possible routing strategies and the resultant end-to-end latency properties of this constellation.

In [19], the authors framed the problem of inter-satellite topology design for large LEO constellation, including Starlink and Kuiper. They proposed a new approach for designing satellite networks using regular repetitive patterns and investigated the topology and impact of ISL range and setup time of satellite network based on this approach. Furthermore, they evaluated and compared the network performance of the two satellite constellations.

In [20], the authors studied the link reliability of the software defined satellite network (SDSN). They designed a failure probability model for ISLs that considered queuing and processing delays and the corresponding SDSN reliability analysis was added. The authors simulated the satellite network with certain constellation and latency parameters, including queuing delay and processing delay, and analyzed the network performance.

In [21], the author investigated the latency performance of large LEO constellations. The author studied the use of ground-based relays for the LEO constellations to provide low-latency wide area networking. For upcoming mega-constellations like Starlink, the author mentioned two inter-continental connections and the results indicated that implementing ISLs could reduce latency and supplement capacity in heavy traffic in network.

In [22], the authors investigated the time delay and outage performance for inter-satellite transmissions, satellite-terrestrial transmissions, terrestrial relay assisted transmissions, and satellite relay assisted transmissions for LEO mage-constellations. They proposed analytical models to investigate the end-to-end time delay and outage performance while considering different transmission scenarios. They presented numerical results to validate the proposed analytical models and compare the performances in terms of time delay and outage probability. They further concluded that the size of coverage area of LEO satellites exhibited a weak negative influence on time delay.

In [23], the authors focused on the reliability performance of downlink transmission in LEO satellite communication. They derived a novel expression for the outage probability of LEO satellite transmission. They simulated the LEO satellite communication scenarios with severe Doppler effects and analyzed the performance of the system based on orthogonal time frequency space (OTFS) scheme and traditional orthogonal frequency division multiplexing (OFDM) scheme.

The tradeoff between average satellite transmission power and network latency in FSOSNs arises from different LISL ranges. Figs. 1 and 2 show different shortest paths for different LISL ranges for the same inter-continental connection from Toronto to Istanbul. In these figures, uplinks and downlinks are highlighted in green, while LISLs are shown in yellow. In Fig. 1, satellites have a longer LISL range and there are only three satellites on the shortest path. In this way, each LISL is longer, and each satellite needs higher transmission power to an LISL. However, only three node delays are added to the network latency in addition to propagation delays of the four links. In Fig. 2, satellites have a shorter LISL range, and so there are more satellites on the shortest path, which leads to higher network latency but shorter LISLs and therefore lower satellite transmission power.

In an earlier work in [7], we studied system models for transmission power and link margin for the LISL and laser uplink/downlink to examine the impact of link distance and link margin on LISL transmission power, and the impact of slant distance, elevation angle, and link margin on laser uplink/downlink transmission power. In another earlier work in [24], we proposed a binary integer linear program formulation of the total network latency over multiple LISL ranges and multi-pairs of inter-continental connection and simulated the Starlink Phase 1 Version 3 constellation. Furthermore, we considered several satellite transmission power constraints to investigate the tradeoff. In this work, we use these models from [7] and [24] in calculating the average satellite transmission power for satellites on the shortest path between a source and destination ground station over an FSOSN for studying the tradeoff between average satellite transmission power and network latency in FSOSNs. Investigations on satellite transmission power and latency have been conducted separately in the literature, but a tradeoff between satellite transmission power and network latency and analysis of the outage probability has never been investigated before. We consider both transmission delay and node delay while analyzing the tradeoff, and we obtain the results based on simulation of the on-going Starlink Phase 1 Version 3 and Kuiper Shell 2 constellation. To the best of our knowledge, this is the first study, which focuses on the tradeoff between average satellite transmission power and network latency for FSOSNs with the Starlink Phase 1 Version 3 and Kuiper Shell 2 constellations, and investigate the outage probability due to atmosphere turbulence. This is the first work of its kind, and no similar previous works exist.

Refer to caption

Figure 1: An illustration of a shortest path for a longer LISL range for an inter-continental connection from Toronto to Istanbul.

Refer to caption

Figure 2: An illustration of a shortest path for a shorter LISL range for an inter-continental connection from Toronto to Istanbul.

3 System Model

Here, we introduce the system model for FSOSNs, which includes models for LISL transmission power, laser uplink/downlink transmission power, link margin, turbulence induced fading channel, link latency, and node latency.

3.1 Transmission Power in FSOSNs

In the following, we introduce the laser link transmission power model and the link margin model.

3.1.1 Laser inter-satellite link transmission power model

For the inter-satellite link, the transmission power PT\textit{P}_{T} is given by [25]

PT=PR/(GTGRLTLRLPSηTηR),\textit{P}_{T}=\textit{P}_{R}/(\textit{G}_{T}\textit{G}_{R}\textit{L}_{T}\textit{L}_{R}\textit{L}_{PS}\textit{$\eta$}_{T}\textit{$\eta$}_{R})\text{,} (1)

where PT\textit{P}_{T} is the transmitted power for LISL in Watts, PR\textit{P}_{R} is the received power for the link in Watts, GT\textit{G}_{T} is the transmitter gain, GR\textit{G}_{R} is the receiver gain, LT\textit{L}_{T} is the transmitter pointing loss, LR\textit{L}_{R} is the receiver pointing loss, LPS\textit{L}_{PS} is the free-space path loss for LISL, ηT\textit{$\eta$}_{T} is the transmitter optics efficiency, and ηR\textit{$\eta$}_{R} is the receiver optics efficiency. The transmitter gain GT\textit{G}_{T} in (1) is given by [26]

GT=16/(ΘT)2,\textit{G}_{T}=16/({\textit{$\Theta$}_{T}})^{2}\text{,} (2)

where ΘT\textit{$\Theta$}_{T} is the full transmitting divergence angle in radians. The receiver gain GR\textit{G}_{R} in (1) is given by [25]

GR=(πDR/λ)2,\textit{G}_{R}=(\textit{$\pi$}\textit{D}_{R}/\textit{$\lambda$})^{2}\text{,} (3)

where DR\textit{D}_{R} is the receiver’s telescope diameter in mm, and λ\lambda is the operating wavelength in nm. The transmitter pointing loss LT\textit{L}_{T} in (1) is given by [25]

LT=exp(GT(θT)2),\textit{L}_{T}=\exp{(-\textit{G}_{T}(\textit{$\theta$}_{T})^{2})}\text{,} (4)

where θT\textit{$\theta$}_{T} is the transmitter pointing error in radians. The receiver pointing loss LR\textit{L}_{R} in (1) is given by [25]

LR=exp(GR(θR)2),\textit{L}_{R}=\exp{(-\textit{G}_{R}(\textit{$\theta$}_{R})^{2})}\text{,} (5)

where θR\textit{$\theta$}_{R} is the receiver pointing error in radians. The free-space path loss for LISL LPS\textit{L}_{PS} is given by [25]

LPS=(λ/4πdSS)2,\textit{L}_{PS}=(\textit{$\lambda$}/4\textit{$\pi$}\textit{d}_{SS})^{2}\text{,} (6)

where λ\lambda is the operating wavelength in nm, and dSS\textit{d}_{SS} is the propagation distance between satellites in km.

3.1.2 Laser uplink/downlink transmission power model

Uplink and downlink are the links between a ground station and satellite. In this work, we consider Mie scattering and geometrical scattering in terms of atmospheric attenuation for optical uplink/downlink. Mie scattering can be caused by the particles in the medium and it can redirect the signal from its intended direction, while geometrical scattering causes reflection, refraction, and scattering of the optical signal.

Atmospheric attenuation must be considered to model laser uplink and downlink, and the transmission power for laser uplink/downlink PT\textit{P}_{T} is given by [27]

PT=PR/(GTGRLTLRLALPGηTηR),\textit{P}_{T}=\textit{P}_{R}/(\textit{G}_{T}\textit{G}_{R}\textit{L}_{T}\textit{L}_{R}\textit{L}_{A}\textit{L}_{PG}\textit{$\eta$}_{T}\textit{$\eta$}_{R})\text{,} (7)

where LA\textit{L}_{A} is the atmospheric attenuation loss, LPG\textit{L}_{PG} is the free-space path loss for laser links between ground stations and satellites, and the rest of the parameters are same as the LISL transmission power model in (1). The slant distance (i.e., the distance between a ground station and a satellite) for laser uplink/downlink dGS\textit{d}_{GS} is given by [32]

dGS=R(((R+H)/R)2(cos(θE))2sin(θE)),\textit{d}_{GS}=\textit{R}(\sqrt{((\textit{R}+\textit{H})/\textit{R})^{2}-(\cos({\textit{$\theta$}_{E}}))^{2}}-\sin({\textit{$\theta$}_{E}}))\text{,} (8)

where R == r ++ hE\textit{h}_{E}, H == hS\textit{h}_{S} - hE\textit{h}_{E}, r is the radius of the Earth, hE\textit{h}_{E} is the height of the ground station, hS\textit{h}_{S} is the altitude of the satellite, and θE\textit{$\theta$}_{E} is the elevation angle of the ground station in degrees. The free-space path loss LPG\textit{L}_{PG} can be expressed on the basis of slant distance accordingly:

LPG=(λ/4πdGS)2.\textit{L}_{PG}=(\textit{$\lambda$}/4\textit{$\pi$}\textit{d}_{GS})^{2}\text{.} (9)

a) Atmospheric attenuation due to Mie scattering
    When the diameter of atmosphere particles equals or exceeds the wavelength of the optical beam, it results in Mie scattering. This occurs because the phase of the optical beam is non-uniform across the particle and leads to the deviation of the transmitted signal from its intended direction. It is mainly caused by microscopic particles of water, and it happens in the lower part of the atmosphere where these particles are more abundant. Mie scattering can be modeled by the following expression [33]:

ρ=a(hE)3+b(hE)2+chE+d,\textit{$\rho$}=\textit{a}(\textit{h}_{E})^{3}+\textit{b}(\textit{h}_{E})^{2}+\textit{c}\textit{h}_{E}+\textit{d}\text{,} (10)

where ρ\rho is the extinction ratio, hE\textit{h}_{E} is the height of the ground station in km, and a, b, c and d are empirical coefficients dependent on wavelength, which can be expressed as follows [33]:

a=0.000545λ2+0.002λ0.0038,\textit{a}=-0.000545\textit{$\lambda$}^{2}+0.002\textit{$\lambda$}-0.0038\text{,} (11)
b=0.00628λ20.0232λ+0.00439,\textit{b}=0.00628\textit{$\lambda$}^{2}-0.0232\textit{$\lambda$}+0.00439\text{,} (12)
c=0.028λ2+0.101λ0.18,\textit{c}=-0.028\textit{$\lambda$}^{2}+0.101\textit{$\lambda$}-0.18\text{,} (13)
d=0.228λ3+0.922λ21.26λ+0.719.\textit{d}=-0.228\textit{$\lambda$}^{3}+0.922\textit{$\lambda$}^{2}-1.26\textit{$\lambda$}+0.719\text{.} (14)

The atmospheric attenuation due to Mie scattering can be calculated as

Im=exp(ρ/sin(θE)),\textit{I}_{m}=\exp{(-\textit{$\rho$}/\sin({\textit{$\theta$}_{E}}))}\text{,} (15)

where θE\textit{$\theta$}_{E} is the elevation angle of the ground station in degrees.

b) Atmospheric attenuation due to geometrical scattering
    Geometrical scattering is used to model the attenuation that takes place near the surface of the Earth and is produced by fog, dense clouds, water molecules, rain or other atmospheric phenomena. Geometrical scattering can be modeled by the following expression [11]:

V=1.002/(NLW)0.6473,\textit{V}=1.002/(\textit{N}\textit{L}_{W})^{0.6473}\text{,} (16)

where V is the visibility in km, N is the cloud number concentration in cm-3, and LW\textit{L}_{W} is the liquid water content in g/m-3. The attenuation coefficient θA\textit{$\theta$}_{A} can be expressed as [11]

θA=(3.91/V)(λ/550)φ,\textit{$\theta$}_{A}=(3.91/\textit{V})(\textit{$\lambda$}/550)^{-\textit{$\varphi$}}\text{,} (17)

where φ\varphi is the particle size related coefficient based on Kim’s model. The Beer-Lambert law is given by [34]

I(z)=exp(μz),\textit{I}(\textit{z})=\exp{(-\textit{$\mu$}\textit{z})}\text{,} (18)

where μ\mu is the wavelength dependent attenuation coefficient, and z is the distance of the transmission path. The atmospheric attenuation caused by geometrical scattering can be expressed using the Beer-Lambert law as

Ig=exp(θAdA),\textit{I}_{g}=\exp{(-\textit{$\theta$}_{A}\textit{d}_{A})}\text{,} (19)

where dA\textit{d}_{A} is the distance of the laser beam through the troposphere and can be calculated as [7]

dA=(hAhE)csc(θE),\textit{d}_{A}=(\textit{h}_{A}-\textit{h}_{E})\csc({\textit{$\theta$}_{E}})\text{,} (20)

where hE\textit{h}_{E} is the height of the ground station in km, hA\textit{h}_{A} is the height of the troposphere in km, and θE\textit{$\theta$}_{E} is the elevation angle of the ground station in degrees.
    In this work, we consider geometrical scattering for uplink, and Mie scattering and geometrical scattering for downlink. The atmospheric attenuation loss for uplink considering only geometrical scattering can be expressed as

LAup=Ig=exp(θAdA),\textit{L}_{Aup}=\textit{I}_{g}=\exp{(-\textit{$\theta$}_{A}\textit{d}_{A})}\text{,} (21)

and the atmospheric attenuation loss for downlink considering both Mie scattering and geometrical scattering can then be expressed as

LAdown=ImIg=exp(ρ/sin(θE))exp(θAdA).\textit{L}_{Adown}=\textit{I}_{m}\textit{I}_{g}=\exp{(-\textit{$\rho$}/\sin({\textit{$\theta$}_{E}}))}\exp{(-\textit{$\theta$}_{A}\textit{d}_{A})}\text{.} (22)

3.1.3 Link Margin Model

The link margin is defined as the ratio of the received signal power PR\textit{P}_{R} to the required signal power Preq\textit{P}_{req} to achieve a specific bit error rate (BER) at a given data rate, and it can be represented as [27]

LM=PR/Preq.\textit{LM}=\textit{P}_{R}/\textit{P}_{req}\text{.} (23)

Since we are interested in computing the link transmission power in this work, we need to know the received power based on link margin. The required signal power is generally referred to as receiver sensitivity, and the received power can be calculated as

PR=LM×Preq.\textit{P}_{R}=\textit{LM}\texttimes\textit{P}_{req}\text{.} (24)

3.2 Atmosphere turbulence in FSOSNs

In the following, we introduce the an atmosphere channel model which is under the influence of atmosphere turbulence resulting in intensity fluctuations.

In this work, we consider that the turbulence induced channel experiences an exponentiated Weibull fading channel.The probability density function (PDF) is given by [28]

f(I)=(αβ/η)(I/η)β1exp((I/η)β)(1exp((I/η)β))α1,\textit{f}(\textit{I})=(\textit{$\alpha$}\textit{$\beta$}/\textit{$\eta$})(\textit{I}/\textit{$\eta$})^{\textit{$\beta$}-1}\exp{(-(\textit{I}/\textit{$\eta$})^{\textit{$\beta$}})}\\ (1-\exp{(-(\textit{I}/\textit{$\eta$})^{\textit{$\beta$}})})^{\textit{$\alpha$}-1}\text{,} (25)

and the cumulative distribution function (CDF) is given by [28]

F(I)=(1exp(I/η)β)α,\textit{F}(\textit{I})=(1-\exp{-(\textit{I}/\textit{$\eta$})^{\textit{$\beta$}}})^{\textit{$\alpha$}}\text{,} (26)

where α\alpha, β\beta are the shape parameters and η\eta is the scale parameter of a ground station. The three parameters can be expressed as [29]

α=7.220α2/3/Γ(2.487σI1/30.104),\textit{$\alpha$}=7.220\textit{$\alpha$}^{2/3}/\Gamma{(2.487\textit{$\sigma$}_{I}^{1/3}-0.104)}\text{,} (27)
β=1.012(ασI2)0.52+0.142,\textit{$\beta$}=1.012(\textit{$\alpha$}\textit{$\sigma$}_{I}^{2})^{-0.52}+0.142\text{,} (28)
η=1/(αΓ(1+1/β)g(α,β)),\textit{$\eta$}=1/(\textit{$\alpha$}\Gamma{(1+1/\textit{$\beta$})}\textit{g}(\textit{$\alpha$},\textit{$\beta$}))\text{,} (29)

where g(α,β)\textit{g}(\textit{$\alpha$},\textit{$\beta$}) is the α\alpha and β\beta dependent constant variable, and σI\textit{$\sigma$}_{I} is the scintillation index of the ground station, which is expressed as [30]

σI=exp(0.49σR2/(1+1.11σR2.4)7/6+0.51σR/(1+0.69σR2.4)5/6)1,\sigma_{I}=\exp\left(0.49\sigma_{R}^{2}/(1+1.11\sigma_{R}^{2.4})^{7/6}\right.\\ +\left.0.51\sigma_{R}/(1+0.69\sigma_{R}^{2.4})^{5/6}\right)-1,

where the Rytov variance σR\textit{$\sigma$}_{R} is given by [30]

σR=2.25k7/6sec11/6(ζ)hEhSCn2(h)(hhE)5/6𝑑h,\textit{$\sigma$}_{R}=2.25\textit{k}^{7/6}\sec^{11/6}{(\textit{$\zeta$})}\int_{\textit{h}_{E}}^{\textit{h}_{S}}\textit{C}_{n}^{2}(\textit{h})(\textit{h}-\textit{h}_{E})^{5/6}\,dh\text{,} (30)

where k is the optical wave number, ζ\zeta is the zenith angle of the ground station, hS\textit{h}_{S} is the altitude of the satellite, and Cn2\textit{C}_{n}^{2} is the altitude dependent refractive index constant, which is expressed as [31]

Cn2(h)=8.148×1056vr2h10eh/1000+2.7×1016eh/1500+C0eh/100m2/3,\textit{C}_{n}^{2}(\textit{h})=8.148\times 10^{-56}\textit{v}_{r}^{2}\textit{h}^{10}\textit{e}^{-\textit{h}/1000}\\ +2.7\times 10^{-16}\textit{e}^{-\textit{h}/1500}+\textit{C}_{0}\textit{e}^{-\textit{h}/100}\textit{m}^{-2/3}\text{,} (31)

where vr\textit{v}_{r} is the r.m.s. ground wind speed, and C0\textit{C}_{0} is the nominal value of the refractive index constant at ground.

3.3 Latency in FSOSNs

In the following, we introduce the modeling of latency in FSOSNs. In optical communication networks, four types of delay can be considered as part of the end-to-end latency. These delays consist of propagation delay, transmission delay, queueing delay, and processing delay. Propagation delay of an optical link depends on its propagation distance, which we also refer to as link latency. Node delay, also referred to as node latency in this work, comprises transmission, queueing, and processing delays at a satellite.

3.3.1 Link Latency Model

The propagation delay of a link (or link latency) is the delay the transmitted signal takes to be received by the receiver through the channel. It depends on the propagation medium and distance. In FSOSNs, the propagation delay of a laser link Lp\textit{L}_{p} or link latency Tlink\textit{T}_{link} can be calculated as

Tlink=Lp=dlink/cs,\textit{T}_{link}=\textit{L}_{p}=\textit{d}_{link}/\textit{c}_{s}\text{,} (32)

where dlink\textit{d}_{link} is the length of the laser link (or propagation distance between transmitter and receiver satellites), and cs\textit{c}_{s} is the speed of light, which is 299,792,458 m/s [35].

3.3.2 Node Latency Model

According to [1], satellite required a precise acquisition, tracking, and pointing (ATP) system for effective and correct connection with other satellites due to the narrow beam divergence of the laser beam and the different satellite motion velocities. Based on the parameters of the using optical terminals and routers, current setup delays for satellites to complete ATP process to establish temporary LISLs for nodes switching are usually as huge as a few seconds and prohibitive for LISLs. Therefore, in this work we consider the node delay without setup delay.
    Transmission delay refers to the delay it takes to transmit a packet on the channel on the transmitter side. It depends on the size of the data packet and the data rate of the link. In FSOSNs, transmission delay Lt\textit{L}_{t} can be calculated as

Lt=l/Rdata,\textit{L}_{t}=\textit{l}/\textit{R}_{data}\text{,} (33)

where l is the length of a packet, and Rdata\textit{R}_{data} is the data rate of the laser link. Transmission delay can be minimal if the link data rate is very high. For example, if the link data rate is 10 Gbps and the packet size is 1,500 bytes, Lt\textit{L}_{t} is 1.2 μ\mus, which is negligible.
    Queuing delay Lq\textit{L}_{q} refers to the delay that a packet incurs while waiting in the buffer of a router. If a new packet arrives while an old one is still waiting or being transmitted, it will stay in the buffer and wait. For an M/M/1 queuing system [36], where the first M denotes the packet arrival rate, the second M denotes the packet service rate, and 1 denotes the number of servers in the queuing system, the queuing delay Lq\textit{L}_{q} is discussed in [37]. The authors assumed a service rate of 740 kpps (where kpps stands for kilo packets per second), an arrival rate of 726 kpps, a link data rate of 9 Gbps, and a packet size of 1,500 bytes with an extra 25 bytes of header and footer and calculated Lq\textit{L}_{q} to be 4 ms. Based on this calculation, we assume the queuing delay as 4 ms for 10 Gbps data rate communication.
   The processing delay Lproc\textit{L}_{proc} is the time taken by a router to process a packet to make an appropriate routing and switching decision before sending the packet to the next appropriate hop. According to [20], when a router has enough processing capability, it can be assumed that the router has approximately equal queuing and processing delays for simplification, and the authors assumed queuing delays and processing of 4 ms and 6 ms, respectively, in their simulations. For LISLs, it’s reasonable to assume that the processing as 6 ms under the assumptions of a congestion-free FSOSN having LISLs with very high data rates as 10 Gbps.
    In this work, we consider the node delay or node latency Tnode\textit{T}_{node} as the sum of queuing and processing delays, since the transmission delay is negligible for a 10 Gbps link data rate, and we assume a node delay of 10 ms for each satellite in the FSOSN, i.e.,

Tnode=Lq+Lproc=10 ms.\textit{T}_{node}=\textit{L}_{q}+\textit{L}_{proc}=10\text{ ms}\text{.} (34)

The LEO satellites in upcoming satellite constellations, such as Starlink, Kuiper, and Telesat, will have the packet switching and routing capability [3, 4, 5, 6], and the queuing and processing delays will indeed be investigated in the node latency model.

4 Methodology

In this section, we describe the Starlink Phase 1 Version 3 and Kuiper Shell 2 constellations and present the methodology for the calculation of the LISL range, ground station range, shortest path for an inter-continental connection over the FSOSN, network latency, and average satellite transmission power.

4.1 Overview of Starlink and Kuiper Constellations

We introduce the two satellite constellations, i.e., Starlink Phase 1 Version 3 [4] and Kuiper Shell 2 [5] constellations, as described in their FCC filings as of 2020 and 2019, respectively. For a satellite constellation, its constellation orbit parameters including number of orbital planes, number of satellites on each plane, and inclination angle can also affect network latency. Greater number of orbital planes and satellites on each plane lead to denser distribution of satellites within the constellation, thus more satellites can be connected for the same LISL range, thereby more optimal links established to achieve the shortest path over an inter-continental connection and results less propagation delay added to the network latency. Similarly, less inclination angle leads to denser satellite distribution in the space above the areas away from the pole, and less propagation delay for inter-continental connections across the low latitude areas.

4.1.1 Starlink Phase 1 Version 3 Constellation

The latest version of the Starlink Phase 1 constellation contains 1,584 satellites and has 22 orbital planes with 72 satellites in each plane. Its altitude is 550 km, and its inclination is 53°. We assume the phasing parameter as 17 to prevent intra-constellation satellite collisions as per [38]. The Walker notation for this constellation is 53°:1,584/22/17, and Fig. 3 shows the constellation pattern. Table 1 provides parameters for different versions of the Starlink Phase 1 constellation since its original proposal in 2016.

Refer to caption

Figure 3: Starlink phase 1 version 3 constellation.
Table 1: Summary of Parameters for Different Versions of the Starlink Phase 1 Constellation.
  Parameter
Version 1
[37]
Version 2
[3]
Version 3
[4]
  Number of Orbital Planes 32 24 22
   Number of Satellites per Plane 50 66 72
  Total Satellites 1,600 1,584 1,584
  Altitude (km) 1,150 550 550
  Inclination (°) 53 53 53

4.1.2 Kuiper Shell 2 Constellation

The second shell of the Kuiper constellation—Kuiper Shell 2—is the biggest of its three shells. This constellation contains 1,296 satellites in total with 36 satellites distributed uniformly on each of the 36 orbital planes. The altitude of this constellation is 610 km, and the inclination is 42°. Based on the analysis in [38], the most suitable value of the phasing parameter for intra-constellation satellite collision avoidance is 11 for Kuiper Shell 2 constellation. Its Walker constellation is 42°:1,296/36/11, and Fig. 4 shows the constellation pattern. Table 2 presents the main parameters of the different shells of the Kuiper constellation [5].

Refer to caption

Figure 4: Kuiper shell 2 constellation.
Table 2: Summary of Parameters for Different Shells of Kuiper Constellation.
  Parameter Shell 1 Shell 2 Shell 3
  Number of Orbital Planes 28 36 34
  Number of Satellites per Plane 28 36 34
  Total Satellites per Shell 784 1,296 1,156
  Altitude (km) 590 610 630
  Inclination (°) 33 42 51.9

4.2 LISL Range Calculation

To study the tradeoff between the average satellite transmission power and network latency in FSOSNs, the LISL range is a critical factor to consider. The LISL range affects the connectivity of satellites, since a satellite can only establish LISLs with other satellites that are within this range. The maximum LISL range for satellites in a constellation can be regarded as the range which is constrained only by visibility, and it can be calculated as [2]

DS=2(r+hS)2(r+ha)2,\textit{D}_{S}=2{\sqrt{(\textit{r}+\textit{h}_{S})^{2}-(\textit{r}+\textit{h}_{a})^{2}}}\text{,} (35)

where r is the radius of the Earth, which is generally assumed to be 6,378 km, hS\textit{h}_{S} is the altitude of the satellite, ha\textit{h}_{a} is the height of the Earth’s atmosphere, which is generally considered to be 80 km, and DS\textit{D}_{S} is the maximum LISL range. For the Starlink Phase 1 Version 3 constellation, hS\textit{h}_{S} is 550 km, and the maximum LISL range can be calculated as 5,016 km. Similarly, for the Kuiper Shell 2 constellation with an altitude of 610 km, the maximum LISL range is 5,339 km.
    To study the tradeoff between network latency and average satellite transmission power, we consider several LISL ranges for satellites in FSOSNs for the two constellations. For the Starlink Phase 1 Version 3 constellation, we assume different LISL range values, specifically 1,575 km, 1,731 km, 2,000 km, 3,000 km, 4,000 km, and 5,016 km. The 1,575 km LISL range is a reasonable minimum feasible LISL range for a satellite in this constellation, as a satellite can establish six permanent LISLs with other satellites at this range, including four with neighbors in the same orbital plane and two with its nearest left and right neighbors in adjacent orbital planes. Shorter LISL ranges are not sufficient for a satellite to establish the two permanent LISLs with neighbors in adjacent orbital planes. At the 1,731 km LISL range, a satellite can establish ten permanent LISLs with other satellites, including four with neighbors in the same orbital plane and six with its nearest neighbors in adjacent orbital planes.
    For the Kuiper Shell 2 constellation, we also consider different LISL ranges, in this case 1,515 km, 2,000 km, 3,000 km, 4,000 km, and 5,339 km. Due to its different design, the Kuiper Shell 2 constellation has a different orbital plane distribution and a different satellite distribution per orbital plane than the Starlink Phase 1 Version 3 constellation; a satellite in the former constellation needs a longer LISL range to establish permanent LISLs with neighbors in same orbital plane as compared to the LISLs with nearest neighbors in the adjacent orbital planes. At the 1,515 km LISL range, a satellite in the Kuiper Shell 2 constellation can set up eight permanent LISLs, including two with satellites in the same orbital plane and six with satellites in adjacent orbital planes. For the Kuiper Shell 2 constellation, 1,515 km is a reasonable minimum LISL range for satellites to establish a reasonable minimum connectivity, since shorter LISL ranges are not long enough for a satellite to establish permanent LISLs with neighbors in same orbital plane.

4.3 Ground Station Range Calculation

The connectivity for the uplink and downlink is affected by the range of the ground station, as the greater the ground station range, the longer the propagation distance for uplink and downlink. The distance of the link between a ground station and satellite can be calculated using (8). The uplink/downlink is longest when the elevation angle between the ground station and an ingress/egress satellite is at its minimum, and the maximum range of the ground station can be calculated for the minimum elevation angle for each constellation. For the Starlink Phase 1 Version 3 constellation, the minimum elevation angle is 25°, and the corresponding maximum range of a ground station is 1,123 km. For the Kuiper Shell 2 constellation, the minimum elevation angle is 35°, and the maximum ground station range is 1,412 km.

4.4 Satellite Orbital Period Calculation

Network latency and average satellite transmission power should be studied based on the entire orbital period of satellites in an FSOSN, which is the time a satellite takes to complete one full orbit around the Earth. The orbital period of a satellite in an FSOSN can be calculated using [40]

T=2π(r+hS)3/(GME),\textit{T}=2\textit{$\pi$}{\sqrt{(\textit{r}+\textit{h}_{S})^{3}/(\textit{G}\textit{M}_{E})}}\text{,} (36)

where r is the radius of the Earth, which is 6,378 km, hS\textit{h}_{S} is the altitude of the satellite in the FSOSN, G is the gravitational constant, which is 6.673×10116.673\times 10^{-11} Nm2/kg2, and ME\textit{M}_{E} is the mass of the Earth, which is 5.98×10245.98\times 10^{24} kg. Based on the altitude of the Starlink Phase 1 Version 3 and Kuiper Shell 2 constellations, which are 550 km and 610 km, respectively, the orbital period for a satellite in these constellations is 5,736 seconds and 5,810 seconds, respectively.

4.5 Outage Probability Calculation

The outage probability (OP) is defined as the probability that the signal-to-noise (SNR) γth\textit{$\gamma$}_{th} can fall below a certain threshold with acceptable communication quality. The OP in this case can be defined as

Pout=Pr{γγth},\textit{P}_{out}=Pr\{\textit{$\gamma$}\leq\textit{$\gamma$}_{th}\}\text{,} (37)

Substituting (26) into (38), the OP can be expressed as [49]

Pout=(1exp((γth/(ηLa)2)γ¯)β/2)α,\textit{P}_{out}=(1-\exp{(-(\textit{$\gamma$}_{th}/(\textit{$\eta$}\textit{L}_{a})^{2})\bar{\textit{$\gamma$}})^{\textit{$\beta$}/2}})^{\textit{$\alpha$}}\text{,} (38)

where γ¯\bar{\textit{$\gamma$}} is the average SNR, and La\textit{L}_{a} is the atmosphere attenuation loss for uplink/downlink. In this work, we investigate the impact of cloud condition on the atmosphere attenuation and turbulence, and therefore the outage probability of optical uplink/downlink communication. We assume three cloud condition scenarios, the first one is clear weather condition with a thin cirrus cloud concentration of 0.5 cm-3 and liquid water content as 3.128×1043.128\times 10^{-4} g/m-3, while the second one is clear weather condition with a cirrus cloud concentration of 0.025 cm-3 and liquid water content as 0.06405 g/m-3, and the third one is cloudy weather condition with a Cumulus cloud concentration of 250 cm-3 and liquid water content as 1.0 g/m-3, according to [10]. Based on these three different weather conditions, we calculated the OP using practical parameters and we hope these deployment scenarios can be used to provide practical insights about the system architecture, so that the reader can get a quick idea about the overall system performance.

4.6 Shortest Path Calculation

To study the tradeoff between the average satellite transmission power and network latency in an FSOSN at different LISL ranges, we need to calculate the shortest path between the source ground station and the destination ground station over the FSOSN at an LISL range at each time slot to calculate the network latency and average satellite transmission power for that shortest path. The shortest path can be determined based on adjacency matrix containing propagation distances of laser links between satellites and between satellites and ground stations using Dijkstra’s algorithm [41]. This algorithm uses weights of the links to determine the shortest path, and these weights are the propagation distances of the optical links in our case. The following steps are used in Dijkstra’s algorithm to find the shortest path for an inter-continental connection from a source ground station to a destination ground station over an FSOSN.
1) Create a shortest path set S that keeps track of satellites that are on the shortest path, another set X for satellites that are not on the shortest path, and a vector D to track the shortest path distance for each satellite to the source ground station gs\textit{g}_{s}.
2) Assign an initial weight to all the links as their propagation distance, add gs\textit{g}_{s} to S, and update the distance for gs\textit{g}_{s} in D as 0, while the distances for other satellites to gs\textit{g}_{s} are infinite in D.
3) While S does not include the destination ground station gd\textit{g}_{d}, go through the following loop.
a) Pick the satellite q with the minimum distance to gs\textit{g}_{s}, which is the closest satellite to gs\textit{g}_{s} based on its distance to gs\textit{g}_{s}.
b) Mark q as visited and include q to S. Update the distance for q in D.
c) For every adjacent satellite pi\textit{p}_{i} of q, check its distance to q. Since Dijkstra’s algorithm keeps tracking the currently known shortest path for each satellite from the source ground station and updates the distance if it finds a shorter one, check the new path to gs\textit{g}_{s} through q from pi\textit{p}_{i}. For each pi\textit{p}_{i}, if the currently known distance of pi\textit{p}_{i} to gs\textit{g}_{s} is greater than the sum of the distance of q to gs\textit{g}_{s} and the distance between q and pi\textit{p}_{i}, then update the distance for pi\textit{p}_{i} to gs\textit{g}_{s}. After iterating through all the adjacent satellites pi\textit{p}_{i} of q, go back to step a to find the next satellite to add to S.

4.7 Network Latency Calculation

We calculate the network latency Tnet\textit{T}_{net} of a shortest path at a time slot as the end-to-end latency of the shortest path between the source ground station and the destination ground station by adding the link latency (or propagation delay) of each laser link, including the uplink, downlink, and ISLs, and the node delay (or node latency) for each satellite on the shortest path, which can be expressed as

Tnet=(Tup+k=1n1Tk+Tdown)+n×Tnode,\textit{T}_{net}=(\textit{T}_{up}+\sum_{k=1}^{n-1}\textit{T}_{k}+\textit{T}_{down})+\textit{n}\texttimes\textit{T}_{node}\text{,} (39)

where Tup\textit{T}_{up} and Tdown\textit{T}_{down} are the link latencies of the laser uplink and downlink, respectively, and Tk\textit{T}_{k} is the link latency (or propagation delay) for kth ISL. Since there are n satellites on the shortest path for the inter-continental connection, n1\textit{n}-1 ISLs exist in total on the shortest path. As mentioned in Section III, Tnode\textit{T}_{node} is the node delay for each satellite and is assumed to be 10 ms.

The propagation delay for each link on the shortest paths for first five time slots for the FSOSN with the Starlink Phase 1 Version 3 constellation at the 3,000 km LISL range for the Toronto–Sydney inter-continental connection is shown in Table 3. Tup\textit{T}_{up} and Tdown\textit{T}_{down} represent the propagation delay for the uplink and downlink, and Tk\textit{T}_{k} is the propagation delay for the kth ISL. For example, the network latency of the shortest path at the first time slot is equal to the sum of Tup\textit{T}_{up}, Tdown\textit{T}_{down}, the propagation delays of all seven LISLs, and Tnode\textit{T}_{node} for each satellite on the shortest path (i.e., 20 ms × 8), and is equal to 137.22 ms.

Table 3: Propagation Delay for Each Link on the Shortest Paths at the First Five Time Slots for the FSOSN with the Starlink Phase 1 Version 3 Constellation at 3,000 km LISL Range for the Toronto–Sydney Inter-Continental Connection.
   Time Slot
Tup\textit{T}_{up}
(ms)
T1\textit{T}_{1}
(ms)
T2\textit{T}_{2}
(ms)
T3\textit{T}_{3}
(ms)
T4\textit{T}_{4}
(ms)
T5\textit{T}_{5}
(ms)
T6\textit{T}_{6}
(ms)
T7\textit{T}_{7}
(ms)
Tdown\textit{T}_{down}
(ms)
Tnet\textit{T}_{net}
(ms)
  1 3.23 8.04 3.94 9.40 9.35 5.27 5.13 9.28 3.53 137.22
  2 3.22 8.05 3.95 9.38 9.38 5.27 5.13 9.28 3.53 137.22
  3 3.21 8.06 3.95 9.36 9.41 5.27 5.13 9.27 3.53 137.22
  4 3.20 8.06 3.95 9.34 9.44 5.27 5.13 9.27 3.53 137.23
  5 3.19 8.07 3.96 9.33 9.47 5.27 5.13 9.26 3.52 137.23

The mean network latency T¯net\textit{$\overline{T}$}_{net} is the mean of the network latency of the shortest paths at all time slots, and it can be expressed as

T¯net=i=1jTnet,i/j,\textit{$\overline{T}$}_{net}=\sum_{i=1}^{j}\textit{T}_{net,i}/\textit{j}\text{,} (40)

where j is the total number of time slots, and Tnet,i\textit{T}_{net,i} is the network latency of the shortest path at the ith time slot. For example, T¯net\textit{$\overline{T}$}_{net} for the first five time slots in Table 3 is equal to 137.22 ms.

4.8 Average Satellite Transmission Power Calculation

The transmission power of each satellite is calculated by adding the transmission power required to establish the incoming link and outgoing link for that satellite. For the ingress (or first) and egress (or last) satellites on a shortest path, the incoming link and outgoing link are the uplink and downlink, respectively, and the transmission power for such satellites can be calculated as

PTS,ingress/egress=PR/(ηTηRGTGRLTLRLALPG,up/down)+PR/(ηTηRGTGRLTLRLPS,out/in),\textit{P}_{TS,ingress/egress}=\\ \textit{P}_{R}/(\textit{$\eta$}_{T}\textit{$\eta$}_{R}\textit{G}_{T}\textit{G}_{R}\textit{L}_{T}\textit{L}_{R}\textit{L}_{A}\textit{L}_{PG,up/down})\\ +\textit{P}_{R}/(\textit{$\eta$}_{T}\textit{$\eta$}_{R}\textit{G}_{T}\textit{G}_{R}\textit{L}_{T}\textit{L}_{R}\textit{L}_{PS,out/in})\text{,} (41)

where LPG,up\textit{L}_{PG,up} and LPG,down\textit{L}_{PG,down} is the free-space path loss for the uplink and downlink, respectively, and LPS,in\textit{L}_{PS,in} and LPS,out\textit{L}_{PS,out} is the free-space path loss for the incoming and outgoing ISLs, respectively. The optical links between satellites, i.e., LISLs, and optical links between ground stations and satellites (uplink and downlink) are bidirectional [42],[43]. That is why, we are considering the power from the satellite to the ground station for uplink and downlink for modelling transmission power of ingress and egress satellites. For other (or intermediate) satellites on the shortest path, the incoming and outgoing links are ISLs, and the transmission power for such satellites can be calculated as

PTS,intermediate=PR/(ηTηRGTGRLTLRLPS,in)+PR/(ηTηRGTGRLTLRLPS,out).\textit{P}_{TS,intermediate}=\\ \textit{P}_{R}/(\textit{$\eta$}_{T}\textit{$\eta$}_{R}\textit{G}_{T}\textit{G}_{R}\textit{L}_{T}\textit{L}_{R}\textit{L}_{PS,in})\\ +\textit{P}_{R}/(\textit{$\eta$}_{T}\textit{$\eta$}_{R}\textit{G}_{T}\textit{G}_{R}\textit{L}_{T}\textit{L}_{R}\textit{L}_{PS,out})\text{.} (42)

The transmission power for each laser link on the shortest paths for the first five time slots for the FSOSN with the Starlink Phase 1 Version 3 constellation at an LISL range of 3,000 km for the Toronto–Sydney inter-continental connection is shown in Table 4. The satellite transmission power and average satellite transmission power for these shortest paths are given in Table 5. In Table 4, Pup\textit{P}_{up} and Pdown\textit{P}_{down} represent the transmission power for the uplink and downlink, respectively, and Pi\textit{P}_{i} is the transmission power for the ith ISL. For example, P1\textit{P}_{1} is the transmission power for the first ISL on the shortest path at a time slot. In Table 5, the transmission power for the first satellite (or ingress satellite) PTS,ingress\textit{P}_{TS,ingress} at the first time slot can be calculated by adding Pup\textit{P}_{up} and P1\textit{P}_{1} (i.e., 89.34 mW and 198.26 mW, respectively) and is equal to 287.60 mW.

Table 4: Transmission Power for Uplink/Downlink and Inter-Satellite Links on the Shortest Paths at the First Five Time Slots for the FSOSN with the Starlink Phase 1 Version 3 Constellation at 3,000 km LISL Range for the Toronto–Sydney Inter-Continental Connection.
   Time Slot
Pup\textit{P}_{up}
(mW)
P1\textit{P}_{1}
(mW)
P2\textit{P}_{2}
(mW)
P3\textit{P}_{3}
(mW)
P4\textit{P}_{4}
(mW)
P5\textit{P}_{5}
(mW)
P6\textit{P}_{6}
(mW)
P7\textit{P}_{7}
(mW)
Pdown\textit{P}_{down}
(mW)
  1 70.42 198.26 47.67 270.82 268.28 85.14 80.80 264.18 111.49
  2 69.93 198.61 47.76 269.78 269.90 85.14 80.75 263.88 111.31
  3 69.45 198.96 47.86 268.74 271.53 85.13 80.71 263.59 111.13
  4 68.99 199.31 47.95 267.71 273.17 85.12 80.67 263.29 110.97
  5 68.53 199.66 48.04 266.69 274.81 85.12 80.63 262.99 110.82
Table 5: Transmission Power for Each Satellite and Average Satellite Transmission Power on the Shortest Paths at the First Five Time Slots for the FSOSN with the Starlink Phase 1 Version 3 Constellation at 3,000 km LISL Range for the Toronto–Sydney Inter-Continental Connection.
   Time Slot
PTS,ingress\textit{P}_{TS,ingress}
(mW)
PTS,1\textit{P}_{TS,1}
(mW)
PTS,2\textit{P}_{TS,2}
(mW)
PTS,3\textit{P}_{TS,3}
(mW)
PTS,4\textit{P}_{TS,4}
(mW)
PTS,5\textit{P}_{TS,5}
(mW)
PTS,6\textit{P}_{TS,6}
(mW)
PTS,egress\textit{P}_{TS,egress}
(mW)
PTS,avg\textit{P}_{TS,avg}
(mW)
  1 268.67 245.93 318.49 539.10 353.42 165.94 344.98 375.67 326.53
  2 268.54 246.37 317.54 539.68 355.04 165.89 344.64 375.19 326.61
  3 268.41 246.81 316.60 540.27 356.66 165.84 344.30 374.72 326.70
  4 268.29 247.26 315.66 540.88 358.29 165.79 343.96 374.26 326.80
  5 268.18 247.70 314.73 541.50 359.93 165.74 343.62 373.81 326.90

The average satellite transmission power for satellites on a shortest path at a time slot PTS,avg\textit{P}_{TS,avg} is the average value of the transmission power of all satellites on the shortest path at that time slot, and it can be expressed as

PTS,avg=(PTS,ingress+m=1n2PTS,m+PTS,egress)/n,\textit{P}_{TS,avg}=(\textit{P}_{TS,ingress}+\sum_{m=1}^{n-2}\textit{P}_{TS,m}+\textit{P}_{TS,egress})/\textit{n}\text{,} (43)

where PTS,m\textit{P}_{TS,m} is the transmission power for the mth intermediate satellite on the shortest path and there are n satellites on the shortest path, and thereby n2\textit{n}-2 intermediate satellites exist in total on the shortest path. The mean of the average satellite transmission power P¯TS,avg\textit{$\overline{P}$}_{TS,avg}, i.e., the mean of the average satellite transmission power of the shortest paths at all time slots, can be expressed as

P¯TS,avg=i=1jPTS,avg,i/j,\textit{$\overline{P}$}_{TS,avg}=\sum_{i=1}^{j}\textit{P}_{TS,avg,i}/\textit{j}\text{,} (44)

where j is the total number of time slots, and PTS,avg,i\textit{P}_{TS,avg,i} is average satellite transmission power at the ith time slot. For instance, PTS,1\textit{P}_{TS,1} is the satellite transmission power for the first intermediate satellite on the shortest path at a time slot in Table 5. The average satellite transmission power PTS,avg\textit{P}_{TS,avg} of the shortest path at the first time slot can be calculated by taking the average value of the transmission powers for the eight satellites on this shortest path, which is equal to 328.89 mW in Table 5. Finally, P¯TS,avg\textit{$\overline{P}$}_{TS,avg} for the first five time slots in Table 5 is equal to 329.02 mW.

5 Power versus Latency Tradeoff

In this work, we assume there is one source ground station in Toronto, and a destination ground station in each of the other cities: Sydney, Istanbul, and London. We then determine the shortest paths for the three inter-continental connections from the source ground station to the destination ground station for the FSOSNs with the two constellations at an LISL range at each time slot. Subsequently, we calculate the network latency and average satellite transmission power for the shortest path at each time slot. In doing so, we are interested in studying the tradeoff between network latency and average satellite transmission power, and we consider several LISL ranges for satellites in each FSOSN as described in Section IV. Finally, we investigate the effect of atmosphere attenuation and turbulence on outage probability of optical uplink/downlink.

5.1 Simulation Methodology

For the computation of the average satellite transmission power and network latency of the shortest path for an inter-continental connection, we first simulate the Starlink Phase 1 Version 3 and Kuiper Shell 2 constellations using the well-known satellite constellation simulator Systems Tool Kit (STK) version 12.1 [44]. We generate a satellite constellation, ground stations, and all possible ISLs and links between satellites and ground stations based on an LISL range within STK, and then we extract the data containing the links between satellites as well as the links between satellites and ground stations and the positions of satellites and ground stations from STK. We process this data in Python to compute the distance for each link at a time slot to find the shortest path between ground stations for an inter-continental connection over the FSOSN at each time slot. To this end, we first discretize the extracted data, as the data obtained from STK is continuous, and for a certain link, it contains the duration of the link for the entire orbital period. We discretize this data by dividing it into time slots and reorganize it so that it shows all links that exist at a particular time slot. The position data for each satellite is also discretized into time slots to match the discretized link data. Next, using the positions of satellites and links between them at each time slot, we compute the distance for a link at each time slot. Finally, we find the shortest path between the source and destination ground stations for an inter-continental connection over the FSOSN using the NetworkX library in Python [45].

5.2 Simulation Setup

To study the tradeoff between network latency and satellite transmission power in FSOSNs with the Starlink Phase 1 Version 3 and Kuiper Shell 2 constellations, the period of the simulation is taken as 6,000 time slots for both constellations. This encompasses the entire orbital period of a satellite in these FSOSNs, and the duration of a time slot is set to one second.

To calculate satellite transmission power in the two FSOSNs, we consider the laser link model parameters summarized in Table 6. These parameters are used in existing or practical laser satellite communication systems. We consider On-Off Keying (OOK) as the laser link’s modulation scheme. We assume the wavelength as 1,550 nm for both LISL and laser uplink/downlink, since it is widely used for optical satellite communication due to its better performance. We take both transmitter and receiver optical efficiency as 0.8, as this is the common value of these parameters for laser communication terminals (LCTs) and is used in many simulation studies as well. We assume the data rate as 10 Gbps since this is a practical link data rate for optical satellite communications and the maximum data rate that is offered by Mynaric’s LCT [46]. The transmitting divergence angle is set to 15 μ\murad and the receiver telescope diameter is 80 mm according to Mynaric’s LCT. We assume both transmitter and receiver pointing errors as 1 μ\murad, since these values are widely used in LCTs. We use a curve-fitting technique to find the receiver sensitivity of -35.5 dBm for the OOK modulation with a 10 Gbps data rate and 101210^{-12} BER from [48]. We assume the ground station altitude as 0.1 km, since this is the general height of tall buildings and towers. We assume clear weather conditions with a thin cirrus cloud concentration of 0.5 cm-3 and set liquid water content as 3.128×1043.128\times 10^{-4} g/m-3 according to [10]. We take 20 km as the height of the troposphere layer [50], since this value is practical in laser uplink/downlink. We consider no special weather condition for all the five inter-continental connection pairs. We consider the LM for inter-satellite links as 3 dB and the LM for uplink/downlink as 6 dB, since there is more turbulence and attenuation in uplink/downlink. This yields a received power PR\textit{P}_{R} of -32.5 dBm for ISLs and -29.5 dBm for uplink/downlink according to (23). To calculate the OP for uplink/downlink, we assume the given SNR threshold as 7 dB, and elevation angle as 30 degrees. As mentioned earlier, we consider three cloud scenarios, each characterized by different parameters: thin cirrus clouds, cirrus clouds, and cumulus cloud concentrations.

Table 6: Simulation Parameters.
Parameter Symbol Units Value
Laser wavelength [46] λ\lambda nm 1550
Transmitter optical efficiency [47] ηT\textit{$\eta$}_{T} 0.8
Receiver optical efficiency [47] ηR\textit{$\eta$}_{R} 0.8
Data rate [46] Rdata\textit{R}_{data} Gbps 10
Receiver telescope diameter [46] DR\textit{D}_{R} mm 80
Transmitter pointing error [47] θT\textit{$\theta$}_{T} μ\murad 1
Receiver pointing error [47] θR\textit{$\theta$}_{R} μ\murad 1
Full transmitting divergence
angle [48]
ΘT\textit{$\Theta$}_{T} μ\murad 1.5
Receiver sensitivity [48] Preq\textit{P}_{req} dBm -35.5
Bit error rate [48] 10-12
Ground station attitude hE\textit{h}_{E} km 0.1
Thin cirrus cloud concentration [10] LW\textit{L}_{W} cm-3 0.5
Thin cirrus liquid water content [10] N g/m-3 3.128×10-4
Cirrus cloud concentration [10] LW\textit{L}_{W} cm-3 0.0255
Cirrus liquid water content [10] N g/m-3 0.06405
Cumulus cloud concentration [10] LW\textit{L}_{W} cm-3 250
Cumulus liquid water content [10] N g/m-3 1.0
SNR threshold [49] γth\textit{$\gamma$}_{th} dB 7
Elevation angle θE\textit{$\theta$}_{E} degree 30
Partial size coefficient [49] Φ\Phi 1.6
Troposphere layer height [50] hA\textit{h}_{A} km 20
Atmosphere layer height [2] ha\textit{h}_{a} km 80
Location of source
ground station
Toronto
Location of destination
ground station
{Sydney,
Istanbul,
London}
Node delay Tnode\textit{T}_{node} ms 10
LISL range for Kuiper Shell 2 km
{1515, 2000,
3000, 4000,
5339}
LISL range for Starlink Phase 1
Version 3
km
{1575, 1731,
2000, 3000,
4000, 5016}
Speed of light cs\textit{c}_{s} m/s 299792458
Duration of a time slot s 1
Total number of time slots j 6000

5.3 T¯net\textit{$\overline{T}$}_{net} versus P¯TS,avg\textit{$\overline{P}$}_{TS,avg} for Different LISL Ranges

To show the relationship and tradeoff between network latency and satellite transmission power, we superimpose the two curves representing T¯net\textit{$\overline{T}$}_{net} and P¯TS,avg\textit{$\overline{P}$}_{TS,avg} for different LISL ranges on a unified plot. For different inter-continental connections, the maximum value of T¯net\textit{$\overline{T}$}_{net} (with the minimum corresponding to the P¯TS,avg\textit{$\overline{P}$}_{TS,avg} case) is realized when the LISL range is at its minimum feasible value for each FSOSN. Conversely, the minimum value of T¯net\textit{$\overline{T}$}_{net} (with the maximum corresponding to the P¯TS,avg\textit{$\overline{P}$}_{TS,avg} case) is observed at the maximum feasible LISL range. Within the plot, the curve of T¯net\textit{$\overline{T}$}_{net} gradually descends over LISL ranges, whereas the curve of P¯TS,avg\textit{$\overline{P}$}_{TS,avg} rises with the increasing LISL range. The point of the intersection between these two curves signifies the juncture of relatively optimal performance between network latency and satellite transmission power can be reached at that LISL range for a given inter-continental connection. Any deviation from this identified LISL range could lead to higher energy consumption or worse latency performance, underscoring the critical importance of maintaining this delicate balance in a FSOSN.

Figs. 5 and 6 show the relationship between T¯net\textit{$\overline{T}$}_{net} and P¯TS,avg\textit{$\overline{P}$}_{TS,avg} for different LISL ranges for the two FSOSNs for the Toronto–Sydney inter-continental connection. For the FSOSN with the Starlink Phase 1 Version 3 constellation, the curves of T¯net\textit{$\overline{T}$}_{net} and P¯TS,avg\textit{$\overline{P}$}_{TS,avg} intersect when the LISL range is approximately 2,900 km, which indicates that T¯net\textit{$\overline{T}$}_{net} and P¯TS,avg\textit{$\overline{P}$}_{TS,avg} are balanced when the LISL range is 2,900 km for satellites in this FSOSN for this inter-continental connection, and T¯net\textit{$\overline{T}$}_{net} and P¯TS,avg\textit{$\overline{P}$}_{TS,avg} are around 135 ms and 380 mW, respectively, at this intersection. For the FSOSN with the Kuiper Shell 2 constellation, T¯net\textit{$\overline{T}$}_{net} and P¯TS,avg\textit{$\overline{P}$}_{TS,avg} intersect when the LISL range is around 3,800 km for this inter-continental connection, and these two parameters are approximately 120 ms and 700 mW, respectively.

In Figs. 7 and 8, we can see the tradeoff between T¯net\textit{$\overline{T}$}_{net} and P¯TS,avg\textit{$\overline{P}$}_{TS,avg} at different LISL ranges in the two FSOSNs for the Toronto–Istanbul inter-continental connection, and Figs. 9 and 10 show this tradeoff in the two FSOSNs for the Toronto–London inter-continental connection. For the Toronto–Istanbul inter-continental connection, T¯net\textit{$\overline{T}$}_{net} and P¯TS,avg\textit{$\overline{P}$}_{TS,avg} intersect when the LISL range is about 2,600 km for the FSOSN with the Starlink Phase 1 Version 3 constellation and 2,900 km for the FSOSN with the Kuiper Shell 2 constellation. For the Toronto to London inter-continental connection, these two parameters intersect when the LISL ranges are 3,400 km and 3,000 km for the FSOSNs with the Starlink Phase 1 Version 3 and Kuiper Shell 2 constellations, respectively. For both FSOSNs in all inter-continental connections, as the LISL range increases, higher transmission power is needed for a satellite but lower network latency is obtained. For T¯net\textit{$\overline{T}$}_{net} curve in Fig. 9, it changes sightly between 3000 km to 4000 km LISL range, which is different from T¯net\textit{$\overline{T}$}_{net} at the same LISL range for Kuiper Shell 2 constellation as shown in Fig. 10. This is because that the distribution of the satellites for Starlink Phase 1 Version 3 constellation is denser at high latitude comparing to Kuiper Shell 2 constellation, in this way for short inter-continental connections like Toronto to London T¯net\textit{$\overline{T}$}_{net} can change unobtrusively at certain LISL ranges.

Refer to caption

Figure 5: T¯net\textit{$\overline{T}$}_{net} vs. P¯TS,avg\textit{$\overline{P}$}_{TS,avg} at different LISL ranges for the FSOSN with the Starlink Phase 1 Version 3 constellation for the Toronto–Sydney inter-continental connection.

Refer to caption

Figure 6: T¯net\textit{$\overline{T}$}_{net} vs. P¯TS,avg\textit{$\overline{P}$}_{TS,avg} at different LISL ranges for the FSOSN with the Kuiper Shell 2 constellation for the Toronto–Sydney inter-continental connection.

Refer to caption

Figure 7: T¯net\textit{$\overline{T}$}_{net} vs. P¯TS,avg\textit{$\overline{P}$}_{TS,avg} at different LISL ranges for the FSOSN with the Starlink Phase 1 Version 3 constellation for the Toronto–Istanbul inter-continental connection.

Refer to caption

Figure 8: T¯net\textit{$\overline{T}$}_{net} vs. P¯TS,avg\textit{$\overline{P}$}_{TS,avg} at different LISL ranges for the FSOSN with the Kuiper Shell 2 constellation for Toronto–Istanbul inter-continental connection.

Refer to caption

Figure 9: T¯net\textit{$\overline{T}$}_{net} vs. P¯TS,avg\textit{$\overline{P}$}_{TS,avg} at different LISL ranges for the FSOSN with the Starlink Phase 1 Version 3 constellation for the Toronto–London inter-continental connection.

Refer to caption

Figure 10: T¯net\textit{$\overline{T}$}_{net} vs. P¯TS,avg\textit{$\overline{P}$}_{TS,avg} at different LISL ranges for the FSOSN with the Kuiper Shell 2 constellation for the Toronto–London inter-continental connection.

5.4 Outage Probability Analysis

Figs. 11 and 12 show the relationship between OP and SNR for optical uplink and downlink, respectively. Based on Table 6, we consider three different weather conditions as thin cirrus, cirrus and cumulus clouds. We assume the average SNR as 30 dB, the SNR threshold as 7 dB, the elevation angle of the ground station as 30 degrees and employing Starlink Phase 1 Version 3 constellation, we plot the two figures. As depicted in Figs. 11 and 12, the OP decreases as the SNR increases under cirrus and thin cirrus cloud conditions. Notably, the thin cirrus cloud condition exhibits better performance compared to the cirrus cloud condition in terms of OP. The OP under cumulus cloud condition is 1.0, which indicates an unapproachable optical communication. Comparing the curves of the two figures under same cloud condition, we find that the OP performance for uplink is better than downlink since less atmosphere attenuation ocurrs through uplink. For upink communication, SNR should reach 45 dB to obtain 10810^{-8} OP performance under thin cirrus cloud condition; while for downlink communication, the requirement of SNR to achieve the same OP performance is 57 dB for thin cirrus cloud condition.

Refer to caption

Figure 11: Outage probability vs. P/N0\textit{N}_{0} under different weather conditions for optical uplink communications.

Refer to caption

Figure 12: Outage probability vs. P/N0\textit{N}_{0} under different weather conditions for optical downlink communications.

5.5 Insights

Here we provide important practical insights for FSOSNs based on the analysis in this work.

\bullet For LISL and laser uplink/downlink, the network latency for the shortest path for an inter-continental connection between two ground stations and the average satellite transmission power for the shortest path have an inverse relationship. As the LISL range increases, the links on the shortest path become longer, the propagation distance of the laser links increases, and the transmission power for the laser links increases along with the average satellite transmission power. At the same time, there are fewer direct links and nodes on the shortest path, and accordingly, the total propagation delay of the path decreases, the total node delay of the path decreases, and the network latency of the path decreases. Therefore, according to this reverse relationship, in order to reach less network latency to achieve better network performance, satellites needed to be designed for high energy payload to sustain higher transmission power consumption. This can be accomplished by implementing larger capacity battery, applying dynamic ISLs instead of establishing the links throughout the entire working period.

\bullet Based on the results of the tradeoff analysis between T¯net\textit{$\overline{T}$}_{net} and P¯TS,avg\textit{$\overline{P}$}_{TS,avg} for different LISL ranges for the two FSOSNs for different inter-continental connections, we find that there is no clear relationship between the LISL range when T¯net\textit{$\overline{T}$}_{net} and P¯TS,avg\textit{$\overline{P}$}_{TS,avg} intersect and the length of the inter-continental connection (i.e., the end-to-end distance between the source and destination ground stations for an inter-continental connection over an FSOSN) within the same FSOSN. A shorter inter-continental connection may have a longer LISL range when T¯net\textit{$\overline{T}$}_{net} and P¯TS,avg\textit{$\overline{P}$}_{TS,avg} intersect. As we can observe that there are different LISL ranges where T¯net\textit{$\overline{T}$}_{net} and P¯TS,avg\textit{$\overline{P}$}_{TS,avg} reach a balance for different inter-continental connections over an FSOSN, all inter-continental connections over an FSOSN must be considered simultaneously while determining an appropriate LISL range for the FSOSN to balance the average satellite transmission power and network latency.

\bullet In laser uplink/downlink communications, the performance of the outage probability is inevitably influenced by cloud conditions. In the case of adverse cloud conditions, such as cumulus clouds, the connectivity of the laser link becomes unavailable, irrespective of the SNR. Under these circumstances, it becomes necessary to transmit through alternative laser links between the ground station and satellites, circumventing the area affected by poor cloud conditions, even though this may cost more network latency and higher satellite transmission power. For laser downlink communication, higher SNR is required to sustain an acceptable OP comparing to uplink communication, and more satellite transmission power will be consumed.

6 Conclusion

In this work, we investigated the tradeoff between satellite transmission power and network latency in FSOSNs resulting from the Starlink Phase 1 Version 3 and Kuiper Shell 2 constellations. We used appropriate system models for laser link transmission power and latency and studied T¯net\textit{$\overline{T}$}_{net} and P¯TS,avg\textit{$\overline{P}$}_{TS,avg} in the two FSOSNs for different inter-continental connections and different LISL ranges. We further studied the effect of atmosphere attenuation and turbulence on outage probability of laser uplink/downlink.

We found that as the LISL range increases, P¯TS,avg\textit{$\overline{P}$}_{TS,avg} increases while T¯net\textit{$\overline{T}$}_{net} decreases for both FSOSNs in all inter-continental connections. The results show that for the inter-continental connection from Toronto to Sydney in the FSOSN based on the Starlink Phase 1 Version 3 constellation, T¯net\textit{$\overline{T}$}_{net} and P¯TS,avg\textit{$\overline{P}$}_{TS,avg} intersect when the LISL range is approximately 2,900 km, with these parameters being approximately 135 ms and 380 mW, respectively. For the FSOSN based on the Kuiper Shell 2 constellation in this inter-continental connection, this LISL range is around 3,800 km, and the two parameters are approximately 120 ms and 700 mW, respectively. For the Toronto–Istanbul inter-continental connection, these parameters intersect when LISL ranges are about 2,600 km for the FSOSN based on the Starlink Phase 1 Version 3 constellation and 2,900 km for the FSOSN based on the Kuiper Shell 2 constellation. For the Toronto to London inter-continental connection, T¯net\textit{$\overline{T}$}_{net} and P¯TS,avg\textit{$\overline{P}$}_{TS,avg} intersect when LISL ranges are 3,400 km and 3,000 km for the FSOSNs for the Starlink Phase 1 Version 3 and Kuiper Shell 2 constellation, respectively. The OP decreases as the SNR increases under cirrus and thin cirrus cloud conditions, while the thin cirrus cloud condition performs better compared to the cirrus cloud condition. For cumulus cloud condition the connectivity of the laser link is unreachable. Furthermore, the OP performance for uplink is better than downlink since less atmosphere attenuation ocurrs through uplink.

In a future work, we will investigate the problem of determining an appropriate LISL range for balancing network latency and average satellite transmission power in an FSOSN. This will require the problem to be formulated as a multi-objective mathematical programming problem. We will further examine the problem of determining an appropriate LISL range for an FSOSN for balancing average satellite transmission power and network latency over multiple inter-continental connections (multiple source-destination pairs) simultaneously for that FSOSN. A LEO satellite constellation at a higher altitude requires fewer satellites to achieve the same coverage, which can affect the network latency and average satellite transmission power, and we will investigate the effect of satellite constellation altitude on the tradeoff between network latency and satellite transmission power. The power consumption for optical amplification and signal regeneration is also important in terms of the total power consumption, we will examine the tradeoff between network latency and the satellite total power consumption in future. In this work, we haven’t considered the mobility of the satellites and its impact on the tradeoff between network latency and average satellite transmission power, and we will examine this important factor in future. In this work, we don’t consider the effect of different satellite constellation parameters, such as number of orbital planes, number of satellites per plane, inclination angle, etc., on the tradeoff between network latency and average satellite transmission power, and we will examine the effect of these parameters in future.

References

  • [1] A.U. Chaudhry and H. Yanikomeroglu, “Free Space Optics for Next-Generation Satellite Networks,” IEEE Consumer Electronics Magazine, vol. 10(6), pp. 21–31, Nov. 2021.
  • [2] A.U. Chaudhry and H. Yanikomeroglu, “Laser Intersatellite Links in a Starlink Constellation: A Classification and Analysis,” IEEE Vehicular Technology Magazine, vol. 16(2), pp. 48–56, Jun. 2021.
  • [3] SpaceX FCC update, Nov. 2018, “SpaceX Non-Geostationary Satellite System, Attachment A, Technical Information to Supplement Schedule S,” [Online]. Available: https://licensing.fcc.gov/myibfs/download.do?attachment_key=1569860, Accessed: December 6, 2023.
  • [4] SpaceX FCC update, Apr. 2020, “SpaceX Non-Geostationary Satellite System, Attachment A, Technical Information to Supplement Schedule S,” [Online]. Available: https://licensing.fcc.gov/myibfs/download.do?attachment_key=2274316, Accessed: December 6, 2023.
  • [5] Kuiper Systems LLC FCC update, Jul. 2019, “Technical Appendix, Application of Kuiper Systems LLC for Authority to Launch and Operate a Non-Geostationary Satellite Orbit System in Ka-band Frequencies,” [Online]. Available: https://fcc.report/IBFS/SAT-LOA-20190704-00057/1773885, Accessed: December 6, 2023.
  • [6] Telesat Canada FCC update, May 2020, “Application for Modification of Market Access Authorization,” [Online]. Available: https://fcc.report/IBFS/SATMPL-20200526-00053/2378318.pdf, Accessed: December 6, 2023.
  • [7] J. Liang, A.U. Chaudhry, E. Erdogan, and H. Yanikomeroglu, “Link Budget Analysis for Free-Space Optical Satellite Networks,” in Proc. 2022 IEEE 23rd International Symposium on a World of Wireless, Mobile and Multimedia Networks, Belfast, UK, 2022, pp. 471–476.
  • [8] V. Knap, L. Vestergaard, and D. Stroe, “A Review of Battery Technology in CubeSats and Small Satellite Solutions,” Energies, vol.13(16), article id. 4097, Aug. 2020.
  • [9] A. Dutta and Y. Yechiam, “Power Management of LEOs under Bursty Broadband Traffic,” in Proc. 17th AIAA International Communications Satellite Systems Conference and Exhibit, Yokohama, Japan, 1998, pp. 494–514.
  • [10] A.U. Chaudhry and H. Yanikomeroglu, “Temporary Laser Inter-Satellite Links in Free-Space Optical Satellite Networks,” IEEE Open Journal of the Communications Society, vol. 3, pp. 1413–1427, Aug. 2022.
  • [11] M.S. Awan, E. Leitgeb, B. Hillbrand, F. Nadeem, and M. Khan, “Cloud Attenuations for Free-Space Optical Links,” in Proc. International Workshop on Satellite and Space Communications, Siena, Italy, 2009, pp. 274–278.
  • [12] D. Bhattacherjee, W. Aqeel, G. Laughlin, B.M. Maggs, and A. Singla, “A Bird’s Eye View of the World’s Fastest Networks,” in Proc. ACM Internet Measurement Conference, New York, NY, USA, 2020, pp. 521–527.
  • [13] DARPA, “Blackjack Focuses on Risk Reduction Flights and Simulations to Prepare for Full Demonstration,” May 2020, [Online]. Available: https://www.darpa.mil/news-events/2020-05-11, Accessed: December 6, 2023.
  • [14] J.F. Kurose and K.W. Ross, Computer Networking – A Top Down Approach Featuring the Internet, Addison-Wesley, Boston, 2009.
  • [15] T.S. Abdu, S. Kisseleff, E. Lagunas, and S. Chatzinotas, “Flexible Resource Optimization for GEO Multibeam Satellite Communication System,” IEEE Transactions on Wireless Communications, vol. 20(12), pp. 7888–7902, Dec. 2021.
  • [16] H. Tsuchida, Y. Kawamoto, N. Kato, K. Kaneko, S. Tani, M. Hangai, and H. Aruga, “Improvement of Battery Lifetime Based on Communication Resource Control in Low-Earth-Orbit Satellite Constellations,” IEEE Transactions on Emerging Topics in Computing, vol, 10(3), pp. 1388–1398, July–Sept. 2022.
  • [17] H. Tsuchida, Y. Kawamoto, N. Kato, K. Kaneko, S. Tani, S. Uchida, and H. Aruga, “Efficient Power Control for Satellite-Borne Batteries Using Q-Learning in Low-Earth-Orbit Satellite Constellations,” IEEE Wireless Communications Letters, vol. 9(6), pp. 809–812, Jun. 2020.
  • [18] M. Handley, “Delay is Not an Option: Low Latency Routing in Space,” in Proc. 17th ACM Workshop on Hot Topics in Networks, Redmond, WA, USA, 2018, pp. 85–91.
  • [19] D. Bhattacherjee and A. Singla, “Network Topology Design at 27,000 km/hour,” in Proc. 15th International Conference on Emerging Networking Experiments and Technologies, Orlando, FL, USA, 2019, pp. 341–354.
  • [20] Z. Guo, J. Liang, N. Xiao, J. Chen, and B. Xie, “Reliability Analysis of Software-Defined Satellite Network Based on Control Delay,” Journal of Physics: Conference Series, vol. 1693(1), article id. 012145, Dec. 2020.
  • [21] M. Handley, “Using Ground Relays for Low-Latency Wide-Area Routing in Mega-constellations,” in Proc. 18th ACM Workshop on Hot Topics in Networks, New York, NY, USA, 2019, pp. 125–132.
  • [22] G. Pan, J. Ye, J. An, and S. Alouini, “Latency versus Reliability in LEO Mega-Constellations: Terrestrial, Aerial, or Space Relay,” in IEEE Transactions on Mobile Computing, doi: 10.1109/TMC.2022.3168081.
  • [23] J. Shi, J. Hu, Y. Yue, X. Xue, W. Liang and Z. Li, "Outage Probability for OTFS Based Downlink LEO Satellite Communication," in IEEE Transactions on Vehicular Technology, vol. 71(3), pp. 3355-3360, March 2022, doi: 10.1109/TVT.2022.3144466.
  • [24] J. Liang, A. U. Chaudhry, J. W. Chinneck, H. Yanikomeroglu, G. K. Kurt,P. Hu, K. Ahmed, S. Martel, “Latency Versus Transmission Power Trade-Off in Free-Space Optical (FSO) Satellite Networks With Multiple Inter-Continental Connections,” in IEEE Open Journal of the Communications Society, vol. PP(99), pp. 1–1, Oct. 2023.
  • [25] S. Arnon, “Performance of a μ\muSatellite Network with an Optical Preamplifier,” Journal of Optical Society of America, vol. 22(4), pp. 708–715, Apr. 2005.
  • [26] A. Polishuk and S. Arnon, “Optimization of a Laser Satellite Communication System with an Optical Preamplifier,” Journal of Optical Society of America, vol. 21(7), pp. 1307–1315, Jul. 2004.
  • [27] H. Lim, J. Park, M. Choi, C. Choi, J. Choi, and J. Kim, “Performance Analysis of DPSK Optical Communication for LEO-to-Ground Relay Link Via a GEO Satellite,” Journal of Astronomy Space Sciences, vol. 37(1), pp. 11–18, Mar. 2020.
  • [28] R. Barrios and F. Dios, “Exponentiated Weibull distribution family under aperture averaging for Gaussian beam waves,” Optical Express, vol. 20(12), pp. 13 055–13 064, May 2012.
  • [29] R. A. Barrios, “Exponentiated Weibull fading channel model in free-space optical communications under atmospheric turbulence,” Ph.D. dissertation, Universitat Politècnica de Catalunya (UPC), May 2013.
  • [30] L. Andrews and R. Phillips, Laser Beam Propagation Through Random Media, ser. SPIE Press monograph. SPIE Press, 2005.
  • [31] ITU-R, “Propagation data required for the design of Earth-space systems operating between 20 THz and 375 THz,” International Telecommunication Union, Recommendation P.1621-1.
  • [32] S. Cakaj, B. Kamo, V. Koliçi, and O. Shurdi, “The Range and Horizon Plane Simulation for Ground Stations of Low Earth Orbiting (LEO) Satellites,” International Journal of Communication, Network and System Science, vol. 4(9), pp. 585–589, Sep. 2011.
  • [33] ITU-R, “Prediction Methods Required for the Design of Earth-Space Systems Operating between 20 THz and 375 THz,” International Telecommunication Union Radiocommunication Recommendation, P.1622, Apr. 2003, [Online]. Available: https://www.itu.int/dms_pubrec/itu-r/rec/p/R-REC-P.1622-1-202208-I!!PDF-E.pdf, Accessed: December 6, 2023.
  • [34] Z. Ghassemlooy, W. Popoola, and S. Rajbhandari, Optical Wireless Communications System and Channel Modelling with Matlab®, CRC Press, Boca Raton, 2019.
  • [35] R.J. MacKay and R.W. Oldford, “Scientific Method, Statistical Method and the Speed of Light,” Statistical Science, vol. 15(3), pp. 254–278, Aug. 2000.
  • [36] A. Roy, J. Pachuau, and A. Saha, “An Overview of Queuing Delay and Various Delay Based Algorithms in Networks,” Computing, vol. 103(10), pp. 2361-2399, Oct. 2021.
  • [37] D. Freedman, T. Marian, J. Lee, K. Birman, H. Weatherspoon, and C. Xu, “Exact Temporal Characterization of 10 Gbps Optical Wide-Area Network,” in Proc. 10th ACM SIGCOMM Conference on Internet Measurement, New York, NY, USA, 2010, pp. 342–355.
  • [38] J. Liang, A.U. Chaudhry, and H. Yanikomeroglu. “Phasing Parameter Analysis for Satellite Collision Avoidance in Starlink and Kuiper Constellations,” in Proc. 2021 IEEE 4th 5G World Forum, Montreal, Canada, 2021, pp. 493–498.
  • [39] SpaceX FCC update, Nov. 2016, “SpaceX Non-Geostationary Satellite System, Attachment A, Technical Information to Supplement Schedule S,” [Online]. Available: https://licensing.fcc.gov/myibfs/download.do?attachment_key=1158350, Accessed: December 6, 2023.
  • [40] O. Montenbruck and E. Gill, Satellite Orbits: Models, Methods, and Applications, Springer-Verlag, Berlin, 2000.
  • [41] E.W. Dijkstra, “A Note on Two Problems in Connexion with Graphs,” Numerische Mathematik, vol. 1, pp. 269–271, Dec. 1959.
  • [42] R. Fields, D. Kozlowski, H. Yura, R. Wong, J. Wicker, C. Lunde, and F. Heine, “5.625 Gbps Bidirectional Laser Communications Measurements Between the NFIRE Satellite and an Optical Ground Station,” in Proc. 2011 International Conference on Space Optical Systems and Applications, Santa Monica, CA, USA, 2011, pp. 44–53.
  • [43] M. Gregory, F. Heine, H. Kämpfner, R. Meyer, R. Fields, and C. Lunde, “TESAT Laser Communication Terminal Performance Results on 5.6 Gbit Coherent Inter Satellite and Satellite to Ground Links,” in Proc. Society of Photo-Optical Instrumentation Engineers, vol. 10565, article id. 105651F, 2017, pp. 324–329.
  • [44] AGI, “Systems Tool Kit (STK),” [Online]. Available: https://www.agi.com/products/stk, Accessed: December 6, 2023.
  • [45] NetworkX Developers, “NetworkX – Network Analysis in Python,” [Online]. Available: https://networkx.org/, Accessed: December 6, 2023.
  • [46] C. Carrizo, M. Knapek, J. Horwath, D. Gonzalez, and P. Cornwell, “Optical Inter-Satellite Link Terminals for Next Generation Satellite Constellations,” in Proc. Society of Photo-Optical Instrumentation Engineers, vol. 11272, article id. 1127203, 2020, pp. 1–11.
  • [47] B. Patnaik and P.K. Sahu, “Inter-Satellite Optical Wireless Communication System Design and Simulation,” IET Communication, vol. 6(16), pp. 2561–2567, Nov. 2012.
  • [48] A. Maho et al., “Assessment of the Effective Performance of DPSK vs. OOK in Satellite-Based Optical Communications,” in Proc. Society of Photo-Optical Instrumentation Engineers, vol. 11180, article id. 111805T, 2019, pp. 1–6.
  • [49] E. Erdogan, I. Altunbas, G. Kurt, M. Bellemare, G. Lamontagne, and H. Yanikomeroglu, “Site Diversity in Downlink Optical Satellite Networks Through Ground Station Selection,” IEEE Access, vol. 9, pp. 31179–31190, Feb. 2021.
  • [50] H. Guo, B. Luo, Y. Ren, S. Zhao, and A. Dang, “Influence of Beam Wander on Uplink of Ground-to-Satellite Laser Communication and Optimization for Transmitter Beam Radius,” Optics Letters, vol. 35(12), pp.1977–1979, Jun. 2010.

Appendix A Dijkstra’s algorithm

The Dijkstra’s algorithm we discissed in section III are also described in the pseudocode in Algorithm 1, where Y is the set for the satellites and ground stations; dist(xi\textit{x}_{i}) is the distance for a node xi\textit{x}_{i} (i.e., a satellite or gd\textit{g}_{d}) to the source ground station gs\textit{g}_{s}; and dq,pi\textit{d}_{q,p_{i}} is the distance of the laser link between nodes q and pi\textit{p}_{i}.

Algorithm 1 Dijkstra’s algorithm to find a shortest path over an FSOSN.

Inputs:
\bullet Adjacency matrix containing propagation distances of laser links between satellites and between satellites and ground stations
\bullet Y, i.e., the set for the satellites and ground stations
Output:
\bullet Shortest path from gs\textit{g}_{s} to gd\textit{g}_{d}
BEGIN

1:S = {}
2:X = Y
3:D(xi\textit{x}_{i}) = \infty for each node xi\textit{x}_{i} in X
4:D(gs\textit{g}_{s}) = 0
5:move gs\textit{g}_{s} from X to S
6:While gd\textit{g}_{d} \notin S
7:    pick q in X that is closest to gs\textit{g}_{s}
8:    move q from X to S
9:    update D(q) = dist(q)
10:    For each pi\textit{p}_{i} adjacent to q in X
11:        If dist(pi\textit{p}_{i}) > dist(q) + dq,pi\textit{d}_{q,p_{i}}
12:            dist(pi\textit{p}_{i}) = dist(q) + dq,pi\textit{d}_{q,p_{i}}
13:        End If
14:    End For
15:End While
16:Output S, i.e., the set containing the shortest path from gs\textit{g}_{s} to gd\textit{g}_{d}

END

{IEEEbiography}

[[Uncaptioned image]]Jintao Liang  (Member, IEEE) received his Bachelor of Engineering degree in Automation Engineering from Shandong University (China) in 2019 and his M.A.Sc degree in Electrical and Computer Engineering from Carleton University (Canada) in 2022.

He is currently a Ph.D. student in the Department of Systems and Computer Engineering in Carleton University (Canada) since 2023.

Mr. Liang’s research interest includes free-space optical satellite networks.

{IEEEbiography}

[[Uncaptioned image]]Aizaz U. Chaudhry  (Senior Member, IEEE) received his B.Sc. degree in Electrical Engineering from the University of Engineering and Technology Lahore (Pakistan) in 1999, and his M.A.Sc. and Ph.D. degrees in Electrical and Computer Engineering from Carleton University (Canada) in 2010 and 2015, respectively.

He is currently a Senior Research Associate with the Department of Systems and Computer Engineering at Carleton University. Previously, he worked as an NSERC Postdoctoral Research Fellow at Communications Research Centre Canada, Ottawa. His research work has been published in refereed venues, and has received several citations. He has authored and co-authored more than thirty-five publications. His research interests include the application of machine learning and optimization in wireless networks.

Dr. Chaudhry is a licensed Professional Engineer in the Province of Ontario, a Senior Member of IEEE, and a Member of IEEE ComSoc’s Technical Committee on Satellite and Space Communications. He serves as a technical reviewer for conferences and journals on a regular basis. He has also served as a TPC Member for conferences, such as IEEE ICC 2021 Workshop 6GSatComNet, IEEE ICC 2022 Workshop 6GSatComNet, IEEE VTC 2022-Spring, IEEE VTC 2022-Fall, the twenty-first international conference on networks (ICN 2022), eleventh advanced satellite multimedia systems conference (ASMS 2022), seventeenth signal processing for space communications workshop (SPSC 2022), and the first international symposium on satellite communication systems and services (SCSS 2022).

{IEEEbiography}

[[Uncaptioned image]]Eylem Erdogan  (Senior Member, IEEE) received B.Sc. and M.Sc. degree from Işık University, Istanbul, Turkey and the Ph.D. degree from Kadir Has University, Istanbul, Turkey in 2014 all in electronics engineering.

He is currently an Associate Professor in the Department of Electrical and Electronics Engineering, Istanbul Medeniyet University. He was a Post-Doctoral Fellow in Electrical Engineering department, Lakehead University, Thunder Bay, ON, Canada from March 2015 to September 2016 and a visiting professor in Carleton University, Ottawa, Canada during summer 2019. He has authored or coauthored more than 40+ papers in peer reviewed SCI/SCI-E journals and international conferences.

His research interests are in the broad areas of wireless communications, including signal processing for wireless communications, the performance analysis of cooperative relaying in cognitive radio networks, unmanned aerial vehicle communications and networks and free space optical communications. He is a Senior Member of IEEE.

{IEEEbiography}

[[Uncaptioned image]]Halim Yanikomeroglu  (Fellow, IEEE) received his Ph.D. degree in Electrical and Computer Engineering from the University of Toronto (Canada) in 1998.

He is currently a Full Professor at Carleton University. He contributed to 4G/5G technologies and non-terrestrial networks. His industrial collaborations have resulted in 39 granted patents. He has supervised or hosted in his lab a total of 140 postgraduate researchers. He has co-authored IEEE papers with faculty members in 80+ universities in 25 countries.

Dr. Yanikomeroglu is also a Fellow of IEEE, Engineering Institute of Canada, and Canadian Academy of Engineering, and an IEEE Distinguished Speaker for ComSoc and VTS. He is currently chairing the WCNC Steering Committee, and he is a Member of PIMRC Steering Committee and ComSoc Emerging Technologies Committee. He served as the General Chair of two VTCs and TP Chair of three WCNCs. He chaired ComSoc’s Technical Committee on Personal Communications. He received several awards, including ComSoc Wireless Communications TC Recognition Award (2018), VTS Stuart Meyer Memorial Award (2020), and ComSoc Fred W. Ellersick Prize (2021).

{IEEEbiography}

[[Uncaptioned image]]Gunes Karabulut Kurt  (Senior Member, IEEE) received the B.S. degree with high honors in electronics and electrical engineering from Bogazici University, Istanbul, Turkey, in 2000 and the M.A.Sc. and the Ph.D. degrees in electrical engineering from the University of Ottawa, ON, Canada, in 2002 and 2006, respectively.

From 2000 to 2005, she was a Research Assistant with the CASP Group, University of Ottawa. Between 2005 and 2006, she was with TenXc Wireless, Canada. From 2006 to 2008, Dr. Karabulut Kurt was with Edgewater Computer Systems Inc., Canada. From 2008 to 2010, she was with Turkcell Research and Development Applied Research and Technology, Istanbul. Between 2010 and 2021, she was with Istanbul Technical University. She is currently an Associate Professor of Electrical Engineering at Polytechnique Montréal, Montréal, QC, Canada. She is a Marie Curie Fellow and has received the Turkish Academy of Sciences Outstanding Young Scientist (TÜBA-GEBIP) Award in 2019. In addition, she is an adjunct research professor at Carleton University.

She is currently serving as an associate technical editor of the IEEE Communications Magazine, an associate editor of IEEE Communication Letters, an associate editor of IEEE Wireless Communications Letters, and an area editor of IEEE Transactions on Machine Learning in Communications and Networking. She is a member of the IEEE WCNC Steering Board. She is serving as the secretary of IEEE Satellite and Space Communications Technical Committee and also the chair of the IEEE special interest group entitled “Satellite Mega-constellations: Communications and Networking”. She is a Distinguished Lecturer of Vehicular Technology Society Class of 2022.

{IEEEbiography}

[[Uncaptioned image]]Peng Hu  (Senior Member, IEEE) received his Ph.D. degree in Electrical Engineering from Queen’s University, Canada.

He is currently a Research Officer at the National Research Council Canada and an Adjunct Professor at the Cheriton School of Computer Science at the University of Waterloo. He has served as an Associate Editor of the IEEE Canadian Journal of Electrical and Computer Engineering, a voting member of the IEEE Sensors Standards committee, and on the organizing/technical committees of industry consortia and international conferences/workshops at IEEE ICC’23, IEEE PIMRC’17, IEEE AINA’15, etc.

His current research interests include satellite-terrestrial integrated networks, autonomous networking, and industrial Internet of Things systems.

{IEEEbiography}

[[Uncaptioned image]]Khaled Ahmed  received his B.S. and MSc degree on wireless communications from Cairo university, Egypt, in 2010 and 2015. He did a PhD and a postdoc in McMaster university, Canada, in 2019 and 2021. His PhD degree was on optical communications while his postdoc was on machine learning for resource-limited devices.

He is currently working as a Member Technical Staff II for MDA, Sainte-Anne-de Bellevue, QC, Canada.

His current research interests include RF and optical communications for terrestrial and satellite networks, optical beamforming networks, and applied machine learning and deep learning.

{IEEEbiography}

[[Uncaptioned image]]Stéphane Martel  received his Bachelor’s degree in Electrical Engineering from McGill University in 2001.

He worked until 2021 in the broadcast television industry where he occupied multiple positions in software and management. He specializes in network communications and embedded software development.

He joined MDA Satellite systems in 2021 as R&D product development manager where he contributes to the advancement of space communication technologies.