Next Article in Journal
Improving Ensemble Forecasting Using Total Least Squares and Lead-Time Dependent Bias Correction
Next Article in Special Issue
Wind Speed Forecast Based on Post-Processing of Numerical Weather Predictions Using a Gradient Boosting Decision Tree Algorithm
Previous Article in Journal
Tourists’ Perceptions on Climate Change in Lisbon Region
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Warm Rain in Southern West Africa: A Case Study at Savè

by
Irene Reinares Martínez
1,*,†,
Jean-Pierre Chaboureau
1 and
Jan Handwerker
2
1
Laboratoire d’Aérologie, Université de Toulouse, CNRS, UPS, 31000 Toulouse, France
2
Institute of Meteorology and Climate Research, Karlsruhe Institute of Technology (KIT), 76344 Karlsruhe, Germany
*
Author to whom correspondence should be addressed.
Current address: LACy, Laboratoire de l’Atmosphère et des Cyclones, CNRS, Université de la Réunion, Météo-France, 97490 Saint-Denis de La Réunion, France.
Atmosphere 2020, 11(3), 298; https://doi.org/10.3390/atmos11030298
Submission received: 24 February 2020 / Revised: 11 March 2020 / Accepted: 14 March 2020 / Published: 19 March 2020

Abstract

:
A warm-rain episode over southern West Africa is analyzed using unprecedented X-band radar observations from Savè, Benin and a Large-Eddy Simulation (LES) over a 240 × 240 km 2 domain. While warm rain contributes to 1% of the total rainfall in the LES, its spatial extent accounts for 24% of the area covered by rainfall. Almost all the warm-rain cells tracked in the observation and the LES have a size between 2 and 10 km and a lifetime varying from 5 to 60 min. During the nighttime, warm-rain cells are caused by the dissipation of large deep-convection systems while during the daytime they are formed by the boundary-layer thermals. The vertical extension of the warm-rain cells is limited by vertical wind shear at their top. In the simulation, their top is 1.6 km higher with respect to the radar observations due to the large-scale environment given by wrong initial conditions. This study shows the challenge of simulating warm rain in southern West Africa, a key phenomenon during the little dry season.

1. Introduction

Rainfall types over West Africa range from shallow, warm-rain-producing cloud elements to mesoscale convective systems (MCSs). While it is now commonly known that MCSs account for the vast majority of rainfall in the Sahel (e.g., [1,2]), there is a lack of robust rainfall type statistics for the wetter and more populous region of southern West Africa (SWA) lying to the south of the Sahel. Among the rainfall types, warm rain formed by clouds composed only of liquid hydrometeors has been understudied in the past.
One reason is the small size of these shallow clouds, on the order of 1 km [3], which makes their observation from space a challenge. Generally, satellites do not allow an accurate estimation of warm-rain amounts. Indeed, rainfall estimations highly rely on passive microwave radiometers and infrared imagers, which are mainly sensitive to ice particles and high-cloud tops. Thus, warm rain associated with shallow clouds may be misclassified when the latter are attached to MCSs or in the case of multi-layer cloud systems (e.g., [3,4]). Furthermore, instantaneous snapshots are not sufficient to determine if the shallow clouds are in their terminal stage or if they are a transient stage of evolving deep convection [5]. This can lead to uncertainties in the warm-rain amount calculations. Using a 9-year climatology of Tropical Rainfall Measuring Mission (TRMM) precipitation radar measurements [3] estimated the contribution of warm rain to total rainfall. They found a contribution of 20% over tropical oceans where shallow clouds occur predominantly (e.g., [6,7]) and of 7.5% over land, where the warm-rain clouds are deeper [3,8]. Note that the probability of rain from warm clouds over land increases towards coastal areas as shown by [7] using 5-year Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) data.
Over SWA, some recent efforts have been undertaken to shed more light on the distribution and regional importance of warm rain. Based on a 16-year TRMM precipitation radar observation, [9] estimated that warm rain contributes at least 2% to total annual rainfall at the immediate coastal region of SWA and that such warm rain events predominantly occur in July and August. [10] developed a SWA-calibrated method to delineate precipitating liquid clouds using the spatio-temporally high-resolution dataset of the Spinning Enhanced Visible and Infrared Imager (SEVIRI). They found a unimodal seasonal cycle of the occurrence frequency with a distinct peak in August, which corresponds well with the results in [9]. They further estimated warm rain frequency over SWA, and found it to be larger over the coastlines, from 5% in Nigeria to 15% in Liberia. Despite the likely small contribution to annual amounts, the larger precipitation frequency yielded by warm rain may be important for agriculture, in contrast to the intense MCS-like rainfall that can generate floods and runoff over dry soils. The latter is particularly important during the “little dry season”, observed from mid-July to mid-September along the SWA coast, sometimes associated with droughts [11]. This motivates a better understanding of warm rain.
Current operational numerical-weather prediction models and convection-permitting research models run with horizontal grid meshes in the range of 1–20 km. The development of the boundary layer and the transition from stratus to convective clouds is at best partly resolved by these models, but usually fully parameterized. Their ability to represent warm rain is thus limited. Large-Eddy Simulations (LESs) can explicitly resolve the most energetic eddies in the boundary layer. Therefore, they are expected to better represent the warm-rain clouds. Among the modeling studies devoted to shallow convective clouds (e.g., [12,13,14,15]), only the work by [14] gives estimates of warm rain contributions from shallow (3 km top) and congestus (7 km top) clouds; the contributions are 3.6% and 34%, respectively. Their analysis is based on a 10-day case transitioning from suppressed to active phases of the Madden–Julian Oscillation.
In this paper, we investigate warm rain in SWA during the field campaign of the Dynamics–Aerosol–Chemistry–Cloud Interactions in West Africa (DACCIWA) project [16]. This field campaign took place in summer 2016. Intensive ground and radiosonde measurements were performed from mid-June to the end of July 2016 at three supersites in SWA: Kumasi (Ghana), Savè (Benin), and Ile-Ife (Nigeria). Savè was instrumented with an X-band rain radar deployed by the Karlsruhe Institute of Technology (KIT) [17]. Also, a radiosonde campaign run parallel in Savè and Kumasi (see [18] for details). Here, we use observations gathered during the campaign and a LES to explore one warm-rain episode. The objectives of the paper are twofold. First, to present a set of unprecedented observations made at the Savè supersite. There, the most remarkable event of warm rain was found on 24 July 2016, it was observed during the night and early morning. This case yielded 1.6 mm, as measured by the distrometer and the rain gauge, roughly half of the warm rain in Savè during the whole campaign. Other notable warm-rain cases were the 21 and 25 June, and the 27 July 2016, responsible for 0.31, 0.42 and 0.36 mm, respectively. In total, 3 mm were yielded by warm rain during the campaign, which represented 1.35% of total rainfall in Savè. An added value of this study is the use of the X-band rain radar, as it collects information on the vertical structure of the warm-rain cells and allows a characterization of warm rain at high spatial and temporal resolution. Both the rain radar and the LES are used to characterize warm rain. The second objective is to determine the mechanisms involved in the development of precipitating ice-free clouds in SWA. The case is analyzed based on a 200-m resolution LES performed with the research model Meso-NH [19]. This grid spacing allows a reasonable representation of turbulent fluxes in the boundary layer and meets our need to simulate a large domain (240 × 240 km 2 ) covering the X-band radar range. The study by [14] has proven to resolve shallow convection with this resolution.
Our methodology to explore this event relies on a comparison of the simulated warm-rain cells and those observed by the X-band radar. The detection of warm-rain cells that lie, by definition, below 4-km altitude, is based on the clustering of reflectivity fields to delimit a volume of precipitating hydrometeors. This allows the characterization of the warm-rain cells in terms of their size, spatial distribution and vertical development. Then, the warm-rain cells are tracked in time, and the lifetime of each track is computed. To the best of the authors’ knowledge, this is the first time that the size or lifetime of warm-rain cells over SWA are statistically quantified.
The paper is organized as follows: Section 2 describes the observational datasets employed to characterize warm rain, particularly the X-band radar and the Meso-NH simulation, and the analysis approach. The case study is presented in Section 3, and the contribution of warm rain to total rainfall in the model is calculated. In Section 4 we analyse the properties of the warm-rain cells for both the radar and the simulation. The vertical profiles of the atmosphere and how they explain the differences between the observed and simulated warm-rain cells are discussed in Section 5. The water content, equivalent potential temperature, turbulence kinetic energy and wind shear inside the warm-rain-producing clouds are also explored. Section 6 concludes the paper.

2. Data and Methods

2.1. Observations from Savè Supersite

Warm-rain cells are defined in terms of reflectivities measured by the X-band radar at Savè (8.00 N, 2.43 E), in Benin. The radar is located at 183 m altitude. It is a Meteor50DX radar emitting at a frequency of 9.37 GHz, or equivalently a wavelength of ∼3 cm, and its horizontal range is set to 100 km. The reflectivity data is available in a polar grid with 12 elevation angles ranging from 0.7 to 90 and a resolution of 500 m and 1 in the radial and azimuth coordinates, respectively. To illustrate the radar scan strategy, Figure 1 shows the altitudes and distances scanned for each elevation angle at a constant azimuth. For example, the reflectivity at 4-km altitude is only measured at a certain distance for each elevation angle; the regions which are not scanned are shaded. The uncertainty in the estimation of the echo heights is inherent to the radar observation, and becomes larger as the distance from the radar increases. The whole volume comprising all azimuths is scanned every 5 min. A calibration factor (7 dBZ) is applied to the logarithmic reflectivities due to a systematic underestimation.
The Savè supersite also took advantage of radiosonde observations daily at 06:00 UTC providing profiles of the atmosphere up to the stratosphere; and of additional high-frequency boundary-layer soundings during the Intensive Observation Periods; Modem M10 meteorological radiosondes were used. For this study, we use the standard observation on 24 July 2016 at 06:00 UTC and two frequent observations at 03:30 and 08:00 UTC. These soundings serve to evaluate the simulated vertical structure of the lower atmosphere in Savè. The rain rate derived from a Joss-Waldvogel-Distrometer and the measurements of reflectivity at vertical incidence from a cloud radar are also used. Further details on the instrumentation of the Savè supersite and the measurements done are given in [17,18].

2.2. Meso-NH Large-Eddy Simulation

The non-hydrostatic mesoscale atmospheric Meso-NH model [19], version 5-4-0, is run with two nested grids with domains centred over the X-band radar location (Figure 2), using a one-way nesting method [20]. Meso-NH is run for 18 h on a large domain covering 1024 × 1024 km 2 with a 1-km horizontal grid. It starts from ECMWF (European Centre for Medium-Range Weather Forecast) operational analysis at 18:00 UTC 23 July 2016. This initialization time allows one to simulate the nocturnal low-level jet (see [21] for a conceptual model of the nocturnal boundary layer in Savè). Boundary conditions are provided by ECMWF operational analysis at 00:00, 06:00 and 12:00 UTC 24 July 2016. Because the period of interest covers 03:00 to 08:00 UTC (see Section 3), the model is run for 10 h with a 200-m grid spacing over a 240 × 240 km 2 domain, starting from the first run fields at 02:00 UTC 24 July 2016. Such a high resolution allows for the terrain, surface roughness and boundary-layer processes to be fairly represented. In this study, only the results from the 200-m simulation are analyzed and discussed. The vertical grid counts 64 levels up to 24-km altitude. A fine vertical spacing is set, varying from 40 to 200 m up to 5 km (approximately the height of the melting layer) in order to well resolve the processes in the edges of the warm-rain clouds. The vertical spacing stretches then up to ∼1 km above the tropopause. Output from the 200-m simulation is generated every 5 min, which is the X-band radar data frequency.
Both runs use the Surface Externalisée (SURFEX) scheme for surface fluxes [22], the Rapid Radiative Transfer Model [23] for longwave radiation and the two-stream scheme [24] for shortwave radiation. The 3-D turbulence scheme employed is a 1.5-order closure scheme [25] with the mixing length being the cubic root of the grid volume. A parameterization of dry thermals and shallow cumuli [26] and a subgrid statistical cloud scheme [27] are activated for the 1-km run while no shallow convection or subgrid cloud parameterizations are used for the 200-m run. We use a microphysical scheme for mixed-phase clouds [28]. In this scheme, the warm processes rely on the 1-moment Kessler scheme in which autoconversion of cloud water into rain is initiated beyond the critical cloud water content value of 0.5 g m - 3 . Finally, from the output of the 200-m simulation, reflectivities are computed in the Rayleigh approximation on each grid point of the model following [29].

2.3. Detection and Tracking of Warm-Rain Cells

Warm rain is detected every 5 min from 03:00 to 08:00 UTC. To ease the processing, the native polar grid of the radar reflectivities ( Δ r = 500 m) is transformed into a Cartesian grid with Δ x = Δ y = Δ z = 500 m. For both the X-band radar observations and the simulation, the identification is done on the domain covering the radar range with a shared horizontal grid mesh of 1 km and the native vertical grid of each dataset. The method relies on the clustering algorithm developed by [30]. It consists of locating the grid points where a condition is verified, and grouping the contiguous points into a three-dimensional object, herein named warm-rain cell. To apply this algorithm, the 3-D reflectivities are first projected into the horizontal, such that the resulting fields correspond to the maximum values of reflectivity over the column. Then, a condition on the intensity of the maximum reflectivity is used to delimit the horizontal spread of the rain cells: it must be larger than 20 dBZ. Thresholds ranging from 18 to 20 dBZ are commonly used to identify light to moderate rain above 0.5 mm h - 1 (e.g., [3,31,32,33]). Figure 3 shows the cumulative frequency distribution of the maximum reflectivity values to which the 20 dBZ threshold is applied, for the radar and model data. The distribution is shifted to lower values in the model with respect to the radar. While 16.9% of the observed reflectivity values exceed 20 dBZ, only 8.6% do in the simulation. A top is then associated with each identified cluster, defined as the highest altitude where a reflectivity of 20 dBZ is found inside the cluster for both datasets. A base is also computed as the lowest altitude where a reflectivity of 20 dBZ is found inside the cluster; but it is only done for the simulated rain cells. This is an added value of the simulation; it is not possible to compute the base height for the radar. Notice that, for the latter, the three first elevation angles are partially blocked by obstacles (e.g., bushes, buildings). The most severe blocking occurs for the first elevation angle (0.7 ), with roughly 1/3 of the total area blocked, mainly covering the northeastern quadrant of the radar range. Each cluster with an effective diameter larger than 2 km (smaller clusters are considered as noise) and with the top located below 4 km is classified as a warm-rain cell. The effective diameter is defined as the diameter of a circle whose area equals the area covered by the cluster. We consider that this 4-km altitude corresponds to the freezing point temperature of 273 K in the clouds. It is, for instance, the same height that was used within TRMM to identify warm rain in the tropics by [6]. All the clouds with higher, thus colder tops cannot produce warm rain. Notwithstanding, in the clouds, it is very common to find supercooled water. In the absence of sufficiently precise data on this phenomenon, we decided to keep 4-km altitude threshold for warm rain, being aware of the limitation of this choice. To illustrate how the detection method is applied, Figure 4 shows a map of the reflectivity obtained from the X-band radar and the warm-rain cells identified at 06:00 UTC. The high-frequency of the detection method allows one to monitor the evolution of the number and characteristics of the cells, such as their size and lifetime, as will be shown in Section 4. Note that the lifetime of warm-rain cells should not be confused with that of warm clouds.
The simulated warm rain at ground is estimated for each cell based on the 5-min total accumulated rainfall maps. The rainfall gridpoints matching the individual cells are attributed to warm rain. As shown in Section 3, warm rain is accumulated over a 6-h period, from 03:00 to 09:00 UTC and over the entire model domain; the clustering algorithm for warm-rain cells is thus extended accordingly.
Finally, the tracking of warm-rain cells over time (developed by [30] to track deep-convective clouds) is done in a way that two warm-rain cells detected at two consecutive time steps are the same if the surface overlap between them is above a certain threshold; different thresholds are tested (see Section 4) and applied for both the observation and the LES.

3. Clouds and Rainfall

The vertical structure of clouds over Savè on 24 July 2016 is shown in Figure 5. During the night and early morning values of reflectivity larger than −10 dBZ are measured by the cloud radar below 2–3 km, corresponding to the presence of clouds. These clouds having tops below the freezing level are in liquid phase. Slightly before 04:00 UTC, from 05:00 to 06:00 UTC and at 07:30 UTC warm rain from these shallow clouds is measured by the distrometer, with intensities ranging from 1 to 9 mm h - 1 (Figure 5, bottom). The warm rain amount is 1.6 mm, approximately half of the warm rain measured in Savè during the whole campaign from 19 June to 30 July 2016. This case and the three other warm rain cases observed during the field campaign yielded a total of 3 mm, which represents 1.35% of total rainfall in Savè. Deep-convective clouds develop in the afternoon and evening, reaching 8–10 km altitude before 21:00 UTC and up to 12 km after 21:00 UTC. In the following, our statistical analysis of the warm-rain cells based on the X-band radar is done from 03:00 to 08:00 UTC, when warm rain is reported at Savè.
We inspect how total rainfall is spatially distributed in the simulation, and compare that with two gridded satellite observations, the TRMM 3B42V7 product (3 h, 0.25 ) and the GPM IMERGV06B product (30 min, 0.1 ). Figure 6 shows maps of the total accumulated rainfall from 03:00 to 09:00 UTC for the three datasets, as well as warm rain for the simulation.
The analysis of rainfall at ground is done from 03:00 to 09:00 UTC (instead of 08:00 UTC) due to the 3-hourly temporal resolution of the TRMM product. Notice that although the colorbars are the same in Figure 6, their values differ between the maps for total and the map for warm rain. Both satellite products retrieve rainfall with values of 5 mm or more to the west of the domain. It corresponds to the signature of a MCS over western Benin and Togo. While TRMM only shows precipitation from the MCS, GPM retrieves rainfall with values up to 7 mm over the eastern half of Benin and over Nigeria. Inside the X-band radar range, the small patterns oriented in the southwest-northeast direction suggest they correspond to warm-rain features or isolated storms. Moreover, GPM retrieves 0.2–0.4 mm over the gridpoint adjacent to Savè, although it is a lower value than the 1.6 mm measured at the supersite. The discrepancy between TRMM and GPM points out the challenge that satellites currently face to capture warm rain and shows the potential of GPM with respect to TRMM. In the simulation, rainfall over Togo is missed, because the MCS is simulated further west (in the 1-km resolution simulation) and is smaller than observed (not shown). Over the eastern boundaries of the domain, large precipitating systems are simulated, in partial agreement with GPM. Inside the radar range, rainfall greater than 5 mm is unevenly distributed and mainly corresponds to deep convection. Smaller precipitation features of less than 3 mm may correspond to shallow clouds. Rainfall below 3 mm mimics the southwest-northeast orientation seen in GPM. The simulated warm rain is sparsely scattered. It sometimes appears as strips or structures with the GPM-like orientation.
The total rainfall observed by TRMM and GPM is 0.8 and 1.5 mm, respectively, and it is 1.0 mm for the simulation. This is encouraging given the short frame of the reporting period. We also estimated the contribution of warm rain to total rainfall and to the area covered by rain. These contributions cannot be calculated for satellite products due to the lack of observed vertical distribution of precipitating hydrometeors. For the simulation, the contribution of warm rain to total rainfall is 1.2%. Although this contribution is only calculated for a 6-h period (from 03:00 to 09:00 UTC), its order of magnitude compares well with the 1.35% warm rain recorded at Savè during the whole campaign and the climatological value of 2% obtained by [9]. While warm rain contributes little to rainfall, its contribution to the area covered by rain is much larger: almost one fourth (24.4%) of the area covered by rainfall is covered by warm rain.

4. Characteristics of the Warm-Rain Cells

The relative distributions of the effective diameter of the warm-rain cells and the lifetime of the tracked warm-rain cells are represented in Figure 7 for the X-band radar and the simulation. The lifetime distribution is shown for a 20% geographical overlap. In both datasets, the number of warm-rain cells decreases with increasing diameter. The total number of 1057 observed warm-rain cells is in agreement with the number detected by the X-band radar, which is in average 20 per 5 min scan. Only 485 warm-rain cells are simulated. Although the simulated warm-rain cell population counts about half the observed number, the relative distribution is well captured. The large majority of warm-rain cells have diameters below 5 km, approximately 78% in the observations and 81% in the simulation. 17% and 16% warm-rain cells in the observations and the simulation, respectively, have diameters comprised between 5 and 10 km. The remaining 5% and 3% correspond to cells with diameters larger than 10 km. The largest diameters are 17.3 km in the radar and 15.9 km in the simulation.
The number (225) of observed tracks (warm-rain cells detected at least at two successive steps) approximately doubles the number of the ones that are simulated (107). This is due to an underestimation of the number of simulated warm-rain cells during the nighttime, as will be shown later. Among the total number of tracks, 67% last up to 10 min for both the radar and the model and 31% and 32% of the observed and simulated tracks, respectively, last between 15 min and 1 h. Lifetime rarely exceeds 1 h. The distribution was calculated using a 20% overlap value. Values from 5 to 40% with an increase of 5% were tested and the sensitivity of the distribution to the threshold has been found to be generally low. Overlapping thresholds in the range from 5 to 30% have little to no impact on the duration of the tracks, with roughly 2/3 of the tracked warm-rain cells that exist for 5 to 10 min, 1/3 that last between 15 and 60 min and 1 to 3% that last longer. Adopting a threshold of 35% or more is, however, more restrictive, such that the frequency of 5–10 min duration reaches about 3/4 of the tracked warm-rain cells. This results in a diminution of the mean lifetime. This is the first time that the lifetime of warm-rain cells is quantified, and has been possible because of their individual identification.
The time evolution of the number and effective diameter of the tracked warm-rain cells is presented in Figure 8 for the radar and the simulation. The simulated depth (calculated as the difference between the top and the base) is also shown. The number of observed tracks increases from around 10 detections per 5 min scan at 03:00 UTC to 20–30 detections at 06:00 UTC. Later on, their number diminishes. Very few warm-rain cells are simulated during the nighttime, then the population grows fast from about 5 detections at 06:00 UTC, which corresponds to sunrise in Savè [17], to more than 20 towards the end of the analysis period. The triggering of warm-rain clouds is enhanced due to the daytime increase of the surface fluxes and the development of the boundary layer. Overall, the mean effective diameter slightly increases for both the observed and simulated tracks. The simulated warm-rain cells also become deeper with time as the surface fluxes increase, from 1-km depth at 05:00 UTC to 3-km depth by 08:00 UTC. In average, the observed tracked warm-rain cells have their tops at 2315 m high. In the LES, the average base and top are 1709 and 3876 m, respectively.
Figure 9 and Figure 10 show the spatial distribution of the tracked warm-rain cells; the colors display the hour of detection for the radar and the LES in Figure 9 and the simulated depth of the warm-rain cells in Figure 10. The solid dots depict the first detection of each track, the empty dots depict the consequent detections. The majority of the observed nighttime warm-rain cells (from 03:00 to 06:00 UTC) are located inside the outer 50-km radar range, in the western half of the domain, whereas the daytime ones (fewer in number) are mainly in the inner 50-km range, some of them lying slightly south. The simulated tracked warm-rain cells, less numerous, are generally distributed over the northeastern quadrant, as seen from Figure 8. The vast majority appear after 06:00 UTC. The observed nighttime warm-rain cells present higher tops. An analysis of MSG-derived brightness temperature indicates that during the nighttime, the observed warm-rain cells tend to develop near cold cloud tops above the 0 C altitude (not shown), and probably correspond to transient stages of deep convection. Some of them are located in the vicinity of the MCS. In the simulation, the vertical extension of the clouds covers the full range of depths.
For the radar, the trajectories of the warm-rain cells show two predominant directions, probably corresponding to two different wind regimes. The nighttime warm-rain cells in the northwestern part of the outer 50-km range, where they generally have higher tops, present a southeasterly propagation. Another group of nighttime and daytime warm-rain cells show preferentially a southwesterly propagation. Finally, a few tracks are rather stationary and appear near orographic features. In the simulation, they mainly show a northeasterly propagation. The long tracks seem to span over longer distances than for the radar. This is possible if they are found within a larger wind speed regime.
In summary, almost every tracked warm-rain cell (98%) lasts 60 min at most, according to the X-band radar. The tracked warm-rain cells present mean effective diameter and top height of 4.9 and 2.3 km, respectively. The LES captures remarkably well the size and lifetime relative distributions. However, the total number of simulated warm-rain cells is underestimated, mainly during the nighttime. They also present tops too high by ∼1.6 km in average. This can affect the amount of warm rain by changing the height in the atmosphere and the path within the cloud in which the collision-coalescence processes occur. The LES has trouble representing the spatial distribution of warm-rain cells partly because of the absence of the MCS, which might also explain the lack of nocturnal warm-rain cells. The LES compares better with the observations during the day, but the main drawback is the simulated southwestward propagation compared to the observed northeastward motion.

5. Vertical Structure of the Atmosphere and the Warm-Rain Cells

To better understand the main processes driving the development of warm-rain-producing clouds in the observation and the LES, the vertical structure of the atmosphere at Savè on 24 July 2016 is now examined as well as the vertical properties of the simulated tracked warm-rain cells.
The vertical profiles of potential temperature, water vapor mixing ratio and wind speed and direction from the radiosoundings at Savè and from the simulation are shown in Figure 11 at 03:30, 06:00 and 08:00 UTC. Profiles of the same variables are also displayed for the ECMWF operational analysis at 06:00 UTC (i.e., when available). The observed potential temperature shows the static stability of the atmosphere at 03:30 and 06:00 UTC. At 08:00 UTC, the daytime development of the boundary layer causes a neutral stratification in the first 200 m. The observed water vapor mixing ratio at the surface takes values between 16 and 17 g kg - 1 , then it decreases with altitude. These near-surface values of moisture content lie within typical values at the Savè site [17]. The profiles at 03:30 and 06:00 UTC are moister than at 08:00 UTC in the lowest 2 km, probably due to daytime mixing. At 03:30 UTC, the observed wind speed shows a maximum of more than 6 m s - 1 at around 700 m altitude, with a southwesterly direction. This is the signature of the nocturnal low-level jet, with an intensity and core height typical of Savè [34,35]. At 06:00 UTC, the low-level jet has disappeared and southwesterlies of 3 to 7 m s - 1 are observed up to around 2 km, corresponding to the monsoon flow. The mean depth and wind speed of the monsoon layer at Savè are 1.9 km and 6 m s - 1 [17]. At 08:00 UTC, the wind speed shows two maxima of 3–4 m s - 1 .
The LES presents a warm bias, of 1 to 2 K, between the surface and around 2 km with respect to the observations. This bias is the largest at 03:30 and then decreases. The profiles of water vapor mixing ratio at 03:30 and 06:00 UTC capture well the values from the soundings up to 1.5 km. Then, they depart from the observations and become drier. At 08:00 UTC, the simulation overall captures the observed values of water vapor mixing ratio. The nocturnal low-level jet at 03:30 UTC is simulated, but with a wind speed of 4 m s - 1 lower than observed. The monsoon flow is reproduced by the simulation with a correct wind direction up to approximately 1.5-km high. Then, the simulated wind exhibits a wrong northwesterly direction and does not capture the observed veering of the wind. Above 3.5-km altitude the simulated wind presents a northeasterly direction, whereas it is southeasterly according to the observations. The wind in the operational analysis of the ECMWF also turns wrongly at the same altitude, indicating that the wind bias in the simulation is largely influenced by the synoptic-scale winds imposed at the initial and boundary conditions. At 08:00 UTC, the simulation misses the wind direction in the first 500 m, probably because of some local event absent in the simulation. Above, the same wind pattern from 03:30 and 06:00 UTC prevails due to the wrong large-scale winds.
Overall, the stable atmosphere at Savè is well captured by the LES. The latter shows, however, a warm bias below 2 km and a dry bias above. The height of the lifting condensation level ( z L C L ) for an air parcel departing from the surface is also represented in Figure 11. It is computed using Equation (24) from [36]:
z L C L z s f c + 20 + T 5 ( 100 - R H ) ,
where z L C L and z s f c are given in m, T is in C and R H in %. z L C L has been calculated for an air parcel departing from the surface and is representative of the lower atmosphere. For the three soundings, z L C L is larger for the LES compared with the soundings (∼730 m average difference). This is due to the warm bias up to 2-km altitude, that partly explains why the simulated warm-rain clouds develop higher in the atmosphere.
Despite the punctual character of these wind profiles, the incorrect wind direction above ∼1.5-km altitude is persistent. Thus, the simulated warm-rain cells, with mean tops at 3.9 km, are wrongly advected from the northeast, as shown from Figure 9. However, in the observations they generally lie below, within the monsoon flow, and exhibit a southwesterly propagation. The observed and simulated mean top height of the warm-rain cells (2.3 and 3.9 km, respectively) correspond to the height where the wind shear is the largest (veering). The observed warm-rain cells develop throughout the depth of the moist monsoon flow whereas in the simulation they grow in the dry northeasterly wind sector. The latter is, at the same time, responsible for the simulated dry bias as air is advected from the drier Sahelian middle atmosphere. This notable discrepancy on the wind direction comes from the initial conditions as these soundings launched during 24 July 2016 were not assimilated in the Integrated Forecasting System of the ECMWF. Furthermore, the denial experiment carried by [37] for the DACCIWA campaign shows that, in the absence of radiosonde assimilation, the main biases in the ECMWF forecasts occurred in the lower atmosphere, especially on wind and temperature during the night and early morning.
We further investigate the properties inside the simulated warm-rain clouds. The mean vertical profiles (and variability) of cloud and rain water, equivalent potential temperature ( θ e ), turbulence kinetic energy (TKE) and vertical wind shear are shown in Figure 12. The individual profiles are calculated at the center of mass of the tracked warm-rain cells, but only with depths lower than 2 km (which have mean base and top at 2.5 and 3.9 km). The reason for using this reduced sample is that the warm-rain clouds observed by the cloud radar at Savè show tops around 2-km height (Figure 5).
The mean rain mixing ratio is positive at the surface, meaning that, in average, rainfall reaches the ground, although the warm-rain base (defined at 20 dBZ) is higher. It then increases until it reaches its maximum value of 0.07 g kg - 1 at 4-km altitude, and then decreases sharply and vanishes at 4.5 km. Note that the actual top of the warm-rain clouds can lie above 4 km. Indeed, there is a peak in cloud water content between 4 and 4.5-km high. Below, the cloud water mixing ratio diminishes drastically whereas rain is present, probably due to too fast autoconversion and accretion of cloud droplets. The profile of θ e shows a maximum value of 346 K at the surface, and then decreases down to a value around 336 K at 4-km high. This minimum value of θ e at 4-km high explains the too high tops found in the simulation with respect to the observations. The simulated clouds grew through the depth of the unstable layer. The largest values of TKE and wind shear are found in the boundary layer, with a peak at around 800 m. TKE is minimum between ∼1.2 and 2.4 km, above the boundary layer and below the warm-rain cell base. Then, larger values of TKE are found up to 4.5 km, due to in-cloud eddies and turbulence. TKE shows a peak of 0.035 m 2 s - 2 at around 3.5 km. Also at that height, a local maximum of 2.5 × 10 - 3 s - 1 is reported for the vertical wind shear. It is found within a 1-km thick layer just below the warm-rain cells top. This indicates, in agreement with Figure 11, that the warm-rain clouds are capped by a high-sheared layer.

6. Conclusions

Warm rain around Savè, Benin, in southern West Africa (SWA) is explored during the case of 24 July 2016, as part of the DACCIWA project. The event yielded 1.6 mm warm rain in Savè, which was one of the three supersites of the field campaign during summer 2016. This amount represents roughly half of the warm rain recorded at the supersite. Overall, total warm rain contributed to 1.35% of total rainfall during the two-month campaign. Several instruments deployed in Savè served to document this warm-rain episode, particularly the X-band rain radar deployed by the KIT. This radar can identify warm-rain features because it delineates the volume of precipitating hydrometeors and determines the top heights, allowing the characterization of the warm-rain cells. A 200-m grid spacing LES run with the Meso-NH model is also used to investigate warm rain, its contribution and its drivers.
The analysis of warm rain is performed from 03:00 to 08:00 UTC, when it was reported at Savè. The detection of warm-rain cells is based on reflectivity retrieved by the X-band radar and simulated by the model at high-temporal frequency. We apply a clustering algorithm to masked reflectivities (>20 dBZ) and determine the vertical development of each cluster or rain cell. If the clusters have tops below 4 km and diameters larger than 2 km they are defined as warm-rain cells. These are also tracked in time, which allows the quantification of their lifetime.
The 5-h analysis of warm-rain cells shows that 78% of the 1057 observed ones have diameters between 2 and 5 km, and only 5% have diameters larger than 10 km. Among the 225 warm-rain cells that were tracked at least twice, approximately two-thirds exist for 5 or 10 min, and approximately one-third exist for 15 min to 1 h. Only a residual percentage of tracks last longer than 1 h. Although the number of both instantaneous and tracked warm-rain cells is roughly half in the simulation, the model reasonably captures the size and lifetime relative frequency distributions.
The observed tracked warm-rain cells lie mainly in the western half of the radar domain. While most of the nighttime warm-rain cells (detected before 06:00 UTC) lie in the outer 50-km radar range and have higher tops, the daytime cells (after 06:00 UTC) are mostly found in the inner 50-km radar range and are shallower. During the night, warm rain is mainly linked to evolving phases of deep convection. During the day, it comes from pop-corn-like shallow convection. Most of the simulated tracked warm-rain cells are located in the eastern half of the domain. They trigger only after sunrise and extend 1.6 km higher in the atmosphere. The observed (simulated) top is 2.3 km (3.9 km), and the simulated base is 1.7 km. While in the observations the predominant direction of propagation is southwesterly, with most of the tracked warm-rain cells embedded in the 2-km thick monsoon flow, the simulated cells exhibit a northeasterly propagation. In other words, the model is able to produce warm rain during this case. It also reproduces, in a satisfactory manner, the size and lifetime of warm-rain cells, but has trouble simulating their spatial distribution, direction of propagation and vertical development.
The main drawbacks in the propagation and vertical extension in the model come from biases in the meteorological fields provided by the initial and boundary conditions. On the one hand, a 1–2 K warm bias in the two first kilometers with respect to the radiosounde measurements at Savè renders the simulated lower atmosphere more unstable and favors the vertical growth of the warm-rain clouds. On the other hand, the simulated wind direction above approximately 1.5-km altitude, where the simulated warm-rain cells propagate, is wrong. This failure is attributed to the wrong large-scale forcing from the ECMWF operational analysis.
The amount of warm rain is further computed in the simulation from the ground rainfall covered by the warm-rain cells. It is found that 1.2% of total rainfall is due to warm rain. Interestingly, the contribution to the area covered by warm rain is much larger (24.4%). This means that warm rain has the potential to distribute rainfall over large areas and suggests that warm rain might be of high importance during the drier periods that SWA experiences during the monsoon season [11]. Although the simulated 1.2% contribution to warm rain is only calculated for a 6-h period, it compares well in order of magnitude with the 1.35% contribution in Savè during the DACCIWA campaign, and with the 2% estimation by [9] in its climatology over the coastal regions of SWA.
This is the first time that warm-rain amounts are quantified from a LES over SWA. The results shown here remain, however, limited, as only one single warm-rain event was studied. Other cases should be investigated. Currently unknown tendencies and seasonal and interannual variability of warm rain in SWA could be addressed by applying the reliable yet inexpensive warm-rain cell detection method used here to longer simulations or radar databases. This study reflects that uncertainties are intrinsic to any data source, which contributes to the known issue of model-observation comparison. The comparison done here indicates that, while a 200-m grid permits shallow convection, proper large-scale conditions remain fundamental to entirely reproduce the warm-rain properties in limited-area models. Further observational efforts are needed to improve today’s operational analysis, probably in the form of continuous observations rather than intensive field operations.
This paper demonstrates the challenge that high-resolution simulations face to capture fine-scale processes leading to warm rain formation in southern West Africa. Because the large-scale forcing was not reliable enough, the representation of the physical processes in the model has not been tested here. Addressing the sensitivity of warm rain to warm-rain microphysics, in particular to the critical cloud water content beyond which autoconversion is initiated, could be a pathway for future research.

7. Data Availability

The DACCIWA data from the Savè supersite are available on the SEDOO database (http://baobab.sedoo.fr/DACCIWA/). The Meso-NH output data are available from I.R.M. upon request.

Author Contributions

I.R.M. designed and performed the simulation and analyses. J.-P.C. participated in the design of the simulation. J.H. provided the X-band radar and disdrometer data. I.R.M. prepared the manuscript with contributions from all co-authors. All authors have read and agreed to the published version of the manuscript.

Funding

The research leading to these results has received funding from the European Union Seventh Framework Program (FP7/2007-2013) under Grant 603502 (EU project DACCIWA: Dynamics-aerosol-chemistry-cloud interactions in West Africa).

Acknowledgments

Computer resources were allocated by GENCI through Project 90569. We thank Norbert Kalthoff (Karlsruhe Institute of Technology) for the ground-based campaign organization at Savè, and Fabienne Lohou (Laboratoire d’Aérologie) for the radiosonde measurements. We thank the anonymous referees for their comments, which helped to improve the overall quality of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mathon, V.; Laurent, H.; Lebel, T. Mesoscale convective system rainfall in the Sahel. J. Appl. Meteorol. 2002, 41, 1081–1092. [Google Scholar] [CrossRef]
  2. Fink, A.; Vincent, D.; Ermert, V. Rainfall types in the West African Sudanian zone during the summer monsoon 2002. Mon. Weather Rev. 2006, 134, 2143–2164. [Google Scholar] [CrossRef] [Green Version]
  3. Liu, C.; Zipser, E.J. “Warm rain” in the tropics: Seasonal and regional distributions based on 9 yr of TRMM data. J. Clim. 2009, 22, 767–779. [Google Scholar] [CrossRef] [Green Version]
  4. Chen, R.; Li, Z.; Kuligowski, R.J.; Ferraro, R.; Weng, F. A study of warm rain detection using A-Train satellite data. Geophys. Res. Lett. 2011, 38. [Google Scholar] [CrossRef] [Green Version]
  5. Luo, Z.; Liu, G.; Stephens, G.L.; Johnson, R.H. Terminal versus transient cumulus congestus: A CloudSat perspective. Geophys. Res. Lett. 2009, 36. [Google Scholar] [CrossRef] [Green Version]
  6. Houze, R.A.; Rasmussen, K.L.; Zuluaga, M.D.; Brodzik, S.R. The variable nature of convection in the tropics and subtropics: A legacy of 16 years of the Tropical Rainfall Measuring Mission satellite. Rev. Geophys. 2015, 53, 994–1021. [Google Scholar] [CrossRef]
  7. Mülmenstädt, J.; Sourdeval, O.; Delanoë, J.; Quaas, J. Frequency of occurrence of rain from liquid-, mixed-, and ice-phase clouds derived from A-Train satellite retrievals. Geophys. Res. Lett. 2015, 42, 6502–6509. [Google Scholar] [CrossRef] [Green Version]
  8. Takahashi, H.; Suzuki, K.; Stephens, G. Land–ocean differences in the warm-rain formation process in satellite and ground-based observations and model simulations. Q. J. R. Meteorol. Soc. 2017, 143, 1804–1815. [Google Scholar] [CrossRef]
  9. Maranan, M.; Fink, A.H.; Knippertz, P. Rainfall types over southern West Africa: Objective identification, climatology and synoptic environment. Q. J. R. Meteorol. Soc. 2018, 144, 1628–1648. [Google Scholar] [CrossRef] [Green Version]
  10. Young, M.P.; Chiu, J.C.; Williams, C.J.; Stein, T.H.; Stengel, M.; Fielding, M.D.; Black, E. Spatio-temporal variability of warm rain events over southern West Africa from geostationary satellite observations for climate monitoring and model evaluation. Q. J. R. Meteorol. Soc. 2018, 144, 2311–2330. [Google Scholar] [CrossRef] [Green Version]
  11. Adejuwon, J.O.; Odekunle, T.O. Variability and the Severity of the “Little Dry Season” in Southwestern Nigeria. J. Clim. 2006, 19, 483–493. [Google Scholar] [CrossRef]
  12. Zhu, P.; Albrecht, B. Large eddy simulations of continental shallow cumulus convection. J. Geophys. Res. Atmos. 2003, 108. [Google Scholar] [CrossRef]
  13. Seifert, A.; Heus, T. Large-eddy simulation of organized precipitating trade wind cumulus clouds. Atmos. Chem. Phys 2013, 13, 5631–5645. [Google Scholar] [CrossRef] [Green Version]
  14. Mechem, D.B.; Oberthaler, A. Numerical simulation of tropical cumulus congestus during TOGA COARE. J. Adv. Model. Earth Syst. 2013, 5, 623–637. [Google Scholar] [CrossRef]
  15. Endo, S.; Fridlind, A.M.; Lin, W.; Vogelmann, A.M.; Toto, T.; Ackerman, A.S.; McFarquhar, G.M.; Jackson, R.C.; Jonsson, H.H.; Liu, Y. RACORO continental boundary layer cloud investigations: 2. Large-eddy simulations of cumulus clouds and evaluation with in situ and ground-based observations. J. Geophys. Res. Atmos. 2015, 120, 5993–6014. [Google Scholar] [CrossRef] [Green Version]
  16. Knippertz, P.; Coe, H.; Chiu, J.C.; Evans, M.J.; Fink, A.H.; Kalthoff, N.; Liousse, C.; Mari, C.; Allan, R.P.; Brooks, B.; et al. The DACCIWA project: Dynamics–aerosol–chemistry–cloud interactions in West Africa. Bull. Am. Meteorol. Soc. 2015, 96, 1451–1460. [Google Scholar] [CrossRef] [Green Version]
  17. Kalthoff, N.; Lohou, F.; Brooks, B.; Jegede, G.; Adler, B.; Babić, K.; Dione, C.; Ajao, A.; Amekudzi, L.K.; Aryee, J.N.A.; et al. An overview of the diurnal cycle of the atmospheric boundary layer during the West African monsoon season: Results from the 2016 observational campaign. Atmos. Chem. Phys. 2018, 18, 2913–2928. [Google Scholar] [CrossRef] [Green Version]
  18. Flamant, C.; Knippertz, P.; Fink, A.H.; Akpo, A.; Brooks, B.; Chiu, C.J.; Coe, H.; Danuor, S.; Evans, M.; Jegede, O.; et al. The Dynamics–Aerosol–Chemistry–Cloud Interactions in West Africa Field Campaign: Overview and Research Highlights. Bull. Am. Meteorol. Soc. 2018, 99, 83–104. [Google Scholar] [CrossRef]
  19. Lac, C.; Chaboureau, P.; Masson, V.; Pinty, P.; Tulet, P.; Escobar, J.; Leriche, M.; Barthe, C.; Aouizerats, B.; Augros, C.; et al. Overview of the Meso-NH model version 5.4 and its applications. Geosci. Model Dev. 2018, 11, 1929–1969. [Google Scholar] [CrossRef] [Green Version]
  20. Stein, J.; Richard, E.; Lafore, J.P.; Pinty, J.; Asencio, N.; Cosma, S. High-resolution non-hydrostatic simulations of flash-flood episodes with grid-nesting and ice-phase parameterization. Meteorol. Atmos. Phys. 2000, 72, 203–221. [Google Scholar] [CrossRef]
  21. Lohou, F.; Kalthoff, N.; Adler, B.; Babić, K.; Dione, C.; Lothon, M.; Pedruzo-Bagazgoitia, X.; Zouzoua, M. Conceptual model of diurnal cycle of stratiform low-level clouds over southern West Africa. Atmos. Chem. Phys. Discuss. 2019, 2019, 1–25. [Google Scholar] [CrossRef] [Green Version]
  22. Masson, V.; Le Moigne, P.; Martin, E.; Faroux, S.; Alias, A.; Alkama, R.; Belamari, S.; Barbu, A.; Boone, A.; Bouyssel, F.; et al. The SURFEXv7.2 land and ocean surface platform for coupled or offline simulation of earth surface variables and fluxes. Geosci. Model Dev. 2013, 6, 929–960. [Google Scholar] [CrossRef] [Green Version]
  23. Mlawer, E.J.; Taubman, S.J.; Brown, P.D.; Iacono, M.J.; Clough, S.A. Radiative transfer for inhomogeneous atmospheres: RRTM, a validated correlated-k model for the longwave. J. Geophys. Res. Atmos. 1997, 102, 16663–16682. [Google Scholar] [CrossRef] [Green Version]
  24. Fouquart, Y.; Bonnel, B. Computations of solar heating of the earth’s atmosphere- A new parameterization. Beitraege Zur Phys. Der Atmos. 1980, 53, 35–62. [Google Scholar]
  25. Cuxart, J.; Bougeault, P.; Redelsperger, J.L. A turbulence scheme allowing for mesoscale and large-eddy simulations. Q. J. R. Meteorol. Soc. 2000, 126, 1–30. [Google Scholar] [CrossRef]
  26. Pergaud, J.; Masson, V.; Malardel, S.; Couvreux, F. A parameterization of dry thermals and shallow cumuli for mesoscale numerical weather prediction. Bound.-Layer Meteorol. 2009, 132, 83. [Google Scholar] [CrossRef]
  27. Chaboureau, J.P.; Bechtold, P. Statistical representation of clouds in a regional model and the impact on the diurnal cycle of convection during Tropical Convection, Cirrus and Nitrogen Oxides (TROCCINOX). J. Geophys. Res. Atmos. 2005, 110. [Google Scholar] [CrossRef] [Green Version]
  28. Pinty, J.; Jabouille, P. A mixed-phase cloud parameterization for use in mesoscale non-hydrostatic model: Simulations of a squall line and of orographic precipitations. In Proceedings of the Conference Cloud Physics, Everett, WA, USA, 17–21 August 1999; American Meteorological Society: Everett, WA, USA, 1999; pp. 217–220. [Google Scholar]
  29. Richard, E.; Cosma, S.; Tabary, P.; Pinty, J.P.; Hagen, M. High-resolution numerical simulations of the convective system observed in the Lago Maggiore area on 17 September 1999 (MAP IOP 2a). Q. J. R. Meteorol. Soc. 2003, 129, 543–563. [Google Scholar] [CrossRef]
  30. Reinares Martínez, I.; Chaboureau, J.P. Precipitation and mesoscale convective systems: Explicit versus parameterized convection over northern Africa. Mon. Weather Rev. 2018, 146, 797–812. [Google Scholar] [CrossRef]
  31. Short, D.A.; Nakamura, K. TRMM Radar Observations of Shallow Precipitation over the Tropical Oceans. J. Clim. 2000, 13, 4107–4124. [Google Scholar] [CrossRef]
  32. Machado, L.A.; Chaboureau, J.P. Effect of turbulence parameterization on assessment of cloud organization. Mon. Weather Rev. 2015, 143, 3246–3262. [Google Scholar] [CrossRef]
  33. Carr, N.; Kirstetter, P.; Gourley, J.; Hong, Y. Polarimetric signatures of midlatitude warm-rain precipitation events. J. Appl. Meteorol. Climatol. 2017, 56, 697–711. [Google Scholar] [CrossRef]
  34. Dione, C.; Lohou, F.; Lothon, M.; Adler, B.; Babić, K.; Kalthoff, N.; Pedruzo-Bagazgoitia, X.; Bezombes, Y.; Gabella, O. Low-level stratiform clouds and dynamical features observed within the southern West African monsoon. Atmos. Chem. Phys. 2019, 19, 8979–8997. [Google Scholar] [CrossRef] [Green Version]
  35. Adler, B.; Babić, K.; Kalthoff, N.; Lohou, F.; Lothon, M.; Dione, C.; Pedruzo-Bagazgoitia, X.; Andersen, H. Nocturnal low-level clouds in the atmospheric boundary layer over southern West Africa: An observation-based analysis of conditions and processes. Atmos. Chem. Phys. 2019, 19, 663–681. [Google Scholar] [CrossRef] [Green Version]
  36. Lawrence, M.G. The relationship between relative humidity and the dewpoint temperature in moist air: A simple conversion and applications. Bull. Am. Meteorol. Soc. 2005, 86, 225–234. [Google Scholar] [CrossRef]
  37. Van der Linden, R.; Knippertz, P.; Fink, A.H.; Ingleby, B.; Maranan, M.; Benedetti, A. The influence of DACCIWA radiosonde data on the quality of ECMWF analyses and forecasts. Q. J. R. Meteorol. Soc. 2020. in review. [Google Scholar]
Figure 1. Radar scan strategy for a single azimuth. The altitudes and distances scanned for each of the 12 elevation angles are shown. The shaded regions illustrate the uncertainty of the radar observation at 4-km altitude.
Figure 1. Radar scan strategy for a single azimuth. The altitudes and distances scanned for each of the 12 elevation angles are shown. The shaded regions illustrate the uncertainty of the radar observation at 4-km altitude.
Atmosphere 11 00298 g001
Figure 2. Meso-NH domain for the 200-m simulation, colour-shaded is the orography. The domain for the 1-km simulation is shown in the insert. The circles indicate the 100-km range of the X-band radar located at Savè.
Figure 2. Meso-NH domain for the 200-m simulation, colour-shaded is the orography. The domain for the 1-km simulation is shown in the insert. The circles indicate the 100-km range of the X-band radar located at Savè.
Atmosphere 11 00298 g002
Figure 3. Cumulative distribution function of the maximum 2-D reflectivities from 03:00 to 08:00 UTC 24 July 2016 at 5-min temporal resolution on a 1-km grid for the X-band radar and Meso-NH computed inside the radar range. For the radar, 16.9% of the reflectivity values lie above 20 dBZ, whereas only 8.6% of the simulated reflectivity values exceed 20 dBZ.
Figure 3. Cumulative distribution function of the maximum 2-D reflectivities from 03:00 to 08:00 UTC 24 July 2016 at 5-min temporal resolution on a 1-km grid for the X-band radar and Meso-NH computed inside the radar range. For the radar, 16.9% of the reflectivity values lie above 20 dBZ, whereas only 8.6% of the simulated reflectivity values exceed 20 dBZ.
Atmosphere 11 00298 g003
Figure 4. Left: horizontal projection of maximum reflectivity as retrieved from the X-band radar. Right: warm-rain cells detected with the clustering algorithm described in Section 2.3, labeled from 1 to 27. The inner and outer circles have 50- and 100-km radius, respectively.
Figure 4. Left: horizontal projection of maximum reflectivity as retrieved from the X-band radar. Right: warm-rain cells detected with the clustering algorithm described in Section 2.3, labeled from 1 to 27. The inner and outer circles have 50- and 100-km radius, respectively.
Atmosphere 11 00298 g004
Figure 5. Time evolution of (top) the reflectivity in the atmospheric column above Savè retrieved by the cloud radar and (bottom) the rain rate derived from a Joss-Waldvogel-Distrometer.
Figure 5. Time evolution of (top) the reflectivity in the atmospheric column above Savè retrieved by the cloud radar and (bottom) the rain rate derived from a Joss-Waldvogel-Distrometer.
Atmosphere 11 00298 g005
Figure 6. Maps of total accumulated rainfall from 03:00 to 09:00 UTC, for TRMM, product 3B42V7 (top left); for GPM, product IMERGV06B (top right); for Meso-NH (bottom left) and for Meso-NH but only the accumulated warm rain (bottom right). Note that the labels of the colorbars differ between the maps of total and warm rain.
Figure 6. Maps of total accumulated rainfall from 03:00 to 09:00 UTC, for TRMM, product 3B42V7 (top left); for GPM, product IMERGV06B (top right); for Meso-NH (bottom left) and for Meso-NH but only the accumulated warm rain (bottom right). Note that the labels of the colorbars differ between the maps of total and warm rain.
Atmosphere 11 00298 g006
Figure 7. Left: histogram of the effective diameter of the warm-rain cells detected by the X-band radar (black) and simulated by Meso-NH (red) on 24 July 2016 between 03:00 and 08:00 UTC. The effective diameter corresponds to the diameter of a circle covering the same area than the warm-rain cell. As a recall, warm-rain cells must have a diameter larger than 2 km. There is a total of 1057 observed and 485 simulated warm-rain cells. Right: histogram of the lifetime of the observed (black) and simulated (red) tracked warm-rain cells, i.e., those detected at least twice, using an overlapping threshold of 20%. There is a total of 225 observed and 107 simulated tracks.
Figure 7. Left: histogram of the effective diameter of the warm-rain cells detected by the X-band radar (black) and simulated by Meso-NH (red) on 24 July 2016 between 03:00 and 08:00 UTC. The effective diameter corresponds to the diameter of a circle covering the same area than the warm-rain cell. As a recall, warm-rain cells must have a diameter larger than 2 km. There is a total of 1057 observed and 485 simulated warm-rain cells. Right: histogram of the lifetime of the observed (black) and simulated (red) tracked warm-rain cells, i.e., those detected at least twice, using an overlapping threshold of 20%. There is a total of 225 observed and 107 simulated tracks.
Atmosphere 11 00298 g007
Figure 8. Time series of the number, effective diameter and depth (only for the simulation) of the tracked warm-rain cells, corresponding to those represented in Figure 9 and Figure 10. The dots show the mean value and the error bars represent the standard deviation for the effective diameter and depth.
Figure 8. Time series of the number, effective diameter and depth (only for the simulation) of the tracked warm-rain cells, corresponding to those represented in Figure 9 and Figure 10. The dots show the mean value and the error bars represent the standard deviation for the effective diameter and depth.
Atmosphere 11 00298 g008
Figure 9. Spatial distribution of the tracked warm-rain cells detected by the X-band rain radar (left) and simulated (right) from 03:00 to 08:00 UTC, at a 5-min temporal resolution. The dots represent the location of the centre of mass of the warm-rain cells; the solid dots show the first detection and the open circles the rest of the trajectory for each track. Colors indicate the time when the warm-rain cells are identified.
Figure 9. Spatial distribution of the tracked warm-rain cells detected by the X-band rain radar (left) and simulated (right) from 03:00 to 08:00 UTC, at a 5-min temporal resolution. The dots represent the location of the centre of mass of the warm-rain cells; the solid dots show the first detection and the open circles the rest of the trajectory for each track. Colors indicate the time when the warm-rain cells are identified.
Atmosphere 11 00298 g009
Figure 10. As in Figure 9 (only for the simulation), but colors represent the depth of the tracked warm-rain cells.
Figure 10. As in Figure 9 (only for the simulation), but colors represent the depth of the tracked warm-rain cells.
Atmosphere 11 00298 g010
Figure 11. Vertical profiles of, from top to bottom, potential temperature, water vapor mixing ratio, wind speed and direction on 24 July 2016 at 03:30 (left), 06:00 (centre) and 08:00 UTC (right), for the soundings at Savè, the Meso-NH LES and the ECMWF operational analysis (the latter only at 06:00 UTC). For the simulation and the ECMWF analysis the profiles are taken at the nearest gridpoint to the Savè location. The horizontal dashed lines in the top panels represent the height of the lifting condensation level (LCL).
Figure 11. Vertical profiles of, from top to bottom, potential temperature, water vapor mixing ratio, wind speed and direction on 24 July 2016 at 03:30 (left), 06:00 (centre) and 08:00 UTC (right), for the soundings at Savè, the Meso-NH LES and the ECMWF operational analysis (the latter only at 06:00 UTC). For the simulation and the ECMWF analysis the profiles are taken at the nearest gridpoint to the Savè location. The horizontal dashed lines in the top panels represent the height of the lifting condensation level (LCL).
Atmosphere 11 00298 g011
Figure 12. Vertical profiles of mean cloud and rain water, equivalent potential temperature ( θ e ), turbulence kinetic energy (TKE), and vertical wind shear inside the simulated tracked warm-rain cells with depths smaller than 2 km. The shaded contours correspond to the standard deviation among the individual profiles of the warm-rain cells. Individual profiles are taken at their center of mass.
Figure 12. Vertical profiles of mean cloud and rain water, equivalent potential temperature ( θ e ), turbulence kinetic energy (TKE), and vertical wind shear inside the simulated tracked warm-rain cells with depths smaller than 2 km. The shaded contours correspond to the standard deviation among the individual profiles of the warm-rain cells. Individual profiles are taken at their center of mass.
Atmosphere 11 00298 g012

Share and Cite

MDPI and ACS Style

Reinares Martínez, I.; Chaboureau, J.-P.; Handwerker, J. Warm Rain in Southern West Africa: A Case Study at Savè. Atmosphere 2020, 11, 298. https://doi.org/10.3390/atmos11030298

AMA Style

Reinares Martínez I, Chaboureau J-P, Handwerker J. Warm Rain in Southern West Africa: A Case Study at Savè. Atmosphere. 2020; 11(3):298. https://doi.org/10.3390/atmos11030298

Chicago/Turabian Style

Reinares Martínez, Irene, Jean-Pierre Chaboureau, and Jan Handwerker. 2020. "Warm Rain in Southern West Africa: A Case Study at Savè" Atmosphere 11, no. 3: 298. https://doi.org/10.3390/atmos11030298

APA Style

Reinares Martínez, I., Chaboureau, J. -P., & Handwerker, J. (2020). Warm Rain in Southern West Africa: A Case Study at Savè. Atmosphere, 11(3), 298. https://doi.org/10.3390/atmos11030298

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop