Effects of Horizontal and Vertical Grid Spacing on Mixing in Simulated Squall Lines and Implications for Convective Strength and Structure

2015 ◽  
Vol 143 (11) ◽  
pp. 4355-4375 ◽  
Author(s):  
Z. J. Lebo ◽  
H. Morrison

Abstract The sensitivity of an idealized squall line to horizontal and vertical grid spacing is investigated using a new approach. Simulations are first performed at a horizontal grid spacing of 1 km until the storm reaches its mature stage. The model output is then interpolated to smaller (and larger) grid spacings, and the model is restarted using the interpolated state plus small thermodynamic perturbations to spin up small-scale motions. This framework allows an investigation of the sensitivity of the storm to changes in without complications from differences in storm initiation and early evolution. The restarted simulations reach a quasi steady state within approximately 1 h. Results demonstrate that there are two -dependent regimes with the transition between regimes occurring for between 250 and 500 m. Some storm characteristics, such as the mean convective core area, change substantially for 250 m but show limited sensitivity as is decreased below 250 m, despite better resolving smaller-scale turbulent motions. This transition is found to be independent of the chosen . Mixing in the context of varying and is also investigated via passive tracers that are initialized 1 h after restarting the simulations (i.e., after the spin up of small-scale motions). The tracer field at the end of the simulations reveals that entrainment and detrainment are suppressed in the simulations with 500 m. For decreasing , entrainment and detrainment are substantially more important, limiting the flux of low-level tracer to the upper troposphere, which has important implications for modeling studies of convective transport from the boundary layer through the troposphere.

2015 ◽  
Vol 143 (8) ◽  
pp. 2998-3024 ◽  
Author(s):  
Corey K. Potvin ◽  
Montgomery L. Flora

Abstract The Warn-on-Forecast (WoF) program aims to deploy real-time, convection-allowing, ensemble data assimilation and prediction systems to improve short-term forecasts of tornadoes, flooding, lightning, damaging wind, and large hail. Until convection-resolving (horizontal grid spacing Δx < 100 m) systems become available, however, resolution errors will limit the accuracy of ensemble model output. Improved understanding of grid spacing dependence of simulated convection is therefore needed to properly calibrate and interpret ensemble output, and to optimize trade-offs between model resolution and other computationally constrained parameters like ensemble size and forecast lead time. Toward this end, the authors examine grid spacing sensitivities of simulated supercells over Δx of 333 m–4 km. Storm environment and physics parameterization are varied among the simulations. The results suggest that 4-km grid spacing is too coarse to reliably simulate supercells, occasionally leading to premature storm demise, whereas 3-km simulations more often capture operationally important features, including low-level rotation tracks. Further decreasing Δx to 1 km enables useful forecasts of rapid changes in low-level rotation intensity, though significant errors remain (e.g., in timing). Grid spacing dependencies vary substantially among the experiments, suggesting that accurate calibration of ensemble output requires better understanding of how storm characteristics, environment, and parameterization schemes modulate grid spacing sensitivity. Much of the sensitivity arises from poorly resolving small-scale processes that impact larger (well resolved) scales. Repeating some of the 333-m simulations with coarsened initial conditions reveals that supercell forecasts can substantially benefit from reduced grid spacing even when limited observational density precludes finescale initialization.


2014 ◽  
Vol 14 (2) ◽  
pp. 609-627 ◽  
Author(s):  
C. A. Klich ◽  
H. E. Fuelberg

Abstract. We use the Weather Research and Forecasting with Chemistry (WRF-Chem) online chemical transport model to simulate a middle latitude cyclone in East Asia at three different horizontal resolutions (45, 15, and 5 km grid spacing). The cyclone contains a typical warm conveyor belt (WCB) with an embedded squall line that passes through an area having large surface concentrations (> 400 ppbv) of carbon monoxide (CO). Model output from WRF-Chem is used to compare differences between the large-scale CO vertical transport by the WCB (the 45 km simulation) with the smaller-scale transport due to its convection (the 5 km simulation). Forward trajectories are calculated from WRF-Chem output using the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model. At 45 km grid spacing, the WCB exhibits gradual ascent, lofting surface CO to 6–7 km. Upon reaching the warm front, the WCB and associated CO ascend more rapidly and later turn eastward over the Pacific Ocean. Convective transport at 5 km resolution with explicitly resolved convection occurs much more rapidly, with surface CO lofted to altitudes greater than 10 km in 1 h or less. We also compute CO vertical mass fluxes over specified areas and times to compare differences in transport due to the different grid spacings. Upward CO flux exceeds 110 000 t h−1 in the domain with explicit convection when the squall line is at peak intensity, while fluxes from the two coarser resolutions are an order of magnitude smaller. Specific areas of interest within the 5 km domain are defined to compare the magnitude of convective transport to that within the entire 5 km region. Although convection encompasses only a small portion of the 5 km domain, it is responsible for ~40% of the upward CO transport. We also examine the vertical transport due to a short wave trough and its associated area of convection, not related to the cyclone, that lofts CO to the upper troposphere. Results indicate that fine-scale resolution with explicitly resolved convection is important when assessing the vertical transport of surface emissions in areas of deep convection.


2016 ◽  
Vol 97 (2) ◽  
pp. 237-243 ◽  
Author(s):  
Dale R. Durran ◽  
Jonathan A. Weyn

Abstract One important limitation on the accuracy of weather forecasts is imposed by unavoidable errors in the specification of the atmosphere’s initial state. Much theoretical concern has been focused on the limits to predictability imposed by small-scale errors, potentially even those on the scale of a butterfly. Very modest errors at much larger scales may nevertheless pose a more important practical limitation. We demonstrate the importance of large-scale uncertainty by analyzing ensembles of idealized squall-line simulations. Our results imply that minimizing initial errors on scales around 100 km is more likely to extend the accuracy of forecasts at lead times longer than 3–4 h than efforts to minimize initial errors on much smaller scales. These simulations also demonstrate that squall lines, triggered in a horizontally homogeneous environment with no initial background circulations, can generate a background mesoscale kinetic energy spectrum roughly similar to that observed in the atmosphere.


2018 ◽  
Vol 75 (4) ◽  
pp. 1243-1269 ◽  
Author(s):  
Kelly Lombardo ◽  
Tristan Kading

Abstract Inland squall lines respond to the stable marine atmospheric boundary layer (MABL) as they move toward a coastline and offshore. As a storm’s cold pool collides with the marine layer, characteristics of both determine the resulting convective forcing mechanism over the stable layer and storm characteristics. Idealized numerical experiments exploring a parameter space of MABL characteristics show that the postcollision forcing mechanism is determined by the buoyancy of the cold pool relative to the MABL. When the outflow is less buoyant, storms are forced by a cold pool within the marine environment. When the buoyancies are equivalent, a hybrid cold pool–internal gravity wave develops after the collision. The collision between a cold pool and less buoyant MABL initiates internal waves along the stable layer, regardless of MABL depth. These waves are inefficient at lifting air into the storm, and ascent from the trailing cold pool is needed to support deep convection. Storm intensity decreases with deeper and less buoyant MABLs, in part due to the reduction in elevated instability. Precipitation is enhanced just prior to the collision between a storm and the deepest marine layers. Storms modify their environment downstream, leading to the development of a moist adiabatic unstable layer and a lowering of the level of free convection (LFC) to below the top of the deepest marine layer. An MABL moving as a sea breeze into the storm-modified air successfully lifts parcels to the new LFC, generating convective towers ahead of the squall line. This mechanism may contribute to increased coastal flash flooding risks during observed events.


10.1175/811.1 ◽  
2004 ◽  
Vol 19 (6) ◽  
pp. 1061-1074 ◽  
Author(s):  
Patrick C. Burke ◽  
David M. Schultz

Abstract A search of radar mosaics and level-II Weather Surveillance Radar-1988 Doppler (WSR-88D) data revealed 51 cold-season (October–April) bow echoes that occurred in the contiguous United States from 1997–98 to 2000–01. Proximity soundings indicated mean 0–2.5-, 0–5-, and 5–10-km shear values of 14, 23, and 19 m s−1, respectively. Mean CAPE was 1366 J kg−1. Most bow echoes developed from squall lines, groups of cells, or squall lines overtaking cells that originated in the path of the squall line. Overall, cell mergers occurred just prior to the development of 34 (67%) of the 51 bow echoes, and embedded supercells were present in the mature stage of 22 (43%) bow echoes. Nine severe, long-lived bow echoes (LBEs) were identified, and seven of these had damage paths that met derecho criteria. LBEs developed in strongly forced, dynamic synoptic patterns with low to moderate instability. As in previous observational studies, proximity soundings suggested that LBEs are possible within much wider ranges of sampled CAPE and shear than idealized numerical modeling studies have indicated. Cold-season bow echoes formed overwhelmingly (47 of 51) in southwesterly 500-mb flow. Twenty (39%) bow echoes formed in a Gulf coast synoptic pattern that produced strong shear and moderate instability over the southeastern United States. Nineteen (37%) and seven (14%) bow echoes, respectively, formed in the plains and east synoptic patterns, which resemble classic severe weather outbreak patterns. Four (8%) bow echoes developed in a northwest flow synoptic pattern that produced strong shear and moderate instability over the southern plains.


2013 ◽  
Vol 13 (6) ◽  
pp. 14871-14925
Author(s):  
C. A. Klich ◽  
H. E. Fuelberg

Abstract. We use the Weather Research and Forecasting with Chemistry (WRF-Chem) online chemical transport model to simulate a middle latitude cyclone in East Asia at three different horizontal resolutions (45, 15, and 5 km grid spacing). The cyclone contains a typical warm conveyor belt (WCB) with an embedded squall line that passes through an area having large surface concentrations (>400 ppbv) of carbon monoxide (CO). Model output from WRF-Chem is used to compare differences between the large-scale CO vertical transport by the WCB (the 45 km simulation) with the smaller-scale transport due to its convection (the 5 km simulation). Forward trajectories are calculated from WRF-Chem output using HYSPLIT. At 45 km grid spacing, the WCB exhibits gradual ascent, lofting surface CO to 6–7 km. Upon reaching the warm front, the WCB and associated CO ascend more rapidly and later turn eastward over the Pacific Ocean. Convective transport at 5 km resolution with explicitly resolved convection occurs much more rapidly, with surface CO lofted to altitudes greater than 10 km in 1 h or less. We also compute CO vertical mass fluxes to compare differences in transport due to the different grid spacings. Upward CO flux exceeds 110 000 t h−1 in the domain with explicit convection when the squall line is at peak intensity, while fluxes from the two coarser resolutions are an order of magnitude smaller. Specific areas of interest within the 5 km domain are defined to compare the magnitude of convective transport to that within the entire 5 km region. Although convection encompasses only a small portion of the 5 km domain, it is responsible for ~40% of the upward CO transport. We also examine the vertical transport due to a short wave trough and its associated area of convection, not related to the cyclone, that lofts CO to the upper troposphere. Results indicate that fine-scale resolution with explicitly resolved convection is important when assessing the vertical transport of surface emissions in areas of deep convection.


2022 ◽  
Vol 9 ◽  
Author(s):  
Zhen Gao ◽  
Liguang Wu ◽  
Xingyang Zhou

It has been numerically demonstrated that the turbulence above the boundary is important to tropical cyclone intensification and rapid intensification, but the three-dimensional structures of the sub-grid-scale (SGS) eddy have not been revealed due to the lack of observational data. In this study, two numerical simulations of Super Typhoon Rammasun (2014) were conducted with the Advanced Weather Research and Forecast (WRF) model by incorporating the large-eddy simulation (LES) technique, in which the enhanced eyewall convection and the process of rapid intensification are captured. Consistent with previous observational studies, the strong turbulent kinetic energy (TKE) is found throughout the whole eyewall inside of the radius of maximum wind in both experiments. The simulations indicate that the strong TKE is associated with horizontal rolls with the horizontal extent of 2–4 km, which are aligned azimuthally in the intense eyewall convection. It is indicated that the three-dimensional structures of the SGS eddy can be simulated with the vertical grid spacing of ∼100 m when the horizontal grid spacing is 74 m. It is suggested that there is considerable turbulence associated with azimuthally-aligned horizontal rolls in the mid-level eyewall of tropical cyclone.


2019 ◽  
Vol 76 (10) ◽  
pp. 3267-3283 ◽  
Author(s):  
Cheng-Ku Yu ◽  
Che-Yu Lin ◽  
Jhang-Shuo Luo

Abstract This study used radar and surface observations to track a long-lasting outer tropical cyclone rainband (TCR) of Typhoon Jangmi (2008) over a considerable period of time (~10 h) from its formative to mature stage. Detailed analyses of these unique observations indicate that the TCR was initiated on the eastern side of the typhoon at a radial distance of ~190 km as it detached from the upwind segment of a stratiform rainband located close to the inner-core boundary. The outer rainband, as it propagated cyclonically outward, underwent a prominent convective transformation from generally stratiform precipitation during the earlier period to highly organized, convective precipitation during its mature stage. The transformation was accompanied by a clear trend of surface kinematics and thermodynamics toward squall-line-like features. The observed intensification of the rainband was not simply related to the spatial variation of the ambient CAPE or potential instability; instead, the dynamical interaction between the prerainband vertical shear and cold pools, with progression toward increasingly optimal conditions over time, provides a reasonable explanation for the temporal alternation of the precipitation intensity. The increasing intensity of cold pools was suggested to play an essential role in the convective transformation for the rainband. The propagation characteristics of the studied TCR were distinctly different from those of wave disturbances frequently documented within the cores of tropical cyclones; however, they were consistent with the theoretically predicted propagation of convectively generated cold pools. The convective transformation, as documented in the present case, is anticipated to be one of the fundamental processes determining the evolving and structural nature of outer TCRs.


1986 ◽  
Vol 1 (20) ◽  
pp. 193 ◽  
Author(s):  
Shiao-Kung Liu ◽  
Jan J. Leendertse

This paper presents the development of a three dimensional model of the Gulf of Alaska. The model extends between the Vancouver Island and the Aleutian Islands covering approximatedly 1.5 million square kilometers over the northern Pacific Ocean. Formulated on an ellipsoidal horizontal grid and variable vertical grid, the model is schematized over a 81 x 53 x 10 grid structure. The solution scheme is implicit over the vertical and is programmed using one-dimensional dynamic array for the efficient use of machine storage. The turbulence closure scheme for the non-homogeneous vertical shear is formulated so that the potential and kinetic energetics are monitored and transferred in a closed form. The hydrodynamic model is coupled to a two-dimensional stochastic weather model and an oil-spill trajectory/weathering model. The former also simulates stochastically the cyclogenetic/cyclolytic processes within the modeled area. The paper also compares the computed results with the available field data. Good agreements are found in tidal amplitude and phases as well as currents.


2016 ◽  
Vol 73 (11) ◽  
pp. 4289-4309 ◽  
Author(s):  
Tomoki Ohno ◽  
Masaki Satoh ◽  
Yohei Yamada

Abstract Based on the data of a 1-yr simulation by a global nonhydrostatic model with 7-km horizontal grid spacing, the relationships among warm-core structures, eyewall slopes, and the intensities of tropical cyclones (TCs) were investigated. The results showed that stronger TCs generally have warm-core maxima at higher levels as their intensities increase. It was also found that the height of a warm-core maximum ascends (descends) as the TC intensifies (decays). To clarify how the height and amplitude of warm-core maxima are related to TC intensity, the vortex structures of TCs were investigated. By gradually introducing simplifications of the thermal wind balance, it was established that warm-core structures can be reconstructed using only the tangential wind field within the inner-core region and the ambient temperature profile. A relationship between TC intensity and eyewall slope was investigated by introducing a parameter that characterizes the shape of eyewalls and can be evaluated from satellite measurements. The authors found that the eyewall slope becomes steeper (shallower) as the TC intensity increases (decreases). Based on a balanced model, the authors proposed a relationship between TC intensity and eyewall slope. The result of the proposed model is consistent with that of the analysis using the simulation data. Furthermore, for sufficiently strong TCs, the authors found that the height of the warm-core maximum increases as the slope becomes steeper, which is consistent with previous observational studies. These results suggest that eyewall slopes can be used to diagnose the intensities and structures of TCs.


Sign in / Sign up

Export Citation Format

Share Document