Abstract
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 files during off-peak hours and push them during on-peak hours. However, for the enormous number of end devices in ISTNs, the long-time waiting for off/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. Specifically, data transmission of this model consists of two phases. For the first phase, named the file-push-and-placement (FPAP), the satellite employs the NOMA protocol to send information to both the relay and first-tier user. While for the second phase, which is called the file-push-and-delivery (FPAD), second-tier users are served by the relay employing NOMA. Performance analysis of the proposed configuration is carried out, focusing on the exact and asymptotic outage probability (OP), diversity order, and hit probability. Moreover, the influence of key factors on the system performance is also investigated. Compared with the traditional configuration, it is shown analytically and numerically that the proposed scheme achieves lower OP and higher diversity order.
1. 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 superimpose 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–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 NOMA-assisted 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 cache-enabled 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 cache-enabled 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-push-and-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’ closed-form 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
2. 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 (), one active user (), one cache-enabled relay (), and two proactive users ( ) in the considered model. All nodes are equipped with a single antenna and operate in half-duplex mode. Moreover, and in the first tier, as well as and 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, pushes the active files to , and at the same time, pushes the proactive popular contents that will be acquired by to . In the FPAD phase, and receive their corresponding files from 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 denote the FPAP and FPAD as the first and second phases, respectively.
2.1. 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 , and proactive files are collected in another one , where and 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 is given by [17] where 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].
2.2. Signal Transmission
2.2.1. The First Phase (FPAP Phase)
During this phase, transmits the superimposed active signals to and . After channel propagation, the received signals at or are where and denotes the channel coefficient and composite fading distribution describing the links from to or , respectively. The satellite transmit power is , and and are the power control coefficients with constraints and . and are the information-bearing symbols intended for and , respectively. Without loss of generality, it is assumed that . Moreover, refers to the additive noise at the receiver, where is the noise variance.
According to the principle of NOMA, decodes first by treating as interference; thus, the signal-to-interfere-plus-noise ratio (SINR) of decoding at is given by where the signal-to-noise ratio (SNR) is defined as .
As for , it first decodes and then removes to decode . Hence, the SINR of decoding and at can be separately expressed by
2.2.2. The Second Phase (FPAD Phase)
In this phase, and are served by with the NOMA protocol. In the meanwhile, continuously sends data symbols to . However, the signal reception at is subject to the interference from . Assume that the file needed by has already been cached in , this can be achieved if possesses storage capacity. After the channel propagation, the received signals of and can be where is intended for with , indicates the channel coefficient from to or , and we assume that the channel condition of is better than that of without the loss of generality. denotes the transmit power of , and and are the power control coefficients, where and . At last, follows the same distribution as .
According to the NOMA principle, first decodes by treating as interference and then removes to decode . As for , it decodes directly. Hence, the SINR of decoding and at can be expressed as where . Besides, the SINR of decoding at is
As for , the signal reception consists of desired signal from satellite and interference from relay ; thus, its received signal is written as where is the terrestrial channel link from to , which suffers more severe fading than . is the transmit signal from in the second phase, with . In addition, is the transmit power coefficient, and follows the same distribution as . Hence, decodes 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 denotes the speed of light, is the carrier frequency, is the distance between satellite and the first-tier user , is the Boltzman constant, is the receiver noise temperature, and is the carrier bandwidth. Besides, denotes the antenna gain at user , whereas gives the satellite beam gain based on both the satellite beam pattern and position of user and is approximated by [19] With , represents the Bessel function of order , while denotes the maximal beam gain. is the angle between user and beam center with respect to the satellite.
In (2), 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 () obey Shadowed-Rician fading, and the probability density function (PDF) is given by [20] where , , and , is the average power of the multipath, is the line-of-sight (LoS) component, denotes the Nakagami-m fading parameter, and represents the confluent hypergeometric function [21]. According to [21], the cumulative distribution function (CDF) of the satellite link gain is given by where is the Pochhammer symbol defined as .
When it comes to the terrestrial link, which follows Nakagami- fading, the PDF of the channel gain is represented by where is the integer fading severity parameter, 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 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 to UE1 and the relay directly with NOMA protocol. In the second phase, the relay retransmits the decoded signal (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., serves and in the first slot, and serves two users at the second slot, at the same time, serves with newly signal, can be described as follows. In the first slot, transmits the superimposed signal to and . The received signal can be where the power coefficients satisfy , , and . Thus, the SINR of to decode is
Besides, the SINR of decoding , , and at is
Meanwhile, the second time slot is the same as that of the proposed scheme; thus, the description is omitted for brevity.
3. Performance Analysis
3.1. Outage Performance
The OP is defined as the probability that the instantaneous SINR is less than a predefined threshold , i.e., where with being the required data rate, and 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.).
3.1.1. OP of
In our model, 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 is failed, i.e., where . is the required SINR of decoding , which satisfies .
In the second phase, the OP is expressed as
According to (21), we can obtain as follows: where and denote the first and second integrals, respectively. Afterwards, it is easy to find that , and then is expanded as
According to [21], one can arrive at
Then, the close-formed expression of is shown as where
Hence, the derivation of the OP of concludes.
3.1.2. OP of
As for , outage occurs when the detection of is failed, or the detection of fails provided the successful detection of . Therefore, the resultant OP can be derived by where
Besides, and are the corresponding SINRs of decoding and , which satisfies .
3.1.3. OP of
The OP of can be derived straightforward as
3.2. Diversity Order
Due to their complex forms, it is difficult to have a deep understanding of the closed-form OP expressions derived above. On the other hand, the diversity order, defined as , is another key parameter to measure the system performance. To facilitate our analysis of diversity order, the asymptotic outage behaviour of OP at high SNR regimes should be investigated. First, the approximation of CDFs of terrestrial and satellite channels is addressed. For terrestrial channels, by expressing the exponential function in (15) in terms of the Maclaurin series, the PDF is approximated 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
According to (20) and (34), the asymptotic OP of in the first phase is
By assuming , it is derived that the diversity order of is 1.
In the second phase, assuming and , then is approximated as . Therefore, the asymptotic OP is calculated by
According to [21], the OP of in the second phase approximates to
Hence, the diversity order of in the second phase is 0.
3.2.2. Diversity Order of
According to (27) and (31), the asymptotic OP of is shown as
After straightforward mathematical manipulations, the diversity order is obtained as .
3.2.3. Diversity Order of
First, the asymptotic OP of is
At last, one can obtain that the diversity order is .
3.3. Hit Probability
When a user requests a certain file cached in , 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 [15]. The hit probability can be expressed as where denotes the OP of to decode file . In our considered system model, the OP of to decode is derived as where . By substituting (41) into (40), the hit probability of is acquired.
3.4. Comparison to the Scheme without Caching
For comparison, this subsection gives the outage performance of users in the traditional CDRT configuration.
3.4.1. OP of
In the first time slot, the OP of detecting is expressed as
While for the second slot, outage occurs in three cases. First, when fails to decode or given the success decoding of ; thus, the SINR of detecting at is denoted by . Second, detection of and is successful at , but decoding encounters with failure, and then the SINR is . At last, all the signals are decoded correctly in the first slot; hence, is equal to . To sum up, the OP can be derived as where , , and denote the OP of decoding , , and at , respectively. Besides, , , and indicate the OP of decoding in three cases described above, respectively. Those outage probabilities can be calculated by
Moreover, by using the same approach in (25), is obtained as
With (43) and (44), is finally obtained.
3.4.2. OP of
As for , the outage event occurs in situations as follows. First, the decoding of at fails. Second, the decoding of fails given the success decoding of . Third, under the assumption of successfully detecting and at , the decoding of at is failed. Finally, the detection at is successful, but the decoding of or fails at . Therefore, the OP can be derived as where and denote the OP of decoding and at , respectively. The detailed expressions are
3.4.3. OP of
As for , an outage first occurs if the decoding of at fails. Then, the decoding of fails given the successful decoding of . Moreover, the decoding of is failed when the detection of both and succeeds. At last, the detection at succeeds, but fails to detect . Hence, the overall OP is where denotes the OP of detecting at , which is
4. 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 () are shown in Table 1. In addition, we consider different Shadowed-Rician fading conditions for satellite links (), i.e., frequent heavy shadowing (FHS), average shadowing (AS), and infrequent light shadowing (ILS) [17]. The sets of channel coefficients are , , and for FHS, AS, and ILS, respectively. Moreover, and are both modeled as ILS. Besides, . Furthermore, the transmit power of is twice that of , namely, . A fixed power allocation scheme is applied that the power coefficients are set by , and . Moreover, the content popularity parameter is , and the number of files is . At last, the transmission rate requirement for , , and 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 , 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 . The reason behind this is that the decoding at fails more frequently with small SNRs (), causing less interference to ; thus, the OP will drop when SNR increases. On the other hand, when the SNR becomes large, the decoding at encounters much less failure, and then the OP will increase because the interference at is more severe. For comparison, it is clear that the OP of proactive users ( and ) 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 , , and in the scheme without caching.


To evaluate the diversity order, we concentrate on OP values at and . Specifically, the comparison of the diversity order between the two schemes is shown in Table 2. Note that the diversity order of shows no difference in two schemes, which is not the case when considering and . In particular, if and , then the diversity orders achieved are the same for two schemes. While for larger and 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 or for the proposed scheme, while it is impacted by all the channel links for the scheme without caching.
Figure 5 illustrates the influence of channel parameters on the OP of . In phase one, the OP is affected by , which is set as ILS or AS. It is clearly observed that the OP reduces when changes from AS to ILS, since the channel condition becomes better. In contrast, the OP of in the second phase is decided by both and . The channel quality concerning is tuned by modifying parameter . It is demonstrated that the OP almost keeps unchanged when varies from 1 to 2, indicating that has little impact on the OP of . Furthermore, the impacts of on and that of on 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 , where different content popularity parameters are included. The number of content files 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 in phase one, where the CSI uncertainty is modeled as a Gaussian noise independent of the actual channel 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.

5. 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.
Acknowledgments
This study is supported by the Chongqing Key R&D Project (cstc2021jscx-gksbX0057) and National Mountain Highway Engineering Technology Research Center (GSGZJ-2022-08).