Open access peer-reviewed chapter

Atmospheric Attenuation due to Humidity

Written By

Milda Tamošiunaite, Stasys Tamošiunas, Mindaugas Žilinskas and Milda Tamošiuniene

Submitted: November 13th, 2010 Reviewed: March 3rd, 2011 Published: June 21st, 2011

DOI: 10.5772/21430

Chapter metrics overview

7,888 Chapter Downloads

View Full Metrics

1. Introduction

Humidity remains in the atmosphere even on bright days. Water of all three states can be found naturally in the atmosphere: liquid (rain, fog, and clouds), solid (snowflakes, ice crystals), and gas (water vapour). Water in any state is an obstacle in the link of the electromagnetic wave. When the wave passes through the water particles, a part of its energy is absorbed and a part is scattered. Therefore the electromagnetic wave is attenuated. Prediction of the influence of these factors is very important in radio system design. Attenuation due to rain, fog, and clouds can lead to the perturbations of the wireless, mobile, satellite and other communications. Another problem is the refractive index of the atmosphere, which affects the curvature of the electromagnetic wave path and gives some insight into the fading phenomenon. The anomalous electromagnetic wave propagation can cause disturbances to radar work, because variation of the refractive index of the atmosphere can induce loss of radar coverage. Accurate prediction of losses due to these factors can ensure a reliability of the radio system, decrease an equipment cost, furthermore, the radio systems can become less injurious to health of people.

When there are no possibilities to gather data for calculations of the specific attenuation due to rain, clouds and fog, and atmospheric refractive index, the values recommended by the International Communication Union’s Radiocommunication sector (ITU-R) can be used. But the recommended values are not always exact. In design of the radio links, the most desirable operating frequencies are below 10 GHz, because in such cases atmospheric absorption and rainfall loss may generally be neglected (Freeman, 2007). However, in most countries, the frequency-band below 10 GHz is highly congested. In addition, high frequencies provide larger bandwidth, narrower beam width, good resolution and smaller component size (Bhattacharyya et al., 2000). Therefore, the operating frequencies of 10 GHz and above are often used in design of radio systems. The higher the operating frequency, the greater attenuation due to hydrometeors (rain, cloud, fog, snow, and etc.) is observed (Tamošiūnaitė et al., 2010a).

In (Ishimaru, 1978), it was mentioned that the electromagnetic wave attenuation due to snow is less than attenuation due to rain, and that the attenuation due to dry snow may be neglected in microwave band. However, the attenuation due to wet snow is higher. Some results of attenuation due to hail are presented in (Ishimaru, 1978). In this chapter, our attention would be concentrated on the attenuation due to rain, clouds, and fog. The variation of the radio refractivity will be the object of our investigation presented there as well.


2. Attenuation due to rain

The electromagnetic wave attenuation due to rain (the rain attenuation) is one of the most noticeable components of excess losses, especially at frequencies of 10 GHz and above (Freeman, 2007). The methods of prediction of the rain attenuation can be grouped into two groups: the physical (exact) models and the empirical models. The physical models attempt to reproduce the physical behaviour involved in the attenuation processes while the empirical methodologies are based on measurement databases from stations in different climatic zones within a given region. The empirical methods are used widely and frequently with the best success (Emiliani et al., 2004). Two main causes of attenuation are scattering and absorption. When the wavelength is large compared to the size of raindrop, scattering is predominant. Conversely, when the wavelength is small compared to the raindrop’s size, attenuation due to absorption is predominant (Ivanovs & Serdega 2006). Water molecules are dipoles. The raindrop’s dipoles have the same time variation as the electromagnetic waves and therefore act as an antenna, which re-radiates the electromagnetic wave energy. Hence, a raindrop becomes an “antenna” with low directivity. Consequently, some energy is reradiated in arbitrary directions giving a net loss of energy in the direction towards the receiver (Ivanovs & Serdega 2006). Water is a loss-making dielectric medium. The relative dielectric constant of water is high, compared to the dielectric constant of the surrounding air. It depends on temperature and the operating frequency of the radio system. The specific heat of the water is high. Therefore, water absorbs a large amount of warmth, while warms itself. The surface tension of water is high. This is the reason why the molecules of water are holding together. One of the problems in prediction of electromagnetic wave power losses is description of shape of the raindrop. It depends on the size of droplet. It is known, that only very small droplets are like spheres. Such droplets form in clouds, as water vapour condenses on the nuclei of condensation. Further, these droplets grow by coalescence. Shape of the raindrops, that are larger than 1 mm in diameter, is no more spherical. They are not tear-shaped, as it commonly presented in pictures. The shape of falling large raindrops is more like a hamburger shape. Therefore, horizontally polarized waves suffer greater attenuation than vertically polarized waves (Freeman, 2007).

As mentioned above, the water molecules are polar ones. Those molecules rotate in such way that positive part of one molecule would be as near as possible to the negative part of another molecule. Therefore, molecules are rotating, hammering one on another and heating (Tamošiūnaitė et al., 2010a). The water molecule also rotates when a negative charge is brought near to it. The fields of electromagnetic wave vary up as time goes and force the water molecules to rotate respectively to the variation of fields.

2.1. Specific rain attenuation

One of the most widely used rain attenuation prediction methods is an empirical relationship between the specific rain attenuation α[dB km-1] and the rain rate R[mm·h-1] (Freeman, 2007, Rec. ITU-R P.838-3, 2005):


where aand bare functions of operating frequency fand rain temperature t; the value of R[mm·h-1] is for an exceedance of 0.01% of the time for point rainfall rates with an integration time of one minute. The coefficients aand b(coefficients ah and bh to be used for horizontal polarized waves; coefficients av and bv to be used for vertical polarized waves) are presented in (Freeman, 2007; Recommendation ITU-R P. 838-3, 2005).

2.1.1. Rain rate

In determination of the rain attenuation, the main parameter is rain rate R, which is expressed in [mm·h-1]. Gauges at the surface measure the accumulation of rain–water (flux) in a known time interval and report the result as a rain rate (accumulation per unit time) averaged over some measurement or aggregation interval (Crane, 1996). The rain rate can be described as the thickness of the precipitation layer, which felled down over the time period of one hour in the case when the precipitation is not evaporated, not soaked into the soil, and is not blown away by the wind (Tamošiūnaitė et al., 2010a). The evaluation of R–value is the first step in the rain attenuation prediction. The rain attenuation depends on the meteorological conditions in the considered localities. This is the reason to analyze the rain attenuation in particular locations (eg. country, city, climatic region).

First attempts to predict the rain attenuation under Baltic region climate conditions are described in (Tamošiūnas et al., 2005, 2006; Ivanovs & Serdega, 2006; Zilinskas et al., 2006, 2008). It was mentioned in (Ivanovs & Serdega, 2006), that rain events produce unavailability of microwave link, which sometimes lead operators to economical losses or even license loosing.

The significant differences in annual, seasonal, monthly, and daily amounts of rainfall are observed in localities of Lithuania. The noticeable local differences of rainfall amounts are characteristic of Lithuania as well. The precipitation amount is probably the most changeable meteorological index on Lithuania’s territory. It varies from 901 mm in Šilalė district to 520 mm in Pakruojis district (Bukantis, 2001). No month of a year could be described as “an average month” in Lithuania. This is the reason to revise the suitability of the models that derived under climatic conditions other than Lithuanian ones. The models using only annual amount of rainfall was analyzed in (Tamošiūnas et al., 2005). Considering the peculiarities of Lithuania’s climate, the change in (Chebil et al., 1999) model was made. This new model for the electromagnetic wave attenuation due to rain medium in atmosphere for the first time has been presented in (Tamošiūnas et al., 2006). Calculation of radio wave attenuation due to rain using annual precipitation and heavy rainfall data is described in (Zilinskas et al., 2006). The heavy rainfall events and showers with thunderstorms occur during the warm season (from May to September) in Lithuania.

2.1.2. Integration time

As was mentioned above, the R-values are expressed in [mm·h-1]. However, time intervals between the readings of rainfall amount in many cases must be much shorter. Those intervals are called the integration time τ. In (Ivanovs & Serdega, 2006; Tamošiūnas et al, 2007; Tamošiūnaitė et al., 2010a) it was mentioned, that the period of time between the readings of the rainfall amount values is a very important parameter, because it can significantly change the R-value. High R-values “hides” when τis long.

Consider an example. There were raining. The duration of the rain was 5 minutes. The total amount of the precipitation was 5 mm. It did not rain during remaining 55 minutes of one hour. Thereby, if we would count the average R-value for that hour (τ=60min.), it would be equal to 5 mm·h-1. But if we would count the R-value for every minute of that hour, we would find that R-values are much higher. Consider that in every of those 5 rainy minutes the amount of the precipitation was 1 mm. Consequently, for each of those 5 minutes the R-value would be 60 mm·h-1. That is why the average R-values are unreliable.

In Lithuania, the τvalues must be as small as possible (Tamošiūnaitė et al., 2010a).

2.1.3. “One-minute” rain rate

Almost all rain attenuation methods require “one–minute” rain rate value. The “one-minute” rain rate value R(1 min.) is expressed in [mm·h-1]. R(1 min.)-value can be defined as the R–value for 0.01% of time of the year, obtained using the rainfall amount value, which was measured inτ=1min and multiplied by 60 (Karasawa & Matsudo, 1991).

However, in many instances data collection is oriented toward agricultural and hydrological purposes, for which annual, monthly, daily, and less commonly, 3– and 6–hourly totals are collected. Therefore the models for conversion of R(τ min.)-values into R(1 min.)-values are used. A review of models for estimation of 1 min rainfall rates for microwave attenuation calculations are presented in (Tattelman & Grantham, 1985).

One of such conversion models was presented in (Moupfouma &Martin, 1995):


where R(1min.)is the “one-minute“ rain rate value, R(τmin.)is the rain rate value measured in τminutes (τ1min.).

In (Zilinskas et al, 2008) another model (4) for calculation of the R(1min.)-value was presented. That model was derived on the basis of model presented in (Rice & Holmberg, 1973) in accordance with the peculiarities of Lithuanian climate.


where MVIXis amount of rainfall which precipitated in May-September, tis the number of hours in a year when the value of rain rate could be equal or exceed the R(1min.)-value.

According to data that was collected in Lithuanian weather stations and (4) formula, the average R(1min.)-value for Lithuanian territory was calculated. That value is 60.23 mm·h-1. This value is double the value, which is suggested by ITU-R (Tamošiūnaitė et al., 2010a).

According to (1) formula, the values of coefficients aand b(presented in Freeman, 2007), and the value of R(1min.)= 60.23 mm·h-1, the dependency of the average specific electromagnetic wave attenuation due to rain, α,on the operating frequency fwas estimated. The results are shown in Fig. 1.

Figure 1.

The dependency of the average specific electromagnetic wave attenuation due to rainαon the operating frequencyf,in Lithuania.

2.1.4. Worst month statistics

The “Worst-month” model was proposed by ITU-R in (Rec. ITU-R P.481-4, 2005). This model is a supplement of the “One-minute” models, which were explained above. In “One-minute” models a lot of precipitation data must be collected and calculated. Furthermore, majority of those models are appropriate only in cases when the reliability of the radio wave system must be equal 99.99%. The main advantage of the “Worst-month” model is that only the worst-month statistics must be collected. Furthermore, the “Worst-month” model is appropriate in cases when the required reliability of the radio system is other than 99.99%.

The worst-month is the month (or 30 days period) from a year (or twelve consecutive calendar months), during which the threshold is exceeded for the longest time. This month is not necessarily the same month in different year. The fraction of time when the threshold value of rain rate (so, and rain attenuation value) was exceeded is identical to probability that the threshold value of rain rate would be exceeded (Crane, 1996).

The average annual worst-month time percentage of excess, pm, is proportional to the average annual time percentage of excess, p, in such relation:


where Qis the conversion factor; pm [%] and p[%] must refer to the same threshold levels (the same rain rate value).

The conversion factor Qis a two-parameters (Q1, β) function of p. In most cases a high reliability of the radio system is required (p3%). Then Qcan be expressed as (Rec. ITU-R P.481-4, 2005):


For global planning purposes the following values of the parameters Q1 and βmay be used: Q1=2.85and β=0.13(Rec. ITU-R P.481-4, 2005).

For global rain rate applications, the following values for the parameters Q1 and βshould be used: Q1=2.82andβ=0.15, for tropical, subtropical and temperate climate regions with frequent rain; Q1=4.48andβ=0.11, for dry temperate, polar and desert regions. Yet ITU-R recommends that more precise values of Q1 and βshould be used where possible.



and (6), consequently:


Mark 1Q1=qand11β=ξ, then:


According to (2), (3) and annual data, the relation between pand R(1min.)can be found. This relation could be compared to the relation calculated according to (8) and ITU-R suggested Q1 and βvalues. According to Lithuanian climate, the values Q1=2.82and β=0.15should be appropriate.

For example, we evaluated the “Worst-month” model in Vilnius, the capital of Lithuania. The results are shown in Fig. 2. As can be seen, the values Q1=2.82and β=0.15are

Figure 2.

The correlation between the real, calculated and corrected values ofp(in Vilnius).

appropriate only in cases whenR(1min.)38mm·h-1. When R(1min.)38mm·h-1, the calculated values are apparently distant from the real values. Therefore, the values of Q1 and βmust be corrected. The best correlation is when in (6) there are q=0.5andξ=1.03. Consequently, the corrected Q1 and βvalues should be Q1=2andβ=0.03. But still, as can be seen in Fig. 2, the corrected values are only correct when R(1min.)30mm·h-1. Furthermore, whenR(1min.)34mm·h-1, the values Q1=2.82and β=0.15are more proper than Q1=2andβ=0.03. As a result, in cases whenR(1min.)34, the values Q1=2and β=0.03should be used, and in cases when R(1min.)34mm·h-1, the ITU-R suggested values Q1=2.82and β=0.15may be used.


3. Attenuation due to clouds

The effect of rain attenuation is greater than that of clouds in many cases, but clouds occur more often than rain. In clouds, water droplets are generally less than 0.01 cm in diameter (Freeman, 2007). In (Altshuler & Mart, 1989), it was mentioned that cloud attenuation was primarily due to absorption by the cloud droplets, and scattering losses were secondary. With increase in operating frequency the attenuation due to clouds also increases, but as the temperature of the clouds decreases the attenuation value increases (Sarkar et. al., 2005).

On average, the clouds cover more than 50% of the territory of Lithuania. According to the data of its weather stations, November and December are the cloudiest months. The clearest sky is in May and June. There are about 100 overcast days in the year.

3.1. Liquid water content

The liquid water content Mis one of the most important parameters of the clouds. Mdescribes the mass of water drops in the volume units of the cloud. It has been mentioned in (Freeman, 2007) that the specific cloud attenuation αC[dB/km] is a function of the liquid water content M[g/m3], the frequency f, and the temperature within the cloud T. The measurements of Mat a point in space or averaged over a radio wave path are very complicated. Direct methods for measuring Mconsists of extracting a known volume through a cotton pad or of rotating cups in an impeller apparatus, both to be weighed; also, resistance changes can be measured with a hot wire probe attached to an aircraft flying through clouds (Liebe et al., 1989). The liquid water content in the cloud varies in a wide range. In most of the cloud attenuation models, it is required to know the value of M.

The climate conditions (humidity, temperature, etc.) and cloud morphology are different over various localities of several regions; accordingly, the liquid water contents differ within the clouds as well. This factor must be considered when analyzing rain attenuation and cloud attenuation. Our first attempt to determine the specific cloud attenuation under the Lithuanian climatic conditions is presented in (Tamošiūnaitė et al., 2008; Zilinskas et al., 2008). The humid weather predominates over the year in Lithuania.

3.2. Calculation of the specific cloud attenuation

The specific cloud attenuation is a function of clouds’ liquid water content and a coefficient, which is a function of frequency and temperature. In this case, the main problem is the value of clouds’ water content, because the direct measurements at a point in space are problematic. In cases when such data is unavailable, models that require only the meteorological parameters, measured at ground level, can be used. These models are based on the fact that the condensation is possible when the water vapour concentration exceeds the saturation density at the temperature, which is prevailing at that height. The water vapour density can be estimated from the humidity measurements carried out at ground level. The cloud’s water content value can be estimated as the difference between water vapour concentration and saturation density at cloud temperature. The specific cloud attenuation is, unlike the case of rain, independent of drop–size distribution (Freeman, 2007). Several cloud attenuation models were developed. In (Freeman, 2007), the specific cloud attenuation was expressed as the function of liquid water content M:


where KC is the attenuation constant.

The attenuation constant KC is the function of frequency fand temperature T. The values of KC for pure water droplets are presented in (Freeman, 2007). The values of KC for salt-water droplets (over the sea and ocean surfaces) are higher. The necessity to know Mvalue is limiting the direct use of relationship (10).

Often there are no possibility to measure the liquid water content and temperature within the clouds. In such cases the methods that require only meteorological parameters measured at the ground level may be used. The basic idea of such models (Dintelmann &Ortgies, 1989) is that the water vapour in the atmosphere would lead to the formation of clouds whenever there would be a possibility for condensation at some height habove ground level. There is also mentioned that the condensation is possible when the water vapour density ρexceeds the saturation density ρs at temperature Tprevailing at that height. It is assumed that the water vapour density ρcan be estimated from humidity measurements carried out at ground level.

The height at which cloud exists is very important for accurate determination of results of attenuation due to clouds (Sarkar et al., 2005). It was assumed in (Ito, 1989, as cited in Dintelmann &Ortgies, 1989) that clouds are created starting in the vicinity of the height h, and h[km] follows ground temperature T0 [K] as:


Relation (11) is based on analysis of temperature profiles in rain and on the Aerological Data of Japan and we have specified the applicability of this relation in the territory of Lithuania.

The condensed water content Mis estimated as the difference between ρand saturation density ρsat cloud temperature (Dintelmann & Ortgies, 1989):


where ρs [g/m3] is the saturated vapour density.

It is assumed that clouds are formed when M>0. As mentioned above, the determination of the water content value Mis complicated. Its values differ in each group of the clouds (the clouds are grouped according to their shape, height, and structure). In our calculations, the main problem was determination of M. According to (Dintelmann & Ortgies, 1989), the values of water vapour density ρat the height hcan be estimated from the equation of state, assuming an adiabatic process:


where ρ0 is the water vapour density at the ground level, T0 is the ground level temperature, Tis the absolute temperature in the vicinity of h, denotes the specific heat ratio which is 4/3 for the water vapour molecule, μis the water molar mass, gis the acceleration due to gravity, his the height, and Ris the fundamental gas constant. The values of ρ0 can be determined by using known relations (Freeman, 2007).

We assume that the clouds are created starting in the vicinity of the height h. We determine the values of hby using relation (11) or the data of the dew point temperature, temperature at the ground level, and the temperature gradient of 6.5˚C/km (Rec. ITU-R P.P.835-3, 2004).

The values of hobtained here we compared to the cloud base height values measured at the weather stations (see Table 1). The analysis of the cloud cover over the localities of Lithuania data shows that the relationship (11) can be used only in the cases when the middle or high clouds are formed over those localities.

T0 [K]Cloud base height (data of weather stations)Cloud base height (equation 11)

Table 1.

Temperature at the ground level and the values of the cloud base heights (data of weather station) in Vilnius in April 2007, as well as the height hdetermined using equation (8) (Tamošiūnaitė et al., 2008).


4. Attenuation due to fog

The influence of the fog on the attenuation of the electromagnetic waves can to lead to the perturbation of the wireless communication. In (Chen et al., 2004), it was mentioned that fog may be one of dominant factors in determination of the reliability of millimeter wave systems, especially in coastal areas, where dense moist fog with high liquid water content happen frequently. Fog results from the condensation of atmospheric water vapour into water droplets that remain suspended in air (Freeman, 2007). Moist fog frequently appears over the localities of Lithuania (Tamosiunas et al., 2009). There are several meteorological mechanisms for determination whether fog will form and of degree of its intensity. The physical mechanism of the formation of the fog can be reduced to three processes: cooling, moistening, and vertical mixing of air parcels with different temperatures and humidity (Duynkerke et al., 1991). All three processes can occur, although one meteorological mechanism may dominate. This circumstance leads to the different types of the fog. In (Galati et. al., 2006), the fog is classified in four types: strong advection fog, light advection fog, strong radiation fog, and light radiation fog.

The calculation methods for determination of fog attenuation are used in many cases. The propagation properties for microwave and millimeter–wave frequencies at the foggy air conditions were examined in (Liebe et. al, 1989). The values of the specific attenuation were derived from a complex refractivity based on the Rayleigh absorption approximation of Mie’s scattering theory. In (Liebe et. al, 1989), the particle mass content and permittivity, which depends on the frequency and the temperature, were key variables. Attenuation due to fog is a complex function of the particle size distribution, density, extent, index of refraction, and wavelength (Altshuler, 1984). Normalized fog attenuation directly, given only the wavelength and fog temperature is presented in (Altshuler, 1984):


where Ais attenuation in [(dB/km)/(g/m3)], λis wavelength in [mm], tis temperature in [ C]; the relation (14) is valid only for 3 mm< λ<3 cm and –8 C< T< 25 C.

It was mentioned in (Altshuler, 1984], that the total fog attenuation could be obtained by multiplying the normalized attenuation by the fog density in [g/m3] and the fog extent in [km]. In (Zhao &Wu, 2000), it was mentioned that fog is often characterized by the visibility and the visibility is defined as the greatest distance at which it is just possible for an observer to see a prominent dark object against the sky at the horizon.

Attenuation due to fog can be expressed in terms of the water content M, and the microstructure of the fog can be ignored (Galati et al., 2000). In (Altshuler, 1984), the empirical formula for fog visibility as a function of fog density was derived:


where Vis the visibility in [km] and Mis the liquid water content in [g/m3].

It was mentioned in (Altshuler, 1984), that the empirical formula (15) is valid for drop diameter between 0.3 μm and 10 μm. For the case of dense haze or other special type fogs, it is recommended to replace the coefficient 0.024 with 0.017 (Altshuler, 1984). If the visibility data are available, but the fog density data are not available, the following expression may be used (Altshuler, 1984):


In (Chen et al., 2004; Galati et al., 2006; Recommendation ITU-R PN 840-4, 2009), based on the Rayleigh approximation, the specific attenuation due to the fog αfog has been written as:


where Kis specific attenuation coefficient.


where θ=300/T, fis frequency, and Tis temperature [K].

V, kmM, g/m3

Table 2.

The values of visibility Vmeasured in the localities of Lithuania and the values of fog water content M(Tamosiunas et al., 2009).

The values of the visibility measured in the localities of Lithuania and the values of fog water content Mdetermined using (16) are presented in Table 2. The highest value of the specific fog attenuation determined using M-data presented in Table 2 was 0.59 dB/km.

In (Naveen Kumar Chaudhary et al., 2011), it was concluded, that the link reliability can be improved by increasing the transmission power or using high gain directional antennas in the cases when the foggy conditions occur and the visibility is less than 500 meters. For the same value of visibility, the fog attenuation decreases when the temperature increases (Naveen Kumar Chaudhary et al., 2011).


5. Radio refractive index and its variability

The atmospheric refractive index is the ratio of the velocity of propagating electromagnetic wave in free space and its velocity in a specific medium (Freeman, 2007). The value of the atmosphere’s refractive index is very close to the unit. Furthermore, changes of the refractive index value are very small in time and space. In the aim to make those changes more noticeable, the term of refractivity is used. It is a function of temperature, atmospheric pressure and partial vapour pressure. The value of the refractivity is about million times greater than the value of refractive index.

In design of the radio communication networks, it is important to know the atmospheric radio refractive index. The path of a radio ray becomes curved when the radio wave propagates through the Earth’s atmosphere due to the variations in the atmospheric refractivity index along its trajectory (Freeman, 2007). Refractivity of the atmosphere affects not only the curvature of the radio ray path but also gives some insight into the fading phenomenon. The anomalous electromagnetic wave propagation can be a problem for radars because the variation of the refractive index can induce loss of radar coverage (Norland, 2006). In practice, the propagation conditions are more complicated in comparison with the conditions predictable in design of radio system in most cases.

The anomalous propagation is due to the variations of the humidity, temperature and pressure at the atmosphere that cause variations in the refractive index (Norland, 2006). The climatic conditions are very changeable and unstable in Lithuania (Pankauskas & Bukantis, 2006). The territory of Lithuania belongs to the area where there is the excess of moisture. The relative humidity is about 70% in spring and summer while in winter it is as high as 85 – 90% (Bagdonas & Karalevičienė, 1987). Lithuanian climate is also characterized by large temperature fluctuations. Difference between the warmest and coldest months is 21.8 C (Pankauskas & Bukantis, 2006). It was noted in (Priestley & Hill, 1985; Kablak, 2007) that even small changes of temperature, humidity and partial water vapour pressure lead to changes in the atmospheric refractive index. In (Zilinskas et al., 2008), the measurements of these meteorological parameters were analyzed in the different time of year and different time of day. The values of the refractive index have been determined by using measured meteorological data. In (Žilinskas et al., 2010), it was mentioned that seasonal variation of refractivity gradient could cause microwave systems unavailability.

5.1. Calculation of radio refractivity

As mentioned above, the value of the radio refractive index, n,is very close to the unit and changes in this value are very small in the time and in the space. With the aim to make those changes more noticeable, the term of radio refractivity, N,is used (Freeman, 2007; Rec. ITU-R P. 453-9, 2003):


According to the recommendation of ITU –R (Rec. ITU-R P. 453-9, 2003):


where T[K] is a temperature; p[hPa] is the atmospheric pressure; e[hPa] is partial water vapour pressure. The refractivity is expressed in N– units.

It was mentioned in (Freeman, 2007; Rec. ITU-R P. 453-9, 2003), that expression (21) may be used for all radio frequencies; for frequencies up to 100 GHz, the error is less than 0.5%. There are two terms (the “dry term” and the “wet term”) in relationship (21).

The values of the refractivity Nin Lithuania were determined by using (21). The data of temperature, humidity, and atmospheric pressure were taken from a meteorological data website (

Figure 3.

The dependences of averageN– values on the time of day in cities of Lithuania: Vilnius (curve 1), Mažeikiai (curve 2), Kaunas (curve 3), and Klaipėda (curve 4) in July 2008 (Valma, et al., 2010).

The dependences of average N–values on the time of day in cities of Lithuania are presented in Fig. 3. As can be seen, the behaviours of those dependences at the diurnal time are similar in all localities that are situated in the Continental part of Lithuania (Vilnius, Kaunas and Mažeikiai) and slightly different in Seacoast (Klaipėda). The climate of Klaipėda is moderate and warm (Pankauskas &Bukantis, 2006; Zilinskas et al., 2008). The climate of Continental part of Lithuania is typical climate of the middle part of the Eastern Europe. This may explain the difference between the daily variations of Nin Klaipėda and in other localities analyzed here. In Lithuania, the highest N-values were in July.


6. Conclusions

The main models for calculation of electromagnetic wave attenuation due to atmosphere humidity were revised. In Lithuania, when the reliability of the radio system of 99,99% is required, the R(1min.)-value is R(1min.)=60.23mm/h. It is twice the ITU-R recommended value. The dependency of the average specific electromagnetic wave attenuation due to rain on the operating frequency (0-100 GHz) was determined. The attenuation of horizontally polarized electromagnetic waves is greater than the attenuation of vertically polarized electromagnetic waves. In cases when the required reliability of the radio system is other than 99,99%, the “Worst-month” model can be used. However, for small R(1min.)-values the parameters of that model should be corrected. In Vilnius, the city of Lithuania, when R(1min.)34mm/h, ITU-R recommended values Q1=2.82and β=0.15could be used. In cases when R(1min.)34mm/h, the corrected values Q1=2and β=0.03are more appropriate.

The main problem of models for calculation of electromagnetic wave attenuation due to clouds and fog is the required value of liquid water content. In Lithuania it is impossible to gather such meteorological information. Therefore, models excluding or calculating the liquid water content were revised. The variations of the atmospheric humidity, temperature and pressure can cause the fluctuations of the atmospheric refractive index. In Lithuania, the atmosphere refractive index fluctuates most in July. The variations of N in diurnal time are similar in all localities that are situated in the Continental part of Lithuania and slightly different in Seacoast.


  1. 1. AltshulerE. E.1984A simple expression for estimating attenuation by fog at millimeter wavelengths. IEEE Transactions on Antennas and Propagation,327(July 1984),757758,0001-8926X
  2. 2. AltshulerE. E.MarrR. A.1989Cloud attenuation at millimeter wavelengths. IEEE Transactions on Antennas and Propagation,3711(November 1989),14731479,0001-8926X
  3. 3. BagdonasA.KaralevičienėR.1987The Reference Book of Agrometeorologist,Mokslas, Lithuania (in Lithuanian).
  4. 4. BhattacharyyaS.DanM.SenA. K.2000Modelling of drop size distribution of rain from rain rate and attenuation measurements at millimetre and optical wavelengths. International Journal of Infrared and Millimeter Waves,2112(2000),20652075,1572-9559
  5. 5. BukantisA.2001Climate fluctuations in Lithuania against a background of global warming. Acta Zool. Lituanica,1122001),113120,1392-1657
  6. 6. Characteristics of precipitation for propagation modeling.20032003, Recommendation ITU-R837834
  7. 7. ChebilJ.RaihmanT. A.1999Development of the one-minute Rain Rate Contour Maps for Microwave Applications in Malaysia Peninsula. Electronics Letters,35201999)177217740013-5194
  8. 8. ChenH.DaiJ.LiuY.2004Effect of fog and clouds on the image quality in millimeter communications. International Journal of Infrared and Millimeter Waves,255(May 2004),749757,1572-9559
  9. 9. Conversion of annual statistics to worst-month statistics.,2005Recommendation ITU-R481484, 2005
  10. 10. CraneR. K.Electromagnetic Wave Propagation Through Rain,1996John Wiley&Sons, Inc., New York, Chichester, Weinheim, Brisbane, Singapore, Toronto,100471613762
  11. 11. DintelmannF.OrtgiesG.1982Semiempirical model for cloud attenuation prediction. Electronics Letters,2522(October 1982),14871488,0013-5194
  12. 12. DissanayakeA.AllnuttJ.HaidaraF.1997A prediction model that combines rain attenuation and other propagation impairments along Earth-satellite paths. IEEE Transactions on Antennas and Propagation,4510(October 1997),15461558,0001-8926X
  13. 13. EmilianiL. D.AgudeloJ.GutierrezE.RestrepoJ.Fradique-MendezC.2004Development of rain-attenuation and rain-rate maps for satellite system design in the Ku and Ka bands in Colombia. IEEE Transactions on Antennas and Propagation Magazine,4665468,1045-9243
  14. 14. FreemanR. L.2007Radio System Design for Telecommunications,Third Edition, John Wiley&Sons,978-0-47175-713-9New York
  15. 15. GalatiG.DalmassoI.PavanG.BrogiG.2006Fog detection using airport radar.Proceedings of IRS 2006 International Radar Symposium,209212, Krakow, Poland May 24-26, 2006,
  16. 16. IshimaruA.1978Wave propagation and Scattering in Random Media,Academic Press,100123747015York, San Francisco, London
  17. 17. ItoS.1989Dependence of 0º C isotherm heoght on temperature at ground level in rain. IEICE Transactions,E722February),98100
  18. 18. IvanovsG.SerdegaD.2006Rain intensity influence on to microwave line payback terms. Electronics and Electrical Engineering,670),6064,0000-1392- 1215
  19. 19. KablakN. I.2007Refractive index and atmospheric correction to the distance to the Earth’s artificial satellites, Kinematics and Physics of Celestial Bodies,2328488,1934-8401
  20. 20. KarasawaY.MatsudoT.1991One-minute rain rate distributions in Japan derived from A Me DAS one-hour rain rate data, IEEE Transactions on Geoscience and Remote Sensing,296890894,0196-2892
  21. 21. LiebeH. J.ManabeT.HuffordG. A.1989Millimeter-wave attenuation and delay rates due to fog/cloud conditions. IEEE Transactions on Antennas Propagation,3712(December),16171623,1001-8926X
  22. 22. MoupfoumaF.MartinL.1995Modelling of the rainfall rate cumulative distribution for the design of sattelite and terrestrial communication systems. International Journal of Satellite Communications and Networking,32(March/April 1995),105115,1542-0973
  23. 23. NaveenKumar.ChaudharyTrivedi. D. K.RoopamGupta.2011Radio link reliability in Indian semi-desert terrain under foggy conditions. International Journal of Latest Trends in Computing,21(March 2011),4750, E-2045-5364
  24. 24. NorlandN.2006Temporal Variation of the Refractive Index in Coastal waters.Proceedings of IRS 2006, International Radar Symposium,221224, Krakow, Poland, 24- 26 May 2006,8-37207-621-9
  25. 25. PankauskasM.BukantisA.2006The dynamics of the Baltic Sea Region climate humidity. Annales Geographicae,3912006),514,1822-6701
  26. 26. PriestleyJ. T.HillR. J.(1985), Measuring High-Frequency Refractive Index in the Surface Layer.Journal of Atmospheric and Oceanic Technology,222233251,0739-0572
  27. 27. Reference Standard Atmospheres.2004Draft Revision to Recommendation ITU-R8358338353
  28. 28. RiceP.HolmbergN.1973Cumulative time statistics of surface-point rainfall rates.IEEE Transactions on Communications,212110(October 1973),11311136,0090-6778
  29. 29. SarkarS. K.KumarA.AhmadI.GuptaM. M.2006Cloud morphology over three Indian tropical stations for Earth space communication.International Journal of Infrared and Millimeter Waves,27277(July 2006),10051017,1572-9559
  30. 30. Specific attenuation model for rain use in prediction methods,2005Recommendation ITU-R838833, 2005
  31. 31. TamošiūnaitėM.TamošiūnasS.TamošiūnienėM.ŽilinskasM.2008Influence of clouds on attenuation of electromagnetic waves. Lithuanian Journal of Physics,4816572,1648-8504
  32. 32. TamošiūnaitėM.TamošiūnienėM.GruodisA.TamošiūnasS.2010Prediction of electromagnetic wave attenuation due to water in atmosphere. 1. Attenuation due to rain. Innovative Infotechnologies for Science, Business and Education,29(2010),310,2029-1035
  33. 33. TamošiūnaitėM.TamošiūnasS.DaukšasV.TamošiūnienėM.ŽilinskasM.2010Prediction of electromagnetic waves attenuation due to rain in the localities of Lithuania.Electronics and Electrical Engineering,9105), ,912,1392-1215
  34. 34. TamošiūnasS.ŽilinskasM.NekrošiusA.TamošiūnienėM.2005Calculation of radio signal attenuation using local precipitation data. Lithuanian Journal of.Physics,455353357,1648-8504
  35. 35. TamošiūnasS.ŽilinskasM.ŠileikaM.TamošiūnienėM.2006Revised model of attenuation of electromagnetic waves due to rain. Lithuanian Journal of Physics,464433436,1648-8504
  36. 36. TamosiunasS.TamosiunaiteM.ZilinskasM.TamosiunieneM.2009The influence of fog on the propagation of the electromagnetic waves under Lithuanian climate conditions. PIERS Online,56576580, 321-325,1931-7360
  37. 37. The radio refractive index: its formula and refractivity data.,2003Recommendation ITU-R453459
  38. 38. ValmaE.TamošiūnaitėM.TamošiūnasS.TamošiūnienėM.ŽilinskasM.2010Determination of radio refractive index using meteorological data.Electronics and Electrical Engineering,10106),125128,1392-1215
  39. 39. ZhaoZ.WuZ.2000Millimeter- wave attenuation due to fog and clouds. International Journal of Infrared and Millimeter Waves,2110(2000),16071615,1572-9559
  40. 40. ZilinskasM.TamosiunasS.TamosiunieneM.2006Calculation of radio signal attenuation using annual precipitation and heavy rainfall data.Proceedings of EMC 2006 18 th International Wroclaw Symposium and Exhibition on Electromagnetic Compatibility,490493,837025947Wroclaw, Poland, June 28-30, 2006
  41. 41. ZilinskasM.TamosiunaiteM.TamosiunasS.TamosiunieneM.2008The influence of the climatic peculiarities on the electromagnetic waves attenuation in the Baltic Sea region. PIERS Online,43321-325,1931-7360
  42. 42. ZilinskasM.TamosiunaiteM.TamosiunasS.TamosiunieneM.2009Calculation of electromagnetic waves attenuation due to rain for various percentages of time.Proceedings of PIERS 2009Progress in Electromagnetics Research Symposium,541545, Beijing, China, March 23-27 2009,1559-9450978-1-93414-208-0

Written By

Milda Tamošiunaite, Stasys Tamošiunas, Mindaugas Žilinskas and Milda Tamošiuniene

Submitted: November 13th, 2010 Reviewed: March 3rd, 2011 Published: June 21st, 2011