Fading Analysis for the High Speed Railway Viaduct and Terrain Cutting Scenarios

A good understanding of the fading characteristics in high speed railway environment is essential for the design of a high reliable railway wireless network. In this paper, measurements have been taken in both high speed viaduct and terrain cutting scenarios using track side base stations (BSs) of the railway wireless network in China. The measurement sites have been chosen with special care; thus the whole measured route can be characterized either by viaduct or terrain cutting. Kolmogorov-Smirnov (K-S) test has been first introduced in the statistical analysis to find out which is the most appropriate model for the small scale fading envelope. Though both Rice and Nakagami distributions provide a good fit to the first-order envelope data in both scenarios, only the Rice model generally fits the second-order statistics data accurately. For the viaduct scenario, higher Rice K factor can be observed. The change tendency of the K factor as a function of distance in the two scenarios is completely different. It can be concluded that the over bridges which span the railway tracks in terrain cutting scenario would affect the Rice K value severely.


Introduction
As an international standard for railway communication and applications, GSM-R (GSM for Railway) has been adopted by China as a wireless communication platform to transfer security data of train control.With the rapid development of high speed railways in China, GSM-R networks have been widely deployed along the railway lines.Since GSM-R aims to transmit security data and information for train control, failure or obstruction of the wireless network will inevitably affect the normal running of the railway system.To provide a safe and reliable network, the design of a railway wireless network is of vital importance.In the meanwhile, appropriate wireless network design and optimization rely heavily on accurate prediction of radio wave propagation.
The propagation environment of the high speed railways is a lot different from the common public wireless network.As indicated in [1], typically a railway track is full of terrain cuttings, tunnels, bridges, and so forth.Track side base stations (BSs) are often used to provide a seamless coverage.Besides, the BS antenna is always much higher than the surrounding environment to insure that there is a lineof-sight (LOS) component between the BS antenna and the train antenna.So far, channel analysis in these special railway environments has attracted more and more research interests.The measured path loss values in viaduct scenarios have been analyzed in [2][3][4], respectively, and the path loss exponents have been found between 2 and 4. A more general path loss model for viaduct scenarios has been proposed in [5], considering the influences of viaduct height and BS antenna relative height.The Rice K factor in the viaduct and terrain cutting scenarios along the high speed railway was first estimated by a moment-based estimator, which employed the second and the fourth moments of the signal envelope [6].However, statistical analysis of two momentbased estimators shows that moment-based estimator using the first and second moments of the signal envelope can be a better estimator, as it has better asymptotic performance [7].And Rice K factor as a function of distance has been studied in [8], which shows the K factor decreases linearly in the distance.However, the estimation interval of the Rice K factors was set to 100 m, which we believe is too large to reflect the fast variation of the K factors.
Besides, though Nakagami distribution has been widely accepted to be a very good fit to mobile radio channel characteristics [9,10], it has never been studied for the small scale envelope characteristics in railway scenarios.The objective of this work is to fill this gap by investigating the suitability of Nakagami distribution for describing the small scale fading envelope.Moreover, second-order statistics of the received envelope, such as the level crossing rate (LCR) and average fade duration (AFD), have not been fully investigated.Although recent work compared the measured second-order statistics with the theoretical values of Rayleigh, Rice, and Nakagami models in viaduct scenario [11], no paper has reported the results of terrain cutting scenario.These statistics provide useful information for the statistic analysis of burst errors and the design of errorcorrecting codes.
This rest of the paper is organized as follows.A brief description of the measurement system and scenarios is given in Section 2. Theoretical models for the small scale fading channel are described in Section 3. Section 4 shows the measurement results of the two scenarios, including the envelope distribution and second-order statistics.Conclusions are drawn in Section 5.

Measurement Setup
The narrowband measurements were done along the Zhengzhou-Xi'an high-speed railway of China with a special test system, which can collect signal level data from GSM-R track side BSs.During the measurements, the train was moving at a high speed of about 277-300 km per hour.The measurement system in the train periodically recorded the signal power value every 10 centimeters.Two typical viaduct and terrain cutting scenarios have been selected for analysis.And the measurements have been done twice in two different days to eliminate the effect of the measurement error.The following sections describe the measurement system and the measurement scenarios in greater detail.
2.1.Measurement System.The measurement system shown in Figure 1 consists of a receive antenna mounted on the top of the train, a Willteck 8300 Griffin fast measurement receiver, a train odometer, a global positioning system (GPS) device, and a laptop to record the test data.The measurement receiver operated in the distance trigger mode.So the train odometer was able to send a trigger signal to the receiver at intervals of 10 centimeters.A software installed in the laptop can store the measured data.And geographic location data of the measured samples can be provided by the GPS device.Table 1 shows some measurement parameters of the two A satellite image of the measured terrain cutting is also displayed in Figure 2(d).Five over bridges which span the railway tracks can be observed in this figure.These bridges are built to facilitate traffic on the two sides of the terrain cutting.Later the effect of the over bridges on radio propagation will be illustrated.BS ZX1714D was transmitting signal during the test.The whole measurement area from BS ZX1714D to BS 1713D can be characterized by terrain cutting scenario.In both scenarios, the railway tracks are almost straight and the BS antennas are situated 17 and 37 m away from the railway tracks.

Theoretical Models
Various theoretical models have been presented to describe the small scale fading behavior of the mobile channel.In this section, both the envelope distribution and second-order statistics of the Rayleigh, Rice, and Nakagami fading models are discussed.

Envelope Distribution.
The Rayleigh probability density function (PDF) is given by [12] p(r  And the maximum likelihood estimate of the parameter σ is where E(r 2 ) is the second moments of the measured samples.The Rice PDF is expressed as where function I 0 is the modified Bessel function of the first kind and zero-order.The Rice distribution is often described in terms of a parameter K which is defined as the ratio between the direct or strong component power of the signal and the variance of the multipath.It is given by The parameter K is known as the Rice factor and completely specifies the Rice distribution.Different methods have been proposed for estimating the Rice K factor.Though the maximum likelihood estimator (MLE) can yield optimal results, it is relatively cumbersome and time consuming [13].A simple K estimation method based on the moments of the receiver envelope has also been proposed.It has been validated by [7] that moment-based estimator using the first and second moments of the signal envelope has a better asymptotic performance, compared to the one using the second and fourth moments of the envelope.The estimator based on the first and second moments of the signal envelope is used in this paper, though it is a little complex.
The moments of the Rice distribution can be expressed as [7] where 1 F 1 (•; •; •) is the confluent hypergeometric function, and Γ(•) is the gamma function.To derive K, define the following functions of K Since f n,m (K) depends only on K, K can be derived by inverting the corresponding f n,m (K).When n = 1 and m = 2, (6) can be calculated using (5) as From ( 6), f 1,2 (K) can also be expressed as where μ 1 and μ 2 can be calculated according to the moments of the measured samples, which can be given by where R(lT s ) is the measured envelope value, T s is the sampling period, and N is the number of available samples.The corresponding estimator K 1,2 can be obtained by combining International Journal of Antennas and Propagation (7) and (8).However, this involves the complex numerical procedure of inverting (7).A simple way to solve the problem is to implement a lookup table.To be specific, first create a lookup table by calculating all the f 1,2 values for a predefined range of K values according to (7).Then compare the f 1,2 results from (8) with the values stored in the lookup table.At last the K value which gives the least deviation of f 1,2 is chosen as the estimated K factor.The Nakagami PDF of the envelope r is given by where m is the Nakagami parameter defined as the ratio between Ω squared and the variance of the envelope squared and Ω is the average power, given by The Gamma function is given by As is generally well accepted, the parameter Ω can be reasonably estimated by where N is the number of available samples R i of the envelope.On the contrary, the shape parameter m can be estimated using different methods.In [14], the performance of the inverse normalized variance, Tolparev-Polyakov, and the Lorenz estimators have been compared through Monte Carlo simulation.And the inverse normalized variance estimator has been proved to be superior to the other two estimators.Thus, it is adopted as the estimator in our paper, which is given by 3.2.Second Order Statistics.The LCR of the signal envelope reveals information about how fast the received signal changes with time.It is defined as the average number of times the signal envelope crosses a certain threshold level in a positive-going direction per second.And the AFD is defined as the average time that the fading signal envelope remains below a certain threshold level [15].Exact formulas for second-order statistics of the Rayleigh, Rice, and Nakagami fading models have been derived in the literature and their fit to the empirical data has also been investigated in various scenarios.It is widely accepted that the goodness of fit for second-order statistics do not appear to be dependent on the accuracy of fit for first-order statistics [16,17].As a result, it is essential to compare our measurement results to all of the three theoretical models, which are listed below.For a fading signal, the LCR expressions for Rayleigh, Rice [18], and Nakagami [19] models are where is the value of the specified level R, normalized to the local rms amplitude of the fading envelope, f d is the maximum Doppler frequency, which can be calculated by where υ is the average speed of the mobile, and λ is the wavelength of the carrier signal.And K is the Rice factor, m is the Nakagami parameter.Their values can be estimated from the measured envelope samples.
The AFD formulas for Rayleigh, Rice, and Nakagami models can be expressed as

Data Processing.
To determine the small scale fading statistics of the received signal envelope, the effects of the path loss and shadowing have to be removed first.In order to extract the fading envelope, the received signal is normalized to its local mean value.So, for the received sample r(x i ), the local mean value is given by where W represents the number of samples in a sliding window (or bin) for the computation.This local mean value is computed for each individual sample in a bin and the normalized samples are then used for distribution fitting.As suggested by [20], a bin size of 40 λ is used for analysis, which corresponds to an average distance of approximately 13 m.For the sampling rate used for these measurements, this distance translates to approximately 128 samples per bin.For the first-order envelope, Kolmogorov-Smirnov (K-S) test is implemented to investigate the fitting performance of the predefined distributions for all the bins in both scenarios.It has been applied in [21][22][23] as a goodness-of-fit test method to verify the suitability of a hypothesized distribution and identify the distribution which best represents the experimental results.
The K-S statistic D N , which measures the maximum deviation between the hypothesized distribution and the empirical distribution derived from the measured data, is given by [24] where S N (x) is the cumulative distribution function (CDF) of the empirical data, F(x) is the CDF of the hypothesized distribution, and N is the number of samples.The value D N is then compared with a critical value J, which is a function of significance level and the sample size N [25].If D N < J, then the hypothesized distribution is equal to the empirical distribution.In other words, the hypothesized distribution has passed the K-S test.

Distribution Fitting Results
. Table 2 shows the K-S testing results.The numerical values represent the percentages of the bins which have passed the test.It is clear that both Rice and Nakagami distribution fit the measured data very well in the two scenarios, as they both have fitted more than 74% of the bins.The differences between the fitting results of the Rice and Nakagami distribution are fairly small.On the contrary, Rayleigh distribution has only fitted a minority of tested bins.Additionally, sample plots of the empirical CDF for a particular bin and theoretical model fits are given in Figures 3 and 4  to the randomness and complexity of the multipath signal components, signal envelope samples of some bins may show a complete different distribution.Thus, in order to obtain more accurate results, the estimated K values which haven't passed the K-S test are discarded in the statistical analysis.Table 3 presents some statistics of the estimated Rice factors (as a linear power ratio) in both scenarios, including the minimum, maximum, mean, and standard deviation values.It is seen that viaduct scenario appears to be a less severe fading channel compared to terrain cutting scenario since the mean K value of the viaduct scenario is larger than that of the terrain cutting scenario.The standard deviation values of both scenarios are approximately the same.
The CDF of the obtained K factors in both scenarios is also shown in Figure 5.It can be seen that for any given probability value, the K value of the viaduct scenario is larger than that of the terrain cutting scenario, further confirming the relative severe propagation environment of the terrain cutting scenario.In most of the cases (about 90%), K value is within the interval 0-4 for viaduct scenario, while the value is 0-2.9 for terrain cutting scenario.In addition, approximately 22% of the K values in terrain cutting scenario are zero, which implies a kind of fading as severe as Rayleigh fading.However, the proportion of zero K values in viaduct scenario is only 9%.
Figure 6 shows the plot of the Rice K factors as a function of BS-train separation distance for both scenarios.In viaduct scenario, the Rice K value tends to increase in the first few hundred meters.This could be explained by the fact that directional antennas are used in the railway wireless networks.The effect of the directional antenna is also reflected in the received signal power, which is low near the BS and gradually increases with the distance.Then the K value shows a general decreasing trend.In terrain cutting scenario, the increasing trend of K value for the first few hundred meters is not obvious.The reason for this is that existence of slopes may weaken the LOS component.The rest part of the Rice K values does not have a regular change tendency as in the viaduct scenario.Instead, it fluctuates sharply.Locations of the five over bridges are also marked in the figure by five vertical dotted lines.It is observed that the Rice K value tends to drop rapidly after passing the over bridges and then rise again.This phenomenon is extremely obvious for the last four over bridges.Thus the K value  appears to change periodically according to the existence of the bridge, while this is not observed in viaduct scenario.

Second-Order Statistics.
All the obtained envelope data are first used to calculate the LCR and AFD values.After that, the measured envelope values are fitted to the theoretical distributions and then theoretical LCR and AFD values can be derived.For viaduct, the obtained Rice K factor is 2.32 and Nakagami parameter m is 1.99.For terrain cutting, the Rice K factor is 1.29 and Nakagami parameter m is 1.50.Both the empirical and theoretical LCR and AFD values are computed for threshold levels from −20 dB to 10 dB.
The empirical and theoretical LCR results in viaduct and terrain cutting scenario for the first measurement are shown in Figure 7. Fairly large LCR values per second can be observed in the figure.This is due to the fact that the train was moving at a very high speed.According to the theoretical formulas (16), the LCR per second increases linearly with the maximum Doppler frequency, which is a monotonic increasing function of the average speed of the mobile station.For viaduct scenario, a good agreement between the empirical LCR and Rice model can be seen for any threshold values, while a bad agreement between empirical data and the Rayleigh model can also be found.Although the Nakagami model does not give good fit for most of the threshold values, it gives excellent fit for the threshold values smaller than −15 dB.It is found that good match between the Nakagami and empirical CDFs does not guarantee a good match for the corresponding LCR curves.The same behavior is also obtained for the second measurement, confirming that close match for the first-order statistics may yield dissimilar second-order statistics [16].For the terrain cutting scenario, as shown in Figure 7  than the other two models.But the deviation between the empirical results and the Rice model is larger than that of the viaduct scenario, especially for the threshold value near 0 dB.The Rayleigh model still gives the worst fit, while the Nakagami model gives a relatively better fit.Still, a relative independence of the goodness of fit between the first-order and second-order statistics can be observed.The empirical and theoretical AFD results in viaduct and terrain cutting scenario for the first measurement are shown in Figure 8.The AFD values are plotted on a logarithmic scale to highlight the slight difference (not visible on a linear scale) between the empirical and theoretical AFDs.For both scenarios, Rice model offers the best fit for most threshold values.However, it deviates larger for the lower and upper

Conclusion
This paper has presented empirical fading characteristic results in typical high speed railway viaduct and terrain cutting scenarios.Measurements were done by a special test system for GSM-R.The small scale fading envelope distribution has been fitted to theoretical Rayleigh, Rice, and Nakagami distributions.K-S test has been specially introduced as a goodness-of-fit test method to verify the suitability of a hypothesized distribution and the results show that both Rice and Nakagami distribution can describe the empirical data very well.Statistical analysis of the Rice K factors indicates that in most cases the value is within the interval 0-4 for viaduct scenario, while the value is 0-2.9 for terrain cutting scenario.And about 22% of the Rice K values in terrain cutting scenario are zero, which reveals that radio wave experiences more severe fading in terrain cutting scenario than in viaduct scenario.Further study of the Rice K factor as a function of distance suggests that in viaduct scenario the K value tends to increase in the first few hundred meters due to the directional antenna.Then it shows a general decreasing trend for the rest part of the distance.On the contrary, the K value exhibits an irregular change tendency with distance in terrain cutting scenario.It is found that the K value appears to drop and rise periodically by the effect of the over bridges which span the railway tracks.Though both Rice and Nakagami distributions provide a good fit to the first-order envelope data in both scenarios, only the Rice model generally fits the second-order statistics data accurately.This result verifies that the goodness of fit between the first-order envelope statistics and second-order statistics is relatively independent.These results will certainly enable better comprehension of the propagation channel in typical high speed railway scenarios and thus have important implications for the system design of the network.
Satellite image of measured viaduct (c) Terrain cutting scenario ZX1714D ZX1713D (d) Satellite image of measured terrain cutting

Figure 3 :Figure 4 :
Figure 3: Sample empirical CDF of the small scale signal envelope and theoretical model fits for viaduct scenario.The Rice K factor is 5.56 and Nakagami parameter m is 4.37.

6 InternationalFigure 5 :
Figure 5: CDF of the K factors in both scenarios.

Figure 6 :
Figure 6: Rice K factor as a function of distance.
(b), the Rice model still gives better fit International Journal of Antennas and Propagation 7

Figure 7 :
Figure 7: Empirical LCR results in both scenarios for the first measurement, together with the theoretical results of Rayleigh, Rice, and Nakagami fading models.

Figure 8 :
Figure 8: Empirical AFD results in both scenarios for the first measurement, together with the theoretical results of Rayleigh, Rice, and Nakagami fading models.

Table 2 :
K-S testing results at 5% significance level.