1. Introduction
The West African monsoon (WAM) brings seasonal rainfall to the Sahel, one of the most densely populated areas on the African continent. The summer monsoon is dominated by mesoscale convective systems (MCSs) which contribute 90% of the annual rainfall in the Sahel (Mathon et al. 2002). During the monsoon period, most regions of West Africa have a single diurnal peak of rainfall either in the afternoon or night, with a minimum of rainfall in the morning (Zhang et al. 2016). Afternoon rainfall peaks are associated with locally initiated storms, some of which sustain into the night, developing into MCSs that propagate westward (Taylor et al. 2010) to give nocturnal rainfall peaks (Duvel 1989; Zhang et al. 2016).
Despite the drastic improvements in global numerical weather prediction over the last 50 years (Bauer et al. 2015), forecasts remain poor over northern Africa with weather models often struggling to provide any real skill for rainfall forecasts, even for lead times of less than a day (Vogel et al. 2018). A patchy observational network hampers successful forecasting (Kniffka et al. 2020), but even during the 2006 African Monsoon Multidisciplinary Analysis campaign (Parker et al. 2008) (when models were initialized with a larger and more consistent set of observations), the gain in skill was lost after 24 h of forecast time (Agustí-Panareda et al. 2010). Accurately capturing the diurnal cycle in numerical weather and climate models is important for correctly representing the entire monsoon system including pressure gradients, winds, and the water budget (Marsham et al. 2013; Birch et al. 2014). In this paper, we explore how aspects of this diurnal cycle support the decay or strengthening of mature squall lines.
The convection parameterizations used in global models struggle to capture MCSs leading to systematic errors in both global numerical weather prediction and climate models (Lee et al. 2007; Moncrieff 2013). Over land surfaces, convection parameterization schemes typically 1) trigger too easily, with rainfall peaks too early in the day, 2) produce rain that is both too weak and too widespread, and 3) lack any representation of organization (Dai et al. 1999; Yang and Slingo 2001; Collier and Bowman 2004; Stephens et al. 2010). This problem is particularly clear in the Sahel due to the dominance of organized convective systems (Peters et al. 2019; Kendon et al. 2019; Pante and Knippertz 2019; Fitzpatrick et al. 2020).
Sahelian storms in convection-permitting models have more intense rainfall, diurnal cycles and lifetimes that are more realistic, and improved propagation (Crook et al. 2019). Berthou et al. (2019a,b) compared convection-permitting climate models against those with parameterized convection on a pan-African domain with 4 km horizontal grid spacing. They found the former drastically improved the distribution of precipitation over the Sahel, particularly through simulating more short-lasting intense rainfall events linked to MCSs, which increases the mean precipitation and gives a better representation of wet and dry spells. This intensification of rainfall occurs in regions of both increased and decreased mean rainfall.
In the summer months, the differential sensible heating produces a low-level pressure gradient between the Saharan heat low and the cooler Gulf of Guinea (see Fig. 1 in Lafore et al. 2010). The resultant regional monsoon circulation (lower-tropospheric southwesterly winds) brings cool, humid air into the continent. This added moisture increases the low-level equivalent potential temperature (Parker et al. 2005b). The shallow monsoon layer undercuts the drier and hotter Saharan air layer, creating a transition zone associated with significant convective inhibition (CIN) and a buildup of convective available potential energy (CAPE). The lower-tropospheric temperature gradient from south to north (or baroclinicity) produces the midtropospheric African easterly jet (AEJ; see Fig. 10 in Parker et al. 2005a). To summarize, the WAM generates an environment with high CAPE and strong vertical wind shear between the low-level southwesterlies and the midlevel AEJ, while the large convective inhibition CIN barrier delays the triggering of convection.
Several studies have indicated how the diurnal cycle is intrinsic to the WAM. Throughout the night, radiative cooling of the surface results in the formation of a shallow, stable boundary layer such that air aloft is decoupled from the surface drag. Consequently, a nocturnal low-level jet (LLJ) forms via an inertial oscillation, which attains its peak intensity between 0500 and 0700 UTC, and is centered near 925 hPa. This southwesterly LLJ produces strong vertical shear with respect to the AEJ at around 600 hPa (Abdou et al. 2010). During the morning, strong solar heating of the ground causes dry boundary layer (BL) convection and erosion of the stably stratified layers formed overnight. Dry air is entrained from above into the BL, and there is downward mixing of momentum from the LLJ and strong near-surface gusts (Knippertz 2008). By midday, deep BL mixing has produced drag on this deep layer of air, and the vertical gradients of momentum and potential temperature are reduced. The LLJ and AEJ both weaken throughout the day, with the LLJ weakening rapidly after the morning and reaching a minimum in the late afternoon and early evening (Abdou et al. 2010). The AEJ has smaller diurnal variations, but weakens to a minimum around 1800 UTC on average (Kalapureddy et al. 2010). Thus, the diurnal cycles of the LLJ and AEJ combine to produce a diurnal minimum in vertical wind shear in the late afternoon (Parker et al. 2005b). The buildup of CAPE throughout the day, when combined with a daily minimum in the CIN barrier in the early evening, is conducive for the triggering of storms which then move westward overnight with the general direction of the midlevel AEJ (Zhang et al. 2016). Storms can be triggered through low-level convergence caused by orography, land surface gradients (e.g., from soil moisture differences), and African easterly waves.
It is well known that line-perpendicular vertical wind shear exerts a controlling influence on convective organization (Browning and Ludlam 1962; Ludlam 1963; Rotunno et al. 1988). However, it is still debated to what extent the ambient wind shear plays a dynamic role by counteracting the horizontal vorticity that is generated baroclinically at the cold pool’s outflow boundary (Weisman 1992; Fovell and Dailey 1995; Weisman and Rotunno 2004; Parker and Johnson 2004; Bryan et al. 2006), rather than a thermodynamic role by modulating the inflow of moisture and CAPE into the storm (Alfaro 2017). In the latter thermodynamic effect—the layer lifting model of convection (LLMC)—the composition of the air flowing into the storm is considered. A strong LLJ would favor the storms, by increasing the fraction of low-level high-θe air compared to the total system-relative inflow into the storm. Similarly, diurnal variations in vertical wind shear as well as changes in cold pool intensity can, according to Rotunno et al. (1988) or RKW theory, alter the structure of the storm including the tilt and intensity of updrafts.
Parker (2008) used idealized simulations to investigate the effects of low-level nocturnal cooling on mature surface-based convective systems, including on where inflowing parcels originate. He found that as the boundary layer cooled and convective systems transitioned from surface based to elevated (inflowing air comes from the lowest 500 m initially then from above that altitude), the mechanism responsible for lifting inflowing parcels evolved from a cold pool to a trapped internal gravity wave. Interestingly, another distinct behavior emerges as the cold pool temperature deficit decreases toward zero. As the system slows, it remains surface based, but produces more intense lifting. Parker (2008) attributes this to the vorticity at the outflow boundary being better balanced by the vertical wind shear, and so largely offsetting the increasing CIN and decreasing CAPE of the parcels that are being lifted. French and Parker (2010) expanded on this study to discuss the relative roles of shear and inflow in idealized MCS simulations with a nocturnal LLJ. They note that changes in the system relative inflow correlate well with changes in total upward mass flux, while changes in shear seem to better describe the changes in updraft intensity.
To summarize, there is a clear diurnal cycle to the WAM, including within the boundary layer flow. MCSs trigger in the evening and persist overnight, with an overall Sahel rainfall minimum in the late morning (Birch et al. 2014). Past studies have shown that the LLJ produces vertical wind shear with the AEJ, and we expect this shear to support the nocturnal squall-line MCSs which characterize much of the Sahel. The morning erosion of the LLJ is similarly expected to reduce shear, and we hypothesize that this provides a less favorable environment for MCSs.
In this study, we consider representative profiles of the Sahel at different times of the day, and use idealized simulations of propagating mature squall lines to study whether the diurnal cycle of shear alone is significant for MCSs and how this compares with the role of the diurnal cycle in temperature and humidity (thermodynamics). In this idealized framework we test 1) whether the nocturnal LLJ supports nocturnal storms and 2) whether the LLMC can be used to explain the modeled characteristics of the MCSs. Finally, we consider what insights these tools and analysis can give us into the actual diurnal cycle of the WAM, which includes both the wind and thermodynamic changes. It is important to note that the simulations used in this study are idealized, and so we do not consider the diurnal cycle of synoptically driven convergence, radiation, or mesoscale flows from inhomogeneous surfaces. Our experiments allow us to isolate the role of the diurnally changing environmental profiles from these factors.
In section 2, we describe our methods, including the control profile and how it is adapted to represent different times in the diurnal cycle and the numerical model that we use. In section 3, we outline how we determine the height, strength, and speed of the simulated cold pool, and introduce the LLMC theory that we will be testing. Numerical results are presented in section 4: results from varying wind profiles alone are in section 4a, results varying thermodynamics alone are in section 4b, and those from varying both wind profiles and thermodynamics together are in section 4c. Throughout we use these results to test the LLMC. We discuss our results and conclude in section 5.
2. Method
a. Numerical model: “Cloud Model 1”
Three-dimensional numerical simulations of idealized squall lines were performed with the Cloud Model 1, version 18 (Bryan and Fritsch 2002). The advection of velocities and scalars is integrated with fifth-order horizontal and vertical advection schemes with implicit diffusion, with a Klemp–Wilhelmson time-splitting, vertically implicit, horizontally explicit pressure solver. The domain extends 1600 km × 60 km in the horizontal (across-squall-line, x, and along-squall-line, y, directions, respectively) with 1 km grid spacing. Horizontal boundaries are open in the across-line (west to east) direction and periodic in the along-line direction. The Durran–Klemp formulation (Durran and Klemp 1983) is employed for the open-flow boundaries, and allows gravity waves to exit the computational domain with minimal reflection. The periodic boundary conditions in the along-line direction allow the squall line to extend across the entire domain in this direction. In the vertical, the domain is 23.4 km deep. The vertical grid stretches from 25 m spacing at the surface to 500 m at 8.4 km and above, resulting in 62 vertical levels. The lower boundary condition is set to a flat, rigid no-slip surface, while a free-slip boundary condition is applied to the top of the domain, with a Rayleigh damping sponge layer above 20 km. Microphysics processes are parameterized with a version of the Morrison double-moment scheme (Morrison et al. 2005), which predicts number concentrations and mixing ratios of cloud droplets, cloud ice, rain, snow, and hail. No radiation and surface heat and moisture fluxes are included in simulations shown here, which aids the isolation of the role of shear.
To trigger squall-line convection, a negative temperature perturbation ΔT(x, z) is applied locally around x = x0 = 912 km and z = z0 = 3 km, but extending uniformly across the entire y domain. We use a locally sinusoidal perturbation given by ΔT(x, z) = −ΔT(x, z)cos2(πβ/2) when β < 1, where
b. Environmental profiles
The model is initialized with vertical profiles of potential temperature, water vapor mixing ratio, and the horizontal wind components (i.e., θ, qυ, u, υ) created from 6-hourly radiosonde measurements made during the AMMA field campaign (Parker et al. 2008) in July–August 2006, at Niamey in Niger (location 13°29′N, 2°10′E). The control profile shown in Fig. 1 is a pre-storm composite profile (Bickle et al. 2020). To produce an atmospheric profile representative of those in which squall lines are observed to occur, only soundings from the evening, night, and early morning (1800, 0000, and 0600 UTC)1 are considered, since squall lines usually form in the evening and persist overnight. These profiles are then filtered to include only those within a 12-h period prior to the passing of a cold pool over Niamey, as identified by Provod et al. (2016). The control profile is a composite of the top 25% pre-storm profiles (13 in total) with the highest near-surface CAPE. The profile allows persistent MCSs in this idealized environment with no imposed large-scale convergence.
The control profile in Fig. 1 has a surface CAPE of 1862 J kg−1 and CIN barrier of −52 J kg−1. There is a decrease in humidity from 750 to 600 hPa indicating the Saharan air layer and the profile is almost dry adiabatic over these levels. The hodograph shows southwesterly winds above the surface, with the nocturnal LLJ spanning from 960 to 900 hPa with a peak in westerlies at 950 hPa. Above 900 hPa there is a continuous decrease in westerly winds, while southerly winds sustain around 3.8 m s−1 until 800 hPa when they start to weaken. The AEJ spans from 725 hPa to a maximum at 600 hPa. It is also worth noting the easterly winds above 250 hPa, as these control the spread of the cloud anvil. If compared to the overall mean July–August profile for Niamey, the pre-storm control profile has higher near-surface CAPE, a lower near-surface CIN barrier, and a stronger westerly component of the LLJ, but a weaker AEJ at 600 hPa. The simulation produced with this profile is discussed in Bickle et al. (2020).
Diurnal cycle from AMMA profiles
Radiosonde measurements at 0600, 1200, 1800, and 0000 UTC have been used to produce both an overall mean profile as well as mean profiles at each of the 6-hourly intervals to generate diurnal anomalies. The diurnal perturbations attained for each time were added to the pre-storm control profile (Fig. 1) to yield four new kinematic (Fig. 2) and thermodynamic (Figs. 3 and 4) environments. The diurnal variations in wind and thermodynamic profiles are shown in Figs. 2–4. This approach was selected as it provided a pre-storm profile to be used for the control, which allowed an MCS to be sustained in our simulations while producing unbiased anomalies of the mean diurnal cycle in wind and thermodynamics.
Nocturnal radiative cooling and the nocturnal LLJ result in near-surface air that is cool and moist at 0600 UTC (Fig. 3a), as well as a stable BL with both a maximum in the surface CIN barrier of −90 J kg−1 and a daily minimum in CAPE for air at all levels (Fig. 4a). The LLJ is at its sampled diurnal maximum of 10 m s−1 (Fig. 2a) and the vertical shear between the level of the LLJ and the AEJ is also at a maximum. By 1200 UTC, surface heating has increased temperatures at the surface to around 30°C, resulting in a well-mixed, turbulent convective BL (Fig. 3b). Additionally, CAPE values from the surface to 850 hPa have increased, while the CIN magnitudes have decreased from the surface to 925 hPa (Figs. 4a,b). Entrainment has deepened and dried the BL, and zonal wind speeds within the LLJ have decreased to less than 5 m s−1 but are at a maximum both near the surface and within the AEJ (Figs. 2b and 4d).
The BL continues to deepen throughout the day, such that there is a minimum CIN barrier at all heights at 1800 UTC, while CAPE values above 950 hPa are at a maximum (Figs. 4a,b). The decrease in BL humidity produces a higher lifting condensation level of 845 hPa and results in lower levels of surface CAPE than at 1200 and 0000 UTC (Figs. 3c and 4c). The LLJ has continued to weaken to a diurnal minimum, leading to a diurnal minimum in the vertical shear between the LLJ and the AEJ (Fig. 4d). After sunset, the surface cools again so that by 0000 UTC a stable nocturnal boundary layer has formed. The increased stratification allows the LLJ, which brings cool moist air to the region, to strengthen and so increase the near-surface qυ (Fig. 3d). The increased qυ and lower temperature of near-surface air results in a lower lifting condensation level at 923 hPa, and an increase in CAPE near the surface to a daily maximum of 2177 J kg−1 (Fig. 4). However, the stratified BL also results in large (around −70 J kg−1) CIN values.
c. Experiments
Three sets of simulations were run to isolate and explore the contributions of the diurnal cycle in the wind and thermodynamic profiles to the development of squall-line MCSs. These included 1) the “Wind experiments,” where the thermodynamic profile was kept constant (pre-storm control) while the wind profile was varied to represent the 6-hourly intervals of 0600, 1200, 1800, and 0000 UTC, 2) the “Thermodynamic experiments,” where the wind profile was constant (control) and the thermodynamic profile varied to represent the 6-hourly intervals, and 3) the “Combined experiments,” where both the thermodynamic and wind profiles were varied to represent each of the 6-hourly intervals. Each experiment was repeated five times with different random noise applied to the initial temperature perturbation, and the average across these five simulations was used in quantitative analysis.
3. Diagnostics
In this section, we describe how we will calculate the height, strength, and front of the simulated cold pool from the buoyancy difference. The LLMC is also introduced, along with the implied definitions of scalings for particular quantities linked to the storm, which will be tested throughout the rest of the paper by comparison with model output.
a. Cold pool
In this study, the cold pool front is identified where the virtual θ perturbation from the original profile first becomes less than −0.1 K (working in the increasing x direction). The propagation speed (PS) of the system is defined as the speed of the cold pool front. This is consistent with Alfaro (2017); however, we consider the speed at 93 m, rather than at the surface. Due to the nocturnal stable boundary layer, there is a clearer signal-to-noise ratio at this height.
b. Layer-lifting model of convection
In this section, we describe the LLMC developed in Alfaro (2017). The LLMC provides a scaling for the strength of convective ascent and rainfall given initial (pre-storm) thermodynamic and dynamic conditions, and the squall-line propagation speed. We will explore whether such scalings, which are primarily based on an environmental profile ahead of a squall line, could be used to predict its relative intensity.
1) Vertical velocities
2) Precipitation
The LLMC calculates the inflow of moisture and air with CAPE into the squall line to produce scalings describing the intensity of an MCS. However, the LLMC does not predict the storm structure, area, differences in microphysics, or indeed the propagation speed (upon which all the LLMC scalings depend, via PS in the integrals). PRll will be compared against 1) the mean rainfall over the storm area and 2) the total rainfall across the domain.
4. Results
In section 4a we first detail experiments investigating how the diurnal wind cycle of the Sahel region can impact the severity and longevity of squall lines in the WAM. In section 4b, we consider the impact of the diurnal thermodynamic cycle alone, and finally in section 4c we consider the impact of the combined (wind plus thermodynamic) diurnal cycle. In each case, we consider whether the LLMC explains variations in storm updrafts and rainfall.
Line-averaged maximum updraft rates (
LLMC variables and ΔU calculated from the initial environmental profile of simulations.
a. Effects of varying the wind profile with constant thermodynamic profile
Vertical cross sections of
Layer lifting indices
The LLMC was tested by plotting wll against
Diagnostics calculated from simulation results.
The maximum of the zonal component of the AEJ (m s−1), the zonal component of the LLJ (m s−1), and the zonal wind shear ΔU = UAEJ − ULLJ (m s−1) for each environment. The average PS of squall lines (m s−1) as defined in section 2a, cmax (m s−1) the maximum of the theoretical cold pool speed as defined in Table 2 and the PS that maximizes wll [
The indices from Alfaro (2017) correlated well with
b. Effects of varying the thermodynamic profile, with a constant wind profile
Simulations were run where the thermodynamic profile was varied with the wind profile held constant. In all cases, the storms were sustained, although the 1800 UTC case reintensifies from hour 5 to hour 6 and the 0600 UTC case begins to die out at 6.25 h (Fig. 5).
Layer lifting indices
The indices from the LLMC correlated with bulk measures, particularly in-cloud upward mass flux (r = 0.74, Fig. 9c) but also Rtot (r = 0.53, Fig. 9d), for all except the 1800 UTC initialized profile. However, the correlation was poorer than for the wind experiments (Fig. 7). When Fig. 9c was replotted for the total domain upward mass flux then the correlation improved (r = 0.88, not shown). This improved correlation is predominantly due to the simulations initialized with the 1800 UTC profile which produces strong upward velocities below the cloud base (Fig. 10c). Thus, Rtot appears to scale on the in-cloud upward mass flux, which is consistent with a quasi-steady system.
The LLMC did not correlate well with
The 0600 UTC thermodynamic initialized ensemble mean has unexpectedly high
The rate of rainfall evaporation in the 0600 UTC case is just above 80% of that of the other profiles (Fig. 11b) resulting in a warmer and so less dense cold pool. Furthermore, the 0600 UTC low-level environment is the coolest of the profiles with the potential temperature near the surface about 10°C lower than for the 1800 UTC case (Figs. 3 and 12). Thus, the cold pool produced with the 0600 UTC thermodynamic profile has a reduced buoyancy deficit compared to the other cases and is therefore slower (Fig. 12) and shallower. The near-surface stable layer that exists at 0600 UTC (Figs. 3 and 4) combined with the weak temperature deficit between the cold pool and its environment produces a system that acts differently to the other cases. These differences include the much smoother wave like shape of the 0600 UTC cold pool (Fig. 10a). Additionally, the ascent at the leading edge of the cold pool is weak, with a broad region of weak ascent feeding the main region of ascent above the yellow dashed line where B0 > 0 and so CAPE is realized.
For the simulations initialized with the 1200, 1800, and 0000 UTC profiles there is a weak system-relative rear-to-front flow near the surface of the cold pool in the 10 km behind the gust front (Figs. 10b–d) which follows the definition of a density current in Crook and Moncrieff (1988). There is also clear dynamical forcing of air at the edge of the cold pool front from levels beneath the yellow dashed line (where it is still negatively buoyant relative to the initial profile) that enables it to reach its level of free convection. For the simulation initialized with the 0600 UTC profile the cold pool is weaker and has a different structure. The system-relative flow is front to back and so has some resemblance to the definition of a gravity wave in Crook and Moncrieff (1988) and is similar to what is described in Parker (2008). Notably a temperature deficit still exists between the environment and the cold pool (Fig. 12a). To summarize, although the LLMC explains Rtot in its current form it does not predict storm structure nor storm area and thus cannot predict average ascent and rainfall rates across the storm. The 0600 UTC case, which has a cold pool with a weaker temperature deficit and a more stable environmental profile, gives more intense rain and ascent due to a narrower storm. The weaker cold pool is also structurally different appearing more wave based rather than density current based.
Even considering Rtot, the 1800 UTC case has lower rainfall than expected from its value of PRll. This profile also produces the fastest moving cold pool (Table 3). To gain insight into differences between the storms, we considered their condensation and evaporation rates and found that for the 1800 UTC case the total condensation of vapor to cloud (Fig. 11a) is around 1.5 × 106 kg s−1 lower than expected from PRll when compared to the other cases. The lifting condensation level of near-surface air is at its highest at this time of the day. Additionally, the rate of cloud which evaporates in the 1800 UTC case is around 50% of the rate at which cloud mass becomes precipitation while it is nearer 40% for the other cases (Fig. 11b). A larger fraction of this precipitation then evaporates in the 1800 UTC case, about 0.75 × 106 kg s−1 more than would be expected from the rate of conversion of cloud to precipitation (Fig. 11c). We conclude that the high rates of rain and cloud evaporation and low rate of condensation explain the low values of Rtot in the 1800 UTC case, compared to the value of PRll. The inflow from above 800 hPa is larger for this case (Fig. 13) and so although entrainment is not well resolved in these simulations more convectively stable and drier air is expected to enter the storm which explains the higher cloud evaporation.
To summarize, the LLMC correlates well with total upward mass flux while deviations from Rtot are explained by the LLMC not considering rain evaporation. The LLMC does not correlate well with average values across the storm, including
We now consider what controls the variations in wll and PRll shown in Fig. 9. The 1800 and 0000 UTC thermodynamic cases have the greatest wll and upward mass flux (Fig. 9c) with the wll values explained by the high levels of CAPE in these profiles. The 0000 UTC case has the largest system-relative inflow of moisture (Fig. 13) which translates into high rainfall. Despite the 0600 UTC thermodynamic profile having enhanced moisture due to the presence of the LLJ overnight, the small values of PS (around 9.7 m s−1) ensure that the system-relative inflow of moisture is not as great as for other times.
c. Combined thermodynamic and wind experiments
Time-averaged properties for the three sets of ensemble-mean simulations (Wind, Thermodynamic, and Combined) are shown in Fig. 14. For both upward mass flux and rainfall, the thermodynamic change tends to dominate the wind change for the runs with the combined profile applied. The 1800, 0000, and 1200 UTC storms have greater upward mass flux and rainfall than the 0600 UTC storms (Fig. 14) and rainfall is at a maximum for the 1200 UTC case. The values of PS of the storms are similar to their thermodynamic experiment counterparts (Table 3).
Layer lifting indices
The in-cloud upward mass flux correlates with wll (Fig. 14a) but not as well as if the total upward mass flux is considered (not shown but with correlations of 0.81 for Wind runs, 0.90 for Thermodynamic runs, 0.81 for Combined runs, and 0.82 when all runs are used together with consistent lines of best fit between the three sets of simulations). The Rtot results are accordant with the role of evaporation discussed in section 4b for cases initialized with the 1800 UTC thermodynamics (red triangles and circles in Fig. 14). Similarly to the thermodynamic experiments, it was found that the indices from Alfaro (2017) did not correlate well with
The 0600 UTC combined case has low values for wll despite the strong nocturnal jet. If the total inflow is considered it can be seen that although the inflow of moisture is high (Fig. 15b) the total inflow of mass (Fig. 15a) is also high compared to the other profiles. Thus, the minimum inflow of CAPE at 0600 UTC (Fig. 15c) produces a low value of wll as the fraction of convectively unstable air as a proportion of total inflow is smaller. The case also has low bulk measures of in-cloud upward mass flux and Rtot compared to the other storms (Fig. 14). This could be a result of the low equivalent potential temperatures (Fig. 13) and consequently the low CAPE values which occur despite the higher levels of qυ produced by the overnight advection of moisture.
Recent observational studies have linked decadal trends in shear to decadal trends in extreme storm rainfall (Taylor et al. 2017). According to the LLMC it is not simply shear, but the system-relative inflows of mass, CAPE, and water vapor that determine the rainfall, which are themselves a function of shear. For the mean values of all runs,
The PS of storms is not predicted by the LLMC and the impact of varying PS on the scaling wll is explored in Fig. 16. The value of PS that maximizes wll appears mainly dependent on the wind profile and is generally slower but within 1 m s−1 of the speed of the AEJ (Table 3). Traveling close to the speed of the AEJ provides an enhanced inflow of low-level air compared to a minimum inflow from midlevels (i.e., the height of the AEJ). Any further increases in PS, i.e., greater than the AEJ, will increase the system-relative inflow into the front of the system of both moist low-level air with high CAPE and dry midlevel air with no CAPE. Thus, the optimum PS is close to the AEJ as further increases will cause a greater fraction of the total inflow to be from midlevels. From Fig. 16 and Table 3 it can be seen that storms do move within 2 m s−1 of this optimum PS but whether this is because the AEJ simply “steers” the storm or whether the storms tend toward this optimal speed is unclear and requires further study with more varied environmental profiles. Both the PS and the y-integrated maximum of the theoretical cold pool intensity cmax of the 0600 UTC combined case are closest to this optimal PS (Table 3).
5. Conclusions
Previous studies have shown that there are strong diurnal cycles in winds, vertical wind shear, and MCSs in the WAM. Idealized simulations were run with environmental profiles representative of different points in the day to isolate and study the effects of both combined and isolated diurnal variations in thermodynamics and shear. A particular motivation was to test whether the nocturnal LLJ supports MCSs overnight in the Sahel. Additionally, the LLMC and whether it provides useful indices for storm updraft speeds and rainfall was investigated.
We conclude that the LLJ winds do support nocturnal storms. Increasing the magnitude of the LLJ also increases the system-relative inflow of CAPE and water vapor. However, the variation in thermodynamics is large diurnally and so dominates when combined with that of the wind.
The LLMC was tested to see whether the relative inflow of moisture and convective instability could be used to explain rates of rainfall and updraft strength. The LLMC correlated well with bulk measures of total upward mass flux and the total condensation. However, as the LLMC neither predicts storm structure nor the area of the storm, LLMC indices correlated poorly with average ascent and rainfall rates averaged over the storm (
Rtot is not particularly well correlated with PRll and our analysis shows that this can largely be explained by variations in thermodynamics affecting microphysical processes. The LLMC does not account for the rate of evaporation or the fraction of cloud that becomes rain. The simulations initialized with the 1800 UTC profile had a deep, dry, and warm BL and much greater fractional evaporation of rainfall, reducing Rtot and producing a strong and fast cold pool, which led to low condensation for its PRll. Despite the LLMC accounting for moisture moving into the storm this case shows that more consideration of the fraction that will condense and reach the surface as precipitation is needed.
Storms where the cold pool intensity was comparatively weak were more compact resulting in a more efficient storm with higher values of
We conclude that the LLMC does provide useful scalings for squall-line upward mass fluxes and rainfall, with results showing these properties are strongly controlled by shear. However, differences in microphysics and storm organization are not accounted for in the LLMC and this is particularly clear when the thermodynamic profile is varied, rather than the wind profile alone.
Our idealized simulations allow us to isolate how key features of the diurnal cycle in thermodynamics and wind shear can control the intensity of squall lines. However, they are artificially triggered storms which exist in an unchanging environment. In reality, storms tend to form in regions of synoptic and mesoscale convergence (Birch et al. 2014) which vary diurnally (Vizy and Cook 2018), as does radiation, which both interacts with clouds directly and produces the surface fluxes which drive BL circulations. Nevertheless, our simulations provide insight into the real world indicating that the most supportive profiles for convection are 1800 and 1200 UTC and then 0000 UTC. The 0600 UTC profile is the least favorable, consistent with observed storms dying through the late night and morning. This is explained by the maximum in surface CAPE and the minimum in the CIN barrier at 1800 UTC and the reverse at 0600 UTC.
Our 1200 UTC profile supports high rainfall due to environmental wind features which include strong surface winds and a maximum in the AEJ, the latter of which results in stronger shear despite the LLJ weakening as the surface heats up and the boundary layer increases. This is in apparent disagreement with the observed minimum in Sahel-wide rain at this time. However, analysis of observational data used in Crook et al. (2019) shows that there is a daily minimum of about 10% in the probability of an active storm dissipating at 1100–1200 UTC, compared with a probability of 22% from 0500 to 0700 UTC (not shown). Thus, our results are consistent with the observed decline in rainfall from 0000 to 1100 UTC, with storms triggering in the evening, persisting into the night, and dying out through the night and early morning, with the fraction of storms which die out at 1200 UTC being much lower than 0600 UTC.
We have also shown that the LLJ does help support storms overnight. Thus, for convection parameterization schemes to account for the diurnal cycle, a feature they currently struggle with, they must capture the shear effects. Future work should consider a temporally varying environment which would reveal how the life cycle of storms depends on different environments.
All times refer to UTC with sunrise and sunset in Niamey (where local time is UTC + 1 h) occurring at approximately 0530 and 1815 UTC, respectively.
Acknowledgments.
Thanks are due to Douglas Parker, whose advice and suggestions were invaluable. We thank George Bryan for making the CM1 model freely available (https://www2.mmm.ucar.edu/people/bryan/cm1/). This work also relied on efforts by others including the AMMA field campaign (Parker et al. 2008). This work was undertaken on ARC3, part of the High Performance Computing facilities at the University of Leeds, United Kingdom. We thank our funders who include EPSRC, University of Leeds CDT in Fluid Dynamics (EP/L01615X/1); Future Climate for Africa AMMA2050 and IMPALA projects (NE/M020126/1); GCRF African SWIFT (NE/P021077/1); National Centre for Atmospheric Science via the NERC/GCRF program, Atmospheric hazards in developing countries: risk assessment and early warning; NERC GENESIS and ParaCon programme (NE/N013840/1, NE/T003898/1). We thank the three anonymous reviewers, whose detailed and constructive comments have significantly improved the manuscript.
Data availability statement.
Radiosonde data used in this research were generated in Parker et al. (2008).
REFERENCES
Abdou, K., D. J. Parker, B. Brooks, N. Kalthoff, and T. Lebel, 2010: The diurnal cycle of lower boundary-layer wind in the West African monsoon. Quart. J. Roy. Meteor. Soc., 136 (Suppl. 1), 66–76, https://doi.org/10.1002/qj.536.
Agustí-Panareda, A., A. Beljaars, C. Cardinali, I. Genkova, and C. Thorncroft, 2010: Impacts of assimilating AMMA soundings on ECMWF analyses and forecasts. Wea. Forecasting, 25, 1142–1160, https://doi.org/10.1175/2010WAF2222370.1.
Alfaro, D. A., 2017: Low-tropospheric shear in the structure of squall lines: Impacts on latent heating under layer-lifting ascent. J. Atmos. Sci., 74, 229–248, https://doi.org/10.1175/JAS-D-16-0168.1.
Alfaro, D. A., and M. Khairoutdinov, 2015: Thermodynamic constraints on the morphology of simulated midlatitude squall lines. J. Atmos. Sci., 72, 3116–3137, https://doi.org/10.1175/JAS-D-14-0295.1.
Bauer, P., A. Thorpe, and G. Brunet, 2015: The quiet revolution of numerical weather prediction. Nature, 525, 47–55, https://doi.org/10.1038/nature14956.
Benjamin, T. B., 1968: Gravity currents and related phenomena. J. Fluid Mech., 31, 209–248, https://doi.org/10.1017/S0022112068000133.
Berthou, S., E. J. Kendon, D. P. Rowell, M. J. Roberts, S. Tucker, and R. A. Stratton, 2019a: Larger future intensification of rainfall in the West African Sahel in a convection-permitting model. Geophys. Res. Lett., 46, 13 299–13 307, https://doi.org/10.1029/2019GL083544.
Berthou, S., D. P. Rowell, E. J. Kendon, M. J. Roberts, R. A. Stratton, J. A. Crook, and C. Wilcox, 2019b: Improved climatological precipitation characteristics over West Africa at convection-permitting scales. Climate Dyn., 53, 1991–2011, https://doi.org/10.1007/s00382-019-04759-4.
Bickle, M. E., J. H. Marsham, A. N. Ross, D. P. Rowell, D. J. Parker, and C. M. Taylor, 2020: Understanding mechanisms for trends in Sahelian squall lines: Roles of thermodynamics and shear. Quart. J. Roy. Meteor. Soc., 147, 983–1006, https://doi.org/10.1002/qj.3955.
Birch, C. E., D. J. Parker, J. H. Marsham, D. Copsey, and L. Garcia-Carreras, 2014: A seamless assessment of the role of convection in the water cycle of the West African monsoon. J. Geophys. Res. Atmos., 119, 2890–2912, https://doi.org/10.1002/2013JD020887.
Browning, K. A., and F. H. Ludlam, 1962: Airflow in convective storms. Quart. J. Roy. Meteor. Soc., 88, 117–135, https://doi.org/10.1002/qj.49708837602.
Bryan, G. H., and J. M. Fritsch, 2002: A benchmark simulation for moist nonhydrostatic numerical models. Mon. Wea. Rev., 130, 2917–2928, https://doi.org/10.1175/1520-0493(2002)130<2917:ABSFMN>2.0.CO;2.
Bryan, G. H., J. C. Knievel, and M. D. Parker, 2006: A multimodel assessment of RKW theory’s relevance to squall-line characteristics. Mon. Wea. Rev., 134, 2772–2792, https://doi.org/10.1175/MWR3226.1.
Collier, J. C., and K. P. Bowman, 2004: Diurnal cycle of tropical precipitation in a general circulation model. J. Geophys. Res., 109, D17105, https://doi.org/10.1029/2004JD004818.
Crook, J., C. Klein, S. Folwell, C. M. Taylor, D. J. Parker, R. Stratton, and T. Stein, 2019: Assessment of the representation of West African storm lifecycles in convection-permitting simulations. Earth Space Sci., 6, 818–835, https://doi.org/10.1029/2018EA000491.
Crook, N. A., and M. W. Moncrieff, 1988: The effect of large-scale convergence on the generation and maintenance of deep moist convection. J. Atmos. Sci., 45, 3606–3624, https://doi.org10.1175/1520-0469(1988)045<3606:TEOLSC>2.0.CO;2.
Dai, A., F. Giorgi, and K. E. Trenberth, 1999: Observed and model-simulated diurnal cycles of precipitation over the contiguous United States. J. Geophys. Res., 104, 6377–6402, https://doi.org/10.1029/98JD02720.
Durran, D. R., and J. B. Klemp, 1983: A compressible model for the simulation of moist mountain waves. Mon. Wea. Rev., 111, 2341–2361, https://doi.org/10.1175/1520-0493(1983)111<2341:ACMFTS>2.0.CO;2.
Duvel, J. P., 1989: Convection over tropical Africa and the Atlantic Ocean during northern summer. Part I: Interannual and diurnal variations. Mon. Wea. Rev., 117, 2782–2799, https://doi.org/10.1175/1520-0493(1989)117<2782:COTAAT>2.0.CO;2.
Fitzpatrick, R. G., and Coauthors, 2020: What drives the intensification of mesoscale convective systems over the West African Sahel under climate change? J. Climate, 33, 3151–3172, https://doi.org/10.1175/JCLI-D-19-0380.1.
Fovell, R. G., and P. S. Dailey, 1995: The temporal behavior of numerically simulated multicell-type storms. Part I: Models of behavior. J. Atmos. Sci., 52, 2073–2095, https://doi.org/10.1175/1520-0469(1995)052<2073:ttbons>2.0.co;2.
French, A. J., and M. D. Parker, 2010: The response of simulated nocturnal convective systems to a developing low-level jet. J. Atmos. Sci., 67, 3384–3408, https://doi.org/10.1175/2010JAS3329.1.
Kalapureddy, M. C., M. Lothon, B. Campistron, F. Lohou, and F. Saïd, 2010: Wind profiler analysis of the African easterly jet in relation with the boundary layer and the Saharan heat-low. Quart. J. Roy. Meteor. Soc., 136 (Suppl. 1), 77–91, https://doi.org/10.1002/qj.494.
Kendon, E. J., R. A. Stratton, S. Tucker, J. H. Marsham, S. Berthou, D. P. Rowell, and C. A. Senior, 2019: Enhanced future changes in wet and dry extremes over Africa at convection-permitting scale. Nat. Commun., 10, 1794, https://doi.org/10.1038/s41467-019-09776-9.
Kniffka, A., and Coauthors, 2020: An evaluation of operational and research weather forecasts for southern West Africa using observations from the DACCIWA field campaign in June–July 2016. Quart. J. Roy. Meteor. Soc., 146, 1121–1148, https://doi.org/10.1002/qj.3729.
Knippertz, P., 2008: Dust emissions in the West African heat trough—The role of the diurnal cycle and of extratropical disturbances. Meteor. Z., 17, 553–563, https://doi.org/10.1127/0941-2948/2008/0315.
Lafore, J. P., C. Flamant, V. Giraud, F. Guichard, P. Knippertz, J. F. Mahfouf, P. Mascart, and E. R. Williams, 2010: Introduction to the AMMA Special Issue on ‘Advances in understanding atmospheric processes over West Africa through the AMMA field campaign’. Quart. J. Roy. Meteor. Soc., 136 (Suppl. 1), 2–7, https://doi.org/10.1002/qj.583.
Lee, M.-I., and Coauthors, 2007: An analysis of the warm-season diurnal cycle over the continental United States and northern Mexico in general circulation models. J. Hydrometeor., 8, 344–366, https://doi.org/10.1175/JHM581.1.
Ludlam, F. H., 1963: Severe local storms: A review. Severe Local Storms, Meteor. Monogr., No. 27, Amer. Meteor. Soc., 1–32, https://doi.org/10.1007/978-1-940033-56-3_1.
Marsham, J. H., N. S. Dixon, L. Garcia-Carreras, G. M. Lister, D. J. Parker, P. Knippertz, and C. E. Birch, 2013: The role of moist convection in the West African monsoon system: Insights from continental-scale convection-permitting simulations. Geophys. Res. Lett., 40, 1843–1849, https://doi.org/10.1002/grl.50347.
Mathon, V., H. Laurent, and T. Lebel, 2002: Mesoscale convective system rainfall in the Sahel. J. Appl. Meteor., 41, 1081–1092, https://doi.org/10.1175/1520-0450(2002)041<1081:MCSRIT>2.0.CO;2.
Moncrieff, M. W., 2013: The multiscale organization of moist convection and the intersection of weather and climate. Climate Dynamics: Why Does Climate Vary, Geophys. Monogr., Vol. 189, Amer. Geophys. Union, 3–26, https://doi.org/10.1029/2008GM000838.
Morrison, H., J. A. Curry, and V. I. Khvorostyanov, 2005: A new double-moment microphysics parameterization for application in cloud and climate models. Part I: Description. J. Atmos. Sci., 62, 1665–1677, https://doi.org/10.1175/JAS3446.1.
Pante, G., and P. Knippertz, 2019: Resolving Sahelian thunderstorms improves mid-latitude weather forecasts. Nat. Commun., 10, 3487, https://doi.org/10.1038/s41467-019-11081-4.
Parker, D. J., C. D. Thorncroft, R. R. Burton, and A. Diongue-Niang, 2005a: Analysis of the African easterly jet, using aircraft observations from the JET2000 experiment. Quart. J. Roy. Meteor. Soc., 131, 1461–1482, https://doi.org/10.1256/qj.03.189.
Parker, D. J., and Coauthors, 2005b: The diurnal cycle of the West African monsoon circulation. Quart. J. Roy. Meteor. Soc., 131, 2839–2860, https://doi.org/10.1256/qj.04.52.
Parker, D. J., and Coauthors, 2008: The AMMA radiosonde program and its implications for the future of atmospheric monitoring over Africa. Bull. Amer. Meteor. Soc., 89, 1015–1027, https://doi.org/10.1175/2008BAMS2436.1.
Parker, M. D., 2008: Response of simulated squall lines to low-level cooling. J. Atmos. Sci., 65, 1323–1341, https://doi.org/10.1175/2007JAS2507.1.
Parker, M. D., and R. H. Johnson, 2004: Structures and dynamics of quasi-2D mesoscale convective systems. J. Atmos. Sci., 61, 545–567, https://doi.org/10.1175/1520-0469(2004)061<0545:SADOQM>2.0.CO;2.
Peters, K., C. Hohenegger, and D. Klocke, 2019: Different representation of mesoscale convective systems in convection-permitting and convection-parameterizing NWP models and its implications for large-scale forecast evolution. Atmosphere, 10, 503, https://doi.org/10.3390/atmos10090503.
Provod, M., J. H. Marsham, D. J. Parker, and C. E. Birch, 2016: A characterization of cold pools in the West African Sahel. Mon. Wea. Rev., 144, 1923–1934, https://doi.org/10.1175/MWR-D-15-0023.1.
Rotunno, R., J. B. Klemp, and M. L. Weisman, 1988: A theory for strong, long-lived squall lines. J. Atmos. Sci., 45, 463–485, https://doi.org/10.1175/1520-0469(1988)045<0463:ATFSLL>2.0.CO;2.
Stephens, G. L., and Coauthors, 2010: Dreary state of precipitation in global models. J. Geophys. Res., 115, D24211, https://doi.org/10.1029/2010JD014532.
Taylor, C. M., P. P. Harrisa, and D. J. Parker, 2010: Impact of soil moisture on the development of a Sahelian mesoscale convective system: A case-study from the AMMA special observing period. Quart. J. Roy. Meteor. Soc., 136 (Suppl. 1), 456–470, https://doi.org/10.1002/qj.465.
Taylor, C. M., and Coauthors, 2017: Frequency of extreme Sahelian storms tripled since 1982 in satellite observations. Nature, 544, 475–478, https://doi.org/10.1038/nature22069.
Vizy, E. K., and K. H. Cook, 2018: Mesoscale convective systems and nocturnal rainfall over the West African Sahel: Role of the inter-tropical front. Climate Dyn., 50, 587–614, https://doi.org/10.1007/s00382-017-3628-7.
Vogel, P., P. Knippertz, A. H. Fink, A. Schlueter, and T. Gneiting, 2018: Skill of global raw and postprocessed ensemble predictions of rainfall over northern tropical Africa. Wea. Forecasting, 33, 369–388, https://doi.org/10.1175/WAF-D-17-0127.1.
Weisman, M. L., 1992: The role of convectively generated rear-inflow jets in the evolution of long-lived mesoconvective systems. J. Atmos. Sci., 49, 1826–1847, https://doi.org/10.1175/1520-0469(1992)049<1826:TROCGR>2.0.CO;2.
Weisman, M. L., and R. Rotunno, 2004: “A theory for strong long-lived squall lines” revisited. J. Atmos. Sci., 61, 361–382, https://doi.org/10.1175/1520-0469(2004)061<0361:ATFSLS>2.0.CO;2.
Yang, G.-Y., and J. Slingo, 2001: The diurnal cycle in the tropics. Mon. Wea. Rev., 129, 784–801, https://doi.org/10.1175/1520-0493(2001)129<0784:TDCITT>2.0.CO;2.
Zhang, G., K. H. Cook, and E. K. Vizy, 2016: The diurnal cycle of warm season rainfall over West Africa. Part I: Observational analysis. J. Climate, 29, 8423–8437, https://doi.org/10.1175/JCLI-D-15-0874.1.