Assessing Spatio-Temporal Variability of Wildfires and their Impact on Sub-Saharan Ecosystems and Air Quality Using Multisource Remotely Sensed Data and Trend Analysis

Globally, wildfires are considered the most commonly occurring disasters, resulting from natural and anthropogenic ignition sources. Wildfires consist of burning standing biomass at erratic degrees of intensity, severity, and frequency. Consequently, wildfires generate large amounts of smoke and other toxic pollutants that have devastating impacts on ambient air quality and human health. There is, therefore, a need for a comprehensive study that characterizes land–atmosphere interactions with regard to wildfires, critical for understanding the interrelated and multidimensional impacts of wildfires. Current studies have a limited scope and a narrow focus, usually only focusing on one aspect of wildfire impacts, such as air quality without simultaneously considering the impacts on land surface changes and vice versa. In this study, we use several multisource data to determine the spatial distribution, frequency, disturbance characteristics of and variability and distribution of pollutants emitted by wildfires. The specific objectives were to (1) study the sources of wildfires and the period they are prevalent in sub-Saharan Africa over a 9 year period, i.e., 2007–2016, (2) estimate the seasonal disturbance of wildfires on various vegetation types, (3) determine the spatial distribution of black carbon (BC), carbon monoxide (CO) and smoke, and (4) determine the vertical height distribution of smoke. The results show largest burned areas in December–January–February (DJF), June–July–August (JJA) and September–October–November (SON) seasons, and reciprocal high emissions of BC, CO, and smoke, as observed by Modern-Era Retrospective Analysis for Research and Applications version 2 (MERRA-2) and Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO). In addition, the results reveal an increasing trend in the magnitude of BC, and CO concentration driven by meteorological conditions such as low precipitation, low relative humidity, and low latent heat flux. Overall, this study demonstrates the value of multisource remotely sensed data in characterising long-term wildfire patterns and associated emissions. The results in this study are critical for informing better regional fire management and air quality control strategies to preserve endangered species and habitats, promote sustainable land management, and reduce greenhouse gases (GHG) emissions.


Introduction
Wildfires are common in many parts of the world, igniting either naturally by lightning or intentionally by anthropogenic activities such as grassland management, deforestation and "slash and burn" agricultural practices. Although wildfires are an integral part of many terrestrial ecosystems, their frequency, intensity, and extent may be detrimental to endangered and threatened ecosystems and their resilience. For example, studies at a global scale show that wildfires led to approximately 400 million hectares (ha) of burned land area, causing destruction of property, biodiversity, ecosystems and associated services in their path [1,2]. Thus, information on fire frequency, intensity and burned area (BA) is critical for understanding and alleviating the environmental, atmospheric and socio-economic impacts of wildfires [3]. Specifically, BA is considered the single most important information for addressing global mandates that seek to reduce fire emissions, such as the Kyoto Protocol, the Paris Agreement, United Nations Sendai Framework on Disaster Risk Reduction 2015-2030, and the Sustainable Development Goals (SDGs), and is a critical input to atmospheric and climate change modelling [4].
Fortunately, BA detection and mapping is well established and has been studied since the dawn of satellite imagery [5] with most recent studies focusing on (a) the development and improvement of detection and mapping techniques [6][7][8], (b) enhancement of existing global products both in detection accuracy and spatial detail [9,10] and (c) the inter-comparison and validation in different environmental settings and regions [11]. Additionally, the availability of operational satellite-based products, such as land cover, temperature, rainfall, tree cover, etc., provide prospects for assessing and quantifying the impact of wildfires on the ecosystems and biodiversity. However, only a few studies [12][13][14] have exploited these products for this purpose. Using GLOBCOVER, the Vegetation Continuous Field (VCF), Moderate Resolution Imaging Spectroradiometer (MODIS) Fire Radiative Power (FRP) and the Tropical Rainfall Measuring Mission (TRMM) datasets, Cardozo, Shimabukuro, Pereira and Silva [12] found that biomass of broadleaved evergreen or semi-deciduous forest were mostly burned, attributed to grassland management practices, agricultural land clearing and the deforestation of the Amazon tropical rainforest.
Alarmingly and despite devastating ecological impacts, wildfires also emit approximately 2 Pg of carbon into the atmosphere, thus affecting the atmosphere's chemical composition, visibility, climate and air quality [15][16][17]. The burning and combustion of vegetation biomass significantly increase the global atmospheric load of trace gases and aerosols such as smoke, which impacts on air quality, the earth's radiation budget and visibility [18,19]. In particular, smoke aerosols perturb regional and global radiation budgets through their light-scattering effects and influences on cloud microphysical processes [20]. In addition, smoke aerosols have been linked to illnesses in human populations [21]. This is not surprising since smoke emissions contain mainly carbon monoxide (CO) and black carbon (BC), which can be detrimental when inhaled by exposed populations depending on particle sizes. For example, BC particle sizes that are between 0.1 µm to 0.25 µm (i.e., ultrafine to fine) can be more harmful to human health than larger particles since they tend to accumulate in the alveolar region of the lungs. On the other hand, particles larger than 4 µm normally deposit in the mouth and throat [22]. Essentially, the heights at which smoke is observed is critical as it determines the lifetime of the pollutants, their downwind transport dispersion pathways and the magnitude of their environmental and health impacts [23,24]. Val Martin et al. [25] further showed that observed altitudes of smoke vary with vegetation type, geographic location and season.
CO and BC particles can be observed at the lower parts of the stratosphere and are able to cause strong radiative forcing [26]. An estimate by Wang et al. [27] showed that the global scale BC radiative forcing can reach +0.53 W m −2 (after removing the preindustrial background) indicating a greater contribution to warming influence. As a result, BC is considered a significant climate warming agent [28]. It can increase the absorption of incoming solar radiation (primarily visible wavelengths), it can reduce the surface albedo of snow and ice, thus accelerating melting, and it can alter the number and composition of small particles on which water vapour condenses, thus affecting the lifetime, reflectivity and stability of clouds [29]. Therefore, studying the sources of BC, its spatial and vertical height distribution has become the key focus in the recent literature due to devastating climate and health implications [30][31][32].
Datasets from sensors such as the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) and reanalysis data such as the Modern-Era Retrospective analysis for Research and Applications-Version 2 (MERRA-2) have been instrumental in characterising the Spatio-temporal and vertical patterns of smoke, BC and CO. For example, using CALIPSO data, Shikwambana [33] showed that biomass burning is responsible for a substantial amount of smoke aerosols in central Africa during the June-July-August (JJA) season. On the other hand, Lasko et al. [34] used the MERRA-2 reanalysis BC concentration data from 2012-2016 to explore the pollution trends over time. Their results show that forest biomass burning from Northwest Vietnam and Laos impacts on the air quality of Hanoi during the peak Ultraviolet Aerosol Index (UVAI) months of March and April. However, these studies are limited in scope and have a specific focus on one aspect of wildfire impacts, such as air quality, without simultaneously considering the impacts on land surface changes and vice versa. The effects of wildfires are interrelated, i.e., have implications for both the surface and atmospheric changes. Therefore, a comprehensive study that characterises land-atmosphere interactions with regard to wildfires is critical in understanding the spatial distribution, frequency, disturbance characteristics, variability and distribution, of pollutants (such as smoke, BC and CO) from wildfires. This paper addresses this gap using multisource remotely sensed data, reanalysis data and trend analysis over a period of 9 years, i.e., 2007-2016. The objectives of the study are to (1) study the sources of wildfires and the period they are prevalent in sub-Saharan Africa, (2) estimate the seasonal disturbance of wildfires on various vegetation types and (3) determine the spatial distribution of black carbon, carbon monoxide and smoke, and (4) determine the vertical height distribution of smoke.

Study Area
The study area consists of the Sub-Saharan African (SSA) region covering latitudes of between 25 • N and 35 • S ( Figure 1). This area is characterised by a number of biomes that are habitats to thousands of terrestrial mammals, birds, herpetofauna and vascular plants. The area around the equator consists of mainly tropical forests, whilst subtropical areas consist of tropical and subtropical Savanna. Most of southern Africa consists of temperate Grasslands, Savanna and Shrubland, as well as Mediterranean forests. Parts of the SSA region also consists of deserts and Xeric Shrublands [35]. A detailed land cover map for 2015 is given in Figure 1. The study area was chosen because it is mostly affected by wildfires resulting from anthropogenic and natural activities and has been found to constitute about 70% of the total burned area across the world [16].

Burned Area
The Moderate Resolution Imaging Spectroradiometer (MODIS) Burned Area (BA) product (MCD64A1) version 6 contains monthly burned-area and quality information per-pixel on global scales, at a 500 m spatial resolution. The MCD64A1 was derived using the Collection 6 algorithm that identifies, at the time of satellite overpass, pixels that contain one or more actively burning fires using brightness temperatures [7,36]. Collection 6 BA algorithm is an improvement of Collection 5.1 MCD64A1 and Collection 5.1 MCD45A1 mapping algorithms. The algorithm has reportedly better detection rate for small burns, has reduced temporal uncertainty and extent of unmapped areas that the product detects [11,37]. The annual total burn area (BA) and pixel-wise frequency of fire incidence per season (i.e., DJF, MAM, JJA, SON) were derived by aggregating MCD64A1 Burn Date product in Google Earth Engine (GEE). From our knowledge, no similar study has been found yet that used MCD64A1.v6 to benefit from its improvements.

Precipitation
Tropical Rainfall Measuring Mission (TRMM) dataset is produced based on Algorithm 3B43 that merges several precipitation estimates from multiple microwave orbital sensors and intercalibrated to the TRMM Combined Instrument. The algorithm produces monthly scale precipitation rate (in mm/month) by combining the 3-hourly merged high-quality/IR estimates (3B42) with the monthly surface precipitation accumulated Global Precipitation Climatology Centre (GPCC) rain gauge analysis. The Monthly Precipitation Estimates (TRMM 3B43) dataset has a spatial resolution of 0.25° × 0.25°.

Burned Area
The Moderate Resolution Imaging Spectroradiometer (MODIS) Burned Area (BA) product (MCD64A1) version 6 contains monthly burned-area and quality information per-pixel on global scales, at a 500 m spatial resolution. The MCD64A1 was derived using the Collection 6 algorithm that identifies, at the time of satellite overpass, pixels that contain one or more actively burning fires using brightness temperatures [7,36]. Collection 6 BA algorithm is an improvement of Collection 5.1 MCD64A1 and Collection 5.1 MCD45A1 mapping algorithms. The algorithm has reportedly better detection rate for small burns, has reduced temporal uncertainty and extent of unmapped areas that the product detects [11,37]. The annual total burn area (BA) and pixel-wise frequency of fire incidence per season (i.e., DJF, MAM, JJA, SON) were derived by aggregating MCD64A1 Burn Date product in Google Earth Engine (GEE). From our knowledge, no similar study has been found yet that used MCD64A1.v6 to benefit from its improvements.

Precipitation
Tropical Rainfall Measuring Mission (TRMM) dataset is produced based on Algorithm 3B43 that merges several precipitation estimates from multiple microwave orbital sensors and intercalibrated to the TRMM Combined Instrument. The algorithm produces monthly scale precipitation rate (in mm/month) by combining the 3-hourly merged high-quality/IR estimates (3B42) with the monthly surface precipitation accumulated Global Precipitation Climatology Centre (GPCC) rain gauge analysis. The Monthly Precipitation Estimates (TRMM 3B43) dataset has a spatial resolution of 0.25 • × 0.25 • .

Land Cover
The land cover data utilised in this study, i.e., ESA Climate Change Initiative Land Cover (CCI-LC)  was generated from global daily reflectance measurements from 1992 to 2015 acquired by multiple  Earth observation systems, namely, MEdium Resolution Imaging Spectrometer (MERIS, 2003-2012), Advanced Very High Resolution Radiometer (AVHRR, 1992(AVHRR, -1999, Satellite Pour l'Observation de la Terre Vegetation (SPOT-VGT, 1999-2013) and PROBA-Vegetation (PROBA-V, 2014-2015). The land cover classification was based on both unsupervised and supervised algorithms, time-series data and a United Nations Land Cover Classification System (UN-LCSS) containing 37 classes that describe Earth's terrestrial surface [38]. The use of multi-temporal, multi-sensor data allowed backwards and forward cross-checking to detect and verify land cover changes and eliminate false changes as a result of inter-annual variability. The CCI-LC database consists of 24 annual land cover maps spanning from 1992 to 2015 at a spatial resolution of 300 m [39,40]. The CCI-LC data from 2007 to 2015 were used in this study to determine burned land cover per season. Detailed Level 2 classes were reclassified to Level 1 classes [38] and non-vegetation classes were removed.

Modern-Era Retrospective Analysis for Research and Applications Version 2 (MERRA-2)
The Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) supplies data from 1980 to the present. MERRA-2 was developed to improve and extend the original MERRA dataset [41] which finished in February 2016. MERRA uses 3-dimensional variational data assimilation (3DVAR) analysis algorithm based on the GridPoint Statistical Interpolation (GPS) scheme with a 6-h update cycle. MERRA makes comprehensive use of satellite radiance data from both operational and research instruments, including data from hyperspectral instruments. The assimilation of radiance data need a forward radiative transfer model (RTM) as the observation operator, to calculate the model-equivalent radiances, and the corresponding Jacobian to calculate the influences in model space of the radiance increments calculated from the analysis. The successful implementation of radiance data needs cautious quality control and bias correction procedures that are channel-specific. The bias in a given satellite channel can vary considerably in space and time, depending on the systematic errors in the radiative transfer model, atmospheric conditions and the quality and age of the instrument. In most data assimilation schemes, the bias in each satellite radiance measurement is represented by a linear predictor model with a relatively small number (~10) of parameters. In MERRA-2, aerosol and meteorological observations are jointly assimilated within Goddard Earth Observing System Model, Version 5 (GEOS-5) [42,43]. The system assimilated MODIS, Advanced very high-resolution radiometer (AVHRR), Multi-angle Imaging SpectroRadiometer (MISR) and AErosol RObotic NETwork (AERONET) data. More details on MERRA can be found in Rienecker et al. [41]. The datasets used in this study are the black carbon extinction 550 nm, carbon monoxide concentration and latent heat flux overland products.

CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations)
The CALIPSO satellite was launched on 28 April 2006 and is an integral part of NASA's A-Train satellite constellation [44]. Its main function is to study the role that aerosols and clouds play in regulating the earth's climate, weather and air quality. The main objectives of the CALIPSO are to (1) provide statistics on the vertical structure of clouds around the globe, (2) detect sub-visible clouds in the upper troposphere and Polar Stratospheric Clouds (PSC), and (3) provide statistics on the geographic and vertical distribution of aerosols and clouds around the globe. The primary instrument aboard CALIPSO is the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP), a two-wavelength laser (532 nm and 1064 nm) operating at a pulse repetition rate of 20.16 Hz [45]. More information on the technical specifications of CALIPSO can be found in Winker et al. [46] and Winker et al. [47]. Furthermore, Hunt, Winker, Vaughan, Powell, Lucker and Weimer [45] gives a detailed discussion on the algorithm that was developed to identify aerosol and cloud layers and to retrieve a variety of optical and microphysical properties. The data used in this study are smoke AOD 532 nm and smoke extinction coefficients 532 nm.

AIRS (Atmospheric Infrared Sounder)
The Atmospheric Infrared Sounder (AIRS) on the EOS/Aqua satellite was launched in polar orbit in May 2002. AIRS is a high spectral resolution spectrometer with 2378 IR channels between 3.74 and 4.61, 6.20-8.22 and 8.8-15.5 µm [48], which is essential for atmospheric temperature and relative humidity soundings. AIRS also has four visible and near-IR channels between 0.40 and 0.94 µm, which are mainly used for the detection of clouds in the IR FOV. The AIRS retrieval algorithm has been developed and validated gradually with a clear sky and clear/cloudy conditions over a non-frozen ocean, and then the non-polar land and polar cases [49]. Retrievals yield cloud fraction, cloud top pressure and temperature, surface temperature, and vertical profiles of temperature and water vapour. The measurements' vertical resolution and uncertainties are described by Aumann et al. [48]. In this study, relative humidity data were used.

Mann-Kendall Test
The Mann-Kendall Trend Test (MK test) [50,51] is used to analyse data collected over time for characterising in a consistent manner, increasing or decreasing trends in Y values. This is a non-parametric rank-based test method which is commonly used to identify monotonic trends in a time series [52] of climate data, environmental data and hydrological data. The MK test is known for its robustness to outliers, because its statistic is based on the sign of differences, rather than actual values of the parameter (α = 0.05).
The test is based on the statistic S defined as follows: where n is the number of data points, x j and x i are data values at time j and i( j > i), respectively.
For large samples (n > 10), the sampling distribution S is assumed to be normally distributed with zero mean and variance as follows: where n is the number of tied (zero difference between compared values) groups and t i is the number of data points in the ith tied group. The Z-statistic or standard normal deviate is then computed by using the equation: If the computed value of Z > Z α 2 , then the null hypothesis of no trend is rejected at α level of significance in a two-sided test (i.e., the trend is significant). In this study, the null hypothesis was tested at a 5% significance level.

Sequential Mann-Kendall Test
The Sequential Mann-Kendall Test (SQMK) proposed by Sneyers [53] is important for determining changes of a trend over time. This test sets up two series, a progressive u(t) and a backward (retrograde) series u (t). If they cross each other and diverge beyond the specific threshold value, then there is a statistically significant trend. The point where they cross each other indicates the approximate year at which the trend begins [54]. Whenever the progressive (u(t)) or the retrograde u (t) row exceeds certain limits, before and after the crossing point, the null hypothesis of the sampled time series has no change points and must be rejected. This trend turning point may be significant at a particular significance level (i.e., a 5% significance level). If the intersection of u(t) and u (t) occur within ±1. 96 (that corresponds to the bounds at 5% significance level) of the standardized statistic, a detectable change has occurred at that point along the temporal time series. The change can be detected by Kendall coefficient (t i ) [55] and when a time series shows a significant trend, the period from which the trend is noticeable can be obtained effectively by this test.

Seasonal Effects of Wildfires on Land Surface Dynamics
The impacts of wildfires on the environment are multifaceted and complex, varying by extent, burning frequency, distribution and types of emissions [1,4,33]. Therefore, in order to characterise the ecological impacts of wildfires, this study employed multisource remotely sensed data from MODIS, MERIS, SPOT, and Proba-V. Generally, the results (see Figure [56]. In addition, the year-to-year fluctuations observed can be explained by the highest burn frequency (see Figure 3) during these seasons and newly burned areas shown by a lower burn frequency (< 20%). During SON, only a few areas show a high burn frequency in the region  [56]. In addition, the year-to-year fluctuations observed can be explained by the highest burn frequency (see Figure 3) during these seasons and newly burned areas shown by a lower burn frequency (< 20%). During SON, only a few areas show a high burn frequency in the region between 0 • and 20 • N. On the other hand, BA during MAM is relatively the lowest across all the years, the highest is recorded in 2007, 2012 and 2013 and generally has a decreasing trend. This is in agreement with the results in Figure 3 that show the lowest burn frequency, i.e., < 20% distributed between the region 0 • to 25 • S and 5 • to 15 • N, respectively.
In general, the results of burn frequency and BA imply that most biomass burning in sub-Saharan Africa occur during the DJF season, i.e., mainly in the northern hemisphere. On the other hand, most biomass burning during JJA and SON occurs mainly in the southern hemisphere. The frequently drier conditions as a result of prevailing meteorological conditions, compounded by poor fire management and the susceptibility of certain vegetation types provided suitable conditions for widespread fires during these seasons. For example, warmer and drier conditions increase fuel load and flammability [57,58]. The higher burn frequency, i.e., > 40% during DJF, JJA and SON may have devastating implications for the sustainability of natural forests. This claim is supported by previous studies [59] that show that wildfires may cause tree mortality lasting up to 3 years after the fire, thus, higher burn frequencies provide little time for vegetation recovery. Higher burn frequency (> 40% over 9-years period) during DJF at latitudes between 0 • and~12 • N, and during JJA and SON between latitudes 0 • and 25 • S may result in permanent loss of habitats for millions of plants, birds and mammal species, forest fragmentation, deforestation, as well as substantial economic losses.  In order to confirm the vegetation types that were burned on an annual basis per season, we extracted vegetation cover types from the annual CCI-LC data from 2007 to 2015 based on seasonally aggregated BA. As shown in Figure 4, the BA varies by vegetation cover types, year and season. The results show that the greatest proportion of BA across the entire period of study and seasons is dominated by the burning of Broadleaved deciduous trees, Shrublands, Rainfed croplands, Grasslands and Mosaic cropland/natural vegetation. This study found that the burning of Broadleaved deciduous trees is the highest in JJA (i.e., total BA of ~ 1 716,855 km 2 ), followed by DJF In order to confirm the vegetation types that were burned on an annual basis per season, we extracted vegetation cover types from the annual CCI-LC data from 2007 to 2015 based on seasonally aggregated BA. As shown in Figure 4, the BA varies by vegetation cover types, year and season. The results show that the greatest proportion of BA across the entire period of study and seasons is dominated by the burning of Broadleaved deciduous trees, Shrublands, Rainfed croplands, Grasslands and Mosaic cropland/natural vegetation. This study found that the burning of Broadleaved deciduous trees is the highest in JJA (i.e., total BA of~1,716,855 km 2 ), followed by DJF (i.e., total BA of~1,394,991 km 2 ) and SON (total BA of~1,124,044 km 2 ). On the contrary, the results show that BA for Shrublands is the highest in SON (i.e. total BA of~843,180 km 2 ), followed by DJF (i.e., total BA of~699,533 km 2 ) and JJA (i.e., total BA of~580,472 km 2 ). Burning of woody vegetation covers such as forests and Shrublands is a common traditional land management practice in Africa for satisfying the requirement of land for food production, particularly by smallholder farmers who often do not have equipment and machinery for land clearing [60]. Under favourable conditions (i.e., higher fuel loads supported by low precipitation, humidity, high air temperatures and strong winds), such fires may spread rapidly, extensively destroying other vegetation cover types in their path [56]. The combustion is also sustained by understory herbaceous vegetation cover in both open and closed forest canopies and thin bark characteristic of closed-canopy forest trees [61,62]. This claim is supported by the results in this study that show that predominant burning of land covers such as Mosaic herbaceous cover/tree and shrubs and Mosaic tree and shrub/herbaceous cover across all seasons.  Furthermore, the results show that Grasslands are largely burned during JJA and SON, evidenced by higher BA with a total BA of ~305,757 km 2 and ~208,469 km 2 , respectively. Moreover, Rainfed croplands exhibits higher BA during the DJF and SON with a total BA of ~216,754 km 2 and ~179,953 km 2 , respectively. The burning of Grasslands is often used as a range and pasture management tool [63]. However, the observed higher frequency compounded by prolonged dry conditions during the JJA and SON may result in reduced livestock production due to forage shortages, the proliferation of alien invasive species, bush encroachment, soil erosion and eventually, land degradation [64][65][66]. Studies by [67] showed that for particular grasslands such as montane grasslands, regular (i.e., annual) burning is beneficial for long-term productivity, while other studies [66,68,69] indicated that the burning of semi-arid grasslands may have drastic impacts on grazing capacity, basal cover, productivity and increases the intensity of seasonal droughts. On the other hand, extensive and frequent burning of Rainfed croplands during DJF and SON can be attributed to the burning of post-harvest crop residuals in preparation for the next season [70,71]. The results also suggest that the burning of croplands often extends to other land cover types such as natural Furthermore, the results show that Grasslands are largely burned during JJA and SON, evidenced by higher BA with a total BA of~305,757 km 2 and~208,469 km 2 , respectively. Moreover, Rainfed croplands exhibits higher BA during the DJF and SON with a total BA of~216,754 km 2 and~179,953 km 2 , respectively. The burning of Grasslands is often used as a range and pasture management tool [63]. However, the observed higher frequency compounded by prolonged dry conditions during the JJA and SON may result in reduced livestock production due to forage shortages, the proliferation of alien invasive species, bush encroachment, soil erosion and eventually, land degradation [64][65][66]. Studies by [67] showed that for particular grasslands such as montane grasslands, regular (i.e., annual) burning is beneficial for long-term productivity, while other studies [66,68,69] indicated that the burning of semi-arid grasslands may have drastic impacts on grazing capacity, basal cover, productivity and increases the intensity of seasonal droughts. On the other hand, extensive and frequent burning of Rainfed croplands during DJF and SON can be attributed to the burning of post-harvest crop residuals in preparation for the next season [70,71]. The results also suggest that the burning of croplands often extends to other land cover types such as natural vegetation, as shown by predominant burning of Mosaic cropland/natural vegetation and Mosaic natural vegetation/cropland across all seasons. Overall, the results have implications for regional agricultural practices, as well as informing wildfire management and conservation strategies.
MAM (see Figure 4b) exhibits the lowest proportions of burned vegetation cover types compared to other seasons, however, it has similar profiles to other seasons. The least burned vegetation cover types include Broadleaved evergreen forest, Irrigated croplands and Sparse vegetation across all seasons.

Seasonal Distribution of BC, CO and Smoke
In addition to devastating ecological impacts, studies have reported the impacts of wildfires on the chemical composition of the atmosphere, visibility, air quality and climate [15][16][17]. Contrarily to previous studies [33], this study provides a complementary and detailed analysis of BC, CO and smoke, their sources and variations across seasons. The results (see Figure 5) show an intense plume of BC, smoke and CO that extends between the 0 • and 18 • S region over the Atlantic Ocean during JJA season. This is as a result of a large distribution and frequency of wildfires in this region during JJA season, as seen in Figure 3. Vast quantities of agricultural land, forest and grass, in this region, are burned on an annual basis (Figure 4), leading to observed high emissions of BC, smoke and CO. However, other activities are present throughout the year that and contribute to the emissions of BC, smoke and CO such as domestic fires used for cooking, motor vehicle traffic and industrial activities [72]. The results in Figure 5 further show that JJA is the only season where BC, smoke and CO have equal dominance in the 0 • to 18 • S region. On the contrary, smoke is more dominant in the same region than BC and CO during SON. This is consistent with the relatively smaller distribution of burned areas during the SON season, as seen in Figure 3. Relatively reduced biomass burning during SON causes reciprocal lower emission of BC and CO. On the other hand, intense plumes of BC and CO, and a less intense plume of smoke are observed during the DJF in the 8 • N to 12 • N region. This can be attributed to the main biomass burning activities in this region that include forest and maquis fires, burning of agricultural wastes, and the domestic use of firewood [73]. In agreement, our results, in Figure 4, indicate that the predominantly burning vegetation cover types are Broadleaved deciduous trees, Shrublands and Croplands. It is known that smoke properties vary between fires depending on fuel type and moisture, combustion phase, wind conditions, and several other variables [74]. This might be one of the reasons why lower emissions of smoke are observed during DJF season in the 8 • N to 12 • N region. As an example, in Figure 1, the region between 0 • and 18 • S consists of mostly sparse trees, needled and broad-leaved trees, whereas the region between 8 • N and 12 • N mainly consists of grasslands, mosaic trees and shrubs. It is anticipated that the different fuel types in these regions produce different amounts and properties of smoke, hence the difference in the detection of smoke in these regions [75]. Overall, MAM has the lowest emission of BC, smoke and CO. Lower emissions of BC, smoke and CO coincides with scattered distribution and relatively smaller burned areas as seen in Figure 3. Other emissions probably come from daily activities such as domestic burning of wood for cooking and motor vehicle emissions. It has been shown that vehicles contribute significantly to the total anthropogenic BC emissions and urban BC concentrations [76]. For example, Efe [77] noted sources such as traffic in highly populated countries such as Nigeria, industry, and oil exploration activities as the major sources of high aerosols in the northern SSA region. However, these emissions sources do not have a strong seasonal pattern at a regional SSA scale.   (2) transport of smoke aerosols from neighbouring areas contributing to the magnitude of aerosol loading [78]. These results have implications for air quality and subsequently human health and visibility in the SSA region [79]. When smoke aerosols move to higher altitudes, they have a high potential for long-range transport by winds and being transported to downwind areas [80]. Another notable smoke aerosol plume is observed at altitudes between ~5.6 and 6.3 km for all seasons except MAM. The source of these smoke plumes is likely attributed to biomass burning. A strong peak at ~3 km is also observed seasonally. This peak is probably due to the high biomass burning during the drier conditions [81].   Figure 6 shows the seasonal averaged vertical profiles of smoke extinction over the study area for the period 2007-2016. Thick smoke aerosol plumes (~2 km) are observed for every season at altitudes between~2.4 and 4.5 km. The peaks between 2.4 and 4.5 km could be due to (1) the development of the planetary boundary layer and (2) transport of smoke aerosols from neighbouring areas contributing to the magnitude of aerosol loading [78]. These results have implications for air quality and subsequently human health and visibility in the SSA region [79]. When smoke aerosols move to higher altitudes, they have a high potential for long-range transport by winds and being transported to downwind areas [80]. Another notable smoke aerosol plume is observed at altitudes between~5.6 and 6.3 km for all seasons except MAM. The source of these smoke plumes is likely attributed to biomass burning. A strong peak at~3 km is also observed seasonally. This peak is probably due to the high biomass burning during the drier conditions [81].

The Effect of Meteorological Conditions on the Spatial and Seasonal Changes of Wildfires
Meteorological conditions such as temperature and humidity show remarkable variability throughout the year, thus, to the patterns of the fire. It is known that the flame over combustible materials is easy to spread at low humidity and difficult at high humidity [82]. In addition, the temperature of the fire will affect the humidity in the air and moisture content within combustible materials [82]. Fires are able to increase surface air temperature predominantly by reducing the latent heat flux due to fire-induced damage to the vegetation [83]. We used precipitation, latent heat flux and relative humidity data to explain the observed seasonal variability of burned area and emissions. The results (see Figure 7) show the seasonal averaged latent heat flux, precipitation and relative humidity over the study area. Latent heat flux is the energy component that best illustrates how vegetation dynamics impact the energy balance of a forest following a fire. The latent heat flux is driven by evaporation of liquid water out of the intercellular spaces within the leaves through the stomata, which are channels from the leaf interior to its surface. When there is no vegetation, the latent heat flux is low, conversely, vegetated surfaces will result in high latent heat flux. Therefore, the recovery of vegetation will result in an increase in latent heat. During the JJA season in the 0-18° S region a low latent heat flux (~40 Wm −2 ), a low precipitation (~50 mm/month) and low relative humidity (~40%) is observed. This is the season with the maximum emissions of BC, smoke and CO (as seen in Figure 4), and widespread BA as well as Burn frequency (as seen in Figure 3). The low latent heat implies an increase in surface air temperature which indicates an increase in the likelihood of intense fires. The intense fire will likely burn more biomass material and inject more BC, smoke and CO aerosols into the atmosphere. The low precipitation and relative humidity imply drier conditions which favour the spread of wildfires. On the other hand, higher values of latent heat (~140 Wm −2 ), precipitation (~250 mm/month) and relative humidity (60%) in the 0-18° S region are observed in the DJF season. This corresponds to the lowest emissions of BC, smoke and CO (as seen in Figure  4). The wetter conditions in the DJF would make it slightly difficult for the burning of any biomass. Hence, the observed lower BA in this region (0-18° S). On the other hand, the 8° N to 12° N region has lower values of latent heat (~80 Wm −2 ), precipitation (>~40 mm/month) and relative humidity (~30%) during the DJF season. This is the period with the maximum biomass burning (as seen in Figure 3), mainly of deciduous, broadleaved trees, shrubland and rainfed crops (Figure 4), as well as

The Effect of Meteorological Conditions on the Spatial and Seasonal Changes of Wildfires
Meteorological conditions such as temperature and humidity show remarkable variability throughout the year, thus, to the patterns of the fire. It is known that the flame over combustible materials is easy to spread at low humidity and difficult at high humidity [82]. In addition, the temperature of the fire will affect the humidity in the air and moisture content within combustible materials [82]. Fires are able to increase surface air temperature predominantly by reducing the latent heat flux due to fire-induced damage to the vegetation [83]. We used precipitation, latent heat flux and relative humidity data to explain the observed seasonal variability of burned area and emissions. The results (see Figure 7) show the seasonal averaged latent heat flux, precipitation and relative humidity over the study area. Latent heat flux is the energy component that best illustrates how vegetation dynamics impact the energy balance of a forest following a fire. The latent heat flux is driven by evaporation of liquid water out of the intercellular spaces within the leaves through the stomata, which are channels from the leaf interior to its surface. When there is no vegetation, the latent heat flux is low, conversely, vegetated surfaces will result in high latent heat flux. Therefore, the recovery of vegetation will result in an increase in latent heat. During the JJA season in the 0-18 • S region a low latent heat flux (~40 Wm −2 ), a low precipitation (~50 mm/month) and low relative humidity (~40%) is observed. This is the season with the maximum emissions of BC, smoke and CO (as seen in Figure 4), and widespread BA as well as Burn frequency (as seen in Figure 3). The low latent heat implies an increase in surface air temperature which indicates an increase in the likelihood of intense fires. The intense fire will likely burn more biomass material and inject more BC, smoke and CO aerosols into the atmosphere. The low precipitation and relative humidity imply drier conditions which favour the spread of wildfires. On the other hand, higher values of latent heat (~140 Wm −2 ), precipitation (~250 mm/month) and relative humidity (60%) in the 0-18 • S region are observed in the DJF season. This corresponds to the lowest emissions of BC, smoke and CO (as seen in Figure 4). The wetter conditions in the DJF would make it slightly difficult for the burning of any biomass. Hence, the observed lower BA in this region (0-18 • S). On the other hand, the 8 • N to 12 • N region has lower values of latent heat (~80 Wm −2 ), precipitation (>~40 mm/month) and relative humidity (~30%) during the DJF season. This is the period with the maximum biomass burning (as seen in Figure 3), mainly of deciduous, broadleaved trees, shrubland and rainfed crops (Figure 4), as well as emission of BC, smoke and CO. DJF has the lowest seasonal precipitation in the 8-12 • N region, implying that this is the driest season. Under warm and dry conditions with low relative humidity, wildfires are easily started. The higher values of latent heat (~140 Wm −2 ), precipitation (up to~350 mm/month) and relative humidity (~60%) are observed in MAM and SON seasons. This corresponds to the low BA and emissions of BC, smoke and CO in this region. Turco et al. [84] found that the fires activity during the wet conditions is not driven by flammability, rather, it is mainly driven by fuel availability and connectivity. Although the factors that affect wildfires considered in the current study were not exhaustive, the results clearly demonstrate and confirm that meteorological conditions are important in the fire regime, which directly affects the frequency of burning, burned area, and emissions of BC, smoke and CO. Future studies should incorporate various variables in addition to meteorological conditions in order to explain complex relationships between climate, vegetation patterns, wildfires and emissions. Drier conditions result in very intensive fires, rapid spreading of fires and burning of extensive biomass material due to availability of fuel loads [85]. Conversely, wetter conditions result in less intense fires, thus, lower burn area, burn frequency, and lower emissions of BC, smoke and CO.  [84] found that the fires activity during the wet conditions is not driven by flammability, rather, it is mainly driven by fuel availability and connectivity. Although the factors that affect wildfires considered in the current study were not exhaustive, the results clearly demonstrate and confirm that meteorological conditions are important in the fire regime, which directly affects the frequency of burning, burned area, and emissions of BC, smoke and CO. Future studies should incorporate various variables in addition to meteorological conditions in order to explain complex relationships between climate, vegetation patterns, wildfires and emissions. Drier conditions result in very intensive fires, rapid spreading of fires and burning of extensive biomass material due to availability of fuel loads [85].
Conversely, wetter conditions result in less intense fires, thus, lower burn area, burn frequency, and lower emissions of BC, smoke and CO.

Seasonal Trend Analysis over the Entire SSA Region
In order to further understand the temporal variability of wildfires, we employed a linear regression analysis and a Mann-Kendall Trend test in BC, CO, precipitation, latent heat flux, relative humidity and BA.

Seasonal Trend Analysis over the Entire SSA Region
In order to further understand the temporal variability of wildfires, we employed a linear regression analysis and a Mann-Kendall Trend test in BC, CO, precipitation, latent heat flux, relative humidity and BA.

Linear Regression over the Entire Region
Linear regression was applied for a trend analysis of BC concentration, CO concentration, latent heat flux, precipitation, relative humidity and BA in sub-Saharan Africa. The annual time series of these parameters averaged over the whole study period are illustrated in Figure 8. The corresponding interpolated regression lines (in red colour) are also plotted. The variability around a mean value of~4.5 × 10 −10 kg.m −3 for BC concentration is noticeable, and an increase in the annual average BC concentration is evident, given the slope of the regression line. Marked increases in CO concentration and latent heat flux are also observed from the slope of the regression lines, respectively. The increase in BC and CO concentrations are inconsistent with linear regression results for BA, which show no change over time. Therefore, the increase in BC and CO could be due to an increase from other emission sources such as biomass cooking stoves, industry, biomass heating stoves, and transport [77]. In fact, the results showing high concentrations of BC, and CO are consistent with Mehta et al. [86] that found high AOD during the JJA and DJF seasons. The increase in latent heat flux indicates drier conditions leading to a favourable environment for fires. Precipitation, on the other hand, showed inter-annual variability around a mean value of~80 mm/month but exhibits a decreasing trend. The decrease in precipitation implies drier conditions which also favours the spread of fires. Relative humidity shows inter-annual variability around a mean value of~41% with a slightly increasing trend. Generally, the lower the relative humidity, the more likely that a fire will start and burn, whereas when humidity is high, evapotranspiration is lower resulting in high moisture preservation by vegetation. Similar to relative humidity, burned area (BA) shows a slightly decreasing trend over time. The decrease can be attributed to the limitation of MODIS BA in detecting fires less than its pixel size, i.e., 500 m. Linear regression was applied for a trend analysis of BC concentration, CO concentration, latent heat flux, precipitation, relative humidity and BA in sub-Saharan Africa. The annual time series of these parameters averaged over the whole study period are illustrated in Figure 8. The corresponding interpolated regression lines (in red colour) are also plotted. The variability around a mean value of ~4.5 × 10 −10 kg.m −3 for BC concentration is noticeable, and an increase in the annual average BC concentration is evident, given the slope of the regression line. Marked increases in CO concentration and latent heat flux are also observed from the slope of the regression lines, respectively. The increase in BC and CO concentrations are inconsistent with linear regression results for BA, which show no change over time. Therefore, the increase in BC and CO could be due to an increase from other emission sources such as biomass cooking stoves, industry, biomass heating stoves, and transport [77]. In fact, the results showing high concentrations of BC, and CO are consistent with Mehta et al. [86] that found high AOD during the JJA and DJF seasons. The increase in latent heat flux indicates drier conditions leading to a favourable environment for fires. Precipitation, on the other hand, showed inter-annual variability around a mean value of ~80 mm/month but exhibits a decreasing trend. The decrease in precipitation implies drier conditions which also favours the spread of fires. Relative humidity shows inter-annual variability around a mean value of ~41% with a slightly increasing trend. Generally, the lower the relative humidity, the more likely that a fire will start and burn, whereas when humidity is high, evapotranspiration is lower resulting in high moisture preservation by vegetation. Similar to relative humidity, burned area (BA) shows a slightly decreasing trend over time. The decrease can be attributed to the limitation of MODIS BA in detecting fires less than its pixel size, i.e., 500m.

Mann-Kendall Trend Analysis over the Entire Region
The MK trend test was run on the annual BC, CO, latent heat flux, precipitation, relative humidity and BA parameters, to test the null hypothesis (H0) that there is no significant trend in the time series at a significance level α = 0.05. The results are shown in Table 1. If the p-value is less than α, H0 is rejected. Rejecting H0 indicates that there is a trend in the time series, while failure to reject H0 (i.e., p-value > 0.05) indicates that no trend was detected. On rejecting the null hypothesis, the time series is said to be statistically significant. The results of the MK test (Table 1) reveal that there is no evidence to support the alternative hypothesis, hence, H0 cannot be rejected (p-value > 0.05) for the time series of BC, relative humidity and burned area. On the other hand, H0 is rejected (p-value < 0.05) for CO, latent heat flux and precipitation. These trends were further analysed using the sequential MK test. The results are consistent with the linear regression results that indicate an increasing linear trend in CO, latent heat flux and precipitation, with the exception of BC, which indicate no significant trend despite increasing linear trend. This may be due to the lower inter-annual variability observed.

Mann-Kendall Trend Analysis over the Entire Region
The MK trend test was run on the annual BC, CO, latent heat flux, precipitation, relative humidity and BA parameters, to test the null hypothesis (H 0 ) that there is no significant trend in the time series at a significance level α = 0.05. The results are shown in Table 1. If the p-value is less than α, H 0 is rejected. Rejecting H 0 indicates that there is a trend in the time series, while failure to reject H 0 (i.e., p-value > 0.05) indicates that no trend was detected. On rejecting the null hypothesis, the time series is said to be statistically significant. The results of the MK test (Table 1) reveal that there is no evidence to support the alternative hypothesis, hence, H 0 cannot be rejected (p-value > 0.05) for the time series of BC, relative humidity and burned area. On the other hand, H 0 is rejected (p-value < 0.05) for CO, latent heat flux and precipitation. These trends were further analysed using the sequential MK test. The results are consistent with the linear regression results that indicate an increasing linear trend in CO, latent heat flux and precipitation, with the exception of BC, which indicate no significant trend despite increasing linear trend. This may be due to the lower inter-annual variability observed.
The seasonal MK test for the parameters was also investigated and the results are shown in Table 2. For relative humidity, the MK test failed to reject the null hypothesis (p-value > 0.05) for all seasons, which is in agreement with the annual MK test in Table 1. For the latent heat flux, the null hypothesis is rejected (p-value < 0.05) for all seasons. Unlike relative humidity and latent heat flux, BC shows no significant trend in DJF and SON seasons (p-value > 0.05) while during MAM and JJA seasons, it showed a significant trend (p-value < 0.05). The BC data shows high seasonal variability in term of trends but no significant annual trend (as can be seen in Table 1). The CO data show a significant annual trend (p-value < 0.05) but no significant seasonal trend except during MAM and JJA seasons. On the other hand, precipitation showed no seasonal significant trend with the exception of DJF season.

Sequential Mann-Kendall Analysis over the Entire Region
The SQMK test was applied for annual and seasonal parameters to identify how the trend gets fluctuated over the study period. The graphical illustration of the results of the sequential MK test for BC concentration is shown in Figure 9. There was no significant trend detected in BC concentration in the annual data, and the DJF and SON seasons. The multiple turning points on these datasets is one of the reasons for the lack of trend. Another reason for the lack of trend could be what seems to be a convergence of the progressive and retrograde series in 2015. The MAM and JJA seasons show increasing trends of the BC concentration. The relative humidity dataset (see Figure 10) did not show any trend in the annual dataset and for all seasons. The reason that there is no trend for the annual and DJF datasets is that there is no intersection point between the progressive and retrograde series. For the MAM, JJA and SON seasons, there are multiple intersections points between the progressive and retrograde series and these curve do not diverge to show a trend. The latent heat flux dataset (see Figure 11) shows a significant increasing trend for the annual dataset and all the seasons. The increase in latent heat flux implies drier environmental conditions.

Seasonal Trend Analysis over Northern and Southern Sub-regions of the SSA
In order to study, in detail, the sub-regions and seasons where large plumes of CO, smoke and BC were observed, a trend analysis was carried out over northern SSA region (0-11° N) and the

Seasonal Trend Analysis over Northern and Southern Sub-regions of the SSA
In order to study, in detail, the sub-regions and seasons where large plumes of CO, smoke and BC were observed, a trend analysis was carried out over northern SSA region (0-11° N) and the

Seasonal Trend Analysis over Northern and Southern Sub-regions of the SSA
In order to study, in detail, the sub-regions and seasons where large plumes of CO, smoke and BC were observed, a trend analysis was carried out over northern SSA region (0-11° N) and the

Seasonal Trend Analysis over Northern and Southern Sub-regions of the SSA
In order to study, in detail, the sub-regions and seasons where large plumes of CO, smoke and BC were observed, a trend analysis was carried out over northern SSA region (0-11 • N) and the southern SSA region (0-17 • S). This was done to exclude the influence to the trend and patterns that may be coming from areas and seasons where there is limited fire activity when the entire SSA region is analysed.
The results in Table 3 show the MK test for BC, latent heat flux and relative surface humidity over the northern and southern SSA regions respectively. On the other hand, Figures 12 and 13 show the results of linear regression and SQMK test, performed to determine the direction and magnitude of the trend since this information could not be obtained from MK test alone. Generally, the results during DJF for the northern SSA region are consistent with the results over the entire region (see Tables 1 and 2) that show no significant trend in BC and relative surface humidity.  On the other hand, MK test results over the southern SSA region are different from the result over the entire region. In particular, the results over the southern SSA region show no significant trend in BC, latent heat flux and relative humidity (Table 3), as opposed to SSA wide results that indicate significant trends for BC and latent heat flux. Moreover, the results of the trend analysis using linear regression and the SQMK test (Figure 13a-c) over the southern SSA region show an increasing trend in BC and latent heat flux, attributable to vast burning activities during the dry season (i.e., JJA) [94]. However, SQMK shows that the observed trend is not significant, attributed to better implementations of environmental policies in this sub-region [86]. On the other hand, the linear regression for relative humidity showed a decreasing trend over time (see Figure 13c). But the SQMK test showed no trend for this parameter (see Figure 13f). These results are consistent with results over the entire region (Figures 8, 9, 10, and 11). On the other hand, MK test results over the southern SSA region are different from the result over the entire region. In particular, the results over the southern SSA region show no significant trend in BC, latent heat flux and relative humidity (Table 3), as opposed to SSA wide results that indicate significant trends for BC and latent heat flux. Moreover, the results of the trend analysis using linear regression and the SQMK test (Figure 13a-c) over the southern SSA region show an increasing trend in BC and latent heat flux, attributable to vast burning activities during the dry season (i.e., JJA) [94]. However, SQMK shows that the observed trend is not significant, attributed to better implementations of environmental policies in this sub-region [86]. On the other hand, the linear regression for relative humidity showed a decreasing trend over time (see Figure 13c). But the SQMK test showed no trend for this parameter (see Figure 13f). These results are consistent with results over the entire region (Figures 8, 9, 10, and 11). Overall, the results are relevant for addressing long-term wildfire information scarcity in Africa, developing fire and emissions inventories for SSA and informing regional and national policies in line with the Sustainable Development Goals (SDGs) and other international mandates to reduce and control emissions. Similarly, the results of linear regression show a slight increase in BC and relative humidity given the slope of linear regression (Figure 12a-c). This increase in BC in the northern SSA region is due to biomass burning from north to south of the Sahel region that occurs in DJF [87]. However, multiple turning points in the SQMK test results indicate that there is no trend. The lack of trend in the northern SSA region can be attributed to the greening of the Sahel region, as reported by various studies [88][89][90]. The trend analysis of latent heat flux using both linear regression and the SQMK test shows a significantly increasing trend, evidenced by the single turning point in 2006. In general, it is known that the latent heat flux increase is largely due to global warming which is driven by the increase in temperature. An overview by Collins [91] shows an increasing trend in temperature in the DJF season for the period 2001-2010 in the SSA region. Furthermore, Serdeczny et al. [92] used the high-emission scenario (RCP8.5) and low-emission scenario (RCP2.6) models to show the increasing trend of temperature in the SSA over time. This steady increase of temperature overtime is one factor responsible for the increase in the latent heat flux. Another factor that could contribute to the increase in latent heat flux is the increase in carbon dioxide (CO 2 ). It is also known that CO 2 is an important gas responsible for regulating the Earth's temperature, where considerable amounts of CO 2 result in an increase in temperatures. Contrarily to the global CO 2 emission change over the period of 1751 to 2017 observed by Ritchie and Roser [93], this study shows an increase of CO 2 emission (1.08 million ton) starting from 2006 over the SSA region for the period 2000 to 2016. This coincides with the period when the increasing trend of latent heat flux is observed. Therefore, the observed increase in CO 2 emission have devastating implications for global warming and wildfire frequency, and has caused increase in latent heat flux. Lastly, the El Niño climatic events that occurred in 2006-2007, 2009-2010 and 2015-2016 contributed to a significantly increasing trend in latent heat [56]. This is in agreement with an annual and seasonal analysis over the entire SSA region that shows no significant trend (Figures 8-11).
On the other hand, MK test results over the southern SSA region are different from the result over the entire region. In particular, the results over the southern SSA region show no significant trend in BC, latent heat flux and relative humidity (Table 3), as opposed to SSA wide results that indicate significant trends for BC and latent heat flux. Moreover, the results of the trend analysis using linear regression and the SQMK test (Figure 13a-c) over the southern SSA region show an increasing trend in BC and latent heat flux, attributable to vast burning activities during the dry season (i.e., JJA) [94]. However, SQMK shows that the observed trend is not significant, attributed to better implementations of environmental policies in this sub-region [86]. On the other hand, the linear regression for relative humidity showed a decreasing trend over time (see Figure 13c). But the SQMK test showed no trend for this parameter (see Figure 13f). These results are consistent with results over the entire region (Figures 8-11).
Overall, the results are relevant for addressing long-term wildfire information scarcity in Africa, developing fire and emissions inventories for SSA and informing regional and national policies in line with the Sustainable Development Goals (SDGs) and other international mandates to reduce and control emissions.

Conclusions
Although wildfires are crucial for promoting diversity and natural vegetation regeneration [95], their frequency and severity may be devastating for particular vegetation types, causing land degradation, loss of habitats, deforestation and supporting the proliferation of alien invasive plant species. At the same time, wildfires increase the global atmospheric load of trace gases and aerosols, thus affecting air quality, visibility and Earth's radiation balance. This study sought to analyse the surface-atmospheric effects of wildfires in sub-Saharan Africa using multisource remotely sensed, reanalysis data and a trend analysis. The results of this study are in agreement with previous studies [33,[96][97][98], however, this study provided a more comprehensive analysis of simultaneous effects of wildfires on the land surface and atmosphere than most previous studies. For example, a study by Jeong and Wang [96] only considered BC seasonal variations at a global scale, whilst [97] showed seasonality and distributions of active fires at a global scale. Additionally, multisource remotely sensed data used in this study were essential in explaining observed long-term and seasonal patterns of wildfires and associated emissions.
Firstly, the results show that frequent wildfires, especially during prolonged dry conditions such as during El Niño events, may result in the disturbance of forest and savannah ecosystems. Additionally, Croplands and Grasslands were amongst the largest burned vegetation cover types after Broadleaved deciduous trees and Shrublands, attributed to agricultural land management practices such as "slash and burn" and pasture management. An analysis of a 9-year period shows that the highest burn frequency, i.e., > 40% in the regions 0-20 • N and 0-25 • S, thus indicates a threat to the sustainability of natural forests and rangelands in Africa. However, the burned area might have been underestimated due to the limitations of MODIS spatial resolution in detecting small fires. Future studies should consider fusing various datasets such as MODIS, Landsat and Sentinel-2 to benefit from the improved spatial and temporal resolutions.
Secondly, large emissions and distributions of smoke, BC and CO in the 0-18 • S region are observed during the JJA season as a result of agricultural land, forest and grass burning. Drier conditions influenced by the low latent heat flux, low precipitation and low relative humidity aid in the intensity of the fire. Conversely large emissions and distributions of smoke, BC and CO are observed in the 8-12 • N region during the DJF season. This is consistent with Jeong and Wang [96]. The MAM season shows the least loading of smoke, BC and CO both in the 0-18 • S and 8-12 • N regions. This is the season that corresponds to the least burned vegetation cover types. This season also comprises of moderate latent heat flux (~20 Wm −2 ), moderate relative humidity (~60%) and low precipitation (~200 mm/month) which are ideal conditions for intensive fires to occur.
The long-term trend analysis showed linear increasing trends in BC, CO, latent heat flux and relative humidity datasets. However, in the northern SSA region, MK and SQMK (95% confidence interval) revealed no significant increasing trend in BC attributed to the greening of the Sahel region studies [88][89][90]. Similarly, there was no significant trend in the southern SSA region attributed to better implementations of environmental policies [86].
Overall, the results in this study necessitate better regional fire management and air quality control strategies and enforcement to combat the devastating effects of wildfires, especially during DJF, JJA and SON, to preserve endangered species and habitats, promote sustainable land management, and reduce greenhouse gases (GHG) emissions.