Evaluation of Refractivity Gradient and k-factor within the Lower
Troposphere of Maiduguri and Enugu under Two Climatic Zones in
Nigeria
FASHADE O. O.1,2, AKINWUMI S. A.1*, OMOTOSHO T. V.1, ARIJAJE T. E.1,
AYO-AKANBI O. A.1, OGUNDOLIE O. I.2, OMETAN O. O.3, ADEWUSI O. M.3
1Department of Physics,
Covenant University, Ota,
Km 10, Idiroko Road, Canaan Land, P.M.B 1023, Ota, Ogun State.
NIGERIA
2Department of Space Engineering,
African Regional Centre for Space Science and
Technology Education in English Ile-Ife, Osun State,
NIGERIA
3Department of Physics,
Lagos State University, Ojo,
NIGERIA
*Corresponding Author
Abstract: - Estimation of radio refractivity is important in the planning and design of terrestrial radio
communication links for the availability and accessibility of strong networks and signals. This paper
investigates the refractivity gradient, effective earth radius factor (k), and the geo-climatic factor K in the first
1km of the troposphere of two selected stations (Maiduguri and Enugu) under different climatic zones in
Nigeria. The indirect method of measuring radio refractivity was employed in this study to take measurements
over the two selected stations. Vertical profile values of pressure (hPa), Temperature (°C), and Relative
humidity (%) within the first 1 km were extracted from MERRA MAIMCPASM V5.20 database profile
obtained from a satellite sounding instrument by NASA in the United States. MatLab programming language
was used to evaluate the refractivity gradient, k-factor, and geo-climatic factor using the equations
recommended by ITU. The results showed that Enugu was predominantly sub-refractive due to the tropical
savannah climate while Maiduguri encountered both sub-refractive and normal refractive conditions due to the
hot semi-arid climate, and unstable and extreme weather conditions in the region.
Key-Words: - Refractivity, Meteorology, MERRA, MatLab, Troposphere, Geo-climatic and k-factor.
Received: August 19, 2022. Revised: February 14, 2023. Accepted: March 14, 2023. Published: April 5, 2023.
1 Introduction
The propagation of radio waves through the earth's
atmosphere encounters some path bending due to
diverse spatial distribution of the refractive index of
air, thereby causing severe effects such as multipath
fading, interference, and attenuation due to
diffraction. These effects significantly damage radio
communication, navigation, environmental
monitoring, and disaster forecasting systems, [1]. In
practice, the change of refractive index n is very
low, and can be difficult to have an accurate
measure of its variation, [2]. Nevertheless, this
variation impedes the transmission of radio wave.
For this reason, the refractive index n is increased
and denoted with N, [3], [4]. The anomalous
propagation that occurs during radio waves
propagation is called bending. Bending waves are
surface waves that appear in thin media, the width
of the media is small compared to the wavelength.
The refractive index N helps to characterize the
curving of the electromagnetic waves, but not the
bending phenomena.
The gradient of the radio refractive index, dN/dH is
the index used to characterize the bending of
electromagnetic waves, [5]. As electromagnetic
waves continually bend as they pass through the
atmosphere, the gradient of the radio refractivity
index can be categorised as normal, sub-refraction,
International Journal on Applied Physics and Engineering
DOI: 10.37394/232030.2023.2.1
Fashade O. O., Akinwumi S. A.,
Omotosho T. V., Arijaje T. E.,
Ayo-Akanbi O. A., Ogundolie O. I.,
Ometan O. O., Adewusi O. M.
E-ISSN: 2945-0489
1
Volume 2, 2023
super-refraction, and trapping/ducting. Figure 1.
Illustrates the propagation paths for different
refractivity conditions.
Fig. 1: Propagation paths for different refractivity
conditions, [6]
1.1 Effective Earth Radius Factor 󰇛󰇜
and Geo-climatic Factor󰇛󰇜
The effective earth radius factor is derived using
the gradient of the refractive index N as radio waves
rise through the altitude of the troposphere. It is also
used to characterize refractive conditions as normal
refraction or standard atmosphere, sub-refraction,
super-refraction, and trapping/ducting. The
variability of the k-factor is dependent on the value
endorsed by the ITU standard which is 1.33. This
value will either underrate or overrate the k-factor
value of the study locations, [7]. According to the
ITU-R-P-530 recommendation, the geo-climatic
factor is used to determine the worst month outage
probability. The estimation of this factor is based on
the point refractivity gradient󰇛󰇜, which is a
function of temperature, pressure, humidity, and
vapour pressure which enhances the seasonal
variation of radio refractive gradient, [8], [9], [10].
Many research studies had been carried out in this
field in West Africa, amongst these studies are
found in the journals, [11], [12], [13], [14], [15],
[16].
1.2 Study Area
Nigeria is a country in West Africa located between
the Equator and the Tropic of Cancer. Nigeria lies
amid latitude 4N and 14N and longitude 2E and
15E respectively with a total area of 923,768
square kilometres. The climatic condition differs in
most parts of the country, in the north the climatic
condition is arid and to the south, there is an
equatorial type of climate or tropical in nature. The
weather condition can be generally categorized into
two; In the North, from April to October is the wet
season, and from November to March is the dry
season, while in the South, from March to October
is the wet season and from November to February
dry season. The two selected areas (Enugu and
Maiduguri) for this study are located under two
different climatic zones in the country. Table 1
presents the geographical information of the two
selected stations in Nigeria.
Table 1. The Geographical Information of the Two
Selected Stations in Nigeria
2 Research Methodology
The in-situ method of measuring refractivity was
employed in this study to take measurements over
two selected stations in Nigeria namely: Enugu and
Maiduguri. Vertical profile values of pressure (hPa),
Temperature (), Geopotential height (m), and
Relative humidity () were extracted from the
measurements obtained from the Satellite Sounding
instrument by NASA in the USA via MERRA
MAIMCPASM V5.20 database profile.
A programming language was written using MatLab
to determine the radio refractivity (N) and
refractivity gradient󰇡
󰇢 applying equations (1) and
(2) recommended by [3].
 


(1)

 

󰇡

󰇢
 

 (2)
The effective earth radius factor (k) for each
selected state was computed on a Microsoft Excel
spreadsheet using the formula in equation (3):
󰇩󰇡
󰇢
󰇪 (3)
The geo-climatic factor (K) was determined using
the point refractivity gradient 󰇛󰇜in the lowest
100 m of the atmosphere not exceeded for 1% of the
average year. The 󰇛󰇜 was deduced from the
probability distribution frequency curve of an
International Journal on Applied Physics and Engineering
DOI: 10.37394/232030.2023.2.1
Fashade O. O., Akinwumi S. A.,
Omotosho T. V., Arijaje T. E.,
Ayo-Akanbi O. A., Ogundolie O. I.,
Ometan O. O., Adewusi O. M.
E-ISSN: 2945-0489
2
Volume 2, 2023
average refractive gradient of each of the two
stations over the period of years (2010 - 2014) under
consideration. The point refractivity gradient dN1
values were computed on a Microsoft Excel
spreadsheet to estimate the geo-climatic factor (K)
using equation (4):
 (4)
3 Results and Discussion
This section presents the results and discussion of
the estimation and analysis of meteorological
parameters (pressure, temperature, relative
humidity) at an altitude of 1 km of the troposphere
in two selected stations from different climatic
zones in Nigeria.
3.1 Monthly Variations of Radio Refractivity
Gradient
The refractivity gradient is strongly dependent on
the radio refractivity 󰇛󰇜 and geopotential height. It
was derived from equation (2). The results obtained
of the mean values of monthly variations of radio
refractivity gradient of the first 1 km for Enugu and
Maiduguri are presented in Table 2 and Table 3. The
result shows that Enugu had positive values both in
the wet season and dry season months and it falls
within the range of 1 N-unit/km and above. The
bending or refractive gradient that occurs within this
range is classified as sub-refractive. On the other
hand, Maiduguri had a high positive value of 54 N-
units/km observed in July (wet season) and negative
values of -3 N-units/km and -4 N-units/km observed
from November to February (dry season). The range
of these negative values is classified as normal
refractive which is within -24 N-units/km to 0 N-
unit/km. This implies that Maiduguri is
predominantly sub-refractive in the wet season and
normal refractive in the dry season. Figure 2 and
Figure 3 show the graphical bar chart of the values
of the refractivity gradient obtained in Table 2 and
Table 3.
Table 2. The mean values of monthly variations of
radio refractivity gradient of the first 1 km for
Enugu
Fig. 2: The mean values of monthly seasonal
variations of radio refractivity gradient of the first
1km for Enugu
3.2 The Result of the Effective Earth Radius
󰇛󰇜
The result of the k-factor was reliant on the type of
season. The mean values of the k-factor during dry
and wet seasons for the two selected stations in
Nigeria from (January 2010 December 2014) are
presented in Table 4. The value of 1.33
recommended by the ITU will either underrate or
overrate the variability of k-factor values obtained
in the two stations under consideration. The results
showed that Enugu had low values of k-factor
compared to the Northern station, Maiduguri which
had higher values of k-factor due to the extreme
climatic conditions that occur in this region. Figure
4 shows the graphical bar chart representation of the
results obtained in Table 4.
International Journal on Applied Physics and Engineering
DOI: 10.37394/232030.2023.2.1
Fashade O. O., Akinwumi S. A.,
Omotosho T. V., Arijaje T. E.,
Ayo-Akanbi O. A., Ogundolie O. I.,
Ometan O. O., Adewusi O. M.
E-ISSN: 2945-0489
3
Volume 2, 2023
Table 3. Mean values of monthly seasonal variations
of radio refractivity gradient of the first 1 km for
Maiduguri.
Fig. 3: Mean values of monthly seasonal variations
of radio refractivity gradient of the first 1km for
Maiduguri
Table 4. Mean effective earth radius factor of
monthly variations of radio refractivity gradient of
the first 1km for Enugu and Maiduguri from (2010-
2014).
Fig. 4: Mean effective earth radius factor of monthly
variations of radio refractivity gradient of the first
1km for Enugu and Maiduguri from (2010-2014).
3.3 The Result of the Geo-climatic Factor ()
The geo-climatic factor (K) determines the worst
month outage probability. The higher the values of
the geo-climatic factor, the more the radio wave
signals in these regions fade away. Table 5 presents
the mean point refractivity gradient 󰇛󰇜and geo-
climatic factor (K) for the period of five years (2010
- 2014) in these two selected stations. It was
observed that in Enugu and Maiduguri, the worst
months outage are found in March and June which
had geo-climatic factors K of (4.370E-05 and
3.506E-05) respectively. The geo-climatic factor
signifies the path fade depth and this also implies
that radio wave propagating signal could encounter
ducting condition in these months. The mean
monthly variations of geo-climatic factor (K) of the
first 1km for Enugu and Maiduguri from (2010-
2014) are presented in Figure 5.
-10
0
10
20
30
40
50
60
2010 2011 2012 2013 2014
Mean refractivity gradient (N-
units/km)
Year
MAIDUGURI
JAN FEB MAR APR MAY JUN
JUL AUG SEP OCT NOV DEC
0,00000
0,50000
1,00000
JAN MAR MAY JUL SEP NOV
Mean effective earth
radius factor (k)
Months
ENUGU MAIDUGURI
Month
Effective
earth
radius
factor (k)
Enugu
Average
refractivity
gradient
(2010 -
2014)
Maiduguri
Effective
earth
radius
factor (k)
Maiduguri
JAN
0.74762
2
0.98742
FEB
0.76214
7
0.95732
MAR
0.76214
13
0.92353
APR
0.76961
19
0.89205
MAY
0.78500
31
0.83511
JUN
0.80928
34
0.82199
JUL
0.83069
37
0.80928
AUG
0.82632
40
0.79695
SEP
0.80928
38
0.80513
OCT
0.82199
36
0.81347
NOV
0.78894
8
0.95152
DEC
0.78894
13
0.92353
International Journal on Applied Physics and Engineering
DOI: 10.37394/232030.2023.2.1
Fashade O. O., Akinwumi S. A.,
Omotosho T. V., Arijaje T. E.,
Ayo-Akanbi O. A., Ogundolie O. I.,
Ometan O. O., Adewusi O. M.
E-ISSN: 2945-0489
4
Volume 2, 2023
Table 5. Mean values of monthly variations of point
refractivity gradient and geo-climatic factor (K) of
the first 1km for Enugu and Maiduguri from (2010-
2014).
Fig. 5: Mean monthly variations of geo-climatic
factor (K) of the first 1km for Enugu and Maiduguri
from (2010-2014).
4 Conclusion
In this study, the evaluation of radio propagation
parameters in the first 1 km of the troposphere has
been examined in two selected stations (Enugu and
Maiduguri) under two climatic regions in Nigeria
from (January 2010 December 2014). The method
recommended by the International
Telecommunication Union (ITU) was adopted in
this study. Radio refractivity gradient, effective
earth radius (k-factor), point refractivity gradient
󰇛󰇜and Geoclimatic factor (K) was estimated and
analysed from the measurements of air temperature,
relative humidity, and atmospheric pressure
obtained from MERRA MAIMCPASM V5.20
database software using a satellite sounding
instrument by NASA, USA. In conclusion, Enugu
had positive values of refractivity gradient in both
wet and dry season throughout the period of
consideration which implies that this region was
predominantly sub-refractive, while, Maiduguri had
both positive and negative values in wet and dry
seasons respectively which indicates that this region
was predominantly sub-refractive in the wet season
and normal refractive in the dry season. The values
obtained from the refractive gradient were used to
characterize the bend in the radio wave signals in
the region if it’s either normal refractive, sub-
refractive, super-refractive, or ducting.
The application of the effective earth radius (k-
factor) in the radio link design will help to calculate
the antenna height required to estimate the
diffraction fading or multipath fading, while the
geo-climatic factor (K) will be applicable in the
calculation of fade depth A (dB). Therefore, it is
important to estimate the correct value of the radio
refractivity gradient, k-factor, and geo-climatic
factor in other to ascertain an adequate fade margin
for a reliable radio link performance.
Acknowledgements:
The authors are grateful to the University’s
Research Centre (CUCRID) of Covenant
University, Ota, Nigeria for sponsoring the
publication of this article.
References:
[1] Serdega, D. and Ivanovs, G. (2007). Refraction
seasonal variation and that influence onto GHz
range microwaves availability. Electronics and
Electrical Engineering. 6(78):39-42.
[2] M. Zilinskas, M. Tamosiunaite, M.
Tamosiuniene, E. Valma, and S. Tamosiunas
(2012). Gradient of Radio Refractivity in
Troposphere. Progress in Electromagnetics
Research Symposium Proceedings, Moscow,
Russia, pp 603-607.
[3] IUT-Radio Communication Assembly:
“Propagation Data and Prediction Methods for
4,370E-05
4,341E-05
4,370E-05
4,227E-05
4,429E-05
4,641E-05
4,831E-05
5,304E-05
4,579E-05
4,831E-05
4,519E-05
5,375E-05
7,607E-05
5,746E-05
4,370E-05
4,227E-05
4,399E-05
3,506E-05
3,529E-05
4,519E-05
4,143E-05
4,143E-05
5,708E-05
4,370E-05
Geo-climatic Factor (K)
Months
ENUGU
MAIDU…
Month
s
Point
Refractivit
y Gradient
󰇛󰇜
(2010 -
2014)
ENUGU
Geo-
climati
c factor
(K)
ENUGU
Point
Refractivity
Gradient
󰇛󰇜
(2010 -
2014)
MAIDUGUR
I
Geo-
climatic
factor (K)
MAIDUGUR
I
JAN
55
4.370E-
05
-28
7.607E-05
FEB
56
4.341E-
05
14
5.746E-05
MAR
55
4.370E-
05
55
4.370E-05
APR
60
4.227E-
05
60
4.227E-05
MAY
53
4.429E-
05
54
4.399E-05
JUN
46
4.641E-
05
88
3.506E-05
JUL
40
4.831E-
05
87
3.529E-05
AUG
26
5.304E-
05
50
4.519E-05
SEP
48
4.579E-
05
63
4.143E-05
OCT
40
4.831E-
05
63
4.143E-05
NOV
50
4.519E-
05
15
5.708E-05
DEC
24
5.375E-
05
55
4.370E-05
International Journal on Applied Physics and Engineering
DOI: 10.37394/232030.2023.2.1
Fashade O. O., Akinwumi S. A.,
Omotosho T. V., Arijaje T. E.,
Ayo-Akanbi O. A., Ogundolie O. I.,
Ometan O. O., Adewusi O. M.
E-ISSN: 2945-0489
5
Volume 2, 2023
the Design of Terrestrial Line of Sight Systems,”
Geneva, Switzerland, IUT-R P.530-15, 2013.
[4] IUT-Radio Communication Assembly: “The
Radio Refractive Index: Its Formula and
Refractivity Data,” Geneva, Switzerland, IUT-R
P.453-9, 2003.
[5] Mangum, J G. and Wallace, P. (2015).
Atmospheric Refractive Electromagnetic Wave
Bending and Propagation Delay. Publications of
the Astronomical Society of the Pacific, 127(947):
74-91
[6] Grady, R. (2014). Introduction to electromagnetic
and electro-optic propagation. Comet Program.
Retrieved on March 2, 2019 from
http://www.meted.ucar.edu
[7] Ojo, L O., Ojo, J S. and Akinyemi, P. (2017).
Characterization of secondary radio climatic
variables for microwave and millimetre wave link
design in Nigeria. Indian Journal of Radio &
Space Physics. Vol. 46: 83-90.
[8] Etokebe, I J., (2016). Determination of
refractivity gradient and geo climatic factor using
radiosonde data and inverse distance weighting
spatial interpolation for missing data.
International Journal of Systems Science and
Applied Mathematics. Volume 1, Issue 4, PP 76-
81.
[9] Bohumil Brtník, David Matoušek, Miroslav
Stehlík, Vojtěch Stejskal (2020). Comparison of
Huelsmann Basic Biquad with View to the
Decrease of the Attenuation over the Transient
Frequency of the Operational Amplifier. WSEAS
Transactions on Circuits and Systems, Volume
19, Pages: 245-249.
[10] Nadia Fezai, Abdessattar Ben Amor (2019).
Traceability Chain for the Improvement of the
Attenuation High Frequency. WSEAS
Transactions on Systems and Control. Volume
14, Pages: 213-219.
[11] Omotosho, T V., Akinwumi, S A., Usikalu, M R.,
Ometan, O O., Adewusi, M O and Abdullah, M.
(2017). Analysis of non-rainy attenuation on
earth-space path in Ota, Southwest, Nigeria. IOP
Conf. Series: Journal of Physics: Conf. Series
852 (2017) 012039 doi: 10.1088/1742-
6596/852/1/012039.
[12] Akinwunmi, S A., Omotosho, T V., Odetunmibi,
O A. (2018). Dataset of surface refractivity in
southeast, Nigeria, Data in brief 16, 470-477.
[13] Ezenugu, I A., Anthony, U I. and Colman, O A.
(2017). Estimation of Clear-Air Atmospheric
Effective Earth Radius (K-Factor) in Calabar.
American Journal of Software Engineering and
Applications, 2(3):,35-37.
[14] Adediji, A T and Ajewole, M O (2008). Vertical
profile of radio refractivity gradient in Akure
SouthWest Nigeria. Progress In
Electromagnetic Research C,Vol. 4. P.157168.
[15] Adagunodo T.A., Akinwumi S.A., Omotosho
T.V. and Akinyemi M.L. (2017). Estimation of
Specific Attenuation of Radio Signal in
Southwest Nigeria. Radio and Antenna Days of
the Indian Ocean (RADIO), September 25 28,
2017, Cape Town, South Africa, DOI:
10.23919/RADIO.2017.8242225.
http://ieeexplore.ieee.org/document/8242225/
[16] Akinwumi S. A., Omotosho T. V., Usikalu M. R.,
Ometan O. O. Study of Tropospheric Scintillation
Effect in West Africa (2018). IEEE Antennas and
Propagation Society International Symposium
and USNC/URSI National Radio Science
Meeting, APSURSI 2018 Proceedings 2018,
Article number 8608261, Pages 2487-2488.
Contribution of Individual Authors to the
Creation of a Scientific Article (Ghostwriting
Policy)
-Fashade O. O, carried out the data acquisition,
analysis and evaluation.
-Omotosho T. V supervision of research activity and
mentorship.
-Ometan O.O, Ayo-Akanbi O. A and Adewusi O. M
implemented the Algorithm used in the MatLab
programming language.
-Arijaje T. E and Ogundolie O. I was responsible for
the Statistics.
-Akinwumi S. A responsible for preparation of the
work for publication from the original research
group.
Sources of Funding for Research Presented in a
Scientific Article or Scientific Article Itself
The authors are grateful to the University’s
Research Centre (CUCRID) of Covenant
University, Ota, Nigeria for sponsoring the
publication of this article.
Conflict of Interest
The authors have no conflict of interest to declare
that is relevant to the content of this article.
Creative Commons Attribution License 4.0
(Attribution 4.0 International, CC BY 4.0)
This article is published under the terms of the
Creative Commons Attribution License 4.0
https://creativecommons.org/licenses/by/4.0/deed.en
_US
International Journal on Applied Physics and Engineering
DOI: 10.37394/232030.2023.2.1
Fashade O. O., Akinwumi S. A.,
Omotosho T. V., Arijaje T. E.,
Ayo-Akanbi O. A., Ogundolie O. I.,
Ometan O. O., Adewusi O. M.
E-ISSN: 2945-0489
6
Volume 2, 2023