Microsoft Word - 233 completed IIUM Engineering Journal, Vol. 12, No. 5, 2011: Special Issue -1 on Science and Ethics in Engineering Dao et al. 53 MODIFICATION OF ITU-R RAIN FADE SLOPE PREDICTION MODEL BASED ON SATELLITE DATA MEASURED AT HIGH ELEVATION ANGLE HASSAN DAO, MD. RAFIQUL ISLAM AND KHALID AL-KHATEEB Electrical and Computer Engineering Department, Faculty of Engineering, International Islamic University Malaysia, Jalan Gombak, 53100 Kuala Lumpur, Malaysia. rafiq@iium.edu.my ABSTRACT: Rain fade slope is one of fade dynamics behaviour used by system engineers to design fade mitigation techniques (FMT) for space-earth microwave links. Recent measurements found that fade slope prediction model proposed by ITU-R is unable to predict fade slope distribution accurately in tropical regions. Rain fade measurement was conducted in Kuala Lumpur (3.3� N, 101.7� E) where located in heavy rain zone by receiving signal at 10.982 GHz (Ku-band) from MEASAT3 (91.5� E) on 77.4� elevation angle. The measurement has been carried out for one year period. Fade slope S parameter on ITU-R prediction model has been investigated. New parameter is proposed for the fade slope prediction modeling based on measured data at high elevation angle, Ku-band. ABSTRAK: Cerun hujan pudar adalah salah satu dinamik tingkah laku pudar yang digunakan oleh jurutera sistem untuk mereka bentuk teknik-teknik pengurangan pudar (FMT) bagi link gelombang mikro ruang bumi. Ukuran baru-baru ini mendapati bahawa cerun pudar ramalan model yang dicadangkan oleh ITU-R tidak mampu untuk meramalkan pengagihan cerun pudar tepat di kawasan tropika. Pengukuran hujan pudar telah dijalankan di Kuala Lumpur (3.3� N, 101.7� E) yang terletak di kawasan hujan lebat dengan menerima isyarat pada 10,982 GHz (Ku-band) dari MEASAT3 (91.5� E) pada sudut ketinggian 77.4�. Pengukuran telah dijalankan untuk tempoh satu tahun. Parameter cerun pudar S pada model ramalan ITU-R telah disiasat. Parameter baru adalah dicadangkan untuk pemodelan cerun pudar ramalan berdasarkan data yang diukur pada sudut paras ketinggian, Ku-band. KEYWORDS: fade slope; ITU-R; fade mitigation techniques; sampling time interval 1. INTRODUCTION In satellite communication on frequencies above 10 GHz, rain attenuation is great propagation phenomenon that influence large variations in received signal levels. Rain attenuaton is inordinately degraded received signal. That is why many mitigating approachs have been developed to reduce fading outage of communication links. Rain attenuation behaviours are extensively studied and measured but dynamic aspect of rain attenuation still has sparse study. Knowledge on rain fade dynamics is noticeable important in the design of new satellite communication systems. Adaptive power control (APC) and Uplink power control (ULPC) are ones of rain countermeasure techniques relies on rain fade dynamics parameters by increasing transmitted power to compensate rain fading on propagation path. fade slope is one behaviour of fade dynamics IIUM Engineering Journal, Vol. 12, No. 5, 2011: Special Issue characteristics that assesses tracking speed contributing factor to fade mitigation techniques (FMT) [1, 2] Various researchers have studied and proposed parameters and prediction models for rain fade slope. Most studies have based on temperate region measurements [1, 3 R P.1623 [8] has recommended a prediction model for fade slope base model [3, 4] measured at Netherlands receiving signal from Olympus satellite ( 26.78˚) on 12.5 - 30 GHz. The model was examined using data from other sites in the UK, France, Belgium, Italy, and the US [1]. However, some measurements [ have indicated the model is required to zones. This paper analyses rain fade slopes comparing between available model slope S parameters has also been focused and proposed on high elevation angle as a function of different fade levels 2. EXPERIMENTAL SETUP Experimental system has been monitored signal from MEASAT 3 (geostationar E) at 10.982 GHz frequency (Ku modulation provided 15 channels for Malaysian direct broadcast satellite (DBS) pay television service (Astro). The received signal level has been collected for 12 months perio Table 1 have summarized experimental setup and system parameters, respectively at Satellite Communication lab, IIUM. Fig. 1: Experimental setup for measurement at Ku Satellite Communications Lab, IIUM. Table Ground Station Location Beacon Frequency Satellite Position Polarization Azimuth Elevation Angle Antenna Diameter NF of LNB g Journal, Vol. 12, No. 5, 2011: Special Issue -1 on Science and Ethics 54 characteristics that assesses tracking speed contributing factor to fade mitigation Various researchers have studied and proposed parameters and prediction models for rain fade slope. Most studies have based on temperate region measurements [1, 3 R P.1623 [8] has recommended a prediction model for fade slope based on Van De Kamp model [3, 4] measured at Netherlands receiving signal from Olympus satellite ( 30 GHz. The model was examined using data from other sites in the UK, France, Belgium, Italy, and the US [1]. However, some measurements [9] in tropical zones have indicated the model is required to improve in order to correspond to different climate rain fade slopes characteristic on different sampling times and available model and measurement in Kuala Lumpur, Malaysia. parameters has also been focused and proposed on high elevation angle as a different fade levels. EXPERIMENTAL SETUP Experimental system has been monitored signal from MEASAT 3 (geostationar E) at 10.982 GHz frequency (Ku-Band), vertical polarization with QPSK modulation provided 15 channels for Malaysian direct broadcast satellite (DBS) pay television service (Astro). The received signal level has been collected for 12 months have summarized experimental setup and system parameters, respectively at te Communication lab, IIUM. Experimental setup for measurement at Ku-band signal over earth- Satellite Communications Lab, IIUM. Table 1: System parameters of the link. Ground Station Location 3.3°N, 101.7°E Beacon Frequency 10. 982 GHz (IF = 1.232GHz ) Satellite Position 91.5° E (MEASAT 3) Vertical 253° Elevation Angle 77.4° Antenna Diameter 2.4 m 0.3 dB 1 on Science and Ethics in Engineering Dao et al. characteristics that assesses tracking speed contributing factor to fade mitigation Various researchers have studied and proposed parameters and prediction models for rain fade slope. Most studies have based on temperate region measurements [1, 3-7]. ITU- d on Van De Kamp model [3, 4] measured at Netherlands receiving signal from Olympus satellite (El = 30 GHz. The model was examined using data from other sites in the UK, 9] in tropical zones improve in order to correspond to different climate characteristic on different sampling times and measurement in Kuala Lumpur, Malaysia. Fade parameters has also been focused and proposed on high elevation angle as a Experimental system has been monitored signal from MEASAT 3 (geostationar Band), vertical polarization with QPSK modulation provided 15 channels for Malaysian direct broadcast satellite (DBS) pay television service (Astro). The received signal level has been collected for 12 months . Fig. 1 and have summarized experimental setup and system parameters, respectively at -space path at IIUM Engineering Journal, Vol. 12, No. 5, 2011: Special Issue -1 on Science and Ethics in Engineering Dao et al. 55 Signal from MEASAT3 has been received by 2.4 m diameter of parabolic dish antenna with Low noise block (LNB), 0.3 dB noise figure. The signal is down converted to Intermediate Frequency (IF) as L-band with Local Oscillator (LO) 975 MHz. The subsequent IF is sent to spectrum analyzer setting post-detection bandwidth (video bandwidth) of 10 Hz [10] to monitor carrier signal and remove spurious noise. Consequent Output signal from spectrum analyzer is logged and saved in databank by using software developed with sampling time interval of 1s. In addition, rainfall rate is also measured synchronously by 0.2mm-Tipping bucket rain gauge every 10s time interval. 3. FADE SLOPE MODEL PROPOSED BY ITU-R ITU-R has proposed statistical distribution model of fade slope developed by Van De Kamp [1, 11] that depends on cutoff frequency, attenuation and sampling time. The proposed prediction model of probability density function (PDF) is shown Eq. 1, ( ) ( )( ) ( ) 1 22 / 1 2 | − + = sdBAp ξς σξπσ ξ (1) Where ξσ is the standard deviation of the conditional fade slope at a given attenuation level given by, dB/s),( AtfFs B ∆=ζσ (2) Where ( ) bbbB B tf tfF /1 2 )2(/1 2 ),( ∆+ π =∆ (3) With b = 2.3, fB is the 3 dB cutoff frequency of the low pass filter (Hz), and S parameter is from experiment as derived above. 4. DATA PRE-PROCESSING 4.1 Clear sky reference level Finding out zero dB reference level from clear sky on measurement of satellite signal attenuation is a fundamental trouble requirement. Defining clear sky reference level by hourly averaging after and before raining days [12]. Raining events were determined on 10s interval time for tipping bucket type of rain gauge. The number of data points (second) during raining events used to generate fade slope for each attenuation values are shown in Fig.3. It is obvious that the number of occurrences decreases with attenuation levels. Fig. 3: Number of time series (second) extracted from data measured from July 1, 2010 to June 30, 2011. IIUM Engineering Journal, Vol. 12, No. 5, 2011: Special Issue -1 on Science and Ethics in Engineering Dao et al. 56 4.2 Rain Attenuation and Scintillation Extraction Attenuation time series during raining events are contaminated by fast fluctuations due to tropospheric scintillation. Therefore, tropospheric scintillation is essentialy decimated prior to proceed next processes [13]. Power spectral density (PSD) is efficient technique in order to separate propagation due to scintillation and rain attenuation by specifying appropriate cutoff frequency (fc) [14, 15]. High pass filter (HPF) order 6 th of Digital Butterworth filter with cutoff frequency of 0.025 Hz [16] obtained by the PSD function and has been implemented to remove fast fluctuation due to scintillation. 5. DATA ANALYSIS AND RESULTS Fade slopes are calculated for each attenuation level thresholds after tropospheric scintillation has been decimated. Each attenuation level (dB), A is calculated when dBAtAA 5.0)(5.0 +≤<− , then fade slope is calculated by Eq. (4) ( ) ( ) ( ) ( )sdB t tiAtiA t ∆ ∆−−∆+ = 2 δ (4) Where ∆t is sampling time, 1s and i is sample number. 5.1 Comparison Between Measured Data and Prediction Model Once fast fluctuation due to tropospheric scintillation has been decimated from the signal, fade slope has been calculated from Eq. 4. The prediction model has been calculated by Eq. 1 using fB 25 mHz obtained from the experiment, ∆t 1s and S = 0.01 (ITU-R proposed). Comparison of measurement and model has been depicted in Fig. 4. Fig. 4: PDF comparison between fitting of measured data (solid line) and ITU-R (dash line) model for 1s sampling time. 5.2 Investigation of Standard Deviation Rain fade slope characteristic has been investigated to evaluate the parameters for ITU-R [8] based on the measured data. Principal parameters analysis was focused on standard deviation (σζ) with different attenuation thresholds. The measured standard deviation and the model have been presented as a function of attenuation thresholds in Fig. 5 (a). Measured σζ indicated that keeps increasing with attenuations from 0.0038 dB/s and reaching maximum 0.018 dB/s at 6dB attenuation and subsequently reducing until 0.0079 IIUM Engineering Journal, Vol. 12, No. 5, 2011: Special Issue -1 on Science and Ethics in Engineering Dao et al. 57 dB/s at 10 dB. Standard deviation, σζ was estimated by ITU-R prediction model indicated that σζ increased linearly with attenuation levels. (a) (b) Fig. 5: (a)Comparison of standard deviations between measured data (solid line) and predicted by ITU-R prediction models (dash line), (b) The measured value of S as a function of attenuations between measurement (solid line) and Zero Offset fit (ZOF) (dash line). 5.3 Proposed Standard Deviation ITU-R prediction model has been proposed σζ which depends on fB, ∆t and parameter S. Parameter S is a constant which depends on climate and elevation angle which affect to σζ of ITU-R model. ITU-R has proposed the constant averagely at 0.01 in Europe and USA at elevation angle between 10�-60�. But there is no model proposed in tropical region and on high elevation angle. Therefore, S is needed to derive in order to propose for the model on this region and high elevation angle. Parameter S is obtained by Eq. (5) using σζ from the experiment and has been presented by Fig. 5(b). AtfF S B ),( ∆ = δ σ (5) Thus, determination of empirical S is implemented by Zero offset Fit (ZOF) [2, 3, 4] yielded 0.0023 to propose in the model for this characteristics. 5.4 Comparison with the Model Statistical distribution comparisons between measurement and model with different parameter S have been depicted in Fig. 6(a) and 6(b) as a function of attenuations. Figure 6 illustrates fade slope distribution comparison between measured data and the propoaed model. The models are drawn S parameters as applying Zero offset fit (ZOF) method. PDF’s predicted by ZOF is closer fit with the measurement for attenuation between 1dB- 6dB. For attenuation higher than 6dB, the model shows wider with attenuations while the measurement showd narrower with attenuations corresponding to standard deviations. Table 2: Descriptive PDF Statistics comparison between measurement and models. Attenuation (dB) 1 3 6 8 10 Mean Measured 3.311 3.308 3.306 3.31 3.311 ZOF 3.357 3.311 3.31 3.309 3.307 SD Measured 17.42 13.181 7.785 8.747 12.548 ZOF 25.123 13.626 9.345 7.921 6.929 Skewness Measured 6.403 5.24 2.859 3.19 4.702 ZOF 9.864 5.105 3.419 2.85 2.449 Kurtosis Measured 43.185 28.818 7.413 9.573 22.645 ZOF 105.662 26.963 11.22 7.316 4.996 IIUM Engineering Journal, Vol. 12, No. 5, 2011: Special Issue -1 on Science and Ethics in Engineering Dao et al. 58 (a) (b) Fig. 6: PDF comparison among Measurement and model with different values of s- parameters of (a)1dB to 5dB attenuations and (b) 6dB to 10dB attenuations levels. 6. CONCLUSION Rain fade was measured in Kuala Lumpur (3.3 0 N, 101.7 0 E) at 10.982 GHz (Ku- band) from MEASAT3 (91.5 0 E) on 77.4 0 elevation angle. The measurement has been carried out for one year period. Fade slope S parameter on ITU-R prediction model has been investigated based on data processed with 1 sec integration time. The probability distribution functions proposed by ITU-R with parameter S = 0.01 are much wider than the measurement for all attenuation levels. The value of S parameter is evaluated using implemention of Zero offset Fit (ZOF) with measured data. The new value of S = 0.0023 is proposed in the fade slope prediction model for high elevation angle, Ku-band in tropical regions. ACKNOWLEDGEMENT The authors are very grateful to the Research Management Centre (RMC), International Islamic University Malaysia (IIUM) for research grant supporting this project. REFERENCES [1] L. Castanet and M.V. de Kamp, “Modelling the Dynamic properties of the propagation channel”, 5th Management Committee Meeting of the COST 280 Action, May 2003 [2] Chambers, A.P., Callaghan S.A. and Otung, I.E., “Analysis of rain fade slope for Ka and V- Band satellite links in southern England”, IEEE Transac. On Antenna and Propagation, Vol.54, No.5, May 2006. [3] M.J.L. Van de Kamp, “Statistical Analysis of Rain Fade Slope”, IEEE Transac. On Antenna and Propagation, Vol.51, No.8, August 2003 [4] M.M.J.L. Van de Kamp, “Climatic Radiowave Propagation Models for the design of Satellite Communication Systems”, Ph.D. Thesis, Eindhoven University of Technology, Eindhoven, the Netherlands, ISBN 90-386-1700-3, Section 4, 1999. [5] E. Matricciani, “Rate of change of signal attenuation from Sirio at 11.6 GHz”, Electron. lett., Vol.17, No.3 ,1981 [6] P.D. Baxter, G.J.G. Upton and D. Eden, “Revise method for calculation of rain-fade-slope”, Electron. Lett. Vol. 37, No. 10, 2001 [7] Stutzman W.L., Pratt. T., Safaai-Jazi A., Laster J., Nelson B. and Ajaz H. “Results from the Virginia tech propagation experiment using the Olympus satellite 12, 20 and 30 GHz beacons”, IEEE Transac. On Antenna and Propagation, Vol.43, No.1, January 1995. [8] ITU-R P.1623, “Prediction method of fade dynamics on Earth-space paths”, ITU, Geneva, Switzerland, 2003. IIUM Engineering Journal, Vol. 12, No. 5, 2011: Special Issue -1 on Science and Ethics in Engineering Dao et al. 59 [9] F.F. Franklin, K. Fujisaki, and M. Tateiba, “Fade Dynamics on Earth-Space Paths at Ku- Band in Fukuoka, Japan Fade-Slope Evaluation, Comparison, and Model”, IEEE Antenna and Wireless Propagation Letters, Vol. 5, 2006. [10] Adhikari A. and Maitara A., “Studies on the inter-relation of Ku-band scintillations and rain attenuation over an Earth–space path on the basis of their static and dynamic spectral analysis”, Journal of Atmospheric and Solar-Terrestrial Physics (JASTP), Volume 73, Issue 4, November 2010. [11] H. Dao, MD. R. Islam, K. Al-Khateeb, “Fade Dynamics review of microwave signals on Earth-Space paths at Ku-Band”, International Conference on Computer and Communication Engineering Proceeding, May 2008. [12] Y. Maekawa, T. Fujiwara, Y. Shibagaki, T. Sato, M. Yamamoto, H. Hashiguchi, and S. Fukao “Effects of Tropical Rainfall to the Ku-Band Satellite Communications Links at the Equatorial Atmosphere Radar Observatory”, Journal of the Meteorological Society of Japan, Volume 84A, 2006. [13] MR Islam, H. Dao and AH Zayuod, “Rain Fade Slope Analysis for Terrestrial Microwave Link in Malaysia”, 9 th Malaysia International Conference on Communications (MICC) Proceeding, December 2009. [14] Y. Karasawa , T. Matsudo, “Characteristics of Fading on Low-Elevation Angle Earth-Space Paths with Concurrent Rain Attenuation and Scintillation”, Vol.39, No.5, May 1991. [15] I.E. Otung, M.O. Al-Nuaimi, and B.G. Evans, “Extracting Scintillations from satellite beacon propagation data”, IEEE Transac. On Antenna & Propagation, Vol.46, No.10, October 1998. [16] H. Dao, MR. Islam, K. Al-Khateeb, “Preliminary Analysis of Ku-Band Rain Fade Data for Earth-to-Satellite Path Measured in Malaysia”, 10 th Malaysia International Conference on Communications (MICC) Proceeding, October 2011