1. Introduction
In Africa, floods and droughts have become a serious issue. The number of flood events is increasing considerably year-to-year, and the occurrence of drought is causing both substantial economic damage and harm to human health (source: Munich Re NatCatSERVICE:
https://www.iii.org/graph-archive/96134 accesed on 30 September 2021 and [
1]). Furthermore, it is predicted that climate change will cause severe floods and droughts to occur more frequently in the future in land areas within the monsoon domains of North Africa [
2]. Additionally, the proportion of the population with access to at least basic drinking water services in 2015 was much lower in sub-Saharan Africa than in the other regions of the world, and 58% of the regional population had no alternative to collecting untreated and often contaminated drinking water directly from surface water sources [
3]. Overall, the level of development of basic sanitation infrastructure is <50% in almost all countries within this region [
3]. Moreover, this region was the only region during 1990–2013, that registered an increase in the absolute number of people living in extreme poverty [
3]. In these circumstances, the occurrence of droughts or floods can be highly detrimental to food security. Africa has an issue with water-related disasters (particularly flood and drought disasters) and their consequences regarding socioeconomic development. Resolving this dire situation will require development of a system for data integration, information fusion, synthesis, information sharing and communication promotion. With such a system, it would be possible to ensure the maximum use of data and information from observation, monitoring, prediction, and socioeconomic surveys and statistics, which can assist African countries in overcoming such problems. In West Africa, the impact of agricultural droughts is becoming increasingly evident because of the Charny effect that links surface albedo and precipitation [
4], i.e., an increase in albedo due to reduced vegetation coverage causes precipitation to decrease, which results in a further increase in albedo because more of the land surface will be exposed owing to diminished growth of vegetation.
In West Africa, many researchers have been studying droughts, and the findings of some have identified the importance of (1) the contrast between the pre-monsoon and peak monsoon seasons, (2) two preferred modes of interannual variability (a latitudinal displacement of the tropical rain belt and changes in its intensity, and (3) the tropical easterly jet [
5]. Moreover, other research identified notable trends of decrease in rainfall in the Sahel region (10–20° N, 18°W–20° E) from the late 1950s to the late 1980s. Although, Sahel rainfall recovered somewhat through to 2003, drought conditions within the region did not end [
6]. These earlier approaches to drought assessment focused on the investigation of monsoon and rainfall trends using conventional drought indexes. Drought prediction over West Africa using the Standardized Precipitation Evapotranspiration Index (SPEI) and the Standardized Precipitation Index (SPI) has been implemented under the RCP4.5 and RCP8.5 scenarios [
7]. In a study of Niger River Basin in West Africa, the performances of three drought indexes, i.e., the Standardized Rainfall Anomaly Index, the Bhalme and Mooley Drought Index, and SPI, were evaluated and compared [
8]. Additionally, the question of how rising global temperatures might affect the spatial pattern of rainfall and the resultant droughts in West Africa was also investigated. Furthermore, precipitation and potential evapotranspiration variables have been simulated using the Rossby Centre RCA4 regional atmospheric model driven by 10 global climate levels under the RCP8.5 scenario (CanESM2, CNRM-CM5, CSIRO-Mk3, EC-EARTH-r12, GFDL-ESM2M, HadGEM2-ES, IPSL-CM5A-MR, MIROC5, MPI-ESM-LR, and NorESM1-M) [
9]. This approach to drought assessment is based on the use of atmospheric and land surface models.
The mainstream approach to the assessment of drought in West Africa has concentrated on monsoons and rainfall trends [
5,
6] using conventional drought indexes [
7,
8], e.g., the SPI, SPEI and the self-calibrating Palmer Drought Severity Index. Although drought assessment using conventional drought indexes can be effective, the emphasis of such an approach is placed on precipitation, with limited consideration of land surface hydrology and energy circulation. In studies using satellite remote sensing, the vegetation condition (e.g., the normalized difference vegetation index (NDVI), leaf area index (LAI), and vegetation optical depth (VOD)) is monitored using visible and near-infrared sensors such as the Moderate Resolution Imaging Spectroradiometer. Furhermore, the near-surface soil moisture content can be monitored using passive and active microwave sensors such as the Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E), Advanced Microwave Scanning Radiometer 2 (AMSR2), the Soil Moisture Active Passive satellite and the Soil Moisture and Ocean Salinity satellite. Thus, although satellite remote sensing cannot be used to assess the root-zone soil moisture content, which is important for vegetation growth dynamics, it can be used to evaluate the near-surface soil moisture content. Therefore, water absorption from the root-zone layer and vegetation growth dynamics are not considered in the NDVI, LAI and VOD. The Global Land Data Assimilation System (GLDAS) [
10] is the land surface data assimilation system integrated between a land surface model and a data assimilation scheme, in which the skin temperature is assimilated. In the GLDAS, not only the near-surface soil moisture content but also the root-zone soil moisture content is calculated by assimilating the skin temperature which is a variable of land surface hydrology and energy circulation used in calculating the flux by the land surface model. For vegetation, NDVI and LAI are used on the basis of visible and near infrared remote sensing. The Coupled Land and Vegetation Data Assimilation System (CLVDAS) [
11,
12,
13] integrates passive microwave remote sensing techniques, a land surface model, a dynamic vegetation model (DVM) and a data assimilation scheme. CLVDAS assimilates microwave brightness temperatures (6.925 GHz; vertical polarization, 6.925 GHz; horizental polarization, 10.7 GHz; vertical polarization and 10.7 GHz; horizental polarization) that represent the product between the skin temperature and microwave emissivity. Thus, CLVDAS can estimate not only the optimized near-surface soil moisture content but also the optimized root-zone soil moisture content through data assimilation, because the microwave brightness temperature is sensitive to moisture. Additionally, the integrated DVM can be used to estimate the optimized LAI and vegetation water content, which can provide evaluation of vegetation growth dynamics, through data assimilation of microwave brightness temperatures. Thus, the gap between conventional study and CLVDAS-based studies is as follows: (1) evaluation of optimized root-zone soil moisture content through data assimilation, and (2) evaluation of optimized LAI and vegetation water content based on vegetation growth dynamics.
To fill this gap, land surface hydrology and energy circulation were evaluated using CLVDAS [
11,
12,
13] in this study. Because this approach provides the possibility of drought monitoring and application to agricultural support [
14], we evaluated the relationship between the pearl millet yield as a major crop and the simulated vegetation water content as a drought index, and analyzed the applicability of the approach to the assessment of agricultural droughts in the Sahel-inland region of West Africa during 2003–2018.
4. Results and Discussion
In passive microwave remote sensing, the low-frequency band with the longest wavelength is used by AMSR-E and AMSR2 because the microwaves emitted from the soil must be scanned through the atmosphere and vegetation. However, a passive microwave sensor can generally detect only the near-surface soil moisture content because even microwaves in the low-frequency band are absorbed by soil moisture [
23]. This is the reason why the AMSR-E and AMSR2 soil moisture product targets only the near-surface soil moisture content [
23]. To overcome this shortcoming, as described in
Section 1, the following processing is implemented in CLVDAS. Eco-HydroSiB evaluates the ecohydrological water cycle, which expresses the penetration of precipitation, water storage in the root-zone, and water absorption by roots and vegetation growth, and simulates the ecohydrological variables (e.g., soil moisture profile, evapotranspiration, and biomass). Furthermore, the RTM calculates the microwave brightness temperatures emitted from the land surface using the simulated ecohydrological variables. Additionally, changes in the ecohydrological variables and assimilation of the microwave brightness temperatures are repeated until the difference between the simulated and satellite-observed microwave brightness temperatures of the ground surface is minimized. Hence, an accurate estimation of the ecohydrological water cycle is derived using this process. This represents the major advantage of this system. Hitherto, the CLVDAS outputs of estimated soil moisture content and LAI were validated by comparison with the following observations: ground-based observed soil moisture content and LAI observed at the Yanco Flux Tower site located in New South Wales (Australia) [
12], ground-based observed soil moisture content from the African Monsoon Multidisciplinary Analyses, Vaira Ranch (USA), Bayantsagaan (Mongolia) [
11,
13], and the Moderate Resolution Imaging Spectroradiometer LAI in West Africa and Northeast Brazil [
14]. The following estimation accuracy was achieved: (1) root mean square error (RMSE) of 0.05 m
3/m
3 or less and bias of 0.045 m
3/m
3 or less in terms of soil moisture content, and (2) RMSE of 0.12 m
2/m
2 or less and bias of 0.14 m
2/m
2 or less in terms of LAI.
For the period from 1 January 2003to 31 December 2018, the near-surface soil moisture content (m
3/m
3), root-zone soil moisture content (m
3/m
3), and vegetation water content (m
3/m
3) were simulated using CLVDAS and used to create a gridded dataset (temporal resolution: daily, spatial resolution: 0.25° × 0.25°).
Figure 3 shows the spatial distribution of the averaged ecohydrological variable in September. In the Sahel-inland region, vegetation water content has shown a trend of decrease from the north since 2011 (
Figure 3c). This trend has also shown in the root-zone soil moisture content (
Figure 3b) but not so clearly in the near-surface soil moisture content (
Figure 3a). Furthermore,
for each ecohydrological variable was calculated for each grid using Equation (2).
Figure 4 shows the spatial distribution of each normalized index of the averaged ecohydrological variables in September. Although variation in the vegetation water content itself is slightly unclear (
Figure 3c), its normalized index (
) based on the z-score theory clearly shows a decrease since 2011 (
Figure 4c). Furthermore, this trend has also shown in the near-surface and root-zone soil moisture content (
Figure 3a,b).
Figure 5 shows the annual variation in the normalized index (
) based on the z-score theory for vegetation of water content in the Sahel-inland region, calculated using Equation (2) (September average in the agricultural drought assessment period from 2003 to 2018, regional spatial average). We calculated the annual total yield of pearl millet, in the Sahel-inland region consisting of Chad, Niger, Nigeria, Benin, Burkina Faso, and Mali from 2003 to 2018 (
Figure 6a), which revealed that the annual yield since 2011 has been approximately half that in 2008. Furthermore, we calculated the normalized index (
) based on the z-score theory of the total yields of pearl millet using Equation (2) (
Figure 6b).
Figure 7 shows the annual variation in
for both pearl millet yield and vegetation water content in the Sahel-inland region (temporal average in September and regional spatial average). The
values for the pearl millet yield are negative during 2003–2004, whereas the concurrent
values for vegetation water content are positive, and the difference between the two sets of
values is large. We investigated external factors other than droughts by considering previous studies [
24,
25] and the FAO Locust watch (
http://www.fao.org/ag/locusts/en/archives/briefs/index.html accessed on 30 September 2021), which revealed that serious locust outbreaks occurred during 2003–2004 in West Africa. The number of days (
) in each year (2003–2018) was counted when more than 10 locust plagues occurred in the Sahel-inland region using the FAO Locust watch. The average (
) and standard deviation (
) were calculated using these the numbers of days (
) in each year. Using
,
, and
, the normalized index (
) based on the z-score theory is calculated by using Equation (2) (brown line). The
values for locust outbreaks (brown line in
Figure 7; calculated as described in
Section 3.3) in 2003 and 2004 are 1.31 and 3.37, respectively; although, all
values after 2005 are negative. We recognize that crop yields of the Sahel-inland region were likely to be adversely affected by the external impact of locust plagues in 2003–2004. Therefore, we assessed agricultural drought in the period from 2005 to 2018. For this period, the RMSE between the
for pearl millet yield (green line in
Figure 7a) and the
for vegetation water content (red line in
Figure 7a) was 0.16 and the correlation coefficient was 0.89, indicating a strong agreement. As an aside, the RSME and correlation coefficient values when including the period 2003–2004 were 0.25 and 0.73, respectively).
The variation in
for precipitation and the simulated ecohydrological variables (near-surface soil moisture content, root-zone soil moisture content and vegetation water content) were investigated for the Sahel-inland region during the agricultural drought assessment period (2005–2018) (
Figure 8 and
Figure 9). In the first half of the agricultural drought assessment period (
Figure 8), the
of vegetation water content was mostly positive, except in 2005, 2008 and 2009, which indicates that the effect of drought on vegetation water content is small. We also found that each peak had a time lag when focusing on the negative peak for precipitation, soil moisture content, and vegetation water content in 2005, 2008, and 2009 (yellow marks and lines in
Figure 8). This is attributable to the following process: (i) the land surface soon dries because of the shortage of precipitation, (ii) the root-zone soil moisture decreases after a further amount of time because water is not supplied from the land surface, and (iii) vegetation water content declines after an even further amount of time because of the lack of root-zone soil moisture available for absorption by roots. In the second half of the agricultural drought assessment period (
Figure 9), the vegetation water content gradually became negative after 2014. As in the first half of the agricultural drought assessment period (
Figure 8), the shortage of precipitation gradually propagated to the lack of vegetation water content (yellow marks and lines in
Figure 9). Although the
values of precipitation became positive in April 2016, negative values remained in the growth and early harvest seasons (May–October), leading to negative
values of soil moisture and vegetation water content simultaneously (green line in
Figure 9). The
values of precipitation became positive in the second half of December 2016 owing to the occurrence of heavy rainfall; however, they reverted to negative values in the beginning of January 2017. The surface soil moisture content showed the same behavior. In contrast, the root-zone soil moisture content was positive in the period from January–April because the water associated with the heavy rainfall was stored in the rootzone. The vegetation grew by absorbing the stored root-zone soil moisture, as indicated by the positive peak of vegetation water content at the beginning of April (blue line in
Figure 9). This indicates that storage of root-zone soil moisture is important for vegetation growth, and that CLVDAS can evaluate such a mechanism. In 2017, soil moisture was not stored in the near-surface soil because precipitation amounts were low. Following subsequent rainfall events, the
values of near-surface soil moisture recovered (became positive) in March 2018; however, the
values of root-zone soil moisture and vegetation water content did not recover (remained negative). Thus, both the root-zone soil moisture content and vegetation water content remained negative in the long term because they have a long retention period of the memory of past water shortages (red line in
Figure 9). By investigating the daily variation in precipitation, soil moisture content, and vegetation water content using CLVDAS, we were able to evaluate the following land surface hydrological water cycle and vegetation growth dynamics mechanism. (i) Reduced precipitation causes aridity of the land surface, which affects the condition of the root-zone soil moisture. (ii) Plants cannot retain sufficient water within their structures because of the lack of soil moisture available for absorption in the root-zone layer. (iii) Near-surface soil moisture content can change rapidly in response to temporary rainfall events. In contrast, both the root-zone soil moisture and the vegetation water content tend to remain in long-term drought conditions because they have a long retention period of the memory of past water shortage. The major finding of this study was establishing that CLVDAS can be used to evaluate the hydrological water cycle (penetration of precipitation to the root-zone soil layer and its absorption by roots) and vegetation growth dynamics. Additionally, we confirmed that vegetation water content output by CLVDAS can be used to assess agricultural drought through comparison with major crop yields and investigation of external factors in the target region.