It is often assumed in parcel theory calculations, numerical models, and cumulus parameterizations that moist static energy (MSE) is adiabatically conserved. However, the adiabatic conservation of MSE is only approximate because of the assumption of hydrostatic balance. Two alternative variables are evaluated here: MSE − IB and MSE + KE, wherein IB is the path integral of buoyancy (
Common assumptions in temperature lapse rate formulas for lifted air parcels include neglecting mixing, hydrostatic balance, the removal of all condensate once it forms (pseudoadiabatic), and/or the retention of all condensate within the parcel (adiabatic). These formulas are commonly derived from the conservation of entropy, which leads to errors when nonequilibrium mixed-phase condensate is present. To evaluate these assumptions, a new general lapse rate formula is derived from an expression for energy conservation, rather than entropy conservation. This new formula incorporates mixing of the parcel with its surroundings, relaxes the hydrostatic assumption, allows for nonequilibrium mixed-phase condensate, and can be formulated for pseudoadiabatic or adiabatic ascent. The new formula is shown to exactly conserve entropy for reversible ascent. Predictions by the new formula are compared to that of older and less general formulas. The errors in previous formulas arise from the assumption of hydrostatic balance, which results in considerable warm biases due to the neglect of the energy sink from buoyancy. Predictions of ascent with entrainment using the new formula are then compared to parcel properties along trajectories in large eddy simulations. Simulated parcel properties are better predicted by the formula using a diluted analogy to adiabatic ascent, wherein condensate is diluted at the same rate as other parcel properties, than by the diluted analogy to pseudoadiabatic ascent, wherein all condensate is removed. These results suggest that CAPE should be computed with adiabatic, rather than pseudoadiabatic, parcel ascent.
more » « less- Award ID(s):
- 1648681
- NSF-PAR ID:
- 10364310
- Publisher / Repository:
- American Meteorological Society
- Date Published:
- Journal Name:
- Journal of the Atmospheric Sciences
- Volume:
- 79
- Issue:
- 3
- ISSN:
- 0022-4928
- Format(s):
- Medium: X Size: p. 815-836
- Size(s):
- p. 815-836
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract B ) and KE is kinetic energy. Both of these variables relax the hydrostatic assumption and are more precisely conserved than MSE. This article quantifies the errors that result from assuming that the aforementioned variables are conserved in large-eddy simulations (LES) of both disorganized and organized deep convection. Results show that both MSE − IB and MSE + KE better predict quantities along trajectories than MSE alone. MSE − IB is better conserved in isolated deep convection, whereas MSE − IB and MSE + KE perform comparably in squall-line simulations. These results are explained by differences between the pressure perturbation behavior of squall lines and isolated convection. Errors in updraftB diagnoses are universally minimized when MSE − IB is assumed to be adiabatically conserved, but only when moisture dependencies of heat capacity and temperature dependency of latent heating are accounted for. When less accurate latent heat and heat capacity formulae were used, MSE − IB yielded poorerB predictions than MSE due to compensating errors. Our results suggest that various applications would benefit from using either MSE − IB or MSE + KE instead of MSE with properly formulated heat capacities and latent heats. -
more » « less
Abstract This article introduces an analytic formula for entraining convective available potential energy (ECAPE) with an entrainment rate that is determined directly from an environmental sounding, rather than prescribed by the formula user. Entrainment is connected to the background environment using an eddy diffusivity approximation for lateral mixing, updraft geometry assumptions, and mass continuity. These approximations result in a direct correspondence between the storm-relative flow and the updraft radius and an inverse scaling between the updraft radius squared and entrainment rate. The aforementioned concepts, combined with the assumption of adiabatic conservation of moist static energy, yield an explicit analytic equation for ECAPE that depends entirely on state variables in an atmospheric profile and a few constant parameters with values that are established in past literature. Using a simplified Bernoulli-like equation, the ECAPE formula is modified to account for updraft enhancement via kinetic energy extracted from the cloud’s background environment. CAPE and ECAPE can be viewed as predictors of the maximum vertical velocity
w maxin an updraft. Hence, these formulas are evaluated usingw maxfrom past numerical modeling studies. Both of the new formulas improve predictions ofw maxsubstantially over commonly used diagnostic parameters, including undiluted CAPE and ECAPE with a constant prescribed entrainment rate. The formula that incorporates environmental kinetic energy contribution to the updraft correctly predicts instances of exceedance ofby w max, and provides a conceptual explanation for why such exceedance is rare among past simulations. These formulas are potentially useful in nowcasting and forecasting thunderstorms and as thunderstorm proxies in climate change studies.Significance Statement Substantial mixing occurs between the upward-moving air currents in thunderstorms (updrafts) and the surrounding comparatively dry environmental air, through a process called entrainment. Entrainment controls thunderstorm intensity via its diluting effect on the buoyancy of air within updrafts. A challenge to representing entrainment in forecasting and predictions of the intensity of updrafts in future climates is to determine how much entrainment will occur in a given thunderstorm environment without a computationally expensive high-resolution simulation. To address this gap, this article derives a new formula that computes entrainment from the properties of a single environmental profile. This formula is shown to predict updraft vertical velocity more accurately than past diagnostics, and can be used in forecasting and climate prediction to improve predictions of thunderstorm behavior and impacts.
-
more » « less
Abstract Are the results of aerosol invigoration studies that neglect entrainment valid for diluted deep convective clouds? We address this question by applying an entraining parcel model to soundings from tropical and midlatitude convective environments, wherein pollution is assumed to increase parcel condensate retention. Invigoration of 5%–10% and <2% is possible in undiluted tropical and midlatitude parcels respectively when freezing is rapid. This occurs because the positive buoyancy contribution from freezing is larger than the negative buoyancy contribution from condensate loading, leading to positive net condensate contribution to buoyancy. However, aerosol‐induced weakening is more likely when realistic entrainment rates occur because water losses from entrainment more substantially reduce the latent heating relative to the loading contribution. This leads to larger net negative buoyancy contribution from condensates in polluted than in clean entraining parcels. Our results demonstrate that accounting for entrainment is critical in conceptual models of aerosol indirect effects in deep convection.
-
more » « less
The one‐dimensional steady state analytical solution of the energy conservation equation obtained by Robin (1955, https://doi.org/10.3189/002214355793702028) is frequently used in glaciology. This solution assumes a linear change in surface velocity from a minimum value equal to minus the mass balance at the surface to zero at the bed. Here we show that this assumption of a linear velocity profile leads to large errors in the calculated temperature profile and especially in basal temperature. By prescribing a nonlinear power function of elevation above the bed for the vertical velocity profile arising from use of the Shallow Ice Approximation, we derive a new analytical solution for temperature. We show that the solution produces temperature profiles identical to numerical temperature solutions with the Shallow Ice Approximation vertical velocity near ice divides. We quantify the importance of strain heating and demonstrate that integrating the strain heating and adding it to the geothermal heat flux at the bed is a reasonable approximation for the interior regions. Our analytical solution does not include horizontal advection components, so we compare our solution with numerical solutions of a two‐dimensional advection‐diffusion model and assess the applicability and errors of the analytical solution away from the ice divide. We show that several parameters and assumptions impact the spatial extent of applicability of the new solution including surface mass balance rate and surface temperature lapse rate. We delineate regions of Greenland and Antarctica within which the analytical solution at any depth is likely within 2 K of the actual temperatures with horizontal advection.
-
The time-dependent freezing rate (TDFR) model here described represents the formation of ice particles by immersion freezing within an air parcel. The air parcel trajectory follows an adiabatic ascent and includes a period in time when the parcel remains stationary at the top of its ascent. The description of the ice nucleating particles (INPs) in the air parcel is taken from laboratory experiments with cloud and precipitation samples and is assumed to represent the INP content of the cloud droplets in the parcel. Time dependence is included to account for variations in updraft velocity and for the continued formation of ice particles under isothermal conditions. The magnitudes of these factors are assessed on the basis of laboratory measurements. Results show that both factors give rise to three-fold variations in ice concentration for a realistic range of the input parameters. Refinements of the parameters specifying time dependence and INP concentrations are needed to make the results more specific to different atmospheric aerosol types. The simple model framework described in this paper can be adapted to more elaborate cloud models. The results here presented can help guide decisions on whether to include a time-dependent ice nucleation scheme or a simpler singular description in models.more » « less
