Sensitivity of a Simulated Squall Line to Horizontal Resolution and Parameterization of Microphysics

2012 ◽  
Vol 140 (1) ◽  
pp. 202-225 ◽  
Author(s):  
George H. Bryan ◽  
Hugh Morrison

Abstract Idealized simulations of the 15 May 2009 squall line from the Second Verification of the Origins of Rotation in Tornadoes Experiment (VORTEX2) are evaluated in this study. Four different microphysical setups are used, with either single-moment (1M) or double-moment (2M) microphysics, and either hail or graupel as the dense (rimed) ice species. Three different horizontal grid spacings are used: Δx = 4, 1, or 0.25 km (with identical vertical grids). Overall, results show that simulated squall lines are sensitive to both microphysical setup and horizontal resolution, although some quantities (i.e., surface rainfall) are more sensitive to Δx in this study. Simulations with larger Δx are slower to develop, produce more precipitation, and have higher cloud tops, all of which are attributable to larger convective cells that do not entrain midlevel air. The highest-resolution simulations have substantially more cloud water evaporation, which is partly attributable to the development of resolved turbulence. For a given Δx, the 1M simulations produce less rain, more intense cold pools, and do not have trailing stratiform precipitation at the surface, owing to excessive rainwater evaporation. The simulations with graupel as the dense ice species have unrealistically wide convective regions. Comparison against analyses from VORTEX2 data shows that the 2M setup with hail and Δx = 0.25 km produces the most realistic simulation because (i) this simulation produces realistic distributions of reflectivity associated with convective, transition, and trailing stratiform regions, (ii) the cold pool properties are reasonably close to analyses from VORTEX2, and (iii) relative humidity in the cold pool is closest to observations.

2020 ◽  
Author(s):  
Ludovic Touzé-Peiffer ◽  
Nicolas Rochetin ◽  
Raphaela Vogel

<p>A considerable amount of literature has been devoted to the study of strong convective squall line. In particular, many studies have noted the role of cold pools on the persistence of these squall lines. Observations and simulations have shown that squall lines are often associated with pools of air cooled by partial rain evaporation. Such cold pools spread at the surface and may initiate new convective cells at their edges, thus contributing to the maintenance of a squall line. Under which environmental conditions the lifting at the edges of cold pools is most efficient has been subject to many debates. Yet, it is generally acknowledged that the environmental wind shear is a critical factor in this process. </p><p>Recent observations and realistic simulations over the trade-wind region have revealed persistent structures of shallow cumuli associated with surface cold pools. We will call these structures shallow convective squall lines, due to their similarity with strong convective squall lines. Based on simulations from the German model ICON and on recent observations from the field campaign EUREC4A, we will study the characteristics of these shallow convective squall lines and their lifecycle. Similarly to strong convective squall lines, shallow convective squall lines organized around a leading edge composed by many updrafts and downdrafts feeding the surface cold pools. We will see that the environmental wind shear plays a key role in the persistence of these shallow convective squall line, and we will compare our findings with classical theories for strong convective squall lines.</p>


2014 ◽  
Vol 71 (7) ◽  
pp. 2733-2746 ◽  
Author(s):  
Xiaowen Tang ◽  
Wen-Chau Lee ◽  
Michael Bell

Abstract This study examines the structure and dynamics of Typhoon Hagupit’s (2008) principal rainband using airborne radar and dropsonde observations. The convection in Hagupit’s principal rainband was organized into a well-defined line with trailing stratiform precipitation on the inner side. Individual convective cells had intense updrafts and downdrafts and were aligned in a wavelike pattern along the line. The line-averaged vertical cross section possessed a slightly inward-tilting convective core and two branches of low-level inflow feeding the convection. The result of a thermodynamic retrieval showed a pronounced cold pool behind the convective line. The horizontal and vertical structures of this principal rainband show characteristics that are different than the existing conceptual model and are more similar to squall lines and outer rainbands. The unique convective structure of Hagupit’s principal rainband was associated with veering low-level vertical wind shear and large convective instability in the environment. A quantitative assessment of the cold pool strength showed that it was quasi balanced with that of the low-level vertical wind shear. The balanced state and the structural characteristics of convection in Hagupit’s principal rainband were dynamically consistent with the theory of cold pool dynamics widely applied to strong and long-lived squall lines. The analyses suggest that cold pool dynamics played a role in determining the principal rainband structure in addition to storm-scale vortex dynamics.


2020 ◽  
Vol 77 (12) ◽  
pp. 4143-4170
Author(s):  
Kelly Lombardo

AbstractIdealized 3D numerical simulations are used to quantify the impact of moving marine atmospheric boundary layers (MABLs) on squall lines in an environment representative of the U.S. mid-Atlantic coastal plain. Characteristics of the MABL, including depth and potential temperature, are varied. Squall lines are most intense while moving over the deepest MABLs, while the storm encountering no MABL is the weakest. Storm intensity is only sensitive to MABL temperature when the MABL is sufficiently deep. Collisions between the storm cold pools and MABLs transition storm lift from surface-based cold pools to wavelike features, with the resulting ascent mechanism dependent on MABL density, not depth. Bores form when the MABL is denser than the cold pool and hybrid cold pool–bores form when the densities are similar. While these features support storms over the MABL, the type of lifting mechanism does not control storm intensity alone. Storm intensity depends on the amplification and maintenance of these features, which is determined by the ambient conditions. Isolated convective cells form ahead of squall lines prior to the cold pool–MABL collision, resulting in a rain peak and the eventual discrete propagation of the storms. Cells form as storm-generated high-frequency gravity waves interact with gravity waves generated by the moving marine layers, in the presence of reduced stability by the squall line itself. No cells form in the presence of the storm or the MABL alone.


2014 ◽  
Vol 53 (8) ◽  
pp. 2072-2090 ◽  
Author(s):  
Charlotte E. Wainwright ◽  
Daniel T. Dawson ◽  
Ming Xue ◽  
Guifu Zhang

AbstractIn this study, power-law relations are developed between the intercept parameter N0 of the exponential particle size distribution and the water content for the rain, hail, graupel, and snow hydrometeor categories within the Milbrandt and Yau microphysics scheme. Simulations of the 3 May 1999 Oklahoma tornadic supercell are performed using the diagnostic relations for rain only and alternately for all four precipitating species, and results are compared with those from the original fixed-N0 single- and double-moment versions of the scheme. Diagnosing N0 for rain is found to improve the results of the simulation in terms of reproducing the key features of the double-moment simulation while still retaining the computational efficiency of a single-moment scheme. Results more consistent with the double-moment scheme are seen in the general storm structure, the cold-pool structure and intensity, and the number concentration fields. Diagnosing the intercept parameters for all four species, including those for the ice species, within the single-moment scheme yields even closer agreement with the double-moment simulation results. The decreased cold-pool intensity is very similar to that produced by the double-moment simulation, as is the areal extent of the simulated storm. The diagnostic relations are also tested on a simulated squall line, with similar promising results. This study suggests that, when compared with traditional fixed intercept parameters used in typical single-moment microphysics schemes, results closer to a double-moment scheme can be obtained through the use of diagnostic relations for the parameters of the particle size distribution, with little extra computational cost.


2021 ◽  
Author(s):  
Sophie Abramian ◽  
Caroline Muller ◽  
Camille Risi

<p>Investigating tropical squall lines with a cloud resolving model</p><p>Using a cloud resolving model, we attempt to clarify the physical processes responsible for the organization of deep clouds into squall lines in the tropics. To do so, we impose a vertical wind shear, and investigate the response of deep convection to different shear strengths in radiative convective equilibrium. As the magnitude of the shear increases, the convection becomes more and more organized into a line, perpendicular to the shear. It is due to the interaction of the low-level shear with the cold pools associated with convective downdrafts. Beyond a certain shear, called optimal shear, the line tends to orient at an angle to the shear. The existing literature suggests that this angle conserves the projection of the shear on the direction perpendicular to the squall line near the optimal value, a hypothesis that we further investigate here.</p><p>In this work, we propose a systematic method, based on image auto-correlation, to determine the angle of the squall line with respect to the shear. We highlight the existence of the sub-critical and super-critical regime, as predicted by earlier studies. In the sub-critical regime, squall lines are indeed perpendicular to the shear. Yet, angles of squall lines in the super-critical regime do not clearly correspond to the conservation of the projected component of the shear near the optimal value. In particular, squall lines often remain more perpendicular to the shear than expected.</p><p>We thus investigate the balance between shear and cold pool winds to explain this difference. Using statistical methods on extreme events, we find that this difference is due to an intensification of cold pool potential energy with shear. Cold pool intensification allows the squall line to better resist to the shear, and thus reduces its angle of orientation. This new feature leads us to conclude that two mechanisms maintain a squall line in wind shear : the orientation of clouds and the intensification of cold pools.</p>


2020 ◽  
Vol 148 (12) ◽  
pp. 4971-4994
Author(s):  
McKenna W. Stanford ◽  
Hugh Morrison ◽  
Adam Varble

AbstractThis study investigates impacts of altering subgrid-scale mixing in “convection-permitting” kilometer-scale horizontal-grid-spacing (Δh) simulations by applying either constant or stochastic multiplicative factors to the horizontal mixing coefficients within the Weather Research and Forecasting Model. In quasi-idealized 1-km Δh simulations of two observationally based squall-line cases, constant enhanced mixing produces larger updraft cores that are more dilute at upper levels, weakens the cold pool, rear-inflow jet, and front-to-rear flow of the squall line, and degrades the model’s effective resolution. Reducing mixing by a constant multiplicative factor has the opposite effect on all metrics. Completely turning off parameterized horizontal mixing produces bulk updraft statistics and squall-line mesoscale structure closest to an LES “benchmark” among all 1-km simulations, although the updraft cores are too undilute. The stochastic mixing scheme, which applies a multiplicative factor to the mixing coefficients that varies stochastically in time and space, is employed at 0.5-, 1-, and 2-km Δh. It generally reduces midlevel vertical velocities and enhances upper-level vertical velocities compared to simulations using the standard mixing scheme, with more substantial impacts at 1- and 2-km Δh compared to 0.5-km Δh. The stochastic scheme also increases updraft dilution to better agree with the LES for one case, but has less impact on the other case. Stochastic mixing acts to weaken the cold pool but without a significant impact on squall-line propagation. It also does not affect the model’s overall effective resolution unlike applying constant multiplicative factors to the mixing coefficients.


Author(s):  
Jake P. Mulholland ◽  
John M. Peters ◽  
Hugh Morrison

AbstractThe influence of vertical wind shear on updraft entrainment in squall lines is not well understood. To address this knowledge gap, a suite of high-resolution idealized numerical model simulations of squall lines were run in various vertical wind shear (hereafter “shear”) environments to study the effects of shear on entrainment in deep convective updrafts. Low-level horizontal mass flux into the leading edge of the cold pool was strongest in the simulations with the strongest low-level shear. These simulations consequently displayed wider updrafts, less entrainment-driven dilution, and larger buoyancy than the simulations with comparatively weak low-level shear. An analysis of vertical accelerations along trajectories that passed through updrafts showed larger net accelerations from buoyancy in the simulations with stronger low-level shear, which demonstrates how less entrainment-driven dilution equated to stronger updrafts. The effects of upper-level shear on entrainment and updraft vertical velocities were generally less pronounced than the effects of low-level shear. We argue that in addition to the outflow boundary-shear interactions and their effect on updraft tilt established by previous authors, decreased entrainment-driven dilution is yet another beneficial effect of strong low-level shear on squall line updraft intensity.


2017 ◽  
Vol 17 (17) ◽  
pp. 10195-10221 ◽  
Author(s):  
Constantino Listowski ◽  
Tom Lachlan-Cope

Abstract. The first intercomparisons of cloud microphysics schemes implemented in the Weather Research and Forecasting (WRF) mesoscale atmospheric model (version 3.5.1) are performed on the Antarctic Peninsula using the polar version of WRF (Polar WRF) at 5 km resolution, along with comparisons to the British Antarctic Survey's aircraft measurements (presented in part 1 of this work; Lachlan-Cope et al., 2016). This study follows previous works suggesting the misrepresentation of the cloud thermodynamic phase in order to explain large radiative biases derived at the surface in Polar WRF continent-wide (at 15 km or coarser horizontal resolution) and in the Polar WRF-based operational forecast model Antarctic Mesoscale Prediction System (AMPS) over the Larsen C Ice Shelf at 5 km horizontal resolution. Five cloud microphysics schemes are investigated: the WRF single-moment five-class scheme (WSM5), the WRF double-moment six-class scheme (WDM6), the Morrison double-moment scheme, the Thompson scheme, and the Milbrandt–Yau double-moment seven-class scheme. WSM5 (used in AMPS) and WDM6 (an upgrade version of WSM5) lead to the largest biases in observed supercooled liquid phase and surface radiative biases. The schemes simulating clouds in closest agreement to the observations are the Morrison, Thompson, and Milbrandt schemes for their better average prediction of occurrences of clouds and cloud phase. Interestingly, those three schemes are also the ones allowing for significant reduction of the longwave surface radiative bias over the Larsen C Ice Shelf (eastern side of the peninsula). This is important for surface energy budget consideration with Polar WRF since the cloud radiative effect is more pronounced in the infrared over icy surfaces. Overall, the Morrison scheme compares better to the cloud observation and radiation measurements. The fact that WSM5 and WDM6 are single-moment parameterizations for the ice crystals is responsible for their lesser ability to model the supercooled liquid clouds compared to the other schemes. However, our investigation shows that all the schemes fail at simulating the supercooled liquid mass at some temperatures (altitudes) where observations show evidence of its persistence. An ice nuclei parameterization relying on both temperature and aerosol content like DeMott et al. (2010) (not currently used in WRF cloud schemes) is in best agreement with the observations, at temperatures and aerosol concentration characteristic of the Antarctic Peninsula where the primary ice production occurs (part 1), compared to parameterization only relying on the atmospheric temperature (used by the WRF cloud schemes). Overall, a realistic double-moment ice microphysics implementation is needed for the correct representation of the supercooled liquid phase in Antarctic clouds. Moreover, a more realistic ice-nucleating particle alone is not enough to improve the cloud modelling, and water vapour and temperature biases also need to be further investigated and reduced.


2006 ◽  
Vol 134 (7) ◽  
pp. 1919-1941 ◽  
Author(s):  
Jeffrey Frame ◽  
Paul Markowski

Abstract Numerical simulations of squall lines traversing sinusoidal mountain ridges are performed using the Advanced Regional Prediction System cloud-resolving model. Precipitation and updraft strength are enhanced through orographic ascent as a squall line approaches a ridge. The simulated squall line then weakens as it descends the ridge because some of the cold pool is blocked by the terrain, resulting in less lift along the gust front and weaker convective cells. The flow within the cold pool accelerates slightly and the depth of the cold air decreases owing to upstream blocking, transitioning the flow in the cold pool head from subcritical to supercritical, then back to subcritical at the bottom of the ridge. A hydraulic jump forms when the flow transitions the second time, enabling the development of a new convective line downwind of the mountain. These new updrafts grow and eventually replace the older updrafts that weakened during descent. This process results in the discrete propagation of a squall line just downstream of a ridge, resulting in the formation of rain shadows downstream from topographic features. Discrete propagation only occurs if a ridge is of sufficient height, however. This replacement process repeats itself if a squall line encounters multiple ridges. The risk of damaging winds from a squall line is greater on the lee side of ridges and on the top of high ridges. These terrain-forced intensity fluctuations increase with mountain height, because the higher terrain permits even less cold air to flow over it. A wider ridge results in a more gradual orographic enhancement and downslope-induced weakening.


2011 ◽  
Vol 139 (10) ◽  
pp. 3163-3183 ◽  
Author(s):  
Casey E. Letkewicz ◽  
Matthew D. Parker

Abstract The complex evolution of convective systems crossing (or attempting to cross) mountainous terrain represents a substantial forecasting challenge. This study examines the processes associated with environments of “crossing” squall lines (which were able to redevelop strong convection in the lee of a mountain barrier) and “noncrossing” squall lines (which were not able to redevelop strong convection downstream of the barrier). In particular, numerical simulations of mature convective systems crossing idealized terrain roughly approximating the Appalachian Mountains were used to test the first-order impact of variations in the vertical wind profile upon system maintenance. By itself, the wind profile showed no ability to uniquely discriminate between simulated crossing and noncrossing squall lines; each test revealed a similar pattern of orographic enhancement, suppression, and lee reinvigoration in which a hydraulic jump deepened the system’s cold pool and renewed the low-level lifting. Increasing the mean wind led to greater enhancement of vertical velocities on the windward side of the barrier and greater suppression on the lee side. Variations in the low-level shear influenced the temperature and depth of the outflow, which in turn altered the lifting along the system’s gust front. However, in all of the wind profile tests, convection redeveloped in the lee. Additional simulations explored more marginal environments in which idealized low-level cooling or drying stabilized the downstream environment. In most such tests, the systems weakened but the presence of CAPE aloft still enabled the systems to survive in the lee. However, the combination of a stronger mean wind with diminished CAPE and increased convective inhibition (CIN) was ultimately found to eliminate downstream redevelopment and produce a noncrossing mesoscale convective system (MCS). Within these experiments, the ability of a squall line to cross a barrier similar to the Appalachians is primarily tied to the characteristics of the downstream thermodynamic environment; however, as the lee thermodynamic environment becomes less favorable, the mean wind exerts a greater influence on system intensity and redevelopment.


Sign in / Sign up

Export Citation Format

Share Document