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