Biomass Burning and Water Balance Dynamics in the Lake Chad Basin in Africa

: The present study investigated the effect of biomass burning on the water cycle using a case study of the Chari–Logone Catchment of the Lake Chad Basin (LCB). The Chari–Logone catchment was selected because it supplies over 90% of the water input to the lake, which is the largest basin in central Africa. Two water balance simulations, one considering burning and one without, were compared from the years 2003 to 2011. For a more comprehensive assessment of the effects of burning, albedo change, which has been shown to have a signiﬁcant impact on a number of environmental factors, was used as a model input for calculating potential evapotranspiration (ET). Analysis of the burning scenario showed that burning grassland, which comprises almost 75% of the total Chari–Logone land cover, causes increased ET and runoff during the dry season (November– March). Recent studies have demonstrated that there is an increasing trend in the LCB of converting shrubland, grassland, and wetlands to cropland. This change from grassland to cropland has the potential to decrease the amount of water available to water bodies during the winter. All vegetative classes in a burning scenario showed a decrease in ET during the wet season. Although a decrease in annual precipitation in global circulation processes such as the El Niño Southern Oscillation would cause droughts and induce wildﬁres in the Sahel, the present study shows that a decrease in ET by the human-induced burning would cause a severe decrease in precipitation as well.


Introduction
The African Sahel experienced a dramatic decrease in precipitation in the 1960s to the 1980s, which has drastically affected the size of Lake Chad [1]. Lake Chad gained widespread attention from the scientific community and the general public due to its rapid shrinkage from 25,000 km 2 in 1963 to less than 3000 km 2 in 2008 [2]. Lake Chad's shrinkage is significantly detrimental to regional sustainability, because it is an economically important water resource upon which agricultural and fishing activities depend. In a series of surveys interviewing 25,000 households in the southwestern portion of the lake, 59% of the households earned three-quarters of their income from farming, 36% earned income from a combination of fishing and farming, and 5% relied entirely on fishing [3]. A number of factors have been attributed to shrinkage of the lake, including a decrease in precipitation [4], poor water management practices [5], and land use changes [6]. Most of the rainfall in the LCB is immediately lost to evapotranspiration (ET) [7], leaving the region sensitive to precipitation anomalies. Precipitation in the region is bi-seasonal, with wet season monsoonal rains lasting from April to October, and the dry season the rest of the year. The monsoonal rains are driven by the Intertropical Convergence Zone (ITCZ), created by the convergence of trade winds from the northern and southern hemispheres and forming the ascending branch of the Hadley circulation [8]. Precipitation in the Sahel has been shown to be affected by the El Niño Southern Oscillation (ENSO) system, which is a cycle of ocean and atmospheric interaction that produces an irregularly periodic variation in winds and sea surface temperatures over the tropical eastern Pacific Ocean. During El Niño years, the wind pattern associated with the West African Monsoon has been shown to create a dry condition across the Sahel region, leading to precipitation deficits [9]. High variability of rainfall due to these meteorological factors and other regional phenomena contributed to a series of droughts in the 1970s and 1980s, with a gradual trend of annual rainfall recovery occurring since the 1990s contributing to the recovery of Lake Chad [10].
Another factor which has received little attention is the occurrence of biomass burning, whether through wildfires or human-induced burning. Biomass burning has been shown to have significant impacts on the climate by altering land cover and vegetation [11][12][13], introducing aerosols and trace gases [14,15], and changing land surface characteristics, especially in regard to surface albedo [16,17]. Dintwe et al. [18] examined the effect of fire-induced albedo change on the surface energy balance in sub-Saharan Africa by using Moderate Resolution Imaging Spectrometer (MODIS) albedo data combined with vegetation cover classes and burn products. The short-term effects of albedo change due to burning contributed to an increase in global radiative forcings, of which the magnitudes were much greater than previously thought for regional climate variability in sub-Saharan Africa. Atchley et al. [19] examined the effect of forest fires on surface and subsurface water balance using the hydrological model ParFlow-CLM. The CLM component of the model solved ET using surface energy balance formulation. The study found that fire disturbance in high severity burn sites increased ET and surface runoff, whereas in low severity areas, reduced transpiration from burning was more significant than the increased runoff. Hodnebrog et al. [20] investigated the effect of biomass burning on precipitation in southern Africa using Global Precipitation Climatology Centre (GPCC) precipitation data in conjunction with black carbon and organic carbon measurements from 1850 to 2000. They showed the strong impact of anthropogenic biomass burning particles during the dry season and a reduction in precipitation by 20-30% over much of Africa from 1950 to 2000.
Ichoku's study [21] of biomass burning in northern sub-Saharan Africa (NSSA), showed an increased conversion of shrubland, grassland, and wetlands to cropland from 2003 to 2013. This in turn, affected the amount of biomass burning that occurred, because prescribed burning is one of the primary sources of biomass burning in NSSA. The study also found that areas with relatively high burn activity experienced much higher evapotranspiration and posited a relationship between biomass burning and decreased precipitation. However, there have been no studies specifically in the LCB examining the relationship between human-induced burning and the water cycle.
A large contributor to biomass burning within the LCB is the use of fire to remove vegetative growth at the end of a growing season and the clearing of grassland and savanna for agricultural use. In the early to late dry season (November-March), a patch mosaic burning regime to clear unwanted vegetation is used to clear savanna and woody areas and facilitate clearing land for planting crops. Burning primarily occurs during the dry season when rainfall is at a minimum to clear the wet season's growth or to convert the land for future agricultural use. Smoke byproducts of biomass burning have been known to produce significant amounts of black carbon, carbon monoxide, and carbon dioxide [14]. Furthermore, the resultant ash and charcoal deposition and change in land cover have been known to have a significant impact on surface albedo [16,22,23]. Surface albedo darkening from burning is caused by char deposition that increases the absorption of solar Earth 2021, 2 342 radiation [24]. Desales [17] found that changes in albedo coupled with a decrease in leaf area index (LAI) and vegetative land cover fraction from burning was linked to convective instability. The cause of this convective instability and precipitation loss is a cooling and drying of the atmosphere, leading to the weakening of upward atmospheric motion during the onset and mature stages of the monsoon.
Biomass burning effects on the water cycle in the Sahel have been addressed to some degree in regard to ET, soil moisture, and precipitation, but conclusive evidence is limited by issues associated with burning and vegetative seasonality [21]. Other studies associating water cycle changes with burning include its effect on soil infiltration rates, surface runoff [25], and soil erosion [26,27]. Previous studies on the water cycle of Sahelian regions have used models focused on terrestrial water storage (TWS) using data from satellites such as Gravity Recovery and Climate Experiment (GRACE) [28][29][30][31]. Although such models have helped to bring about great advances in understanding hydrologic dynamics on a climatic scale, their spatial resolution is too coarse for detailed characterizations of burned areas.
Lack of in situ data is largely due to the sparseness of meteorological stations and the fact that existing stations are inconsistent at recording data on a daily or even monthly basis [32][33][34]. Additionally, installing new stations is difficult in the LCB due to insecurity in the region. WetSpass, which stands for Water and Energy Transfer between Soil, Plants and Atmosphere under quasi-Steady State, is a physically based water balance model that is not as data-intensive as other water balance models [35,36]. Babamaaji [6] demonstrated the effective use of WetSpass in the LCB. One drawback of WetSpass is its inability to calculate water balance more frequently than seasonally. This, however, has been addressed in the newer version, WetSpass-M [36], which was designed specifically to calculate monthly water balance.
There has been little research on using a water balance model to calculate the effect of biomass burning on a hydrological system, especially in a data-scarce region in Africa. In order to accomplish this, the present study investigated burning effects on the water balance in the LCB by incorporating burning parameters such as surface albedo and land cover characteristics of burned areas into a water balance model. The present water balance modeling approach determines how biomass burning affects regional scale water cycles and interprets whether burning would affect the shrinkage of the lake. The time period chosen for this study was from 2003 to 2011 in order to use albedo change from burning, as suggested by [16]. This study investigates the correlation between increased evapotranspiration or increased precipitation due to burning. If there is a correlation, what vegetative land types are most greatly affected? Supposing a continued conversion of other land use types to cropland, what does this imply for water balance in the LCB? To answer these questions, it is necessary to use the model to simulate how burning affects water balance within the LCB, particularly examining water balance changes in vegetative areas. Determining which parameters are affected by burning and the percentage of land use within the study area are also important steps to building a comprehensive model.

Study Area
The Chari-Logone catchment (CLC), a sub-basin of the LCB located south of Lake Chad, was chosen for this study because it provides more than 90% of the water input to Lake Chad and is the main driver for the entire lake system [37]. The Chari and Logone rivers provide an average of 1946 mm/year to the lake, whereas direct rainfall contributes a mere 329 mm annually [38]. The spatial extent of the study area is 601,350 km 2 , and it is located within the borders of Cameroon, Chad, Central African Republic, and Sudan. The Chari River originates in the Central African Republic, while the Logone river does so on the Adamawa plateau in Cameroon. The catchment lies between 5 • N and 15 • N at the eastern part of the Sahel (Figure 1). The topography of the CLC is a plain of clayey soils in the northern region and laterite deposits in the south. The CLC is located in the Chad Graben, which formed prior to the upper Jurassic. The stratigraphic sequence of the Chad Graben is composed of fine and coarse sands along with reddish and gray clays from the Cretaceous period [39].
Ecologically, the area is in the Sudano-Sahelian savanna zone, covered by short grass savannah with a mixture of grasses, shrubs and trees [40]. Since the 1950s, a large proportion of the land cover of the Sahel region had evolved from savanna to agricultural land use [41]. The most significant land use/land cover transformations have taken place after 1975, particularly the dramatic increase in the areal extent of croplands [42]. An overall increase in runoff and erosion rates in the catchment of Lake Chad has been observed in conjunction with decreased vegetation cover, land clearing, and cropland emergence [41,43]. Land clearing has decreased the content of organic matter in the topsoil and enhanced surface crusting, which lead to an increase in runoff [44]. These changes are related to the demographic expansion of cropland and the heightened demand for agricultural production.
Earth2021, 2, FOR PEER REVIEW 4 eastern part of the Sahel (Figure 1). The topography of the CLC is a plain of clayey soils in the northern region and laterite deposits in the south. The CLC is located in the Chad Graben, which formed prior to the upper Jurassic. The stratigraphic sequence of the Chad Graben is composed of fine and coarse sands along with reddish and gray clays from the Cretaceous period [39]. Ecologically, the area is in the Sudano-Sahelian savanna zone, covered by short grass savannah with a mixture of grasses, shrubs and trees [40]. Since the 1950s, a large proportion of the land cover of the Sahel region had evolved from savanna to agricultural land use [41]. The most significant land use/land cover transformations have taken place after 1975, particularly the dramatic increase in the areal extent of croplands [42]. An overall increase in runoff and erosion rates in the catchment of Lake Chad has been observed in conjunction with decreased vegetation cover, land clearing, and cropland emergence [41,43]. Land clearing has decreased the content of organic matter in the topsoil and enhanced surface crusting, which lead to an increase in runoff [44]. These changes are related to the demographic expansion of cropland and the heightened demand for agricultural production.

The WetSpass-M: The Water Balance Model
WetSpass-M utilizes raster data, and calculations are performed on a cell-by-cell basis where the input data are geospatially overlayed and sub-divided into vegetated, bare soil, open water, and impervious surface fractions. Water balance calculations of each pixel are processed in the following order: interception (process 1); surface runoff (process 2); evapotranspiration (process 3); and recharge (process 4) ( Figure 2). Additionally, landuse/land-cover fraction values for each pixel are used as weighting factors for the calculation of the water balance [36].

The WetSpass-M: The Water Balance Model
WetSpass-M utilizes raster data, and calculations are performed on a cell-by-cell basis where the input data are geospatially overlayed and sub-divided into vegetated, bare soil, open water, and impervious surface fractions. Water balance calculations of each pixel are processed in the following order: interception (process 1); surface runoff (process 2); evapotranspiration (process 3); and recharge (process 4) ( Figure 2). Additionally, land-use/land-cover fraction values for each pixel are used as weighting factors for the calculation of the water balance [36].
WetSpass-M has demonstrated favorable results for modeling water balance in the western Sahel in the Black Volta Basin [36] and in northern Africa [45]. Calculations for the resolution of the input land cover and soil raster imagery are important to influence sub-cell land cover heterogeneity, because their classification is compared to a series of look-up tables for parameterization during water balance calculations [35]. The look-up tables include land use and land cover, soil type, and rainy days per month. For this study, additional land cover parameters were included in the land use look-up table to account for burning. These land cover parameters include LAI, vegetated area, burned area, impervious area, open water area, root depth, minimum stomata opening, vegetative height, Manning coefficient, land factor, and aerodynamic resistance. LAI change was calculated using the same method employed by [17], with original unburned vegetation cover, the fractional burned area, and vegetative survival rate as input parameters. Raster datasets required for input into WetSpass-M include precipitation, temperature, wind, soil, land cover, elevation, slope, and potential evapotranspiration (PET). The water balance model outputs are interception, actual evapotranspiration, surface runoff, and recharge, the summation of which are equal to precipitation. WetSpass-M calculates impervious or actual evapotranspiration (ET) using vegetation coefficients in conjunction with PET. WetSpass-M has demonstrated favorable results for modeling water balance in the western Sahel in the Black Volta Basin [36] and in northern Africa [45]. Calculations for the resolution of the input land cover and soil raster imagery are important to influence sub-cell land cover heterogeneity, because their classification is compared to a series of look-up tables for parameterization during water balance calculations [35]. The look-up tables include land use and land cover, soil type, and rainy days per month. For this study, additional land cover parameters were included in the land use look-up table to account for burning. These land cover parameters include LAI, vegetated area, burned area, impervious area, open water area, root depth, minimum stomata opening, vegetative height, Manning coefficient, land factor, and aerodynamic resistance. LAI change was calculated using the same method employed by [17], with original unburned vegetation cover, the fractional burned area, and vegetative survival rate as input parameters. Raster datasets required for input into WetSpass-M include precipitation, temperature, wind, soil, land cover, elevation, slope, and potential evapotranspiration (PET). The water balance model outputs are interception, actual evapotranspiration, surface runoff, and recharge, the summation of which are equal to precipitation. WetSpass-M calculates impervious or actual evapotranspiration (ET) using vegetation coefficients in conjunction with PET.

Datasets and Parameters for Burning
Ground station observations from the National Oceanic and Atmospheric Administration's (NOAA) National Center for Environmental Information (NCEI) database are available for sub-Saharan Africa. However, the stations in the study area are extremely sparse when compared to the United States or Europe, and monthly or even yearly data are not always recorded for many stations. African stations are also not always accurate, either due to instrument error or poor data management.
With such few ground stations, even if data are accurate, when the data are interpolated across large areas within the study area, they are skewed based on spatial distribution. Therefore, we adopted the precipitation data from the NASA's Tropical Rainfall Measurement Mission (TRMM) https://disc.gsfc.nasa.gov/datasets/TRMM_3B31_7/summary, accessed 17 April 2021. TRMM level 3 3B31 was used for the study, with a 0.5 degree × 0.5 degree resolution at a monthly time scale and units in mm/month. The data cover the area from 40° N to 40° S latitude, and were collected from January 1998 to April 2015, shortly before the satellite re-entered Earth's atmosphere. A number of studies have examined the accuracy of TRMM data compared to station data in different regions [46,47]. One such study by Ojo and Omotosho [48] looked specifically at the accuracy of TRMM in Nigeria. The study found that TRMM performed very well when compared to the few

Datasets and Parameters for Burning
Ground station observations from the National Oceanic and Atmospheric Administration's (NOAA) National Center for Environmental Information (NCEI) database are available for sub-Saharan Africa. However, the stations in the study area are extremely sparse when compared to the United States or Europe, and monthly or even yearly data are not always recorded for many stations. African stations are also not always accurate, either due to instrument error or poor data management.
With such few ground stations, even if data are accurate, when the data are interpolated across large areas within the study area, they are skewed based on spatial distribution. Therefore, we adopted the precipitation data from the NASA's Tropical Rainfall Measurement Mission (TRMM) https://disc.gsfc.nasa.gov/datasets/TRMM_3B31_7/summary, accessed 17 April 2021. TRMM level 3 3B31 was used for the study, with a 0.5 degree × 0.5 degree resolution at a monthly time scale and units in mm/month. The data cover the area from 40 • N to 40 • S latitude, and were collected from January 1998 to April 2015, shortly before the satellite re-entered Earth's atmosphere. A number of studies have examined the accuracy of TRMM data compared to station data in different regions [46,47]. One such study by Ojo and Omotosho [48] looked specifically at the accuracy of TRMM in Nigeria. The study found that TRMM performed very well when compared to the few meteorological stations available. TRMM did, however, have a tendency to marginally overestimate precipitation, especially toward southwestern Nigeria, which is far from the study area.
Soil rasters were extracted from the United Nation's Food and Agriculture Organization (FAO) Harmonized World Soil Database (HWSD). The HWSD was released in 2008, and is a 30 arc-second raster database with over 16,000 different soil mapping units that combine existing regional and national updates of soil information worldwide. Land cover data were derived from MODIS MCD12Q1 Collection 6. This Level 3 dataset is a land classification raster created from an algorithm combining data from the Moderate Resolution Imaging Spectrometer (MODIS) instrument onboard NASA's Terra and Aqua satellites. These satellites were launched December 1999 and May 2002, respectively. An advantage of MODIS is its ability to cover the entire globe every 1-2 days with a 2330 km wide swath. MODIS has thirty-six discrete spectral bands, ranging in wavelength from 0.4 µm to 14.4 µm, thereby allowing the retrieval of a wide variety of geophysical prod-ucts [49,50]. MODIS data have widely been used and validated in a number of studies involving Africa [21,[51][52][53]. MCD12Q1 is available yearly at a 1 km × 1 km resolution. The yearly availability of MCD12Q1 makes it preferable in this study over Globcover, which has shown better performance in Africa although only provides data from December 2004 to June 2006 and January to December 2009 [54].
The MCD12Q1 data are divided into eighteen distinct land classifications, which were converted into another set of classifications required by WetSpass-M. Of the eighteen MODIS classifications, only ten were found in the study area. A model was built in ArcGIS to convert the numbers signifying those ten MODIS classes into eight WetSpass-M classes ( Table 1). Two different land cover datasets were created, one considering burning and one without burning. The dataset considering burning combined the land classification raster with fire counts derived from MOD14/MYD14 Collection 5 fire data. MOD14 data are derived from a fire detection algorithm which takes the two 4 µm channels, numbered 21 and 22. Classification of a pixel as containing fire is only considered if it satisfies a variety of conditions which indicate a strong probability that it contains fire [55,56]. For each grid cell, Terra/Day or Terra/Night settings were selected based on which had the most fire detected; this value was then used for the day. Each daily total was then summed in order to obtain a monthly value per grid cell (Figure 3a,b). Although the monthly value represents one overpass per day, the data are not quantitatively representative of the diurnal burn pattern, however do provide statistically qualitative information. The fire data were then overlaid onto the land cover map using GIS to simulate which vegetated areas were likely affected by burning.
A GIS model was built, which first extracted the fire count for the desired month, and then a raster of only pixels detecting fire was created using the Raster Calculator tool. Specific land cover types were extracted from the previously created yearly land cover data used in the non-burning WetSpass-M model. The extracted land cover classes were cropland, savanna, mixed forest, deciduous forest, coniferous forest, and grassland ( Figure 4a). Each land cover area was then used as a mask to extract only those areas which experienced burning. For example, a new raster for savanna was created, only showing those areas in the savanna in which fire was detected for a particular month. Each of the newly created land cover rasters was then assigned a number corresponding to a WetSpass code used in the look-up table. Savanna and grassland were distinguished due to savanna containing more woody vegetation. The difference in vegetative cover affected the parameters used in the look-up table [Table S1]. For example, land cover roughness and vegetative height are higher in savanna areas. The burned areas were merged with the other land cover types such as urban, bare soil, etc., and new classes were assigned to vegetated areas indicated as burned: burned coniferous, burned deciduous, burned mixed forest, burned savanna, and burned grassland (Figure 4b). The land cover values were then used by WetSpass-M to assign land use parameters by means of a look-up table used in water balance calculation. A GIS model was built, which first extracted the fire count for the desired month, and then a raster of only pixels detecting fire was created using the Raster Calculator tool. Specific land cover types were extracted from the previously created yearly land cover data used in the non-burning WetSpass-M model. The extracted land cover classes were cropland, savanna, mixed forest, deciduous forest, coniferous forest, and grassland (Figure 4a). Each land cover area was then used as a mask to extract only those areas which experienced burning. For example, a new raster for savanna was created, only showing those areas in the savanna in which fire was detected for a particular month. Each of the newly created land cover rasters was then assigned a number corresponding to a WetSpass code used in the look-up table. Savanna and grassland were distinguished due to savanna containing more woody vegetation. The difference in vegetative cover affected the parameters used in the look-up table [Table S1]. For example, land cover roughness and vegetative height are higher in savanna areas. The burned areas were merged with the other land cover types such as urban, bare soil, etc., and new classes were assigned to vegetated areas indicated as burned: burned coniferous, burned deciduous, burned mixed forest, burned savanna, and burned grassland (Figure 4b). The land cover values were then used by WetSpass-M to assign land use parameters by means of a look-up table used in water balance calculation.  For this study, additional land cover types were included in the land use look-up table to account for burning. For example, LAI (the quantity of leaf coverage compared to ground surface area) changes in the look-up table parameters were calculated using the same method employed by [46], with modification of the equation: For this study, additional land cover types were included in the land use look-up table to account for burning. For example, LAI (the quantity of leaf coverage compared to ground surface area) changes in the look-up table parameters were calculated using the same method employed by [46], with modification of the equation: where VC is the burned vegetation cover, VC U is the original (unburned) vegetation cover, FBA is the fractional burned area, and SR is the survival rate. Relevant land cover parameters include LAI, vegetated area, burned area, impervious area, open water area, root depth, minimum stomata opening, vegetative height, Manning coefficient, land factor, and aerodynamic resistance. LAI is a quantity, measuring the leaf area per unit of ground surface area. LAI modifies the amount of water from precipitation and heat fluxes, affecting ET and runoff. Minimum stomata opening refers to the minimum size of stomata for that land cover type. Stomata, which are the pores on leaves used for gas exchange, have been shown to change size with temperature and depending on water availability [57].
Burning also reduces canopy height as land is cleared of vegetation, and thus increases wind speed. Increased wind speed reduces sensible heat, increasing leaf temperature and stomatal resistance, thus increasing transpiration. Conversely, decreased vegetative height tends to decrease aerodynamic resistance, also known as drag. Aerodynamic resistance is calculated for the WetSpass-M table using: where r a is the aerodynamic resistance, K is the von Karman constant (0.41), U a is the wind speed at elevation Z a , Z d is the zero-displacement elevation, and Z 0 is the aerodynamic roughness height of the surface [35]. Decreased aerodynamic resistance has the effect of decreasing ET [58].
Burning within the dominant land cover types showed grassland to have the highest fire detection, with an average of 5-20% of grassland showing fire during the dry season. Grassland fires may be due to controlled burning or wildfires, but that they make up such a large portion of the catchment is significant. Peak burning appears to occur from November to January for all classes, with little to no burning occurring from May to July ( Figure 5). REVIEW 9 where ra is the aerodynamic resistance, K is the von Karman constant (0.41), Ua is the wind speed at elevation Za, Zd is the zero-displacement elevation, and Z0 is the aerodynamic roughness height of the surface [35]. Decreased aerodynamic resistance has the effect of decreasing ET [58].

Earth2021, 2, FOR PEER
Analysis of land use model inputs for the CLC show the dominant classes to be grassland (74.88%), cropland (21.81%), and savanna (1.90%). Coniferous, deciduous, and mixed forest combined only make up 0.12%, with wetland covering 0.2%. Non-vegetative cover are open water (0.28%), urban build-up (0.18%), and bare soil (0.93%) (Figure 4b). Burning within the dominant land cover types showed grassland to have the highest fire detection, with an average of 5-20% of grassland showing fire during the dry season. Grassland fires may be due to controlled burning or wildfires, but that they make up such a large portion of the catchment is significant. Peak burning appears to occur from November to January for all classes, with little to no burning occurring from May to July ( Figure 5). Temperature data were extracted from the MOD11C3 Version 5 monthly daytime CMG land surface temperature product from MODIS Terra. MOD11C3 is derived from the MOD11C1 daily global product, taking daily clear sky land surface temperature values and averaging them for one month at a 0.05 degree resolution (~5.56 km). Validation of MODIS land surface temperature has been well researched [59][60][61]. Hulley and Hook [62] showed the high degree of MODllC1 accuracy in Africa by comparing the data to in situ measurements taken in Namibia, Botswana, and South Africa. MODllC1 Version 5 showed a mere 1.93% combined mean temperature difference from ground measurements. Analysis of land surface temperature showed a decrease in temperature during the wet season due to the cooling effect of precipitation over land.
Global climate models and datasets, such as the NOAA Global Forecast System (GFS) or the NASA Modern-Era Retrospective analysis for Research and Applications (MERRA), provide the wind speed at various atmospheric layers. Unfortunately, GFS has a relatively low resolution of 28 km, which is far too coarse for this study. Mesoscale models such as the Weather Research and Forecasting (WRF) model could be configured for the study area. However, this was beyond the scope of the study and, even then, such models have difficulty accounting for the complex terrain when calculating near-surface winds at a desired resolution [63,64]. For these reasons, station data were used for the wind speed, despite issues with data continuity. NOAA does not provide monthly wind speed averages, therefore daily records from NCEI Global Summaries of the Day were used and averaged monthly. The average wind speed values were also converted from 0.1 knots to m/s. Once the wind speed data were averaged, they were imported into ArcGIS as a shapefile and converted to rasters using kriging [65]. The data were then processed and converted into the proper format for input into WetSpass-M.
PET is the maximum evaporation and transpiration that could occur, assuming sufficient water. As noted earlier, burning has a significant impact on surface albedo change. This change in albedo affects the energy balance and therefore should not be overlooked when calculating PET. Two PET datasets were created for use with the two WetSpass-M models. Albedo values considering burned and non-burned vegetative classes were derived from [16], and those of non-vegetative areas were taken from Yu and Lu [66]. The spatial distribution of albedo was assigned using the modified MDC12Q1 (Type1) yearly land cover data showing burned areas. To consider the effect of albedo change on PET, a simplified Penman equation was used, following the method outlined by Valiantzas [67]: where E PEN is the Penman potential evapotranspiration, α is the albedo, R s is the solar radiation, T is the surface temperature, R A is the extraterrestrial radiation, RH is the relative humidity, a u is a wind function coefficient of 1, and u is the wind speed. The wind coefficient of 1 was originally used by Penman [68]. No solar radiation information was available from regional station data; therefore, it was estimated from measured sunshine hours using the equation: where n is the measured bright sunshine hours per day, and N is the maximum possible duration of daylight which is derived from the latitude of the site and the number of the Julian months [69,70]. R A is calculated using: where d r is the relative distance between the sun and the earth, φ is the latitude in radians, and δ is the solar declination. The maximum possible duration of daylight was calculated using the equation: where ω S is the sunset hour angle in radians. δ is the solar declination and is determined using: where J is the Julian day corresponding to the respective month [69,70]. Relative humidity was calculated using the equation: where e • (T) is the saturation vapor pressure (kPa) and e a is the actual vapor pressure. Climate station data from NCEI Global Summaries of the Month were used to calculate solar radiation, extraterrestrial radiation, and relative humidity; these values were used to create raster datasets. The data were imported into GIS and a model was used to calculate the two PET datasets and export them for use as WetSpass-M inputs. PET was validated by comparing the calculated PET to the Consultative Group for International Agricultural Research (CGIAR-CSI) PET, which averaged PET observations from 1950 to 2000 [71,72]. Calculated PET with albedo matched well against the CGIAR PET [ Figure S1].
Elevation from a digital elevation model (DEM) was provided by the Shuttle Radar Topography Mission (SRTM). The Shuttle Radar Topography Mission (SRTM) is an international project coordinated by the National Imagery and Mapping Agency (NIMA) and NASA [73]. The primary objective of the mission was the acquisition of a complete, high-resolution, digital topographic database of the Earth. SRTM, which flew onboard the Space Shuttle Endeavour, consisted of a specially modified radar system and had an 11-day mission in February 2000. The digital topographic map products were designed to meet Interferometric Terrain Height Data (ITHD) specifications, with 30 m × 30 m spatial sampling and absolute vertical height accuracy (90% linear error) of 16 m. The absolute horizontal accuracy (90% circular error) was 20 m. The SRTM resolution used in this study was 90 m × 90 m [74]. The slope was calculated using the elevation input by the WetSpass-M program.
Finally, a groundwater raster dataset was generated using 2009 measurements taken from a collaborative study in the Lake Chad Basin by the University of Missouri, Kansas City and NASA Goddard Space Flight Center and [75]. The data were imported into ArcGIS and converted to shapefile points and converted to raster using kriging. These same data were used by Babamaaji [6].

Validation
We performed model validation by comparing runoff outputs from WetSpass-M to the Logone river discharge measurements acquired at the Bongor station from 2003 to 2007 [76]. A Nash-Sutcliffe efficiency (NSE) value of 0.57 was found for the datasets, indicating that the model did well in calculating hydrological parameters (Figure S2a), where a value between 0 and 1 is considered good with 1 being the perfect match [77]. An r 2 value of 0.68 ( Figure S2b) found between the measured and calculated value was satisfactory.

Analysis
Analysis of the model results considering albedo change from burning showed the majority of precipitation being distributed to ET and recharge. The majority of available water from precipitation was lost to ET (45.19%) and groundwater recharge (45.35%), leaving a smaller portion available for surface runoff. We determined the burning impact by subtracting the monthly average outputs without burning from those with burning for 2003 to 2011.
To obtain a more detailed understanding of fire effect in the CLC, an investigation of the effect of fire by land cover type was performed. Particular attention was paid to the three primary land cover classes in the study area: grassland, savanna, and cropland. Average wet season ET decreased 6.44 + 03 mm per month, and decreased 2.01 + 02 mm per month for the dry season. This is in agreement with other studies [78,79], which found that ET was especially reduced by burning during the early monsoonal season. Increases in ET, recharge, and runoff due to fire were observed in the dry season and can be attributed primarily to grassland burning (Figure 6a,b). Considering burning, grassland showed an increase in ET and recharge during the dry season, but showed a decreasing pattern in the wet season. These parameters have increased runoff during the wet season (Figure 6b).
Considering no burning, there was more ET observed during the wet season in the savanna, whereas there was less ET, more runoff, and more recharge when considering burning (Figure 6c). Cropland showed more runoff, but less recharge and ET during the wet season (Figure 6d). To take the amount of monthly fire distribution into account, MODIS MOD14/MYD14 fire count data were compared to water cycle indicators using general linear least squares regression analysis (Figure 7). Correlations between annual burning changes to water balance parameters were negative. This matches findings by Ichoku [21], and supports the idea that with increased severity of fire, there is a severe decrease in hydrologic parameters including precipitation.
To minimize the effects of quantitative biases and uncertainties, correlations between fire count and water cycle parameters were calculated for the dry season and wet season. Both seasons showed a consistent negative correlation between fire count and water cycle parameters. The dry season findings differed from what was found in the Ichoku [21] study. However, it should be noted that the study area of their entire North Central block included significant wetland converted to agriculture, especially near the LCB. The CLC has very little wetland, making up a mere 0.2% of total land coverage; the majority of wetlands in the LCB lie to the northwest of this area. Therefore, it is reasonable to presume that changes due to wetland conversion to cropland mentioned by Ichoku [21] would not be observed, because the land cover is different in this regard. Results from a scatter plot of fire count against precipitation for the CLC show that burning has an inverse relationship with precipitation [ Figure 8]. The results look remarkably similar to those found by To take the amount of monthly fire distribution into account, MODIS MOD14/MYD14 fire count data were compared to water cycle indicators using general linear least squares regression analysis (Figure 7). Correlations between annual burning changes to water balance parameters were negative. This matches findings by Ichoku [21], and supports the idea that with increased severity of fire, there is a severe decrease in hydrologic parameters including precipitation.   To minimize the effects of quantitative biases and uncertainties, correlations between fire count and water cycle parameters were calculated for the dry season and wet season. Both seasons showed a consistent negative correlation between fire count and water cycle parameters. The dry season findings differed from what was found in the Ichoku [21] study. However, it should be noted that the study area of their entire North Central block included significant wetland converted to agriculture, especially near the LCB. The CLC has very little wetland, making up a mere 0.2% of total land coverage; the majority of wetlands in the LCB lie to the northwest of this area. Therefore, it is reasonable to presume that changes due to wetland conversion to cropland mentioned by Ichoku [21] would not be observed, because the land cover is different in this regard. Results from a scatter plot of fire count against precipitation for the CLC show that burning has an inverse relationship with precipitation [ Figure 8]. The results look remarkably similar to those found by Ichoku [21], and further emphasize the inverse relationship between burning and monsoonal rainfall.

Discussion
When comparing model simulations with and without fire activity considerations, a clear pattern of decreased ET was observed in the burning model. However, lower ET may limit the amount of water available for precipitation and could have significant meteorological impacts on agriculture, although more study in this area is needed. Burning seems to also have impact on the primary land cover types of the CLC, which, in turn, significantly impacted the hydrological parameters. Analysis of the burning scenario showed that burning grassland, which comprises almost 75% of Chari-Logone total land cover, causes increased ET and runoff during the winter months. Recent studies have demonstrated an increasing trend in the LCB of converting shrubland, grassland, and wetlands to cropland. This change from grassland to cropland has the potential of decreasing water availability for water bodies during the winter.
It is known that surface runoff tends to increase on burned bare surfaces after fire due to a decrease in infiltration into soils caused by a reduction in vegetation-induced roughness and an increase in the water repellency of soils [25,26]. Vadlonga et al. [80] showed that plots with increased prescribed burning had a tendency to decrease storage capacity with no significant increase in surface runoff. However, Onda et al. [81] found that ash layers from burning clogged the preferential flow paths of rainwater, leading to a significant increase in surface runoff. With an increase in accumulated ash towards the end of the dry season burning regimen, it is reasonable to expect that precipitation would likewise not infiltrate well into the soil instead running off the surface and contributing to river flow. Therefore, fires can contribute to the intra-annual irregularity of river flow at the end of the dry season and even during the wet season. Figure 6 shows the increase in runoff with burning from all three classes at the end of the dry season and the beginning of the wet season. To fully address the intra-annual irregularity of river flow, more studies are needed, using river flow models with more discharge data.
The ability of WetSpass-M to approximate discharge found in situ shows great promise for using it to calculate hydrological parameters despite a severe lack of data. By the use of satellite data in lieu of meteorological station records, temporal and spatial limitations may be mitigated. Development of a system to incorporate albedo into PET calculations at a relatively high resolution may prove useful for future studies examining the impact of burning on albedo and land cover type, as well as on the hydrological cycle.

Conclusions
Burning and no-burning models were used in WetSpass-M to determine the effect of burning on land classes. The effect of burning varied by vegetative type and whether the burning occurred during the dry or wet season in the LCB. A negative correlation of fire to water cycle was clearly observed. Water balance is affected by land cover, with each land cover responding differently to burning. If the pattern of converting shrubland to cropland continues, the primary difference will be less groundwater recharge. Grassland constitutes the majority of land cover in the CLC; therefore, it is worth noting that burning increases ET and recharge during the dry season but decreases ET and recharge during the wet season. These parameters cause an increase in runoff during the wet season. This research provides insight into how burning affects different land use types within a critical region of the LCB, which provides a significant amount of water to Lake Chad. All vegetative classes in burning scenarios showed a decrease in ET during the wet season.
Further research detailing the groundwater recharge along the Chari and Logone rivers would be beneficial to understanding the hydrology of the CLC in relation to fire, but more detailed field data for rainfall and wind are probably needed, as well as groundwater measurements. A greater distribution of these measurements increasing in number near the confluence of the Chari and Logone rivers and near their combined discharge to Lake Chad would be beneficial. If more data for burning to non-burning albedo comparisons become available, a longer study period could allow for better analysis of the effect of land cover change on burning in the CLC.

Data Availability Statement:
The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest:
The authors declare no conflict of interest.