Refractivity
Refractivity
The variability and structure of tropospheric refractivity (N) over Nigeria has been studied using monthly data of
pressure, temperature and relative humidity for the period 2004-2006 for twenty six stations, obtained from the archives of
the Department of Satellite Applications Facility on Climate Monitoring (CM-SAF) DWD, Germany. The data was
retrieved from ATOVS observations onboard polar orbiting NOAA-15, 16 and 18 satellites. The results have shown that
variations in each region and at different atmospheric levels are influenced by the North–South movement of inter–tropical
discontinuity (ITD). Using the analysis of variance (ANOVA) technique on the climatological data, an empirical
relationship of the form N = ω + γNs (ω, γ constants) have been established between columnar refractivity and its surface
value, Ns, for the four regions (coastal, Guinea Savannah, midland and sub-Sahelian) in Nigeria. The initial tests carried out
on the model using data obtained over the regions for the year 2007 yielded an encouraging result as have been established
by the use of Kolmogorov-Smirnov tests.
Keywords: Inter-tropical discontinuity (ITD), Columnar refractivity, Tropospheric refractivity; Surface refractivity,
Analysis of variance (ANOVA)
obtained for the period 2004-2006 was used in Meteorological factors (pressure, P; air
computations of tropospheric refractivity for twenty temperature, T; and relative humidity, H) are obtained
six stations in Nigeria (Fig. 1), classified into four directly from the retrieved data. Water vapour
regions based on their climatic conditions19. pressure, e, on the other hand, was determined using23
As in Adeyemi13, Aro20 and Adedokun21,22, radio
refractivity between surface and 925 mb level has HeS
been labelled as low-level refractivity (NL); between e= …(6)
775 mbar and 600 mbar as mid-level refractivity 100
(Nm); and between 400 mbar and 250 mbar as upper-
where, H, is relative humidity (%); es, saturation
level refractivity (Nu).
vapour pressure (hPa). It should be noted that Nwet is
2.1 Data analysis techniques mainly responsible for the variability in N within the
The radio refractivity, N, of air for frequencies up troposphere24,25. Nwet contributes about 34% to the
to 100 GHz is generally expressed as8,13: total value of Ns while Ndry contributes about 66%.
Radio refractivity was determined using Eq. (4) for
77.6 4810 e the five different layers.
N= P + …(3) Then for each region, correlations of (NS and NL),
T T
(NS and Nm) and (NS and Nu) were obtained. High
P e correlation values were found to be predominant at
N = 77.6 + 3.75 × 105 …(4) most of the regions, most especially, at the lower and
T T2 middle levels of the atmosphere (Table 1), therefore,
using analysis of variance (ANOVA) technique (table
N dry and N wet term components from Eq. (4) are
not shown), a linear regression relation of the form
defined as:
N = ω + γN S …(7)
P e
N dry = 77.6 ; N wet = 3.75 × 105 …(5)
T T2 was developed taking into consideration the
associated F-ratios and the probability values.
where, P, is the atmospheric pressure in hPa; T, air
temperature in Kelvin; and e, water vapour pressure in 3 Result and Discussion
hPa. 3.1 Space – time distribution of tropospheric radio refractivity
In the coastal region, Figs 2(a-d), both the surface
refractivity, NS, and low-level refractivity, NL, display
double peaks with a dip in between them around
July/August. Their mean values during the dry season
are 373.3±2.96 and 327.2±2.57, respectively, while
during the rainy season they were 374.7±1.37 and
328.4±1.19, respectively. The mid-level refractivity,
Nm, on the other hand displays seemingly double
peaks with the first in May/June and the other in
September. The observed double peaks here are not as
conspicuous as in NL. Nm values are generally lower
than those of NL. This is expected to be so as water
vapour decreases with increasing height from the
surface. Its average value during the dry season is
217.5±0.95 while during the rainy season, it is
222.9±0.26. Upper level refractivity, Nu, has its values
slightly increasing from January to July, after which,
it then decreases for the rest of the year. It lies
between 109.4±0.085 and 110.3±0.095, with its mean
Fig. 1 — Map of Nigeria showing study locations value being 109.8±0.09.
304 INDIAN J RADIO & SPACE PHYS, DECEMBER 2011
Table 1 — Values of the best fit parameters ω and γ in the regression equation of NL on NS, Nm on NS and Nu on NS
(a) Low level k NL + ω + γNs
Region Ω Error in ω Γ Error in γ r CD, % F-ratio p-value
Fig. 2 — Monthly mean variations of: (a) surface refractivity, Ns; (b) low level refractivity, NL; (c) mid level refractivity, Nm; (d) upper
level refractivity, Nu in the coastal region
ADEYEMI & EMMANUEL: MONITORING TROPOSPHERIC RADIO REFRACTIVITY OVER NIGERIA 305
In the Guinea Savannah region, average NS and values of 294.5±5.19 and 291.5±3.64, respectively.
average low-level refractivity display similar pattern During the rainy months (April – October), high
with those of the coastal region [Figs 3(a-d)]. Here values of NS and NL are discernible. Their mean values
double peaks with a dip in between them in June- are 349.5±3.26 and 329.3±2.20, respectively. Nm
August are also discernible. The June-August dip is displays a single peak during May and August with
characterized by a slight increase in July. Their mean low values prevailing during the dry season and high
values during the dry season are 370.4±4.48 and values characterizing the rainy season. Mean Nm
324.7±3.89, respectively. During the rainy season, values during the dry and rainy periods are
they are 376.2±3.32 and 329.7±2.88, respectively. Nm, 211.0±0.65 and 222.4±1.13, respectively. In the case
on the other hand, display double peaks with a dip in of Nu, its values gradually decrease from January to
between them lasting between June and August. This April. The month of April seems to be where the
depression is more conspicuous than in Nm for coastal minimum for Nu occurs. A gently increasing trend
region. Its mean value is 219.7±1.09. Upper-level was then noticeable between April and July. July is
refractivity, Nu, produces a single peak in June. The the month where the annual peak occurred. Its mean
values are uniformly low, and lie between 109.4± value is 109.819 ± 0.0642 and its range lies between
0.05 and 110.2± 0.06 throughout the year. They gave 109.4±0.0643 and 110.2±0.0641 showing that in both
a mean value of 109.9± 0.065 the coastal and guinea savannah regions, Nu is
In the midland region, the structure of NS and NL uniformly low.
are well related to those of the Guinea savannah In the sub-Sahelian region, the variations observed
region [Figs 4(a-d)]. They are comparably low during in the refractivity parameters are almost similar to
the dry months (November – March), with mean those of the midland region [Fig 5(a-d)]. NS and NL
Fig. 3 — Monthly mean variations of: (a) surface refractivity, Ns; (b) low level refractivity, NL; (c) mid level refractivity, Nm; (d) upper
level refractivity, Nu at Guinea Savanna region
306 INDIAN J RADIO & SPACE PHYS, DECEMBER 2011
are low during the dry months with mean two seasons. They are the Harmattan and the rainy
refractivities of 251.7±3.47 and 269.4±2.52, season periods. During the Harmattan period, the
respectively. Their values are high during the rainy Nigeria troposphere is characterized by dry dust
season with mean refractivities of 327.2±11.57 and particles transported from Azores’s sub-tropical high
321.9±7.90, respectively. The variations observed in pressure system in the Sahara desert by the north-
Nm over this region resemble those of NS and NL. But easterly tropical continental (cT) air mass blowing
the mean refractivity values of Nm are lower than inland. This coincides with the period of no rain. The
those of NS and NL; they are 205.2±1.56 and rainy season, on the other hand, is characterized by
218.2±2.29 during the dry and rainy seasons, high humidity (heavy rainfall) brought about by the
respectively. Nu, on the other hand, has almost south-westerly tropical maritime (mT) air originating
uniform values throughout the year. However, its from the southern hemisphere. Wedged in between
value gently decreases during the late dry season the air masses, is a front known as inter-tropical
(January – April) before it rises, gently, to its highest discontinuity (ITD) (refs 26,27). The ITD migrates in
value in December. Its mean value is 109.9 ± 0.064 north-south and south-north direction reaching its
and it lies between 109.3±0.065 and 110.1±0.063. maximum northward extent 22-25°N in August and
From the foregoing, the observed variations in the its southward extent 4-6°N in January. While
refractivity parameters are attributable to the migrating, the ITD oscillates backward and forward
prevailing weather conditions pervading the country within a few latitudinal points. In January, all regions
(ref. 3). The Nigeria troposphere is characterized by north of 4°N are under the influence of the cT air.
Fig. 4 — Monthly mean variations of: (a) surface refractivity, Ns; (b) low level refractivity, NL; (c) mid level refractivity, Nm; (d) upper
level refractivity, Nu at the midland region
ADEYEMI & EMMANUEL: MONITORING TROPOSPHERIC RADIO REFRACTIVITY OVER NIGERIA 307
Conditions are such that little or no precipitation is The values of the parameters ω and γ for NL, Nm, and
experienced and refractivity parameter values are low, Nu together with the correlation coefficient, r,
and show marked decrease from the coast inland. coefficient of determination, CD and probability
During March and April, the ITD position would be (p–value) at which the null hypothesis was either
between 10°N and 15°N meaning that the regions are accepted or rejected are as shown in Table 1(a-c). The
now beginning to experience little amount of rainfall degree of association evident in Table 1 shows that NS
making variability in the refractivity parameters to be is better associated with low-level refractivity (NL)
higher. On reaching its most northern position in than those of Nm and Nu at all the regions. The degree
August, all areas in Nigeria would have been of association between NS and Nm has also been found
subjected to widespread rainfall and increased cloud to be higher at both the midland and sub-Sahelian
cover. Nm and Nu values increase and become less regions than at the coastal and Guinea Savannah
variable. NS and NL on the other hand, in the coastal regions. Regression relations obtained for all the
and guinea savannah regions, experience lower values refractivity parameters are statistically significant at
than they are in the midland and sub-sahelian regions. all the levels and at the regions except Nm in the
This is as a result of the existence of the little dry coastal region which is statistically insignificant
season that often pervades the area around the (p-value = 0.060).
West-African coast in August. To verify the reliability of the linear relations
obtained in Table 1, they were then applied to
3.2 Regression analysis evaluate refractivity parameters at the different levels
Using Eq. (7), N can be replaced by NL, Nm, or Nu and at the regions using satellite data for the year
and the corresponding parameters ω and γ evaluated. 2007. This was first done annually (January –
Fig. 5 — Monthly mean variations of: (a) surface refractivity, Ns; (b) low level refractivity, NL; (c) mid-level refractivity, Nm; (d) upper-
level refractivity, Nu at the sub-Sahelian region
308 INDIAN J RADIO & SPACE PHYS, DECEMBER 2011
December); and then for the different seasons: dry is the actual distribution function. Using observed
season (November – March); and rainy season values, Fc (xj), j = 1,2,3,……. k, are the values of
(April – October). The results of these, as compared Fc(x) evaluated at x1≤ x2 ≤ x3 ≤…..≤ xk. It has been
with the actual and calculated values, are shown in chosen to compare the Dk values with the critical
Table 2. The agreement between them is remarkable. values of Dk for the 5 percent significant level
(ref 14,28).
3.3 Kolmogorov – Smirnov (KS) test
In order to determine whether any agreement exists 1.3581
D k (α =0.05) = …(9)
between the actual distribution function Fk(x) and k1/2
each generated series Fc(x), the Kolmogorov-Smirnov
(KS) test was applied (refs 13,22,24,28). The KS test Hence, if the value of Dk obtained exceeds that of
is used to determine whether agreement exists Dk(α = 0.05), the two functions are not close enough to
between the actual and calculated values. The be related. This criterion is used to judge the degree of
maximum deviation between Fk (x) and Fc (x) is: fit or relationship, if any, of the two functions under
consideration. On applying this test to the refractivity
Dk = max Fk (x) – Fc (x) …(8) parameters at different levels and for regions, the
results obtained are shown in Tables 3 and 4.
where, x, is ordered in ascending order such that 0 ≤ Hence, as shown in Table 3, the actual and model
x1 ≤ x2 ≤ x3 ≤……..≤ xk. functions are related at all the regions at the lower and
upper levels. At the middle level and for coastal and
For k observations, guinea savanna regions, their values showed
insignificance. This may be because, as is common to
Number of values ≤ x all surface based models, perturbations due to
Fk ( x) =
k disturbances, which originate (ref. 22) from within the
Table 2 — Application of proposed model for each of the regions annually and seasonally
middle troposphere and whose effects are not quickly south movement of the intertropical discontinuity
communicated to the surface, may not be effectively (ITD). This is in consonant with the findings of
accounted for by the model. In addition, over the Balogun30 and Garbutt et al.31, who have used
hinterland of West Africa, the ITD surface is known radiosonde data obtained over Nigeria to investigate
to slope southwards with the moisture-laden south precipitation and precipitable water variations and
westerlies forming a wedge under a drier north- have concluded that differences in precipitation
easterlies29. Therefore, a station’s position relative to climatologies of the different regions in Nigeria
ITD is an important factor determining the structure account for the observed differences in the variations
of the moisture profile over the station. When a of their weather parameters; and that the variations
shallow moist layer is overlain by a deep dry layer, a observed in this weather parameters are dependent
model based on surface parameter may overestimate upon the north-south transect of the ITD. Applying
refractivity parameters aloft, and on the other hand, the analysis of variance (ANOVA) technique, a linear
where a dry layer is interspersed by a deep moist relation of the form
layer, an underestimation may result. These N = ω + γN S ( ω , γ , are constants)
limitations may then affect to a greater extent, the connecting the level refractivity parameters, N-units
effectiveness of the model. The actual and model (low-level refractivity, NL; mid-level refractivity, Nm;
functions for the mid-level refractivity during the dry and upper-level refractivity, Nu) with surface
and rainy seasons showed the same trend with it being refractivity, NS, have also been established from the
extended to the midland region (Table 4). satellite data for each region in Nigeria. The
difference observed in precipitation climatologies of
4 Conclusions the different regions has made it impossible to obtain
Using CM-SAF water vapour and temperature a single model that fits adequately the entire Nigeria
products retrieved from Advanced TIROS region. These relations when used to evaluate
Operational Vertical Sounder (ATOVS) onboard refractivity aloft obtained over each of the regions
National Oceanic and Atmospheric Administration during the year 2007 yielded an encouraging result.
(NOAA) satellites NOAA-15, 16, 18, it was possible On applying the KS test to the evaluated refractivity
to establish that seasonal variations of refractivity parameters, the actual and model functions are related
parameters over Nigeria is dependent upon the north- at almost all the levels considered at all the regions.
Table 3 — Result of KS test on actual and model results on
an annual basis Acknowledgements
The authors wish to express their gratitude to the
Region Dk (α=0.05) NL Nm Nu officials and management of the Department of
DL Dm Du Climate Monitoring SAF, Deutscher Wetterdienst,
Coastal 0.392 0.001 0.662* 0.130 Offenbach, Germany for providing the data. The
Guinea savanna 0.392 0.002 0.643* 0.034 efforts of Dr Jorg Schulz and Mr Markus Jonas of the
Midland 0.392 0.027 0.325 0.043 same Department, in making sure that the data were
Sub - Sahelian 0.392 0.069 0.292 0.048 available in ASCI format during the visit of the
*Association does not exist between actual and model values at corresponding author to CM-SAF is also greatly
5 percent significant level appreciated.
Table 4 — Result of KS test on actual and model NL, Nm and Nu values during the dry and the rainy seasons
Region Dk (α=0.05) Value of Dk for various refractivity Dk (α=0.05) Values of Dk for various refractivity
parameters during the dry season parameters during the rainy season
Dk for NL Dk for Nm Dk for Nu Dk for NL Dk for Nm Dk for Nu
*Association does not exist between actual and model values at 5 percent significant level
310 INDIAN J RADIO & SPACE PHYS, DECEMBER 2011