Microsoft Word - vol11_no1_final IIUM Engineering Journal, Vol. 11, No. 1, 2010 Islam et al. 71 PREDICTION OF SIGNAL ATTENUATION DUE TO DUSTSTORMS USING MIE SCATTERING MD. RAFIQUL ISLAM, ZAIN ELABDIN, OMER ELSHAIKH, OTHMAN O. KHALIFA, AHM ZAHIRUL ALAM, SHEROZ KHAN AND A. W. NAJI Wireless Communication and Signal Processing Research Group, Kulliyyah of Engineering, International Islamic University Malaysia P.O. Box 10, Kuala Lumpur, 50728, Malaysia. E-mail: rafiq@iiu.edu.my ABSTRACT: The present trend in radio design calls for the use of frequencies above 40 GHz for short links carrying wide-band digital communication signals. In order to utilize the new frequency band efficiently, signal attenuation studies due to duststorms is needed urgently for desert areas. This paper presents a mathematical model which has been developed to predict the signal attenuation due to duststorm. The proposed model enables the convenient calculation of the signal path attenuation based on Mie solution of Maxwell's equations for the scattering of electromagnetic wave by dust particles. The predicted values from the proposed mathematical model are compared with the measured values observed in Saudi Arabia and Sudan and show relatively close agreement. Keywords: Signal Attenuation, Duststorm, Mathematical Model, Duststorm Visibility. 1. ABSORPTION AND SCATTERING BY THE ATMOSPHERE When microwaves pass through the medium containing precipitations like rain, snow or dust particles, the signals get attenuated through two phenomena [1]. a) Absorption of energy by these particles. b) Scattering of energy out of the beam by these particles. Microwaves suffer absorption and scattering by the atmosphere especially at higher frequencies where the scattering effects become more severe [2]. The knowledge of these characteristics is essential to design reliable communications. 1.1 Theory of single-particle scattering The basic theory underling mathematical model for attenuation is the theory for single- particle scattering [3]. The characteristics of these two phenomena, scattering and absorption, can be expressed most conveniently by assuming an incident plane wave. The total electromagnetic field around the scattering object is split into an incident and a scattered field as shown in Fig. 1 [4]: IIUM Engineering Journal, Vol. 11, No. 1, 2010 Islam et al. 72 Fig. 1: Scattering configuration. The total amplitude of the fields at any point on the surface of the sphere can be given as [5]: i s H H H  (1) i sE E E  (2) where, Ei and Hi denote the incident electric and magnet field, and Es and Hs denote corresponding scattered fields. 1.2 The total cross-section efficiency factors The power density of the incident wave is uniform in space. Therefore, the power loss of the incident wave can be calculated by multiplying the uniform power density with a cross-section, having the dimension of the area. The total power removed from the incident field by absorption and scattering can be represented by the total or extinction cross-section [4]. When the problem is concerned with the flow of energy, poynting vector must be used. So the power density is represented by the averaged poynting vector S . This poynting vector can be splitted into separate terms using the separation of the total field into the incident and scattered field (Eq. (1) and Eq. (2)), resulting in [4]: i s dS S S S   (3) where, 1 Re( *) 2i i i S E H  (4) 1 Re( *) 2s s s S E H  (5) 1 Re( * *) 2d i s s i S E H E H    (6) Ei (r) E s (r) ε (r) u i r u s IIUM Engineering Journal, Vol. 11, No. 1, 2010 Islam et al. 73 and where, Ei and Hi denote the incident electric and magnet field, whereas Es and Hs denote corresponding scattered fields. The net power flow through the surface of this sphere can be calculated by integrating the radial component of S over the sphere. If the scatterer is a lossless dielectric, the power flow will be zero. If it is lossy (permittivity is complex), the net power flow will represent the absorption rate Pa of energy by the scatterer as shown in Fig. 2 and can be written as: a i s dP P P P    (7) where Pa represents the integral of the radial component of S . u n u n r P s P a P i u i Fig. 2: Scattering points. The medium surrounding the particle is assumed to be lossless, so the net flow of the energy of the incident field is zero for every closed surface (Pi = 0), and thus: a s dP P P   (8) and from Eq. 6; 1 Re ( * *) 2a s d i s s i n s P P P E H E H u dS         (9) where un is the unit vector normal to the surface S and directed outwards and the dot represent the inner product as shown in Fig. 2. Using Maxwell equation the sum of absorbed and scattered energy can be expressed as [4]:   1 22 Im ( , )a s i e s iP P E u f u u k            (10) where, k : wave number ( (2 )k      ). IIUM Engineering Journal, Vol. 11, No. 1, 2010 Islam et al. 74 1 2        : characteristic impedance of the medium. ue : unit vector defining the polarization of the wave. f (us, ui): the scattering amplitude function; it is a function of the direction of propagation ui of the incident field and of the direction us from the object to the observation point. The ratio of the rate of removal of energy (Pa + Ps) to the rate of the incident energy (Pi) on a unit cross-sectional area of the scatterer is defined as the extinction or total cross- section efficiency factors [4]:  4 Im ( , )a st a s e s i i P P u f u u P k           (11) 2. ANALYTICAL MODELS FOR SCATTERING Two models have been presented for an analytical solution of Maxwell's equations for the scattering of electromagnetic wave by dielectric spherical particles, Rayleigh approximation and Mie solution. 2.1 Rayleigh Approximation The assumption for Rayleigh approximation is that (ka < 1), meaning that the radius a of the particle is small compared with the wavelength. For a dielectric sphere that is small compared with wavelength, the field inside the particle is chosen to be the solution of the equivalent dielectric electrostatic problem. In the Rayleigh approximation, the field inside the scattering particle is not modeled properly. For an exact solution the scattered field reflects all the properties of the object which is no longer true for an approximation solution [5]. The extinction or total cross-section efficiency factors by Rayleigh approximation expressed as [6]:     2 42 2 2 2 8 '' 1 12 3 ( ' 1) ( '') 2 t s a a ka a ka                      (12) where, a : the particle radius k : the wavenumber (k =2π / λ) ε' , ε'' : the real and imaginary contributions of the relative dielectric constant of the particles IIUM Engineering Journal, Vol. 11, No. 1, 2010 Islam et al. 75 2.2 Mie Solution Mie solution is a complete analytical solution of Maxwell's equations for the scattering of electromagnetic wave by dielectric spherical particles. In contrast to Rayleigh scattering Mie solutions to scattering embraces all possible ratios of diameter to wavelength [7]. Mie solutions unlike Rayleigh scattering do not depend upon any such limitation and can be utilized to predict attenuation in microwave wave band with high reliability especially at higher frequencies used by new systems with larger bandwidth. The extinction or total cross-section efficiency factors by Mie solutions can be expressed as [6]:        2 3 2 3 1 2 3 2 t ka c c ka c ka       (13) where, a : the particle radius k : the wavenumber (k =2π / λ) λ : the wavelength, and C1, C2 and C3: constants whose values depend on real (ε') and imaginary part (ε'') of the dielectric constant of the particles as:   1 2 2 6 2 C         (14)     2 2 2 22 22 2 6 7 7 4 20 1 5 5 152 3 2 3 4 C                                    (15)          2 4 3 22 2 1 2 2 1 2 94 3 2 C                                  (16) 3. SIGNAL ATTENUATION DUE TO DUSTSTORM The theory of attenuation due to sand and duststorm can be explained in terms of scattering and absorption cross section of a single particle. The methods to predict the signal attenuation due to rain effects can be applied for duststorm because the general model for scattering in sand and dust particle populations is essentially the same as that for a population of hydrometeors; both of them are discrete random medium [3]. IIUM Engineering Journal, Vol. 11, No. 1, 2010 Islam et al. 76 The signal attenuation due to duststorm is estimated generally by solving the forward scattering amplitude function of a single particle [8,9,10]. The solution may be carried out using the Rayleigh approximation or Mie solutions. The method depends largely on the wave number and particle radius. The attenuation of electromagnetic radiation (AT) over a path of extent L through precipitating particles may be written as [11]: 0 L T pA A dx  [dB] (17) where Ap(dB/km) is the specific attenuation characterizing the precipitating particles. The following expression has been used to calculate the attenuation due to rain [11,12,13]:    max min 34.343 10 a p ta A a N a da   [dB/km]. (18) where, N(a)da : the particles number per unit volume of air with particles radius between a and a+da, σt : the total attenuation cross section efficiency factors of particle of radius a. 3.1 Duststorm Predicted Model Starting from the above, we can expressed the attenuation of electromagnetic radiation (AT) over a path of extent L through duststorm as: 0 L T dA A dx  [dB] (19) where Ad is the s pecific attenuation characterizing the duststorm which can be expressed as:    max min 34.343 10 a d ta A a N a da   [dB/km]. (20) and where, N(a)da : the dust particles number per unit volume of air with dust particles radius between a and a+da, σt : the total attenuation cross section efficiency factors of dust particle of radius a. IIUM Engineering Journal, Vol. 11, No. 1, 2010 Islam et al. 77 a. Using Rayleigh Approximation Goldhirsh [14], in using the Rayleigh approximation expression for the total cross- section efficiency factor (σt), expressed the specific attenuation due to duststorm Ad as:   3 2 2 2.317 10 1 2 A V                  [dB/km] (21) where, V : the visibility in kilometer, λ : the wavelength in meter, γ : constants that depend on the type of land from which the storm originated as well as the climatic conditions, ε' , ε'' : the real and imaginary part of the dielectric constant of the dust particles. Goldhirsh formulation is applicable at wavelengths for which the Rayleigh condition is applicable which is equal to 48 GHz as a maximum frequency [14]. b. Using Mie solution In order to derive a mathematical model that can predict the signal attenuation due to duststorm at higher frequencies, the Mie expression for total cross-section efficiency factors (σt) that appear in Eq. (13) was selected for substitution in Eq. (20). So the duststorm attenuation (Ad) can be written as:         max min 2 3 2 33 1 2 34.343 10 2 a d a A ka c c ka c ka N a da        [dB/km] (22) i. Dependence on Visibility To calculate the attenuation by Eq. (22) requires data for the number of particles of dust N, which is difficult to measure accurately. On the other hand, statistical information on duststorm visibility is available. Goldhirsh [14] expresses the visibility in term of the particle density and the radius as: 4 2 5.5 10 ( ) e V km Na   (23) where, N : the dust particles number per unit volume, ae : the equivalent particle radius in meters. IIUM Engineering Journal, Vol. 11, No. 1, 2010 Islam et al. 78 By solving N in the above formula we can express the particle density in term of the visibility and the radius as: 4 2 5 .5 1 0 e N V a   (24) ii. General Formula By substituting the particle density expression which appear in Eq. (24) into Eq. (22), the specific attenuation characterizing the duststorm (Ad ) may alternately be expressed as:       max min 2 4 3 2 33 1 2 3 2 5.5 10 4.343 10 2 a d a e A ka c c ka c ka da V a              [dB/km]. (23) A further approximation can be made in these calculations, assuming that every dust particle in a real storm may be replaced by an equivalent particle (ae) whose radius is the mean radius for all dust particle. By this assumption the value of equivalent particle radius (ae) considered as constant value and Eq. (23) may alternately be expressed as an algebraic expression:        2 4 3 2 33 1 2 3 2 5.5 10 4.343 10 2 d e e e e A ka c c ka c ka V a              [dB/km]. (24) and can be simplified to:   2 3 5 3 6 4 1 2 3 0.38 d e e eA c k a c k a c k a V     [dB/km]. (25) By substituting k = 2π / λ in Eq. (25): 3 5 6 3 4 1 2 33 4 0.38 8 32 64 d e e eA c a c a c a V               [dB/km]. (26) which can be simplified to: 3 4 1 2 33 4 94.3 3721.2 23381e e ed a a a A c c c V V V         [dB/km]. (27) Finally, we can express the specific attenuation due to duststorm Ad (dB/km) as: 3 4 3 4 e e e d X a Y a Za A V V V      [dB/km]. (28) where, ae : the equivalent particle radius in meters, IIUM Engineering Journal, Vol. 11, No. 1, 2010 Islam et al. 79 V : the visibility in kilometer, λ : the wavelength in meter and; X, Y and Z: constants whose values depend on real (ε') and imaginary part (ε'') of the dielectric constant of the dust particles as:  2 2 565.8 2 X         (29)     2 2 3 22 22 2 6 7 7 4 20 1 5 3.7 10 5 152 3 2 3 4 Y                                      (30)          2 4 4 22 2 1 2 2 1 2 9 3.12 10 2 Z                                   (31) iii. The Formula in Term of Frequencies (f) in GHz For more simplification of Eq. (31) it is better to express the wavelength (λ) in term of frequencies (f) in GHz which is easier to use by microwave network engineers. So the specific attenuation of duststorm Ad (dB/km) can be expressed as: 3 3 4 4 3 4(0.3) (0.3) (0.3) e e e d X a f Y a f Za f A V V V    [dB/km]. (32) and can be simplified to:  2 2 3 3ed e e a f A x ya f za f V    [dB/km]. (33) where, ae : the equivalent particle radius in meters, V : the visibility in kilometer, f : the frequency in GHz, and x, y and z: constants whose values depend on real (ε') and imaginary part (ε'') of the dielectric constant of the dust particles as: IIUM Engineering Journal, Vol. 11, No. 1, 2010 Islam et al. 80  2 2 1886 2 x         (34)     2 2 3 22 22 2 6 7 7 4 20 1 5 137 10 5 152 3 2 3 4 y                                      (35)          2 4 4 22 2 1 2 2 1 2 9 379 10 2 z                                   (36) 3.2 Determination of Attenuation at Different Bands A further simplification can be made in this general formula (Eq. (33)) to more effective and become easier for use by microwave network designers and engineers. This can be achieved by subdividing the general formula into several formulas, each one can deal with specific microwave band suitable to calculate the signal attenuation due to dust particles at every microwave bands. It is clear from the general formula that the signal attenuation due to dust particles depends on visibility, frequency, dust particles radius and dielectric constant. The dielectric constant values of dust particles are important components in the determination of attenuation. A number of investigators have measured the refractive index of sand and dust samples [9, 15, 16]. A further approximation can make in the new subdivided formulas, assuming that a single value of the refractive index can be fitted for all frequencies range in specific microwave band. i. Attenuation at S-band For S-band which is in the range between 2 – 4 GHz, the dielectric constants measured were ε'= 4.56, ε''= 0.251 [15]. The specific attenuation of duststorm Ad (dB/km) can be expressed as: 3 3 4 4 511 433 2.46 10e e edS a f a f a f A V V V       [dB/km]. (37) Figure 3 showed the signal attenuation (dB/km) versus frequency plot at S-bands for four different values of visibility at equivalent particle radius of 50 μm. IIUM Engineering Journal, Vol. 11, No. 1, 2010 Islam et al. 81 Signal Atte nuation due to Dus tstorm with Visibility (10 -500 m) at 50 μm dus t particle radius in S-band. 0 0.05 0.1 0.15 0.2 0.25 3 3.1 3.2 3.3 3.4 3.5 3.6 3.7 3.8 3.9 4 S-band Frequency (2 - 4 GHz). S ig n a l A tt en u a ti o n (d B /k m ). V=500m V=100m V=50m V=10m Fig. 3: Signal attenuation (dB/km) versus frequency at S-band. ii. Attenuation at X-band For X-band which is in the range between 8 – 12 GHz, the dielectric constants measured were ε' = 5.73, ε''= 0.415 [8]. The specific attenuation of duststorm Ad (dB/km) can be expressed as: 3 3 4 4 513 702 2.5 10e e edX a f a f a f A V V V       [dB/km]. (38) The signal attenuation (dB/km) versus frequency plot at X-bands for four different values of visibility at equivalent particle radius of 50 μm is shown in Fig. 4. Signal Atte nuation due to Dus tstorm with Vis ibility (10 -500 m) at 50 μm dust particle radius in X-band. 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 8 8.5 9 9.5 10 10.5 11 11.5 12 X-band Frequency (8 - 12 GHz). S ig n a l A tt en u a ti o n ( d B /k m ). V=500m V=100m V=50m V=10m Fig. 4: Signal attenuation (dB/km) against frequency at X-band. IIUM Engineering Journal, Vol. 11, No. 1, 2010 Islam et al. 82 iii. Attenuation at Ku-band For Ku-band which is in the range between 12 – 18 GHz, the dielectric constants measured were ε' = 5.5, ε''= 1.3 [16]. The specific attenuation of duststorm Ad (dB/km) can be expressed as: 3 3 4 4 542 2806 2.4 10e e edku a f a f a f A V V V       [dB/km]. (39) Figure 5 showed the signal attenuation (dB/km) versus frequency plot at Ku-bands for four different values of visibility at equivalent particle radius of 50 μm. Signal Atte nuation due to Dus ts torm with Vis ibility (10-500 m) at 50 μm dus t particle radius in Ku-band. 0 0.5 1 1.5 2 2.5 3 3.5 4 12 12.5 13 13.5 14 14.5 15 15.5 16 16.5 17 17.5 18 Ku-band Frequency (12 - 18 GHz). S ig n a l A tt en u a ti o n ( d B /k m ). V=500m V=100m V=50m V=10m Fig. 5: Signal attenuation (dB/km) vs frequency at Ku-band. iv. Attenuation at K-band For K-band which in the range between 18 – 26.5 GHz, the dielectric constants measured were ε' = 5.1, ε''= 1.4 [16]. The specific attenuation of duststorm Ad (dB/km) can be expressed as: 3 3 4 4 550 3368 2.38 10e e edk a f a f a f A V V V       [dB/km]. (40) The signal attenuation (dB/km) versus frequency plot at K-bands for four different values of visibility at equivalent particle radius of 50 μm is shown in Fig. 6. IIUM Engineering Journal, Vol. 11, No. 1, 2010 Islam et al. 83 Signal Atte nuation due to Dus tstorm with Visibility (10 -500 m) at 50 μm dust particle radius in K-band. 0 1 2 3 4 5 6 7 8 18 18.5 19 19.5 20 20.5 21 21.5 22 22.5 23 23.5 24 24.5 25 25.5 26 26.5 K-band Frequency (18 - 26.5 GHz). S ig n a l A tt en u a ti o n ( d B /k m ) V=500m V=100m V=50m V=10m Fig. 6: Signal attenuation (dB/km) vs frequency at K-band. v. Attenuation at Ka-band For Ka-band which in the range between 26.5 – 40 GHz, the dielectric constants measured were ε' = 4, ε''= 1.325 [16]. The specific attenuation of duststorm Ad (dB/km) expression can be written as: 3 3 4 4 566 3712 2.2 10e e edka a f a f a f A V V V       [dB/km]. (41) Figure 7 showed the signal attenuation (dB/km) versus frequency plot at Ka-bands for four different values of visibility at equivalent particle radius of 50 μm. Signal Atte nuation due to Dus ts torm with Visibility (10-500 m) at 50 μm dus t particle radius in Ka-band. 0 2 4 6 8 10 12 14 26.5 28 29.5 31 32.5 34 35.5 37 38.5 40 Ka-band Frequency (26.5 -40 GHz). S ig n a l A tt en u a ti o n ( d B /k m ). V=500m V=100m V=50m V=10m Fig. 7: Signal attenuation (dB/km) against frequency at Ka-band. IIUM Engineering Journal, Vol. 11, No. 1, 2010 Islam et al. 84 vi. Attenuation at W-band For W-band which in the range between 56 – 100 GHz, the dielectric constants measured were ε' = 3.5, ε''= 1.64 [16]. The specific attenuation of duststorm Ad (dB/km) expression is; 3 3 4 4 594 6232 2.1 10e e edW a f a f a f A V V V       [dB/km]. (42) The signal attenuation (dB/km) versus frequency plot at W-bands for four different values of visibility at equivalent particle radius of 50 μm is shown in Fig. 8. Signal Atte nuation due to Dus ts torm with Vis ibility (10 -500 m) at 50 μm dus t particle radius in W-band. 0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 W -ban d Fre qu e n cy (56 - 100 GHz ). S ig n a l A tt en u a ti o n ( d B /k m ). V=500m V=100m V=50m V=10m Fig. 8: Signal attenuation (dB/km) versus frequency at W-band. 4. VALIDITY OF THE MATHEMATICAL MODEL In order to verify the validity of the proposed mathematical model for prediction of signal attenuation due to duststorm, a comparison run between the results obtained from the mathematical model and field measurements were made. At the same time a comparison run between the proposed model and one of the current models was also investigated. 4.1 Saudi Arabia Measurements During the course of measurements, Alhaider and Ali recorded several duststorms during 1987 in the city of Riyadh, Saudi Arabia [17]. The duststorm was observed both by the meteorological sensors and by the radio links. Visibility reduction was measured by using the known marked distance method. For most of the events, the cell size of the storm exceeded 20 km. The measurements run over 14 km microwave link at 40 GHz (Ka-band). The millimeter wave transmitters and receivers are placed at 100 m and 25 m above ground, respectively. IIUM Engineering Journal, Vol. 11, No. 1, 2010 Islam et al. 85 4.2 Sudan Measurements Measurement in Khartoum-Sudan on September 1, 2007 shows attenuation 0.67 dB per kilometer observed by the author on 15 km link at 13 GHz (Ku band) [18]. The duststorm produced a visibility smaller than 5 m. The predicted value found by proposed model using the Ku-band specific attenuation formula is 0.55 dB per kilometer for the same visibility which is quite close to measured value of 0.67 dB/km. 4.3 Comparison of Measured and Calculated Values Calculated attenuation values for different values of visibility at 40 GHz are given in Table 1 together with attenuation measured values observed in Saudi Arabia. The predicted values calculated by the proposed Ka-band specific attenuation formula (Eq. 41) and Goldhirsh model (Eq. 21). The following values are considered in calculation:  Visibility from less than 1 km to more than 5 km.  Equivalent dust particle radius (ae) = 30 μm.  The dielectric constants (ε'=4, ε''=1.325).  Frequency = 40 GHz. Table 1: Result of Measurements and Calculations of Signal Attenuation due to Duststorm for Different Value of Visibility at 40 GHz. Calculated Value (dB/km) Measured Value (dB/km) A c (Goldhirsh Model) A c (Proposed Model) A m Visibility (km) 0.02 0.13 0.14 0.625 0.01 0.064 0.1 1.25 0.007 0.06 0.071 1.42 0.003 0.021 0.05 3.75 0.002 0.014 0.036 5.56 The signal attenuation (dB/km) versus different values of visibility plot at frequency equal to 40 GHz for measured and calculated values is shown in Fig. 9. Unlike the predicted values from Goldhirsh model the attenuation values calculated by the proposed model show close agreement to measurements. IIUM Engineering Journal, Vol. 11, No. 1, 2010 Islam et al. 86 Comparison of Measured and Calculated Values of Attenuation at 40GHz. 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 5.56 3.75 1.42 1.25 0.625 Visibility (km) A tt e n u a ti o n ( d B /k m ) Measured Value Goldhirsh Model(Ray) Proposed Model (Mie) Fig. 9: Comparison between Measured and Calculated Attenuation. 5. CONCLUSION A model has been developed to predict microwave attenuation due to dust particles using Mie solution of Maxwell's equations for the scattering of electromagnetic wave. In this proposed model the term visibility (V) is applied to denote the degree of duststorm density instead of total number of dust particles (N). The predictions show that the attenuation varies from 13 dB/km to 0.2 dB/km at 40 GHz for dust particle radius equal to 50 μm and the visibility varies from 10 m to 500 m. At higher frequency of 100 GHz, the attenuation varies from 47 dB/km to 2dB/km for dust particle radius equal to 50 μm and the visibility varies from 10 m to 500 m. The predicted values from the mathematical model are compared with the measured values observed in Saudi Arabia and Sudan. The proposed model predicts the measured data more accurately than existing model. The accurate prediction is useful to design higher frequency links at areas affected by duststorms. REFERENCES [1] R. K. Crane, Propagation Phenomena Affecting Satellite Communication Systems Operating in the Centimeter and Millimeter Wavelength Bands, Proc. IEEE, 59(2), 173-88, 1971. [2] K. Williams and R. Greeley, Radar Attenuation by Sand: Laboratory Measurements of Radar Transmission, IEEE Transactions on Geoscience and Remote Sensing, vol. 39, No. 11, Nov 2001. [3] G. Brussaard and P.A. Watson, Atmospheric Modelling and Millimeter Wave Propagation (Ghapman & Hall, UK 1995). IIUM Engineering Journal, Vol. 11, No. 1, 2010 Islam et al. 87 [4] Ishimaru, Wave Propagation and Scattering in random media (IEEE, N.Y. 1997). [5] B.R. Vishvakarma and C.S. Rai, Limitations of Rayleigh Scattering in The Prediction of Millimeter Wave Attenuation in Sand and Dust Storms, Geoscience and Remote Sensing Symposium, IEEE Intr., 1993. [6] Collin R., E., Antenna & radiowave propagation (McGraw-hill International Edition, Singapore, 1985). [7] V. E. Cachorro, New Improvements For Mie Scattering Calculations, eprintarXiv:physics/0103052, March 2001. [8] S. I. Ghobrial and S. M. Sharief, Microwave attenuation and cross polarization in dust storms, IEEE Trans. Antennas Propagat., v. AP-35, pp. 418–425, Apr. 1987. [9] T. S. Chu, "Effects of sandstorms on microwave propagation." Bell Sys. Tech. J.. v. 58, no. 2, Feb. 1979. [10] J. Goldhirsh, A parametric review and assessment of attenuation and backscatter properties associated with dust storms over desert regions in the frequency range of 1 to 10 GHz, IEEE Trans. Antennas Propagation, vol. AP-30, no. 6, pp. 1121-1127, 1982. [11] R. R. Rogers, Statistical Rainstorm Models: Their Theoretical and Physical Foundations, IEEE Trans. On Antennas and Propagation, July 1976. [12] R. L. Olsen et al., The aRb relation in the calculation of rain attenuation, IEEE Trans. On Antennas and Propagation, vol. AP-26, no. 2, Mar. 1978. [13] Oguchi, Electromagnetic Wave propagation and Scattering in hydrometeors, Proceedings of The IEEE, vol.71,No.9, Sep 1983. [14] J. Goldhirsh, “Attenuation and Backscatter from a Derived Two-Dimensional Duststorm Model” IEEE Trans. Antennas Propagation, v. 49, no. 12, pp. 1703-1711, 2001. [15] S. I. Ghobrial, The effect of sand storms on microwave propagation, Proc. Nat. Telecommun. Conf., Houston, TX. vol. 2, Proc. No. CH1539-6/80/0000-0216, pp. 43.5.1-43.5.4, 1980. [16] Y. Ruike, W. Zhensen, and Y. Jinguang, The Study Of MMW and MW Attenuation Considering Multiple Scattering Effect in Sand and Dust Storms at Slant Paths, International Journal of Infrared and Millimeter Waves, Vol. 24, No. 8, Aug 2003. [17] Alhaider, M.A. Ali, A.A., Experimental studies on Millimeterwave and Infrared propagation in Arid Land: The Effect of Sand Storms, Sixth International Conference on Antennas and Propagation ICAP, 1989. [18] Zain Elabdin Omer Elshaikh, Md. Rafiqul Islam, Othman O. Khalifa and Momoh Jimoh E Salami, Duststorm Measurements for the Prediction of Attenuation on Microwave Signals in Sudan, Submitted to International Conference on Computer and Communication Engineering 2008 (ICCCE‘08) Kuala Lumpur, May 2008.