NOMA-Based Integrated Satellite-Terrestrial Networks with Wireless Caching

To decrease the transmission delay and alleviate communication congestion, caching has been applied in the integrated satellite-terrestrial networks (ISTNs). In traditional caching schemes, the cache-enabled relays update ﬁ les during o ﬀ -peak hours and push them during on-peak hours. However, for the enormous number of end devices in ISTNs, the long-time waiting for o ﬀ /on-peak hours would decrease the quality of service. To address this issue, this paper proposes a two-tier nonorthogonal multiple access-(NOMA-) based ISTN with wireless caching. Speci ﬁ cally, data transmission of this model consists of two phases. For the ﬁ rst phase, named the ﬁ le-push-and-placement (FPAP), the satellite employs the NOMA protocol to send information to both the relay and ﬁ rst-tier user. While for the second phase, which is called the ﬁ le-push-and-delivery (FPAD), second-tier users are served by the relay employing NOMA. Performance analysis of the proposed con ﬁ guration is carried out, focusing on the exact and asymptotic outage probability (OP), diversity order, and hit probability. Moreover, the in ﬂ uence of key factors on the system performance is also investigated. Compared with the traditional con ﬁ guration, it is shown analytically and numerically that the proposed scheme achieves lower OP and higher diversity order.


Introduction
Future networks are expected to provide diverse services to cope with the ever-increasing traffic demands of various services. Nevertheless, limited by the network capacity and coverage, only depending on the terrestrial communication systems cannot provide wireless access with high data rate and reliability at every place on the earth, especially in environmentally harsh areas like oceans and mountains [1]. Hence, integrated satellite-terrestrial networks (ISTNs) are deemed to be new network architectures to accommodate diverse services and applications with different quality of service (QoS) requirements [2].
Meanwhile, as the latest member of the multiple access family, nonorthogonal multiple access (NOMA) is regarded as a promising technique in the next-generation wireless communication system [3]. In this paper, we focus on the power-domain NOMA, of which the key idea is to superim-pose multiple signals in the power domain at the transmitter and employ successive interference cancellation (SIC) at the receiver [4]. With NOMA scheme, more users can be served simultaneously; hence, the spectrum efficiency is improved as a result. Therefore, the application of NOMA in ISTNs has attracted great attention in academia [5][6][7][8].
In future wireless networks, not only the spectral efficiency but also the transmission delay should be paid special attention to. Stemming from that, the idea of caching is introduced. Instead of retrieving information from a central server, users can ask for the cache-enabled server that has replicated popular information [9]. As a result, the response time required to fetch a content file can be reduced. The authors in [10] analyzed and optimized the outage performance of a multiple cache-enabled amplify-and-forward relay network. Besides, caching can also help to further enhance the spectral efficiency of NOMA [11]. Cache content placement optimization in NOMA networks was investigated in [12,13] which studied the corresponding coding mode. Furthermore, the authors in [14] proposed a NOMA-based multicast scheme that pushing and multicasting content objects can be accomplished at once, which boosts the spectrum efficiency. On this basis, two NOMAassisted caching strategies were developed, namely, the push-then-deliver and the push-and-deliver [15]. However, these two strategies were studied individually, and the importance to combine them in practice was not investigated.
Recently, the significance of applying caching to the ISTN has become a consensus in the academic community. The authors in [16] compared the outage probability (OP) of a relay-assisted ISTN with two representative cache placement schemes. As for the NOMA-assisted caching schemes applied in ISTNs, [17] investigated the OP of the cacheenabled relays and the hit probability of users in the NOMA-based hybrid satellite-terrestrial content delivery network. Moreover, a satellite-aerial-terrestrial network with cache-enabled aerial relays was introduced in [18], where the NOMA scheme was implemented to deliver content and push other currently most popular content to cacheenabled aerial relays simultaneously.
The combination of NOMA and ISTNs can solve the problem of multiple nodes, but in addition to the requirements of massive nodes, future mobile communications also require low latency data transmission; so, cache is introduced to reduce latency. The combination of NOMA technology and cache technology can reduce the delay and increase spectral efficiency. However, the existing literature rarely considers how to make better use of the system spectrum resources and time resources in an environment with complex channel conditions, where the server needs to push files to users and place files to cache nodes at the same time. Therefore, this paper proposes a NOMA-based ISTN with wireless caching. Specifically, the active users request the active contents and are served by the satellite directly. In contrast, proactive users request proactive contents. However, due to heavy shadowing and masking effects, it is highly probable that proactive users lack direct links to the satellite and thus need help from cache-enabled relays. By exploiting the NOMA protocol, the content file placement of relays, and the files pushed to the active users from the satellite can be accomplished simultaneously.
The main contributions are summarized as follows: (a) The overall communication includes the file-pushand-placement (FPAP) and file-push-and-delivery (FPAD) phases. Specifically, in the FPAP phase, active and proactive files are sent to the active users and cache-enabled relays from the satellite by using the NOMA scheme. In the FPAD phase, the satellite pushes new files to active users, while the relays deliver files to the proactive users by using the NOMA technique (b) The performance of the proposed system is thoroughly analyzed, with emphasis on users' closedform OP expressions. Then, the diversity gain is derived from the asymptotic behavior of OP. In addition, the hit probability of the relay node is studied. At last, the influence of key system parameters on outage performance is investigated (c) A comparison between the proposed scheme and the NOMA-based ISTN without caching is carried out, where the result demonstrates that the system performance is improved by our scheme

System Model
The proposed two-tier heterogeneous cache-enabled NOMA-based ISTN is depicted in Figure 1. Due to the complexity and diversity of channel conditions, users with stable and reliable direct links with satellites are divided into the first tier and defined as active users. Users suffering from severe shadow effect without a direct link are divided into the second tier and defined as proactive users. It is worth noting that the satellite can directly serve the active users, but it needs to use the relay node with good channel conditions to assist the proactive users. Due to the long backhaul link of the satellite communication network, to avoid the high delay of proactive users, the cache-enabled relay nodes are used to help the communication of proactive users. Assuming that there is good channel state quality between the satellite and the cache-enabled relay node, the relay node is divided into the first tier. In addition, content files are divided into active files and proactive files according to different requirements, i.e., the active file is defined as the file currently requested by the user, and the proactive file belongs to the file that the user does not request at present but will request in the future. Therefore, when the satellite directly serves the first-tier users, the pushed files are active; when the satellite places the content file to the cache-enabled relay node, the pushed file is the proactive. Above all, select active users and cache-enabled relays as the first-tier nodes and proactive users as the second-tier nodes. As for the signal transmission, the satellite communicates with the nodes in the first tier directly, while with the nodes in the second tier with the help of relays. To facilitate the analysis, it is supposed that there is one satellite (S), one active user (U 1 ), one cache-enabled relay (R), and two proactive users (U 2 ,U 3 ) in the considered model. All nodes are equipped with a single antenna and operate in half-duplex mode. Moreover, U 1 and R in the first tier, as well as U 2 and U 3 in the second tier, form NOMA clusters. It is worth noting that two-user NOMA-based downlink transmission has been proposed for long-term evolution (LTE) advanced [15].
As aforementioned, the overall communication consists of two phases. In the FPAP phase, by using the NOMA scheme, S pushes the active files to U 1 , and at the same time, pushes the proactive popular contents that will be acquired by U 2 to R. In the FPAD phase, U 2 and U 3 receive their corresponding files from R using the NOMA protocol. Without loss of generality, we assume that a NOMA cluster consists of two users, where the user with good channel condition is denoted by the strong user, while the other one is the weak user. For ease of description, in the following analysis, we 2 Wireless Communications and Mobile Computing denote the FPAP and FPAD as the first and second phases, respectively.

Caching Model.
In the light of [14], the content files are divided into two categories: active files which are requested by the active users currently and proactive files not requested by the proactive users now but will be requested soon. Correspondingly, the relay caches the proactive files in the first phase, and then it can provide service to the proactive users in the second phase. Besides, the active files are pushed to active users directly in the first phase. The active files are collected in a finite content catalog denoted by F a = f f a 1 ,⋯,f a N g, and proactive files are collected in another one F p = f f p 1 ,⋯,f p M g, where N and M are the numbers of active and proactive files, respectively. Besides, the popularity of the requested files obeys a Zipf distribution, e.g., the popularity of f p m ð1 ≤ m ≤ MÞ is given by [17] where ξ > 0 denotes the shape parameter defining the content popularity skewness. The popularity of the file becomes more concentrated when ξ gets larger, while the cache utilization ratio depends on the number of files requested by users [10].

Signal
where h i and ffiffiffiffiffi Q i p denotes the channel coefficient and composite fading distribution describing the links from S to U 1 or R, respectively. The satellite transmit power is P s , and ω 1 and ω 2 are the power control coefficients with constraints ω 1 + ω 2 = 1 and ω 1 > ω 2 . x 1 and x 2 are the information-bearing symbols intended for U 1 and R, respectively. Without loss of generality, it is assumed that

Wireless Communications and Mobile Computing
Efjx 1 j 2 g = Efjx 2 j 2 g = 1. Moreover, n 1 ∼ CN ð0, N 0 Þ refers to the additive noise at the receiver, where N 0 is the noise variance.
According to the principle of NOMA, U 1 decodes x 1 first by treating x 2 as interference; thus, the signal-to-interfereplus-noise ratio (SINR) of decoding x 1 at U 1 is given by where the signal-to-noise ratio (SNR) is defined as ρ s = P s /N 0 . As for R, it first decodes x 1 and then removes x 1 to decode x 2 . Hence, the SINR of decoding x 1 and x 2 at R can be separately expressed by 2.2.2. The Second Phase (FPAD Phase). In this phase, U 2 and U 3 are served by R with the NOMA protocol. In the meanwhile, S continuously sends data symbols to U 1 . However, the signal reception at U 1 is subject to the interference from R. Assume that the file needed by U 3 has already been cached in R, this can be achieved if R possesses storage capacity. After the channel propagation, the received signals of U 2 and U 3 can be where x 3 is intended for U 3 with Efjx 3 j 2 g = 1, h j indicates the channel coefficient from R to U 2 or U 3 , and we assume that the channel condition of h ru 2 is better than that of h ru 3 without the loss of generality. P r denotes the transmit power of R, and τ 2 and τ 3 are the power control coefficients, where τ 2 + τ 3 = 1 and τ 3 > τ 2 . At last, n 2 follows the same distribution as n 1 . According to the NOMA principle, U 2 first decodes x 3 by treating x 2 as interference and then removes x 3 to decode x 2 . As for U 3 , it decodes x 3 directly. Hence, the SINR of decoding x 3 and x 2 at U 2 can be expressed as where ρ r = P r /N 0 . Besides, the SINR of decoding x 3 at U 3 is As for U 1 , the signal reception consists of desired signal from satellite S and interference from relay R; thus, its received signal is written as where h ru 1 is the terrestrial channel link from R to U 1 , which suffers more severe fading than h su 1 . x 1 ′ is the transmit signal from S in the second phase, with Efjx 1 ′ j 2 g = 1.
In addition, τ 1 is the transmit power coefficient, and n 3 follows the same distribution as n 1 . Hence, U 1 decodes x 1 ′ directly, and the corresponding SINR is expressed by 2.3. Channel Model. For the satellite links, effects such as antenna gain, path loss, and link fading should be taken into account. We consider the GEO satellite; thus, the scaling parameter is given as where c denotes the speed of light, f is the carrier frequency, d i ≈ 35786 km is the distance between satellite S and the first-tier user iði ∈ fsu 1 , srgÞ, K B = 1:38 × 10 −23 J/°K is the Boltzman constant, T is the receiver noise temperature, and B c is the carrier bandwidth. Besides, T i denotes the antenna gain at user i, whereas G i ðφ i Þ gives the satellite beam gain based on both the satellite beam pattern and position of user i and is approximated by [19] With u i = 2:07123ðsin φ i /sin φ i3dB Þ, J v ð·Þ represents the Bessel function of order v, while G max denotes the maximal beam gain. φ i is the angle between user i and beam center with respect to the satellite.
In (2), h i denotes the channel fading vector, which is usually modeled by composite fading distributions to describe more accurately the amplitude fluctuation of the signal envelope. Let the satellite link (jh i j 2 , i ∈ fsu 1 , srg) obey Shadowed-Rician fading, and the probability density function (PDF) is given by [20] where denotes the Nakagami-m fading parameter, and 1 F 1 ða ; b ; cÞ represents the confluent hypergeometric function [21].

Wireless Communications and Mobile Computing
According to [21], the cumulative distribution function (CDF) of the satellite link gain is given by where ð•Þ k is the Pochhammer symbol defined as ðxÞ k = xðx + 1Þðx + k − 1Þ.
When it comes to the terrestrial link, which follows Nakagami-m fading, the PDF of the channel gain is represented by where m j is the integer fading severity parameter, Ω j is the average power of each link, and Γð•Þ is the Gamma function [21]. In addition, the CDF of the terrestrial link gain can be expressed as where γða, xÞ = Ð x 0 e −t t a−1 dt denotes the lower incomplete Gamma function [21].
2.4. Traditional Scheme without Caching. The model highly related to the proposed is the NOMA-based coordinated direct and relay transmission (CDRT) [22], which is depicted in Figure 2. In CDRT, a base station (BS) directly communicates with user equipment 1 (UE1) while communicating with user equipment 2 (UE2) via a relay. The communication consists of two phases, i.e., in the first phase, BS transmits the superposed signal ðα 1 x 1 + α 2 x 2 Þ to UE1 and the relay directly with NOMA protocol. In the second phase, the relay retransmits the decoded signal x 2 (required by UE2) to UE2, while the BS transmits a new data symbol to UE1. The main challenge of nonorthogonal CDRT can be solved by using the inherent property of NOMA that allows a receiver to obtain side information such as other UE's data for interference cancellation. Stemming from that, it is reasonable to apply NOMA-based CDRT to achieve high spectrum resource utilization.
The main difference between proposed model and the traditional CDRT lies in the second phase. Specifically, in the proposed model, relay can serve multiple users simultaneously, while CDRT serves only one. To achieve the same time and spectrum utilization as the proposed scheme, the NOMA-based CDRT without caching, i.e., S serves U 1 and R in the first slot, and R serves two users at the second slot, at the same time, S serves U 1 with newly signal, can be described as follows. In the first slot, S transmits the superimposed signal to R and U 1 . The received signal can be where the power coefficients satisfy ς 1 + ς 2 + ς 3 = 1, ς 1 > ς 2 > ς 3 , and ς 1 > ς 2 + ς 3 . Thus, the SINR of U 1 to decode x 1 is Besides, the SINR of decoding x 1 , x 2 , and x 3 at R is Þh sr Q sr j j 2 ρ s + 1 , Meanwhile, the second time slot is the same as that of the proposed scheme; thus, the description is omitted for brevity.

Performance Analysis
3.1. Outage Performance. The OP is defined as the probability that the instantaneous SINR λ is less than a predefined threshold λ th , i.e., where λ th = 2 2R th − 1 with R th being the required data rate, and F λ ð•Þ is the CDF of λ. We emphasize that the analytical derivation is based on the following assumptions. Firstly, the CSI and SIC are crucial to the outage performance, and there have been several studies concerning this issue [23]. The CSI can be obtained with channel estimation methods, which have been adopted in most the existing works [24,25]. As the focus herein is the caching scheme, we assume that the perfect CSI is available, and the SIC is error-free. Secondly, both the terrestrial and satellite channels are independent and identically distributed (i.i.d.).

Wireless Communications and Mobile Computing
3.1.1. OP of U 1 . In our model, U 1 receives signals within two phases, and then its OP needs to be discussed separately. In the first phase, the outage occurs when the decoding of x 1 is failed, i.e., where ϕ 1 = λ th 1 /ðω 1 − ω 2 λ th 1 Þρ s Q su 1 . λ th 1 is the required SINR of decoding x 1 , which satisfies λ th 1 < ω 1 /ω 2 .
In the second phase, the OP is expressed as According to (21), we can obtain 1 − P ð2Þ u 1 ,out as follows: where J 1 and J 2 denote the first and second integrals, respectively. Afterwards, it is easy to find that J 1 = 1, and then J 2 is expanded as According to [21], one can arrive at Then, the close-formed expression of P ð2Þ u 1 ,out is shown as where Hence, the derivation of the OP of U 1 concludes.
3.1.2. OP of U 2 . As for U 2 , outage occurs when the detection of x 3 is failed, or the detection of x 2 fails provided the successful detection of x 3 . Therefore, the resultant OP can be derived by where Besides, λ th 2 and λ th 3 are the corresponding SINRs of decoding x 2 and x 3 , which satisfies λ th 3 < τ 3 /τ 2 .

OP of U 3 . The OP of U 3 can be derived straightforward as
Thus, the corresponding CDF of (30) is obtained straightforward, shown as As for satellite links, via using the series representation in [21], one can arrive at Thus, the approximated CDF expression for the satellite links can be obtained as follows: And the final approximated result is With (31) and (34), diversity orders of users can be derived as follows.
3.2.1. Diversity Order of U 1 . According to (20) and (34), the asymptotic OP of U 1 in the first phase is By assuming ρ s ⟶ ∞, it is derived that the diversity order of U 1 is 1.
In the second phase, assuming ρ s ⟶ ∞ and ρ r ⟶ ∞, then λ u 1 ,x 1 ′ is approximated as τ 1 jh su 1 j 2 Q su 1 /jh ru 1 j 2 . Therefore, the asymptotic OP is calculated by According to [21], the OP of U 1 in the second phase approximates toP Hence, the diversity order of U 1 in the second phase is 0.

Diversity
Order of U 2 . According to (27) and (31), the asymptotic OP of U 2 is shown as After straightforward mathematical manipulations, the diversity order is obtained as m ru 2 .

Diversity Order of
At last, one can obtain that the diversity order is m ru 3 .

Hit Probability.
When a user requests a certain file cached in R, a cache hit event happens. Hence, an effective criterion to evaluate the performance of content pushing is the hit probability, which indicates the probability that, during the content delivery phase, a user finds its requested file in the associated cache-enabled R [15]. The hit probability can be expressed as where P f p n r,out denotes the OP of R to decode file f p n . In our considered system model, the OP of R to decode f p 2 is derived as where ϕ 4 = λ th 2 /ω 2 ρ s Q sr . By substituting (41) into (40), the hit probability of R is acquired.

Comparison to the
While for the second slot, outage occurs in three cases. First, when R fails to decode x 3 or x 2 given the success decoding of x 1 ; thus, the SINR of detecting x 1 ′ at U 1 is denoted by λ no,1 u 1 ,x 1 ′ = τ 1 jh su 1 j 2 ρ s Q su 1 . Second, detection of x 1 and x 2 is successful at R, but decoding x 3 encounters with failure, and then the SINR is λ no,2 u 1 ,x 1 ′ = τ 1 jh su 1 j 2 ρ s Q su 1 /τ 2 jh ru 1 j 2 ρ r + 1 . At last, all the signals are decoded correctly in the first slot; hence, λ no,3 u 1 ,x 1 ′ is equal to λ u 1 ,x 1 ′ . To sum up, the OP can be derived as where P 1 , P 2 , and P 3 denote the OP of decoding x 1 , x 2 , and x 3 at R, respectively. Besides, P 1 x 1 ′ , P 2 x 1 ′ , and P 3 x 1 ′ indicate the OP of decoding x 1 ′ in three cases described above, respectively. Those outage probabilities can be calculated by Moreover, by using the same approach in (25), P 2 x 1 ′ is obtained as With (43) and (44), P no,2 u 1 ,out is finally obtained.

OP of U 2 .
As for U 2 , the outage event occurs in situations as follows. First, the decoding of x 1 at R fails. Second, the decoding of x 2 fails given the success decoding of x 1 . Third, under the assumption of successfully detecting x 1 and x 2 at R, the decoding of x 3 at R is failed. Finally, the detection at R is successful, but the decoding of x 3 or x 2 fails at U 2 . Therefore, the OP can be derived as where P 4 and P 5 denote the OP of decoding x 2 and x 3 at U 2 , respectively. The detailed expressions are 3.4.3. OP of U 3 . As for U 3 , an outage first occurs if the decoding of x 1 at R fails. Then, the decoding of x 2 fails given the successful decoding of x 1 . Moreover, the decoding of x 3 is failed when the detection of both x 1 and x 2 succeeds. At last, the detection at R succeeds, but U 3 fails to detect x 3 . Hence, the overall OP is where P 6 denotes the OP of detecting x 3 at U 3 , which is

Numerical Results
In this part, numerical simulations are performed to showcase the performance of the proposed cache-enabled NOMA-based ISTNs. First, the satellite link parameters (Q i ) are shown in Table 1. In addition, we consider different Shadowed-Rician fading conditions for satellite links (h i ), i.e., frequent heavy shadowing (FHS), average shadowing (AS), and infrequent light shadowing (ILS) [17]. The sets of channel coefficients ðb i , m i , Ω i Þ are ð0:063,0:739,8:97 × 10 −4 Þ, ð0:126,10:1,0:835Þ, and ð0:158,19:4,1:29Þ for FHS, AS, and ILS, respectively. Moreover, h su 1 and h sr are both modeled as ILS. Besides, m ru 1 = 2, Ω ru 1 = 1, m ru 2 = 2, Ω ru 2 = 1, m ru 2 = 1, Ω ru 2 = 1. Furthermore, the transmit power of S is twice that of R, namely, P s = 2P r . A fixed power allocation scheme is applied that the power coefficients are set by ω 1 = 0:6, ω 2 = 0:4, τ 1 = 1, τ 2 = 0:3, τ 3 = 0:7, and ς 1 = 0:6, ς 2 = 0:3, ς 3 = 0:1. Moreover, the content popularity parameter is ξ = 0:7, and the number of files is M = 10. At last, the transmission rate requirement for x 1 , x 2 , and x 3 is 0.5 (bps/Hz), 0.8 (bps/Hz), and 0.8 (bps/Hz), respectively. Figures 3 and 4 illustrate the outage performance of users for different schemes. First of all, the good match between the analytical and numerical results validates our derivations. As for the outage performance of U 1 , due to the same power control and signal transmission, the OP is the same in the first phase for both the proposed and traditional schemes. While in the second phase, it is shown that the OP for the scheme without caching has a turning point at SNR = 10 dB. The reason behind this is that the decoding at R fails more frequently with small SNRs (<10 dB), causing less interference to U 1 ; thus, the OP will drop when SNR increases. On the other hand, when the SNR becomes large, the decoding at R encounters much less failure, and then the OP will increase because the interference at U 1 is more severe. For comparison, it is clear that the OP of proactive users (U 2 and U 3 ) in our scheme is lower than the scheme without caching. This is because the proactive users are affected by two factors, namely, the channel links and files' popularity, in the proposed scheme, while they are influenced by the channel links of h ru 1 , h ru 2 , and h ru 3 in the scheme without caching.

Wireless Communications and Mobile Computing
To evaluate the diversity order, we concentrate on OP values at SNR = 25 dB and SNR = 30 dB. Specifically, the comparison of the diversity order between the two schemes is shown in Table 2. Note that the diversity order of U 1 shows no difference in two schemes, which is not the case when considering U 2 and U 3 . In particular, if m ru 2 < 1 and m ru 3 < 1, then the diversity orders achieved are the same for two schemes. While for larger m ru 2 and m ru 3 values, the proposed scheme is superior to the counterpart. This observation lies in that the signal transmission of proactive users is subject to the channel quality of h ru 2 or h ru 3 for the proposed scheme, while it is

Wireless Communications and Mobile Computing
impacted by all the channel links for the scheme without caching. Figure 5 illustrates the influence of channel parameters on the OP of U 1 . In phase one, the OP is affected by h su 1 , which is set as ILS or AS. It is clearly observed that the OP reduces when h su 1 changes from AS to ILS, since the channel condition becomes better. In contrast, the OP of U 1 in the second phase is decided by both h su 1 and h ru 1 . The channel quality concerning h ru 1 is tuned by modifying parameter m ru 1 . It is demonstrated that the OP almost keeps unchanged when m ru 1 varies from 1 to 2, indicating that h ru 1 has little impact on the OP of U 1 . Furthermore, the impacts of h ru 2 on U 2 and that of h ru 3 on U 3 are depicted in Figure 6, where phase two is considered. As expected, the outage event occurs more frequently if the channel condition gets worse. Figure 7 depicts the hit probability of R, where different content popularity parameters are included. The number of content files M is set as 10, and ξ varies from 5 to 1, and then to 0.1. A large ξ indicates the request for files with high popularity, whereas a small one corresponds to the requests with heavy-tailed popularity. Hence, a higher hit probability is obtained with a larger ξ, which is consistent with the numerical results shown in Figure 7. It is also worth noting that enhancing the SNR value is beneficial to level up the hit probability.
In the simulations above, we assume that the perfect CSI is available to all nodes. However, in practice, the channel estimation can never be error-free; thus, it is necessary to investigate the impact of uncertainty in CSI on the outage performance. To this end, Figure 8 draws the outage performance of U 1 in phase one, where the CSI uncertainty is modeled as a Gaussian noise independent of the actual chan-nel coefficients. Note that the independence between the actual CSI and estimation error holds if the MMSE channel estimation method is employed. In this figure, the case of zero variance indicates the perfect CSI acquisition. From this figure, the CSI uncertainty harms the outage probability, and the larger the uncertainty is, the higher the outage probability will be.

Conclusion
In this paper, a NOMA-based ISTN with wireless caching is proposed to reduce the transmission delay and improve spectrum efficiency. Specifically, the overall signal transmission includes two phases. In the first phase, the active user and cache-enabled relay are served by the satellite using NOMA protocol. In the second phase, the proactive users are served by the relay, and at the same time, the active user is served by the satellite. The exact-form and asymptotic OP, diversity order, and the hit probability are derived. In addition, the impact of key parameters on the system performance is discussed. Both simulation and analytical results are provided to validate the superior performance of the proposed scheme over the conventional one.

Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
The authors declare that they have no conflicts of interest.