Abstract

This paper analyses the impact of 3D multipath in an LTE-Advanced heterogeneous network. The impact of the base station array configuration is considered in both the azimuth and elevation planes. Spatial and temporal multipath statistics are generated for example macro- and picocellular base stations, and these are combined with appropriately oriented complex polarimetric antenna patterns. The resulting wideband channels are then passed to an LTE-Advanced physical layer simulator. The optimal spatial multiplexing mode is determined by computing the throughput using the received bit mutual information rate for all modulation and coding schemes. The fastest link speed for each user is identified given a packet error rate threshold of 10%. Results show that RMS angle spread statistics and the base station MIMO array configuration strongly influence user performance in an LTE-A network. The deployment of macrocellular vertical arrays is shown to significantly degrade network performance. In contrast, the elevation angle spread is found to be far greater in picocells, and this allows vertical arrays to be applied to produce compact high performance picobase stations.

1. Introduction

The growth of smart handsets has led to mobile users demanding higher data rates and better quality mobile services. Traditional wireless telecommunication systems no longer meet the aggregate traffic demands of the user base. In order to meet the needs of next generation networks, such as LTE-Advanced (LTE-A) [1], the spectrum efficiency of future mobile communication systems must be improved. This can be achieved by adopting advanced technologies such as multiple-input multiple-output (MIMO) antenna arrays and heterogeneous networks (HetNets). Developed by 3GPP, LTE-A is an evolved version of LTE. LTE-A meets the requirement of the International Telecommunication Union (ITU) for classification as a fourth generation (4G) radio standard. LTE-A makes use of carrier aggregation, advanced MIMO techniques, wireless relays, intercell interference coordination (ICIC), and coordinated multipoint (CoMP) transmission/reception. LTE-A targets low mobility (pedestrian) users at velocities up to 15 km/h.

Considerable work has been published on MIMO performance in an LTE-A network by Prasad et al. and Li et al. [2, 3]. Lee et al. in [4] discussed the MIMO technologies in 3GPP LTE and LTE-A. The work focused on average cell throughput and cell-edge user throughputs on the downlink. Results showed that for and antenna configurations the use of spatial multiplexing (SM) improved the average cell throughput. SM increases the user data rate by transmitting independent data streams simultaneously over unique spatial channels. Gains in average cell throughput of 8.2% have been reported for a MIMO configuration compared to a reference network. MIMO technologies were also shown to significantly improve cell-edge user throughput.

This paper analyses the performance of the MIMO techniques that enable LTE-A to achieve a peak spectral efficiency of 30 bps/Hz on the downlink (DL) and 15 bps/Hz on the uplink (UL). In particular, this paper quantifies the relationship between the MIMO data rates and 3D multipath angle statistics at the base station and mobile terminals. According to the mobile data traffic statistics reported by Cisco Visual Networking Index in [5], mobile data traffic is expected to grow at a compound annual growth rate (CAGR) of 92% from 2010 to 2015, reaching 6.3 exabytes per month by 2015. In order to enhance the peak data rates as well as the overall network capacity, one interesting solution is to complement the macrolayer with a number of low power picobase stations. The resulting solution is a HetNet [69]. HetNets are a new technique for increasing the capacity of LTE-A networks. The initial deployments of LTE consist of macrobase stations and user equipment (UE) terminals that cater for high mobility users. However, HetNets use a combination of macro- and picobase stations to increase the spectral efficiency per unit area. Khandekar et al. in [6] discussed the use of HetNets to enhance range.

Few papers to date have analysed the throughput performance of LTE-A users in a picocellular scenario. HetNets allow many users to associate with the pico base stations and this enables a more equitable per user distribution of the radio resources, especially in regions with hotspots. In [7] a number of picocells were deployed at the edge of the macrocell to mitigate coverage holes. Previous work by Landström et al. in [8] has studied the position of the picocells relative to the macrosites. It was found that the location of the picocells has a significant impact on capacity. The performance of three different HetNet deployments was compared by Hiltunen in [10]. The results indicated that HetNet deployments are valuable alternatives to traditional macrosites.

The impact of user handset antenna pattern and orientation was investigated by Mansor et al. [11] in a single-input single-output (SISO) HetNet scenario with a single base station and six users. This paper expands on [11] by considering a HetNet deployment for MIMO links with different base station array configurations for a very large number of links (7,283 pico base station-UE links and 4,076 macro base station-UE links). All analyses are performed in 3D space. The urban ray-tracing channels described in this paper have been used as the basis of a proposed 3D extension of the 3GPP/ITU channel model [1214]. Figure 1 illustrates a graphical representation of the simulation methodology used in this paper.

3D ray tracing is used to model all the significant multipath components between the base station and mobile terminal given the 3D geometry of the base station and mobile arrays. Measured 3D complex field antenna patterns are then spatially and polarimetrically convolved with a set of 1000 channels generated from the original ray data. The resulting channel impulse response for each link is then transformed into the frequency domain to determine the wideband channel response. As shown in the left-hand side of Figure 1, a bit accurate LTE-A physical layer simulator is used to determine the single antenna AWGN packet error rate (PER) versus signal-to-noise ratio (SNR) graph for each Modulation and coding scheme (MCS). The AWGN graphs are then passed along with the set of 1000 wideband channels into the received bit mutual information rate (RBIR) simulator to determine the instantaneous packet error rate (PER) for all MIMO MCS modes at the SNR determined from the ray-tracing data. The MCS mode that maximises the throughput for each link while maintaining a PER < 10% is then chosen by the link adaptation algorithm. Finally, statistics are generated based on the MCS modes and data throughputs of the 11,000+ links.

This work exploits the novel combination of MIMO LTE-Advanced base station geometries, detailed 3D ray-tracing data sets, measured 3D antenna patterns, and advanced physical layer simulation and link speed selection algorithms. In particular, for urban 2.6 GHz pico- and macrocellular MIMO deployments, the paper makes the following key contributions.(1)A detailed analysis of the 3D multipath statistics based on more than 11,000 ray-traced links. For each link the multipath data is spatially and polarimetrically convolved with measured 3D antenna patterns to form an integrated antenna/channel model.(2)The end-to-end performance of horizontally and vertically oriented MIMO base station antenna arrays is investigated for LTE-A. Results are compared with the azimuth, and elevation RMS angle spread results in both cell types.(3)An analysis of the eigenvalue spread of the MIMO channel and its relationship to the optimal LTE-A link speed.(4)An investigation into the likelihood (occurrence rate) of selecting a particular MIMO enabled MCS mode as a function of cell type, antenna configuration and inter-element spacing.

The outline of this paper is as follows. Section 2 presents the key parameters and assumptions behind the MIMO channel model. Section 3 describes the propagation channel model and the antenna radiation patterns used in this study. It also includes a statistical analysis of the RMS delay spread, K-factor, and 3D arrival and departure angle statistics for urban pico- and macrocellular deployments. Results focused on LTE-A performance in both cell types using vertical and horizontal base station antenna arrays (for various antenna spacings) are reported in Section 4. Metrics based on PER, throughput, and occurrence rate for specific MCS modes are presented, and these are linked to the RMS angle spread and eigenvalue spread statistics of the channel correlation matrix. Finally, conclusions are drawn in Section 6.

2. MIMO Channel Model

The performance of MIMO is critically dependent on the availability of independent spatial channels. It is well-known that channel correlation degrades the performance of a MIMO system. Spatial correlation is dependent on the 3D multipath structure. Multipath signals depart from the base station with a given Angle of Departure (AoD). At the mobile terminal, the multipath signals arrive via a unique set of Angle of Arrivals (AoA). Spatial correlation in the horizontal plane is related to the azimuth angle spread, while spatial correlation in the vertical plane is related to the elevation angle spread. Current models assume propagation in a 2D horizontal plane. However, recent work has proposed 3D extensions to the 3GPP/ITU channel model [1214]. This paper exploits these new 3D channel models to investigate the viability of vertically spaced MIMO arrays at the base station. Results for macro- and picocells are presented and compared against traditional horizontally spaced MIMO arrays.

3. System Design and Parameters

The LTE-A downlink for MIMO with two transmit and two receive antennas is studied using a baseband link-level simulator. Table 1 shows the key parameters of the LTE-A FDD downlink simulation. We assume a channel bandwidth of 10 MHz and a carrier frequency of 2.6 GHz. All physical layer parameters, unless explicitly stated, are listed in Table 1. A receive noise figure (NF) of 5 dB was assumed at each base station [9]. The PER for each of the MIMO OFDMA PHY layer transmission modes is simulated as a function of SNR using MIMO channel data extracted from our 3D ray tracer (see Section 4). The predicted channel data is combined with appropriately oriented complex polarimetric antenna patterns for the BS and each UE location.

Three data modulation schemes are supported in the LTE-A system (QPSK, 16-QAM, and 64-QAM) along with four different coding rates. The MCS modes considered in our simulator assume MIMO and are listed in Table 2. The achievable throughputs at the PHY layer can be calculated from the error-free data rate and the residual PER. An approximation for throughput is given by , where is the peak error-free transmission rate and PER is the residual packet error rate for a specific MCS mode.

The transmission date rate is defined as , where denotes the number of spatial streams, is number of data subcarriers, is the coding rate, is the number of coded bits per subcarrier, is the number of OFDMA symbols per time slot, and is the duration of a time slot. Based on Table 2 different MIMO modulation and coding schemes are chosen to match the SNR and spatial correlation properties of the channel.

In order to perform link-level analysis in an efficient and scalable manner, a PHY abstraction technique is required. In this paper, the RBIR abstraction technique is used to determine the optimal MCS mode. RBIR abstraction has been used by Kong et al. in [15, 16] and is fully described by Halls et al. in [17]. As described in Figure 1, the mode with the highest throughput (assuming PER < 10%) is selected for each UE location [18, 19]. The selection is optimum since it relies on simulating the throughput for all MCS modes for all UE locations.

Figure 2(a) shows the PER versus mean SNR from the full bit level LTE-A simulator in an ideal AWGN channel. This result was verified using the data reported in [20]. This data is used for link-to-system mapping in the RBIR abstraction engine. As shown in Figure 2(b), the abstraction approach has been thoroughly validated against our LTE-A link-level simulator in order to validate its accuracy.

In this investigation, the impact of MIMO antenna array configuration in an interference-free urban scenario is considered. As illustrated in Figure 3, both horizontal and vertical array configurations for various antenna spacings at the base station and handset are considered. The macrocell and picocell BSs were modelled to employ two horizontally and two vertically spaced antennas. Meanwhile, the UE was modelled with a maximum of two horizontally spaced antennas (to exploit the higher angle spread in the azimuth plane). For each BS-UE link, a set of 1000 uncorrelated channel snap-shots is generated. Based on this specific channel data and antenna structure, the resulting PER and throughput are calculated. A noise floor of −99.4 dBm is assumed at the mobile receiver [19]. The terms handset or user equipment (UE), base station (BS) and eNode B (eNB) or macro eNode B (macro-eNB), and picobase station or pico eNode B (pico-eNB) are used interchangeably in this paper.

4. Channel Modelling and Antennas

The radiation patterns of a macrocell BS antenna (6-element array with patches in slant polarisation, downtilted by 10°), a picocell BS antenna (patch in slant polarisation), and a UE antenna (a common mobile phone) were measured in our anechoic chamber. All patterns are recorded in 3D and include full power, phase, and polarisation information. The 3D patterns provide insight into the behavior of the antenna systems that is simply not possible from single planes (e.g., , , and ). Figure 4 shows the antennas and their total power radiation patterns. Table 3 summarises the radiation pattern statistics. It should be noted that all the antennas in this study were assumed to be 100% efficient; thus their gain is equal to their directivity. Both ends of the link were modelled to employ two antennas of the same type separated by different distances in either the horizontal or vertical plane, resulting in different MIMO system configurations. The propagation channel between each BS and UE was modelled via the spatial convolution of the polarimetric antenna patterns with the spatial and temporal multipath ray components from our 3D outdoor ray tracer [2124].

The ray-tracing engine identifies all possible ray paths between the base station and the mobile receiver and was used previously to generate the models reported in [13, 14, 25]. The geographic database includes terrain, buildings, and foliage. The ray-tracing model has been validated for cellular (and picocellular) applications, where the base station was located above (and well below) the rooftop level at frequencies from 450 MHz to 6 GHz, with a Root-mean-square-error (RMSE) of approximately 7 dB [2124]. The urban environment database of a 17.6 km² area in the centre of Bristol, UK, is investigated at a carrier frequency of 2.6 GHz.

For the macrocell study, 23 BSs were placed on rooftop locations (rooftop heights varied from 7 m to 122 m as shown in Figure 1). Each BS was modelled to cover a 3-sector cell, with a cell radius of 1 km. 300 UEs were randomly scattered in each sector within the BS antenna 3 dB beam width. The BS transmit power was 43 dBm, the BS antenna height was 3 m above rooftop level, and the UE antenna height was 1.5 m above ground level. A minimum distance of 50 m was assumed between the BSs and the UEs. For the picocell study, 600 BSs were randomly scattered at a height of 5 m above ground level. Each pico-BS was mounted on the external wall of a building and covered a 1-sector cell with a radius of 150 m. 50 random UEs were considered per BS at a height of 1.5 m above ground level and within the BS antenna 3 dB beam-width.

A pico-BS transmit power of 30 dBm and a minimum BS-UE separation distance of 5 m were assumed. Figures 5(a) and 5(b) show macrocell and selection of picocell examples. Figures 5(c) and 5(d) show examples of a macrocell and picocell coverage maps based on the total average received power at the user terminal. Base stations are denoted by red circles, while white crosses represent UE locations. A detailed statistical analysis of the propagation parameters can be found in [25, 26].

In order to compute a statistically valid set of wideband channel matrices suitable for Orthogonal Frequency Division Multiplexing (OFDM) modelling, the procedure reported in [27] was followed. Point-source ray tracing was conducted from each BS to each UE location. This provided information on the amplitude, phase, time delay, Angle of Departure (AoD), and Angle of Arrival (AoA) of each multipath component (MPC). The complex gain of each MPC was adjusted according to each transmitting/receiving antenna electric field pattern response for the corresponding AoD/AoA and polarisation. The double-directional time-invariant channel impulse response for the transmitting and receiving antenna links is given by [28] where where represents the Dirac delta function, is the time of flight, represents the departure/arrival angle, and is the total number of MPCs. The MPC between the transmitting and the BS antennas has a double-directional time-invariant channel impulse response , a complex amplitude (a matrix for all four polarisation combinations), a time of flight , and departure/arrival solid angles , and represents the vertical/horizontal polarisation component of the transmitting/ BS receiving antenna electric field radiation patterns. Time binning was then applied with a bin resolution of 32.55 ns. The wideband channel frequency response , where represents the frequency-domain channel for the -transmitting antenna, the -receiving antenna, and the subcarrier (the total number of subcarriers is ), was computed using a 2048-point Discrete Fourier Transform (DFT) as shown in where represents the DFT.

An example of a set of instantaneous wideband channel matrices can be seen in Figure 6 for a picocell link. To compute the average channel performance for a given link, this procedure was repeated for 1000 independent channel realizations, applying a uniformly distributed random phase to each MPC.

The previous process resulted in 7,283 pico-UE links and 4,076 macro-eNB-UE links. The process is summarised by the following 8 steps:(1)associate UEs to BSs;(2)ray trace each UE-BS link to determine relevant multipath components (in 3D space);(3)compute a set of 1000 -matrices for each UE-BS link based on isotropic antennas;(4)spatially and polarimetrically convolve the measured BS and UE antenna patterns with the 3D channel data;(5)compute a second set of 1000 -matrices for each BS-UE link;(6)perform RBIR abstraction to determine the LTE-A link performance based on the BS-UE SNR and scenario specific MIMO channel -matrices;(7)determine PER and link throughput for all MCS and MIMO modes;(8)select the MCS mode that maximises link throughput given PER < 10%.

Repeat the process for all BS-UE links.

The macrocell and picocell propagation statistics are presented in Figure 7 in terms of CDFs for K-factor, RMS delay spread, and RMS azimuth and elevation spreads (for the base station and user terminal). For the macrocells the values of the RMS azimuth spread at the base station (AoD) far exceed the values of the RMS elevation spread. It can also be seen that for picocellular base stations there is a considerable increase (relative to the macrocells) in both the RMS azimuth and RMS elevation spreads. Figure 7 also shows that the RMS delay spreads are statistically higher in the macrocells (since they illuminate a larger volume of the city), as are the values of K-factor (higher mounted macrocells are more likely to generate a dominant multipath component).

5. PHY Layer Performance Evaluation

This section presents results from our LTE-A PHY layer simulator [11]. In particular the performance of the DL SM MIMO OFDMA physical layer is investigated for NLoS UE terminals connected to horizontal and vertical MIMO arrays at the eNodeB. The 3D ray-tracing model introduced in Section 4 is used to compute the MIMO channel response matrices for each user. Results are presented in terms of the PER, mean throughput, and eigenvalue spread. Link selection is performed for each user, and a UE is assumed to be in outage if none of the MIMO-enabled MCS modes are able to meet the 10% PER threshold.

5.1. Throughput Performance Analysis

Figure 8 shows the usage (number of occurrences) of each modulation and coding scheme (assuming an SNR of 15 dB) from QPSK 1/2 (lowest data rate) to 64QAM 3/4 (highest data rate) when the eNB uses either a horizontal or vertical array. Antenna spacings of and are explored for picocells and and for macrocells. The graphs show that the higher-order MIMO modes are more likely to be selected in the picocellular environment. There are significantly less outage cases (denoted by “No MIMO Connection”) for the picocellular scenarios compared to the macrocellular scenarios.

Furthermore, for picocellular eNodeBs, the use of vertical arrays is now seen as a viable option. In the macrocells the usage of MIMO modes is seriously compromised by low angular spread in both the azimuth and elevation planes. In addition, the throughput occurrence rates for both the vertical and horizontal eNB antenna configurations are shown in Figures 9 and 10. It is interesting to note that the picocells offer significantly higher data rates compared to the macrocells. In both scenarios increasing the eNB antenna spacing has little overall effect in the picocells; however, this has a significant impact in the macrocells.

5.2. Packet Error Rate Performance Analysis

The PER performance for QPSK is compared for NLoS locations in a pico- and macrocell for different antenna spacings in Figure 11. The UE antenna separation in this analysis is maintained at . The PER performance is also shown for horizontal and vertical eNB antenna configurations. Perfect channel estimation is assumed along with a packet size of 54 bytes. All results are averaged over 1000 channel realisations unless otherwise stated. From Figure 11(a) it can be seen that that macroperformance using a vertical eNB array is very poor (around 23 dB worse than the equivalent picocellular channel). Figure 11(b) confirms that the vertical and horizontal eNB arrays offer similar performance in the picocell. However, as shown in Figure 11(c), it is clear that the vertical eNB configuration for a macrocell results in significantly degraded performance (relative to a traditional horizontal configuration). In the picocell, the horizontal configuration offers a 0.3 dB gain over the vertical configuration. However, for the same PER target in the macrocell, the vertical eNB antenna configuration is 6.5 dB worse than the horizontal case. Figures 11(b) and 11(c) present results for different eNB antenna separations. It can be seen that larger antenna separations improve the PER performance in both environments. As shown in Figure 11(b), the results for spacing in the picocell are around 2 dB worse than the use of spacings. Figure 11(c) shows the same scenario for the macrocell. The macrocell with eNB antenna spacing has an 8 dB gain over the case for horizontal orientations and 3 dB for vertical configurations.

5.3. Signal Correlation Performance Analysis

Figure 12 shows the relationship between the throughput (in Mbps) and the eigenvalue spread of the MIMO channel correlation matrix at a fixed SNR of 15 dB for the picocell and macrocell scenarios, respectively. High eigenvalue spreads indicate high spatial correlation in the eNB-UE links. Spatial correlation is directly linked to angle spread at the base station and UE terminal. In the picocells, where the elevation angle spread at the eNB is much closer to the azimuth values, the performance of vertically spaced arrays is similar to this of horizontally spaced arrays. However, in the macrocell the throughput performance for vertical element spacings at the eNB is much worse than for traditional horizontal spacings. Overall, we see that the higher throughput modes in LTE-A require lower values of eigenvalue spread.

Figure 13 shows the probability of the eigenvalue spread being less than the abscissa for horizontal and vertical eNB antenna configurations for various eNB antenna separations. It can be seen that vertical configurations degrade the eigenvalue spread in macrocells far more than in picocells. In all cases horizontal eNB antennas achieve better levels of spatial decorrelation (i.e., lower eigenvalue spreads). This result is particularly noticeable for the macrocell base stations since the elevation spreads in this scenario are significantly lower than the azimuth values.

6. Conclusions

This paper has combined polarimetric base station and UE 3D antenna patterns with state-of-the-art ray tracing to analyse the performance of horizontal and vertical MIMO arrays in an urban 2.6 GHz LTE-A network. Our propagation models were built on those used to guide the development of 3D channel model extensions in 3GPP. Results showed that for picocells the elevation angular spread at the base station was significantly increased relative to traditional microcells. A small increase was also observed here in the azimuth angle spread. Furthermore, our picocell studies demonstrated reduced values of RMS delay spread and Rician K-factor.

By combining the output of our detailed ray model with an LTE-A PHY layer simulator, it was possible to relate antenna and propagation statistics to PHY layer data rates and MIMO modulation and coding schemes. Results showed that vertical arrays at the eNodeB performed very poorly in macrocells but offered viable performance enhancement in picocells due to their increased levels of elevation angular spread.

Results showed that 3D angular statistics are necessary at the eNodeB and UE in order to adequately quantify the performance of LTE-A links. This paper has shown that picocells, which play a key capacity enhancing role in LTE HetNet configurations, are able to exploit vertical element spacings to deliver more compact antenna arrays at the base station.

Acknowledgments

Zuhanis Mansor would like to thank the Majlis Amanah Rakyat (MARA) and the Universiti Kuala Lumpur (UniKL) for her postgraduate scholarship and the Communication Systems and Networks Group, University of Bristol, for the provision of the simulation and laboratory facilities.