Tracking-and Scintillation-Aware Channel Model for GEO Satellite to Land Mobile Terminals at Ku-Band

Recent advances in satellite to land mobile terminal services and technologies, which utilize high frequencies with directional antennas, have made the design of an appropriate model for land mobile satellite (LMS) channels a necessity. This paper presents LMS channel model at Ku-band with features that enhance accuracy, comprehensiveness, and reliability. The effect of satellite tracking loss at different mobile terminal speeds is considered for directional mobile antenna systems, a reliable tropospheric scintillation model for an LMS scenario at tropical and temperate regions is presented, and finally a new quality indicator module for different modulation and coding schemes is included. The proposed extended LMS channel (ELMSC)model is designed based on actual experimental measurements and can be applied to narrowand wide-band signals at different regions and at different speeds and multichannel states. The proposed model exhibits lower root mean square error (RMSE) and significant performance observation compared with the conventional model in terms of the signal fluctuations, fade depth, signal-to-noise ratio (SNR), and quality indicators accompanied for several transmission schemes.


Introduction
The recent applications and services in the broadcasting satellite communication to land mobile terminal have increased the demand for higher data rate and thus higher transmission frequency.Therefore, land mobile satellite (LMS) channel estimation at high frequency has become a necessity to develop efficient adaptive transmission models and techniques as solutions for channel impairments.Considerable interest has been directed recently toward LMS communication at Ku-band [1], which is the main concern of this study, owing to the high number of existing transponders and the large amount of available bandwidth.Nevertheless, Ku-band receivers require a high-gain directional antenna [1].A new generation of statistical LMS channel models capable of producing time series generative models was initially introduced for lower frequencies, particularly for L and S bands.These models have replaced the ones that solely provide cumulative distribution functions (CDFs), which are insufficient for appropriate investigations [2].
The LMS channel condition at Ku-band depends on mobility impairments and tropospheric scintillation.The latter, which causes rapid fluctuations in satellite signal level, occurs due to the irregularities in radio refractivity as the wave travels along different medium densities in the troposphere [3,4].Several researchers [1,2,[5][6][7][8][9][10] have created well-designed LMS channel models.The accuracy of the estimated models has increased notably over time through the addition of several features for approaching the real-world environment along with recent LMS technologies and services.This condition has motivated researchers, as well as the authors of this paper, to design more reliable LMS channel models.
However, existing channel models do not consider tropospheric scintillation under nonrainy conditions, which significantly affect the signal performance at Ku-band, particularly at high humidity regions, such as tropical areas [4].Moreover, 2 International Journal of Antennas and Propagation these models do not consider the impairments caused by different vehicle speeds at Ku-band for systems utilizing mobile directional antenna.Therefore, a comprehensive LMS channel modeling that considers these significant impairments is needed to be designed with lower root mean square error (RMSE).
This paper presents a channel model for satellite to land mobile terminals based on actual measured channel conditions.The proposed model is referred to as extended LMS channel (ELMSC) model.The term "extended" refers to four new features included in the model design.Firstly, the improvement is based on actual signal measurements to enhance the accuracy and reliability of the previously developed multistate LMS model at Ku-band.Secondly, the impairment attributed to the variable vehicular speeds is modeled concerning the clear line-of-sight (LOS) and shadowing scenarios.Thirdly, an LMS tropospheric scintillation model is developed.Lastly, a quality indicator module is developed and added to the ELMSC model.The model considers tropospheric scintillation and vehicular environments as well as its application to narrow-and wide-band signals worldwide because the LMS environment varies with respect to different regions in the world, particularly in temperate and tropical regions.
The LMS channel characteristics and impairments regarding the mobility and tropospheric scintillation effects are presented in the next section.The experimental setup is enlightened in Section 3 taking into account the mobility impairments for different states accompanied with tropospheric scintillation.The proposed ELMSC model is presented in Section 4. Section 5 provides a discussion of the results obtained from the measurements and proposed model, along with comparisons made concerning the existing models.The conclusion is presented in Section 6.

LMS Channel Characterizations and Impairments
Transmission parameters considerably affect the channel characteristics, the amount of signal attenuation, and the quality of service (QoS) of the LMS signal.In particular, the transmitted frequency, bandwidth, and elevation angle exert effects that should be identified by satellite system designers prior to the design process.The amount of signal attenuation is directly proportional to the carrier frequency and inversely proportional to the elevation angle; thus, the signal quality indication metrics, such as the signal-to-noise ratio (SNR) and energy of bit to noise ratio (  /  ), are affected.For signal frequencies below 3 GHz, ionospheric scintillation has a paramount effect on signal performance.This effect begins to disappear as the frequency increases above the said value [11].The bandwidth effect will be explained in Section 2.1.
The mobility impairments and channel states are presented in Sections 2.2 and 2.3, respectively.The tropospheric scintillation effect is discussed in Section 2.4.

Bandwidth and Directional
Antenna.The bandwidth of the signal is considered narrow-band when it is smaller than the coherence bandwidth; otherwise, it is wide-band.The coherence bandwidth, which is inversely proportional to the delay spread, is typically in the range of 7 MHz to 11 MHz at Lband and approximately 30 MHz at EHF band [1,12].Under the narrow-band condition, which is the most probable condition in the LMS channel [13], the channel causes signal amplitude variations, whereas the time dispersion (time delay of the received echoes) is insignificant.Under the wide-band condition, the channel not only causes variations in signal amplitude but also experiences significant time dispersion in the received signal.This condition will cause distortion effects, such as frequency selectivity or intersymbol interference; this distortion effect is attributed to the arrival of the signal echoes at the receiver at different excess delays with respect to the direct signal [2].Practically, the radiation pattern of a directional antenna, which is commonly used in LMS scenarios at Ku-band, filters out the echoes with significant delays; therefore, narrowband models are typically employed for Ku-band signals for narrow-and wide-band signals [1,12,14].

Mobility.
Advanced technologies for satellite communication services have resulted in a significant increase in mobile satellite terminals.This condition gave rise to the demand for the estimation of the mobility effect on the satellite link channel.Mobility impairments are typically produced by the shadowing or blockage of signal energy, multipath, and antenna tracking error.Shadowing is usually caused by roadside trees, whereas blockage is normally caused by surrounding tall buildings or bridges.The multipath effect, which can be modeled by Rayleigh distribution [15], is primarily attributed to nearby buildings in urban areas or mountains.
These major factors affect the performance and quality of the received signal and have motivated LMS system designers to build well-formulated mobility channel models to predict the signal performance of mobile scenarios.For accurate modeling, the model design should be based on actual channel condition measurements.Thus, a number of experiments are conducted to compare the measured results with that of our estimated model (discussed in Section 3).The received signal strength in a tropical region is measured in consideration of the scenarios of signal quality degradation using a mobile receiver system mounted at the roof of a car.To the best of our knowledge, this is the first experiment conducted to determine the mobility effects on satellite link channel conditions in tropical regions.

Channel States and Conventional Model.
According to the type of impairments encountered during the movements, LMS channel models can be categorized into single-state and multistate models [16,17].The single-state model describes the signal level distribution of a specific propagation environment over time, such as clear LOS, direct signal blockage, or shadowing states.Typically, in a single-state LMS channel, the unblocked LOS signal is modeled by Rice distribution.The Rice factor () is defined as the ratio of carrier to multipath normalized average power.Consider where  is the amplitude of the direct signal and  is the standard deviation.If the signal is totally blocked, it is usually represented by Rayleigh distribution.Meanwhile, several researchers have established singlestate channel models that focus on the shadowing condition; they defined this state as a combination of two distributions.The most used models for the LMS scenario are Suzuki and Loo models, which have gained global attention [13].Loo model assumes that, in the case of shadowing, the Riciandistributed direct signal amplitude varies according to a log-normal distribution, and the received complex envelope consists of the sum of the two phasors.The model considers the change in the phase caused by the Doppler effect, which is practically detected by the omnidirectional antenna.Loo model is applied to multistate LMS channel models at L and S bands [10,16,18].Suzuki model is a combination of Rayleigh and log-normal distributions and is conventionally applied to the modeling of multistate LMS channels at Ku-band using directional antenna [1].
However, in a real-world environment, two or more LMS channel states may occur during normal satellite terminal movement.Therefore, the signal performance exhibits several significant transformations in its characteristics over time unlike the terrestrial scenario whose behavior is usually described by a single distribution [2].In such case, multistate statistical models should be used to design LMS channel models [2,19].
Based on data obtained from measurement campaigns in Europe, Lutz et al. [8] designed a two-state model that estimates the LMS signal performance under clear sky (Rician distribution) and signal blockage (Rayleigh distribution) conditions.A high-state channel model is typically based on the Markov chain approach [15,19].The Markov chain is a stochastic model in which a system takes on discrete states.This approach is commonly used to model the variations in signal performance attributed to shadowing or blockage effects [10].The correlation between rain impairment and mobility impairment at Ku-band is presented in [9]; they modified the Rice factor to include the rainy sky scenario.
The first stochastic multistate model based on received signal measurements at Ku-band was supported by the European Space Agency (ESA) [1].A three-state channel model that reflects the proper approach for the real-world environment, including clear sky, shadowing, and blockage scenarios with multipath effect, was employed.The model was built based on the measurement campaign in Germany with a constant high speed of 120 Km/h in a highway environment and lower than 60 Km/h in other environments.However, the receiver antenna was mounted on a mechanically steerable platform during measurements; this condition may cause pointing misalignments and hence a high tracking accuracy error.
The International Telecommunication Union (ITU) [20] recommended the use of a three-state model for statisticalapproach LMS channel modeling.Moreover, [2] stated that the three-state model provides a tradeoff between complexity and accuracy.For the aforementioned reasons, a three-state channel is utilized in the proposed model.Thus, the model is compared with the conventional model [1] that also implemented the three-state scenarios at Ku-band, as explained in the previous paragraph.Nevertheless, other states can be added besides the aforementioned three states; however, this condition represents only special case modeling, such as that presented in [10,18], where a fourth state was added.The former described the fourth state as signal unavailability attributed to the failure of the satellite tracking process for very low elevation angles (lower than 10 ∘ ), whereas the fourth state of the latter represented multisatellite diversity (angle diversity).

Tropospheric Scintillation.
Tropospheric scintillation occurs inherently as the wave travels along different medium densities in the troposphere and causes rapid fluctuations in satellite signal level due to irregularities in radio refractivity [3,21].The standard deviation of fluctuations in tropical regions differs from that in temperate regions owing to their different related weather parameters, particularly the temperature and humidity of the medium.Accordingly, this condition results in different satellite channel assessments and signal quality performance.However, tropical regions are warmer and have higher relative humidity (  ) than temperate regions for most time of the year [22].ITU [23] claimed that the fade depth of tropospheric scintillation may reach a few dBs based on its probability of exceedance time.

Experimental Setup
To achieve accurate channel estimation, the model design should be based on the measured values of the LMS channel conditions.Several experiments were therefore conducted to measure the received signal level at the receiver (installed in a vehicle) at different mobility impairments.The impairments that occur typically in everyday life were considered (Figure 1), namely, effects attributed to different vehicle speeds (up to 150 Km/h), shadowing by roadside trees, and blockage by crossroad bridges.
The experimental measurements obtained in the southernmost region of Kuala Lumpur in Malaysia are shown in Figure 2. The measurements are represented by the three states listed in Table 1.
A mobile antenna system (upper right image of Figure 2) was used in the measurements setup mounted at the roof of a car (upper left image of Figure 2).The system consists of a 0.44 m mobile autosteerable antenna pointed at MEASAT 3/3A broadcasting satellite (91.5 ∘ E) at 76.5 ∘ elevation angle.A high tracking rate (75 ∘ /sec) antenna was employed to reduce the tracking error as much as possible.The antenna system also includes a built-in motor and gyroscope to point to the satellite automatically.The detection, filtering, and amplification processes were performed by the antenna system

State 2 Shadowing
The vehicle moved at constant speed (40 km/h) under link shadowing by tall trees (approximately 12 m tall and 7 m away from the vehicle on the average) on the roadsides.

State 3 Blockage
The vehicle moved at constant speed (40 km/h) under blockage by road bridges.

Scenari o 1 S c e n a r i o 3
S c e n a r i o 2  As the signals pass through the mobile antenna system, they are split into two directions as shown in Figure 3; the first direction is toward the decoder, and the second is toward the spectrum analyzer for data evaluation and analysis.
Meanwhile, the relative humidity was obtained on the day of the measurement campaign from the Malaysian Meteorological Department (MMD).The satellite and receiver systems parameters are listed in Table 2.
Typically, when the directional antenna is used, the reflections from buildings come from angles that are different from those in the satellite direction.This condition causes the reflections from the buildings to be rarely measured [1].The measurement campaign was conducted in highway and suburban areas as shown in Figure 2. The areas contain many crossroad bridges, tunnels, and road toll counters (whose sheds cause signal blockages) as well as roadside trees that cause shadowing because of their dense leaves and branches.During the continuous measurements, 837 data recordings were obtained at a clear sky state at various vehicle speeds, whereas 21 and 17 recordings were obtained during signal shadowing and blockage, respectively.The probability of accordance of each state is obtained depending on the aforementioned numbers of recordings and shown in Table 3.
The large number of measured samples is used in the comparison with the measured data to improve the accuracy of the ELMSC model.However, the probability of each state can be assumed in the design of the LMS channel model according to the environment of the specified region of interest.

ELMSC Model
Fading channel model and characterization gain the outmost importance for systems utilizing LMS link adaptive techniques and channel-aware modeling performance evaluation and assessment.The proposed model consists of several sections as shown in Figure 4.These sections include the   In an actual LMS environment, channel output () can be defined as where () is the transmitted signal, () is the LMS channel fading, and () is the signal noise.The proposed ELMSC model imposes a comprehensive impairments estimation effect which includes the impairments of mobility () and tropospheric scintillation ().Consider where   () is the total mobility impairments.The mobility impairments of the LMS channel include the attenuation attributed to the satellite terminal speed, shadowing or blockage attributed to physical obstacles in the link, and the multipath effect.In practice, both impairments are modeled separately to produce the time series for correlated scintillation and mobility impairments.Afterwards, their values are correlated according to the mobility state environment.The International Journal of Antennas and Propagation details of the proposed model blocks are explained in the subsequent subsections.

Multipath Generator.
The diffused multipath fading model can be characterized by Rayleigh distribution [15].The envelope is a result of the reflected signals from number of paths  as illustrated in (4).Consider where   () is the reflected power of the th signals and Γ  () is a coefficient that depends on angular Doppler frequency   as well as angle of arrival  and phase  as expressed in (5).Consider Therefore, (4) can be rewritten as ( 6) or (7).Consider where   () and   () represent the inphase and quadrature components of the Rayleigh distribution, respectively.These two components, as predicted by the central limit theorem, can be produced by first generating two, zero-mean white Gaussian distributions [13,24] as shown in Figure 4.The standard deviation  varies with respect to the clear sky, shadowing, and blockage states because the fluctuation level attributed to the multipath effect differs with respect to each state [13]. reaches the maximum in the blockage state and minimum in the state of clear LOS.Therefore,  is included in the multistate model (explained in Section 4.2).Power spectral shaping is achieved by coloring "filtering" the resultant complex signal.The inphase and quadrature trails are filtered with a Doppler filter to shape "reproduce" the Doppler effects for the fading channel attributed to the motion of the terminal.The filtering process is necessary in simulating the channel.Without this process, the resultant signal variations would be unnaturally fast for a given satellite terminal speed and would thus produce error bursts and distribution of outage durations [25].Therefore, a Butterworth filter is used in the proposed model; this filter is commonly used in LMS channel modeling because it provides a more realistic approach [13,24,26,27].The filtering process is designed based on [13,27]  The cumulative distribution of the resultant signal should be the weighted combination of the three states [10,23].
Therefore, the probability behavior for the three states is expressed in (8).Consider where  is the sample number from the total number of samples .  is the probability of occurrence of the th state ( 1 +  2 +  3 = 1).The probabilities were obtained from the experimental measurements as mentioned in Section 3 and were used to generate random state series.The total time (or number of samples) of being in a specific state depends on the probability of occurrence.The output discrete signal is a result of switching between the generated signal samples of the aforementioned states according to the state series.(1), where  represents the ratio between the direct signal level and the noise floor level.Therefore, our assumption on the Rice factor is based on the actual measurements of the signal.Typically, the noise floor level is assumed based on the measurement setup [1,2].Moreover, to estimate the signal losses that may occur because of the antenna tracking error attributed to terminal movements, a new formula is proposed (see (9)) to estimate the signal losses for satellite terminal movements at different speeds in dB based on the actual measurement explained in Section 3. Consider where  is the satellite terminal speed in km/h and the constant values   ,   , and   are −4.667× 10 −5 , 1.468 × 10 −2 , and −1.364 × 10 −3 , respectively.These values were obtained with nonlinear curve-fitting optimization technique; the proposed formula provides accurate mobility impairment estimation with a RMSE of 0.02407.Terminal speed is used as an input to the module.This input series can be generated using sinusoidal distribution random numbers with a minimum value of 0. The maximum value can be assumed based on the nature of the road of interest and its highest permitted speed.Moreover, the probability of speed values can be obtained through statistical observations.

Shadowed Signal Model.
The shadowing (state 2) attributed to roadside trees causes additional attenuation to the received signal and thus results in a slow variation in the direct signal simultaneous with the diffuse multipath components.Typically, the shadowed signal at Ku-band is modeled with the Suzuki model [1].This model involves the combination of log-normal and Rayleigh distributions.
Normal distribution signal   is generated with its specific mean  and standard deviation Σ obtained from [1] in dB before being converted to a log-normal distribution using the relation 10   /20 .The resultant signal passes through two stages: rate conversation (to identify the slope of the envelop and number of the added samples based on   of the fast multipath fluctuations) and interpolation (adding samples between two different leveled samples corresponding to the rate conversation and then convolution for envelope shaping of the interpolated samples).In summary, these steps were implemented to add an interpolated sloped envelop sample to approach the number of multipath samples and to avoid sharp multipath fluctuations.For periods with no state transitions, this rate is utilized to generate rapid variations with a low sampling period to account for the Doppler spread bandwidth [13].
The multipath fluctuations of this state were generated as discussed in Section 4.1.The resultant combination of the two distributions is then obtained to reflect the performance of the signal in the shadowing environment.The impairments that occur primarily because of the antenna tracking error due to the terminal movements, expressed in (9), affect the shadowed signal as well because the direct signal component still exists despite being shadowed.Therefore, in the case of a partially shadowed environment (state 2), the   module is also involved in the simulated shadowed signal.

Blocked Signal Model.
LMS signal blockage (state 3) is typically caused by bridges across the highway, sheds of the road tolls counters, tunnels, or tall buildings in urban areas.The received signal power vanishes, and the fade depth can be expressed by the noise floor level below the LOS.The LMS signal blockage model is typically designed with Rayleigh distribution with specific standard deviation  BL as discussed earlier.Therefore,  BL is applied to the resultant fluctuated signal from the multipath model discussed in Section 4.1.Considering that the direct signal is totally blocked in this state, the signal would be unaffected by the   module.The status of the current state is obtained from state series block as shown in Figure 4 to switch the   module on during the clear LOS and shadowing states and off during the blockage state.

Tropospheric Scintillation Model.
Tropospheric scintillation depends on its standard deviation (  ) [21,23].The modeling of tropospheric scintillation effect begins with the generation of a zero-mean white Gaussian distribution and with standard deviation conventionally set to 1 [23,28].However,   depends mainly on the relative humidity and the temperature [23], which is different in tropical regions compared with temperate regions according to their weather parameters.For temperate regions, Cioni et al. [28] assumed that the value of   is 0.066.However, it is declared by [22] that the standard deviation of scintillation in tropical regions is larger than that in temperate regions.Therefore,   is adjusted to 0.117 to increase the accuracy of the estimated scintillation standard deviation.The value of   is estimated based on two steps.First, the   value was obtained from the MMD as discussed in Section 3. Second, the mean value is used to estimate   according to (10) [23].Consider where  is the frequency in GHz,  is the angle of elevation, and  is the antenna averaging factor that depends mainly on the frequency and antenna parameters. 0 can be obtained in dB from the wet term of radio refractivity  wet which, in turn, depends on the value of temperature  and   [29].Consider where   is the temperature-dependent saturation vapor pressure that can be calculated from [29].The mean value of   measured by MMD is 82.375%.With   , the tropospheric scintillation fade depth  scint is estimated using the ITU model [23] for different probabilities of occurrence ().
The power spectral density (PSD) of the fading is time variant and can be approximated by low-pass filter (LPF) [28].The LPF in [28] is used, and the resultant signal is then uniformly distributed and combined with  scint to obtain the estimated signal series of tropospheric scintillation.

Complementary Blocks and Quality
Indicator.The total performance of the correlated impairments is a combination of the LMS multistate mobility and tropospheric scintillation models as declared in (2).However, as the signal passes through the channel, several other constraints cause additional losses, such as system loss and free space loss (FSL).The latter depends on the transmitted frequency and the distance difference between the satellite and the mobile terminal as revealed in (12).Consider FSL = 20 log ( 4  ) .
Consequently, the estimation of the received SNR, which is the received signal power   to noise [  = 10 log(  TBW)] ratio is included in the design.Equation ( 13) [30,31] is utilized to estimate the received power   in dBW.Consider where EIRP is the effective isotropic radiated power in dBW that depends on the transmitted power and the transmitter antenna gain;   is the receiver gain;   is the total system loss; and  LMS is the total LMS channel loss that has not been considered in previous studies, which represents the total fading losses that can be obtained from the ELMSC model.After the insertion of the channel fading effect to the transmitted signal, additive white Gaussian noise (AWGN) is added.
The atmospheric impairments have a negative effect on the data after being demodulated in the receiver.The effect appears, for example, as a decrease in   /  or packet error International Journal of Antennas and Propagation rate.These ratios can be utilized to generate the proposed indexing for optimal physical or link layer adaptations for the current instantaneous transmission time interval (TTI).Equation ( 14) is used to obtain   /  (in dB) for widely used modulation schemes and code rates (MODCOD) in satellite communications.Consider where  = log 2 () ×   ,  is the modulation order,   is the code rate,   /  is the symbol energy to noise ratio, BW is the bandwidth, and  is the symbol rate.For fixed  and BW,   /  is directly proportional to SNR and inversely proportional to  and   .Therefore, after substituting (( 12)-( 14)) to ( 15), the estimated   /  is obtained as shown in (16).Consider where / is the figure of merit of the receiver antenna and   is the Boltzmann constant.  /  is used as a quality indicator for the ELMSC model.

Results and Discussion
For realistic channel estimation, the estimation was based on the actual measurements of signal performance under dynamic weather conditions.The model design results are organized based on the arrangements of the blocks in Figure 4.
The generated signal (Figure 5(a)) exhibits rapid variation at a transmission frequency of 11.6 GHz. Figure 5(b) shows that the distribution of the resultant signal did not change after filtering.
The state series shown in Figure 5(c) comprises the three states distributed randomly with respect to the probability of occurrence of each state.Figure 6 provides a zoomed in version of Figure 5 to recognize the generated samples clearly.The succeeding results are presented based on the zoomed in version to identify clearly the effect of each block in the ELMSC model.For accurate comparison with other models and better modeled signal resolution, the comparison was based on extended samples as in Figure 5.
The resultant multipath signal passes through the multistate model with respect to the state series.If the current state is 1, then the multipath signal samples will pass through the clear LOS model, which is composed of three stages.First, the samples are multiplied by  LOS to obtain the Rician multipath signal.Second, the direct signal component () is added.Lastly,   is added at different speeds.The resultant clear LOS fading is shown in Figure 7.
The direct signal component () is calculated based on (1), where Rice factor  = 11.5 dB.The value of  is assumed based on the measurements discussed in Section 3.
The channel parameters are listed in Table 4.
During the measurement campaign, vehicle speed was recorded simultaneously with the measurement data to account for the additional losses from the movement of the satellite terminal at different speeds.These losses are mainly due to the pointing misalignment caused by the antenna  tracking error.The resultant average measurement data of signal attenuation with respect to satellite terminal speed is shown in Figure 8.With respect to terminal speed, the average values of several measured signal samples have been obtained.The speed of approximately 40, 80, 100, 120, 140, and 150 km/h has been considered in the measurement campaign as well as the stationary (zero speed) case.The measurements show that the signal is further attenuated by 0.6 dB when the terminal moves at 60 km/h.The attenuation slope changes when the speed reaches 80 km/h.The attenuation level increases more slightly than when speed is below 80 km/h.The losses reach approximately 1 dB at 100 Km/h, 1.15 dB at 140 Km/h, and 1.18 dB at 150 Km/h.However, as discussed in Section 4.2 and expressed in Figure 8, the proposed equation ( 9) provides an accurate estimation of signal performance at different terminal speeds with RMSE equal to 0.02407.To the best of our knowledge, this model is the first to be designed for the estimation of LMS channel impairments at different terminal speeds based on measurement data at the Ku-band.
Based on the experimental measurements in the clear LOS scenario, the signal model at state 1 provides less RMSE estimation than the conventional model proposed earlier.This finding can be clearly proven by comparing the CDF and PDF of the measured SNR and the conventional model with the proposed model.After adding the effect of FSL from (12) to obtain the channel fade and AWGN to the transmitted signal carrier, the CDF and PDF of the estimated received power were obtained as shown in Figure 9.
To identify the effect of the proposed multistate part of the ELMSC model with respect to the clear LOS measured signal and conventional model, Figure 9 shows the CDF and PDF of the received SNR prior to the addition of the proposed tropospheric scintillation model.The effect of the proposed   module appeared as a change in the slope difference of CDF as shown in Figure 9(a).This leads to a decrease in the RMSE of the proposed model, whereas the change in the  factor shifted the CDF for few dBs towards the measured value, which approaches the CDF slope of the measured SNR compared with the conventional model.Until this point, the PDF performance of the proposed model exhibits a significant enhancement compared with the conventional model such that it approaches the PDF of the measured SNR and will be further enhanced if tropospheric scintillation is included.
In the case of state 2, the shadowing model is applied to the generated multipath signal in Figure 5(b).The model is a combination of two distributions.The slow variation trail is composed of three stages before the addition of rapid variation.Finally, the signal passes through the proposed   block as discussed in Section 4.2.The effect of these steps is shown in Figure 10.
The slow variation trail begins with the generation of random normal distribution samples.One of the objectives of the experimental measurements is to measure the average power absorption for each scenario.13 dB mean value is used as proposed in [32].
The distribution is then converted to a log-normal distribution before rate conversion and interpolation as discussed in Section 4.2 and shown in Figure 10(a).Figure 10(b) shows that the signal variation performance does not change after this step, except for the higher number of samples that accounts for the generated rapid variations.Figure 10(c) shows the simulated shadowing fade.
The signal level in a full blocking environment (state 3) usually reaches the noise floor level, as discussed in Section 4.2, with multipath fluctuations of specific standard deviation  BL that is typically greater than  LOS [13].Therefore, the resultant signal from the multipath generator is multiplied by the value of  BL , which is assumed according to International Journal of Antennas and Propagation [1].However, the simulated signal of the multistate model is a result of the switching process of the three aforementioned states.Switching control is based on the state series generated (shown in Figure 5(c)).
Another remarkable impairment to the LMS channel is tropospheric scintillation which gives rise to signal turbulence because of the rapid change in the refractive index of the medium.The change in the refractive index of the medium is because of the water particles in the troposphere.Therefore, the   percentage has a vital role in the signal performance and standard deviation of fluctuations.
is varied from one region to another according to their different weather parameters.Therefore, two different signals were modeled each with its specific region, one with lower standard deviation of fluctuation as proposed by [28] represents a temperate region and the other proposed in this study for higher humidity regions as in tropics based on the measured   discussed in Section 4.3.Figure 11 shows the  estimated tropospheric scintillation performance according to the proposed model.
Figure 11 shows the difference in the estimated signal samples between temperate and tropical regions.To identify the effect of tropospheric scintillation clearly in both regions, Figure 12 shows the CDF and PDF of the received SNR in temperate and tropical regions for clear LOS state after adding the tropospheric scintillation effect.
It is shown in Figure 12 that the distribution of the samples in tropical regions is different from the temperate regions according to its weather parameters especially the higher relative humidity and hence higher water particles which cause extra fluctuations compared with the temperate regions.Therefore, the fluctuation performance reflects a change in CDF and PDF with respect to the original distribution functions when tropospheric scintillation is excluded.It can be noticed that the peak value of the PDF of the estimated SNR in tropical region is moved to lower SNR because of the different weather parameters and approaches the peak value of the measured SNR (as seen later in Figure 14).As discussed in Section 2.4, fade depth depends on the probability of exceedance time ().Figure 13 shows The probability of exceedance time significantly affects the distribution function. is directly proportional to fade depth; the highest level reached  ≥ 1.
The effect of  was predicted based on the experimental measurements obtained by ITU [23].The proposed ELMSC model at tropical regions in clear LOS is compared with the measurement data (discussed in Section 3) and the conventional model used, as shown in Figure 14.
The CDF and PDF show that the ELMSC model provides less RMSE regarding the estimation of the channel compared with the conventional model.The PDF of the proposed model reached a peak value of approximately 0.3 among 8 dBs.This value approaches the measured value (difference = 0.02).By contrast, the conventional model reached the peak value of approximately 0.59 among 5 dBs, which significantly differs from the measured value (difference = 0.31).Table 5 shows the RMSE and the differences in variance and PDF's peak values of the ELMC and conventional model compared with the measured time series signal.The proposed fade for three states is shown in Figure 15.
The variance in fluctuation during state 1 is equal to 1.8133, which approaches the measured value of 1.6408 and is closer than the value of conventional variance (0.4508).The minimum duration for one state is set to four times the sampling period to reflect the real-world environment, such that, when the mobile terminal passes under trees, it usually stays in this environment for more than one sample period.The quality indicator included in the model is used to identify the channel quality in several widely used modulation schemes and code rates (MODCODs) that are used in video broadcasting protocols such as digital video broadcasting via satellite-second generation (DVB-S2) [33].Figure 16 shows the mean   /  value for the three states.These values were estimated at the input of the decoder.
/  is used as a quality indicator for the ELMSC model.The value is examined at different MODCODs.  /  determines the availability of the received signal as well as its quality.Availability is determined based on the value of a quality indicator.Therefore, to achieve signal availability for systems that use, for example, a minimum   /  value of 4 dB, the signal is available in the clear LOS scenario at all times.In the shadowing scenario, the signal is available when QPSK at   = 3/5 and 2/3 only, whereas no signal is available in the blockage state.Moreover, the values of   /  at different MODCODs and environments can determine the error rates in the data received.Therefore, the quality indicator is essential in the design of the channel model to specify the availability and error rates for the signal under different channel conditions and system parameters.

Conclusion
A model for LMS channel modeling at the Ku-band was presented in this study based on experimental measurements.The design considers the three typical channel states, namely, clear LOS, shadowing, and blockage states.The proposed model has new features and was proven to be comprehensive, reliable, and less RMSE than conventional model.The channel states were modeled and combined with the new module to account for the satellite tracking error attributed to the dynamic terminal speeds.The tropospheric scintillation impairments were presented and modeled for the LMS scenario at temperate and tropical regions.Finally, a quality indicator module was presented and included in the proposed ELMSC model.The measurement data were provided, and the time series synthesizer, PDF, and CDF of the proposed model were presented and compared with the conventional model.A fairly good agreement was observed between the proposed ELMSC model outputs and the measurements.The model can be applied to temperate and tropical regions for narrow-and wide-band signals at various modulation and coding schemes and different satellite terminal speeds.The model and its associated modules can be used to study the signal performance, availability, and error rates of different services, including communications, broadcast, and navigation, as well as to develop fade mitigation techniques for channel-aware strategies.

4. 2 .
Multistate Model.Three states, namely, clear LOS, shadowing, and blockage states, were implemented in the ELMSC model.The model introduces one of the extensions of the ELMSC model that represents the losses of satellite terminal attributed to movements at different speeds (  ); the new extension is applied to the clear LOS and shadowing states.

Figure 7 :
Figure 7: Fading in the clear LOS state.

Figure 8 :
Figure 8: Signal losses with respect to satellite terminal speed.

Figure 11 :
Figure 11: Modeled fluctuations caused by tropospheric scintillation: (a) temperate region and (b) tropical region LMS link scenarios.

Table 1 :
Measurements scenarios.Clear sky/different speedsThe vehicle moved at different speeds (from 0 km/h up to 150 km/h) during clear LOS.Measurements at zero speed (stationary) were obtained as well.
4.2.1.Clear LOS Signal Model.The signal received by the mobile satellite terminal moving at different speeds in a clear (unblocked) channel environment is usually modeled with Rician distribution as discussed in Section 2.3.The model includes the constant direct signal with multipath fluctuations.The latter, Rician multipath, is designed by first multiplying the resultant signal from the multipath model discussed in Section 4.1 with its specific standard deviation ( LOS ).The direct signal component () is then added to the Rician multipath signal. is extracted using

Table 5 :
Comparison with the measured data.