Completo_DEF.qxd 1333 ANNALS OF GEOPHYSICS, VOL. 47, N. 4, August 2004 Ten years analysis of Tropospheric refractivity variations Stergios A. Isaakidis and Thomas D. Xenos Department of Electrical and Computer Engineering, Aristotle University of Thessaloniki, Greece Abstract The refractivity variations of the troposphere are responsible for various effects on radio wave propagation, such as refraction, bending, radio-station interference, etc. In this work, the refractivity variations of the Hel- lenic troposphere are studied using data from Helliniko airport of Athens/Greece. The data were analyzed us- ing various interpolation procedures, i.e. in a day-by-day manner for temperature and relative humidity, to transform the data according to a reference height common for the whole dataset and finally for the refractiv- ity N with respect to time, using piece-wise hermite interpolation polynomials for the low and medium alti- tudes and linear interpolation factors for the high altitudes. Since refractivity varies with time and height, two height independent basic parameters were computed and analyzed: the refractivity at station height N0 and scale height Hs. These parameters can be used to calculate the refractivity profile. Using statistical tools as the mov- ing average, the β0 parameter and the monthly mean values, together with the corresponding standard devia- tions, useful results were obtained for the variations of the refractivity with respect to observation hour, height, season, month, day and level. 1980; Anderson, 1995), refraction and fading of the electromagnetic waves (Flock, 1987), ducting phenomena (Patterson et al., 1994), etc. The index of refraction is defined as /c vr =fn = , where εr is the dielectric constant of the troposphere, c is the speed of light and v is the phase velocity of the wave in the medium. Because n is slightly greater than unity near the earth’s surface (1.00025-1.00040), in practice we use the scaled index of refraction N which is called refractivity and is defined as (Flock, 1987) ( ) . . . N n T p T e T e 1 10 77 6 5 6 3 75 106 2 5 # # = - = - + (1.1) where p is the total pressure in mbar, e is the water vapor pressure and T is the temperature in Kelvin degrees. If the Relative Humidity (RH %) and tem- perature are given, the water vapor pressure e can be computed with the following set of equations (Lear, 1980): 1333 Key words troposphere – refractivity – scale height – β0 parameter 1. Introduction The troposphere is the lower region of At- mosphere and it extends from the earth’s surface to an altitude of about 10 km at the middle lati- tudes (6 km at the poles and up to 18 km at the equator). The basic parameters that describe the tropospheric medium are pressure, temperature and relative humidity. Propagation in the tropo- sphere is affected by the variations of the refrac- tive index with time and height. These refractivi- ty variations are responsible for various phenom- ena in the wave propagation, such as range and elevation angle errors in radar acquisition (Lear, Mailing address: Dr. Stergios A. Isaakidis, Department of Electrical and Computer Engineering, Aristotle Univer- sity of Thessaloniki, GR-54 124 Thessaloniki, Greece; e- mail: sisaakid@otenet.gr 1334 Stergios A. Isaakidis and Thomas D. Xenos 1334 . . ( . ) x T T 35 85 17 2694 273 15 = - - (1.2) . ( )( )expe x mbars6 11s #= (1.3) .e e RH 100 s = (1.4) Theoretically speaking, the refractivity with re- spect to height is expressed by the exponential relation expN N h H s0= -_ i (1.5) where N0 is the refractivity at the station’s height, h is the height and Hs is the scale height. It is obvious that the parameters N0 and Hs are independent of height and since their changes will be reflected in N according to eq. (1.5), they constitute a useful tool for the study of the refractivity variations. The purpose of this work is to study the tro- pospheric refractivity variations over the Hel- lenic region. These variations describe the tro- pospheric conditions over the specific region of interest and can be taken into account as an em- pirical reference for troposheric wave propaga- tion, scale height estimations, bending of rays for terrestrial or satellite paths and so on. 2. Data analysis For the calculation of the refractivity varia- tions, data from Helliniko airport of Athens/ Greece (37°53′N, 23°43′E) were used. The data cover the period from 1991 to 1999 and the measured parameters are temperature, relative humidity and height at constant pressure levels. The measured data were provided by the Hel- lenic National Meteorological Service (HNMS) and come from observations with radio balloons launched twice a day (00 and 12 LT). The measured parameter values were classi- fied according to observation hour, day of the year and pressure levels and after a preliminary day-by-day statistical process according to the standard deviation of each pressure level, the data were cleared from any unrealistic values. Before calculating the refractivity indices, in- terpolation procedures were followed to fill in missing values, enrich with data the low alti- tudes, where the intensity of the tropospheric variations is higher and also to create smoother curves of the initial parameters. The first inter- polation procedure was made in a day-by-day manner according to the reference pressure lev- els and was the most important in the sense that it was applied directly to the initial measured data and provided most of the interpolated val- ues, forming smoothed and enriched curves of temperature, relative humidity and height ver- sus pressure. The next interpolation procedure was also applied each day. Its purpose was to transform the data according to a reference height, common for the whole dataset. This procedure also enriched the low altitude values and also detected and corrected days where ab- normal fluctuations of the original data were present. After the analysis of the data according to the above procedures, the refractivity index N was calculated using eqs. (1.1) to (1.4). Fi- nally, an interpolation procedure for N accord- ing to time was followed for each altitude lev- el, in order to fill in missing values that could not be interpolated in the previous methods. For all the interpolation procedures, piece- wise-cubic-hermite interpolating polynomials were used, except for high altitudes (higher than 10 km), where linear interpolation functions have been applied. This was due to the fact that at high altitudes, when a large number of data were miss- ing, the cubic functions in many cases returned unrealistic results. Spline interpolation functions instead of cubic were also tested, but they result- ed in strong oscillations about the measured data, especially in cases where many consecutive val- ues were missing, as is the case for relative hu- midity. After calculation of refractivity, the scale height Hs was determined for each day by apply- ing a least squares nonlinear fitting. The analysis of the data by means of the pre- viously described procedures and the calcula- tion of refractivity were followed by a statistical evaluation in order to assess their variability. 3. Results and discussion The statistical treatments applied are the second order 31-days moving averages of N0 1335 Ten years analysis of Tropospheric refractivity variations 1335 and Hs versus days and their monthly mean val- ues, together with the corresponding standard deviations. Although calculated N0 and Hs with respect to time form a time series with periodic variations, moving averages eliminate strong fluctuations by low pass filtering the data. Also, the β 0 parameter was determined for each month for the two observation hours of 00 and 12 LT. This parameter is defined as the percent- age of time that the refractive index gradient in the first 100 m above ground exceeds 100 N- units/km and it represents the probability that anomalous propagation may occur. Figure 1 shows the second order 31-days moving average of N0 versus all the days of the period 1991-1999 and for both observation hours (00 and 12 LT). In the horizontal axis the serial number of the day is presented with day 1 corresponding to the 1st January 1991 and Day 3287 to 31st December 1999. These diagrams constitute a de-noised version of the original calculated refractivity and they show clearly the cycle variations of N0 for all years (separated with vertical grid lines). It can be seen that in general, the refractivity increases from January to June, where it appears to reach a first peak. During the summer period, refractivity levels remain high, but a decrement usually occurs during July, reaching a minimum level in Au- gust. Then, the refractivity increases till Sep- tember-October, where it reaches a second peak, higher than the previous one. After October and for the whole winter period, refractivity shows a high rate of attenuation. From the same figure, it can be deduced that only the amplitude of re- fractivity is changing, being lower at 12 LT, while the refractivity variations with respect to time preserve their initial form. Figure 2 shows the second order 31-days moving average of Hs versus the days for the two observation hours of 00 and 12 LT. Comparing this diagram with fig. 1 it can be seen that the scale height follows the form of N0, but its varia- tions are lower. The results corresponding to 12 LT are similar to those of 00 LT since the scale height does not show time variations, except for its level which seems higher compared to 00 LT. Further information on seasonal refractivity variations can be extracted from fig. 3 which shows the monthly mean values of N0 versus months for the whole analyzed period. In this figure, the increment of the refractivity values during the summer period is clear. This diagram also shows comparative variations of No for each year. Although all the years do not follow the same level variations, the decrement from June to July is clear, except for 1996, when this phenomenon started on May. Except for the seasonal variations of refrac- tivity, diurnal variations are also present (fig. 4). This figure shows the standard deviations of the monthly values of N0 versus months. This diagram shows that the increase in refractivity Fig. 2. Second order moving average of Hs versus days (observation hour 00 and 12 LT). Fig. 1. Second order moving average of N0 versus days (observation hour 00 and 12 LT). 1336 Stergios A. Isaakidis and Thomas D. Xenos levels during the summer is followed by strong deviations within each month. In general, these deviations increase till July where they reach their maximum level and they start to decrease again during the autumn and winter periods. The identification of the monthly variations of the scale height is not as clear as the varia- tions of N0 (fig. 5). This figure shows the month- ly mean values of Hs versus month for the peri- od 1991-1999. It can be seen that in general Hs follows an inverse variation pattern compared to N0. It is decreasing during summer and a local peak appears in July but its deviations are very strong and they do not follow the distribution of the standard deviations of N0 (figs. 4 and 6). The results for the β 0 parameter are shown in figs. 7 and 8, for 00 and 12 LT respectively. By studying fig. 7 it can be seen that the prob- ability of anomalous propagation is very high during the summer months reaching a value of about 90%. During the other periods of the year the β 0 values are much lower but strong varia- tions can be observed in April, September and October where the β 0 values show a fluctuation Fig. 3. Monthly mean values of N0 versus months (observation hour 00 LT, years 1991-1999). Fig. 5. Monthly mean values of Hs versus months (observation hour 00 LT, years 1991-1999). Fig. 4. Standard deviations of the monthly values of N0 versus months (observation hour 00 LT, years 1991-1999). Fig. 6. Standard deviations of the monthly values of Hs versus months (observation hour 00 LT, years 1991-1999). 1337 Ten years analysis of Tropospheric refractivity variations of the order of 50%. For the 12 LT observation hour, it can be seen (fig. 8) that the β 0 values are in general higher during the summer period, but the annual variations are stronger with re- spect to 00 LT and thus the prediction of anom- alous propagation conditions is very difficult. The relation between N0 and Hs can be fur- ther revealed by comparing their standardized values, No’ and Hs’ respectively. The standardi- zation can be made using the equations / ( )Nstd= ( )N N Nmean0 0 0 0-l 7 A (3.1) / ( )Hstd= ( )H H Hmeans s s s-l 7 A (3.2) where the moving averages of N0 and Hs are used. In the above equations «mean» stands for the global mean value of the moving averages and «std» for their corresponding standard devi- ation. In fig. 9 the N′0 and – H′s are plotted with respect to days. It can easily be seen that their correlation is very high with respect to time (cor- relation coefficient ≈ 0.77), whereas N′0 shows weaker variations. REFERENCES ANDERSON, K.D. (1995): Radar detection of low-altitude targets in a maritime environment, IEEE Trans. Anten- nas Propag., 43, 609-613. FLOCK L.W. (1987): Propagation effects on satellite sys- tems at frequencies below 10 GHz, NASA Reference Publication Handbook, 1108 (2). LEAR M.W. (1980): Computing Atmospheric Scale Height for Refraction Corrections, (NASA Mission Planning and Analysis Division, Lyndon B. Johnson Space Center). PATTERSON, W.L., C.P. HATTAN, G.E. LINDEM, R.A. PAULUS, H.V. HITNEY, K.D. ANDERSON and A.E. BARRIOS (1994): Engineer’s Refractive Effects Prediction Systems (EREPS) Version 3.0, NRaD Technical Document 2648, May. (received February 02, 2004; accepted July 02, 2004) Fig. 9. Comparison of the N′0 and H′s variations (ob- servation hour 00 LT). Fig. 7. β 0 parameter versus months (observation hour = 00). Fig. 8. β 0 parameter versus months (observation hour = 12).