Assessment of WRF Land Surface Model Performance over West Africa

Simulations with four land surface models (LSMs) (i.e., Noah, Noah-MP, Noah-MP with ground water GW option, and CLM4) using the Weather Research and Forecasting (WRF) model at 12 km horizontal grid resolution were carried out as two sets for 3 months (December–February 2011/2012 and July–September 2012) overWest Africa..e objective is to assess the performance of WRF LSMs in simulating meteorological parameters overWest Africa..emodel precipitation was assessed against TRMMwhile surface temperature was compared with the ERA-Interim reanalysis dataset. Results show that the LSMs performed differently for different variables in different land-surface conditions. Based on precipitation and temperature, Noah-MP GW is overall the best for all the variables and seasons in combination, while Noah came last. Specifically, Noah-MP GW performed best for JAS temperature and precipitation; CLM4 was the best in simulating DJF precipitation, while Noah was the best in simulating DJF temperature. Noah-MP GW has the wettest Sahel while Noah has the driest one. .e strength of the Tropical Easterly Jet (TEJ) is strongest in Noah-MP GW and Noah-MP compared with that in CLM4 and Noah. .e core of the African Easterly Jet (AEJ) lies around 12°N in Noah and 15°N for Noah-MP GW. Noah-MP GW and Noah-MP simulations have stronger influx of moisture advection from the southwesterly monsoonal wind than the CLM4 and Noah with Noah showing the least influx. Also, analysis of the evaporative fraction shows sharp gradient for Noah-MP GW and Noah-MP with wetter Sahel further to the north and further to the south for Noah. Noah-MP-GW has the highest amount of soil moisture, while the CLM4 has the least for both the JAS and DJF seasons. .e CLM4 has the highest LH for both DJF and JAS seasons but however has the least SH for both DJF and JAS seasons. .e principal difference between the LSMs is in the vegetation representation, description, and parameterization of the soil water column; hence, improvement is recommended in this regard.


Introduction
e land surface is the principal constituent located in the space separating the atmosphere and the lithosphere, and it obviously impacts the exchanges of energy and moisture with the boundary layer.
rough the controlling of the surface energy balance and water balance, land-surface processes strongly affect the weather and climate from local to regional and global scales [1]. One of the ways in which the land surface affects the climate system is the fact that due to the direct contact of the land surface with the atmosphere, the land surface reacts as a source and sink of heat and moisture through the sensible heat flux and evaporation. e surface conditions regulate the important feedback cycles in the climate system [2]. Also, the partitioning of the surface net radiation into sensible and latent heat fluxes determines the soil wetness development. e surface energy fluxes to a large extent control the surface weather parameters such as the temperature, humidity, and wind speed and to a lesser extent, low-level cloudiness and precipitation [2]. According to Tiwari et al. [3], it is noteworthy to have a better representation of the surface boundary conditions in a model as much as that of the surface and the atmosphere have an effect on the prognostic variables. e land surface and the atmosphere are inseparable and cannot be isolated as they are firmly coupled systems. e planetary boundary layer (layer where intense turbulence takes place) is the interface that regulates the feedbacks. Despite the significant influence of land-atmosphere exchange on the daytime planetary boundary layer (PBL), uncertainty remains in the parameterization of surface heat and moisture fluxes in numerical weather models [4]. e land surface model (LSM) provides heat and moisture fluxes over land to provide a lower boundary condition for vertical transport in the PBL scheme. ese fluxes of sensible heat and latent energy are dependent on the surface meteorology, radiative forcing, soil properties, and land use type. us, an accurate description of the land surface and vegetation characteristics is needed in any numerical weather prediction model [5]. e surface energy exchange is determined by the terrestrial radiation budget, as shown in the following equation: where R net is the net radiation flux (W/m 2 ), SH is the sensible heat flux (W/m 2 ), LH is the latent heat flux (W/m 2 ), and G is the ground heat flux (W/m 2 ). e surface available fluxes (i.e., SH and LH) are net radiation minus ground heat flux. e net radiation is balanced by outgoing fluxes of LH, SH, and G, whose partitioning strongly depends on the prevailing surface conditions. Some energy is absorbed by the ground G, this is in fact much lower than average value of LH and SH for most plant canopies, as most part of available energy will be transferred back into the atmosphere as sensible and latent heat [6]. LH is the quantity of heat absorbed or released by water undergoing a change of state. LH is most often the heat released by water as it changes from a liquid to gaseous state over a plant canopy through evapotranspiration. e partitioning of the energy available at the surface into latent and sensible heat depends crucially on the soil moisture [7]. Vegetated surfaces have the ability to draw water from a depth of order 1 m (the root layer), while for bare ground, only the water in the top few centimeters of soil contributes to evaporation. According to Zheng and Eltahir [8], vegetation cover and soil moisture content play a comparable role in the concept of land-atmosphere interactions. e main difference is that soil moisture anomaly patterns could last for many days to weeks, while vegetation is capable of mobilizing the root zone soil moisture that would otherwise not be in contact with the atmosphere. is consequently foists a lower boundary condition that could be in force for a longer time scale on the atmosphere. So, vegetation responds much slower on a single precipitation occurrence than the surface soil moisture. e land-atmosphere interactions are modulated by the extent of the associated north-south gradient of heat and moisture in the lower atmosphere [9]. A number of regions in the world, for example, the Sahel [10,11], Amazon [12], and Asian monsoon regions [13], have been identified as hot spots of land-atmosphere interactions, where interactions through feedback loops play a critical role in the surface water and energy balances as well as regional climate. In West Africa (WA) where we have strong land-atmosphere coupling, land surface processes, like soil-precipitation feedback [14], soil moisture initial conditions [14,15], and vegetation feedbacks [16,17], have significant impacts on the dynamic downscaling of regional climate models (RCMs). However, explanation of the results, from any one of such studies, ought to be tempered by the fact that there are considerable discrepancies in African land-atmosphere coupling strength among current state-of-the-art GCMs [18]. Li et al. [19] have demonstrated that land surface processes play a crucial role in the climate system. e land surface has been shown to be an important factor in modulating the West African monsoon (WAM) [20]. Nicholson [21] stressed that the land surface characteristics and processes have been shown to have a significant impact on the interannual variability of rainfall in the Sahel region based on observations. Furthermore, the importance of surface-atmosphere interactions was one of the main principles of the international African Monsoon Multidisciplinary Analysis (AMMA) project [22] which was later investigated in several studies. According to Dirmeyer [23], the WA region typically appears as the one where the soil moisture feedbacks with the atmosphere are among the strongest over the globe. e Sahel region has been identified as one of strong soil moisture-atmosphere coupling [10]. Furthermore, it has been resolved to be the region of the world with the highest impact of biophysical processes on the climate [11]. As reviewed by Xue et al. [24], studies by Steiner et al. [25] and Lavender et al. [26] have shown the importance of the land surface on modulating the West African Monsoon (WAM). For example, previous studies have examined the role of changes in the surface albedo, e.g., [27][28][29], and the vegetation, e.g., [30][31][32][33], on modulating the WAM. All of these studies lead to the general conclusion that reduced vegetation leads to reduced rainfall.
Parameterization is a method of replacing processes that are of a very small scale or very complex to be physically represented in the model through a simplified process. Researchers have focused on the appropriate selection of Weather Research and Forecasting (WRF) parameterizations schemes for varying conditions and applications [34][35][36][37][38][39][40][41][42][43][44][45][46][47]. Some of these studies concluded that the performance of the physics schemes varies according to the region of study including the time of the year. Hence, care should be taken in applying each for a specific study area. However, different climate variables are found to be sensitive to different physical parameterizations [37] making the need for an all-inclusive sensitivity analyses more pressing and the procedure of physics parameterization selection more demanding [48]. WRF has a number of physics options, which includes options for selecting the radiation scheme, surface layer scheme, planetary boundary layer (PBL) scheme, microphysics scheme, convective scheme, and land-surface scheme. is study assesses the performance of recent WRF 2

Advances in Meteorology
LSMs over West Africa for better simulations over the region. Land surface models (LSMs) differ in the number of soil and canopy layers and in the treatment of vegetation-related processes and hence are able to perform differently. Changes in surface conditions are shown to affect the position and intensity of the African Easterly Jet (AEJ) [49]. e AMMA Land Surface Model Intercomparison (ALMIP) was designed as a step towards a better understanding and explanation of surface processes over West Africa [50]. e idea was to develop a forcing database with the best quality and highest spatial and temporal resolution data available and use this database to force state-of-the-art LSMs in order to better understand key processes at different scales [51]. Concerning the LSMs, the AMMA results showed that they have a tendency to underestimate surface sensible heat flux and perhaps baseflow runoff and overestimate evaporation from the vegetation canopy. During the monsoon season, evapotranspiration is generally the most significant component of the water budget, acting to recycle much of the rainfall, in particular, over the Sahel [51]. Nevertheless, qualitative estimates of recycling by the LSMs must be interpreted in view of the fact that the LSMs have a substantial intermodel scatter over the Sahel [51].
Steiner et al. [25] coupled the community land model (CLM3) and the biosphere-atmosphere transfer scheme (BATS) with the ICTP regional climate model (RegCM3) and found out that CLM3 improves the timing of the monsoon advance and retreat across the Guinean Coast and reduces a positive precipitation bias in the Sahel and Northern Africa.
is results in a higher simulated temperature, and by that means, it reduces the negative temperature bias found in the Guinean Coast and Sahel in the BATS. ey also emphasized that the improvement in CLM3 is triggered by drier soil conditions in the simulation.
In spite of the notable impact of the land-atmosphere exchange on the daytime planetary boundary layer, there are still some uncertainties in the parameterization of surface heat and moisture fluxes in numerical weather models (NWM) [4]. WRF has many LSM options that range from simple treatment of subsurface processes like absence of vegetation or snow cover prediction to advanced physical models with sophisticated vegetation and soil models and snow cover prediction. Moreover, WRF LSMs use seasonally adjusted or annually fixed land-use properties read from look-up tables in order to allocate surface variables; yet, some land cover parameters like vegetation can change often within a short time.
Several sensitivity and evaluation studies have been carried out in different regions [5,19,[52][53][54][55][56]. is is necessary in order to help firstly the model developers for further improvements, secondly, the users in the selection of LSM options especially in the WRF model, and thirdly, improving the simulations for the area under study. erefore, considering the number of LSMs currently in the WRF model, it remains a difficult task for researchers to choose a suitable land surface scheme that fits their needs [57] especially over a data sparse region like West Africa. Although, some researchers [58][59][60] have replaced the land-use data by incorporating remote sensing data into the LSMs for more updated and realistic land-use land-cover data in order to improve the models; this was not done in this research because we wanted to evaluate the performance of each WRF LSM in the same standard mode.
Researchers such as Evans et al. [37] examined the performance of various physics scheme combinations on the simulation of a series of rainfall events near the southeast coast of Australia. ey created a thirty-six member multiphysics ensemble such that each member had a unique set of physics parameterizations. ey found out that no single ensemble member was found to perform best for all events, variables, and metrics, and this also reflects the fact that different climate variables are found to be sensitive to different physical parameterizations.
Hagos et al. [61] performed different WRF simulations with NOAH, SSIB, and PLEIM-XIU LSMs and cumulus parameterization schemes over West Africa. ey found out that the model simulations consistently differ from the observations over the Gulf of Guinea, with the model seemingly overestimating precipitation, and the simulations differ in their northward incursion of precipitation over land. eir analysis also showed that the African Easterly Jet (AEJ) position for the simulations was too far south for SSIB and too far north for PLEIM-XIU; hence, SSIB was too dry while PLEIM-XUI was too wet.
Wharton et al. [5] in their study found that the choice of land-surface model led to a ∼10% improvement in simulating hub-height wind speed and that the overall best performing models were Noah and Noah-MP. ey emphasized that the variability of LSM performance across different soil moisture and vegetation canopy conditions suggests LSM representation of surface energy exchange processes in WRF remains a large source of model uncertainty.
Jin et al. [53] studied the sensitivity of four land surface schemes in the WRF, namely, the simple soil thermal diffusion (STD) scheme, Noah scheme, RUC scheme, and community land model version 3 (CLM3). ey ran four simulations with four LSMs over the western United States.
eir results show that land-surface processes strongly affect temperature simulations over the area. ey also found out that WRF CLM3 with the highest complexity level significantly improves temperature simulations, except for the wintertime maximum temperature as compared with the observations. Also, precipitation was substantially overestimated by WRF with the LSMs over the area and does not show a close relationship with land-surface processes.
All these researchers and many more have performed different evaluations for different regions, but this research is focused on examining the role of land-surface processes in the simulation of some variables in West Africa by carrying out a series of WRF runs using four recent LSM options which includes Noah scheme, Noah-MP, Noah-MP GW, and Community Land Model Version 4 (CLM4), each having different complexity levels. ere are only few extensive studies on the assessment of regional climate models sensitivity to multiple land surface models over West Africa. erefore, the impact of different LSMs on West African

Advances in Meteorology
Monsoon (WAM) and land surface energy balance is still a study of importance. Atmospheric simulations over a variety of land-surface conditions that would be affected by the surface energy budget, including wet and dry seasons, were carried out in order to study the impact of LSMs choice on some WRF simulated parameters. In particular, we evaluate how the change in LSMs affect the average simulation of precipitation and 2 m temperature and tried explaining some reasons behind the differences using the mean wind, soil moisture, and surface energy parameters. e simulations were performed with a 12 km horizontal resolution for a domain that encompasses latitude 10°S-30°N and longitude 28°W-28°E as shown in Figure 1(a). Initial and lateral boundary conditions are from European Center for Medium-Range Weather Forecasts (ECMWF) Interim 6 hourly Re-Analysis (ERA-Interim) data, which has a horizontal resolution of 0.75° [62] and National Centers for Environmental Protection (NCEP) Final Analysis (FNL) 6 hourly initial soil data (sea surface temperature, soil moisture, and temperature) from the NCAR's Computational and Information System Laboratory Research Data Archive (CISL RDA) [63] with resolution of 1°. So, the initial and lateral boundary conditions were updated every 6 hours. e first 15 days were used as spin-up, hence 16 th of July to 29 th of September 2012 (JAS), for the rainy season and 16 th of December 2011 to 28 th of February 2012 (DJF), for the dry season were analysed. e variables were analysed spatially within the 20°W-20°E and 0°-20°N (the green box area) as shown in Figure 1, while the analysis involving time averages were within the 10°W-10°E and 5°-15°N boundary as shown in Figure 1(a) so that the domain of interest is moved away from the errors introduced by boundary values. Figure 1(a) also shows the elevation of the study area revealing the highlands such as the Jos Plateau, Cameroon Mountains, and Fouta Djalon highlands. e geogrid program by default will interpolate land use land cover LULC categories from the MODIS IGBP 21-category data in the model. e LULC distribution derived from MODIS data is shown in Figure 1(b). It shows barren or sparsely vegetated surfaces spread across the northern area (from 16°N upwards) and Evergreen Broadleaf Forest around the coastal regions. Most part of the Sahel is covered by grasslands with scattered croplands and open shrublands between the Sahel-Sahara interface. e distribution of the dominant soil type in West Africa is also shown in Figure 1 WRF version 3.9.1.1 has about seven different land surface physics options from the simple 5-layer soil model simply with thermal diffusion in soil layers and no vegetation or snow cover prediction to Noah LSM, Noah-MP, RUC LSM, PX LSM, CLM4, and SSiB LSMs with sophisticated vegetation models and snow cover prediction. e WRF LSMs are driven by surface energy and water fluxes and predict soil temperature and moisture in different layers depending on the LSM option in which Noah and Noah-MP has 4 layers, RUC has 9 layers, PX has 2 layers, and SSiB has 3 and 10 layers for CLM4. Noah-MP, RUC, SSiB, and CLM4 may predict snow water equivalent on ground. All the WRF land surface models utilize the wind and stability information from the atmospheric surface layer scheme.

Materials and Methods
e Noah LSM (as the most commonly used LSM) was used to test for four other different physics combinations (results not shown) in which the best combination was used with other three LSMs, which includes the Noah-MP, CLM4, and Noah-MP GW.
e Noah-MP GW is a new option including a free drainage soil lower boundary condition, variable water table, and 1-dimensional interaction with horizontal aquifer transport but no river routing or overland flow scheme as in hydrological models. All parameterization schemes apart from the LSMs were the same for all simulations. e new version of rapid radiative transfer model (RRTMG) based on Iacono et al. [64] is selected for describing the long-wave and short-wave radiative transfer within the atmosphere and to the surface.
e Mellor-Yamada Nakanishi and Niino Level 2.5 (MYNN2.5) scheme by Nakanishi and Niino [65] was used for boundarylayer processes with a consistent MYNN surface layer scheme.
e new Tiedtke scheme [66] is applied to parameterise the unresolved deep cumulus clouds. e 4 Advances in Meteorology microphysics scheme chosen is the WRF Single-Moment 6class (WSM6) scheme [67]. ese were selected based on some preliminary tests that also agreed with the findings of Gbode et al. [68] who carried out a verification study of many atmospheric physics combinations over West Africa. e Noah LSM [69] was developed jointly by National Center for Atmospheric Research (NCAR) and National Centers for Environmental Prediction (NCEP). It has been validated by many intercomparison studies both in coupled [70,71] and uncoupled [72][73][74] modes. Its moderate complexity and computational efficiency have made it effective in both operational weather and climate models. It has the advantage of being compatible with the time-dependent soil fields made available by NCEP in the global analysis datasets. It has a 4-layer (0.10, 0.30, 0.60, and 1 m thick) soil temperature and moisture model in addition to canopy moisture and snow cover prediction. It has a rooting depth fixed at 1 m including three layers at the top and also has one snow layer and one canopy layer. e mass conservation law and the diffusive form of Richard's law controls the vertical water mass movement between soil layers, while a conceptual parameterization for the subgrid treatment of soil moisture and precipitation controls the infiltration [75]. Total evapotranspiration in Noah is the sum of the canopy intercepted water evaporation, transpiration from vegetation canopies, and evaporation from bare soil as weighted by the respective land surface coverage fractions [76]. e drainage is due to the gravitational percolation  below the soil layers. e surface skin temperature is calculated from a surface energy balance equation. e surface energy fluxes of Noah are calculated over a combined surface layer of vegetation and bare soil surface. Such model hinders its further development as a process-based dynamic leaf model, because it cannot compute photosynthetically active radiation (PAR), canopy temperature, and related energy, water, and carbon fluxes explicitly [77]. Prediction includes a root zone, evapotranspiration, soil drainage, and runoff, taking into account vegetation categories, monthly vegetation fraction, and soil texture. It also predicts soil ice, and fractional snow cover effects, has urban treatment options, and takes surface emissivity and albedo properties into consideration.
Noah-MP (multiphysics) LSM [77] is an augmented version of the Noah LSM. e vital augmentation is the introduction of (a) a vegetation canopy layer for individually computing the canopy and the ground surface temperatures, achieved by introducing a semitile subgrid scheme to represent land surface heterogeneity [77], (b) a modified twostream scheme [78,79] for transfer of radiation through vegetation canopy and considering the canopy gaps, (c) a Ball-Berry scheme for the canopy stomatal resistance [80][81][82] which connects stomatal resistance to the photosynthesis of sunlit and shaded leaves, and (d) a short-term dynamic vegetation model [83] with two options (off and on) in which leaf area index (LAI) and vegetation greenness fraction (GVF) can be predicted from the model when turned on. Also, a simple TOPMODEL runoff model for the computation of surface runoff and groundwater discharge, a 3-layer snow model, is given in [84], and a frozen soil scheme with exceptional soil permeability is discussed by Niu and Yang [85]. e Noah-MP uses multiple options for key landatmosphere interaction processes. It contains a separate vegetation canopy defined by a canopy top and bottom with leaf physical and radiometric properties used in a twostream canopy radiation transfer scheme that includes shading effects [77]. Noah-MP contains a multilayer snow pack with liquid water storage and melt/refreeze capability and a snow-interception model describing loading/ unloading, melt/refreeze, and sublimation of the canopyintercepted snow. Various options are available for surface water infiltration and runoff, and groundwater transfer and storage including water table depth to an unconfined aquifer. Horizontal and vertical vegetation density can be prescribed or predicted using prognostic photosynthesis and dynamic vegetation models that allocate carbon to vegetation (leaf, stem, wood, and root) and soil carbon pools (fast and slow) [77].
CLM4 (Community Land Model Version 4 [86] is a state of the science LSM more often used in climate applications, developed at NCAR in cooperation with many external collaborators and includes sophisticated handling of biogeophysics, hydrology, biogeochemistry, and dynamic vegetation. e vertical structure comprises of a single-layer vegetation canopy, a five-layer snowpack, and a ten-layer soil column [87]. It has 3.8 m soil depth divided into 10 layers approximately at 1.8 cm, 4.5 cm, 9.1 cm, 16.6 cm, 28.9 cm, 49.3 cm, 82.9 cm, 138.3 cm, 229.6 cm, and 342.3 cm below the surface. e soil water is predicted from the modified Richards equation by Zeng and Decker [88]. CLM4 centers on biogeophysics of land surface including vegetation dynamics modules. e overland flow is calculated with a simple conceptual TOPMODEL approach to parameterise the surface runoff. An exchange of water between an unconfined aquifer and the overlying soil column is incorporated in the soil hydrology scheme [85]. e parameterization schemes of Niu and Yang [89] and Wang and Zeng [90] is used to calculate the snow cover and snow burial fraction. Some of the major differences between the LSMs are highlighted in Table 1.

Model Verification and Validation Datasets.
Simulated surface air temperature was assessed using the ERA-Interim dataset from the ECMWF and simulated precipitation using the Tropical Rainfall Measurement Mission (TRMM [92]). e 0.25 o resolution TRMM 3B42 product was used as the standard for evaluating the model outputs because of its reliability and also that it is a merged dataset (i.e., a combination of in situ and satellite products) with high-quality precipitation estimates [92]. e Global Precipitation Climatology Project (GPCP [93,94]) is another genuine source for merged estimates computed from microwave, infrared, and sounder data observed by the international constellation of precipitation-related satellites and precipitation gauge analysis. is was also analysed to see how good it is in comparison with TRMM over West Africa. ese datasets were used for the validation because the region (West Africa) is a data sparse region with a poorly spread synoptic weather station network, which is insufficient to validate the spatial and temporal spread of the model outputs. Hence, the use of several datasets for the model validation was necessary to examine their differences in order to clearly understand the model uncertainties.

Method of Model Evaluation.
e model evaluation was also done using the Taylor diagram that can concisely summarize the degree of agreement between simulated and observed variables [95]. e correspondence amongst different patterns is quantified by their correlation, their centred root-mean-squared difference, and their standard deviations. Houghton et al. [96] stressed that the diagram is useful specifically in evaluating multiple aspects of complex models or in gauging the relative skills of many dissimilar models. e location of each model on the diagram describes how closely the model's output pattern corresponds with the observations. However, each model point in the diagram represents the correlation, the standard deviation, and the centered RMS and is related by the following equations: where R is the correlation coefficient between the forecast and observation fields; E′ is the centered root-mean-squared (RMS) difference between the fields; and σ 2 f and σ 2 o are the variances of the forecast and the observation fields.

Precipitation and Surface Temperature Analysis.
From the plot of the average air temperature for JAS (Figure 2), all the simulations were in close resemblance to ERAinterim, which is used as the reference. FNL seems to have higher temperature than ERA-interim and all the models around 17°N and 20°N. e areas of high and low values in ERA-interim were captured by all the four LSMs. From Figure 3, the average temperature bias between ERA-Interim and the models for JAS shows that the outputs of all the model simulations agree with ERA-interim within the range of −4 to 3.5°C. e maximum positive bias was shown around the EL Djouf basin in Mauritania and Northern Mali, while the maximum negative bias was found around Lake Chad, possibly due to a bias in the lake temperature in the model, for all simulations. But, Noah shows a cold bias relative to ERA-interim in areas above 16°N at the northeastern part which other models did not reveal. Areas below 15°N mostly showed bias value close to zero except for CLM4. Nevertheless, it shows that all the models could fairly predict rainy season air temperature. Figure 4 shows that the spatial distribution of average 2 m surface temperature for DJF decreases with increasing latitude just like the reference from areas above 10°N. All the models and ERA-interim were able to capture reasonably well the cooling around the highlands of Jos Plateau, Guinea, and Cameroon. However, FNL has higher temperature around the coastal areas than ERA-Interim.
From Figure 5, it is seen that the average temperature bias between ERA-Interim and the models for DJF ranges between −4°and 5°C. ere is a maximum warm bias around Fouta Djallon highlands in all the simulations. Noah LSM output tends to have the coldest bias while CLM4 LSM tends to have the warmest. Despite having the most sophisticated snow, soil, and vegetation physics among the LSMs, CLM4 has the strongest positive temperature bias in some major highlands and the El Djouf basin. e deep soil column and water table in CLM4 might however need a considerably long time period for temperature and soil moisture to reach equilibrium, mostly in the cold and dry regions [97]. is could make it require more spin-up time than other LSMs and might be the reason for having high bias around the semiarid/arid areas of the region for both temperature and rainfall in JAS season.
From the time series of JAS daily 2 m surface temperature averaged over 10°W-10°E and 5°-15°N (as shown in Figure 6(a)), the general biases between ERA-Interim and the models can be clearly seen as ERA is predominantly cooler than the models. On average, August is the coolest period while September is the warmest, but the model tended to be too warm in September for all LSMs. From the time series of average daily temperature for DJF ( Figure 6(b)), the model patterns were very similar to the reference. Noah-MP and the ground water option were the closest while Noah was the furthest being 1-2°C colder than ERA-Interim. Figure 7 spatially displays the average JAS daily precipitation. ere seems to be a more intense precipitation core formed around the Fouta Djallons in Guinea both in TRMM and GPCP and all the models. e orographic effect of the highlands is the reason the models deliver a high precipitation amount but it appears to be more than that of the reference in those areas. When air moves over a mountain, it will normally be forced to rise. e forced rising air cools and expands as the temperature and pressure drop and consequently condenses and form clouds, if moist. is will enhance the rainfall activities as seen for the simulations and the reference datasets. Outputs from Noah and CLM4 models showed almost no rainfall from latitude 17°N northwards, whereas TRMM, GPCP, Noah-MP, and Noah-MP GW displayed some amount of precipitation. Also, the models differ in their northward advance of precipitation as compared with the reference. e 1-3 mm/day precipitation belt stopped around 16°N for Noah, 17°N in CLM4, and 19-20°N for Noah-MP and the ground water option, and by implication, it shows that Noah is relatively dry in comparison with other LSMs understudy. Hagos et al. [61] did a similar comparison with Noah, SSIB, and Pleim-Xui LSMs in WRF and also found out that the simulations differ in their northward incursion of precipitation, where the 2-4 mm/day precipitation band reaches Lake Chad for Noah LSM, while Pleim Xui was further north and SSIB was a little to the south. From the spatial bias of JAS precipitation as presented in Figure 8, the difference between the models and TRMM ranges between 5 and 10 mm/day. Positive bias >9 mm/day was noted around the Fouta Djallon highlands. Researchers like Jenkins et al. [98] and Jones et al. [99] have emphasized that orography plays an important role in West African precipitation patterns. Noah-MP GW showed the largest area of high precipitation bias among all the LSMs tested.

Advances in Meteorology
From Figure 9, all the models showed a good pattern of the West African precipitation as areas above 7°N were dry. All the models were able to capture the coastal precipitation during this period. However, Noah-MP and the ground water option (Figures 9(c) and 9(d)) seem to overestimate the DJF precipitation as the entire coastal region has precipitation, whereas the dry coasts of Ghana, Togo, Benin, and some parts of Nigeria as seen in TRMM and GPCP were demonstrated in Noah and a little in CLM4. From the spatial bias of precipitation for DJF in Figure 10, there was a clear overestimation in Noah-MP and the ground water option, while Noah simulation has the least spatial bias.
All the models showed maximum positive temperature bias around the EL Djouf basin in Mauritania and Northern Mali which could be as a result of the effect of the West African heat low prevalent in this area. is is a low surface         pressure region, which develops during the summer, linked with high insolation and seasonal surface temperatures. However, the lower bias prevalent for Noah-MP and Noah-MP GW than Noah in areas between latitude 11°and 14°N (which falls to the Savanna region) shows that the inclusion of a dynamic vegetation and groundwater process (as in the Noah-MP GW) can improve the simulation of temperature over this semiarid region of West Africa to some extent especially during the monsoon season. But, the area within 15°N and 20°N (the Sahel-Sahara interface) showed high positive bias more than any other region. Figure 11(a) shows that the model precipitation patterns agree with the reference in some instances and disagree in others. All the model simulations were completely in disagreement with TRMM on July 20 th and 21 st , August 2 nd and 24 th , and September 1 st . All the series have a similar pattern with the reference aside from the listed points. ere was a clear overestimation from Noah-MP GW option on 21 st of August. From Figure 11(b), the simulated precipitation pattern disagrees with the reference on 29 th of January and 22 nd of February but has a similar pattern in other periods. e very little or no precipitation notable for this period of the year in the region was well captured in all the models.
From the table showing the average values of the reference and all the model outputs (Table 2), it can be seen that most of the simulations have average values close to the observed value both for rainy and dry season except for precipitation with a little wide range in some instances. Also, Noah LSM gives a little different average value from other LSM, which could be related to the fact that it is the least sophisticated LSM among the ones under the study. However, the change in WRF LSM does not significantly change the average output of temperature over the area, and this could mean that the differences in some values might be as a result of atmospheric changes and also the data used for the initial and boundary conditions. e Noah MP and the GW option better simulates the minimum temperature than the other LSMs. According to Niu et al. [77] and Milovac et al. [100], the reason for this could be attributed to the high values of exchange coefficient in the Noah as compared with Noah MP. e high exchange coefficient leads to stronger ventilation from the surface and results in stronger surface cooling in the case of Noah [55]. e high exchange coefficient in Noah can overcome the blocking of ventilation by the moisture and create high temperature during the JAS season. During the DJF season, the temperature simulated by Noah MP was higher than that from Noah. For the JAS season, the Noah MP and Noah-MP GW simulated wetter condition that corresponds to a lower surface temperature than that from Noah. e mechanism behind this will be explained better in Section 3.2.
Using the same frequency and data as used in Figures 6  and 11 (i.e., daily averaged data series), Taylor diagrams [95] were used to evaluate the models performance in comparison with the ERA-interim (for surface temperature) and TRMM (for precipitation). e 2 m temperature and precipitation outputs from the eight different WRF runs at 12 km resolution for rainy (JAS) and dry (DJF) seasons with four different land surface schemes and options are compared in the diagrams as shown in Figures 12(a) and 12(b). For precipitation (Figure 12(a)), the correlation coefficient for dry season rainfall ranges from 0.69 to 0.74 for the simulations, while the GPCP has a correlation coefficient of 0.95. Standard deviations range from 0.62 to 1.48 for the simulations. For the rainy season precipitation over the whole region, the correlation coefficient ranges from 0.09 to 0.22 for the models and 0.94 for the GPCP. e possible reason the correlation coefficients in surface air temperature are generally higher than those of precipitation is because precipitation is more variable in time and harder to simulate precisely [101]. According to Chotamonsak et al. [102] and Chawla et al. [103], the deficiency of rainfall simulation with regional climate models is generally a common occurrence.
is results from the inability of the models to adequately tackle complex biogeochemical and biogeophysical processes [104] and the nonlinear interaction of the microphysics, cumulus, and planetary boundary layer schemes. Also, poor correlations with summer rainfall in major parts of the Guinea coast subregion in West Africa have also been 14 Advances in Meteorology suggested to be due the region's complexity associated with orographic and Atlantic Ocean Sea Surface Temperature (SST) influences [105]. Standard deviations range from 1.6 to 2.13 for the models and 0.8 for GPCP. erefore, there were higher standard deviations for precipitation during the rainy season, whereas there were higher values for surface air temperature during the dry season. Rainy season precipitation is best simulated with the Noah-MP GW, while CLM4 and Noah-MP almost performed equally for dry season precipitation simulations over the study region. For the 2 m temperature (Figure 12(b)), FNL shows a high correlation coefficient of 0.98, while the simulations  Noah LSM however performs best in simulating the dry season surface air temperature while the Noah-MP GW option performs best in simulating the rainy season surface air temperature. On average, from the four compared schemes based on the analysis of 2 m temperature and precipitation, Noah-MP GW best simulates the rainy season events while CLM4 is good for dry season events over the region. Although each LSM performed well for some variable and worse for others, but with a view to get the overall best performing LSM for all parameters as shown in Tables 3 and  4 below, a score table was formulated for each set of variables as named A, B, C, and D. Four points were assigned to each model that performed best for each variable, three were assigned for model that came second, two for third, and one for the worst performing model, i.e., fourth position. Noah got the highest score for DJF 2 m temperature, CLM4 got the highest for DJF precipitation, and Noah-MP GW got the highest score for JAS 2 m temperature and JAS precipitation. Based on the two most important meteorological parameters, Noah-MP GW is the overall best for all the variables and seasons in combination. However, if this is based on the JAS season (the monsoon season), then Noah-MP GW comes first, Noah-MP second, CLM4 came third, and Noah came fourth. is is also the most important season for weather and climate studies over the region as it is the period for crop growing as most farmers depend on the monsoon rainfall for agricultural activities due to poor mechanised farming and irrigation practices.

Atmospheric Circulation.
As a means to understand some of the mechanism behind the differences in the simulations with the LSMs, it is necessary to examine the atmospheric circulation and the main mechanism of moisture transport for the West African Monsoon (WAM). e WAM circulation is forced by the Sahara heat low (SHL) [17,49] and the connected midtropospheric high. e SHL drives the near surface and moistened southwesterly winds migrating from the Gulf of Guinea and the high associated with the heat low in turn drives the African Easterly Jet (AEJ), which transports moisture westward out of the Sahel region into the Atlantic Ocean [106,107]. According to Hagos et al. [61], the meridional temperature gradient together with the dry north-easterlies determines how far north the moisture is transported through the low-level south-westerlies, the latitudinal location of the AEJ, and finally the Sahel precipitation. So, the position and strength of the AEJ could influence the large-scale circulation over West Africa [49]. Also, at about 200 hPa, there exists the Tropical Easterly Jet (TEJ), which is associated with the South Asian monsoon outflow and propagates across West Africa during boreal summer [105]. Grist and Nicholson [108] have observed a weaker AEJ shifted to the north and stronger TEJ when rainfall is above average over the Sahel. Furthermore, the southwesterly monsoon wind flow at about 850 hpa is the main source of moisture that migrates from the Atlantic Ocean into West Africa. e strength and position of these systems are known to affect precipitation over West Africa [109] and also affects other connected parameters. . e RMS error and the standard deviation have been normalized according to [95] before plotting the Taylor diagram. 16 Advances in Meteorology From Figure 13, the core of the TEJ is located approximately at around 200 hPa in all the model runs. However, there are noticeable differences in terms of the strength of the TEJ, which is strongest in Noah-MP GW and Noah-MP (>16 m/s) compared with CLM4 and Noah (about 14 m/s). is suggests that the inclusion of the new option with interactive groundwater introduced into Noah-MP has the tendency to affect the simulation as West Africa and mostly the semiarid area of the region has a strong landatmosphere interaction (i.e., any change in land properties has a strong influence on the atmospheric conditions). Also, this interaction could increase the soil moisture and increase the latent heat fluxes/evapotranspiration that eventually affects the water vapour and cloud development, leading to precipitation increase as stronger TEJ enhances deep convection, which implies intensified precipitation.
For the AEJ, the strength of the jet (∼8 m/s) is almost the same across the four simulations and the core is located around 600 hPa.
ere are, however, slight differences in terms of the latitudinal position of the core as AEJ lies around 12°N in Noah, the southernmost, and Noah-MP GW resides around 15°N, the northernmost. is explains the reason why the model simulates low-level monsoon flow in Noah and a more enhanced flow in Noah-MP GW. is could be a reason for the dryer Sahel in Noah LSM, as the AEJ is known to transport more moisture out of the Sahel region and reduce precipitation. is is in line with the findings of Abiodun et al. [110] and Hagos et al. [61]. Figure 14 shows the low-level (850 hpa) wind vector and geopotential height from the simulations using the four LSMs for JAS period. It shows that Noah-MP GW and Noah-MP simulations have stronger influx of moisture advection from the southwesterly monsoonal wind than the CLM4 and Noah with Noah showing the least influx. is is evident as the flow from Noah-MP GW and Noah-MP goes further into the north than that from CLM4 and Noah. erefore, the weaker influx could also be a good reason why the 1-3 mm/day precipitation belt stopped around 16°N for Noah and 17°N in CLM4, while that of Noah-MP and Noah-MP GW extended up to 19-20°N as shown in Figure 7. Hence, the characteristics of the two easterly jets (AEJ and TEJ) and monsoon flow (at ∼850 hPa) in Noah-MP GW therefore provide the necessary environment for convective activities that favours the production of more rainfall events than other models.

Surface Energy and Soil Moisture.
In order to explain more on the reason or mechanism behind the differences as presented, soil moisture, sensible heat flux, and latent heat flux were also analysed. e LSM affects the energy balance at the land surface, which afterward affects the atmospheric conditions through land-atmosphere interactions [19]. In WRF, surface sensible heat (energy for direct heating of the air from the warmer ground) and latent heat (energy for evapotranspiration) fluxes are computed in the LSM and passed to the planetary boundary layer (PBL) and thus influences the simulation of West African Monsoon (WAM). Also, the LSM modulates the energy partitioning between the SH and LH fluxes which has a crucial impact on the atmospheric circulation as demonstrated in a number of studies [6,17,111,112]. Latent energy has been stressed to be the main contributor of differences between LSMs [25]. According to [113], higher latent heat fluxes increase the moist entropy flux and the amount of convective available potential energy in the boundary layer, likely increasing the frequency and magnitude of convective rainfall events. e spatial pattern of JAS latent heat fluxes from the four simulations ( Figure 15) matches that of precipitation (Figure 7) and soil moisture ( Figure 16). Also, the spatial pattern of JAS sensible heat flux from the four simulations (Figure 17) matches that of 2 m temperature ( Figure 2). Furthermore, the spatial pattern of DJF latent heat fluxes from the four simulations in Figure 18 is in line with that of precipitation in Figure 9. From Figure 19, it is seen that the DJF sensible heat is stronger in Noah-MP and Noah-MP GW than other simulations and also higher as seen in Table 5. is could be one of the reasons we have more DJF precipitation in Noah-MP and Noah-MP GW as increase in sensible heat could lead to increased boundary layer heating, increase convective clouds, and could eventually increase convective precipitation.
From Figures 20(a) and 20(b), it is clear that there were differences in the simulated sensible heat and latent heat flux among all the four simulations. e sensible heat flux on  average increases while the latent heat flux decreases from south to north. is could be due to the effect of the negative gradient in soil moisture from south to north as shown in Figures 16 and 20(d). is is also in line with the findings of Wolters et al. [114]. Also, the meridional variation of sensible heat flux (Figure 20(a)) and temperature (Figure 20(e)) are not very identical. is could be because the temperature is controlled not only by the sensible heat flux but also by advective effects [114]. Figure 20(c) shows the meridional structure of the zonally averaged evaporative fraction (EF). e EF is the ratio of latent heat flux to the sum of the sensible and latent heat fluxes. e southwesterly monsoon flow and local moisture recycling are the principal sources of moisture during the WAM [115]. e latter process contributes about 30% to the precipitation over the West African subcontinent [116]. From a simple one-dimensional perspective, increased value of evaporative fraction can be linked to higher values of PBL equivalent potential temperature and more favorable conditions for deep convection (e.g., [117,118]). Also, according to [61], simulations that show relatively dry conditions of the Sahel tend to shift the steep gradient further to the south and those that are relatively wet shift the gradient further north. is is in line with the sharp gradient from Noah-MP GW and Noah-MP ( Figure 20(c)) with wetter Sahel further to the north and for Noah, which is further to the south. e situation with the CLM4 is not in line with this as it is further to the north in the Sahel despite having dryer JAS Sahel than Noah-MP and Noah-MP GW.
Differences in simulated precipitation are probably to cause substantial differences in soil moisture. Also, the vegetation composition is different for all the LSMs, which affects the surface zone soil moisture. From Table 5, averaging the SMOIS for the entire soil column for all the simulations, Noah-MP GW has the highest amount of SMOIS while the CLM4 has the least for both the JAS and DJF seasons. e high soil moisture in Noah-MP GW could be attributed to the effect of the interaction with the ground water. e CLM4 has the highest LH for both DJF and JAS seasons but however has the least SH for both DJF and JAS

Advances in Meteorology 21
seasons. It could be deciphered that simulations with higher LH have lower SH and vice versa. e Noah scheme has a combined surface layer of vegetation and surface soil, while the Noah MP has a separated ground surface and canopy layer. Also, the bulk layer of snow/surface vegetation and soil associated with the Noah tends to under estimate the heat flux [77]. An important point for a lower sensible heat flux in Noah could be attributed to the underestimation of surface air temperatures during the dry season. e soil moisture considerably influences the latent heat flux, as high part of the surface energy has a tendency to be shared between sensible heat flux and soil heat flux in a dry condition. Just as the soil in most part of West Africa is practically dry in the DJF season in comparison with the JAS season, it leads to a reduction in the amount of latent heat flux in the DJF season. e reduced latent heat flux associated with the dry condition simulated by Noah MP and Noah-MP GW would result in corresponding increment in soil heat flux as against CLM4 with the highest LH in both DJF and JAS but still has the lowest soil moisture. e increased soil moisture content as seen in the Noah-MP GW increases the latent heat flux and decreases sensible heat flux and thereby decreases surface temperature as seen in the JAS period. e transmission of intermodel soil moisture differences to the atmosphere happens through the LH.
e Noah-MP GW simulations had wetter soil, produced more LH, and had upward transport of aquifer soil water (to upward flux of moisture) in JJA relative to the other Noah and Noah-MP simulations [119].
West Africa is one of the regions on earth having a strong coupling between soil moisture and precipitation [114], and this strong coupling is generally found mostly in transition   Advances in Meteorology zones between wet and dry climates like the Sahel [10]. e soil moisture controls the albedo and Bowen ratio. Eltahir [118] suggests that a positive feedback mechanism exists between soil moisture and rainfall. According to Wang et al. [120], higher soil moisture increases the latent heat flux, and this is liable to increase the moist entropy flux per unit mass of air and the amount of convective available potential energy in the boundary layer. ese processes probably increase the frequency and magnitude of convective precipitation events [113]. According to Steiner et al. [25], during the monsoon season, wetter soils and increased moisture advection could lead to an increase in the release of water vapor to the atmosphere, increasing the amount of latent energy from the surface. Wolters et al. [114] also concluded that soil moisture modifications could influence the 850 hPa flow pattern causing deviations in the path of the squall line over West Africa. ese results are in line with the conclusions of Hacker [121], who discussed the influence of soil moisture initialization in NWP on the predictability of winds in the PBL. is is also in line with our findings, as simulation with higher soil moisture tends to have higher precipitation as in the case with Noah, Noah-MP, and Noah-MP GW.
Changes in the distribution of roots can dramatically change the amount of soil moisture available to plants to transpire as the increase in root depth could increase water uptake, increase transpiration, increase water vapour, increase cloud formation, and increase rainfall [122]. While evaporation from bare soil and transpiration from short vegetation is mostly dependent on soil moisture conditions in the surface soil layer, tall vegetation can additionally access water in deeper layers [123]. However, despite having the lowest average amount of SMOIS (Table 5), the CLM4 has high average precipitation corroborated with high average LH. is could be explained as CLM4 has total soil depth of about 3.8 m and thus gives room for deeper penetration of plant roots to access deeper soil moisture, which could influence the LH. is also has the tendency to increase the transpiration and hence could explain why CLM4 has higher LH and EF than other LSMs and higher precipitation than Noah LSM. e difficulty in comparing CLM4 with the other tested LSMs here would be as a result of the differences in mostly the number of soil levels and might bring about some inconsistencies in averaging for better comparison. Also, the high average LH fluxes in CLM4 could lead to a greater evapotranspiration per unit precipitation which tends to return more water to the atmosphere, thereby leading to the drier soils/reduced soil moisture as seen for both DJF and JAS season. However, other factors beyond the scope of this study could be in operation.
Land surface models have become more important in numerical weather prediction models like WRF as more complex and sensitive planetary boundary layer schemes are used. e principal difference between the models is in the vegetation representation, the description, and parameterization of the soil water column in each model. Vegetation substantially influences the exchange of momentum, energy, moisture, and trace gases between the surface and the atmosphere and hence affects all other parameter in one way or the other. erefore, it is recommended that improvement on the treatment of surface parameters in the models should be paramount, as this will go long in improving the simulations especially in areas of strong land-atmosphere interaction.  Figure 1 for the JAS period. 24 Advances in Meteorology

Conclusion
In all, eight simulations with 4 LSMs (Noah, Noah-MP, CLM4, and Noah-MP GW) using the WRF model were carried out for two sets of 3 months' (December-February 2011/2012 and July-September 2012) period over West Africa. ese were done using the 3.9.1.1 version of WRF model to assess the performance of the listed LSMs over West Africa. e models' air temperature was evaluated using ERA-Interim, while precipitation was evaluated using TRMM. In order to provide a succinct statistical summary, the results were plotted in the Taylor diagram.
On average within the red box in Figure 1, CLM4, Noah-MP, and Noah-MP GW all overestimated the precipitation during both DJF and JAS seasons, while Noah underestimated it. e strength of the Tropical Easterly Jet (TEJ) is strongest in Noah-MP GW and Noah-MP compared with CLM4 and Noah. e core of the African Easterly Jet AEJ lies at the southernmost position in Noah and northernmost for Noah-MP GW. e Noah-MP GW and Noah-MP simulations have stronger influx of wind than the CLM4 and Noah with Noah showing the least influx. Also, the number of soil layers and the depth was a considerable factor as regards the cause of the differences between the models as CLM4 with higher soil depth behaved differently by having low soil moisture with high latent heat. e variable of interest is also a factor to be considered in the selection of LSMs. However, each LSM performed differently for different variables as a result of the difference in the model physics, soil layers, and columns, vegetation structure, and level of model sophistication. erefore, based on precipitation and 2 m temperatures, Noah-MP GW is the overall best for all the variables and seasons in combination. If this is based on the JAS season (the monsoon season), then Noah-MP GW will also come first, Noah-MP second, CLM4 third, and Noah will come fourth. e West African Monsoon (WAM) controls the rainy season in West Africa. Researchers like [124] have stressed that it supplies most of the precipitation because it transports rich moisture into the region through the southwest trade wind coming from the Atlantics. However, the WAM is a very complex large-scale circulation system, which entails the interactions of the elements of many scale atmospheric circulations, thereby making the close to accurate simulation of precipitation and other variables in the region a bit more demanding. So long as the improvement of the land surface models can result into a better simulation of the weather parameters, a more in-depth analysis of the impact of the recent WRF land surface models on major features of the general circulation of West Africa which includes the intertropical convergent zone (ITCZ), the monsoon depth and strength, the AEJ, TEJ, and AEW, and the deep convective systems is recommended as a good area for future research.
Precipitation could be better simulated more extensively by testing with several other sets of physics combinations like the radiation schemes, cumulus schemes, microphysics, and planetary boundary layer schemes in the WRF model. Also, a longer period of simulation to accommodate more seasons for better comparison between same seasons across each zone in West Africa could lead to a better simulation. From the study, it can be concluded that differences in land surface have the tendency to alter the weather parameters. Earlier studies have shown that the West African land surface is an important part in the region's climate dynamics. is has been confirmed in this study, as the studied parameters are very sensitive to each LSM used. However, differences between the land use in the model and the actual land use could cause discrepancies in the simulations. Also, the accuracy of initial and boundary conditions is crucial to the successful production of WRF simulations [45] as the difference could propagate through the atmosphere and have an effect on the simulation of the dynamical features in West Africa. Hence, future research would also incorporate a real-time land use data products into the model over West Africa with various initial and boundary conditions and probably compare with an in situ data (if made available) in order to improve the simulations and have a perfect comparison.
Finally, despite having constant roughness length values for all season, Noah-MP still performed better than Noah over West Africa.
is could be due to the other improvements and modifications in the model, but the constant roughness length in the model is a factor to look into.
is could be a good area for future improvement as to reduce uncertainties in the model for better simulations and forecasts.
Data Availability e data used for the research are archived in the NCAR HPSS storage. is could only be made available from the corresponding author upon reasonable request.

Conflicts of Interest
e authors declare no conflicts of interest.