{'fraction_non_alphanumeric': 0.04115100973560286, 'fraction_numerical': 0.016272546107662822, 'mean_word_length': 4.4733916909823535, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
On the Use of NB-IoT over GEO Satellite Systems with Time-Packed Optical Feeder Links for Over-the-Air Firmware/Software Updates of Machine-Type Terminals

The verticals of 5G, such as the automotive, smart grid, and smart cities sectors, will bring new sensors and IoT devices requiring Internet connectivity. Most of these machine-type terminals will be sparsely distributed, covering a very large geographical area and, from time to time, will have to update their software, firmware, and/or other relevant data. Given this situation, one viable solution to implement the “Over-the-Air” update of these IoT terminals can be done with the aid of GEO satellite systems. However, due to the ultra-dense radio frequency reuse factor that contemporary High-Throughput Satellite (HTS) systems implement in the access link to serve the IoT terminals, the use of a time-packed Free Space Optical (FSO) link represents a practical solution to avoid the bottleneck that the satellite gateway experiences in the feeder link. The performance of both Detect-and-Forward and Decode-and-Forward relaying strategies are studied, assuming that the single-carrier M-PAM symbols that are transmitted on the optical feeder link are mapped into M-QAM symbols that modulate the multiple sub-carriers of the OFDM-based radio access link. In addition, the benefits of encapsulating the NB-IoT frames into DVB-S2(X) satellite frames is also analyzed in detail. The effects of the impairments introduced in both the optical feeder and radio access links are characterized in detail, and the end-to-end error correction capabilities of the Modulation and Coding Schemes (MCS) defined in the contemporary releases of the NB-IoT and DVB-S2(X) standards are studied for different working regimes.


Introduction
In the forthcoming years, an increased data rate capacity will be needed to provide enhanced Mobile Broadband (eMBB) and massive Machine-Type Communications (mMTC) fueled, among other reasons, by the large demand of video transmissions and IoT communications that are foreseen in future vertical services [1]. Specifically, according to CISCO, it is expected that Machine-to-Machine (M2M) connections will ammount to 14.7 billion by 2023. Applications such as home automation, home security, video surveillance, connected white goods, and tracking are expected to represent almost half of the total M2M connections by 2023 [2]. Furthermore, the M2M services that are currently experiencing the fastest growth are connected car applications, such as fleet management, in-vehicle entertainment, emergency calling, vehicle diagnostics, and navigation, with a Compound Annual Growing Rate (CAGR) in the order of 30%. All these applications share the same requirement, known as Over-The-Air (OTA) programming, to update the software (e.g., maps for navigation), security keys (e.g., for cryptography in the IoT devices), among other firmware updates.
Lately, 3GPP has started its studies to analyze the feasibility of integrating satellite networks into its mobile communications systems [3]. Therefore, until these studies are completed, OTA services provided over satellites are based on proprietary solutions owned by the satellite operators. Regarding the feeder links, the dominant State-of-Art (SoA) technology for implementing them nowadays is based on multi-beam reflect arrays in the Ka/Ku radio frequency band [4,5], using the DVB-S2 standard. However, there are many satellites in different orbits allocated in the Ka/Ku bands, and limited bandwidth is currently available for new deployments, particularly for 5G networks and beyond [6]. Note that 6G systems will have to support terminals with a mobility speed of up to 1000 km/h, a connectivity density of 10 7 devices/km 2 , and peak data dates in downlink of up to 1 Tera-bit-per-second (Tbps) to name a few of its requirements [7][8][9]. This means that the capacity of a satellite feeder link must be increased, in order to avoid possible bottlenecks. Non-terrestrial networks, particularly satellite ones, will have to play a key role in the next generation of mobile communication systems, integrating airborne, terrestrial, and satellite networks to support in-flight connectivity. Then, multiple technologies, not only radio-based ones, must be used to face the constrains that are foreseen for 6G.
In practical terms, the growth of OTA applications that is expected in the near future implies that the downlink channel (i.e., from the base station to the IoT terminals) will have to be improved to support the forecast data traffic demand. Furthermore, given that machine-type terminals may be sparsely deployed in large geographical areas, then the use of satellite network would be an excellent option to enable reliable IoT connectivity in a global scale. Towards this regard, the advent of Very High Throughput Satellite (VHTS) systems will allow one to achieve a total network capacity of few Tbps [10,11]. In line with this, satellite operators are starting to study the viability of deploying optical feeder links for Beyond 5G applications in GEO [6], where primary research analyses have been already conducted [12]. Recent experimental studies in single optical LEO-to-ground links have also been reported in [13]. Larger communication bandwidths and free access to spectrum are the key points in favor of optical wireless communications. Unfortunately, both radio and optical wireless communications suffer from channel impairments. Therefore, from a longterm perspective, hybrid solutions combining both technologies should be implemented in the forthcoming generations of satellite networks [14] to increase the capacity of the feeder link.
Regarding the potential techniques to increase the spectral efficiency of the optical feeder link, e.g., beamforming [15], NOMA [16], and frequency reuse [17], this paper resorts to the ones based on shrinking the transmitted pulses. The so-called Faster-Than-Nyquist or time-packed techniques. Initially proposed by Mazo in the 70s [18], has been proposed as a a potential technique for increasing the spectral efficiency of Beyond 5G systems [19]. From the satellite point of view, the possibility of reducing the transmission time without augmenting the transmission bandwidth fits very well for overcoming the drawbacks of the optical channel, which heavily depend on atmospheric conditions such as clouds [6]. By doing so, the link layer may adjust the transmission time of the frames without reducing the number of symbols to transmit according to weather forecasts [20]. This paper focuses on using optical feeder links for GEO satellite networks. Since initial studies about the introduction of satellites into the 3GPP landscape have focused on low-latency applications, most of these satellite communication scenarios are dominated by the use of Low Earth Orbit (LEO) satellite constellations [21]. However, due to the longer space-to-ground link distances of GEO satellites, more spectral-efficient techniques should be introduced in their architectures. Furthermore, observing the evolution that 3GPP standards have experienced in the past, first integrating vehicular and railway networks into 5G [8] and then airbone and satellites into 6G [7], it is highly probable that at some point in the future, extra-terrestrial communications will also be integrated into the 3GPP ecosystem. In this latter scenario, GEO satellites are in a good position for relaying data from the Earth to the outer space, as discussed in the Moonlight initiative from ESA [22]. Therefore, this paper studies the approach that GEO satellites should use to forward the information that they receive from the satellite gateway [23].
In the "ideal" case of a GEO satellite that implements a fully regenerative payload, the optical feeder link would be terminated in the satellite, and a robust modulation and coding scheme should be selected to address the bit error bursts that the turbulent optical satellite wireless channel introduce [24]. In the transparent non-regenerative "bent-pipe" solutions, on the other hand, the instantaneous value of the radio signal is used to modulate the intensity of the optical carrier of the feeder link Laser Diode (LD) with the aid of an external Match-Zehnder Modulator (MZM) [25]. Moreover, if time-packing encoding is applied on the real-valued electrical signal that is used to modulate the intensity of the LD beam, the data throughput of the optical feeder link can be increased even further, without the necessity of using a wider communication bandwidth. This effect is obtained by shrinking the separation between adjacent transmitted pulses [26], mitigating part of the Inter-Symbol Interference (ISI) power that time-packing introduces with the aid of a linear equalizer that is placed on-board the GEO satellite before symbol detection [27]. It is important to highlight that the impact of the residual ISI, which remains in the forward link after the GEO satellite relaying, can be further mitigated with the proper selection of the Modulation and Coding Schemes (MCS) to communicate with the NB-IoT terminal [28].
The forward link of a satellite system can be divided into two parts, namely: (i) The optical feeder link (uplink), from the ground station to the GEO satellite, and (ii) the radio access link (downlink), from the GEO satellite to the NB-IoT user terminals. Therefore, and in order to improve the achievable throughput of the forward direction of communication of the GEO satellite system under different working conditions, this paper studies three different relaying architectures, namely: (1) Detect-and-Forward (a non-regenerative strategy), where the GEO satellite only detects the symbols of the NB-IoT frames that modulate the intensity of the optical beam, and forwards them to the NB-IoT devices after M-PAM to M-QAM mapping; (2) Decode-and-Forward with NB-IoT (a regenerative strategy), where the GEO satellite detects and decodes the symbols of the NB-IoT frames that are transported on the optical feeder link, and forwards them to the IoT devices after the NB-IoT frame regeneration for the downlink radio transmission is over; (3) Decode-and-Forward with NB-IoT/DVB-S2(X) (a regenerative strategy), where the NB-IoT frames are encapsulated into DVB-S2(X) satellite frames for uplink transmission and, in the GEO satellite, the DVB-S2(X) decoding is performed to recover the NB-IoT frame that is then transmitted to the NB-IoT terminals. We note that the previously listed relaying architectures are also applicable to LEO and Medium Earth Orbit (MEO) satellites. However, to adapt the study presented in this paper from GEO satellites with fixed positions in the sky to LEO/MEO satellites that change their position on the sky at different speeds, the modeling of the optical channel and its link budget must be adapted accordingly, reducing the (distance-dependent) path loss attenuation that is experienced at lower orbits but adding an additional effect of the atmosphere at low-elevation angles, as well as pointing errors that may be incurred when the satellite is moving.
The remaining part of this article has been structured as follows: Section 2 summarizes the key concepts to model the MCS defined in the NB-IoT standard, the three proposed relaying architectures using GEO satellites with optical feeder links (one non-regenerative and two regenerative), and the details of the time-packing equalization and low-complexity Log Likelihood Ratio (LLR) computation for soft decoding. Section 3 studies the effect of the turbulent atmosphere in the optical feeder link, with emphasis on the beam wander and scintillation that is introduced in the uplink transmission. Section 4 presents the simulation set up and the figures in terms of end-to-end Block Error Rate (BLER) and throughput. Finally, conclusions are drawn in Section 5.

System Model
This section summarizes the key technological concepts that are needed to model the link level of the NB-IoT communication, as well as the different relaying architectures that could be used on-board the GEO satellite to interface the optical wireless signal (feeder link) into the radio wireless signal (access link) that is forwarded by the GEO satellite to the NB-IoT terminals.

NB-IoT Signal Format for the Satellite Forward Link
NB-IoT has been developed by 3GPP to cope with the large demand on IoT connectivity that is foreseen by the designers of the future generations of mobile communication standards (i.e., 5G and beyond). NB-IoT has been conceived to be deployed in three different configurations or typologies, which are compatible with the spectrum allocation (channelization) that is used in contemporary mobile communication standards such as GSM (2G) and LTE (4G). An overview of these deployment typologies, which are known as stand-alone, in-band, and guard-band, can be appreciated in Figure 1. In the stand-alone deployment, the NB-IoT signal occupies the bandwidth that corresponds to one (or few) of the 200-kHz GSM radio channels of the 2G radio spectrum; note that this strategy is suitable for the re-farming process of the GSM bands. In the in-band deployment, the NB-IoT signal is placed on the radio spectrum that corresponds to one (or few) Physical Resource Blocks (PRB) of an LTE carrier, where each PRB spans over 180 kHz of bandwidth and is formed by 12 sub-carriers of a 15-kHz bandwidth each. Finally, in the third typology of deployment known as guard-band deployment, the NB-IoT signal is placed on the guard bands that are reserved to prevent adjacent-channel interference between LTE carriers. Note that these strategies of deployment do not imply any additional cost and time to enter in service, provided the operator owns a licence either in the 2G/4G radio bands.
According to 3GPP standardization, both NB-IoT uplink and downlink transmissions occupy a communication bandwidth of 180 kHz in the radio portion of the electromagnetic spectrum. Moreover, as the downlink of NB-IoT defines the technology to communicate from the base station (eNB) to the user terminal (IoT device), we focus on this direction of communication to design the radio frame that should be used in the forward link of the IoT satellite system. Specifically, the downlink of NB-IoT uses few 15 kHz sub-carriers, providing a sampling time unit of T s = 1/(15000 × 2048) s., which is identical to the one used in the LTE standard. Similarly, the time slot duration in NB-IoT is T slot = 15360 × T s = 0.5 ms [29]. Two consecutive NB-IoT time slots constitute a sub-frame, which spans 1 ms. Similarly to LTE, a group of 10 subframes with total duration T frame = 10 × 2 × T s = 10 ms constitutes a NB-IoT frame.
The NB-IoT standard enables to repeat the transmission of the same information (block data) up to 2048 times, in order to extend the coverage range and increase the reliability of the data communication [29]. However, the higher the number of repetitions that are performed, the lower the spectral efficiency of the data communication that takes place. The NB-IoT link selects the Transport Block Size (TBS) on a MAC layer from a variety of sizes, which range from 2 bytes (16 bits) up to 317 bytes (2536 bits) [30]. The number of Modulation and Coding Schemes (MCS) that NB-IoT supports is equal to 14, and the combination of a number of subframes and MCS to be used for communication determines the code rate of the NB-IoT transmission. Regarding the error control coding, the downlink of NB-IoT uses a 1/3 tail-biting convolutional encoding mother code [31]. This encoding procedure is formed by three generator polynomials, which are known as the G 0 = 133, G 1 = 177, and G 2 = 165 polynomials in the octal notation (see Figure 2 for more details). Then, after channel encoding, data rate matching is utilized to obtain the desired code rate. This rate-matching procedure is a puncturing process to obtain code rates that are higher than the one provided by the mother code (i.e., code rates higher than 1/3). However, in order to obtain code rates that are lower than the one provided by the mother code, the NB-IoT matching procedure combines block data repetition with puncturing [29]. The downlink of NB-IoT is formed by four channels, namely: Narrowband Physical Downlink Control Channel (NPDCCH), Narrowband Physical Downlink Shared Channel (NPDSCH), Narrowband Physical Broadcast Channel (NPBCH), and Narrowband Synchronization Signals (NPSS/NSSS) [32]. The first channel, the NPDCCH, is used for the control plane and provides the scheduling information for the downlink and uplink data channels. The second channel, the NPDSCH, is used for the data plane and for paging, and contains dedicated and common downlink data. The third channel, the NPBCH, contains information for the initial acquisition conveying information about the cell parameters. The NPDCCH, NPDSCH, and NPBCH channels are QPSK modulated [31]. Finally, the NPSS/NSSS signals are used to perform the cell search, time and frequency synchronization, and cell identity detection procedures, which are modulated using Zadoff-Chu sequences. In this paper, however, we focus the attention on the study of the QPSK-modulated downlink channels.
Finally, the forward link of the satellite relaying system is formed by two links, namely: (i) The link from the gateway to the satellite, the so-called feeder link and (ii) the link from the satellite to the corresponding IoT terminal, the so-called access link. The satellite gateway aggregates the NB-IoT downlink channels and send them to the corresponding satellite beam of the access link. At the satellite, the received data is switched to the access beam of the target IoT device. As expected, the satellite architecture enables to increase the throughput of the system. The following section gives further details on the satellite architectural options that are here studied to forward data from the gateway to the NB-IoT devices.

Architectures for Forwarding NB-IoT Frames over a Satellite Relaying Node
In the coming years, the feeder links of satellite systems will start to introduce optical wireless technology to cope with the capacity demand that new 5G/6G services will require, such a the OTA applications [6]. However, the throughput of the optical links can be increased even further by adding more advanced spectral-efficient communication techniques. Specifically, this paper considers that time-packing can be one the enablers to enlarge the end-to-end throughput of HTS systems. Nevertheless, the architecture of the optical receiver placed on the satellite payload is also important to increase the HTS throughput. Some of these architectures have been already considered in the deployment of non-terrestrial networks, resulting in bent-pipe and fully-regenerative satellite configurations [3]. However, this paper extends the possibilities of fully-regenerative satellites (i.e., Decode-and-Forward) by considering the partial regeneration of NB-IoT frames until the modulation level (i.e., Detect-and-Forward), and by considering the encapsulation of the NB-IoT modulated symbols in DVB-S2(X) frames. Specifically, this paper analyzes the throughput of the following satellite architectures, namely: (i) Case 1: NB-IoT modulated symbols detected at the satellite and forwarded to the corresponding IoT terminal; (ii) Case 2: NB-IoT frames decoded until the bit-level, re-encoded, re-mapped at the satellite, and re-sent to the IoT terminal; and (iii) Case 3: NB-IoT modulated symbols are encapsulated in DVB-S2(X) frames and re-sent to the IoT terminal (see Figure 3 for more details).

Case 1: Detect-and-Forward Relaying of NB-IoT Frames
This strategy is similar to the one used in transparent (bent-pipe) satellite architectures, but with the difference that the NB-IoT amplitude-modulated symbols in the optical feeder link are time-packed. The Inter-Symbol Interference (ISI) that time-packing introduces is removed in the satellite (see Section 2.3) and, after that, the 4-PAM modulated symbols of the NB-IoT signal are detected, remapped to QPSK symbols, OFDM modulated, and re-sent to the IoT terminal over the radio access link. Finally, the IoT terminal computes the LLRs of the QPSK symbols (see Section 2.4) and feeds them in the soft convolutional decoder to estimate the transmitted NB-IoT frames from the gateway (see Figure 3). Let s ntp [k] be the k-th NB-IoT modulated symbols, g tx (·) the transmit pulse-shaping square-root raised-cosine filter, T s the Nyquist symbol time of the NB-IoT signal, and δ the time-packing overlapping factor. Then, the time-packed signal that is generated by the gateway can be written as: Next, the Electrical-to-Optical (E/O) conversion of s tp (t) is done with a MZM driven by voltage: v where V B and V π are the bias and half-wavelength voltages of the MZM, β is the intensity modulation index (scaling factor) that is selected for the communication in the optical feeder link, and: where E{·} is the mathematical expectation operator that determines the mean value of the signal. After that, the time-packed signal is transmitted through the optical feeder link. At the satellite, the Optical-to-Electrical (O/E) conversion is performed with the aid of a Photodetector [28], obtaining: where h fl is the equivalent channel gain of the optical feeder link (see Section 3), s tp [n] is the n-th discrete sample of the time-packed signal in the gateway, and η fl [n] is the resulting noise signal of unit power after the O/E conversion in the satellite (see Section 4.1). The value of (E b /N 0 ) fl represents the equivalent bit-to-noise-energy ratio of the optical feeder link. The estimations of the non-time packed NB-IoT symbols at the satellite (i.e., s ntp,sat ) are obtained after the received signal samples in (4) are first matched-filtered and then equalized (see Section 2.3). Then, the 4-PAM non-time packed signals are re-mapped into QPSK symbols, OFDM modulated, and finally transmitted to the corresponding IoT device. Due to that, the k-th sample of the signal that the IoT device receives in the radio access link attains the form: where η al [k] is the unit power noise signal, h al is the gain of the radio access channel, and (E b /N 0 ) al is the bit-to-noise-energy ratio of the radio access link. Finally, the LLR of the received QPSK modulated symbols are computed, and soft-decoding it is conducted to detect the transmitted bits from the gateway, b gw (see Figure 3).

Case 2: Decode-and-Forward Relaying of NB-IoT Frames
Similar to Case 1 (see Section 2.2.1), the optical receiver of the satellite removes the time-packed interference (see Section 2.3) but, instead of remapping the 4-PAM modulation to QPSK directly, it computes the LLRs of the 4-PAM modulation that are used by the soft-Viterbi decoder to detect the message bits transmitted from the gateway (see Section 2.4). After that, the estimated message bits b gw,sat are re-encoded, QPSK mapped, OFDM modulated, and forwarded to the corresponding IoT device. There, the NB-IoT receiver computes the LLRs of the receive QPSK symbols, to use them in the soft-Viterbi algorithm to detect the transmitted bits from the gateway b gw,iot .
Let s ntp,sat [n] be the estimated symbols of the non-time-packed 4-PAM modulated symbols at the satellite. Then, the following step of Case 2 consists in computing the two LLRs per modulated symbol of the 4-PAM constellation, denoted as LLR b 0 and LLR b 1 (see Section 2.3). After that, the soft-Viterbi decoder at the satellite estimates the original transmitted bits b gw,sat , re-encodes them to regenerate the coded NB-IoT bits c ntp,sat , maps them to QPSK symbols s ntp,sat , and follows as in Case 1 (see Figure 3).

Case 3: Detect-and-Forward Relaying of NB-IoT Frames Encapsulated in DVB-S2(X)
The third case consists in encapsulating the NB-IoT encoded bits into DVB-S2(X) frames. Next, the LDPC coded bits from DVB-S2(X) are 4-PAM modulated, pulse-shaped, time-packed, and sent to the satellite. There, the time-packed interference is removed (see Section 2.3), and the LLRs for the 4-PAM modulation are computed and fed into the soft-LDPC decoder of the DVB-S2(X) [33]. After that, the NB-IoT encoded bits are desencapsulated from the DVB-S2(X) frames, mapped to QPSK symbols, OFDM modulated, and forwarded to the IoT device. There, the LLRs of the QPSK symbols are first computed, and then used by its soft-Viterbi decoder to detect the message of the gateway.
Let c iot,gw be the NB-IoT convolutional encoded signal at the terrestrial gateway, which is encapsulated into the DVB-S2(X) frame. Due to the difference in the payload size of a NB-IoT and a DVB-S2(X) frame, several NB-IoT convolutional codewords can be packed together into the input message of the DVB-S2(X) physical layer frames. We recall that at the physical layer, a DVB-S2(X) frame is first BCH encoded and then LDPC encoded. The BCH encoding is used to remove the possible error floors of the LDPC decoding, which aims at compensating the impairments that introduce the communication channel (i.e., the optical feeder link in our case). Let us assume that c dvb,gw denotes the input frame to the DVB-S2(X) encoder at the gateway. Then, it can convey up to P NB-IoT convolutionalcoded frames as c dvb,gw = c iot,gw (0) · · · c iot,gw (P − 1) . Next, the DVB-S2(X)-coded frames are 4-PAM modulated, E/O converted, and sent to the satellite on the optical feeder link. At the satellite, the O/E conversion is carried out, and the received signal samples are matched-filtered and equalized to mitigate the ISI introduced by time-packing. Likewise, in Case 2, the two LLRs are computed per 4-PAM symbol, i.e., LLR b 0 and LLR b 1 (see Section 2.3). However, in Case 3 these LLRs are introduced into the soft LDPC decoder, not in the soft-Viterbi decoder as it was used in Case 2. After that, the NB-IoT encoded bits are de-encapsulated from the decoded DVB-S2(X) frames, QPSK mapped, OFDM modulated, and forwarded to the IoT terminal over the radio access link. Finally, at the IoT device, the processing for detecting the transmitted message is the same as the one that was explained for Case 1.

Equalization of the Time-Packed Signal
The three proposed relaying architectures rely on time-packing signalling in the optical feeder link to increase the throughput even further. Unfortunately, the use of time-packing introduces ISI that must be mitigated in reception at the satellite. The optimal strategy for cancelling this unwanted ISI consists in resorting the Maximum-Likelihood Sequence Decoding (MLSD) [34], which can be implemented efficiently with the Viterbi algorithm [35]. However, to increase the throughput of the optical feeder link, it is necessary use a low roll-off factor and large overlapping factor, two things that increase notably the complexity of the Viterbi Algorithm to be implemented in the satellite [27]. Due to the impracticability of this solution, alternative strategies should be considered. Towards this regard, in this paper we use a two-side Least Mean Square (LMS) filter to equalize the time-packed channel. Thus, the equalization strategy consists of two steps: (1) Compute the weights of the interference canceller by means of a training sequence and (2) apply the pre-computed weights to the received time-packed signal from the gateway. These weights have to be computed for each overlapping factor that can be used.
Let us assume that r is the vector that stacks the received samples of the time-packed signal, T is the vector of L weights in the LMS equalizer, and y[n] is the buffer that contains the received samples that participates in the equalization process of the n-th time-packed received symbol. Then, the equalization of the n-th transmitted symbol attains the form: During the training process, the ideal values of the signal samples at the output of the equalizer z[n] are known, enabling to determine the most convenient equalization weights. Specifically, the weights w have been computed by using the LMS algorithm strategy [26,36], i.e., where w[p] contains the equalizer weights at the p-th training iteration, µ is the forgetting factor, and is the error between the training symbol z t [p] and its estimation from (6).

Computation of the LLRs for the 4-PAM and the QPSK Modulation Schemes
In this paper, the NB-IoT transmitted data is encoded with a convolutional coding scheme. Furthermore, for Case 3, the encapsulated NB-IoT frames are protected using the LDPC code of the DVB-S2(X) standard. In both cases, soft-decoding is used and due to that, the LLRs has to be determined. However, the IoT devices and the satellite node may be limited in their power consumption. Fortunately, for both modulation schemes used in this paper (i.e., 4-PAM and QPSK), it is possible to derive reduced complexity closed-form formulas for computing the optimal LLRs. The 4-PAM modulation scheme is used in the optical feeder link, whereas the QPSK modulation is used in the radio access link. We assume that Gray-mapping is used in both cases, and that their corresponding modulation symbols are Figure 4). Since both modulation schemes transport two bits per modulated symbol, it is necessary to compute two LLRs per modulated symbol. Thus, the closed-form expression for computing these LLRs is given by: where LLR b m is the LLR that corresponds the m-th bit of the modulated symbol, σ 2 n is the noise power, M is the number of constellation symbols, h is the communication channel, z is the received data at the input of the de-mapper, and s p (s q ) symbolizes the constellation symbols in which the m-th bit is 0 (1). Thus, according to (8) and the Gray mapping proposed in Figure 4, the LLR of the first bit b 0 for both 4-PAM and QPSK modulation schemes is given by: whereas LLR of the second bit b 1 attains the form: If (9) and (10) are computed for 4-PAM and QPSK modulation schemes, the closedform expressions for LLR b 0 and LLR b 1 are shown in Table 1. In these closed-form formulas, we have that: where s 1 is the second symbol of the 4-PAM constellation, i.e., s 1 = −1/ √ 5 (see Figure 4) and signal z represents the data after equalizing the time-packed signal (see Section 2.3). For Case 2, this signal z corresponds to the NB-IoT one, and its LLRs are introduced in the soft-convolutional decoder for recovering the transmitted message to regenerate the NB-IoT signal. For Case 3, the signal z corresponds to the DVB-S2(X) after the time-packing equalization. The computed LLRs for this case is introduced to the LDPC decoder of the DVB-S2(X) receiver [33]. After that, the encoded NB-IoT bits are QPSK mapped, OFDM modulated, and forwarded to the IoT-device. Finally, at the IoT receiver, the LLRs of the received QPSK modulation symbols are computed for all cases under study. These LLRs are used by the soft-Viterbi decoder to recover the message transmitted from the satellite gateway. Table 1. Closed-form expressions for the LLR b 0 and LLR b 1 of both 4-PAM and QPSK modulations. Figure 4. Gray mapping between 4-PAM and QPSK symbols performed in the satellite node.

Optical Wireless Satellite Channel Model
The optical signal that is used to transport the data symbols from the satellite gateway to the GEO satellite must go through the different layers of the Earth's atmosphere. Unfortunately, the power loss that the optical feeder link experiences in the uplink direction of communication is larger than the one observed in downlink. This is because, in the ground-to-satellite communication, the optical signal starts to spread and accumulate distortion as soon as it leaves the satellite gateway transmitter.

Atmospheric Power Losses: Absorption and Scattering Modeling
The power loss of the optical feeder link is a function of the atmospheric attenuation, which depends on both absorption and scattering effects that the light signal experiences while propagating [37]. To compute this value, the atmospheric attenuation coefficient: needs to be computed, where α m and α a are the molecular and aerosol absorption coefficients, respectively, whereas β m and β a are the molecular and aerosol scattering coefficients, respectively.

Modelling of absorption:
At Infrared (IR) wavelengths, the principal atmospheric absorbers are the molecules of water, carbon-dioxide, and ozone. As expected, the atmospheric absorption is a wavelength-dependent phenomenon. Therefore, the operating wavelength for optical feeder link transmissions should be chosen to minimize this loss, using the atmospheric transmission windows in which the molecular and aerosol absorption is less than 0.2 dB/km for clear sky conditions [38]. In addition to the low-absorption requirements, most optical feeder links are designed to work in the 780-850 nm and 1520-1600 nm windows because there are off-the-shelf lasers and detectors commercially available to work in these wavelengths.

Modeling of scattering:
Like absorption, scattering is also a phenomenon that is strongly dependent on the operating wavelength. If the size of the atmospheric particles is small in comparison to the optical feeder link wavelength, then Rayleigh scattering is produced. Particles like air molecules and haze cause Rayleigh scattering [39] and affect notably optical wireless transmissions in the Visible Light (VL) and Ultraviolet (UV) regions; on the other hand, Rayleigh scattering can be neglected for optical feeder link wavelengths in the IR range (i.e., when λ 1 µm). Similarly, when the atmospheric particles size is comparable with the operating wavelength, then Mie scattering is produced. Aerosol particles, fog, and haze are the major contributors of Mie scattering, and this phenomenon is dominant for wavelengths in the IR range. Finally, if the atmospheric particles are much larger than the operating wavelength, the scattering is better described by geometrical optical models, which should be used in case of rain, snow, and hail weather conditions [40].
Modeling the transmittance of the Earth's atmosphere: Apart from the previously described λ-dependent effects, the specific value that the atmospheric attenuation coefficient takes depends on the concentration of molecules (and aerosols) of the Earth's atmosphere at different altitudes h. Based on this assumption, the atmospheric transmittance that an optical feeder link with zenith angle ζ experiences is given by: where h 0 is the altitude of the satellite gateway (ground station) over the sea level, H is the vertical height at which the GEO satellite is placed, and γ(λ, h) is the attenuation coefficient at wavelength λ and altitude h. Based on this formula, it is possible to see that atmospheric transmittance is increased at low zenith angles (i.e., at high elevation angles), as the fraction of the incident electromagnetic power that is transferred through the atmosphere layers is increased.
Power losses due to fog: From the common weather conditions, fog is the one that contributes most in the absorption and scattering of the optical signal when it propagates through the Earth's atmosphere. In the presence of fog, the optical feeder link connectivity is seriously put at risk, particularly when the fog layer next to the ground station extends vertically very high, forming a fog layer that can be as thick as 400 m over the Earth's surface. In such critical weather conditions, the use of very high power lasers (1550 nm) with special mitigation techniques is the only option to maximize the chances of optical feeder link connectivity. As an alternative method to the Mie scattering theory, the attenuation due to fog for different wavelengths can be estimated using empirical models that use as the input parameter the visibility in km measured on the VL region (550 nm). For a comparison of the fog attenuation at different wavelengths (850 nm and 950 mn), please refer to [41]. Note that in extreme cases, where the visibility due to fog is reduced to about 50 m, atmospheric attenuation can be as high as 350 dB/km [42].
Power losses due to rain: The impact of rain in the propagation of optical signals is not as pronounced as fog, because the size of the rain droplets are significantly larger in size (from 100 µm to 1000 µm) when compared to operating wavelengths of optical feeder links. For example, the attenuation loss in light rain (2.5 mm/h) to heavy rain (25 mm/h) ranges from 1 dB/km to 10 dB/km for 850 nm and 1550 nm operating wavelengths, respectively [43]. Note that the low clouds, which usually accompany rainy weather, are the source of strong attenuation in most optical feeder links. In order to combat the huge power loss that takes place in such conditions, it is recommended to include a few tens-of-dB margin (e.g., 30 dB) when designing the link budget of the optical wireless link. Moreover, optical feeder link designers can also implement adaptive coding and modulation schemes to address the varying weather conditions in the geographical area around the ground station [26].
Power losses due to snow: Finally, since the size of snow droplets is between the size of rain and fog droplets, the atmospheric attenuation for dry/wet snow conditions is usually stronger than the one in the presence of rain, but not as severe as the one in case of fog. However, during heavy snow storms, the path of the optical feeder link can be completely blocked for the presence of densely-packed snow flakes in the propagation path. In such cases, the attenuation is similar to the one observed in foggy weather (30-350 dB/km) and, as expected, can seriously put at risk the optical feeder link connectivity.

Atmospheric Turbulence: Beam Wander, Beam Spreading, and Beam Scintillation
Atmospheric turbulence is a random phenomenon that is caused by the variation of the temperature and pressure on the atmosphere layers that are in the propagation path of the optical wireless signal. These temperature and pressure inhomogeneities form turbulent cells, known as eddies, which have different sizes and different diffractive indexes. The eddies act as if prisms/lenses were deployed in the propagation path, introducing constructive and destructive interference in the received optical signal. The perturbations that atmospheric turbulence introduces in the wave-front of the optical beam can be physically described by the Kolmogorov model [44]. Depending on the size of the turbulent eddies with respect to the transmitted beam size, three types of atmospheric turbulenceinduced effects can be identified, namely: beam wander, beam spreading, and beam scintillation.

Turbulence-induced beam wander:
This phenomenon takes place when the size of the turbulent eddies is larger than the size of the optical beam. Beam wander results in a random deviation of the optical beam from its planned (rectilinear) propagating path and, in extreme displacement situations, may lead to the failure of the optical wireless link. Beam wander is a major concern in the uplink transmission of an optical feeder link, as the beam size in the ground-to-satellite transmission is often smaller than the size of the turbulent eddies, resulting in a beam displacement at the receiver side that can be as large as several hundred meters.
In case of a collimated beam (plane wave model), the Root Mean Square (RMS) displacement due to beam wander for an uplink path with zenith angle ζ can be written as: where H is the altitude of the GEO satellite (receiver), h 0 is the altitude of ground station (transmitter), W 0 is the initial beam size, and r 0 is the atmospheric coherence width, which is also known now as Fried's coherence length, Fried's parameter, or simply coherence length [45]. The Fried's coherence length is a widely-used descriptor of the level of atmospheric turbulence at a particular site and, for a known structure constant profile C 2 n (h) and plane wave model (collimated beam), it is given by: [46] where k = 2π/λ is the wavenumber of the optical beam. As expected, C 2 n (h) varies with the time of the day, the geographical location, and the altitude. Therefore, for vertical optical links (slant paths), the value of C 2 n (h) has to be integrated over the complete propagation path, starting from the height of the ground station above the sea level and ending at the altitude in which the Earth's atmosphere vanishes (i.e., at about 40 km).
Various empirical models for C 2 n (h) have been proposed in the literature to estimate the turbulence profiles, using as reference the experimental measurements that were carried out at different geographical locations, time of the day, wind speed, terrain types, among others. The most widely-used model to characterize the refractive index structure of the atmosphere for vertical links (slant paths) is the so-called Hufnagel-Valley (H-V) model [45], i.e., where h [m] is the altitude, v [m/s] is the RMS wind-speed, and parameter A 0 = C 2 n (h 0 ) [m −2/3 ] is the nominal value of the refractive index near the ground level. The RMS wind speed in (17) is determined from the Bufton wind model, and can take values that range from v = 10 to 30 m/s for moderate and strong wind speeds, respectively. Similarly, the ground turbulence level can take values between A 0 = 1.7 × 10 −14 and 1.7 × 10 −13 m −2/3 , which depends on the location and day time, among other parameters.
When using A 0 = 1.7 × 10 −14 m −2/3 and v = 21 m/s, this model is commonly referred to as the H-V 5/7 model because, for wavelength λ = 0.5 µm and a transmitter on the ground looking up (i.e., with ζ = 0 deg.), it predicts a value of atmospheric coherence diameter r 0 = 5 cm according to (16) and a value of an isoplanatic angle: of 7 µrad in case of a spherical wave with output-plane beam parameters Θ = Λ = 0. The refractive index profile along the vertical/slant path is shown in Figure 5 for two nominal values of a refractive index at the ground level and three different RMS wind speeds. From this figure, it is possible to observe that the ground turbulence level A 0 has little effect above 1 km, and that the wind speed governs the profile behavior primarily in the vicinity of altitudes in the 10 km range. Similarly, Figure 6 shows the RMS angular displacement due to beam wander (σ 2 BW ) as a function of the beam radius W 0 , when operating wavelength λ = 1.55 µm and refractive index structure C 2 n (h) follows the H-V 5/7 model. As expected, the RMS beam wander displacement is higher for the largest zenith angle, as the section of the atmosphere through which the optical beam needs passes through is thicker, the beam deviation with respect to the straight path grows. Finally, according to (16), the Fried's coherence length is r 0 = 19.25 and 12.70 cm for zenith angle ζ = 0 and 60 deg., respectively.

Turbulence-induced beam spreading:
This phenomenon takes place when the turbulent eddies are smaller than the size of the optical beam. Beam spreading generates a widening of the beam size, beyond the natural broadening due to diffraction that the non-turbulent atmosphere introduces. Beam spreading does not affect the direction of the optical beam but, in contrast, reduces the optical power at the receiver aperture due to the energy dispersion that takes place.
Turbulence-induced beam scintillation: When the size of turbulent eddies is of the same order of the size of the optical beam, then the eddies act as lenses that focus and defocus the incoming beam. In this situation, the eddies lead to a redistribution of the signal energy that generates a temporal and spatial fluctuation of the irradiance at the receiver aperture. This phenomenon, which is known as scintillation, represent one of the major sources of degradation in the performance of an optical feeder link. Atmospheric turbulence also leads to loss of spatial coherence of an initially coherent optical beam, and may also produce depolarization of the light and temporal stretching of the optical pulse.
The atmospheric scintillation is measured in terms of the scintillation index, which is the normalized variance of the intensity fluctuations, i.e., where I is the irradiance (intensity) in the detector plane and · denotes the ensemble average. The Gamma-Gamma distribution has been proposed to describe the turbulence-induce scintillation over a broad range of beam diameters. The Probability Density Function (PDF) of the Gamma-Gamma turbulence model and the scintillation index are given by: respectively, where I m denotes the mean irradiance, Γ(x) is the Gamma function, and K n (x) is the modified Bessel function of the second kind. The parameters α = 1/σ 2 X and β = 1/σ 2 Y of the Gamma-Gamma distribution in (20) are directly related to the atmospheric conditions, and for the untracked beam case are given: The various parameters that appear in Equations (21) and (22) are defined as follows: and are the jitter-induced angular pointing error, the pointing error variance, slant path length, and scaling constant, respectively. Similarly, the diffractive beam radius at the receiver is given by: are the input-plane beam parameters. Note that for a collimated beam, the phase front radius of curvature at the transmitter output aperture F 0 → ∞ and, due to that, Θ 0 ∼ = 1. Finally, the irradiance flux variance in the focal plane of the receiver: where: is the normalized distance for the uplink propagation path, and: are the output-plane beam parameters. In Figure 7 we plot the corresponding Gamma-Gamma PDF for three different beam sizes W 0 , which are equal to 1, 10, and 50 cm. Once again, the wavelength was set to λ = 1.55 µm and the analysis was done for the uplink direction of communication of a perfectly vertical GEO satellite feeder link (i.e., ζ = 0 deg. and H = 36,000 km) when using the H-V 5/7 refractive index model (i.e., v = 21 m/s). Note that in this situation, the Fried's coherence length is r 0 = 19.25 cm. As expected, for small beam sizes in which the 2W 0 /r 0 1 relationship is verified (e.g., similar to W 0 = 1 cm in Figure 7), the longitudinal component of the scintillation index will be much less than 1; due to that, the corresponding PDF of the normalized irradiance will have a shape that resembles the one of a log-normal distribution, but with some differences in the upper and lower tails. On the other hand, for large beams in which the 2W 0 /r 0 1 relationship is observed (e.g., similar to W 0 = 50 cm in Figure 7), the scintillation index becomes larger than 1 and the shape of the PDF starts to resemble a negative exponential distribution.

Evaluation
The error correction capabilities that the Modulation and Coding Schemes (MCS) of the NB-IoT (and DVB-S2(X)) standard have on the end-to-end forward link of the GEO satellite system (i.e., from the satellite gateway to the IoT terminals) is now evaluated in detail. For this purpose, we first present the simulation setup and, after that, we show the different figures of merit that are relevant to characterize the end-to-end performance of the three GEO satellite relaying strategies.

Simulation Setup of the Optical Channel
According to the analysis presented in [28], the mean SNR of the electrical signal that is direct-detected by the PD that is placed on-board the satellite is given by where is the DC component of the time-varying electrical current i D (t) that the PD generates when being excited by the intensity modulated optical signal, β is the intensity modulation index, and includes the contribution of all noise sources in the optical feeder link, namely the shot noise sources, thermal noise, Relative Intensity Noise (RIN) of LD, and beat noise [25]. Note that shot noise term includes the contribution of the received optical signal, the Amplified Spontaneous Emission (ASE) noise, the background optical noise and the dark current noise, whereas the beat noise term accounts the effect of combining the received optical signal with the ASE noise. When the received optical power is between −90 and −20 dBW, it can be shown that the beat noise between received optical signal and ASE noise dominates the SNR of the optical feeder link [47]. In this situation, where B o is the bandwidth of the optical signal at the PD input, B e is the bandwidth of the electrical signal at the PD output, and I ase = µ G o,edfa P ase is the DC component generated by the ASE noise, whose equivalent noise power at the input of the EDFA is given by P ase = ρ ase B o . Table 2 summarizes the parameters of the optical feeder link, taking into account both the optical gains and losses, as well as the different sources of optical noise [25,48]. Note that the values of the parameters that appear in this table have been determined taking into account the different phenomena described in Section 3, when modeling the contribution of each of them in the space-to-ground optical link budget. The effect of any other parameter not listed in this table is considered as negligible. When L o,atm = 0 dB (i.e., clear-sky conditions), the DC current at the PD output is I D = 75.68 mA, whereas the DC current generated by the ASE noise is I ase = 0.125 mA regardless of the weather. Thus, when we set the intensity modulation index β = 0. 5 Note that for larger intensity modulation indexes β, the non-linear distortion introduced by the MZM of the optical transmitter becomes more notable. For those situation, the use of Digital Pre-Distortion (DPD) is necessary to keep non-linear distortion under control [25].  Figure 8 shows the Cumulative Distribution Function (CDF) of the power loss that turbulence-induced beam wander and beam scintillation introduce on the received optical power of the optical uplink transmission in the satellite forward link (see Section 3). Without loss of generalization, it is assumed that the ground station is placed next to the city of Madrid, Spain (i.e., Lat. 40.43 • North, Long. 4.25 • West, and altitude h 0 = 864 m), with an expected system availability due to cloud-free line-of-sight conditions close to 99.8%. Similarly, the satellite is assumed to be placed in the GEO position of 19.2 • East at an altitude of H = 36,000 km, giving as result a zenith angle ζ = 52.61 deg. for the optical feeder link in the forward direction of communication. The nominal value of the refractive index near the ground level is A 0 = 1.7 × 10 −14 m −2/3 , and the CDF is plotted for three different wind speeds (i.e., v = 10, 21 and 30 m/s). Note that when v = 21 m/s, the refractive index model becomes the well-known H-V 5/7 model. Finally, an untracked collimated beam is assumed, with a beam radius of W 0 = 10 cm at a wavelength of λ = 1550 nm. As expected, the stronger is the wind speed, the more variability is introduced in the instantaneous intensity of the received optical signal at the GEO satellite. For example, the 10-th percentile for the turbulence-induced power loss (i.e., power loss margin for a 90% link availability without coding) is equal to 2.4, 3.0 and 3.8 dB for a wind speed of 10, 21 and 30 m/s, respectively. These turbulence-induced power losses grow to 4.3, 5.4 and 7.0 dB when we study the 1-st percentile (i.e., power loss margin for a 99% link availability without coding) for the same three values of wind speeds (i.e., 10, 21 and 30 m/s, respectively).

Simulation Setup of NB-IoT Signal
This paper evaluates the BLER and Throughput of the following three forward link architectures for transmitting NB-IoT over satellite when the feeder link is optical: In all cases, the NB-IoT frame length is set to N = 1032 bits. The pulse-shaping filter for NB-IoT and DVB-S2(X) waveforms is assumed to be Square-Root Raised-Cosine (SRRC) with roll-off factor ρ = 0.15. The overlapping factors that have been considered in the simulations to implement the time-packing signaling are δ = 0, 15, 25 and 40%. Note that when δ = 0% (i.e., no overlapping), the symbol time is equal to the Nyquist symbol time. The modulation scheme used in the optical feeder links is 4-PAM, in order to make it compatible with the Intensity Modulation/Direct Detection approach. On the contrary, the access link uses QPSK, which is one of the modulation schemes considered for the downlink in the NB-IoT standard. So, in both optical feeder and radio access links, the number of modulation states is M = 4 and the number of bits per modulated symbol is In order to make fair comparisons, all the relaying satellite architectures under analysis use the same data rate. In particular, the BLER and throughput of the different satellite relaying strategies have been tested with nine possible code rates, namely: Three low code rates that are more robust than the mother code rate of the convolutional encoding (i.e., 0.2222, 0.259 and 0.3), three medium code rates that provide similar error correction capabilities than the mother code rate (i.e., 0.4444, 0.5319 and 0.6), and three high code rates that are the least robust (i.e., 0.6386, 0.7678 and 0.854). For the first and second use cases, these code rates correspond to the code rate of the NB-IoT waveform. However, for the third case, these code rates represent the joint code rate of the DVB-S2(X) and NB-IoT waveforms. This joint code rate is defined as the product between the code rates of the DVB-S2(X) and NB-IoT signals as In (33), the code rates of NB-IoT and DVB-S2(X) have been selected such that the their joint code rate R c,case 3 is equal to the code rates of NB-IoT for the first and second cases. Thus, three code rates (i.e., low, medium and high) have been considered for each waveform (i.e., NB-IoT and DVB-S2(X)), which enable to obtain the nine code rates for cases 1 and 2.
Specifically, the code rates of the DVB-S2(X) and NB-IoT standards that have been used are: R c,dvb = 0.66667, 0.8 and 0.9 [33] and R c,iot = 0.3333, 0.66667 and 0.96, respectively. Thus, we obtain nine possible combinations of the code rates, which are values that are similar to the ones used in the three relaying architectures (see Table 3 for more details). Therefore, after obtaining a similar code rate for the three cases under study, the architecture that provides the best in terms of BLER and throughput is determined for the different SNR working regimes in the optical feeder link (uplink) and radio access link (downlink), respectively. Regarding the throughput, the following closed form formula has been used to estimate its value, i.e., where the code rate R c for the two first cases under study corresponds to the code rate of the NB-IoT waveform. For the third case, the value of R c represents the joint code rate of the NB-IoT and DVB-S2(X) waveforms, which is given in (33). Without loss of generality, the roll-off factor is kept fixed in all simulations (i.e., ρ = 0.15). Thus, according to the overlapping factors and code rates that have been previously presented, the maximum achievable throughput in each situation are shown in Table 4. This maximum achievable throughput may not be reached due to the presence of a residual time-packing interference, unexpected variations in the optical channel gain, and noise in the optical feeder link may increase the BLER. So, simulation results give information about the transmission structure that provides the best in terms of both BLER and throughput. Note that lowest code rates offer a very low BLER but, in contrast, worst results in terms of throughput. On the contrary, high code rates may provide the highest achievable throughput but, at the same time, the worst BLER. For this reason, a target BLER ≤ 10 −2 has been considered as requirement for a normal service according to what is stated in the 5G specifications. This condition helps to decide the best transmission architecture, such that the case that satisfies the target BLER and provides the best throughput will be categorized as the most convenient transmission architecture. However, in order to compare the three cases under study, additional figures of merit needs to be defined in the following section.

Definition of the Figures of Merit
In the previous section, we defined the throughput in terms of the overlapping factor and code rate. However, in the evaluation process of the three cases, it is important to identify the envelope of the throughput defined in (34). Specifically, we will study the envelope of the throughput when the overlapping factor δ is kept constant. This throughput, denoted as THR(δ) * , is defined as the maximum throughput for all code rates with the same overlapping factor, and is computed as follows: Next, it is also interesting to determine the envelope of the throughput for the no-timepacking and time-packing signaling schemes. For the case of no-time-packing signaling, the envelope of its throughput is obtained by setting δ = 0 in (35), that is, THR ntp = THR(0) * . For the case of time-packing signaling, the envelope of the throughput is determined by computing the maximum of (35) for all overlapping factors that are higher than zero (i.e., δ > 0), according to Then, the overall envelope of the throughput for a particular relaying strategy is obtained from the maximum envelope throughput of the time-packed and non-time packed ones, i.e., THR case q = max THR ntp , THR tp q = 1, 2, 3.
After defining the throughput envelope in different working conditions, we now introduce two relative figures of merit that evaluate the throughput gain of the relaying strategies. The first one determines whether using time-packing is beneficial or not with respect to non-time-packing. This parameter, which provides the percentage gain in throughput, attains the form The second figure of merit aims at determining the relative gain in throughput of the different relaying strategy. From the three cases under analysis, the first one represents the simplest approach, since the satellite does not manipulate the content of the NB-IoT frame beyond the symbol-by-symbol mapping that is needed to adapt the optical wireless transmission (M-PAM) to radio wireless (M-QAM). Therefore, it is considered as the benchmark system. Then, the relative throughput gain that cases 2 and 3 provide with respect to case 1 is given by ∆G thr,(q,1) (%) = THR case q − THR case 1 × 100 THR case 1 q = 2, 3, where THR case q is the throughput envelope for the three relaying strategies under study, computed with the aid of (37). Finally, after defining the figures of merit to carry out the analysis, we are ready to present the obtained simulation results.

Simulation Results
The results provided in this section evaluate the BLER and throughput of the three aforementioned NB-IoT satellite relaying architectures in the following aspects:

1.
Evaluation of the end-to-end BLER and throughput when the E b /N 0 of the feeder link varies from 2 to 22 dB in steps of 1 dB (see Figure 9), the SNR of the access link is 20 dB, the wind speed is 21 m/s, and the NB-IoT code rates are R c,iot = 0.3333, 0.66667, and 0.95556. There is no time-packing (δ = 0%) and configuration is Case 1 (Detect-and-Forward NB-IoT). The throughput of Figure 9 is the THR 0, R c,iot of (34); 2.
Evaluation of the end-to-end BLER and throughput when the E b /N 0 of the access link varies from −8 to 12 dB in steps of 1 dB (see Figure 10), the SNR of the feeder link is 15, 20, and 25 dB, the wind speed is v = 10, 21, and 30 m/s, and the NB-IoT code rates are R c,iot = 0.3333, 0.66667 and 0.95556. There is no time-packing (δ = 0%) and the tested configuration is Case 1 (Detect-and-Forward NB-IoT). The throughput of Figure 10 is THR 0, R c,iot of (34); 3.
Evaluation of the throughput for Case 1 (see Figure 11): Detect-and-Forward NB-IoT when the SNR of the feeder link is 15 dB, the E b /N 0 of the access varies from −8 to 12 dB, and the code rates of the NB-IoT and overlapping factors are as defined in Section 4.2. The throughput of Figure 11 corresponds to THR δ, R c,iot of (34); 4.
Evaluation of the throughput for Case 2 (see Figure 12): Decode-and-Forward NB-IoT when the SNR of the feeder link is 15 dB, the SNR of the access varies from 0 to 20 dB, and the code rates of the NB-IoT and overlapping factors are as defined in Section 4.2. The throughput of Figure 12 corresponds to THR δ, R c,iot of (34); 5.
Evaluation of the throughput for Case 3 (see Figure 13): Detect-and-Forward NB-IoT encapsulated in DVB-S2(X) frames when the SNR of the feeder link is 15 dB, the SNR of the access varies from 0 to 20 dB, and the code rates of the DVB-S2(X) and NB-IoT and overlapping factors are defined as in Section 4.2. The throughput of Figure 13 corresponds to THR δ, R c,iot of (34); 6.
Evaluation of the envelope of the throughput when the overlapping factor δ remains constant for all cases under study (see Figure 14a-c) when the SNR of the feeder link is 15 dB, the SNR of the access varies from 0 to 20 dB, and the code rates of the DVB-S2(X)/NB-IoT frames and the overlapping factors are defined in Section 4.2. The throughput of these figures corresponds to THR(δ) * of (35); 7.
Evaluation of the gain in throughput of time-packing signaling with respect to notime-packing, as shown in Figure 14d. The SNR of feeder link is 15 dB, the SNR of the access varies from 0 to 20 dB. This gain is denoted as G thr (%) and is computed following (38); 8.
Comparison among the envelope of the throughput for all cases under study (see Figure 15a-d) in terms of the overlapping factor δ, when the SNR of feeder link is 15 dB, the SNR of the access link varies from 0 to 20 dB, and the code rates of the DVB-S2(X)/NB-IoT frames and overlapping factors are defined as in Section 4.2. The throughput of these figures corresponds to THR(δ) * of (35); 9.
Comparison of the envelope of the throughput for all cases under study (see Figure 16a) when the SNR of the feeder link is 15 dB, the SNR of the access link varies from 0 to 20 dB. The throughput plotted in this figure is denoted as THR case q and is computed according to (37); 10. Comparison of the gain in throughput offered by Cases 2 and 3 with respect to Case 1 (see Figure 16b) when the SNR of the feeder link is 15 dB, the SNR of the access link varies from 0 to 20 dB. The gain in throughput of this figure is denoted as ∆G thr,(q,1) (%) with q = 2, 3, and is determined following (39).
From the simulations of the end-to-end BLER and throughput in terms of the SNR in the optical feeder link (uplink), assuming a very high SNR in the radio access link (see Figure 9), it is possible to observe the region in which the uplink communication limits the performance of the end-to-end transmission. In all these plots, the dotted lines represents the cases where the BLER is higher than the target BLER (i.e., BLER > 0.01), whereas the continuous lines identify the situations in which the BLER is lower than the target BLER (i.e., BLER ≤ 0.01). For a better indication, the grey area in the figure shows the region of SNRs in which the tested code rates for NB-IoT satisfies the target BLER. Thus, it is possible to observe that R c,iot = 0.3333 provides the best throughput when the SNR of the feeder link is ranged between 13 and 19 dB, R c,iot = 0.6666 is the best option when the feeder link SNR is between 19 and 22 dB, and finally R c,iot = 0.95556 is the most convenient alternative when the SNR in the uplink direction of communication is larger than 22 dB. According to these results, we now consider that the SNR of the feeder link is 15, 20, and 25 dB, for determining the degradation that different wind speeds introduce on the end-to-end link.
(a) BLER for δ = 0 %. SNR al = 20 dB.   According to these results, the higher is the SNR of the optical feeder link, the lower is the degradation that the wind introduces on the end-to-end BLER and throughput of the forward NB-IoT link. Consequently, a SNR of 15 dB is used as reference value in the optical feeder link from now on, to compare the performance of the three satellite relaying architectures under study.      Figure 11 presents the throughput for each overlapping factor and code rate for Case 1 (Detect-and-forward with NB-IoT). Simulation results show that the segment of the NB-IoT code rates that satisfies the target BLER is the lowest one (i.e., R c,iot = 0.2222, 0.259, and 0.3) for all the overlapping factors under study. The medium code rate of R c,iot = 0.4444 could only achieve its maximum throughput with no overlapping (i.e., when δ = 0%). Medium and larger code rates are mainly penalized by the SNR of the optical feeder link (uplink), and they are not able to achieve their maximum throughput (see Table 4). When comparing time-packing versus no-time-packing, results show that time-packing signaling provided a better throughput, and the gain was depended on the E b /N 0 of the access link.      Figure 12 shows the throughput for each code rate and overlapping factor for Case 2 (Decode-and-forward with NB-IoT). From this figure, it is possible to observe that the regeneration of the NB-IoT signal in the satellite permits to use also medium code rates (i.e., R c,iot = 0.4444) with overlapping factor δ = 0 and 15%. Therefore, it means that the largest throughput is not obtained with the largest possible overlapping factor (i.e., δ = 40%) and low code rates (i.e., R c,iot = 0.3) as in the previous case, but with a low-to-moderate overlapping factor (i.e., δ = 15%) and a medium code rate (i.e., R c,iot = 0.4444).     Similarly, Figure 13 shows the throughput for each code rate and overlapping factor for Case 3 (Decode-and-forward with NB-IoT and DVB-S2(X)). Results of this figure show that the best code rate to use in the optical feeder link is R c,dvb = 0.6666. For the NB-IoT system, in contrast, the best option is R c,iot = 0.3333 when the E b /N 0 in the access link varies between −1 and 5 dB, R c,iot = 0.6666 when the E b /N 0 in the access link grows from 5 to 9 dB, and finally R c,iot = 0.96 for E b /N 0 values beyond 9 dB in the access link.     Figure 14a-c show the envelope of the throughput for Cases 1-3 and the different overlapping factors under study. From these figures, it is possible to observe that for Cases 1 and 3, the larger is the overlapping factor, the higher is the maximum throughput that is achieved when δ = 40%. On the contrary, for Case 2, the maximum throughput is achieved when the overlapping factor is δ = 15%. Figure 14d shows the percentage gain in throughput of the time-packed versus no-time-packed signalling. From this subfigure, it is possible observe that for when the access link E b /N 0 is between −2 and 5 dB, the percentage gain of Case 1 is in the order of 50%. For other values of E b /N 0 , this gain falls down to about 10%. For Case 2, the percentage gain on the throughput envelope of time-packed schemes is larger than 40% when the E b /N 0 varies between −4 and 1 dB. However, for values of E b /N 0 larger than 1 dB in the access link, this gain in the throughput envelope reduces to about 20%. For Case 3, we have that for an E b /N 0 larger than 0 dB, the minimum gain in the throughput envelope is 40%. Though Case 3 provides the largest gain, the minimum E b /N 0 in the access link to use time-packing is larger than in Cases 1 and 2.      Figure 15 shows the envelope of the throughput for all cases in terms of the overlapping factor. From these figures, it is possible to observe that for Cases 1 and 3, the larger is the overlapping factor, the larger is the maximum throughput that is feasible when δ = 40%. On the contrary, for Case 2, the maximum throughput is achieved when δ = 15%.    Next, Figure 16 presents the envelope of the throughput for all cases (Figure 16a) and the percentage gain of Cases 2 and 3 with respect to Case 1 (Figure 16b). Studying Figure 16a, it is possible to conclude that Cases 2 and 3 have a better throughput than Case 1 in most cases. Thus, Case 1 is used as benchmark to assess the goodness of Cases 2 and 3. From this figure, it can also be noted that the regeneration of the NB-IoT signal enables to lower the minimum E b /N 0 to about 1.5 and 4 dB for Case 2 and Case 3, respectively. Regarding the gain in throughput for Cases 2 and 3 with respect to Case 1, for E b /N 0 values lower than −2 dB, the best option is to use Case 2; for an E b /N 0 ranging between −2 and 4 dB, both Cases 2 and 3 provide practically the same throughput. Finally, for E b /N 0 larger than 4 dB, the best strategy is to encapsulate the NB-IoT frames into the DVB-S2(X) ones (Case 3). Finally, Figures 17-19 show the BLER of the three relay cases under study.

Conclusions and Futures Extensions
This section presents the main conclusions of the paper and possible future research to extend this work.

Concluding Remarks
Future mobile systems over satellite may suffer a bottleneck in the feeder link when they have to provide connectivity to OTA applications on a global scale, such as software and firmware updates for autonomous driving. To tackle this issue, the use of optical wireless technologies has been seriously considered to implement the feeder link and satisfy the high data rate demand that will be required. In this regard, this paper presented closedform formulas that estimate the impairments that are introduced in the optical uplink transmission in terms of the wind speed, wavelength of the optical signal, beam width and telescope aperture, refractive index structure of the atmosphere, and azimuth angle of transmission, among others. Furthermore, the use of time-packing has been also considered to increase the spectral efficiency of the optical feeder link even further. Three different relaying strategies have been analyzed for the GEO satellite, namely: (1) Detect-and-Forward with NB-IoT; (2) Decode-and-Forward with NB-IoT; and (3) Decode-and-Forward with NB-IoT/DVB-S2(X). From the simulation results that were obtained, it was possible to conclude that the wind speed has little effect on the end-to-end for mean SNR values larger than 20 dB in the optical feeder link. For this reason, the performance of the different relaying strategies has been studied in further detail for an optical feeder link at 15 dB (i.e., up to 10 dB atmospheric loss in the optical link budget).
For fair comparisons among all the cases under study, a similar code rate has been included. For Detect-and-Forward and Decode-and-Forward relaying architectures with NB-IoT, these code rates correspond to the ones defined in the NB-IoT standard. In contrast, for Decode-and-Forward with NB-IoT/DVB-S2(X), the equivalent code rate is equal to the product between the code rates of the DVB-S2(X) and NB-IoT frames, which can be obtained by using two MCS with high code rates. After evaluating these relaying configurations, it was possible to conclude that the Detect-and-Forward with NB-IoT achieved its larger throughput when using the lowest NB-IoT code rate and the largest overlapping factor. In contrast, the relaying architecture based on Decode-and-Forward with NB-IoT obtained its largest throughput when using intermediate code rates and intermediate overlapping factors. Finally, for the Decode-and-Forward with NB-IoT/DVB-S2(X) relaying architecture, the largest throughput was obtained when using an intermediate code rate for DVB-S2(X), with a variable code rate for NB-IoT that was adjusted according to the E b /N 0 in the radio access links. In this latter case, the larger the E b /N 0 , the higher the code rate that NB-IoT should be used. Regarding the most convenient overlapping factor, the highest throughput was observed when using the largest overlapping factors.
All the proposed relaying architectures were evaluated assuming time-packed signalling in the optical feeder link. The ISI introduced by time-packing was mitigated using an adaptive MMSE equalizer. From simulations, it was concluded that the use of timepacking increases the throughput in all proposed regenerative strategies under study, when compared to the cases in which time-packing was not used. However, this gain was depended on the E b /N 0 of the access link. For Detect-and-Forward with NB-IoT, this gain varied from 10% to 50-65% for low and medium E b /N 0 in the access link, respectively. For Decode-and-Forward with NB-IoT, the gain in throughput observed with time-packed signalling with respect to no-time packed going from 45-65% for low E b /N 0 of the access link to 10% for medium-to-large E b /N 0 values. Therefore, the decoding of the information permits to achieve higher gains at lower E b /N 0 . However, at larger E b /N 0 , the behavior of time-packed and no-time-packed signals is quite high and, as consequence, the improvement in throughput when comparing both schemes is reduced. Finally, for the Decode-and-Forward with NB-IoT/DVB-S2(X), the gain in throughput varies between 50-65% for medium E b /N 0 values in the access link to 45% for high E b /N 0 ones. In this scenario, the large error correction capability of LDPC codes enables one to remove the residual ISI that the equalizer is not able to eliminate. Consequently, it can deal successfully with the strong ISI that time-packing with large overlapping factors introduce.
Next, when comparing the aggregate throughput that the three relaying strategies offer, it can be concluded that the ones based on the Decode-and-Forward architecture provide the largest values. Specifically, the Decode-and-Forward with NB-IoT offers a better throughput for low E b /N 0 in the access link, whereas the Decode-and-Forward with NB-IoT/DVB-S2(X) provides a larger throughput for medium-to-high E b /N 0 values. Note that Decode-and-Forward with NB-IoT regenerates the NB-IoT signal in the satellite, whereas the Decode-and-Forward with NB-IoT/DVB-S2(X) protects the NB-IoT symbols by introducing the code rate of DVB-S2(X) in the optical feeder link. In this situation, when the E b /N 0 of the access link was low, the LDPC decoder of DVB-S2(X) did not converge and, due to that, it introduced errors in the NB-IoT frames that were decoded on-board the satellite for re-transmission in the radio access link. In this situation, these NB-IoT frames were protected with a code rate larger than their equivalent for the other cases; therefore, it can be concluded that at a low E b /N 0 in the radio access link, the Case 3 relaying strategy is not the optimum one. However, when the E b /N 0 of the access increased, the LDPC decoder started to converge and, due to that, it managed to remove the erroneous bits that were introduced by the optical feeder link in the NB-IoT frames that were encapsulated in the DVB-S2(X) signal. Here, as the NB-IoT frames had a higher code rate than the code rates that were used for other configurations, it offered the largest throughput.
Finally, it is remarked that the link layer should be prepared for adjusting dynamically the transmission architecture according to the E b /N 0 of the access link. Therefore, a control channel from the NB-IoT terminals to the gateway would be necessary, such that the estimated E b /N 0 of the access link can be known in advance when selecting the MCS of the NB-IoT frame, enabling to improve the end-to-end throughput. In all cases, the IoT terminal would receive the data in NB-IoT signalling format. Note that all required modifications of the Decode-and-Forward of NB-IoT with and without DVB-S2(x) encoding would be transparent to the IoT terminal, since the time-packing signalling, NB-IoT regeneration, and the DVB-S2(X) encapsulation would be performed at the satellite feeder link.

Future Extensions
After presenting the main conclusions, we introduce possible future research extensions of this work. Specifically, the following ones are considered: (i) Adapt the optical channel to LEO and MEO scenarios, in order to evaluate the performance of the proposed relaying strategies at different satellite orbits; (ii) evaluate the benefits of using time-packing with optical feeder link in the reverse flow of information, from the IoT terminals to the satellite gateway; and (iii) study the viability of using time-packing techniques as a potential physical layer security scheme. In the first case, the aim would be to assess the benefits that LEO and MEO satellites could reap from the proposed relaying architectures. For both LEOs and MEOs, the slant range and elevation angle of the feeder link changes continuously with the time and, as consequence, the optical channel modeling has to be adjusted accordingly to consider these effects. The complexity of the proposed schemes should be also analyzed in detail, since LEO and MEO satellites have more energy constraints than GEO ones due to the Earth/Moon possibly blocking the sunlight that reaches their solar panels. In the second case, the reverse link should be evaluated by replacing the time-packing scheme with a frequency-packing one. It is well known that NB-IoT uses SC-FDMA in the uplink and, as consequence, solutions that increase the spectral efficiency by overlapping the subcarriers should be better considered. As a side effect, SC-FDMA can compensate in part the increase of PAPR that overlapping in the time and frequency domains introduce. Finally, in the third case, the ISI that the process of shrinking the pulses/subcarriers introduces could be used as an artificial noise signal, which can mask the desired information from potential eavesdroppers.
Author Contributions: J.B. and A.A.D. conceived, designed, and analyzed the experiments; J.B. was oriented to DVB-S2(X) and NB-IoT signalling, time-packing generation, and its suppression, A.A.D. focused on modeling, and J.B. and A.A.D. wrote the paper. All authors have read and agreed to the published version of the manuscript.
Funding: This research was funded by Agencia Estatal de Investigación (AEI) and the European Regional Development Fund (FEDER) with grant number TERESA-TEC2017-90093-C3-1-R (AEI/FEDER, UE), and has been based upon work from COST Action CA19111 NEWFOCUS, supported by COST (European Cooperation in Science and Technology).

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Acknowledgments:
The authors would like to express their own deepest gratitude to Càndid Rèig for his support in publishing this paper in this Special Issue.

Conflicts of Interest:
The authors declare that there is no personal circumstances or interest that may be perceived as inappropriately influencing the representation or interpretation of reported research results. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.