Hemler, Richard S., F B Lipps, and Bruce B Ross, 1991: A simulation of a squall line using a nonhydrostatic cloud model with a 5-km horizontal grid. Monthly Weather Review, 119(12), 3012-3033. Abstract PDF
A three-dimensional nonhydrostatic cloud model is used to simulate the squall line observed in central Texas on 11 April 1979. The cloud model covers an area 400 x 400 km2 with a 5-km horizontal resolution and is supplied initial and boundary conditions by a larger hydrostatic mesoscale model.
The model produces a back-building squall line ahead of the surface cold front, as would be expected based on an analysis of the pre-squall-line environment. A well-defined gust front and cold pool develop with the squall line. At the end of the 5-h simulation, deep convection is found along a line nearly 400 km long. The simulated squall line compares favorably both with observations and with a higher-resolution model simulation in an environment of similar shear, suggesting that the 5-km horizontal resolution is adequately representing the significant features of the squall line.
The major shortcoming of this study is the failure of the cloud model to produce the observed squall line at the proper time. Without the observed small-scale forcing, which was unresolved in the Severe Environmental Storms and Mesoscale Experiment (SESAME) dataset, the model is unable to generate the squall line until a larger-scale convergence area evolves, some 2-3 h after the appearance of the observed squall line.
Ross, Bruce B., 1987: The role of low-level convergence and latent heating in a simulation of observed squall line formation. Monthly Weather Review, 115(10), 2298-2321. Abstract PDF
The effects of mesoscale forcing and diabatic heating on the development of convective systems have been investigated using a simplified numerical model to simulate the squall line and the convective system preceding it that occurred over Texas and Oklahoma on 10-11 April 1979. A simulation run without including latent heat showed both systems to be initiated and maintained by convergence produced by larger-scale forcing. The first cloud system formed downwind of the convergence zone that was produced by the confluence of airstreams along a dryline. A cloud front approaching from the west then merged with this dryline, destroying its horizontal gradients through diffusive effects and replacing it with a frontal convergence line that was aligned with the low-level flow. This new configuration was then favorable for the formation of the squall line that developed in the simulation.
When latent heat was included, the continuous cloud in the first convective system broke down into isolated cells which moved downstream from the convergence zone. In the non-latent heat case, the primary mechanism for providing moisture to this cloud was vertical diffusion from the moist surface layer. When latent heat was added, vertical advection within cell updrafts provided a more efficient means to supply moisture to the convective system.
In the simulated squall line, latent heat release produced a deeper cloud system while intensifying and maintaining the low-level convergence. However, unlike the earlier system, the squall line did not break into convective cells when latent heat was included in the simulation.
Orlanski, Isidoro, and Bruce B Ross, 1986: Low-level updrafts in stable layers forced by convection. Journal of the Atmospheric Sciences, 43(10), 997-1005. Abstract PDF
An investigation is made of the stability of a convectively unstable atmosphere in the presence of a stably stratified layer beneath, which is moving with a constant velocity relative to the upper air. This work is an extension of the linear model presented as part of the recent study of Orlanski and Ross in which they sought to explain the structure of their simulated squall line. A stability analysis shows that two modes are possible: 1) the gravitational or convective mode due to the unstable stratification in the upper layer which modifies the stable region below and 2) the classical Kelvin-Helmholtz mode due to shear across the interface. The Kelvin-Helmholtz mode is of limited physical interest in this case. On the other hand, the gravitational mode produces an updraft structure similar to updrafts in the stable lower layer of a convective system. Analysis of the horizontal displacement of the surface convergence for this mode relative to the convergence in the convective zone shows this displacement to depend primarily on the wind, stratification, and depth of the stable lower layer. The resulting relationship provides a method for determining whether a dual or single updraft will occur in a convective system.
Ross, Bruce B., 1986: An overview of numerical weather prediction In Mesoscale Meteorology and Forecasting, Boston, MA, American Meteorological Society, 720-751. PDF
Orlanski, Isidoro, and Bruce B Ross, 1984: The evolution of an observed cold front. Part II: mesoscale dynamics. Journal of the Atmospheric Sciences, 41(10), 1669-1703. Abstract PDF
A detailed analysis is made of a three-dimensional numerical simulation of the evolution of an observed moist frontal system over a 48 h period. The simulated front undergoes an initial period of frontogenetic growth, characterized by an alignment of vertical vorticity and horizontal convergence near the surface. The front then evolves to a mature, quasi-steady state as the line of maximum convergence moves ahead of the maximum vorticity. This phase shift is shown to be the result of a negative feedback mechanism which inhibits further vorticity growth while reducing the amount of viscous damping required to achieve a steady state. The influence of viscosity and surface drag upon this mechanism is also assessed.
When moisture is included in the numerical solution, the squall line which develops along the front exhibits a dual updraft structure with low-level convergence near the nose of the front and midlevel convergence located 100 km to the rear at a height of 3 km. This configuration is very similar to that found by Ogura and Liou in their analysis of an Oklahoma squall line not associated with a cold front.
Analysis of the equations of motion within the convective zone of the mature squall line shows the diabatic heating to be closely balanced by adiabatic cooling due to vertical temperature advection. As a result, the only net warming within this region occurs as adiabatic warming in the clear air outside of the cloud zone.
A linear, two-layer, dry model containing stable lower and unstable upper layers is shown to reproduce the dual updraft structure for certain low-level wind intensities without requiring microphysics. Also, for all wind conditions, this simple model produces strong convergence at the interface between the two layers. This suggests that the occurrence of a convergence maximum at the level of free convection should be a common feature of convectively unstable cloud systems.
Ross, Bruce B., and Isidoro Orlanski, 1982: The evolution of an observed cold front. Part I: Numerical simulation. Journal of the Atmospheric Sciences, 39(3), 296-327. Abstract PDF
The 48 h evolution of an observed cold front is simulated by a three-dimensional mesoscale-numerical model for a typical springtime synoptic situation over the southeastern United States. The model used in this study employs anelastic equations of motion on a limited-area domain with locally determined inflow/outflow side boundaries.
Both the observed and simulated characteristics of the weather system indicate a mature front which intensifies and then weakens over the 48 h period. Moist convection occurs in the form of intermittent squall lines in the observed case; in the numerical simulation, convection develops above and somewhat ahead of the surface front after 24 h as an ensemble of convective cells.
An investigation is made of the mesoscale and subsynoptic-scale features of this solution to determine their sensitivity to the inclusion of moisture and to the magnitude of the eddy viscosity used in the numerical simulation. The primary effect of increased eddy viscosity is to reduce somewhat the propagation speed of the front. The major changes due to moisture inclusion occur when convection develops along the cold front; these convective effects, which are apparent in the subsynoptic as well as the mesoscale features of the solution, include increased low-level convergence, reduced surface pressure due to diabatic heating, and the deflection of winds due to upper-level divergence. In addition, small temperature changes occur in the middle troposphere between the jet stream and the surface front when either viscosity or moisture is varied; these disturbances are a clear manifestation of the effect which changes in the cross-stream circulation intensity have upon the frontal system.
A fundamental feature of the mesoscale structure of the front in all cases is the tendency of the line of maximum horizontal convergence at the surface to move ahead of the line of maximum vertical vorticity. This phase shift appears to be related to the propagation characteristics of the frontal system. Also, the mesoscale moist convection develops a cellular structure throughout the convective zone in the low-viscosity solution; the use of higher viscosity tends to suppress these cells, particularly near the surface.
Orlanski, Isidoro, and Bruce B Ross, 1979: 3-D simulation of a frontal squall line: Sensitivity studies In Eleventh Conference on Severe Local Storms, Boston, MA, American Meteorological Society, 557-558.
Ross, Bruce B., and Isidoro Orlanski, 1978: The circulation associated with a cold front. Part II: Moist case. Journal of the Atmospheric Sciences, 35(3), 445-465. Abstract PDF
The effect of moisture upon the dynamics of mature idealized cold front systems is investigated using a two-dimensional numerical model. Lifting produced by the initial cross-stream frontal circulation studied by Orlanski and Ross (1977) is shown to saturate the warm moist air above the nose of the front when initial humidity levels are sufficiently high. If the atmosphere is convectively unstable, this saturated air will develop into deep convection with the convection-induced circulation overwhelming the initial frontal circulation. The initial development of convection is also shown to produce a gravity wave exhibiting similar scales to those of the convective zone. This wave propagates into the warm air at a much faster speed than the moving front-cloud system. Comparisons are made of the intensity of convection for different initial humidity and temperature conditons and when a low-level capping inversion is present. Also a comparison is made of cloud development caused by a combination of frontal lifting and surface heating when temperature inversions of different intensities are present. The stonger inversion is shown to suppress convection produced by surface heating alone with the combined effect of frontal lifting and surface heating required to release the convective instability.
Orlanski, Isidoro, and Bruce B Ross, 1977: The circulation associated with a cold front. Part I: Dry case. Journal of the Atmospheric Sciences, 34(10), 1619-1633. Abstract PDF
The transient behavior of an idealized dry frontal system is investigated using a two-dimensional numerical model. The development of a cross-stream circulation within stationary and moving cold fronts is determined for various frontal and synoptic conditions. In the stationary front, a circulation is generated by symmetric baroclinic instability, but nonlinear effects restrict this circulation to remain very weak. In the moving cold front, the vertical shear of the synoptic wind which advects the front produces an ageostrophic residue as a result of the differential advection of the vertical shear of the frontal jet and the horizontal temperature gradient across the front. This residue, which depends upon the vertical synoptic shear and the thermal wind structure of the frontal system, will generate a cross-stream circulation which maintains the cold front in a quasi-steady state. The resulting motion field is described well by the streamfunction balance equation. The lifting produced by the cross-stream circulation in the moving cold front system may be sufficient to trigger deep convection under favorable conditions in the moisture and synoptic wind fields.
A two-dimensional mesoscale atmospheric model is presented and used to study unsteady dynamic processes occurring in the planetary boundary layer (PBL) driven by diurnal heating at the ground. The model reproduces turbulent fluxes of heat and momentum both by explicitly modeling resolvable eddies and by employing a single parameterization at all levels of the model to represent vertical fluxes caused by subgrid-scale eddies. The unsteady behavior of horizontally-averaged profiles of temperature and velocity respond quite realistically to the diurnally-varying heat flux at the ground, particularly with regard to the time variation of lapse rates and the occurrence times of maximum and minimum temperatures at various levels in the lower boundary layer. The spatial variation of predicted atmospheric quantities shows a great deal of resolved eddy activity during the day with a significant remnant persisting through the night at higher levels of the PBL. These eddies account for the predominant means of vertical heat and momentum transfer away from the surfaces with the model realistically reproducing the unsteady behavior of heat fluxes in the PBL. Temporal variation of vertical heat and momentum profiles shows boundary layer activity to be confined to a few hundred meters at night while extending up to a kilometer during the day. A weak heat flux source at the ground with an amplitude of 10% of the maximum daytime heating produced a nocturnal heat island some 60 m high with a maximum city-country temperature contrast of ~1C.
Orlanski, Isidoro, and Bruce B Ross, 1973: Numerical simulation of the generation and breaking of internal gravity waves. Journal of Geophysical Research, 78(36), 8808-8826. Abstract
A numerical model is presented that permits the simulation of stratified fluid phenomena in which gravitationally unstable regions are present. The influence of subgrid scale turbulence generation due to convective instability is parameterized by relating eddy viscosity coefficients to the local Rayleigh number in unstable regions. The model is used to study three different laboratory scale flow problems involving gravity wave generation, wave breaking, and penetrative convection. The numerical solutions show good agreement with available experimental and analytic results as well as with a numerical solution obtained by other investigators.