Search Results
You are looking at 31 - 40 of 67 items for
- Author or Editor: George H. Bryan x
- Refine by Access: All Content x
Abstract
A combination of idealized numerical simulations and analytical theory is used to investigate the spacing between convective orographic rainbands over the Coastal Range of western Oregon. The simulations, which are idealized from an observed banded precipitation event over the Coastal Range, indicate that the atmospheric response to conditionally unstable flow over the mountain ridge depends strongly on the subridge-scale topographic forcing on the windward side of the ridge. When this small-scale terrain contains only a single scale (λ) of terrain variability, the band spacing is identical to λ, but when a spectrum of terrain scales are simultaneously present, the band spacing ranges between 5 and 10 km, a value that is consistent with observations. Based on the simulations, an inviscid linear model is developed to provide a physical basis for understanding the scale selection of the rainbands. This analytical model, which captures the transition from lee waves upstream of the orographic cloud to moist convection within it, reveals that the spacing of orographic rainbands depends on both the projection of lee-wave energy onto the unstable cap cloud and the growth rate of unstable perturbations within the cloud. The linear model is used in tandem with numerical simulations to determine the sensitivity of the band spacing to a number of environmental and terrain-related parameters.
Abstract
A combination of idealized numerical simulations and analytical theory is used to investigate the spacing between convective orographic rainbands over the Coastal Range of western Oregon. The simulations, which are idealized from an observed banded precipitation event over the Coastal Range, indicate that the atmospheric response to conditionally unstable flow over the mountain ridge depends strongly on the subridge-scale topographic forcing on the windward side of the ridge. When this small-scale terrain contains only a single scale (λ) of terrain variability, the band spacing is identical to λ, but when a spectrum of terrain scales are simultaneously present, the band spacing ranges between 5 and 10 km, a value that is consistent with observations. Based on the simulations, an inviscid linear model is developed to provide a physical basis for understanding the scale selection of the rainbands. This analytical model, which captures the transition from lee waves upstream of the orographic cloud to moist convection within it, reveals that the spacing of orographic rainbands depends on both the projection of lee-wave energy onto the unstable cap cloud and the growth rate of unstable perturbations within the cloud. The linear model is used in tandem with numerical simulations to determine the sensitivity of the band spacing to a number of environmental and terrain-related parameters.
Abstract
In high-resolution numerical simulations (using horizontal grid spacing less than ∼1 km), the convective region of squall lines will sometimes overturn as quasi-horizontal convective rolls. The authors study one case in detail using output from a simulation with 125-m grid spacing. The rolls have an average spacing of 3 km and are aligned parallel to the vertical wind shear. Individual convective cells often have long-lived, undiluted cores that entrain primarily on the sides of the rolls (i.e., between the roll updraft and downdraft). The following set of conditions is proposed for obtaining roll overturning: the formation of a moist absolutely unstable layer (MAUL); vertical shear of the horizontal wind within the MAUL; an environment without large-amplitude perturbations; and quasi-horizontal flow over the squall line’s surface-based cold pool.
Further insight is gained through a series of more idealized simulations wherein a specified MAUL is perturbed by analytic potential temperature perturbations. These simulations confirm classical studies based on linear analysis because the smallest perturbations grow fastest (with the exception of the very smallest scales that are affected by diffusion). The results also confirm that, with shear, updrafts oriented across the shear vector are inhibited by the shear. An explanation for the ∼3-km roll spacing also emerges from these simulations. The argument focuses on the perturbations that exist in the cold pool underneath the MAUL; they induce pressure fields that extend upward into the overlying MAUL. The perturbations with large horizontal scale have pressure fields that extend farther vertically and with a greater amplitude, and thus are more effective at initiating motions in the overlying MAUL. The convective scale that ultimately emerges within the MAUL is somewhere between two scales, whereby comparatively large scales are perturbed more strongly by perturbations in the cold pool, but the comparatively small scales grow faster.
Abstract
In high-resolution numerical simulations (using horizontal grid spacing less than ∼1 km), the convective region of squall lines will sometimes overturn as quasi-horizontal convective rolls. The authors study one case in detail using output from a simulation with 125-m grid spacing. The rolls have an average spacing of 3 km and are aligned parallel to the vertical wind shear. Individual convective cells often have long-lived, undiluted cores that entrain primarily on the sides of the rolls (i.e., between the roll updraft and downdraft). The following set of conditions is proposed for obtaining roll overturning: the formation of a moist absolutely unstable layer (MAUL); vertical shear of the horizontal wind within the MAUL; an environment without large-amplitude perturbations; and quasi-horizontal flow over the squall line’s surface-based cold pool.
Further insight is gained through a series of more idealized simulations wherein a specified MAUL is perturbed by analytic potential temperature perturbations. These simulations confirm classical studies based on linear analysis because the smallest perturbations grow fastest (with the exception of the very smallest scales that are affected by diffusion). The results also confirm that, with shear, updrafts oriented across the shear vector are inhibited by the shear. An explanation for the ∼3-km roll spacing also emerges from these simulations. The argument focuses on the perturbations that exist in the cold pool underneath the MAUL; they induce pressure fields that extend upward into the overlying MAUL. The perturbations with large horizontal scale have pressure fields that extend farther vertically and with a greater amplitude, and thus are more effective at initiating motions in the overlying MAUL. The convective scale that ultimately emerges within the MAUL is somewhere between two scales, whereby comparatively large scales are perturbed more strongly by perturbations in the cold pool, but the comparatively small scales grow faster.
Abstract
Numerical models of supercell thunderstorms produce near-ground rotation about a vertical axis (i.e., vertical vorticity) after the development of rain-cooled outflows and downdrafts. The physical processes involved in the production of near-ground vertical vorticity in simulated supercells have been a subject of discussion in the literature for over 30 years. One cause for this lengthy discussion is the difficulty in applying the principles of inviscid vorticity dynamics in a continuous fluid to the viscous evolution of discrete Eulerian simulations. The present paper reports on a Lagrangian analysis of near-ground vorticity from an idealized-supercell simulation with enhanced vertical resolution near the lower surface. The parcel that enters the low-level maximum of vertical vorticity has a history of descent during which its horizontal vorticity is considerably enhanced. In its final approach to this region, the parcel’s enhanced horizontal vorticity is tilted to produce vertical vorticity, which is then amplified through vertical stretching as the parcel rises. A simplified theoretical model is developed that exhibits these same features. The principal conclusion is that vertical vorticity at the parcel’s nadir (its lowest point), although helpful, does not need to be positive for rapid near-surface amplification of vertical vorticity.
Abstract
Numerical models of supercell thunderstorms produce near-ground rotation about a vertical axis (i.e., vertical vorticity) after the development of rain-cooled outflows and downdrafts. The physical processes involved in the production of near-ground vertical vorticity in simulated supercells have been a subject of discussion in the literature for over 30 years. One cause for this lengthy discussion is the difficulty in applying the principles of inviscid vorticity dynamics in a continuous fluid to the viscous evolution of discrete Eulerian simulations. The present paper reports on a Lagrangian analysis of near-ground vorticity from an idealized-supercell simulation with enhanced vertical resolution near the lower surface. The parcel that enters the low-level maximum of vertical vorticity has a history of descent during which its horizontal vorticity is considerably enhanced. In its final approach to this region, the parcel’s enhanced horizontal vorticity is tilted to produce vertical vorticity, which is then amplified through vertical stretching as the parcel rises. A simplified theoretical model is developed that exhibits these same features. The principal conclusion is that vertical vorticity at the parcel’s nadir (its lowest point), although helpful, does not need to be positive for rapid near-surface amplification of vertical vorticity.
Abstract
Output from idealized simulations of a squall line are compared with in situ aircraft-based observations from the Deep Convective Clouds and Chemistry campaign. Relative humidity distributions around convection are compared between 1-Hz aircraft observations (≈250-m horizontal scale) and simulations using a double-moment bulk microphysics scheme at three horizontal grid spacings: Δx = 0.25, 1, and 4 km. The comparisons focus on the horizontal extent of ice supersaturated regions (ISSRs), the maximum and average relative humidity with respect to ice (RHi) in ISSRs, and the ice microphysical properties during cirrus cloud evolution, with simulations at 0.25 and 1 km providing better results than the 4-km simulation. Within the ISSRs, all the simulations represent the dominant contributions of water vapor horizontal heterogeneities to ISSR formation on average, but with larger variabilities in such contributions than the observations. The best results are produced by a Δx = 0.25-km simulation with the RHi threshold for initiating ice nucleation increased to 130%, which improves almost all the ISSR characteristics and allows for larger magnitude and frequency of ice supersaturation (ISS) > 8%. This simulation also allows more occurrences of clear-sky ISSRs and a higher spatial fraction of ISS for in-cloud conditions, which are consistent with the observations. These improvements are not reproduced by modifying other ice microphysical processes, such as a factor-of-2 reduction in the ice nuclei concentration; a factor-of-10 reduction in the vapor deposition rate; turning off heterogeneous contact and immersion freezing; or turning off homogeneous freezing of liquid water.
Abstract
Output from idealized simulations of a squall line are compared with in situ aircraft-based observations from the Deep Convective Clouds and Chemistry campaign. Relative humidity distributions around convection are compared between 1-Hz aircraft observations (≈250-m horizontal scale) and simulations using a double-moment bulk microphysics scheme at three horizontal grid spacings: Δx = 0.25, 1, and 4 km. The comparisons focus on the horizontal extent of ice supersaturated regions (ISSRs), the maximum and average relative humidity with respect to ice (RHi) in ISSRs, and the ice microphysical properties during cirrus cloud evolution, with simulations at 0.25 and 1 km providing better results than the 4-km simulation. Within the ISSRs, all the simulations represent the dominant contributions of water vapor horizontal heterogeneities to ISSR formation on average, but with larger variabilities in such contributions than the observations. The best results are produced by a Δx = 0.25-km simulation with the RHi threshold for initiating ice nucleation increased to 130%, which improves almost all the ISSR characteristics and allows for larger magnitude and frequency of ice supersaturation (ISS) > 8%. This simulation also allows more occurrences of clear-sky ISSRs and a higher spatial fraction of ISS for in-cloud conditions, which are consistent with the observations. These improvements are not reproduced by modifying other ice microphysical processes, such as a factor-of-2 reduction in the ice nuclei concentration; a factor-of-10 reduction in the vapor deposition rate; turning off heterogeneous contact and immersion freezing; or turning off homogeneous freezing of liquid water.
Abstract
In a previous study, the authors showed that the intensification process of a numerically simulated axisymmetric tropical cyclone (TC) can be divided into two periods denoted by “phase I” and “phase II.” The intensification process in phase II can be qualitatively described by Emanuel’s intensification theory in which the angular momentum (M) and saturated entropy (s*) surfaces are congruent in the TC interior. During phase I, however, the M and s* surfaces evolve from nearly orthogonal to almost congruent, and thus, the intensifying simulated TC has a different physical character as compared to that found in phase II. The present work uses a numerical simulation to investigate the evolution of an axisymmetric TC during phase I. The present results show that sporadic, deep convective annular rings play an important role in the simulated axisymmetric TC evolution in phase I. The convergence in low-level radial (Ekman) inflow in the boundary layer of the TC vortex, together with the increase of near-surface s* produced by sea surface fluxes, leads to episodes of convective rings around the TC center. These convective rings transport larger values of s* and M from the lower troposphere upward to the tropopause; the locally large values of M associated with the convective rings cause a radially outward bias in the upper-level radial velocity and an inward bias in the low-level radial velocity. Through a repetition of this process, the pattern (i.e., phase II) gradually emerges. The role of internal gravity waves related to the episodes of convection and the TC intensification process during phase I is also discussed.
Abstract
In a previous study, the authors showed that the intensification process of a numerically simulated axisymmetric tropical cyclone (TC) can be divided into two periods denoted by “phase I” and “phase II.” The intensification process in phase II can be qualitatively described by Emanuel’s intensification theory in which the angular momentum (M) and saturated entropy (s*) surfaces are congruent in the TC interior. During phase I, however, the M and s* surfaces evolve from nearly orthogonal to almost congruent, and thus, the intensifying simulated TC has a different physical character as compared to that found in phase II. The present work uses a numerical simulation to investigate the evolution of an axisymmetric TC during phase I. The present results show that sporadic, deep convective annular rings play an important role in the simulated axisymmetric TC evolution in phase I. The convergence in low-level radial (Ekman) inflow in the boundary layer of the TC vortex, together with the increase of near-surface s* produced by sea surface fluxes, leads to episodes of convective rings around the TC center. These convective rings transport larger values of s* and M from the lower troposphere upward to the tropopause; the locally large values of M associated with the convective rings cause a radially outward bias in the upper-level radial velocity and an inward bias in the low-level radial velocity. Through a repetition of this process, the pattern (i.e., phase II) gradually emerges. The role of internal gravity waves related to the episodes of convection and the TC intensification process during phase I is also discussed.
Abstract
Unsteadiness and horizontal heterogeneities frequently characterize atmospheric motions, especially within convective storms, which are frequently studied using large-eddy simulations (LES). The models of near-surface turbulence employed by atmospheric LES, however, predominantly assume statistically steady and horizontally homogeneous conditions (known as the equilibrium approach). The primary objective of this work is to investigate the potential consequences of such unrealistic assumptions in simulations of tornadoes. Cloud Model 1 (CM1) LES runs are performed using three approaches to model near-surface turbulence: the “semi-slip” boundary condition (which is the most commonly used equilibrium approach), a recently proposed nonequilibrium approach that accounts for some of the effects of turbulence memory, and a nonequilibrium approach based on thin boundary layer equations (TBLE) originally proposed by the engineering community for smooth-wall boundary layer applications. To be adopted for atmospheric applications, the TBLE approach is modified to account for the surface roughness. The implementation of TBLE into CM1 is evaluated using LES results of an idealized, neutral atmospheric boundary layer. LES runs are then performed for an idealized tornado characterized by rapid evolution, strongly curved air parcel trajectories, and substantial horizontal heterogeneities. The semi-slip boundary condition, by design, always yields a surface shear stress opposite the horizontal wind at the lowest LES grid level. The nonequilibrium approaches of modeling near-surface turbulence allow for a range of surface-shear-stress directions and enhance the resolved turbulence and wind gusts. The TBLE approach even occasionally permits kinetic energy backscatter from unresolved to resolved scales.
Significance Statement
The traditional approach of modeling the near-surface turbulence is not suitable for a tornado characterized by rapid evolution, strongly curved air parcel trajectories, and substantial horizontal heterogeneities. To understand the influence of statistically unsteady and horizontally heterogeneous near-surface conditions on tornadoes, this work adopts a fairly sophisticated approach from the engineering community and implements it into a widely used atmospheric model with necessary modifications. Compared to the traditional approach, the newly implemented approach produces more turbulent near-surface winds, more flexible surface-drag directions, and stronger wind gusts. These findings suggest a simulated tornado is very sensitive to the modeling approach of near-surface turbulence.
Abstract
Unsteadiness and horizontal heterogeneities frequently characterize atmospheric motions, especially within convective storms, which are frequently studied using large-eddy simulations (LES). The models of near-surface turbulence employed by atmospheric LES, however, predominantly assume statistically steady and horizontally homogeneous conditions (known as the equilibrium approach). The primary objective of this work is to investigate the potential consequences of such unrealistic assumptions in simulations of tornadoes. Cloud Model 1 (CM1) LES runs are performed using three approaches to model near-surface turbulence: the “semi-slip” boundary condition (which is the most commonly used equilibrium approach), a recently proposed nonequilibrium approach that accounts for some of the effects of turbulence memory, and a nonequilibrium approach based on thin boundary layer equations (TBLE) originally proposed by the engineering community for smooth-wall boundary layer applications. To be adopted for atmospheric applications, the TBLE approach is modified to account for the surface roughness. The implementation of TBLE into CM1 is evaluated using LES results of an idealized, neutral atmospheric boundary layer. LES runs are then performed for an idealized tornado characterized by rapid evolution, strongly curved air parcel trajectories, and substantial horizontal heterogeneities. The semi-slip boundary condition, by design, always yields a surface shear stress opposite the horizontal wind at the lowest LES grid level. The nonequilibrium approaches of modeling near-surface turbulence allow for a range of surface-shear-stress directions and enhance the resolved turbulence and wind gusts. The TBLE approach even occasionally permits kinetic energy backscatter from unresolved to resolved scales.
Significance Statement
The traditional approach of modeling the near-surface turbulence is not suitable for a tornado characterized by rapid evolution, strongly curved air parcel trajectories, and substantial horizontal heterogeneities. To understand the influence of statistically unsteady and horizontally heterogeneous near-surface conditions on tornadoes, this work adopts a fairly sophisticated approach from the engineering community and implements it into a widely used atmospheric model with necessary modifications. Compared to the traditional approach, the newly implemented approach produces more turbulent near-surface winds, more flexible surface-drag directions, and stronger wind gusts. These findings suggest a simulated tornado is very sensitive to the modeling approach of near-surface turbulence.
Abstract
Understanding momentum exchange at the air–sea interface is important for accurate hurricane predictions and understanding fundamental storm dynamics. One method for estimating air–sea momentum transfer in high winds is the flux-profile method, which infers surface momentum fluxes and the corresponding drag coefficient from mean velocity profiles obtained from either dropsondes or meteorological towers, under the assumption that the boundary layer wind profile at low altitudes exhibits a logarithmic profile with height. In this study, we use dropsonde data from reconnaissance aircraft, as well as “virtual sondes” from a turbulence-resolving simulation of an intense tropical cyclone, to critically analyze the diagnosis of drag coefficient C D at hurricane-force wind speeds. In particular, the “rolloff” of the drag coefficient, where C D decreases at 10-m wind speeds > 35 m s−1, is called into question based on uncertainty due to relatively low sample size and a lack of robustness of the flux-profile method at high winds. In addition, multiple factors appear to favor an underestimate of C D at hurricane-force winds relative to their true values, including uncertainty in the height of recorded dropsonde data, in violation of Monin–Obukhov similarity theory near the eyewall, and the short vertical extent of the logarithmic layer. Due to these and other related sources of uncertainty, it is likely that a quantitative limit has been reached in inferring the specific values of u * and C D using the flux-profile method, while at the same time the potential for underestimation may cast doubt on the C D –U 10 relationship inferred from this method at high winds.
Abstract
Understanding momentum exchange at the air–sea interface is important for accurate hurricane predictions and understanding fundamental storm dynamics. One method for estimating air–sea momentum transfer in high winds is the flux-profile method, which infers surface momentum fluxes and the corresponding drag coefficient from mean velocity profiles obtained from either dropsondes or meteorological towers, under the assumption that the boundary layer wind profile at low altitudes exhibits a logarithmic profile with height. In this study, we use dropsonde data from reconnaissance aircraft, as well as “virtual sondes” from a turbulence-resolving simulation of an intense tropical cyclone, to critically analyze the diagnosis of drag coefficient C D at hurricane-force wind speeds. In particular, the “rolloff” of the drag coefficient, where C D decreases at 10-m wind speeds > 35 m s−1, is called into question based on uncertainty due to relatively low sample size and a lack of robustness of the flux-profile method at high winds. In addition, multiple factors appear to favor an underestimate of C D at hurricane-force winds relative to their true values, including uncertainty in the height of recorded dropsonde data, in violation of Monin–Obukhov similarity theory near the eyewall, and the short vertical extent of the logarithmic layer. Due to these and other related sources of uncertainty, it is likely that a quantitative limit has been reached in inferring the specific values of u * and C D using the flux-profile method, while at the same time the potential for underestimation may cast doubt on the C D –U 10 relationship inferred from this method at high winds.
Abstract
This study examines how midlevel dry air and vertical wind shear (VWS) can modulate tropical cyclone (TC) development via downdraft ventilation. A suite of experiments was conducted with different combinations of initial midlevel moisture and VWS. A strong, positive, linear relationship exists between the low-level vertical mass flux in the inner core and TC intensity. The linear increase in vertical mass flux with intensity is not due to an increased strength of upward motions but, instead, is due to an increased areal extent of strong upward motions (w > 0.5 m s−1). This relationship suggests physical processes that could influence the vertical mass flux, such as downdraft ventilation, influence the intensity of a TC. The azimuthal asymmetry and strength of downdraft ventilation is associated with the vertical tilt of the vortex: downdraft ventilation is located cyclonically downstream from the vertical tilt direction and its strength is associated with the magnitude of the vertical tilt. Importantly, equivalent potential temperature of parcels associated with downdraft ventilation trajectories quickly recovers via surface fluxes in the subcloud layer, but the areal extent of strong upward motions is reduced. Altogether, the modulating effects of downdraft ventilation on TC development are the downward transport of low–equivalent potential temperature, negative-buoyancy air left of shear and into the upshear semicircle, as well as low-level radial outflow upshear, which aid in reducing the areal extent of strong upward motions, thereby reducing the vertical mass flux in the inner core, and stunting TC development.
Abstract
This study examines how midlevel dry air and vertical wind shear (VWS) can modulate tropical cyclone (TC) development via downdraft ventilation. A suite of experiments was conducted with different combinations of initial midlevel moisture and VWS. A strong, positive, linear relationship exists between the low-level vertical mass flux in the inner core and TC intensity. The linear increase in vertical mass flux with intensity is not due to an increased strength of upward motions but, instead, is due to an increased areal extent of strong upward motions (w > 0.5 m s−1). This relationship suggests physical processes that could influence the vertical mass flux, such as downdraft ventilation, influence the intensity of a TC. The azimuthal asymmetry and strength of downdraft ventilation is associated with the vertical tilt of the vortex: downdraft ventilation is located cyclonically downstream from the vertical tilt direction and its strength is associated with the magnitude of the vertical tilt. Importantly, equivalent potential temperature of parcels associated with downdraft ventilation trajectories quickly recovers via surface fluxes in the subcloud layer, but the areal extent of strong upward motions is reduced. Altogether, the modulating effects of downdraft ventilation on TC development are the downward transport of low–equivalent potential temperature, negative-buoyancy air left of shear and into the upshear semicircle, as well as low-level radial outflow upshear, which aid in reducing the areal extent of strong upward motions, thereby reducing the vertical mass flux in the inner core, and stunting TC development.
Abstract
This study demonstrates how midlevel dry air and vertical wind shear (VWS) can modulate tropical cyclone (TC) development via radial ventilation. A suite of experiments was conducted with different combinations of initial midlevel moisture and VWS environments. Two radial ventilation structures are documented. The first structure is positioned in a similar region as rainband activity and downdraft ventilation (documented in Part I) between heights of 0 and 3 km. Parcels associated with this first structure transport low–equivalent potential temperature air inward and downward left of shear and upshear to suppress convection. The second structure is associated with the vertical tilt of the vortex and storm-relative flow between heights of 5 and 9 km. Parcels associated with this second structure transport low–relative humidity air inward upshear and right of shear to suppress convection. Altogether, the modulating effects of radial ventilation on TC development are the inward transport of low–equivalent potential temperature air, as well as low-level radial outflow upshear, which aid in reducing the areal extent of strong upward motions, thereby reducing the vertical mass flux in the inner core, and stunting TC development.
Abstract
This study demonstrates how midlevel dry air and vertical wind shear (VWS) can modulate tropical cyclone (TC) development via radial ventilation. A suite of experiments was conducted with different combinations of initial midlevel moisture and VWS environments. Two radial ventilation structures are documented. The first structure is positioned in a similar region as rainband activity and downdraft ventilation (documented in Part I) between heights of 0 and 3 km. Parcels associated with this first structure transport low–equivalent potential temperature air inward and downward left of shear and upshear to suppress convection. The second structure is associated with the vertical tilt of the vortex and storm-relative flow between heights of 5 and 9 km. Parcels associated with this second structure transport low–relative humidity air inward upshear and right of shear to suppress convection. Altogether, the modulating effects of radial ventilation on TC development are the inward transport of low–equivalent potential temperature air, as well as low-level radial outflow upshear, which aid in reducing the areal extent of strong upward motions, thereby reducing the vertical mass flux in the inner core, and stunting TC development.
Abstract
Nearly all previous numerical simulations of supercell thunderstorms have neglected surface fluxes of heat, moisture, and momentum. This choice precludes horizontal inhomogeneities associated with dry boundary layer convection in the near-storm environment. As part of a broader study on how mature supercell thunderstorms are affected by a convective boundary layer (CBL) with quasi-two-dimensional features (i.e., boundary layer rolls), this paper documents the methods used to develop a realistic CBL in an idealized environment supportive of supercells. The evolution and characteristics of the modeled CBL, including the horizontal variability of thermodynamic and kinematic quantities known to affect supercell evolution, are presented. The simulated rolls result in periodic bands of perturbations in temperature, moisture, convective available potential energy (CAPE), vertical wind shear, and storm-relative helicity (SRH). Vertical vorticity is shown to arise within the boundary layer through the tilting of ambient horizontal vorticity associated with the background shear by vertical velocity perturbations in the turbulent CBL. Sensitivity tests suggest that 200-m horizontal grid spacing is adequate to represent rolls using a large-eddy simulation (LES) approach.
Abstract
Nearly all previous numerical simulations of supercell thunderstorms have neglected surface fluxes of heat, moisture, and momentum. This choice precludes horizontal inhomogeneities associated with dry boundary layer convection in the near-storm environment. As part of a broader study on how mature supercell thunderstorms are affected by a convective boundary layer (CBL) with quasi-two-dimensional features (i.e., boundary layer rolls), this paper documents the methods used to develop a realistic CBL in an idealized environment supportive of supercells. The evolution and characteristics of the modeled CBL, including the horizontal variability of thermodynamic and kinematic quantities known to affect supercell evolution, are presented. The simulated rolls result in periodic bands of perturbations in temperature, moisture, convective available potential energy (CAPE), vertical wind shear, and storm-relative helicity (SRH). Vertical vorticity is shown to arise within the boundary layer through the tilting of ambient horizontal vorticity associated with the background shear by vertical velocity perturbations in the turbulent CBL. Sensitivity tests suggest that 200-m horizontal grid spacing is adequate to represent rolls using a large-eddy simulation (LES) approach.
