Advances in Statistical Channel Modeling for Wireless Communications
View this Special IssueResearch Article  Open Access
Marie Rieche, Alexander Ihlow, Daniel Arndt, Fernando PérezFontán, Giovanni Del Galdo, "Modeling of the Land Mobile Satellite Channel considering the Terminal’s Driving Direction", International Journal of Antennas and Propagation, vol. 2015, Article ID 372124, 21 pages, 2015. https://doi.org/10.1155/2015/372124
Modeling of the Land Mobile Satellite Channel considering the Terminal’s Driving Direction
Abstract
A precise characterization of the Land Mobile Satellite (LMS) channel, that is, the channel between a satellite and a mobile terminal, is of crucial importance while designing a satellitebased communication system. Stateoftheart statistical LMS channel models offer the advantage of requiring only a few input parameters, which include the environment type and the elevation angle of the satellite. However, the azimuth angle relative to the driving direction of the mobile terminal is usually ignored, as its proper modeling requires either an extensive measurement campaign or a significant effort from the user, as a precise geometrical description of the scenario is required. In this contribution we show that the impact of the driving direction on the channel statistics is not negligible and requires to be modeled explicitly. Moreover, we propose a statistical LMS channel model whose parameters are obtained via an imagebased state estimation method. The imagebased method is verified by a comparison with measured radio frequency signal levels. The proposed method allows obtaining a complete statistical description of the channel for arbitrary elevation and azimuth angles.
1. Introduction
Satellite communication services need to ensure a certain quality of the received signal. For stationary applications, this is a relatively undemanding task, as it is normally possible to provide a nearly ideal LineofSight (LoS) connection. However, this is not the case for moving terminals. In order to design appropriate systems for this specific application, a profound knowledge of the propagation channel, the socalled land mobile satellite (LMS) channel, becomes necessary.
The LMS channel is mainly impaired by the environment in the direct vicinity of the moving terminal. In fact, the LoS component may be completely obstructed by solid objects (e.g., a building) or may be more or less shadowed by vegetation depending on the density of the foliage and branches. Moreover, diffraction, reflection, and scattering of the transmitted signal generate multipath components.
Looking at the received signal level over traveled distance, as in Figure 1, it is possible to distinguish different effects: the LoS component (which might be blocked or not) gives rise to “very slow fading”, while in a partially obstructed condition, it gives rise to “slow fading”, for example, due to the vegetation. Furthermore, the multipath propagation introduces “fast fading”, which is superimposed to slow and very slow fading. Besides the impact on the signal level described above, the channel also causes temporal dispersion (induced by the different lengths of the propagation paths) and the Doppler shifts, as terminal, satellite, and possibly some objects in the environment are moving.
The LMS channel can be modeled with different approaches, where statistical channel models are frequently used and are able to generate timeseries of the received fading signal [1–3]. Their main advantage, compared to other approaches, is that they do not require many input parameters and offer the possibility to generate a large quantity of different channels. These models usually feature two different input parameters: the environment type (e.g., urban, suburban) and the elevation angle of the satellite. Normally, they do not consider the driving direction of the mobile terminal, that is, the azimuth angle of the satellite. Ignoring the driving direction in the model results in a timeseries of the signal fading, which is stationary with respect to the driving direction. This is intuitively an oversimplification for urban scenarios, where a satellite oriented in driving direction leads to a higher signal quality compared to a satellite perpendicular to the driving direction. This is exemplified in Figure 2.
(a) Lineofsight (LoS)
(b) Nonlineofsight (NLoS)
Traditionally, the internal model parameters are derived from statistical distributions, whose parametrization is based on Radio Frequency (RF) measurements. The measurement data utilized by the statistical LMS models, available in the literature, is gathered either from one or more satellites at specific azimuth positions or from a helicopter emulating the satellite. If measurements are carried out for several satellites at different azimuth positions, the average over the available satellite positions in azimuth is calculated. Consequently, the extracted model parameters reflect the flight/satellite conditions and the angle between the radio path and the mobile but do not explicitly consider it. To develop an LMS model which explicitly considers the driving direction, it is necessary to collect enough RF measurements to carry out a statistical analysis in dependence of the relative azimuth angle. For instance, by dividing the full circle () into intervals of each, one needs 36 times the measurement data to achieve an equivalent statistical significance.
In this contribution we extend the statistical LMS model for SBand proposed by Arndt et al. [5]. We consider the driving direction with a resolution of 1° over the whole azimuthal range of . To solve the challenging task of gathering the necessary large amount of RF measurement data, we follow an alternative imagebased approach [23], which allows us to extract the very slow fading of the received signal for arbitrary positions of the satellite, even for those which are not yet in orbit. The method enables the derivation of accurate parameters for different elevations and especially different driving directions with an angular resolution of 1°. Based on the high azimuthal resolution we found that the leftright symmetry assumed in [5] is an oversimplification, because the vehicle travels on one side of the street. We show the significant impact of considering the azimuthal satellite position relative to the driving direction, in contrast to neglecting the azimuthal position, on the quality of the received signal in terms of first order statistics.
The remainder of this contribution is organized as follows: Section 2 gives an overview of stateoftheart LMS channel models and closer analyzes those, feasible to reflect the effects of the driving direction. Section 3 illustrates the imagebased state estimation method, which allows the extraction of the very slow fading, Furthermore, the accuracy of this method is shown by comparing with simultaneously recorded RF signal levels. Section 4 introduces the main components of the LMS model. Section 5 presents the impact of the driving direction on the statistics of the satellite signal. Finally, Section 6 summarizes the main results of the contribution.
2. Modeling Approaches of the Land Mobile Satellite Channel
This section concisely reviews different principles of LMS channel modeling and explains the most common LMS channel models.
It is convenient to distinguish between four basic principles [24].(i)Empirical channel models describe the signal attenuation by fitting mathematical expressions to measurements, without including any physical interpretation. An overview of different empirical LMS channel models can be found in [25].(ii)Deterministic channel models describe the propagation effects for specific environments, while approximating Maxwell’s equations. Most of them are based on assumptions of a propagating electromagnetic wave and are socalled raybased. Raytracing approaches applied to the LMS channel are described in [26].(iii)Statistical channel models describe the channel in terms of statistical distributions, which may or may not have physical meaning. The parameters of these models are extracted from a representative number of measurements.(iv)Hybrid channel models are a combination of the approaches mentioned above. In the literature, a large variety of these models can be found, making them difficult to categorize. When deterministic and physical channel models are combined, they are socalled physicalstatistical models. A synthetic propagation environment is generated by considering deterministic/physical techniques, whereas statistical distributions are utilized for a description of the environment or the signal propagation.
Furthermore, we can identify more categories of available LMS channel models depending on other features.(i)One feature for categorizing LMS channel models is the number of reception states supported. The received satellite signal can no longer be assumed stationary for a longer time period. Different reception states are introduced to represent the blockwise stationarity of the received signal, which is caused by the slow varying environmental conditions (like LoS or NLoS due to, e.g., buildings, trees), while the terminal moves. Singlestate LMS channel models are those which only consider the slow and fast fading of the received satellite signal. Over a short time period these two components (slow and fast fading) are modeled by a stationary stochastic processes, for example, Loo [27–29]. Multistate LMS channel models differentiate two states (“good” and “bad”) [30], three states (“lineofsight,” “shadowed,” and “blocked”) [31] or higher (variable) number of states [32]. The sequence of reception states can either be generated by a Markov process or by physical/geometrical considerations of the environment (see Section 2.2).(ii)Besides the number of states we can also distinguish between LMS channel models for different number of satellites. In the literature, several singlesatellite models can be found [1, 2, 5, 30]. Available multisatellite models usually consider dualsatellite systems [3, 4, 6]. Lutz [3] already emphasizes the impact of the correlation between two satellite links, as it provides information about the visibility of satellites and the diversity gain. Propagation models which consider more than two satellites (e.g., [7, 33]) are rather rare in the literature, due to the more complex description of the satellites correlations.
In the following we deal more in detail with the purely statistical and the physicalstatistical approaches, as they provide the possibility to include the driving direction and offer less complexity, compared to purely deterministic models.
2.1. Statistical Land Mobile Satellite Channel Models
Stateoftheart purely statistical models [2], which produce timeseries of the envelope of the received signal, derive their parameters from RF measurements. They implicitly include the driving direction without evaluating it adequately. If the driving direction was considered explicitly in the channel modeling, the resulting statistics were more realistic for specific situations, as shown in Section 5. Therefore, it would be necessary to collect a sufficient amount of RF data to carry out a statistical analysis in dependence of the relative azimuth angle. In [1] the measurement data was conducted for relative azimuth angles of approximately , , and . Arndt et al. [5] utilize measurement data for different environments, elevation angles, and driving directions. The model consists of a state sequence generator implemented as semiMarkov chains and a Loogenerator simulating slow and fast fading effects within each state. The model parameters are derived from extensive measurements of RF signal levels in SBand for different environments, elevation angles, and driving directions. The intervals of driving directions with respect to the satellite are , , , and . This model assumes a symmetry of for the different driving directions and clusters the data accordingly, which is an oversimplification, as we show in Section 5.
2.2. PhysicalStatistical Land Mobile Satellite Channel Models
There are several hybrid LMS channel models available, which consider the driving direction by different geometrybased approaches. The “virtual city” approach represents a special case for the urban environment. Building heights, street widths building block lengths, and so forth are modeled following a statistical distribution. The output of these models are multiple state sequences, which allow to calculate the correlation factors between different satellite links of azimuth and elevation. One of the first physicalstatistical LMS models is introduced by Tzaras et al. [8]. The basis is a physical approach based on synthetic environments, where the heights of obstacles follow certain distributions (or are assumed to be constant). The mobile terminal can also be shifted away from the middle of the street. The result is presented in terms of an analytic expression based on diffraction theory, depending only on the azimuth angle separation.
Oestges et al. [9] develop a prediction of the statistical fade distribution based on a conditional Rice distribution. The parameters of the Rician distribution vary with the physical parameters and depend on the satellite position (elevation and azimuth) and on the environment (building height and street width). The relationship between the Rice parameters and the physical parameters are validated by raytracing, based on the Uniform Theory of Diffraction (UTD). Building height and street width are assumed to be lognormally distributed, based on measurements in London and Guildford. The results are shown for an azimuth angle of relative to the street axis and the vehicle was positioned in the middle of the street. Based on this model, the shadowing correlation coefficient for a dualsatellite system as a function of the azimuth separation is presented in [10]. The correlation model is based on the Rician dominant component of the physicalstatistical model explained in [9]. Here, the vehicle is positioned at one side of the road, which shows the effect of the righthand side and lefthand side driving, respectively. In [11], again the same geometrybased approach is used to model timeseries, combined with a series of roadside buildings according to measured distributions of geometrical parameters such as building height, street width, and height of the nearest building. This model considers the position of the car to be in the middle of the street.
There are several models based on a street mask concept, mask angle (MKA) or mask function (MKF), respectively [12–16]. These graphical representations of the azimuth and elevation angle indicate the visibility of satellites from the mobile terminal’s point of view for different urban scenarios (e.g., street canyons, street crossings). The input parameters are building height, terminal height, and street width, whereby the user is assumed to be in the middle of the street, therefore the correlation between two satellite links is symmetric. To generate the building heights and street widths statistical distributions are used.
A hybrid statisticdeterministic model [17] is based on a Markov chain, considering elevation and street orientation angles (assuming symmetry of ). The visibility of the satellite is also described by the MKA, when the terminal is located in the middle of the street. The resulting “on/off” timeseries is enhanced by multipath and diffraction effects based on a simple raytracing model (statistical distributions for building height and width are used).
A methodology, which uses also fisheye images to derive the state information (clear, shadowed, and blocked) and adding different distributions (i.e., Loo, Rayleigh, and Rice) to simulate the fast fading effects, can be found in [19–21, 34]. Meenan et al. [18] also use fisheye images and apply the approach to the correlated fourstate Lutz model [3]. However, none of these models explicitly takes into account the dependency on azimuth angle relative to the street orientation.
Table 1 gives an overview of the statistical and physicalstatistical models, which are concisely reviewed in this section. These models are frequently used in the literature and differ in their assumptions, constraints, and outputs. Therefore, the models are compared with respect to their modeling technique, the number of satellites, if they consider the driving direction of the terminal, as well as the correlation coefficient (corr). Further evaluation criteria are the cumulative distribution function (CDF) and the bad state probability of a satellite.

In summary, stateoftheart LMS channel models are either uncapable of describing the influence of the azimuth angle of the satellite due to the lack of data, or they necessitate significant effort from the user, as a precise description of the geometry of the environment is required.
3. State Parameter Extraction
As mentioned in the introduction, to carry out a statistical analysis for different driving directions it is necessary to collect an unfeasible amount of measurement data. To solve this problem, we rely on the imagebased method described in [35], which is reviewed in the next Section 3.1. This method is capable of estimating the very slow fading conditions of the LMS channel for an arbitrary position of the satellite using hemispheric images of the environment surrounding the mobile terminal. The very slow fading is modeled with a sequence of the reception states good and bad, corresponding to lineofsight (LoS) and nonlineofsight (NLoS), respectively [2]. Section 3.2 presents the measurement campaign realized to collect the necessary images for the model proposed here. Finally, Section 3.3 addresses the validation of the imagebased method by comparing its estimates with actual RF measurements.
3.1. ImageBased State Estimation Method
In the following, we review the imagebased method from which we derive the LMS channel reception states [35]. Hemispheric images of the environment are obtained from a fisheye camera pointing towards the sky. Such an image of the upper hemisphere is depicted in Figure 3(a). The image is classified in the region sky, where the satellite is visible and the region obstruction, where it is shadowed. The result is shown in Figure 3(b) as an overlay of the binary classification result and the original image. Finally, the classified fisheye image is converted into a rectangular image in landscape panoramic form yielding a resolution of one degree in elevation and azimuth (see Figure 3(c)). Knowing the time, location, and heading information of the vehicle we can determine the position of arbitrary satellites within the image given their azimuth and elevation with respect to a known coordinate system. From the position of the satellite within the panoramic image we can directly extract the reception state, where sky represents the good state and obstruction represents the bad state, respectively. Based on the extracted state sequence we can determine different state parameters, representing an essential part of the LMS channel model parameters, as explained in detail in Section 4.
(a) Exemplary hemispheric image
(b) Overlay of the original and the binary classification into sky (blue) and obstruction (red)
(c) The converted panoramic binary images are of size pixels, where the resolution is one degree in elevation and in azimuth. White represents the reception state good and black bad, respectively
3.2. Measurement Data Acquisition
To validate the proposed method we consider extensive measurements conducted in the context of the project Mobile satellite channeL with Angle DiversitY (MiLADY) [36]. During this campaign the received power levels of two geostationary (GEO) satellites of XM Satellite Radio and two of three highly elliptical orbit (HEO) satellites of Sirius Satellite Radio were recorded at the operational frequency of 2.3 GHz (SBand) with a sampling rate of 2.1 kHz. The precise measuring time and vehicle position was logged via GPS. In addition, a hemispheric camera was mounted on the van to document the environment, pointing towards the sky (5 frames per second with a resolution of 880 768 pixels of the fisheye region). The campaign was carried out along a distance of 3700 km from Jacksonville (Florida) to Portland (Maine) in different environments (urban, suburban, commercial, and highway). In this contribution we focus on data from Portland (Maine) comprising urban and suburban environment. Figure 4 illustrates the measurement route in urban and suburban environments (depicted by the different colors) in Portland, Maine.
3.3. Evaluation of the ImageBased State Estimation Method
In order to validate the imagebased method we compare the estimated channel states of the images with actual radio frequency (RF) Data.
In LMS channel models the sequence of states is commonly modeled by Markov/semiMarkov processes [1, 37]. Therefore, parameters like state probabilities and the distribution of state durations are of interest. The state duration statistics are derived for the satellites Sirius 1, Sirius 2, XM3, and XM4 with an average elevation of 65°, 68°, 38°, and 23°, respectively. The results based on image and RF data are depicted in Figure 5 for bad state and in Figure 6 for good state. In order to evaluate the estimation accuracy, we show the Complementary Cumulative Distribution Function (CCDF).
(a) Sirius 1
(b) Sirius 2
(c) XM3
(d) XM4
(a) Sirius 1
(b) Sirius 2
(c) XM3
(d) XM4
The statistics show very good agreement between the RFbased (dashed lines) and the imagebased results (solid lines). Additionally, the lognormal fit is plotted in black. As a result, the imagebased approach can be seen as a valid technique for state estimation.
4. ImageBased Model of the Land Mobile Satellite Channel
In the following we explain the proposed LMS channel model. Statistical channel models reproduce the timeseries of the received satellite signal comprising two processes, the very slow fading and the slow/fast fading. The nonstationarity of the signal due to the different degrees of shadowing is expressed by states (e.g., LoS, moderate shadowing, and deep shadowing). Well accepted in the literature are LMS models, which describe the very slow fading component by either three states [1], namely, lineofsight, shadowed, and blocked, or by two states [2], good and bad. The slow and fast fading caused by varying shadowing and multipath effects are treated as a second process and is commonly reproduced by a stationary stochastic process.
In Figure 7, the LMS model is illustrated, which is based on the versatile twostate approach presented in [2, 5]. The model distinguishes between two states, good and bad, corresponding to LoS or light shadowing and NLoS or heavy shadowing, respectively. The state modeling is performed by a semiMarkov chain [37], which improves the accuracy of the state duration modeling, compared to the more common first order Markov chains [1]. Based on the measured state durations, the state duration probability density function (SDPDF) is modeled by a lognormal distribution. It is expressed by where is the state duration, is the mean value of , and is the standard deviation of . These parameters and are called the lognormal mean and lognormal standard deviation, respectively, and are estimated for both good and bad state. The mean of the state duration (in meter) can be calculated by ; the standard deviation (in meter) can be calculated by .
The slow and fast signal fading within each state are modeled by a Loo distribution [27]. The Loo distribution considers the received signal as a sum of two signal components. A lognormally distributed direct signal expresses the slow fading component corresponding to varying shadowing conditions of the direct signal. A Rice distribution characterizes the fast fading component due to multipath effects. The Loo parameter triplet consists of the mean , the standard deviation for the lognormally distributed direct signal, and the average multipath power . Depending on the current state interval and on the environment of the terminal, a new random Loo parameter triplet is generated. Compared to an earlier threestate model [1], where the Loo parameter triplet is fixed for each state, the versatile Loo parameter selection enables a more realistic modeling over the full dynamic range of the received signal [2]. The output of the channel model is a timeseries of the received signal in form of a complex envelope. For the purpose of comparability, we convert the output to a timeseries of the carriertonoise ratio (). As the LoS level is dependent on the transmit power, pathloss (distance, weather, and frequency), and the receiver characteristics, we remove it from the signal to obtain a representation of the fading effects exclusively caused by the environment. Thus, we end up with a carriertonoise ratio normalized to the LoS level in dB.
Figure 8 shows how the parameters of the model are preprocessed and how they are applied. First, hemispheric images are binarized into the classes sky and obstruction. Second, a set of parameters for the state duration statistics is extracted from the images, depending on the environment type, satellite elevation, and satellite azimuth. Depending on the user’s needs, representative parameters are selected for the actual LMS channel modeling process. Despite the extracted state duration parameters (very slow fading), the Looparameters (slow and fast fading) are chosen. Finally, a timeseries of the LMS channel in terms of the complex envelope is simulated.
In Figure 9 an overview of the detailed implementation of the LMS channel model is depicted. The model is structured in three different generator types [5].(1)The State Sequence Generator (SSG) reproduces as sequence of good and bad states by a semiMarkov model process. The necessary parameters are the lognormal mean and lognormal standard deviation of the state duration distribution for the good and bad state, respectively (cf. (1)). With this approach, an accurate modeling of the state duration and the state probabilities can be provided. At the th interval, the SSG outputs the current state, which can be either good or bad, and the duration for which the channel remains in the same state. The lognormal state duration parameters (, ), as input parameters of the semiMarkov model, are provided in Tables 2 and 3.(2)The Propagation Parameter Generator (PPG) is invoked at each new state interval producing a new random Loo triplet . As depicted in Figure 9, the PPG features three independent standard normal generators, whose means and variances () are estimated from measured RF data by curvefitting as presented in [38].(3)The SmallScale Fading Generator (SSFG) computes the Loo distributed timeseries including Doppler shaping for every new state interval, which is the output of the proposed LMS channel model. Furthermore, a Doppler shift of the lognormal varying direct signal is included, which is related to the speed of the vehicle. An optional Doppler shift of all signal components is also possible, to include the movement of the satellite.


5. Results
In the following we show the results obtained by the LMS channel model explained in Section 4, where the model parameters have been obtained with the imagebased state estimation presented in Section 3. The first part of this section focuses on the state parameters of the model, which are important for a reliable state sequence generation. Secondly, the output of the whole channel model is evaluated in terms of first order statistics.
Exploiting the imagebased state parameter estimation (instead of RF measurements) offers the advantage to perform a reliable analysis of the impact of the driving direction. Without any further geometrical consideration, a statistical availability prediction depending on the direction of the road is possible.
One crucial parameter of satellite service prediction is the service unavailability, or socalled bad state probability . Figure 10 depicts with respect to different elevation angles, azimuth angles (one degree resolution), and environment types (urban in Figure 10(a) and suburban in Figure 10(b)). The different curves represent different elevation angles, from which an inverse proportionality of against the elevation angle can be seen. Moreover, a nonsymmetric dependency on the azimuth angle relative to the heading (abscissa) is visible. In both plots two maxima can be identified, one global at 90° azimuth and one local at 270° azimuth, which correspond to the right side and left side of the vehicle, respectively. The global maximum is caused by the righthand driving on the street; therefore obstacles at the right side of the vehicle appear to be higher in relation to the position of the terminal on the street. This result shows that an assumption of the symmetric dependency of the azimuth angle is an oversimplification. At 0° and 180°, corresponding to front and back of the vehicle, two minima are found. This clearly reflects the street canyon characteristics, especially in urban scenarios. The influence of the street canyon decreases with an increasing elevation angle. The constant (solid red line) results from a calculation of , which ignores the heading of the vehicle (exemplary for an elevation angle of 30°). Here, the heading is assumed to be uniformly distributed, representing the stateoftheart LMS channel models, which do not include the azimuth angle of the satellite, for example, [2]. Ignoring the driving direction leads to a significant overestimation or underestimation of the bad state probability . The difference between urban (Figure 10(a)) and suburban environment (Figure 10(b)) becomes visible as approaching higher elevations. In contrast to urban environment, the bad state probability for suburban does not become zero due to the trees close to the road, reaching into regions of high elevation.
(a) Urban environment
(b) Suburban environment
The hemispheric images offer the possibility to analyze with a resolution of one degree in azimuth and in elevation. Based on the results of Figure 11 we further focus on eight exemplary satellite positions relative to the street orientation, see Figure 11. Turning clockwise, the positions are front, frontright, right, backright, back, backleft, left, and frontleft and correspond to the azimuthal angles 0°, 45°, 90°, 135°, 180°, 225°, 270°, and 315° relative to the driving direction, respectively.
In order to model the state sequence by a Markov process, the state probability , the lognormal mean of the state duration , and the lognormal standard deviation are of prime interest. In Figures 12 and 13 the three parameters are plotted versus the elevation angle (left for urban and on the right for suburban environment). For the sake of clarity we only show four different main positions of the satellite relative to the street orientation (front, right, back, and left). In Tables 2 and 3, the lognormal state duration parameters (, ) for the urban and suburban environment are provided for eight satellite positions in azimuth and elevation angles ranging from 15° to 85°. Due to the limited building height, the values of the bad state probability, the lognormal mean for bad, and the lognormal standard deviation for bad state decrease with an increasing elevation angle. The lognormal mean good state duration and the good state lognormal standard deviation show an inverse relationship. In particular, for urban scenarios a strong difference between the values of all the parameters for different satellite positions can be noticed; for example, a difference up to 20% in service unavailability can be seen. The estimated lognormal mean and the lognormal standard deviation are crucial parameters to describe the state duration probability density function (SDPDF) of the utilized semiMarkov chain (see Section 4). As we have shown that both parameters are dependent on the driving direction, we present the resulting complementary SDPDFs for the good and bad state in Figures 14 and 15, respectively. Each figure comprises urban (left) and suburban environment (right). From top to bottom, the elevation angle of the satellite is increasing from 20° to 40° and 60°, and again the four different positions in azimuth are analyzed. In urban environment considerable differences in the state duration statistics can be identified, where right represents the worst case and front/back the best case. For suburban scenarios the differences in the statistical behavior decrease. In the first part of this section we reveal significant differences of the state parameters depending on the environment, elevation, and on the azimuth relative to the driving direction.
(a) Urban environment
(b) Subrban environment
(a) Urban environment
(b) Suburban environment
(c) Urban environment
(d) Suburban environment
(e) Urban environment
(f) Suburban environment
(g) Urban environment
(h) Suburban environment
(a) Satellite at 20° elevation angle in urban environment
(b) Satellite at 20° elevation angle in suburban environment
(c) Satellite at 40° elevation angle in urban environment
(d) Satellite at 40° elevation angle in suburban environment
(e) Satellite at 60° elevation angle in urban environment
(f) Satellite at 60° elevation angle in suburban environment
(a) Satellite at 20° elevation angle in urban environment
(b) Satellite at 20° elevation angle in suburban environment
(c) Satellite at 40° elevation angle in urban environment
(d) Satellite at 40° elevation angle in suburban environment
(e) Satellite at 60° elevation angle in urban environment
(f) Satellite at 60° elevation angle in suburban environment
Building on the results presented so far, we now investigate the impact of the driving direction by comparing the performance of the imagebased LMS model output (see Section 4). The evaluation is performed in terms of first order statistics (Cumulative Distribution Functions (CDFs)) of the LoSnormalized carriertonoise ratio ( in dB), for eight exemplary satellite positions relative to the street orientation, see Figure 11. Turning clockwise, the positions are front, frontright, right, backright, back, backleft, left, and frontleft and correspond to the azimuthal angles 0°, 45°, 90°, 135°, 180°, 225°, 270°, and 315° relative to the driving direction, respectively. By extracting the satellite reception states from the images for the particular satellite position (see Section 3.1), the driving direction dependency is considered within the state parameter set. The Loo parameters (, , ), chosen from statistical distributions, depend only on the environment and the elevation angle [5]. The results, corresponding to routes in the urban and suburban environment of Portland (Maine), are depicted in Figures 16 and 17, respectively. Each plot refers to a specific elevation angle (from 20° to 70°) and the colored curves with different markers represent the eight chosen satellite positions with respect to the driving direction. As a reference, the black curve is generated by considering a route in the corresponding environment, while the heading of the vehicle is assumed to be uniformly distributed. This method represents traditional modeling without an explicit consideration of the azimuth angle.
(a) Satellite at 20° elevation angle
(b) Satellite at 30° elevation angle
(c) Satellite at 40° elevation angle
(d) Satellite at 50° elevation angle
(e) Satellite at 60° elevation angle
(f) Satellite at 70° elevation angle
(a) Satellite at 20° elevation angle
(b) Satellite at 30° elevation angle
(c) Satellite at 40° elevation angle
(d) Satellite at 50° elevation angle
(e) Satellite at 60° elevation angle
(f) Satellite at 70° elevation angle
For all plots a clear tendency of an increasing LoSnormalized related to an increasing elevation angle can be seen. The limited height of obstacles in the vicinity of the vehicle leads to better LoS conditions for higher elevations. At low elevation angles of the satellite (20°–30°) the impact of the street canyon becomes obvious (Figures 16(a) and 16(b)). A satellite oriented in the driving direction (front) or in the opposite (back) leads to an almost 17 dB higher LoSnormalized compared to satellites positioned behind the buildings on the left or right side. The reference curve (ignoring the heading) is very close to the worst case.
In Figures 16(c) and 16(d), a third group with approximately 6 dB better than the worst case can be identified, consisting of the satellite positions frontleft, backleft, and left, respectively. Here, the effect of the righthand driving of the vehicle becomes evident; objects at the right side are closer and in relation to the position of the car, higher in elevation. The reference estimates the channel close to the worst case or close to the CDFs of the left positions. Consequently, modeling the satellite signal without any information of the heading of the vehicle leads to an overly optimistic estimation for satellites positioned at the right side.
Figure 16(e) shows the diversity of the satellite channel at an elevation of 60°. The CDFs of frontleft and backleft superpose with front and back. The satellite at the left position behaves slightly worse. The LoSnormalized decreases even more for frontleft and backleft (7 dB) and right is still the worst case. The reference CDF is approximately located in the middle of the set of curves.
For 70° elevation (Figure 16(f)) the positions on the right side represent the worst case with a difference of 15 dB to the other superposed curves. Again, the reference curve represents almost an average of worst case and best case.
In Figure 17 the results for suburban environment are shown. Not surprisingly, the set of CDFs is closer to each other (10 dB); nonetheless the same worst and best cases can be identified.
6. Conclusions
In this contribution, we present a statistical Land Mobile Satellite (LMS) channel model, which includes the azimuth angle of the satellite with respect to the driving direction. To derive the necessary LMS model parameters for different headings and arbitrary satellite positions, an imagebased parameter estimation method is used. This allows us to reduce the measurement effort significantly, compared to a conventional RF measurement campaign. We verify the imagebased method by comparing the statistical state parameters derived from images with the ones derived from the corresponding measured RF signal levels, serving as a ground truth. The results show a good match, demonstrating that the imagebased approach is a reliable method for detecting the very slow fading of the received satellite signal.
The proposed LMS channel model comprises a semiMarkov model using imagebased state duration parameters to generate realistic state sequences for different driving directions. In particular, the extracted state duration statistics are strongly influenced by the driving direction. The model is able to generate timeseries of the received satellite signal for different azimuth positions relative to the heading. We show that varying headings lead to significant differences in the first order statistics of the carriertonoise ratio, demonstrating the considerable nonstationarity of the received satellite signal. Although the proposed approach provides an accurate and realistic model, only a few user input parameters are necessary. This makes it an attractive tool for designing and testing of satellite systems.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
The measurements were carried out in the context of the project MiLADY (Mobile satellite channeL with Angle DiversitY), funded by the European Space Agency (ESA) under contract no. C21150.
References
 F. PérezFontán, M. VázquezCastro, C. E. Cabado, J. P. García, and E. Kubista, “Statistical modeling of the LMS channel,” IEEE Transactions on Vehicular Technology, vol. 50, no. 6, pp. 1549–1567, 2001. View at: Publisher Site  Google Scholar
 R. PrietoCerdeira, F. PérezFontán, P. Burzigotti, A. BoleaAlamaac, and I. SanchezLago, “Versatile twostate land mobile satellite channel model with first application to DVBSH analysis,” International Journal of Satellite Communications and Networking, vol. 28, no. 56, pp. 291–315, 2010. View at: Publisher Site  Google Scholar
 E. Lutz, “A Markov model for correlated land mobile satellite channels,” International Journal of Satellite Communications, vol. 14, no. 4, pp. 333–339, 1996. View at: Publisher Site  Google Scholar
 D. Arndt, T. Heyn, A. Heuberger, R. PrietoCerdeira, and E. Eberlein, “State modeling of the land mobile satellite channel with angle diversity,” in Proceeding of the 6th European Conference on Antennas and Propagation (EuCAP '12), pp. 3140–3144, Prague, Czech Republic, March 2012. View at: Publisher Site  Google Scholar
 D. Arndt, T. Heyn, J. König et al., “Extended twostate narrowband LMS propagation model for SBand,” in Proceedings of the IEEE International Symposium on Broadband Multimedia Systems and Broadcasting (BMSB '12), pp. 1–6, Seoul, Republic of Korea, June 2012. View at: Publisher Site  Google Scholar
 D. Arndt, A. Ihlow, T. Heyn, A. Heuberger, R. PrietoCerdeira, and E. Eberlein, “State modelling of the land mobile propagation channel for dualsatellite systems,” EURASIP Journal on Wireless Communications and Networking, vol. 2012, article 228, 2012. View at: Publisher Site  Google Scholar
 D. Arndt, A. Ihlow, A. Heuberger, and E. Eberlein, “State modelling of the land mobile propagation channel with multiple satellites,” International Journal of Antennas and Propagation, vol. 2012, Article ID 625374, 15 pages, 2012. View at: Publisher Site  Google Scholar
 C. Tzaras, S. R. Saunders, and B. G. Evans, “A physicalstatistical pro pagation model for diversity in mobile satellite PCN,” in Proceedings of the 48th IEEE Vehicular Technology Conference (VTC '98), vol. 1, pp. 525–529, Ottawa, Canada, May 1998. View at: Google Scholar
 C. Oestges, S. R. Saunders, and D. VanhoenackerJanvier, “Physical statistical modelling of the land mobile satellite channel based on ray tracing,” IEE Proceedings: Microwaves, Antennas and Propagation, vol. 146, no. 1, pp. 45–49, 1999. View at: Publisher Site  Google Scholar
 C. Oestges and D. VanhoenackerJanvier, “A physicalstatistical shadowing correlation model and its application to lowearthorbit systems,” IEEE Transactions on Vehicular Technology, vol. 50, no. 2, pp. 416–421, 2001. View at: Publisher Site  Google Scholar
 C. Oestges and D. VanhoenackerJanvier, “Time seriesbased prediction of fade dynamics and shado wing correlation statistics for navigation systems,” in Proceedings of the 2nd European Conference on Antennas and Propagation (EuCAP '07), pp. 1–6, IET, Edinburgh, UK, 2007. View at: Google Scholar
 M. A. VázquezCastro, F. PérezFontán, H. IglesiasSalgueiro, and M. A. BarciaFernandez, “A simple threesegment model for shadowing cross correlation in multisatellite systems in street canyons,” Microwave and Optical Technology Letters, vol. 28, no. 3, pp. 160–164, 2001. View at: Publisher Site  Google Scholar
 M. A. VázquezCastro, F. PérezFontán, and S. R. Saunders, “Shadowing correlation assessment and modeling for satellite diversity in urban environments,” International Journal of Satellite Communications, vol. 20, no. 2, pp. 151–166, 2002. View at: Publisher Site  Google Scholar
 F. PérezFontán, M. A. VázquezCastro, C. Enjamio, and J. Pita, “Modeling the building blockage crosscorrelation in multisatellite systems,” in Proceedings of the 55th Vehicular Technology Conference (VTC '02), vol. 2, pp. 928–932, May 2002. View at: Google Scholar
 F. PérezFontán, M. A. VázquezCastro, C. EnjamioCabado, and J. Pita, “Simplified multisatellite system availability calculations using street masking functions,” in Proceedings of the 55th IEEE Vehicular Technology Conference (VTC '02), vol. 4, pp. 1864–1868, May 2002. View at: Google Scholar
 F. PérezFontán, M. A. VázquezCastro, C. Enjamio, P. Mariño, and J. Pita, “Approximate urban area call dropout probability estimation in mobile multisatellite systems,” International Journal of Satellite Communications and Networking, vol. 21, no. 2, pp. 183–197, 2003. View at: Publisher Site  Google Scholar
 F. PérezFontán, S. Martínez, B. Sanmartin, C. Enjamio, P. Mariño, and F. Machado, “An enhanced Markov chain based model for the narrowband LMS channel in builtup areas,” International Journal of Satellite Communications and Networking, vol. 23, no. 2, pp. 111–128, 2005. View at: Google Scholar
 C. Meenan, M. Parks, R. Tafazolli, and B. Evans, “Availability of 1st generation satellite personal communication network service in urban environments,” in Proceedings of the IEEE 48th Vehicular Technology Conference (VTC Spring '98), vol. 2, pp. 1471–1475, IEEE, Ottawa, Canada, May 1998. View at: Google Scholar
 R. Akturan and W. J. Vogel, “Optically derived elevation angle dep endence of fading for satellite PCS,” Proceedings of the 19th NASA Propagation Experimenters Meeting (NAPEX '95), vol. 1, pp. 127–132, 1995. View at: Google Scholar
 R. Akturan and W. J. Vogel, “Path diversity for LEO satellitePCS in the urban environment,” IEEE Transactions on Antennas and Propagation, vol. 45, no. 7, pp. 1107–1116, 1997. View at: Publisher Site  Google Scholar
 R. Akturan, H.P. Lin, and W. J. Vogel, “Propagation modeling in land mobile satellite systems using photogrammetry,” in Proceedings of the IEEE 46th Vehicular Technology Conference, pp. 1018–1022, Atlanta, Ga, USA, May 1996. View at: Google Scholar
 M. Rieche, D. Arndt, A. Ihlow, F. PérezFontán, and G. Del Daldo, “Impact of driving direction on land mobile satellite channel modeling,” in Proceedings of the 8th European Conference on Antennas and Propagation (EuCAP '14), The Hague, The Netherlands, 2014. View at: Google Scholar
 R. Akturan and W. J. Vogel, “Photogrammetric mobile satellite service prediction,” Electronics Letters, vol. 31, no. 3, pp. 165–166, 1995. View at: Publisher Site  Google Scholar
 L. Castanet, “Influence of the variability of the propagation channel on mobile, fixed multimedia and optical satellite communications,” in Mobile/Indoor Multipath, L. Castanet, Ed., Shaker, Aachen, Germany, 2008. View at: Google Scholar
 A. G. Kanatas and P. Constantinou, “A narrowband land mobile satellite channel software simulator for urban environments,” International Journal of Satellite Communications, vol. 18, no. 1, pp. 17–45, 2000. View at: Google Scholar
 M. Döttling, A. Jahn, D. Didascalou, and W. Wiesbeck, “Twoand threedimensional ray tracing applied to the land mobile satellite (LMS) propagation channel,” IEEE Antennas and Propagation Magazine, vol. 43, no. 6, pp. 27–37, 2001. View at: Publisher Site  Google Scholar
 C. Loo, “A statistical model for a land mobile satellite link,” IEEE Transactions on Vehicular Technology, vol. 34, no. 3, pp. 122–127, 1985. View at: Google Scholar
 G. E. Corazza and F. Vatalaro, “A statistical model for land mobile satellite channels and its application to nongeostationary orbit systems,” IEEE Transactions on Vehicular Technology, vol. 43, no. 3, pp. 738–742, 1994. View at: Publisher Site  Google Scholar
 H. Suzuki, “A statistical model for urban radio propogation,” IEEE Transactions on Communications, vol. 25, no. 7, pp. 673–680, 1977. View at: Publisher Site  Google Scholar
 E. Lutz, D. Cygan, M. Dippold, F. Dolainsky, and W. Papke, “The land mobile satellite communication channel: recording, statistics, and channel model,” IEEE Transactions on Vehicular Technology, vol. 40, no. 2, pp. 375–386, 1991. View at: Publisher Site  Google Scholar
 Y. Karasawa, K. Minamisono, and T. Matsudo, “A propaga tion channel model for personal mobilesatellite services,” in Proceedings of Progress of Electromagnetic Research Symposium of the European Space Agency (ESA '95), pp. 11–15, Noordwijk , The Netherlands, July 1995. View at: Google Scholar
 S. Scalise, C. Alasseur, L. Husson, and H. Ernst, “Accurate and novel modeling of the land mobile satellite channel using reversible jump markov chain monte carlo technique,” in Proceedings of the IEEE Global Telecommunications Conference (GLOBECOM '06), pp. 1–6, San Francisco, Calif, USA, November 2006. View at: Google Scholar
 M. Milojević, M. Haardt, A. Heuberger, and E. Eberlein, “Channel state modeling for single and multiple satellite broadcasting systems,” in Proceedings of the International Workshop on Satellite and Space Communications (IWSSC '08), pp. 102–106, Toulouse, France, October 2008. View at: Publisher Site  Google Scholar
 H.P. Lin, R. Akturan, and W. J. Vogel, “SatellitePCS channel simulation in mobile user environments using photogrammetry and Markov chains,” Wireless Networks, vol. 3, no. 4, pp. 299–308, 1997. View at: Publisher Site  Google Scholar
 M. Rieche, D. Arndt, A. Ihlow, and G. Del Galdo, “State modeling of the land mobile satellite channel by an imagebased approach,” in Proceedings of the 7th European Conference on Antennas and Propagation (EuCAP '13), pp. 672–676, Gothenburg, Sweden, April 2013. View at: Google Scholar
 E. Eberlein, A. Heuberger, and T. Heyn, “Channel models for systems with angle diversity—the MiLADY project,” in Proceedings of the ESA Workshop on Radiowave Propagation Models, Tools and Data for Space Systems, Noordwijk, The Netherlands, December 2008. View at: Google Scholar
 L. E. Bråten and T. Tjelta, “SemiMarkov multistate modeling of the land mobile propagation channel for geostationary satellites,” IEEE Transactions on Antennas and Propagation, vol. 50, no. 12, pp. 1795–1802, 2002. View at: Publisher Site  Google Scholar
 J. König, D. Arndt, A. Ihlow, and A. Heuberger, “Estimation of fading parameter correlation for modeling the land mobile satellite channel,” in Proceedings of the IEEE International Symposium on Broadband Multimedia Systems and Broadcasting (BMSB '12), Seoul, Republic of Korea, June 2012. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Marie Rieche et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.