Olusegun G. Fawole 1,2, Xiaoming Cai2, Rachel T. Pinker4, A. Robert MacKenzie2,3

Department of Physics and Engineering Physics, Obafemi Awolowo University, Ile-Ife 220005, Nigeria
School of Geography, Earth and Environmental Sciences, University of Birmingham, B15 2TT, UK
Birmingham Institute of Forest Research (BIFoR), University of Birmingham, B15 2TT, UK
Department of Atmospheric and Oceanic Science, University of Maryland, College Park, MD 20742, USA


Received: December 26, 2017
Revised: July 5, 2018
Accepted: July 31, 2018

Download Citation: ||https://doi.org/10.4209/aaqr.2017.12.0600  

  • Download: PDF

Cite this article:

Fawole, O.G., Cai, X., Pinker, R.T. and MacKenzie, A.R. (2019). Analysis of Radiative Properties and Direct Radiative Forcing Estimates of Dominant Aerosol Clusters over an Urban-Desert Region in West Africa. Aerosol Air Qual. Res. 19: 38-48. https://doi.org/10.4209/aaqr.2017.12.0600


  • In-depth analysis of radiative properties of aerosols in an urban-desert region.
  • Estimates of the radiative forcing abilities of dominant aerosol classes.
  • Comparison of the radiative forcing abilities of the dominant aerosol classes.
  • Estimate of the DRF and forcing ability of aerosols of gas flaring origin.


The strategic location of the AERONET site in Ilorin, Nigeria, makes it possible to obtain information on several aerosol types and their radiative effects. The strong reversal of wind direction occasioned by the movement of the ITCZ during the West Africa Monsoon (WAM) plays a major role in the variability of aerosol nature at this site, which is confirmed by aerosol optical depth (AOD) (675 nm) and Ångström exponent (AE) (440–870 nm) values with 1st and 99th percentile values of 0.08 and 2.16, and 0.11 and 1.47, respectively. The direct radiative forcing (DRF) and radiative forcing efficiency (RFE) of aerosol, as retrieved from the AERONET sun-photometer measurements, are estimated using radiative transfer calculations for the periods of 2005–2009 and 2011–2015. The DRF and RFE of the dominant aerosol classes—desert dust (DD), biomass burning (BB), urban (UB) and gas flaring (GF)—have been estimated. The median (± standard deviation) values of the DRF at the top of the atmosphere (TOA) for the DD, BB, UB and GF aerosol classes are –27.5 ± 13.2 Wm–2, –27.1 ± 8.3 Wm–2, –11.5 ± 13.2 Wm–2 and –9.6 ± 8.0 Wm–2, respectively, while those of the RFE are –26.2 ± 4.1 Wm–2 δ–1, –35.2 ± 4.6 Wm–2 δ–1, –31.0 ± 8.4 Wm–2 δ–1 and –37.0 ± 10.3 Wm–2 δ–1, respectively. Arguably due to its high SSA and assymetric values, the DD aerosol class shows the largest DRF but the smallest RFE. Its smallest AOD notwithstanding, the GF class can cause greater perturbation of the earth-atmosphere system in the sub-region both directly and indirectly, possibly due to the presence of black carbon and other co-emitted aerosol and the ageing of the GF aerosols. This study presents the first estimate of DRF for aerosols of gas flaring origin and shows that its radiative potential can be similar in magnitude to that of biomass burning and urban aerosol in West Africa.

Keywords: Gas flaring; West African Monsoon; Direct radiative forcing; Radiative forcing efficiency; Assymetric parameter.


 Atmospheric aerosols perturb the earth’s radiative energy balance both indirectly and directly on regional and global scales (Charlson et al., 1992; Haywood and Shine, 1995; Rana et al., 2009). The ability of the aerosols to alter the amount of radiation depends on their concentration, composition, and particle size distribution (Verma et al., 2017). All of these determining factors vary significantly with aerosol sources. Increased concentrations of anthropogenic aerosols in the atmosphere since the pre-industrial times has been suggested to be partly responsible for the onset of global warming (IPCC, 2013).

When perturbation of the radiative budget of the earth-atmosphere system results from the scattering and absorption of incoming solar radiation by atmospheric aerosol, the resulting radiative forcing is termed Direct Radiative Forcing (DRF). When atmospheric aerosols absorb radiation, they eventually dissipate such radiation, thereby altering the microphysical properties and lifetime of clouds, which invariably affect precipitation. Forcing resulting from such alterations is termed Indirect Radiative Forcing (IPCC, 2013). The contribution of aerosol to the total forcing due to well-mixed greenhouse gases is still associated with large uncertainties (Myhre, 2013).

The West African climate has a unique weather pattern due to the West African Monsoon (WAM) which is characterised by large-scale seasonal reversals of wind regimes (Sultan and Janicot, 2000; Barry and Chorley, 2009). The movement of the Intertropical Convergence Zone (ITCZ) and Intertropical front (ITF) are responsible for the seasonal reversal of the prevailing wind pattern in the region. Deep convection occurs in organised systems referred to as Mesoscale Convective System (MCS) (Mathon and Laurent, 2001). MCSs associated with the ITCZ can lead to rapid uplift and large scale redistribution of aerosols (Reeves et al., 2010).

In recent years, in the West African region, anthropogenic emissions of aerosols and gaseous pollutants have increased substantially, largely due to increasing population and industrialisation, a trend expected to continue until 2030 (Liousse et al., 2014). Dominant anthropogenic sources of aerosol at the study site are fossil fuel combustion, vehicular emission, biomass burning, and industrial emission while the dominant natural aerosol source is desert dust. Despite growing evidence in support of the impacts of anthropogenic aerosols on regional radiative budget, strict regulations on emissions are still not available in major African cities and, where they are available, they are very weak (Liousse et al., 2012).

In Nigeria, the various wind patterns and seasons are associated with different dominant aerosol types. While the north-easterly Harmattan (NEH) wind, pre-dominant in the dry season (November–February), brings desert dust and biomass burning aerosols into the region, the South-westerly Monsoon (SWM) wind, associated with the onset of the WAM (April–October), brings predominantly urban and industrial aerosol which are believed to contain more carbonaceous aerosols (Knippertz et al., 2015). The properties and concentrations of these aerosol types vary significantly with the wind pattern. Studies of atmospheric aerosol and their radiative effects are very scarce in Nigeria.

In this study, the radiative properties of key aerosol types at this urban-desert station were analysed and their DRF and RFE at the TOA was estimated to provide the first estimate of climate forcing of gas-flaring aerosols in the region. 


Description of the AERONET Site and Prevailing Climatic Condition

The Ilorin AERONET site (8.32°N, 4.34°E) is located at a site between the densely populated monsoonal forest region of the south and Sahel Savannah region of the north. There is pronounced variation in the climatic conditions of the region governed by the movement of the Intertropical Convergence Zone (ITCZ) and Intertropical Front (ITF), which are responsible for the seasonal reversal of the wind direction (the West African Monsoon (WAM)).

The WAM is a coupled atmosphere-ocean-land system which is characterised by summer rainfall and winter drought (Lafore et al., 2010). The rainfall in the West African sub-region results essentially from the northward movement of the low-level monsoon airflow from March to August and the southward retreat from September to November. At their northern-most position, the humid monsoonal wind from the south meet drier and warmer air to form the ITF (Cornforth, 2012). During the dry season (November–February), the West African sub-region experiences strong emissions of pollutants resulting from extensive biomass burning of vegetation often from land preparation for the incoming planting season. During the wet season (May to October), the region is strongly influenced by mesoscale convective systems (MCSs), which affects the compositions of the atmosphere through several ways including rapid vertical transport of aerosols to the upper troposphere (Law et al., 2010; Mari et al., 2011).

Gas flaring is a prominent and persistent source of atmospheric aerosols which includes soot (black carbon), SO2, CO, NOx (NO + NO2), PAH and VOCs, especially in the oil-rich regions of the world (Fawole et al., 2016a). There are over 300 active flare sites in the region where an estimated 23.7 (44.4 metric tons of CO2 equivalent) and 15.1 (28.3 metric tons of CO2 equivalent) billion cubic meters (bcm) of natural gas was flared in 2006 and 2008, respectively (Elvidge et al., 2011; Fawole et al., 2016a). In 2012, of the 325 active flare sites identified in the Nigerian oil field, 97 (~30%) ranked among the top 1000 largest flares out of the 7467 identified globally (Elvidge et al., 2015). 

Trajectory Calculation and Classification

Seven-day (168 hours) back trajectories were calculated using the UK’s Universities Global Atmospheric Modelling Programme (UGAMP) offline trajectory model. The model is driven by six-hour ERA-Interim (European Centre for Medium-Range Weather Forecasts Interim Re-Analysis) wind analyses data. The trajectories of particles are calculated backward in time by interpolating these wind analyses to the current particle position. The position (latitude, longitude) and pressure were output every trajectory time step of 0.6 hours. The choice of 7-day back trajectory length is due to the atmospheric lifetime of between 5 and 9 days estimated for black carbon (BC) and particulate organic matter (POM), respectively (Cooke and Wilson, 1996; Cooke et al., 1997; Stier et al., 2006; Koch et al., 2009). Both BC and POM are major constituents of aerosol in the study area.

As shown in several studies, for example, Bibi et al. (2016) and Alam et al. (2016), atmospheric aerosols could be clustered using the inter-relationships between different pairs of their microphysical and optical properties. Using similar techniques, prominent aerosol classes were identified at the study site in Fawole et al. (2016b) as: Biomass burning (BB), Desert dust (DD), Urban (UB) and Gas flaring (GF) aerosols. In terms of optical and microphysical properties, these classes vary significantly (Fawole et al., 2016b); mixed classes (DD-BB, DD-UB, GF-UB and GF-DD) were also identified. Using similar clustering technique of analysing aerosol optical and microphysical properties, Bibi et al. (2016) classified aerosol in the Indo-Gangetic Plain into dust, biomass and urban/industrial aerosol classes.

In this study, the properties of the single-source dominant classes were analysed to estimate their direct radiative forcing and forcing efficiency. For details of the trajectory classification and analysis of the variation of the optical and microphysical properties of the identified aerosol classes see Fawole et al. (2016b).

AERONET Data Analysis

The absolute magnitude of aerosol radiative forcing is determined, predominantly, by the values of aerosol optical depth (AOD) and single scattering albedo (SSA), while its sign is dependent on the SSA and surface albedo. Both AOD and SSA vary significantly with the source of the aerosol (Pani et al., 2016). In this study, Version 3 Level 1.5 of AERONET data released in January 2018 to which improved cloud screening and new quality controls have been applied were used to estimate the DRF and RFE of anthropogenic and natural aerosol classes in the West African sub-region. Adequate knowledge of aerosol SSA, hemispheric backscatter fraction (b) and AOD can be used to calculate the mean TOA aerosol radiative forcing for optically thin, partially absorbing aerosol (Haywood and Shine, 1995). For sites like Ilorin, where differences in the diurnal variation of aerosol properties (extensive and intensive) could be highly pronounced, the use of monthly averages of aerosol parameters will only provide highly generalised estimates of the optical and microphysical properties of the aerosol at such a site.

One of the key properties that determine the climate forcing ability of an aerosol is the angular distribution of the light scattered by the aerosol particles (Marshall et al., 1995). The angular distribution of scattered light intensity at a specific wavelength is referred to as the phase function (P). The asymmetry parameter, g, an important intensive parameter of aerosols for estimating its climate forcing ability could be derived from P. Values of g range between –1 for entirely backscattered light to +1 for entirely forward scattered light (Andrews et al., 2006). The fraction of backscattered light is the ratio of the integral of the volume scattering function over the backward half solid angle divided by the integral of the volume scattered function over the full solid angle (Horvath et al., 2016).

For AERONET retrievals, uncertainties in the direct sun measurements are within ±0.01 for longer wavelengths greater than 440 nm and ±0.02 for shorter wavelengths less than 440 nm. AOD estimated uncertainty varies spectrally from ±0.01 to ±0.02 with the highest error in the ultraviolet wavelengths (Holben et al., 1998; Eck et al., 1999). For all sky radiance wavelengths (that is, 440, 675, 870, and 1020 nm), the uncertainty in SSA is expected to be ±0.03 based on Version 1 almucantar retrieval computations (Dubovik et al., 2000; Holben et al., 2006).

Relationship between the Asymmetry Parameter and the Backscatter Fraction

Several studies, e.g., Wiscombe and Grams (1976), Marshall et al. (1995) and Kokhanovsky and Zege (1997), have attempted to parameterise the backscatter fraction (b) in terms of the asymmetry parameter (g). Studies estimating aerosol DRF have either adopted an approximate relation between b and g or look-up tables of parameterisation of aerosol optical properties such as those of Hess et al. (1998) and D’Almeida et al. (1991).

In this study, assuming spherical particles, approximate relations given in Eq. (1) as cited in Horvath et al. (2016) and Eq. (2) according to Delene and Ogren (2002) have been used to estimate backscatter fraction (b) and average upscatter fraction, β, respectively.

Aerosol Radiative Forcing

Estimating Direct Radiative Forcing

The direct radiative forcing (DRF), ∆F, of aerosol at the top of the atmosphere (TOA) is estimated using the expression derived by Charlson et al. (1992). According to Haywood and Shine (1995) the radiative transfer equation proposed by Charlson et al. (1992) is simplified as given in Eq. (3). 

where D is the fractional day length, ω is the spectrally weighted single scattering albedo, So is the solar constant, Tat is the atmospheric transmission, Ac is the fractional cloud amount, Rs is the surface reflectance,  is the spectrally weighted backscattered fraction and  is the spectrally weighted AOD. The critical value of SSA at which the DRF shifts from positive to negative is dependent on the surface albedo and asymmetric parameter, g (Haywood and Boucher, 2000; Kassianov et al., 2007). One advantage of the analytical solution for the radiative transfer equation as stated above (Eq. (3)) over a radiative transfer model is an explicit dependence on individual parameters determining the radiative forcing (Chylek and Wong, 1995). As cloud cover (Ac) is a parameter in Eq. (3), to use the expression, the assumption is that the cloud cover is above the aerosol layer which is a typical atmospheric condition in the region considered in this study.

Schemes of wavelength-dependent aerosol parameters are time-consuming and quite complex to be incorporated into radiative forcing calculations and radiative transfer codes that can produce representative and accurate estimates of radiative forcing with one or two wavelength regions (Blanchet, 1982). In their study to examine the possibility of replacing aerosol parameters by wavelength-independent parameters and the accuracy and representativeness of such average parameters for the complete solar spectrum, Blanchet (1982) found out that results of calculations with average parameter are in close agreement with corresponding terms at a wavelength (λ) of 700 nm. Haywood (1995), using detailed radiative transfer codes, tested the representativeness of average aerosol parameter and found that results at around λ = 700 nm were quite similar to those of using the entire solar spectrum. Hence, the use of aerosol parameters at λ = 675 nm, which is the nearest to 700 nm in the range of wavelengths at which aerosol parameters are measured by AERONET sun-photometers, in our estimations.

Fractional day-length, solar constant and atmospheric transmittance are assumed to be 0.5, 1370 W m–2 and 0.76, respectively (Haywood and Shine, 1995). To estimate DRF and RFE of the different aerosol classes identified in Fawole et al. (2016b), monthly mean values of cloud amount (Ac)were obtained from the ASOS-AWOS-METAR dataset (NOAA, 1998; Yang et al., 2016) for the nearest airport (Cotonou) to the site, where sufficient amount of cloud cover data are available. The surface reflectance data used are model output for albedo simulations for Ilorin (2005–2009) (R.T. Pinker, personal communication, August 2016). Fig. 1 presents the time series for the mean monthly surface albedo for Ilorin during 2005–2009. The yearly pattern of the variation of surface albedo (reflectance) is quite similar for the available five-year model output.

Fig. 1. Time series for monthly mean surface albedo at Ilorin for the period (01/2005–12/2009).Fig. 1. Time series for monthly mean surface albedo at Ilorin for the period (01/2005–12/2009).

Radiative Forcing Efficiency

SSA and backscatter fraction of the particle can be used to calculate the TOA aerosol forcing (ΔF) per unit aerosol optical depth (AOD); this is called aerosol forcing efficiency (Sheridan et al., 2002; Kaufman et al., 2005). In this study, to compare the forcing potential of the various aerosol classes, we estimated the forcing efficiency using Eq. (4). Forcing efficiency, (ΔFδ), is the aerosol radiative forcing per unit AOD. AOD is a major extensive property of the aerosol, which determines the magnitude of its radiative forcing. Forcing efficiency depends only on the nature and composition of the aerosol rather than its amount (Sheridan and Ogren, 1999).

Virkkula et al. (2014), in their study to assess the effect of aerosol from different phases of biomass burning—flaming and smouldering—on the chemical and physical properties of airborne aerosols in the boreal forest, used the expression in Eq. (4) to estimate the radiative forcing efficiency of the biomass burning aerosols. Rizzo et al. (2013), using Eq. (4), estimated aerosol forcing efficiency over a primary forest site in Amazonia.


Climatology of Aerosol Properties

Significant variation of aerosol optical and microphysical properties in the multiyear analysis of aerosol properties at the Ilorin AERONET site is a strong indication of the varying sources of aerosols at the site. The range of values, at λ = 675 nm, for aerosol optical depth (AOD), single scattering albedo (SSA) and asymmetry parameter (g) are 0.04–3.71, 0.68–0.99, and 0.58–0.8, respectively. The variation of aerosol properties is more pronounced between the non-WAM and WAM months due to seasonal reversal of the prevailing wind direction, and hence, a change of dominant sources of aerosols.

Temporal Variability of Aerosol Optical Depth, Fine Mode Fraction and Angström Exponent

During the NEH months (November–February), the values of the aerosol optical depth (AOD675) and Angström Exponent (AE440-870) are 1.22 ± 0.17 and 0.35 ± 0.06, respectively. These are months of intense biomass burning in the West African sub-region as well as intrusion of dust from the Sahara and Sahel regions. Aerosol loading in the SWM months (April–October) are characterised by lower AOD (675 nm) and high AE (440–870 nm) with median values of 0.58 ± 0.23 and 1.02 ± 0.19, respectively (see Fig. 2). The monthly plots in Fig. 2 are obtained from average daily data of aerosol parameters. Compared to AE values for similar dust sites, the relatively high average AE value of the dust aerosol in the NEH months is probably due to contributions of biomass burning aerosol at that time of the year (Fawole et al., 2016b).

The significant seasonal pattern in the AE and fine mode fraction (FMF) of aerosol at the site, as seen in Fig. 2(b) and 2(c), is due to varying aerosol sources and/or changes in atmospheric transport. There is the dominance of fine mode aerosol fraction during the West African monsoon months when the prevailing wind is the moist SWM. The influx of urban-industrial air is expected to predominate during the WAM months, between April and October (Fawole et al., 2016b). The value of AE peaks between July and September and is lowest between February and March. The lower AE values during the peak of the dry season show the strong intrusion of dust in this region at that period of the year. During the WAM months, the peak AE values between July and September coincide with the peak values of backscatter fraction, b (see Figs. 2(b) and 3(b)). This strongly suggests an increase in fine particle fraction, which is attributable to inflow of urban and industrial emissions from the south of the AERONET site. The variation of aerosol parameter with prevailing monsoonal wind as observed in this study has been observed in similar studies over Karachi, Pakistan, during the period 2006–2008 (Bibi et al., 2017) and Ahmedabad, India (Ramachandran and Kedia, 2010).

Fig. 2. Multiyear monthly variations of (a) AOD, (b) Angstrom exponent (AE), and (c) Fine Mode Fraction (FMF) of aerosol during the year 2005–2015.Fig. 2. Multiyear monthly variations of (a) AOD, (b) Angstrom exponent (AE), and (c) Fine Mode Fraction (FMF) of aerosol during the year 2005–2015.

Temporal Variability of Single Scattering Albedo (SSA) and Backscatter Fraction

In the West African sub-region, there are significant differences in the relative amount of scattering and absorption of aerosols at different periods of the year. These differences result in the variation of SSA during the years as shown Fig. 3(a). During the SWM months, as shown in the multiyear mean monthly SSA values in Fig. 3(a), inland flow of south-westerly monsoon winds are rich in partially absorbing aerosols from the urban and industrial site including gas flaring emissions from the intense gas flaring activities in the Niger Delta region. In Fig. 3(a), the WAM months (April–August) exhibit the widest range of SSA (0.66–0.98) which could result in a wider range of DRF. This wide range of SSA could be attributable to the diverse nature of aerosol in the urban and industrial emissions from the south of Ilorin.

Fig. 3. Multiyear average monthly (a) variation of SSA (675 nm) and (b) variation of backscatter fraction during the year 2005–2015.Fig. 3.
 Multiyear average monthly (a) variation of SSA (675 nm) and (b) variation of backscatter fraction during the year 2005–2015.

The boxplots in Fig. 3(b) show a steady increase in the backscatter fraction from the lowest average values of 0.07 ± 0.01 during March (peak of the NEH months) to the highest average values of 0.1 ± 0.02 during the peak of the WAM months. The median values of the backscatter fraction of the non-WAM and WAM months correspond to asymmetric parameter (g) values of 0.72 ± 0.1 and 0.61 ± 0.1, respectively. This arguably suggests a steady increase in the concentration of fine-mode aerosol fraction during the WAM months which is attributable to increased inflow of combustion aerosols from urban and industrial emissions. Mie theory predicts a higher backscatter fraction for fine-mode spherical aerosol particles (Andrews et al., 2011). The wide range of backscatter fraction during the non-WAM months (NDJF) is due to mixture of biomass burning and intrusion of desert dust which are intense during the Harmattan haze period in the region. The mean monthly value of backscatter fraction in Fig. 3(b) shows a bi-modal distribution with peaks during the intense biomass burning season (NDJ) and the peak of the WAM months (July–August) when the ITCZ is northernmost allowing enhanced inland flow of aerosol from south of the AERONET site.

Variability of Angström Exponent and AOD for the Aerosol Classes

The median (± standard deviation) values of the optical and microphysical properties of the identified aerosol classes are presented in Table 1. Unless otherwise stated, average values of aerosol parameters are reported at 675 nm. The Angström exponents (AE) discussed for the aerosol classes were estimated using the 440 nm and 870 nm wavelength pair. The highly varying range of AOD440 and AE values, 0.07–3.87 and 0.01–1.74, respectively, strongly suggests a broad range of contributing sources to the aerosol loading at the study site. As the distribution of most of the aerosol parameters for the classes are non-Gaussian, the median values are reported with the standard deviations given in brackets.

Table 1. Summary of parameter for aerosol classes used to estimate their DRF.

Desert Dust (DD)

The DD aerosol class consists of 209 days of aerosol signals, which are predominant in the NEH months and the early days of the onset of WAM months. The major source of desert dust considered in this class classification is the Sahara and Sahel dust regions (13–18°N; 6–17°E). The median values of AOD440 and AE440-870 for this aerosol class is 1.13 (± 0.54) and 0.3 (± 0.12), respectively. These values agree well with those from studies for similar sites in the Bodélé Depression of northern Chad (Todd et al., 2007), Indo-Gangetic Plain (Bibi et al., 2016) and dust regions of China (Wang et al., 2004). The average value for AOD is highest for the desert dust class while AE is the least. With a median SSA value of 0.95 (± 0.02), this class is the least absorbing.

Urban Aerosol (UB)

Aerosol signature in the urban aerosol class is prominent in the WAM months when the south-westerly moist monsoon wind is prevalent in the region. For this class, the median AOD440 and AE values are 0.53 (± 0.35) and 0.52 (± 0.34), respectively. Even though this AE value is low, it is stillhigher than that for the DD aerosol class. The DD aerosol class is expected to contain a higher fraction of coarse aerosol. This class of aerosol (Urban), with a median value of SSA of 0.93 (± 0.04), is partially absorbing arguably due to increased carbonaceous particle content from anthropogenic sources in the urban area. 

Gas Flaring Aerosol (GF)

This class is similar to the urban class but has a lower median value of AOD440 and an average AE value, which is higher than that of the urban aerosol by a factor of ~2. For this class, the median values of AOD440 and AE values are 0.41 (± 0.26) and 1.16 (± 0.29), respectively. A median value of SSA of 0.9 (± 0.06) makes it more absorbing than the urban aerosol class, which is attributable to it having a relatively higher carbonaceous particulate content. This class is estimated to have an average Absorption Angstrom Exponent (AAE) of 0.98 (± 0.25) in contrast to urban aerosols, which has an AAE value of 1.2 (± 0.38) (Fawole et al., 2016b). The gas-flaring region, south of the AERONET site, contains more than 300 active flares (Elvidge et al., 2015), where it is estimated that more than 25% of the annual natural gas production is flared (Elvidge et al., 2009; Ite and Ibok, 2013; Anejionu et al., 2015).

It should, however, be noted that gas flaring, a prominent source of soot (BC), also emits other aerosol including volatile organic compounds (VOCs), SO2 and NOx, some of which exert a cooling effect on the climate (USEPA, 2012).

Biomass Burning

Similar to findings from the studies by Bibi et al. (2016) and Tiwari et al. (2016), the BB aerosol class is characterised by high AOD and high AE, which is typical of biomass burning sites. Although it is prevalent almost at the same time as the desert dust season during the NEH months, it can be distinguished by its lower SSA and higher AE values. For this class, the median values of AOD440 and AE are 0.93 (± 0.3) and 1.0 (± 0.25), respectively. The range of values for AOD440 and AE are in agreement with values reported by Ogunjobi et al. (2008) and Pace et al. (2006) for similar biomass burning sites in West Africa and around the central Mediterranean, respectively. The regions of biomass burning considered in the classification of this class are (i) 6.5–11.5°N; 3°W–3°E and (ii) 6.5–11.5°N; 13.6°E–22.5°E (see Fawole et al. (2016b)). The choice of these BB regions are based on data obtained from MODIS active fire detection over Africa as reported by Roberts et al. (2009). In agreement with reports from previous studies from similar sites, this class, with median SSA value of 0.86 (± 0.04), is the most absorbing class presumably due to its enhanced organic carbon (OC) content.

Table 2 presents the mean monthly surface reflectance, RS, and cloud amount, AC, used in the DRF estimation. Sources of these parameters are stated in Section 2.4.1.

Table 2. Mean monthly surface reflectance and cloud amount for study site. 

Aerosol Radiative Forcing

The direct radiative forcing (DRF) for each aerosol class at the TOA was estimated using the relationship in Eq. (3). The values of surface reflectance, Rs, and cloud amount, AC, used in the estimations range between 0.12–0.15 and 0.15–0.39, respectively. Table 1 presents, for each cluster, the range of values of AOD (τ), SSA (ω), backscatter fraction (b) used in the DRF estimations. The daily average values of these aerosol parameters were used for the DRF estimations of the four classes.

Fig. 4 presents the variation of the DRF estimates at TOA for each aerosol cluster. The DD aerosol class with median AOD675 value of 0.91 (± 0.44) has the highest DRF of –30.3 ± 13.4 Wm–2 at TOA. Of the four classes identified, this class has the highest mean AOD. The median DRF value for this class, as shown in Fig. 4, is highest. It is believed that this high DRF value is due to the large SSA (average value of 0.97 ± 0.02), which brings about less absorption and more scattering, and the largest asymmetry factor, g, of 0.74 ± 0.03, which causes more forward scattering of incoming radiation, compared to the other aerosol classes. Thus, consistent with the findings of García et al. (2012), this class has the most effective cooling effect on the earth-atmosphere system at the TOA in the region.

Fig. 4. Direct radiative forcing (DRF) at the TOA for the different aerosol classes.Fig. 4. Direct radiative forcing (DRF) at the TOA for the different aerosol classes.

The biomass burning (BB) aerosols class has an estimated DRF of –23.6 ± 8.9 Wm–2 at TOA with average AOD675 value of 0.61 (± 0.26). Compared to the DD class, this class has a relatively shorter range of DRF. This value of DRF is comparable to the mean DRF obtained by García et al.(2012) and Yoon et al. (2005) for similar biomass burning site in South America and South Africa, respectively. Compared to the DD and BB aerosol classes, the urban (UB) class exerts a smaller cooling effect. This class (urban) with average AOD675 of 0.38 (± 0.23), rich in anthropogenic urban aerosol, is estimated to have a DRF value of –11.7 (± 7.5) Wm–2. In a study by Yoon et al. (2005), similar values of DRF were obtained for US East Coast (Goddard Space Flight Center (GSFC)), a heavily populated urban area. With median DRF value of –8.2 ± 5.8 Wm–2 (AOD675 = 0.29 ± 0.21), the GF class has the least cooling effect at TOA. The GF class is believed to be rich in fossil fuel combustion emissions including sulfate and black carbon (Fawole et al., 2016a).

Radiative Forcing Efficiency (RFE)

The absolute magnitude of the DRF is dependent not only on the amount of radiation entering the atmosphere but also on the quantity of aerosol perturbing the atmosphere (Bush and Valero, 2003). For a better understanding of the impact of aerosol optical depth (AOD) on the estimation of aerosol DRF, the radiative forcing efficiency (RFE) of the different classes was estimated using Eq. (4). Since RFE is independent of AOD, it is a useful tool to compare the forcing abilities of different aerosol types. As such, the influences of other variables, such as SSA, absorption and scattering properties, and surface albedo might become more evident (García et al., 2012).

Fig. 5 shows the variation of the RFE for the different aerosol classes. The natural aerosol, desert dust (DD), has the least average RFE of –31.0 ± 3.3 Wm–2δ–1. The BB aerosol class, like the DD class, has a short range of RFE, but a higher mean RFE value of –39.0 ± 4.0 Wm–2δ–1. The UB and GF aerosol classes have average RFE values of –32.4 ± 5.4 and –36.0 ± 7.8, respectively. The BB aerosol class has the highest median RFE value. These two classes (UB and GF), compared to the DD and BB aerosol classes, have relatively wider ranges of RFE. As such, aerosols in these two classes have the ability to perturb the earth-atmosphere system more in this region.

Fig. 5. Radiative forcing efficiency (RFE) for the different classes.Fig. 5.
 Radiative forcing efficiency (RFE) for the different classes.

Fig. 6 presents the relationship between DRF and aerosol optical depth AOD675 for the different aerosol types. The slope of best-fit line, shown in red, gives the average forcing efficiency, as estimated by Eq. (4). The regression equation for the plot and correlation of DRF and AOD values are also shown in red on Fig. 6. In Fig. 6, N is the number of days clustered into each aerosol class.

Fig. 6. The relationship between DRF and AOD (675 nm) for (a) Desert dust cluster, (b) Biomass burning cluster, and (c) Urban cluster and (d) Gas flaring cluster.Fig. 6.
 The relationship between DRF and AOD (675 nm) for (a) Desert dust cluster, (b) Biomass burning cluster, and (c) Urban cluster and (d) Gas flaring cluster.


 The variations of aerosol optical and microphysical parameters—AOD, SSA, asymmetric parameter, Ångström exponent and backscatter fraction—were studied for the West African sub-region using AERONET retrievals from Ilorin, Nigeria. The DD aerosol class is characterised by high AOD and low AE. The BB aerosol class is characterised by high AOD and high AE, while the GF class is characterised by low AOD and high AE. The direct radiative forcing of the various dominant aerosol types was estimated using aerosol parameters from AERONET retrievals as inputs in a simplified radiative transfer equation proposed by Haywood and Shine (1995). Due to differences in methodologies and varying aerosol sources/nature, it is difficult to directly compare our results (average DRF values) with those from literature. Desert dust (DD) and biomass burning aerosols were found to be the most effective cooling aerosol in the region at the TOA. The UB and GF aerosol classes, which have been suggested as being rich in emissions from the combustion of fossil fuel (i.e., black carbon and sulphate), have less of a cooling effect. The more absorbing aerosols (GF and BB) show a higher forcing efficiency, and the GF aerosol class exhibits the largest variability in RFE. These results suggest the need for concerted efforts to adequately characterise and quantify emissions from real-world gas flares, as they make significant contributions to the radiativetransfer in the earth-atmosphere system, particularly in oil-rich regions, where gas flaring is persistent, continuous and substantial. To the best of our knowledge, this is the first estimate of DRF for the gas-flaring dominant aerosol class.

Findings from this study, especially as they relate to the GF cluster, suggest the need for an adequate understanding of the behaviour and transformation of atmospheric aerosol of gas flaring origin. A chemistry transport model with adequate schemes to simulate the behaviour of aerosols will be very appropriate for this proposed study.

Don't forget to share this article 


Subscribe to our Newsletter 

Aerosol and Air Quality Research has published over 2,000 peer-reviewed articles. Enter your email address to receive latest updates and research articles to your inbox every second week.

Latest coronavirus research from Aerosol and Air Quality Research

2018 Impact Factor: 2.735

5-Year Impact Factor: 2.827

SCImago Journal & Country Rank

Aerosol and Air Quality Research (AAQR) is an independently-run non-profit journal, promotes submissions of high-quality research, and strives to be one of the leading aerosol and air quality open-access journals in the world.