Atmospheric Effects on OFDM Wireless Links Operating In

electronics

Article Atmospheric Eﬀects on OFDM Wireless Links Operating in the Millimeter Wave Regime

Yosef Golovachev 1,* , Gad A. Pinhasi 2,* and Yosef Pinhasi 1,*

1 Department of Electrical and Electronic Engineering, Ariel University, Ariel 98603, Israel 2 Department of Chemical Engineering, Ariel University, Ariel 98603, Israel * Correspondence: [email protected] (Y.G.); [email protected] (G.A.P.); [email protected] (Y.P.)

Received: 6 September 2020; Accepted: 27 September 2020; Published: 29 September 2020

Abstract: The development of millimeter wave communication links and the allocation of bands within the Extremely High Frequency (EHF) range for the next generation cellular network present significant challenges due to the unique propagation effects emerging in this regime of frequencies. This includes susceptibility to amplitude and phase distortions caused by weather conditions. In the current paper, the widely used Orthogonal Division Frequency Multiplexing (OFDM) transmission scheme is tested for resilience against weather-induced attenuation and phase shifts, focusing on the effect of rainfall rates. Operating frequency bands, channel bandwidth, and other modulation parameters were selected according to the 3rd Generation Partnership Project (3GPP) Technical Specification. The performance and the quality of the wireless link is analyzed via constellation diagram and BER (Bit Error Rate) performance chart. Simulation results indicate that OFDM channel performance can be significantly improved by consideration of the local atmospheric conditions while decoding the information by the receiver demodulator. It is also demonstrated that monitoring the weather conditions and employing a corresponding phase compensation assist in the correction of signal distortions caused by the atmospheric dispersion, and consequently leads to a lower bit error rate.

Keywords: MMW; OFDM; 5G; 3GPP; atmospheric radio propagation

1. Introduction The rapid growing use of cellular services and the expanding line of products that are based on the wireless communications technologies have led to an increasing interest in the standardization of the 5th Generation (5G) mobile communication technology. The 5G technology is expected to use millimeter waves (MMW) to offer unprecedented spectrum and extremely high data rates to the wireless communications. The implementation of MMW technology has led to the recognition of the necessity of a consistent standardization of networks, channels, and modulation techniques that are relied on corresponding modeling efforts [1]. The US Federal Communications Commission (FCC) has stated their intention to establish the U.S. as leaders in the field of 5G [2]. One of the options considered for application was the allocation of bandwidth in the Extremely High Frequency (EHF) range, 30–300 GHz that can lead to the development of practical consumer applications of MMW technologies. The spectrum available for licensed use was set in three bands within the Ultra High Frequency (UHF) 850 MHz, 2.4 GHz, as well as MMW in the 28 GHz, 37 GHz, and 39 GHz bands. Furthermore, the 64–71 GHz band was recently designated for unlicensed use, which increases the total amount of unlicensed spectrum to 14 GHz as well as 600 MHz in the 37 GHz band for shared access. Different modulation techniques have been proposed for usage in 5G systems in order to provide greater flexibility and thereby increase transmission efficiency and reliability. Currently, the orthogonal

Electronics 2020, 9, 1598; doi:10.3390/electronics9101598 www.mdpi.com/journal/electronics Electronics 2020, 9, 1598 2 of 14 frequency-division multiplexing (OFDM) scheme is used in the 4G as a digital multi-carrier modulation method and is one of the popular techniques considered for the 5G cellular wireless communications [3]. The communication in the EHF range has many advantages; however, it is subjected to atmospheric propagation issues. In general, atmospheric effects are due to absorption and dispersion phenomena. In short wavelengths, the MMW propagation is significantly affected by weather conditions, such as air pressure, temperature, relative humidity, fog, and rain. Almost all of the studies on the effect of atmospheric conditions on the performance of communication links are referring to the attenuation only [4]. The signal phase shift and time delay effects are usually ignored, although they can be significant in extreme conditions [5,6]. Pinhasi et al. presented in their previous studies on MMW atmospheric propagation, theoretical and experimental results of absorptive and dispersive propagation effects and their influence on the signal strength and phase dispersion [7–11]. It was found that apart from attenuation, the dispersive effects like group delay should be considered in the design of MMW and THz radar and wireless communication links. In the present work, the performance of an OFDM wireless communication link operating in the MMW regime is studied under different weather conditions such as humidity, fog, and rain. The atmospheric effects on the OFDM waveform were formulated theoretically and demonstrated by simulation model. For this purpose, an atmospheric channel model for wireless communication was developed, combining our MMW propagation model with 5G Toolbox of MATLAB R2019b (MathWorks, Natick, MA, USA) [12]. MMW propagation model is applied to the transmitted baseband OFDM complex envelope and the performance of the wireless link is presented by using constellation diagram and BER (Bit Error Rate) performance chart. The present work is focused on the dispersive effects like group delay and phase shift and their effect on the OFDM communication link. The current paper presents an integration of several issues as sub-models: the MMW propagation model (Section2), the frequency dependent attenuation and phase dispersion factors for di fferent weather conditions (Section3), transmission of OFDM waveform in a medium (Section4), technical specifications of the next generation cellular network and the simulation parameters (Section5). Section6 presents simulations for MMW communication scenarios. The results of the simulation are presented in the form of a constellation diagram and BER performance graph for different atmospheric conditions, focusing on rainfall rate effect (Section7). The novel approach of predicting phase distortions due to the dispersion of the millimeter waves propagating in the atmospheric medium is presented in Section8. It is shown that by employing an adaptive phase compensator, the bit error rate can be considerably reduced.

2. Channel Model Besides modelling of the transmitter and receiver, analyzing the wireless link operating in the millimeter wave regime also requires accurate description of the medium. The transmitted signal ET(t) is propagating in the medium, which can be treated as a linear system via its impulse response h(t). The received signal is the result of the convolution ER(t) = ET(t) h(t). In the presence of Additive ∗ White Gaussian Noise (AWGN) n(t), the signal and noise obtained at the receiver’s input and fed to the demodulator is the summation Ein(t) = ER(t) + n(t). The link is illustrated in a block diagram in Figure1.

Figure 1. Block diagram of a wireless link operating in the millimeter wavelengths. Electronics 2020, 9, 1598 3 of 14

The atmospheric medium is characterized as a linear system with a transfer function H(j f ) given in the frequency domain. In the time domain, the transmitted signal ET(t) is presented as a carrier wave at a frequency fc modulated by a wideband complex envelope Aein(t) [8]:          j2π fct 1 +j2π fct 1 j2π fct ET(t) = Re[i(t) + jq(t)] e  = Aein(t) e + Aein∗(t) e− (1)  ·  2 · 2 · | {z }     Aein(t) 

where Aein(t) = i(t) + jq(t) is a complex envelope, given in terms of the In-phase i(t) and Quadrature q(t) baseband components. A Fourier transformation of Equation (1) is carried out to present the transmitted signal in the frequency domain:

+ R∞ j2π f t ET(j f ) = ET(t) ET(t) e dt = F ≡ · − (2) 1 −∞1 = Ae [j( f fc)] + Ae [ j( f + fc)] 2 in − 2 in∗ − n o where Ae (j f ) = Ae (t) . Using the transfer function H(j f ) of the medium, the received signal can in F in be now expressed in the frequency domain by:

ER(j f ) = ET(j f ) H(j f ) = 1 1 · (3) = Ae [j( f fc)] H(j f ) + Ae [ j( f + fc)] H(j f ) 2 in − · 2 in∗ − · Inverse Fourier transformation of Equation (3) results in the signal obtained at the receiver site:

+ R∞ +j2π f t ER(t) = ER(j f )e d f =    −∞         Z+    ∞   (4) = Re Ae (j f ) H[j( f + f )] e+j2π f td f  ej2π fct  in c    · ·  ·     −∞  | {z }    Aeout(t)

The expression for the received complex envelope is identiﬁed in Equation (4) to be:

+ Z ∞ +j2π f t Aeout(t) = Ae (j f ) H[j( f + fc)] e d f (5) in · · −∞ 3. Propagation Factors of the Atmosphere The millimeter wave propagation model (MPM) is used for calculation of the atmospheric frequency response taken into consideration of weather conditions. Contributions of dry air, water vapor, suspended water droplets (haze, fog, cloud), and rain are all addressed in the model. For completeness and clarity of notations, the model [7,9] is hereby reviewed. In general, the propagation of millimeter waves through the atmospheric medium along a distance d is described by the transfer function:

q 2π f j n( f ) d H(j f ) = TFS(d)e− c · · (6) Electronics 2020, 9, 1598 4 of 14 where n( f ) is the refractive index of the atmospheric medium and c is the speed of light in vacuum. TFS(d) is the free-space power transmission along a distance d excluding the atmospheric transmittance:

!2 c TFS(d) = GR GT (7) 2π fcd where GR and GT are the directivity gains of the receiving and transmitting antennas, respectively. The refractive index n( f ) is a frequency dependent quantity that can be expressed in terms of dispersive complex refractivity N( f ) (given in ppm):

6 n( f ) = 1 + [N0 + N0( f ) jN00 ( f )] 10 (8) | {z− } × N( f ) where the nondispersive part N0 is real and positive constant, N0( f ) and N00 ( f ) are the frequency dependent real and imaginary refractivity terms, respectively. The millimeter wave propagation model (MPM) is used for calculation of these refractivity terms in the frequency domain [13,14].

Figure 2. Attenuation factor (a,c,e) and incremental phase dispersion factor (b,d,f) in different weather conditions: (a,b) show the relative humidity RH [%] effect (c,d) show the effect of suspended droplets via the liquid water content W [g/m3] and (e,f) are for rainy conditions, where R is the rainfall rate in [mm/h]. Electronics 2020, 9, 1598 5 of 14

Substituting Equations (7) and (8) into Equation (6), the transfer function of the atmospheric medium can be presented as: q [α( f )+jβ( f )]d H(j f ) = TFS(d)e− (9) where α( f ) is the attenuation factor, given in terms of the imaginary part of the atmospheric refractivity:

2π f 6 α( f ) = N00 ( f ) 10− (10) c × where β( f ) is the wavenumber of the propagating ﬁeld, expressed in terms of the real part of the atmospheric refractivity:

2π f 2π f 6 β( f ) = + [N + N0( f )] 10− (11) c c 0 × | {z } ∆β( f )

The term ∆β( f ) is defined as the incremental phase dispersion caused by the dielectric substances composing the air. The frequency dependent attenuation and incremental phase dispersion are illustrated in Figure2, for different weather conditions characterized by relative humidity RH [%], suspended liquid water content W [g/m3], and rainfall rate R [mm/h]. The results reveal absorption peaks at 22, 183, and 325 GHz, where resonance absorption of water (H2O) occurs, as well as absorption peaks at 60 and 119 GHz, due to absorption resonances of oxygen (O2). Between these frequencies, minimum attenuation is obtained at 35 GHz (Ka-band), 94 GHz (W-band), 130 GHz (D-band), and 220 GHz, which are known as atmospheric transmission “windows”. Increasing the water content in the air in different forms (vapor (humidity), suspended water droplets (mist and fog), and drops (rain)) significantly increases the attenuation and incremental phase dispersion. The model was tested experimentally for various foggy conditions in a controlled environment demonstrating good agreement with measured results [9–11].

4. Transmission of OFDM Waveform in a Dispersive Medium The orthogonal frequency division multiplexing (OFDM) modulation technique is based on a digital multi-carrier signaling. A set of symbols are mapped onto adjacent sub-carriers in the frequency domain, employing an Inverse Fast Fourier Transformation [15]. A Cyclic Prefix (CP) is added to preserve orthogonality and reduce Inter-Symbol Interference (ISI) [16]. One of the main benefits of CP-OFDM is its simplicity, and its susceptibility against ISI and fading. One of the notable drawbacks of OFDM is its wide bandwidth. While preventing ISI, the use of a Cyclic Prefix leads to a reduction in spectral efficiency, as sending a redundant information. The use of CP-OFDM is also characterized by having high Out of Band Power, leading to possible Adjacent Channel Interference [17]. The transmitted complex envelope of the OFDM signaling is given by [7]:

KX1 − jk 2π t Aein(t) = [Im(k) + jQm(k)]e T (12) k=0 where K is the number of subcarriers, Im(k) and Qm(k) are the respective in-phase and quadrature components of the mth Quadrature Amplitude Modulation (QAM) symbol transmitted in subcarrier k, and T is symbol temporal duration. The Fourier transform of the transmitted complex envelope is:

K 1 ! X− k Ae (j f ) = [Im(k) + jQm(k)]δ f (13) in − T k=0 Electronics 2020, 9, 1598 6 of 14

Substitution of Equation (12) into Equation (5) results in the received complex envelope:

q KX1 ! − k jk 2π t Aeout(t) = T (d) [Im(k) + jQm(k)] H + fc e T (14) FS · T · k=0

Examination of Equation (14) reveals that the constellation of the received signal is aﬀected by the frequency response H(j f ) of the medium. In the case of a wireless link operating in millimeter waves, the resulted complex envelope is aﬀected by the dispersion of both attenuation and phase of the atmospheric refractivity:

q KX1 − [α( k + f )+jβ( k + f )]d jk 2π t Aeout(t) = T (d) [Im(k) + jQm(k)] e− T 0 T 0 e T (15) FS · · k=0

Inspection of Equation (15) for the complex envelope reveals the inﬂuences of absorptive and dispersive atmospheric propagation eﬀects on the constellation expected at the receiver site. It is noted that not only is attenuation resulted due to the atmospheric absorption, but also a frequency dependent phase shift that leads to a distorted constellation. The expression will be used to study the role of atmospheric refractivity in the determination of the performance of OFDM wireless link operating in millimeter wavelengths. In the following sections, we will focus on EHF bands allocated for the next generation of cellular networks—the 5G.

5. Technical Specifications of the Next Generation Cellular Network The definition of signal frame structure for the next generation cellular network was discussed in detail during the 3GPP and in many regulation forums [18–21]. Requirements throughout the RF specifications defined separately for different frequency ranges (FR). Two different frequency bands are available for the 5G technology, designated as: FR1-410 MHz to 7.125 GHz and FR2-24.25 GHz to 52.3 GHz [19,21,22]. The 5G New Radio (NR) Cellular Network is designed to operate in the FR2 operating bands presented in Table1. Both Uplink (UL) and Downlink (DL) bands operated in Time Division Duplexing (TDD).

Table 1. Frequency operating bands of New Radio (NR) in Ka-Band spectrum [21,22].

Uplink (UL) and Operating Band Carrier Frequency [GHz] Downlink (DL) Duplex Mode Operating Bands [GHz] Time Division n257 28 26.5–29.5 Duplexing (TDD) n258 26 24.25–27.5 TDD n260 39 37–40 TDD n261 28 27.5–28.35 TDD

In addition to FR1 and FR2, the Multiple Gigabit Wireless System (MGWS) was defined to operate in frequencies around 60 GHz. The MGWS network can be used in short-range, line-of-sight, and non-line-of-sight circumstances with traditional Wireless Local Area Network (WLAN) topologies [23]. Flexible Multi-Numerology are supported in an effort to present diverse services and user requirements of NR [24]. This numerology is characterized by different values of OFDM waveform parameters, such as SCS (Subcarrier Spacing) and CP (Cyclic Prefix) overhead. According to the 3GPP [22], the maximum channel bandwidth per NR carrier is 400 MHz and the maximum SCS is 480 kHz. At least for single numerology cases, the maximum number of subcarriers per band is Electronics 2020, 9, 1598 7 of 14

3300 or 6600 [25]. Frame Length is ﬁxed to be 10 ms. The number of OFDM symbols per slot time is depending on subcarrier spacing. For subcarrier spacing less than 60 kHz, the number of OFDM symbols per slot is 7 or 14. For subcarrier spacing greater than 60 kHz, the number of OFDM symbols per slot is 14 [26]. Table2 . presents the OFDM parameters designated for the 5G technology.

Table 2. Orthogonal frequency division multiplexing (OFDM) symbols conﬁguration in 5G NR [26].

Number Subcarrier Spacing (SCS) Number of Slots Number of OFDM SCS ∆f=2µ 15 [kHz] of Symbols Per Conﬁguration µ Per Frame Symbols Per Subframe · Slot 0 15 14 10 1 1 30 14 20 2 2 60 14 40 4 3 120 14 80 8 4 240 14 160 16

6. Channel Simulation Model The channel model was simulated by integrating the OFDM signaling with the medium frequency dependent transfer function resulted from the millimeter propagation model (MPM) for different weather condition scenarios. The received signal is decoded in the presence of AWGN by the OFDM demodulator. The flow chart shown in Figure3 summarizes the procedure employed in the numerical simulation calculating the resulted atmospheric absorptive and dispersive effects on the waveform constellation. Information related to the signaling designated for the next generation cellular network was taken from technical specifications presented above.

Figure 3. Flow chart of a baseband OFDM complex envelope evolution along the channel.

The model was implemented in MATLAB environment using 5G modeling Toolbox [12]. The model consists of the following steps:

1. Transmitted data signal settings: generation of transmission vector containing a given number (Nb = 24,000 or 50,400) of random bits in time domain. 2. Channel properties: the channel characteristics were calculated from the MPM by setting the carrier frequency (28, 39 or 60 GHz) and bandwidths (200 or 400 MHz). 3. Signaling: the modulation method for the current analysis is considered to be OFDM with the parameters presented in Table3: Symbol Rate, Cyclic Preﬁx Length, Guard Band Length, and Modulation Order. The modulation was built by creating complex QAM envelopes for each of the OFDM subcarriers. The modulator transmits a complex envelop waveform in accordance with Equation (12). 4. The medium: the atmospheric medium was modeled by the MPM to generate the transfer function (9) in the frequency domain for a given distance (200 m) and diﬀerent weather condition scenarios. The meteorological parameters were set to: temperature (25 ◦C), pressure (1 atm), relative humidity (RH = 0, 10, 50, 100%) suspended liquid water content (W = 0, 0.1, 1, 5 g/m3), or rain rate (R = 0, 5, 10, 15 mm/h). Note that for the case of suspended liquid water droplets and rain the relative humidity is set to be 100%. Electronics 2020, 9, 1598 8 of 14

5. Noise: Additive White Gaussian Noise (AWGN) with a given Signal-to-noise ratio (SNR) was introduced to the signal. 6. Demodulator: after passing through the atmospheric medium, the received signal is demodulated to decode the digital data.

Table 3. OFDM parameters for the simulation.

Parameter Units Value Channel Bandwidth [MHz] 200 Number of Subcarriers - 3300 OFDM Subcarrier Spacing [kHz] 60 Subcarrier Spacing [kHz] 15 Quadrature Amplitude Modulation (QAM) Modulation Order - 6 OFDM Symbols Rate [kSymbols/sec] 14 Bits Per Subcarrier - 8 Cyclic Preﬁx Length - 64

Figure 4. Simulation results of 64 QAM constellation diagram and eﬀect of Phase Compensator on it at carrier frequencies: (a,b) 28, (c,d) 39, (e,f) 60 GHz for 200 m wireless link in diﬀerent rainfall rate. Electronics 2020, 9, 1598 9 of 14

The channel performance was studied to ﬁnd the constellation diagram at the receiver and the resulted Bit Error Rate (BER). The BER is calculated by comparing the received bit sequence with the transmitted one. The quality of the reception is also measured via the constellation diagram in the signal vector I/Q space. Each point in the constellation diagram represents a given QAM symbol as a vector. The constellation diagram demonstrates amplitude and phase distortions emerging due to the dispersive propagation as well as the additive noise associating with the demodulation.

7. Simulation Results The performance of the OFDM wireless link under various weather conditions was studied using the developed model. A single transmitter and receiver (SISO) transmission link was simulated, within bandwidth of 200 MHz, centered at 28 GHz, 39 GHz, or 60 GHz. Typical transmission distance of 200 m is assumed. Simulations were carried out for reference conditions: an ambient temperature of 25 ◦C and air pressure of 101 kPa for an urban environment. Using meteorological parameters, the MPM predicts frequency-dependent absorptive and dispersive phenomena, causing distortions in amplitude and phase of propagated waveform.

Figure 5. Bit Error Rate (BER) Performance of 64 QAM and eﬀect of Phase Distortion Compensator on it at carrier frequencies of (a,b) 28, (c,d) 39, and (e,f) 60 GHz for 200 m wireless link in diﬀerent rainfall rate. Electronics 2020, 9, 1598 10 of 14

In the current work, the presented results are focused on the effect of rain, where the most significant effect on the link budget was identified. The propagation parameters calculated by the MPM model are summarized in Table3 for the designated 5G bands. The attenuation and phase shift are revealed to increase with the rain rate and with the carrier frequency. These results correspond with experimental climatic wind tunnel results [4]. At the second stage, the effect of the rain intensity on the quality of waveform reception and OFDM demodulation is examined. The bit error ratio (BER) as function of SNR at carrier frequencies 28, 39, and 60 GHz of a wireless link of 200 m is presented in Figure4a,c,e for 64-QAM constellation. For each frequency, the BER is found for different rainy scenarios: no rain, 5, 10, and 15 mm/h. It is important to note that the rain plays a different role at different frequency bands. For instance, setting the carrier frequency to 28 GHz, SNR = 25 dB, the resulted BER is increased from 2 10 5 (no rain) × − to 2 10 4 in rainy conditions (when the rain intensity is R = 15 mm/h). However, for 39 GHz carrier × − frequency, SNR = 25 dB and rain intensity of R = 15 mm/h, the BER further increases to 2 10 3 × − (see Figure5a,c). The received signal quality is analyzed by displaying the constellation diagram for each scenario. A contraction of the constellation diagram demonstrates attenuation of the OFDM signal, while rotation of the diagram is a result of phase shift. Figure4a,c,e shows the e ffect of rain on 64 QAM constellation diagram at carrier frequencies of 28, 39, and 60 GHz of wireless link of 200 m. The attenuation expressed by α( f ) causes the constellation to “shrink”. The revealed counterclockwise rotation of the diagram is a result of dispersion of the propagating OFDM signal due to the rain. The rotation angle increases with frequency corresponding to the incremental phase dispersion, given in term of ∆β( f ) in Equation (11). Both effects cause the symbols to change their position in the constellation I/Q domain resulting in misinterpretation by the demodulators and erroneous symbol receptions. While constellation diagrams are useful for graphically visualizing signal data, it is also useful to quantify the disparity between received signal waveform and the ideal constellation by calculating the error vector magnitude (EVM). EVM is the scalar distance between the end points of the measured and ideal phasors and is a measure of how well a digital communications system is performing [27]. As is presented in Figure6, the higher the QAM signaling order is, the resulted constellation rotation angle (Table4) at the receiver is more expressed. Therefore, the EVM will increase for high modulation orders, resulting in erroneous detected symbols even at relatively low rainfall intensities.

Figure 6. Comparison of 64 QAM and 256 QAM constellation diagrams at a carrier frequency of 28 GHz for 200 m wireless link for diﬀerent rainfall rates. Electronics 2020, 9, 1598 11 of 14

Table 4. Propagation factors of 28, 39, and 60 GHz for 200 m wireless link.

Parameter Units No Rain 5 mm/h 10 mm/h 15 mm/h 28 GHz Attenuation Factor [dB/km] 0.3074 1.0460 1.8501 2.6809 Phase Factor [deg/km] 1.7740 3.9402 6.5050 9.5047 Attenuation for 200 m [dB] 0.0615 0.2092 0.3700 0.5362 Incremental part of phase for 200 m [deg] 0.3548 0.7880 1.3010 1.9009 39 GHz Attenuation Factor [dB/km] 0.3310 1.7524 3.1170 4.4610 Phase Factor [deg/km] 2.9819 2.6178 9.0821 16.4558 Attenuation for 200 m [dB] 0.0662 0.3505 0.6234 0.8922 Incremental part of phase for 200 m [deg] 0.5964 0.5236 1.8164 3.2912 60 GHz Attenuation Factor [dB/km] 14.63 17.5066 19.8636 22.0577 Phase Factor [deg/km] 3.3368 11.9821 28.9276 47.4929 Attenuation for 200 m [dB] 2.926 3.5013 3.9727 4.4115 Incremental part of phase for 200 m [deg] 0.6674 2.3964 5.7855 9.4986

8. Channel Performance Improvement via Phase Distortion Compensation As mentioned in Section3, the e ffect of the atmosphere on the communication channel is expressed by the attenuation coefficient α( f ) and by the incremental, frequency dependent phase dispersion ∆β( f ), given in Equations (10) and (11), respectively. Figure4a,c,e presents the e ffect of these two parameters on the modulated signal, revealing respective “shrink” and “rotation” in the constellation. While the effect of the attenuation α( f ) is usually considered in an appropriate design of the link budget, we suggest to neutralize the phase dispersion ∆β( f ) caused by the atmospheric medium by adding a phase compensator in the pre-demodulating stage in the receiver MMW chain (see Figure7). The frequency dependent phase dispersion coefficient can be estimated by meteorological conditions at the receiver site, including relative humidity RH [%], suspended liquid water content W [g/m3], and rainfall rate R [mm/h].

Figure 7. Flow chart of millimeter wave (MMW) link including Phase Distortion Compensator in the receiver site.

Figure4b,d,f shows a positive e ffect of the Phase Compensator on a distorted 64 QAM constellation diagram for different rainy scenarios: no rain, 5, 10, and 15 mm/h. As it presented, Phase Distortion Compensator completely neutralizes the “rotation’” effect dramatically. As a result, the receiver demodulator becomes able to reproduce the data almost as in no-rain scenarios (see Figure5b,d,f). Table5 summarizes the expected performances of the average EVM and BER when an adaptive phase compensator is installed before the OFDM demodulation. A comparison with the performances attained when no compensating mechanism is applied during rainy conditions demonstrates that the proposed novel approach improves the reliability of the channel, leading to a reduction in the EVM Electronics 2020, 9, 1598 12 of 14 and consequently lowering the error rate. A quantitative measure of the improvement is expressed via the average EVM relative reduction (in [%]):

EVMw/o comp EVMw comp − − − 100% (16) EVMw/o comp · −

Table 5. Average EVM and BER improvements due to Phase Dispersion Compensation for 200 m wireless link operating with 64QAM and Signal-to-noise ratio (SNR) = 25 dB.

Without Compensator With Compensator Frequency Rainfall Improvement Average EVM Average EVM [GHz] Rate [mm/h] BER BER Factor [%] EVMw/o comp EVMw comp − − 5 0.170 2.817 10 5 0.165 1.210 10 5 2.99 × − × − 28 10 0.364 7.044 10 5 0.335 1.548 10 5 7.85 × − × − 15 0.567 2.137 10 4 0.507 1.667 10 5 10.61 × − × − 5 0.346 2.183 10 5 0.301 1.151 10 5 12.79 × − × − 39 10 0.687 1.756 10 4 0.573 1.290 10 5 16.58 × − × − 15 1.028 1.996 10 3 0.839 1.369 10 5 18.39 × − × − 5 0.514 5.335 10 4 0.397 1.349 10 5 22.69 × − × − 60 10 0.953 2.994 10 2 0.694 1.528 10 5 27.19 × − × − 15 1.376 1.183 10 1 0.961 1.587 10 5 30.11 × − × −

9. Conclusions The transmission quality of the 5th generation (5G) millimeter wave (MMW) wireless communication was studied for various atmospheric weather conditions. The current work is focused on the investigation of the attenuation and phase shift introduced by weather effects, especially the influence of different rainy scenarios. The link reliability was examined via constellation diagram and BER in the presence of AWGN for different rainfall rates. Simulation results indicate that OFDM channel performance is strongly affected by rainy conditions. It is shown that absorptive and dispersive behavior of the medium transfer function results in distorted magnitude and phase QAM signaling, leading to an increased error rate. The OFDM MMW channel performance can be significantly improved by considering the effects emerging due to local atmospheric conditions in the receiver demodulator scheme. Online weather reports from local meteorological stations can provide side information on the propagation characteristics that will enable compensation of waveform magnitude and phase distortions at the receiver. This will significantly reduce the error rate leading to an improvement of the link availability and reliability. Furthermore, the on-line meteorological data could also be useful for selecting the optimal parameters for the QAM constellation and the appropriate OFDM modulation numerology to increase network efficiency. It is demonstrated that by employing an adaptive, pre-demodulation phase correction scheme that the atmospheric dispersion phenomena on the millimeter wave can be compensated, resulting in an improvement of the link reliability, as summarized in Table5. It is important to note that by using Phase Dispersion Compensator in the pre-demodulation stage of the receiver, the reliability of the channel during rainy conditions is improved. It is shown that for all examined scenarios, even for extremely bad weather conditions, the compensation revealed a reduction in the error probability, almost to the level expected in clear weather (no rain) scenarios (see Figure5). Electronics 2020, 9, 1598 13 of 14

The necessity of phase compensation depends on the link operational frequency, the signal to noise ratio, as well as on the rain intensity. Major improvements are demonstrated in high rain rates, when the weather conditions cause severe phase distortions. It is also shown that for frequency bands, where highly absorptive and dispersive atmospheric eﬀects are expected (as in 60 GHz), the compensator modify the phase distortions and as a result reduced EVM and lower error probability.

Author Contributions: Conceptualization, Y.G., G.A.P. and Y.P.; Data curation, Y.G.; Formal analysis, Y.G., G.A.P. and Y.P.; Investigation, Y.G., G.A.P. and Y.P.; Methodology, G.A.P. and Y.P.; Project administration, Y.P.; Software, Y.G. and G.A.P.; Supervision, G.A.P. and Y.P.; Validation, Y.G.; Visualization, Y.G.; Writing—original draft, Y.G.; Writing—review & editing, G.A.P. and Y.P. All authors have read and agreed to the published version of the manuscript. Funding: This research received no external funding. Conﬂicts of Interest: The authors declare no conﬂict of interest.

References

1. Rappaport, T.S.; Xing, Y.; MacCartney, G.R.; Molisch, A.F.; Mellios, E.; Zhang, J. Overview of Millimeter Wave Communications for Fifth-Generation (5G) Wireless Networks—With a Focus on Propagation Models. IEEE Trans. Antennas Propag. 2017, 65, 6213–6230. [CrossRef] 2. Federal Communications Commission. Fact Sheet: Spectrum Frontiers Proposal to Identify, Open Up Vast Amounts of New High-Band Spectrum for Next Generation (5G) Wireless Broadband. Available online: https://www.fcc.gov/document/fact-sheet-spectrum-frontiers-item (accessed on 7 June 2020). 3. Lathi, B.P.; Ding, Z. Modern Digital and Analog Communication Systems, 5th ed.; Oxford University Press: New York City, NY, USA, 2019. 4. Hilt, A. Availability and Fade Margin Calculations for 5G Microwave and Millimeter-Wave Anyhaul Links. Appl. Sci. 2019, 9, 5240. [CrossRef] 5. Antes, J.; Boes, F.; Messinger, T.; Lewark, U.J.; Mahler, T.; Tessmann, A.; Henneberger, R.; Zwick, T.; Kallfass, I. Multi-Gigabit Millimeter-Wave Wireless Communication in Realistic Transmission Environments. IEEE Trans. Terahertz Sci. Technol. 2015, 5, 1078–1087. [CrossRef] 6. Matricciani, E. A Relationship Between Phase Delay and Attenuation Due to Rain and Its Applications to Satellite and Deep-Space Tracking. IEEE Trans. Antennas Propag. 2009, 57, 3602–3611. [CrossRef] 7. Pinhasi, Y.; Yahalom, A.; Harpaz, O.; Vilner, G. Study of ultra-wideband transmission in the extremely high frequency (EHF) band. IEEE Trans. Antennas Propag. 2004, 52, 2833–2842. [CrossRef] 8. Pinhasi, Y.; Yahalom, A.; Pinhasi, G.A.; Lutvak, M. Ultra Wideband Wireless Communications in the 60 GHz Band. WSEAS Transact. Comm. 2010, 12, 773–781. 9. Golovachev, Y.; Etinger, A.; Pinhasi, G.A.; Pinhasi, Y. Millimeter Wave High Resolution Radar Accuracy in Fog Conditions—Theory and Experimental Verification. Sensors 2018, 18, 2148. [CrossRef][PubMed] 10. Golovachev, Y.; Etinger, A.; Pinhasi, G.A.; Pinhasi, Y. Propagation properties of sub-millimeter waves in foggy conditions. J. Appl. Phys. 2019, 125, 151612. [CrossRef] 11. Etinger, A.; Golovachev, Y.; Shoshanim, O.; Pinhasi, G.A.; Pinhasi, Y. Experimental Study of Fog and Suspended Water Effects on the 5G Millimeter Wave Communication Channel. Electronics 2020, 9, 720. [CrossRef] 12. 5G Toolbox MATLAB, MathWorks Inc. Available online: https://www.mathworks.com/products/5g.html (accessed on 7 June 2020). 13. Liebe, H.J. MPM? An atmospheric millimeter-wave propagation model. Int. J. Infrared Millim. Waves 1989, 10, 631–650. [CrossRef] 14. Liebe, H.J.; Hufford, G.A.; Manabe, T. A model for the complex permittivity of water at frequencies below 1 THz. Int. J. Infrared Millim. Waves 1991, 12, 659–675. [CrossRef] 15. Gerzaguet, R.; Bartzoudis, N.G.; Baltar, L.G.; Berg, V.; Dore, J.-B.; Kténas, D.; Font-Bach, O.; Mestre, X.; Payaró, M.; Färber, M.; et al. The 5G candidate waveform race: A comparison of complexity and performance. EURASIP J. Wirel. Commun. Netw. 2017, 2017, 1312. [CrossRef] Electronics 2020, 9, 1598 14 of 14

16. An, C.; Kim, B.; Ryu, H.-G. Design and evaluation of spectrum efficient WR-OFDM system for 5G and B5G mobile system. In Proceedings of the 2017 IEEE International Conference on Microwaves, Antennas, Communications and Electronic Systems (COMCAS), Tel Aviv, Israel, 13–15 November 2017; pp. 1–5. 17. Levy, D.; Reichman, A. Filter bank multi carrier modulation performance. In Proceedings of the 2017 IEEE International Conference on Microwaves, Antennas, Communications and Electronic Systems (COMCAS), Tel Aviv, Israel, 13–15 November 2017; pp. 1–6. 18. Mourad, A.; Yang, R.; Lehne, P.H.; De La Oliva, A. A Baseline Roadmap for Advanced Wireless Research Beyond 5G. Electronics 2020, 9, 351. [CrossRef] 19. Series, M. IMT Vision—Framework and Overall Objectives of the Future Development of IMT for 2020 and Beyond. In Recommendation ITU; Electronic Publication: Geneva, Switzerland, 2015. 20. Guey, J.-C.; Liao, P.-K.; Chen, Y.-S.; Hsu, A.; Hwang, C.-H.; Lin, G. On 5G radio access architecture and technology [Industry Perspectives]. IEEE Wirel. Commun. 2015, 22, 2–5. [CrossRef] 21. Gozalvez, J. New 3GPP Standard for IoT [Mobile Radio]. IEEE Veh. Technol. Mag. 2016, 11, 14–20. [CrossRef] 22. 3GPP Technical Specification 38.101-2. NR. Available online: https://www.3gpp.org/dynareport/38101-2.htm (accessed on 7 June 2020). 23. ITU-R. M.2003-2, Multiple Gigabit Wireless Systems in frequencies around 60 GHz. Available online: https://www.itu.int/rec/R-REC-M.2003-2-201801-I/en (accessed on 7 June 2020). 24. Yazar, A. Flexible Multi-Numerology Systems for 5G New Radio. J. Mob. Multimed. 2018, 14, 367–394. [CrossRef] 25. 3GPP Technical Report TR 38.802. Available online: https://www.3gpp.org/dynareport/38802.htm (accessed on 7 June 2020). 26. 3GPP Technical Specification TS 38.211. NR. Available online: https://www.3gpp.org/dynareport/38211.htm (accessed on 7 June 2020). 27. Keysight Technologies, Digital Modulation in Communication Systems—An Introduction. Available online: https://literature.cdn.keysight.com/litweb/pdf/5965-7160E.pdf (accessed on 7 June 2020).

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).