Microsoft Word - ELEC991201.doc Science and Technology, 5 (2000) 47-54 © 2000 Sultan Qaboos University 47 Electromagnetic Modeling of the Propagation Characteristics of Satellite Communications Through Composite Precipitation Layers, Part1: Mathematical Formulation * H.M.Al-Rizzo,** H.T. Al-Hafid, and *** J.M.Tranquilla *College of Information Science and Systems Engineering, University of Arkansas, 2801 South, University Avenue, Little Rock, AR 72204; and Department of Electrical and Electronics Engineering, College of Engineering, Sultan Qaboos University, P.O. Box 33, Al Khod 123, Muscat, Sultanate of Oman;**Electrical Engineering Department, United Arab Emirates University, Al Ain, United Arab Emirates; ***EMR Microwave Technology Corporation, 64 Alison Blvd., Suite 105, Fredericton, NB,Canada, E3C 1N2. خالل طبقات المطر جة انتشار الموجات الكهرومغناطيسية لمنظومات األقمار الصناعيةنمذ النموذج الرياضي : المتراكمة حسين مظفر الريزو ، حافظ طه الحافظ ، و جم ترانكوال ى من في هذا البحث تم تطوير نموذج رياضي لبرنامج تنبؤ انتشار الموجات الكهرومغناطيسية في الطبقة السفل :خالصة هذا وقد أخذت الدراسة في االعتبار التوزيع العشوائي ألبعاد .الغالف الجوي وذلك في إطار الموجات الدقيقة والمللمترية .أشكال وكثافة الدقائق المعلقة في الجو والتي تشمل الثلج وقطرات المطر ABSTRACT: A systematic and general formulation of a Propagation Simulation Program (PSP) is developed for the coherent field of microwave and millimeter wave carrier signals traversing intermediate layered precipitation media, taking into account the random behavior of particle size, orientation, shape and concentration distributions. Based on a rigorous solution of the volumetric multiple-scattering integral equations, the formalism offers the capability of treating the potential transmission impairments on satellite- earth links and radar remote sensing generated by composite atmospheric layers of precipitation in conjunction with the finite polarization isolation of dual-polarized transmitting and receiving antennas. A multi-layered formulation is employed which encompasses an ensemble of discrete particles comprising an arbitrary mixture of ice crystals, melting snow and raindrops that may exist simultaneously along satellite-earth communication paths. KEYWORDS: Microwave, Multiple-Scattering, Remote Sensing. ecent demand for additional communication capacity, combined with the requirement of alleviating the congested electromagnetic (EM) spectrum, has led system planners to utilize increasingly higher carrier frequencies. Many attractive features for earth-space communications, radar, remote sensing and imaging applications such as improved resolution and wide absolute bandwidths essential for civilian and military high data rate transmission needs are feasible in the millimeter wave (MMW) regime. In this paper the analytical formulation of a PSP, based on a rigorous solution of the volumetric multiple-scattering integral equations (Ishimaru, 1977), is presented to evaluate the impairments caused by precipitation along arbitrary slant earth-space paths AL-RIZZO, AL-HAFID AND TRANQUILLA 48 utilizing the Generalized Multipole Technique (GMT) developed by Al-Rizzo and Tranquilla (1997). A concise matrix approach is employed for the treatment of the ensemble-averaged, coherent field of MMW signals traversing a random collection of hydrometeors including ice crystals, melting snow and rainfall. The purpose of the work reported in this paper is twofold. First, to integrate the simulation through the whole spectrum of precipitation using a multi-layered formulation which includes inhomogeneous mixture of ice at top of rain clouds and the transition from melting snow to rain below the melting layer. The second objective, and perhaps the foremost, is to extend the frequency range and improve the accuracy of previously published results, within the limitations of available meteorological data and physical parameters concerning the various forms of hydrometeors in the lower part of the atmosphere. General Medium Modelling An arbitrarily polarized time-harmonic incident plane EM wave, propagating in the positive z- direction, is assumed. The polarization state is expressed in terms of the angle pair, (γ,δ), where γ = tan-1(E2/E1), E1 and E2 are the maximum projections of the instantaneous electric field vector along the x and y axes, respectively, and δ is the time-phase angle by which Ey leads Ex (Kraus and Carver, 1973). A Cartesian coordinate system is employed, the z-axis defines the incident wave direction, pointing from the satellite antenna toward the receiving antenna at the earth-terminal location. The x and y axes coincide with the reference horizontal and vertical directions, respectively, forming a right-handed triplet where the x-axis is parallel to the local horizontal at the receiving end. Complex normalized column vectors; Vt of impressed and Vr of received phasor voltages are defined such that Pti = Vti . V * ti, is the input power to the transmitting antenna and Pri = Vri . V * ri, i = 1, 2 is the power delivered to a matched load by the receiving antenna. The finite polarization isolation, inherently present in the transmitting and receiving antennas, is accounted for via two matrices; [ ]tA and [ ]rA which relate the incident Ei and received Er electric field vectors to Vt and Vr such that [ ] tt VAE i = and [ ] rrr EAV = where (Bostian et al 1982)             =      2t 1t 2tj 2t 1tj 1t 2t1t y i x i V V esinesin coscos E E δδ γγ γγ             =      − − y r x r 2rj 2r2r 1rj 1r1r 2r 1r E E esincos esincos V V δ δ γγ γγ (1) The subscript t (r) refers to the transmitting (receiving) antenna, respectively and the two (nominally) orthogonal polarizations are denoted by the additional subscript 1 (2). The coherent electric field vector exiting a discrete, random anisotropic and lossy dielectric medium is related to the corresponding incident field using the eigenvalues and eigenvectors of the ensemble- averaged, complex propagation constant matrix [], the exponential form of which can be expressed as (Tsoloakis and Stutzman, 1982) [ ] [ ]             =      − − − − y i x i1 2j 1j L0jk y r x r E E e0 0e e E E QQ λ λ (2) where L is the path length along the propagation direction and ko is the free-space wave number. The term e-jkoL characterizes the pure phase-shift properties of the scatterer-free medium and [Q] is the matrix of non-trivial eigenvector solutions of [] which are found from MODELING THE PROPAGATION OF SATELLITE COMMUNICATIONS 49 ( ) ( ) 2 k2kkkk 212 2 22112211 2,1 〉〈+〉〈−〉〈±〉〈+〉〈 =λ (3) [ ]       〉〈−〉〈− −〉〈−〉〈 = 1212 222122 kk kk λλ Q (4) When the matrix of eigenvectors and its inverse are substituted into (2), then (2) can be transformed into the following concise form Er(L) = e-jkoL [] Ei(L) (5) The elements of the transmission matrix [], which describes the amplitude, phase and polarization of the coherent field exiting the medium, can be written as ( ) ( ) 21 1j 122 2j 222 11 e ke k T λλ λλ λλ − −〉〈−−〉〈 =〉〈 −− 21 1j2j 122112 e e kTT λλ λλ − − 〉−〈=〉〈=〉〈 −− ( ) ( ) 21 2j 122 1j 222 22 e ke k T λλ λλ λλ − −〉〈−−〉〈 =〉〈 −− (6) The transfer-matrix equation describing the input-output relationship of the complete orthogonal- polarized, frequency re-use system can then be found by combining (1) and (5) as Vr = e - jkoL [Ar] [] [At] Vt (7) Theoretical Development For The Rainfall Medium A PSP must include empirical or statistical averaging to account for the temporal and spatial randomness of the rain characteristics. Falling raindrops are nonspherical in shape and are influenced by aerodynamic forces. As the size of a raindrop increases, a concave depression develops. Pruppacher and Pitter (1971),(P-P), provided 13 realistic raindrop shapes with spherical equivolume radii from 0.25 to 3.25 mm. With the exception of the limited cases presented by Oguchi (1977), Fang and Lee (1978) and Li et al. (1995), previous models considered only spherical and/or oblate spheroidal raindrops. The GMT is utilized in this paper to compute the forward scattering functions of an oblate and/or P-P raindrop, the symmetry axis of which coincides with the z-axis as shown in Figure 1. A plane EM wave, incident in the ki direction and inclined at an angle β with respect to the horizontal plane, is assumed. The incident electric field vector is given by .rk iE i i je−= (8) where the unit vector i describes the polarization state of the incident field. The propagation direction and the raindrop's axis of rotational symmetry define the incident plane. The incident electric field vector is decomposed into two mutually orthogonal components, parallel C iE and perpendicular ⊥ iE to the plane of incidence, respectively, as ⊥⊥+= iiiii iEiEE CC (9) AL-RIZZO, AL-HAFID AND TRANQUILLA 50 In the far-field region, the forward-scattered wave has the same polarization as the incident field and the scattering process can be modeled by the following diagonal matrix equation             =      ⊥⊥⊥ −⊥ C CC C i i w0, w0, rjk s s E E ),,a(f0 0),,a(f r e E E 0 εβ εβ (10) where f║,║(a0,β,εw) and f┴,┴(a0,β,εw) are the forward scattering amplitude functions for the two principal incident polarizations and εw is the complex relative dielectric constant of water. Figure 1. Geometry for the scattering from an axisymmetric raindrop. At a given elevation angle, β and drop temperature, f║,║(a0,β,εw) and f┴,┴(a0,β,εw) are computed for 13 drops with spherical equivolume radii from 0.25 mm to 3.5 mm within a convergence criterion of 10-4 for a0≤ 2 mm and 10-3 for a0≥ 2 mm for the two principal polarizations. The integrals of the complex-valued scattering functions over the appropriate raindrop size distribution have been numerically evaluated using a 14-point composite Simpson’s rule. The orientation geometry of an oblate and/or P-P falling raindrop is defined in the x-y plane in terms of the apparent canting angle between the projection of the symmetry axis and the reference polarization directions (Oguchi, 1977). For an assemblage of identically oriented raindrops, a diagonal transmission matrix is defined, the diagonal elements describe the coherent attenuation and phase shift encountered by the principal polarizations. If the canting angle is denoted by θ, a matrix transformation is performed to bring the reference polarizations into the principal axes of the raindrop via the rotation matrix MODELING THE PROPAGATION OF SATELLITE COMMUNICATIONS 51       − = θθ θθ θ cossin sincos R (11) The components of the forward-scattered electric field vector are next rotated back into the reference x and y directions by applying the inverse rotation 1R −θ . The ensemble-averaged, complex propagation constant matrix, [ ]〉〈 1k which characterizes the coherent forward scattering effects due to a population of raindrops with a prescribed size, canting angle and shape distributions can then be written as [ ] [ ] ds dda )z,s,,a(N k 2 00 0 1 θθ π ∫∫∫=〉〈 fk (12) where N (a0,θ ,s, z) is the number of raindrops of a given shape, s, within a specified size, a0, and orientation, θ, intervals per unit volume of space which is a function of the spatial variable, z along the propagation direction. The forward-scattering matrix obtained from (10) after applying the rotation matrix of (11) and its inverse, when averaged over the statistical parameters of shape and canting angle distributions, can be written as { }〉〈+〉〈−+=〉〈 ⊥⊥ θθ 2,2,s0ss11 sinfcosf)P1()a(fPf CC { } 2 2sin ff)P1(ff ,,s2112 〉〈 −−=〉〈=〉〈 ⊥⊥ θ CC { }〉〈+〉〈−+=〉〈 ⊥⊥ θθ 2,2,s0ss22 cosfsinf)P1()a(fPf CC (13) where fs is the forward scattering amplitude of a spherical water drop. The common assumption of statistically independent shape and canting angle distributions has been adopted in (13). A discrete bimodal shape distribution is assumed, Ps being the effective fraction of spherical and the remaining (1 - Ps) are oblate and/or P-P drops. For linearly polarized geo-stationary satellite transmission, the horizontal and vertical orientations at the satellite are referenced with respect to the subsatellite point as either in or perpendicular to the plane of the equator, respectively. For earth stations not located on the meridian, the electric field vector appears tilted from the local vertical or horizontal by an angle τ where (Allnutt, 1989)       −= − η ξ τ sin tan tan90 10 (14) and ξ and η are the longitude and latitude of the earth station, respectively, defined with respect to the subsatellite point. If the receiving antenna is aligned to the satellite transmission, then equation (7) must be modified to [ ][ ] [ ][ ][ ] tt 1 r L0jk r e VARTRAV ττ 〉〈= −− (15) where the rotation matrix       − = ττ ττ τ cossin sincos R (16) AL-RIZZO, AL-HAFID AND TRANQUILLA 52 transforms the satellite's vertical and horizontal to the local vertical and horizontal at the earth station. Theoretical Development For The Melting – Layer Medium The melting layer is a prominent feature of widespread rain associated with stratiform precipitation. It contains a mixture of gradually melting snowflake aggregates, which eventually form raindrops, attaining the terminal fall speed of water upon complete melting at the bottom of the layer (Ekpenyong and Srivastava, 1970). The establishment of a simulation model for the melting layer is imperative to accurately predict the total attenuation and phase shift for low- elevation satellites where a portion of the radio path may be contained within the layer boundary particularly in temperate maritime, northerly latitude and tropical climates. The meteorological model of the melting process developed by Ekpenyong and Srivastava (1970) is adopted in this paper with a spectrum of polydispersed snowflake sizes, densities and fall speeds. Due to the lack of definitive information on the actual shape and orientation of partially melted snowflakes, the present model assumes a two-layered spherical scatterer shape. Initially, a dry snowflake of radius rs and density ρs originating at the top of the melting layer is considered. A two-layered, concentric spherical snowflake of radius r0 with the melted water forming a uniform shell around a snow core of radius ri is assumed. If y denotes the fraction of total mass melted, then 3 w 3 i 3 0 r rr y − = (17) where rw is the radius of the resulting water drop at the bottom of the layer. The particle number density and the terminal fall speed, denoted respectively by N0 and v, depend on r0 and the height of fall (z - z0) below freezing level where the 0 0C isotherm starts at z = z0, y = 0 at the top of the layer. The precipitation intensity is considered constant throughout the melting region with aggregation and breakup ignored such that a one-to-one correspondence exists between raindrops at the bottom and dry snowflakes at the top. If N(rw)drw denotes the raindrop size distribution, the corresponding distribution of the melting snow particles at any height within the melting zone can be estimated from (Dissanayake and McEwan, 1978) ( ) www 000 ww 0000 dr )r(N)zz,r(v )r(v dr zz,rN − =− (18) Consideration of mass conservation in a certain air volume implies ( ) 3ww3i30w3is3ss rrrrr ρρρρ =−+= (19) If the density and kinematic viscosity of air and mass of a snowflake are assumed constants during the melting process, it can be shown that the product vr remains constant that is ww00ss rvrvrv == (20) From (19) and (20), the terminal fall speed, vs of a dry snowflake can be obtained as 3/1 w s ws vv       = ρ ρ (21) At any stage of melting, the fall speed, v0, can be obtained from (17), (20) and (21) as MODELING THE PROPAGATION OF SATELLITE COMMUNICATIONS 53 y v v 11 v v 3 3 w s s 0               −− = (22) For a dry snowflake, ρs approaches a constant value of 0.008 g/cm3 for r0 > 0.75 cm. For smaller radii, ρs may be taken as 0.3 g/cm3. In between these two limiting values, the density can be obtained from (Ekpenyong and Srivastava, 1970) 2 s s r 005.0 =ρ (23) where rs is in cm, ρs is in g/cm3 and 3/1 s w ws rr       = ρ ρ (24) The dielectric constant of the snow core, εs is evaluated from the following expression (Dissanayake and McEwan, 1978) 2 1 2 1 i i i s s s + − = + − ε ε ρ ρ ε ε (25) where a snowflake is treated as a two-component mixture of ice and air, εi, ρi are the complex dielectric constant and density of ice, respectively. Once the mass fraction of a melted snowflake has been computed, the core and outside radii, terminal fall speed, dielectric constant of the core and particle-size distributions corresponding to a particular size of an eventually melted drop can be evaluated at any prescribed height within the melting layer. Conclusions A theoretical formulation is introduced for the coherent field of a monochromatic EM signal traversing a discrete medium of composite precipitation layers. The model accounts for the random behavior of particle size, orientation, shape and intensity distributions. Although the degrading effects due to a single tropospheric layer of precipitation, which may hinder reliable space communications, were investigated extensively over the past three decades, the contributions at frequencies in the upper MMW band are still scarce. Furthermore, little has been done to integrate the effects throughout the whole spectrum of precipitation in the lower part of the atmosphere. References ALLNUTT, J.E. 1989. Satellite-to-Ground Radio Wave Propagation, London, U.K.: Peter Peregrinus Ltd. AL-RIZZO, H.M. and TRANQUILLA, J.M. 1997. Application of the generalized multipole technique (GMT) to high- frequency electromagnetic scattering from perfectly conducting and dielectric bodies of revolution, J. Comput. Phys., 136: 1-18. BOSTIAN, C.W, STUTZMAN, W.L. and GRAINES, J.M. 1982. A review of depolarization modeling for earth-space radio paths at frequencies above 10 GHz, Radio Sci., 17: 1231-1241. DISSANAYAKE, A.W. and MCEWAN, N.J. 1978. Radar and attenuation properties of rain and bright band, IEE Conf. Publ. 169-2, 125-129. EKPENYONG, B.E. and SRIVASTAVA, R.C. 1970. Radar characteristics of the melting layer- A theoretical study, Tech. Rep. No. 16, Dept. of the Geophys. Sci., The University of Chicago and Illinois Inst. of Techn., Dept. of Elect. Eng., Chicago, Illinois. FANG, D.J. and LEE, F.J. 1978. Tabulation of raindrop induced forward and backward scattering amplitudes, Comsat. Tech. Rev., 8: 455-486. AL-RIZZO, AL-HAFID AND TRANQUILLA 54 ISHIMARU, A. 1977. Theory and application of wave propagation and scattering in random media, Proc. IEEE, 65: 1030- 1061. KRAUS, J.D. and CARVER, K.R. 1973. Electromagnetics, New York: McGraw-Hill. LI LE-WEI, PANG-SHYAN KOOI, MOOK-SENG LEONG, TAT-SOON YEO, and MIN-ZHAN GAO, 1995. Microwave attenuation by realistically distorted raindrops: Part I- Theory, IEEE Trans. Antennas Propagat., AP-43: 811-821. OGUCHI, T. 1977. Scattering properties of Pruppacher-and-Pitter form raindrops and cross polarization due to rain : calculations at 11, 13, 19.3 and 34.8 GHz, Radio Sci., 12: 41-51. PRUPPACHER, R. and PITTER, R.L. 1971. A semi-empirical determination of the shape of cloud and rain drops, J. Atmos. Sci., 28: 86-94. TSOLOAKIS, A. and STUTZMAN, W.L. 1982. Multiple scattering of electromagnetic waves by rain, Radio Sci.,12: 1495- 1502. Received 1 December 1999 Accepted 19 June 2000