The atmospheric boundary layer (ABL) and, more essentially, boundary layer turbulence, are indispensable elements of weather and climate  . They control fundamental fluxes of motion and heat amid the atmosphere and earth surfaces. The atmospheric boundary layer which begins from the surface layer up to 100 m generates the convective activities that influence circulations within the troposphere and beyond. The boundary layer extends beyond 1km during the day and shrinks to below 500 m during the periods of the night  . The boundary layer exists in two states: laminar and turbulent boundary layers  . During laminar conditions, flow is restricted within a layer while erratic mixing across layers occurs at turbulent periods. Turbulence is responsible for the contraction and expansion of the boundary layer height  . This variation of height in relation to time is subject to the degree of mixing created by the earth surface  . The stronger the mixing, the higher the boundary layer height while the weaker the mixing the thinner the layer. The boundary layer obtains considerably its heat and water through the process of turbulence and is predominantly enriched by both shear and thermal forces   . The degree of turbulence within the boundary layers shows the extent of stability situations prevalent at any location. Laminar conditions mean stable situations while turbulent conditions indicate the existence of unstable conditions. However, the amount of turbulence determines the prevailing levels of stability conditions. Studies have been conducted on atmospheric turbulence over cities and areas  -  . The level of turbulence in the boundary layer varies from one location to another due to geographical dissimilarities. This dissimilarity across locations impacts on the dispersive potentials of any local environments. The effectiveness of emission dispersion hinges on two vital components: wind speed and the strength of turbulence. While average wind speeds are reliant on local flow peculiarities, turbulence intensity relates to the generation of instability either mechanically, thermally or both  .
Boundary layer studies are very important due to the following reasons: creatures cohabit there, energy fluxes mediate creating transfers of air masses and emission dispersion is prevalent. Also, the state of any local weather is predicted from boundary layer processes and clouds within the boundary layer regulates local weather pattern. This study intends to assess the levels of turbulence within the boundary layer in Maiduguri across three vertical layers i.e. 10 - 50 m, 50 - 100 m and 100 - 1300 m. The Gradient Richardson (Rig) number method will be used to assess the levels of turbulence.
2. Study Area
Maiduguri is located at the tip of north-eastern Nigeria and roughly positioned between Latitude 11˚47'N to 11˚54'N and Longitude 13˚05'E to 13˚12'E (Figure 1). It has a climate pattern classified as Bsh i.e. semi-arid, hot and steppe by Koppen and Gigger  . Due to its closeness to the Sahel region, the effect of the hot and dry tropical continental air mass from across the Sahara desert is felt
Figure 1. Map showing study area.
throughout the year. The area has two different climatic periods, which are: the hot dry spell with quite high temperatures and the wet season. The wet season is preceded by the very cold, dry and dusty Harmattan winds  from across the Sahara desert. During this period, temperature could drop below 16˚C. Mean temperature for the area is roughly 32˚C and could extend up to 45˚C - 48˚C  between the months of March to May. Mean daily sunshine hours range between 6 - 10 hours with maximum periods between October-May and minimum periods in July and August  .
Rainfall amount over the years have ranged between 265.5 - 925.7 mm with mean value of 580.5 mm  . This low rainfall amount hinges on the location’s constraint in accessing the buoyant warm and moist tropical maritime air mass from the Ocean. However, the area receives much of this air mass that brings rainfall during July-August periods when the ITD is overhead the northern extreme of Nigeria  . It has noted  that mean relative humidity for the area falls around 25% in January to less than 60% in August. Due to the low amount of relative humidity, cloud cover is very minimal at the area throughout the year. The mean wind velocity pattern as shown on Figures 2-4 reveals range of 2 - 6 m/s for the indicated hours i.e. 0000 and 1200 hours from December-February (DJF) as well as June-August (JJA). The wind direction pattern is north-easterly in DJF and south-westerly in JJA (Figures 2-4).
Figure 2. Wind velocity/direction pattern at 0000Hr in Maiduguri from December-February.
Figure 3. Wind velocity/direction pattern at 1200Hr in Maiduguri from December-February.
Figure 4. Wind velocity/direction pattern at 0000Hr in Maiduguri from June- August.
3. Data and Methodology
The data for used for this study were obtained from the European Centre for Medium Ranged Weather Forecast (ECMWF) Era-Interim Re-analysis data for five years (2011-2015). The Era-Interim Reanalysis Platform is the up-to-date widespread atmospheric reanalysis data and has been found resourceful in the analysis of atmospheric circulation of tropical Africa. The use of the improved Era-Interim reanalysis data has exceeded prospects and expresses positive views about the achievements in the analysis of weather data realised within the last decade. The data for the indicated period of time acquired at 0.125 degree resolution was retrieved 6-hourly at 0000, 0600, 1200 and 1800 synoptic hours. Tem- perature and wind speed data were obtained at pressure level of 1000 mbar which is an approximation of the surface level data.
Many different models of assessing boundary layer turbulence exist. Researchers have generally combined one destabilizing force with one stabilizing force  and expressed these features as dimensionless ratio. Examples of these models Reynolds number, Richardson number, Rayleigh number, Ross by number and Froude number. This study utilizes the gradient Richardson number to assess the level of turbulence across three layers (10 - 50 m, 50 - 100 m and 100 - 1300 m). An approximate height where the destabilising force equals the stabilizing force through Richardson’s number relations to Monin-Obukhov length (L) was estimated. Wind velocity and potential temperature were determined at 50 m, 100 m and 1300 m altitudinal levels. MATLAB software was utilised in analysing the mathematical equations.
3.2.1. The Application of Gradient Richardson Number (Rig)
When meteorologist mention the convectional Richardson number what is being referred to is the gradient Richardson number  . The flux Richardson number is an approximation of the gradient Richardson number. The gradient Richardson number (Rig) is a dimensionless ration of the restoring force to the square of the generated force  . It is given by:
where S (a stability parameter) is the restoring force. The stability parameter is of the equation below.
Therefore, the gradient Richardson number is given by:
g is the gravitational acceleration (m/s2)
T is the mean temperature (˚C)
θ is the potential temperature
(du/dz)2, is the mean wind speed (m/s2)
Z is the vertical height (m).
The mean wind speed can be resolved into
where Zm is the mean vertical height considered. The relationship between Monin-Obukhov Length and Richardson number is given by:
The Equation (5) was used to estimate the approximate vertical distance where the restoring force equals the generated force.
3.2.2. Estimation of Wind Velocity
The analysed wind component of the gradient Richardson formula at 50 m, 100 m and 1300 m was calculated using Equation (6).
where V0 is the referenced surface wind velocity and V is the estimated wind velocity at the specified vertical height (H). The assumed roughness length for Maiduguri was 0.03 i.e. for open flat terrain with scattered settlements  .
3.2.3. Estimation of Atmospheric Pressure and Potential Temperature
The estimated atmospheric pressure at 50 m, 100 m and 1300 m levels was determined with the following equation:
where, P is the atmospheric pressure in bars, h, the height in (km), P0, pressure at height; h = 0 (P0 = 1 bar) and h0 = 7 (an approximate scale height for the atmosphere). The correspondent potential temperature (θ) at the indicated altitudes was calculated with the following equation  .
where, “Tz” is the temperature (K) of the air parcel at reference height (z), “R” is the gas constant of air and “Cp” is the specific heat capacity at constant pressure. The ratio (R/Cp) is given as (0.286) for air. At any level, z, there is a temperature (Tz) and a corresponding potential temperature (θz).
3.2.4. The Richardson Number Level of Turbulence
The Richardson number (Ri) is very important because it can identify the onset and cessation of turbulence. Theoretical and laboratory studies recommend that lamina flow turn out to be turbulent when Ri is smaller than the critical Richardson (Ric) number of 0.25  . Another term is the Richardson Termination (RT) level of 1. This indicates the termination of turbulence. Turbulent flow becomes laminar when Ri > RT. It is noted that Ric is between 0.21 - 0.25. The critical Richardson number only relates to local gradients. The thinner the layer, the closer the Ri will likely be to 0.25  .
It has been disclosed  that as the gradient Richardson number (Rig) approaches 0.25, vertical fluxes become weak. However,  have acknowledged that turbulence could persist for Ri values greater than 1, due to numbers of experimental studies conducted. An upsurge in turbulence and a decline in wind shear is observed when Ri number turns negative, conforming to unstable situations during daytime. Throughout night-time, for positive Ri values, the strength of turbulence reduces and wind shear is high  .
4. Results and Discussion
The results from analysis of the mean gradient Richardson (Rig) number across the layers for the respective hours are shown on Table 1. Findings show that the first layer (10 - 50 m) across the hours is in a turbulent state as all Rig values were
Table 1. Gradient Richardson number values in Maiduguri for the specified hours.
below the Richardson Termination point of 1. However, over 96% of the Rig values in this layer were below the Richardson critical level of 0.25. The Rig values across the second layer (50 - 100 m) have indicated an entirely lamina condition except at 1200H in February when mean Rig value was below the RiT of 1 i.e. 0.88.
Turbulence conditions within the atmospheric boundary layers differs beginning from the surface layer to the outer layer  . It has been noted that the surface layer (up to 50 m) is characterised by continuous intense  and small scale turbulence created by wind shear and thermal activities. The next layer is the less-turbulent layer (laminar boundary layer) in close contact with the earth surface. This layer establishes a buffer between the surface layer and the free atmosphere  . It has been disclosed that top the stable boundary limit in the planetary boundary layer is lower than 400 m  .
Results from the third layer (100 - 1300 m) with more vertical stretch than the other layers have indicated more turbulent conditions across the hours. All Rig values were largely negative. This indicates strongly unstable situations over the layer. It has been emphasised that largely negative Ri values indicates turbulent or unstable situations  and this is due to massive free ascents of air mass compared to very minimal forced ascents.
Results from Table 1 also highlighted the estimated height where free ascent begins to displace forced accent within the study area. The mean maximum estimated height of displacement recorded was 134.42 m at 1200H in February while the lowest was 17.89 m at 1200H in September (Table 1). This shows that thermal ascents of air parcel in the study area are prevalent at altitudinal distances below 140 m.
A correlation of surface layer wind shear values based on Richardson number  with atmospheric stability categories for Inland and Coastal sites in Northeast region of Brazil are shown on Table 2.
Furthermore, Rig values from this study at the first layer (10 - 50 m) compare well with the wind shear outcomes on Table 2. This further confirms the turbulent nature of surface layer in Maiduguri which when related with in-land site values in northeast Brazil on Table 2 shows that the surface layer is unstable for most of the time especially during the 0000H and 0600H. At the 0000H in Maiduguri, there exists terrestrial radiation from the earth surface which is being modified by surface wind speed thereby generating unstable condition at the surface layer. However, the layer above it will be laminar or stable due to the warm air overriding the cool surface.
Additionally, the Rig results for the study area indicated turbulent patterns in two dimensional grade: the first at the surface layer (10 - 50 m) and at the upper layer (100 - 1300 m) during 0000, 0600, 1200 and 1800 hours. Results for the surface layer were minor positives while that at the upper layer were largely negative. This implies that turbulence due to wind shear dominated the surface layer while that resulting from thermal or free ascent prevailed at the upper layer considered in this study. However, turbulence due to free convention is stronger than that due to wind shear. The mid layer was in a state of laminar condition throughout the period considered. Wind velocity pattern at the surface shows a mean range of 2 - 6 m/s for specified night and day hours (Figures 2-5) and this
Table 2. Atmospheric stability on wind parameters in northeast Brazil.
Source: Vieira and Sampaio, 2012.
Figure 5. Wind velocity/direction pattern at 1200Hr in Maiduguri from June-August.
is capable of modifying thermals that may arise at night and during the day. Therefore, turbulence could be regulated by the intervention of varying wind velocity across layers. Massive sensible heat exchange takes place at tropical desert where above 60 kly per annum is transferred to the atmosphere  . The study area which is of close proximity to the Saharan desert will be almost similar in trend. This sensible heat is moderated by wind shear at the surface which results in turbulent events. The stable mid layer couple with the re-radiated terrestrial radiation by dominant high clouds in the study area could have caused the unstable nature of the upper layer due to increased air temperature  . It has been admitted  that across the Saharan desert, warming from surface sensible heat flux dominates the surface layer while radiative cooling is foremost upward from the planetary layer. This radiative cooling is extremely linked with latent heat releases which generate penetrating heat almost throughout the troposphere. This process makes the upper layer turbulent and in a state of unstable conditions.
4.1. Mean Monthly Trend Analysis of Turbulence Pattern in Study Area
The monthly trends of Rig for the various layers in relation to the specified synoptic hours are displayed on Figures 6-8. During the 0000H on Figure 6, the months of April and May are more turbulent due to forced convention at the surface layer (10 - 50 m) while September and October (Figure 7) are most laminar (i.e. very stable) at the mid layer (50 - 100 m). At the upper layer (100 -
Figure 6. The monthly Rig trend at the 10 - 50 m layer for the specified hours.
Figure 7. The monthly Rig trend at the 50 - 100 m layer for the specified hours.
Figure 8. The monthly Rig trend at the 100 - 1300 m layer for the specified hours.
1300 m) during the same hour, September and October reveals the most turbulent periods due to free convection (Figure 8). This coincides with the peak/ cessation of rainy season where the moist and warm tropical maritime air mass from the Atlantic Ocean is at its maximum. The ITD reaches the Sahel region between the months June-July  and starts retreating by September across study area. Also the easterly wind which undercut the tropical maritime air mass during this time is higher in strength leading to very turbulent situations.
During 0600H period, June experiences slightly higher turbulence than the rest months at the 10 - 50 m layer while laminar conditions is stronger during the months of September to November at the 50 - 100 m layer (Figure 6 and Figure 7). At the upper layer (100 - 1300 m) during the hour, turbulent conditions due to free convection are least in June and stronger in October (Figure 8).
The turbulence pattern at 1200H on Figure 6 shows that at the surface layer turbulent is least in the months of May and September and more turbulent during the months of November to February. However, strong laminar conditions prevail between May-June as well as September-October respectively.
While September maintains high turbulence pattern due to free ascent at the upper layer (100 - 1300 m), the months of November-March is least in the fray (Figure 8). At 1800H, the months of December and January show stronger turbulent situations due to wind shear at the surface layer (Figure 6). However, at the mid layer (50 - 100 m), laminar conditions were stronger from August to October (Figure 7). At the upper layer in Figure 8, September retained the strongest turbulence due to free convection with least months from December to March.
Result analysis from the expression of relationship between Richardson number and Monin-Obukhov Length which estimates the mean height where free ascent of air mass starts dominating forced ascents is shown on Figure 9. It is revealed that free ascents destabilises forced ascents below 140 m during 1200H period in February and below 20 m during 1800H period in September. This
Figure 9. The 6-hourly monthly variation of vertical heights where free ascent destabilises forced ascent.
shows that the months of February and September at 1200 and 1800 hours are the periods of maximum and minimum heights were free convection destabilises forced convection in the study area. It has been disclosed  that the stable boundary layer is generated by wind shear and thereafter destabilised by negative buoyancy and viscosity.
4.2. Effect of Turbulence Pattern on Boundary Layer Activities in the Study Area
Rig results indicate that the strength of turbulence within the boundary layer of the study location changes amid layers. The application of Rig as a turbulent indicator also reveals the pattern of atmospheric stability conditions prevalent in the study area. This will be critical to the events that occur at the planetary layer such as emissions from various anthropogenic and natural sources, micrometeorological variations as well as the level of comfort experienced by inhabitants of the boundary layer most especially at the surface layer. The most critical factor in air pollution analysis is the level of turbulence within the boundary layer as it has intense effect on the dispersion of ground level emissions. Result analysis show that pollutant dispersions will be enhanced at 1800H from February to November due to the lower level (below 50 m) at which free ascent destabilises forced ascent (Figure 9). This will allow for convective overturning that disperses pollutants away from ground surface even at moderate wind velocity. Pollutant dispersions will be worse during December to February at 1200H; June to July at 0600H; December to January at 1800H and June at 0000H. This is due to the higher altitude at which free ascent destabilises forced ascent (Figure 9). During inversion periods at night when inversion level might be at the surface layer below an emission source of more than 50 m, pollutants will be trapped at the mid layer. However, the degree of turbulence at the upper layer of 100 - 1300 m may create windows for the escape of trapped emissions from the mid layer (50 - 100 m). When emission sources are below inversion level at the surface layer, pollutant dispersions will be concentrated within the surface layer and distributed effectively by moderate to strong wind velocity. At low wind velocity, emissions will be stagnated hence impacting receptors downwind of the source. The inversion depth during the night is governed by turbulent mixing, controlled strongly by dynamic stability  . The unstable nature at the upper layer (100 - 1300 m) will be effective in transporting emissions across boundaries as dispersions will take place above the mid layer (50 - 100 m) that is stable and slow to vertical movements. Energy transfer within laminar layers is by conduction and vertical exchanges of heat and moisture are very slow  . The turbulent pattern at the upper layer could also engender violent storms like squall line occasionally. This will be enhanced when the easterly winds that intermittently come in contact with either the moist tropical maritime air mass or the high- pressure subsidence air mass that characterises the study area from across the Saharan desert. The dusty Hamattan wind that is prevalent at the area especially during the dry season will impact receptors when transported both at the surface and upper layer. This will have adverse effects on the health of habitats and concentrations will be significant both below and aloft at close and far distances. It has been noted that the dust emissions  in the desert region are related to the occurrence of high wind velocities at the surface. Additionally, in the aviation sector, the Richardson number is utilised as a rough measure of expected air turbulence. The Rig values in the range 0.1 to 10 are typical  . Values below unity indicate significant turbulence. The Rig values obtained in this study suggest that the aviation industry in the study area must continue to monitor the levels of turbulence in the area for efficient air traffic.
Turbulence activities in the atmospheric boundary layer affect both weather processes and other vital activities that take place within the layer. Maiduguri, an urban centre located in northeastern Nigeria responds to such realities. Using the Era-Interim Re-analysis data (2011-2015), examining turbulence pattern across the boundary layers in the area have shown that the surface layer (10 - 50 m) is always in a turbulent state. The analysed gradient Richardson number (Rig) across the synoptic hours: 0000H, 0600H, 1200H and 1800H show that over 95% of values were small positive and below the Richardson critical (Ric) value of 0.25. However, all Rig values at the surface layer were below the Richardson Termination (RT) value of 1. Studies has revealed that turbulence exist due to more of wind shear when Rig values are lesser than RT. On the other hand, Rig findings shows that the mid layer (50 - 100 m) indicates laminar situation as a results of Rig values greater than 1. It was a different pattern at the upper level considered in this study (100 - 1300 m) as Rig values were largely negative indicating a strongly turbulent layer. In the study area, turbulence due to forced and free convention was prevalent. While the former was more prevalent at the surface layer, the later was prevalent at the upper layer. Turbulence due to forced ascent results from surface wind shear. That of free ascent could be as a result of the subsidence of high-pressure stable air mass that descends to the surface from across the Sahara desert. This creates unstable conditions as the air mass gains heat and assumes horizontal flow to low-pressure area. The turbulence pattern at the area suggest that air pollution will be transported to longer distance at the upper layer due to the laminar situations at the mid layer. Also, air pollution will be dispersed within the surface layer if the emission source stack is below 50 m both at low and high wind speeds. However, if the emission source stack is above 50 m, emission dispersion will take place aloft. Results also suggest that wind turbines development for power generation will be appropriate within the surface layer when compared to Ri wind shear values in northeast Brazil.
The authors show profound appreciation to the anonymous reviewers whose observations have improved the content of this study. Also the authors appreciate this Journal platform for the opportunity given to publish this study.
 Mauritsen, T. (2007) On the Arctic Boundary Layer. PhD Thesis at Stockholm University, Stockholm University, Sweden.
 Capman, N. (2011) Turbulence Intensity in Complex Environments and Its Influence on Small Wind Turbines. Department of Earth Sciences, Uppsala University, Uppsala.
 Weiss, A. (2002) Determination of Thermal Stratification and Turbulence of the Atmospheric Surface Layer over Various Types of Terrain by Optical Scintillometry. A Published Edited Version of Dissertation Submitted to the SWISS Federal Institute of Technology, Zurich.
 Friedrick, K., Lundquist, J.K., Aitken, M., Kalina, E.A. and Mashall, R.F. (2012) Stability and Turbulence in the Atmospheric Boundary Layer: A Comparison of Remote Sensing and Tower Observations. Geophysical Research Letters, 39, 1-6.
 Weli, V.E. and Emenike, G.C. (2016) Turbulent Weather Events and Aircraft Operations: Implications for Aviation Safety at the Port Harcourt International Airport, Nigeria. International Journal of Weather, Climate Change and Conservation Research, 2, 11-21.
 Peel, M.C., Finlayson, B.L. and McMahon, T.A. (2007) Updated World Map of the Koppen-Geiger Climate Classification. Hydrology and Earth System Sciences, 11, 1633-1644.
 Ogugbuaja, V.O. and Barsisa, L.Z. (2001) Atmospheric Pollution in North-East Nigeria: Measurement and Analysis of Suspended Particulate Matter. Bulletin of the Chemical Society of Ethiopia, 15, 109-117.
 Lawal, J.R., Bello, A.M., Balami, S.Y., Dauda, J., Malgwi, K.D., Ezema, K.U., Kazim, M. and Biu, A.A. (2016) Prevalence of Haemoparasites in Village Chickens (Gallus gallusdomesticus) Slaughtered at Poultry Markets in Maiduguri, North-Eastern Nigeria. Journal of Animal Science and Veterinary Medicine, 1, 39-45.
 Obot, N.I., Chendo, M.A.C., Udo, S.O. and Ewona, I.O. (2010) Evaluation of Rainfall Trends in Nigeria for 30 years (1978-2007). International Journal of the Physical Sciences, 5, 2217-2222.
 Bailey, D.T. (2000) Meteorological Monitoring Guidance for Regulatory Modelling Applications. U.S. Environmental Protection Agency. Office of Air and Radiation; Office of Air Quality Planning and Standards.
 Dickerson, R.R. and Li, Z.Q. (2015) Physics and Chemistry of the Atmosphere I. Department of Atmospheric & Oceanic Sciences, University of Maryland. AOSC 620 Lecture 2.
 Grachev, A.A., Andreas, E.L., Fairall, C.W., Guest, P.S. and Persson, P.O.G. (2012) The Critical Richardson Number and Limits of Applicability of Local Similarity Theory in the Stable Boundary Layer. Boundary Layer Meteorology, 147, 51-82.
 Texier, O., Bezault, C., Girald, N., Houbart, J.-C. and Pham, S. (2012) Results of Integration of Atmospheric Stability in Wind Power Assessment through CFD Modelling.
 Gounou, A., Guichard, F. and Couvreux, F. (2012) Observations and Diurnal Cycles over a West African Meridional Transect: Pre-Monsoon and Full-Monsoon Season. Boundary Layer Meteorology, 144, 329-357.
 Li, W.P., Chineke, T.C., Liu, X. and Wu, G.X. (2001) Atmospheric Diabatic Heating and Summertime Circulation in Asia-Africa Area. Advances in Atmospheric Sciences, 18, 257-269.
 Karam, D.B., Flamant, C., Knippertz, P., Reitebuch, O., Pelon, J. and Chong, M. (2008) Dust Emissions over the Sahel Associated with the West African Monsoon Inter-Tropical Discontinuity Region: A Representative Case-Study. Quarterly Journal of the Royal Meteorological Society, 134, 621-634.
 Wang, S., Zheng, X. and Jiang, Q. (2012).Strongly Sheared Stratocumulus Convection: An Observationally Based Large-Eddy Simulation Study. Atmospheric Chemistry and Physics, 12, 5223-5235.