On the Influence of Vegetation Cover Changes and Vegetation-Runoff Systems on the Simulated Summer Potential Evapotranspiration of Tropical Africa Using RegCM4

The community land model version 4.5 provides two ways for treating the vegetation cover changes (a static versus an interactive) and two runoff schemes for tracking the soil moisture changes. In this study, we examined the sensitivity of the simulated boreal summer potential evapotranspiration (PET) to the aforementioned options using a regional climate model. Three different experiments with each one covering 16 years have been performed. The two runoff schemes were designated as SIMTOP (TOP) and variable infiltration capacity (VIC). Both runoff schemes were coupled to the carbon–nitrogen (CN) module, thus the vegetation status can be influenced by soil moisture changes. Results show that vegetation cover changes alone affect considerably the simulated 2-m mean air temperature (T2M). However, they do not affect the global incident solar radiation (RSDS) and PET. Conversely to the vegetation cover changes alone, the vegetation-runoff systems affect both the T2M and RSDS. Therefore, they considerably affect the simulated PET. Also, the CN-VIC overestimates the PET more than the CN-TOP compared to the Climatic Research Unit observational dataset. In comparison with the static vegetation case and CN-VIC, the CN-TOP shows the least bias of the simulated PET. Overall, our results show that the vegetation-runoff system is relevant in constraining the PET, though the CN-TOP can be recommended for future studies concerning the PET of tropical Africa.


Introduction
Among different processes within the terrestrial hydrological cycle, evapotranspiration (ET) plays an important role in determining the global water balance and linking the global water, energy and carbon cycles (Pielke et al. 2002;Bonan 2008;Yang 2015). Indeed, it controls the retroactions between land and atmosphere through energy fluxes (Mamadou et al. 2014(Mamadou et al. , 2016 and thus sustains precipitation (Taylor et al. 2012;Nicholson et al. 2013;Koster et al. 2004). The 'potential evapotranspiration' (PET) however, can be seen as the amount of the water vapor that would evaporate if a sufficient water source is available. Estimation of the PET is valuable for many operational applications especially in agriculture and hydrology. Indeed, PET is required for estimating the crop water requirements, supporting irrigation scheduling, drought management, and climate change studies. Furthermore, it is an important input to a variety of hydrological model applications 1 3 Published in partnership with CECCR at King Abdulaziz University (Andréassian et al. 2004;Vicente-Serrano et al. 2014). Numerous formulations of the PET exist in literature. However, the most used and called 'the Penman-Monteith method (PM)', has been proposed by Allen et al. (1994). Indeed, they defined the PET for a hypothetical reference crop of 0.12 m height, bulk surface resistance of 70 s m −1 , and albedo of 0.23 assuming active green grass. In addition, Allen et al. (1998) calculated the PET using a process-based approach. This approach needs many meteorological inputs including the 2-m air temperature (T2M), 2-m relative humidity, and near surface wind speed (see Eq. 1). Kjelgaard and Stokes (2001) calculated the PET (or λE as a source of latent heat flux) as follows: where: λE, R n , and G are expressed in W m −2 , VPD is the vapor pressure deficit (kPa), Δs is the slope of saturation vapor pressure curve (kPa °C −1 ) at 2-m air temperature, ρ a is the density of air (kg m -3 ), c p is the specific heat of the air (J Kg -1 °C −1 ), γ is the psychrometric constant (kPa °C −1 ), r a is aerodynamic resistant (s m -1 ), r s is the surface resistance to vapor transport (s m -1 ).
Besides, the PM method was recommended by the Food and Agriculture Organization of the United Nations (FAO) to accurately estimate the PET from climate data (Allen et al. 1998). Due to the numerous meteorological variables which are not often all available, alternative approaches that consider limited meteorological inputs are needed. Since the 1940s, several empirical methods have been developed. These methods can be grouped into five categories: (1) water budget (Guitjens 1982), (2) mass-transfer (Harbeck 1962), (3) combination (Penman 1948), (4) radiation (Priestley and Taylor 1972), and (5) temperature-based (Thornthwaite 1948). Among these, the temperature-based ones offer an interesting option to study the future projected PET. Hargreaves method (hereafter HG;Hargreaves and Allen 2003) is the best alternative method for quantifying the PET because of its ability to capture information relative to the land surface property, such as soil moisture and albedo through the diurnal temperature range (Almorox et al. 2015;Li et al. 2018). Note that, this approach has been widely applied and evaluated across different climate regions with good performance (Droogers and Allen 2002). According to the HG method, PET is computed as: where PET HG is the daily PET (unit: mm day −1 ); R a is the extra-terrestrial radiation for a given day and latitude; T max , PET HG = 0.0023 × R a × T max − T min 0.5 × T mean + 17.8 , T min , and T mean are the daily values of maximum, minimum, and mean air temperature at 2 m above the surface (°C) respectively. When the global incident solar radiation (RSDS) is readily available, Eq.
(2) can be rewritten as follow: R a /RSDS is expressed in units of mm day −1 to show how much energy is used to evaporate water (Allen et al. 1998). Performance of the HG method with respect to the other approaches has been highlighted in various studies. For instance, Er-Raki et al. (2010) found that the HG method is the best to calculate PET compared to the PM formula. Sperna et al. (2012) recommended calibrating the HG equation to calculate the PET compared to the Climatic Research Unit (CRU) gridded product. Likewise, Sheikh and Mohammadi (2013) showed the superior performance of the HG method over the other empirical methods, with the HG showing the least compared to the PM observations. By comparing the PET obtained by different methods (e.g. HG) using either station data or CRU gridded data, Potop and Boroneant 2014) showed that the HG method gives results closer to the CRU-PET data at the Chisinau station. Finally, Murat et al. (2017) found that the modified HG equation can give a relatively precise estimate of the PET compared to the PM observations. Recently, regional climate models (e.g. RegCM4) were used to study the influence of vegetation cover/fraction changes on the tropical African climate Yu et al. 2015;Erfanian et al. 2016). For instance, Yu et al. (2015) studied the future climate changes of West Africa using the RegCM4 model. They observed that the future precipitation changes are reversed as a result of changes in the vegetation cover. In addition, the RegCM4 model was used to examine the role of runoff parameterization in constraining the surface climate of tropical Africa (Anwar et al. 2019). These authors showed that the Variable Infiltration Capacity (shortly VIC; Liang et al. 1994) runoff scheme is better than the SIMTOP (hereafter TOP; Niu et al. 2005) runoff scheme for simulating the tropical African climate. Adeniyi (2019) utilized the RegCM4 model to examine the potential influence of the total solar irradiance (TSI) on the precipitation, temperature, evapotranspiration, water availability for plant, and heat stress over West Africa. In general, the RegCM4 showed a good capability to capture the diurnal precipitation trend, intra-annual variations of precipitation, and temperature over West Africa. Furthermore, it was found that increase (decrease) of TSI results in increase (decrease) of the aforementioned variables. The influence of the two runoff schemes (TOP and VIC) on the terrestrial carbon cycle of tropical Africa with this model revealed that the runoff scheme has a considerable influence on the terrestrial (3) PET HG = 0.0135 × RSDS × T mean + 17.8 .
Published in partnership with CECCR at King Abdulaziz University carbon cycle of tropical Africa (Anwar and Diallo 2021a). Also, the VIC runoff scheme shows superior performance over the TOP runoff scheme for simulating the net ecosystem exchange (NEE). Anwar and Diallo (2021b) showed that considering changes in the leaf area index (LAI), the RegCM4 model bias is amplified. However, the VIC runoff scheme still outperforms the TOP scheme relative to the static case. In a most recent study, Anwar and Diallo (2021c) found that the gross primary production (GPP) is more sensitive to the vegetation-runoff systems than when the vegetation cover changes alone. In addition, the GPP can be properly simulated when the VIC surface dataset is calibrated with respect to reanalysis products and/or in-situ observations. Finally, Anwar and Diallo (2021d) utilized the RegCM4 model to examine the potential influence of the new LAI formula on the surface energy balance and surface climate of tropical Africa. They found that utilizing the new LAI formula leads to a better performance (than using the old LAI formula) in simulating the surface energy balance and surface climate with respect to reanalysis products. While many studies have focused on the tropical African climate using different physical configurations, none of them focusing on PET and in particular on the influence of the vegetation cover changes alone and vegetation-runoff systems in simulating the boreal summer PET climatology. This paper aims to examine the sensitivity of the simulated boreal summer PET to the vegetation cover changes alone (as represented by LAI) and vegetation-runoff systems using the RegCM4 model. Section 2 describes the study area, experiment design, and validation data. Section 3 shows the results of the study. Section 4 provides the discussion, and while conclusions are presented in Sect. 5.

Model Description
The regional climate model version 4.5.0 (hereafter RegCM4; Giorgi et al. 2012) was used in the present study. RegCM4 is a limited area model and has been used for numerous studies including process studies (Sylla et al. 2015;Diallo et al. 2016;Kébé et al. 2020) and future climate projections (Diallo et al. 2016;Sawadogo et al. 2020;Ashfaq et al. 2020;among others). Performance of the RegCM4 model has been examined in tropical Africa as in Ogwang et al. (2016); Komkoua et al. (2017) andN'Datchoh et al. (2018). The selected physical configuration is documented in Anwar and Diallo (2021a;b, c, d). Recently, the Community Land Model version 4.5 (CLM45; Oleson et al. 2013) is coupled to the RegCM4. The CLM45 simulates various complex biogeophysical processes such as absorption, reflection, and transmittance of solar radiation, absorption, and emission of long wave radiation. Also, it includes complex biogeochemical processes such as photosynthesis and vegetation phenology.
In the CLM45 model, vegetation can be treated in three ways: derived from the Moderate Resolution Imaging Spectroradiometer (MODIS). This mode is known as Satellite Phenology (SP; Lawrence et al. 2011). 2. Prediction of vegetation cover changes (e.g. LAI) using a prognostic carbon-nitrogen model (CN; Thornton and Rosenbloom 2005) to simulate/predict the seasonal variation of vegetation structure with a prescribed vegetation demography. 3. Prediction of both vegetation cover and fractional percentage using the CN model combined with a dynamic vegetation module (DV; Levis et al. 2004).
In addition, the CLM45 model has two runoff schemes: 1. SIMTOP (TOP; Niu et al. 2005). In the TOP runoff scheme, both surface and sub-surface runoff generation are predicted based on knowledge of the fractional saturated area, which in turn is determined by topographic characteristics and soil moisture state of a grid cell. Moreover, the total runoff distributed over the land grid is constituted by surface, sub-surface runoff, and liquid runoff from different sources (glacier, wetlands, and lakes). 2. Variable infiltration capacity (VIC; Liang et al. 1994).
The VIC is a large-scale, semi-distributed hydrologic model. Basically, the VIC model requires the knowledge of four parameters (B inf , D smax , D s and W s ) to calculate the surface, sub-surface, runoff and related surface variables (e.g. energy balance). These four parameters are provided in a global surface dataset. B inf is a parameter that controls the shape of the soil moisture-holding capacity curve; D smax is the maximum sub-surface flow rate (in kg m −2 s −1 ); D s is a fraction of D smax . Finally, W s is a fraction of w m,bot (maximum soil moisture in kg m −2 ).
When the CN module is enabled, the LAI is predicted using the formula of Thornton and Zimmermann (2007) defined as follow: where m (m 2 gC −1 ) is a linear slope coefficient for each PFT, SLA o (m 2 gC −1 ) is the specific leaf area at the top of the Published in partnership with CECCR at King Abdulaziz University canopy, C L (gC m −2 ) is the leaf carbon content and L is the predicted leaf area index (m 2 m −2 ). In the CLM45 land model, LAI is predicted for both sunlit and sun-shaded sides as: where TLAI is the total projected leaf area index.

Experiment Design
The model domain was centered at latitude 5° and longitude 12° using a 60 km horizontal grid-spacing resolution, 18 vertical levels, and 160 horizontal grid points in both zonal and meridional directions (Fig. 1). Three simulations were conducted for 16 years spanning from 01 January 1995 to 31 December 2010. The first three years of each simulation were discarded to allow equilibration of the soil moisture (Steiner et al. 2009). The period from 01 January 1998 to 30 November 2007 was considered for further analysis. Such period was chosen to be consistent with the availability of the Solar Radiation Budget (SRB) reanalysis product. The SP-VIC, CN-TOP, and CN-VIC form the identical-triplet sets (5) TLAI = LAI sunlit + LAI sunshaded , of 16-year simulations, which have the same physical parameterization options and initial conditions, with either the vegetation status or the runoff scheme varying between the experiments (see Table 1 file). In all simulations, the PET was calculated using Eq. (3) (see Sect. 1). The boreal summer (June-July-August; JJA) season was selected because it is the rainy season of tropical Africa. Wang et al. (2015) showed that the RegCM4 performance in the West African domain is similar when driven either by the 1.5° × 1.5° resolution ERA-Interim reanalysis (EIN15; Dee et al. 2011) or the 2.5° × 2.5° National Centre of Environmental prediction/ National Centre of Atmospheric Research version 2 reanalysis (NCEP/NCAR2; Kanamitsu et al. 2002). In the present study, the NCEP/NCAR version 2 reanalysis data was selected as the initial and the lateral boundary conditions (LBCs), while the EIN15 was used as the prescribed Sea Surface Temperature (SST). A long-term spinning up of the CN was used to initialize the CN-TOP/CN-VIC simulations with an equilibrium state of the vegetation carbon and net ecosystem exchange as recommended by Fang et al. (2015) and Wang et al. (2015).

Validation Data
In this study, the Climatic Research Unit gridded Time Series version 4.05 (CRU TS v4.05; Harris et al. 2020; hereafter CRU) dataset was used to evaluate the RegCM4 model performance. The CRU dataset includes various variables: cloud cover, 2-m diurnal temperature range, frost day frequency, PET, precipitation, 2-m mean, maximum and minimum air temperature as well as vapor pressure. These variables are monthly averaged and integrated over the period 1901-2020. The CRU product is available on a 0.5° × 0.5° horizontal grid over land domains. The CRU-PET is considered as the best available reference PET data across the globe (Droogers and Allen 2002;Mitchell and Jones 2005;IPCC 2007;Sperna et al. 2012;Potop and Boroneant 2014). Also, Solar Radiation Budget (SRB) reanalysis product was used to evaluate the global incident solar radiation (RSDS). The SRB is a part of the World Climate Research Programme's (WCRP) Global Energy and Water Exchanges (GEWEX) framework, and it covers the period of July 1983 to November 2007.
SRB is constructed on a 1° × 1° global grid using satellite-derived cloud parameters and ozone fields, reanalysis meteorology, and a few other ancillary data sets. SRB product consists of downward and upward components of shortwave (SW) and longwave (LW) radiation as a major component of the energy exchanges between the atmosphere and land/ocean surfaces and thus affects surface temperature fields, sensible and latent heat fluxes, as well as energy and hydrological cycles. To facilitate the comparison between the model and observation, CRU gridded-observation and SRB reanalysis products were bi-linearly interpolated on the RegCM4-CLM45 horizontal grid. This bi-linear interpolation already used in some previous studies (Diallo et

Results
Equation (3) shows that PET depends on both the 2-m mean air temperature (T2M) and global incident solar radiation (RSDS). Therefore, it is necessary first to assess the simulated T2M and RSDS with respect to reanalysis products. For instance, Anwar (2021) examined the RegCM4 performance for simulating the T2M with respect to the CRU product. In this study, it was found that the CN-VIC overestimates the T2M more than the SP-VIC with respect to the CRU. This overestimation of T2M will likely affect the simulated PET (see Sect. 3.1). Furthermore, uncertainty of the RSDS was quantitatively evaluated with respect to the Solar Radiation Budget (SRB) reanalysis product ( Figure  S1 in supplements). From Figure S1, it can be observed that the RegCM4 overestimates the RSDS in the Sahara region by 20-60 W m −2 with respect to the SRB. Also, the CN-VIC reduces the RSDS negative bias by 20-30 W m −2 relative to the CN-TOP (which can explain the considerable difference between the two configurations in simulating the PET; see Sect. 3.2).

Effects of Vegetation Cover Changes
In this section, the influence of vegetation cover changes alone was examined. The SP-VIC was considered as the control simulation, while the CN-VIC was considered as the sensitivity simulation. Figure 2 shows the pattern of the spatial distribution for the following variables: 2-m relative humidity (RH2M; in %, top panels), T2M (in K, middle panels), and RSDS (in W m −2 , bottom panels). The statistical significance of the difference was calculated using a student t test with α = 0.05. The CN-VIC simulates lower RH2M than the SP-VIC by 15-25 % over all tropical Africa (Fig. 2a-c). This is an indication that the CN-VIC severely decreases the vegetation transpiration (QVEGT) relative to the SP-VIC particularly over the Congo basin as found in Anwar and Diallo (2021b) as a consequence of the simulated decrease of total projected leaf area index (TLAI) (particularly over the Congo basin). An immediate effect is, that the CN-VIC shows higher T2M in the range of 1-3 K compared to the SP-VIC (Fig. 2d-f). Regarding the RSDS, the difference between the two simulations is quite small (± 5 W m −2 ) (see Figs. 2 g-i). Figure 3 shows the spatial pattern of the simulated boreal summer PET with respect to the CRU. It is clear that that both simulations reproduce adequately the spatial pattern of the PET with respect to the CRU (Fig. 3a-c). CN-VIC and the significant difference (Diff) between the two simulations. The significant difference was calculated using student t test with α = 0.05 Published in partnership with CECCR at King Abdulaziz University However, the two simulations overestimate the PET by 0.5-2 mm day −1 particularly over the Northern Savanna and the Congo basin. Such biases of the PET can be partially attributed to the T2M warm biases (Anwar 2021). Qualitatively, the difference between the two simulations is less than 1 mm day −1 and is only significant over the Northern Savanna (Fig. 3d-f). Overall, our results show that the simulated summer PET doesn't change considerably between the two simulations. Therefore, considering vegetation cover changes alone does not cause a notable change in the simulated PET.

Effects of Vegetation-Runoff Systems
The CN-VIC simulates lower RH2M (by 10-30 %) relative to the CN-TOP (Fig. 4a-c). This is due to the simulated severe decrease of the TLAI and QVEGT in CN-VIC compared to the CN-TOP as reported in Anwar and Diallo (2021b). In addition, the CN-VIC shows a warmer T2M than the CN-TOP by 2-4 K (Fig. 4d-f). This is a consequence of a considerable overestimation of the sensible heat flux in the CN-VIC compared to the CN-TOP (Anwar and Diallo 2021b). Finally, the CN-VIC simulates higher RSDS than the CN-TOP by about 15-30 W m −2 (Fig. 4 g-i). Figure 5 shows the spatial pattern of the simulated boreal summer PET with respect to CRU. It can be seen that both Published in partnership with CECCR at King Abdulaziz University simulations can capture the spatial pattern of the PET with respect to the CRU (Fig. 5a-c). However, the difference between the two simulations is quite large. The CN-TOP overestimates PET by 1 mm day −1 over Northern and Southern Savanna regions and by 0.4 mm day −1 over the Congo basin, while the CN-VIC overestimates the PET by 1.2-2.4 mm day −1 over the entire tropical region. It is thus clear that, the CN-VIC overestimates the PET more than the CN-TOP by 0.8-1.4 mm day −1 (0.2-0.6 mm day −1 ) over the Congo basin (Northern Savanna) (Fig. 5d-f). The large difference between the two simulations can be linked to the combined changes of the T2M and RSDS. In addition, it can be understood that vegetation-runoff systems have a considerable influence on the simulated PET more than when vegetation cover changes alone. In comparison with the SP-VIC and the CN-VIC, the CN-TOP shows the TOP, CN-VIC and the significant difference (Diff) between the two simulations. The significant difference was calculated using student t test with α = 0.05 least bias with respect to the observations. However, this bias is still high particularly over both the Northern and Southern Savanna regions. It is worth mentioning that calibrating the HG method (Eq. 3) is necessary to ensure a low bias as suggested by previous studies (Jensen et al. 1997;Vanderlinden et al. 1999;Xu and Singh 2002;Dinpashoh 2006;Murat et al. 2017). Therefore, the constant "17.8" (see Eq. 3) has been adjusted to be "10.8" after conducting various tests. Equation (3) then takes a new form and is rewritten as follow: where the PET HGnew is the calculated new PET. Figure 6 shows the spatial pattern of the simulated PET before and after calibrating the HG equation. When Eq. (6) is used, (6) PET HGnew = 0.0135 × RSDS × T mean + 10.8 , RegCM4 simulates the PET with a lower bias than when Eq. (3) is used (Fig. 6a-c). For instance, RegCM4 simulates the PET HG with a bias in the range of 0.8-1.6 mm day −1 over Northern and Southern Savanna regions. However, the biases for the simulated PET HGnew are in the range of ± 0.4 mm day −1 (see Fig. 6d-e). Qualitatively, RegCM4 underestimates the PET HGnew more than the PET HG by 0.2-0.4 mm day −1 (0.8-1.2 mm day −1 ) over the Northern and Southern Savanna regions (the Congo basin) (Fig. 6d-f). To quantitatively evaluate the RegCM4 performance, we plotted (in Fig. 7) the Box-and-Whisker of the boreal summer potential evapotranspiration (PET; in mm day −1 ) before calibrating the HG equation (RegCM4_HG; in blue) and after calibrating the HG equation (RegCM4_HGnew; in red) with respect to the CRU (in green) for the three regions (Northern Savanna, Evergreen Forest and Southern Savanna). This visualization allows displaying measures of the central tendency of data and spreading free from the assumption of a normal distribution. It can be obviously seen that, the RegCM4_HGnew compares better with the CRU than RegCM4_HG particularly over the Evergreen Forest and Southern Savanna regions. Thus, the proposed calibration of the HG method performs better in simulating PET compared to CRU product.

Discussion
The PET is an essential variable for assessing hydrological and agricultural activities, as well as monitoring drought events. The FAO organization recommends the utilization of the PM method to calculate the PET; but this method requires many meteorological variables. Furthermore, the Averaged potential evapotranspiration (PET; in mm day −1 ) over the period 1998-2007 for the boreal summer season (June-July-August; JJA) as simulated by the coupled RegCM4-CLM45-CN-TOP comparing between the default HG equation (HG) and the calibrated HG equation (HGnew). In the first row, HG is the first panel on the left, HGnew is the second panel from the left and CRU data is on the right; the second row shows the Diff1 (HG minus CRU), Diff2 (HGnew minus CRU) and the difference between HGnew and HG. Significant changes are indicated in black dots using student t test with α = 0.05 uncertainty of the meteorological variables can amplify the error of the calculated PET (since the PM equation is not linear). Therefore, there was an urgent need to compute the PET using a simple empirical method. Among the temperature-based methods, the HG method was used to compute the PET. We compute PET from the variables simulated by RegCM4, and examine the influence of different physical configurations on the simulated boreal summer PET. Finally, we calibrate the HG method with respect to the CRU product using the RegCM outputs. The results of this study can be summarized as follows: 1. The influence of vegetation cover changes alone induces a little impact on the simulated PET. Larger impact is shown in the vegetation-runoff systems changes. 2. In comparison with the SP-VIC and CN-VIC, the CN-TOP shows the least bias with respect to the CRU product.

Fig. 7
Box-and-Whisker plot to show the boreal summer potential evapotranspiration (PET; in mm day −1 ) spread before calibrating the HG equation (RegCM4_HG; in blue) and after calibrating the HG equation (RegCM4_HGnew; in red) with respect to the CRU gridded product (in green) for the three sub-area averages highlighted in Fig. 1 1 3 Published in partnership with CECCR at King Abdulaziz University 3. Using the calibrated HG equation, the regional coupled model simulation RegCM4-CLM45-CN-TOP shows a very good performance for simulating the boreal summer PET with respect to the CRU gridded-based product.

Conclusions
We conducted three 16-year simulations to examine the potential influence of changes in the vegetation cover alone and vegetation-runoff systems on the simulated boreal summer PET of tropical Africa. We observed that vegetation-runoff systems affect the PET more than when vegetation cover changes are considered alone. When the calibrated HG equation is used, the regional coupled RegCM4-CLM45-CN-TOP model can be recommended for future studies concerning the PET and the Standardized Precipitation-Evapotranspiration Index (SPEI) of tropical Africa. In our future work, we will address the following points: 1. Incorporating the role of the nitrogen deposition, carbon dioxide (CO 2 ), land use land cover changes (LULCC), crops, and aerosols with enabling the dynamic vegetation (DV) using a high-resolution grid-spacing and at a domain similar to the one utilized within the framework of the Coordinated Regional Climate Downscaling Experiment-Africa (CORDEX-Africa; Giorgi et al. 2009;Mariotti et al. 2014). 2. Studying the potential influence of different climate forcing using GCMs participating in the CMIP5/CMIP6 simulations (e.g. Erfanian et al. 2016;Mehboob et al. 2020) on the simulated PET during the reference period and under different global warming scenarios.
The above-mentioned future research works will be supportive to solve the meteorological and environmental related issues in the region.