Abstract Ice crystal number is a critical ingredient in the potential climate impact of persistent contrails and contrail-induced cirrus. We perform an extensive set of large-eddy simulations (LES) ...of ice nucleation and growth within aircraft exhaust jets with an emphasis on assessing the importance of detailed plume mixing on the effective ice-number emission index (EI iceno ) produced for different conditions. Parameter variations considered include ambient temperature, pressure, and humidity; initial aerosol origin (exhaust or ambient), number, and properties; and aircraft engine size. The LES are performed in a temporal representation with binned microphysics including the basics of activation of underlying aerosol, droplet growth, and freezing. We find that a box-model approach reproduces EI iceno from LES well for sufficiently low aerosol numbers or when crystal production is predominantly on ambient aerosol. For larger exhaust aerosol number the box model generally overestimates EI iceno and can underestimate the fraction from ultrafine aerosol. The effects of different parameters on EI iceno can largely be understood with simpler analytic models that are formulated in low and high aerosol-number limits. The simulations highlight the potential importance of “cold” contrails, ambient ultrafine aerosols, crystal loss due to competition between different-sized crystals, and limitations on reducing EI iceno . We find EI iceno insensitive to engine size for lower aerosol numbers, but decreasing with increasing engine size for higher aerosol numbers. Temporal versus spatial representations for jet LES are compared in an appendix.
Data from the first research flight (RF01) of the second Dynamics and Chemistry of Marine Stratocumulus (DYCOMS-II) field study are used to evaluate the fidelity with which large-eddy simulations ...(LESs) can represent the turbulent structure of stratocumulus-topped boundary layers. The initial data and forcings for this case placed it in an interesting part of parameter space, near the boundary where cloud-top mixing is thought to render the cloud layer unstable on the one hand, or tending toward a decoupled structure on the other hand. The basis of this evaluation consists of sixteen 4-h simulations from 10 modeling centers over grids whose vertical spacing was 5 m at the cloud-top interface and whose horizontal spacing was 35 m. Extensive sensitivity studies of both the configuration of the case and the numerical setup also enhanced the analysis. Overall it was found that (i) if efforts are made to reduce spurious mixing at cloud top, either by refining the vertical grid or limiting the effects of the subgrid model in this region, then the observed turbulent and thermodynamic structure of the layer can be reproduced with some fidelity; (ii) the base, or native configuration of most simulations greatly overestimated mixing at cloud top, tending toward a decoupled layer in which cloud liquid water path and turbulent intensities were grossly underestimated; (iii) the sensitivity of the simulations to the representation of mixing at cloud top is, to a certain extent, amplified by particulars of this case. Overall the results suggest that the use of LESs to map out the behavior of the stratocumulus-topped boundary layer in this interesting region of parameter space requires a more compelling representation of processes at cloud top. In the absence of significant leaps in the understanding of subgrid-scale (SGS) physics, such a representation can only be achieved by a significant refinement in resolution-a refinement that, while conceivable given existing resources, is probably still beyond the reach of most centers.
Abstract
Mobile, polarimetric radar data were collected on a series of tornadoes that occurred near Dodge City, Kansas. A poststorm survey revealed a series of tornadic debris swaths in several dirt ...fields and high-resolution pictures of the tornado documented the visual characteristics of the tornado and the lofted debris cloud. The main rotational couplet associated with the tornado was identified in the single-Doppler velocities; however, no secondary rotational couplets were resolved in the low-level data performed during two consecutive volume scans. Numerical simulations have suggested that cycloidal damage swaths can result when debris is deposited as the low-level inflow turns upward in the corner region of the updraft annulus of the tornado core. This mechanism can dominate even when suction vortices are present in the simulations and can produce these swaths in the absence of these smaller-scale vortices. It is hypothesized that the observed cycloidal damage swaths were a result of the low-level inflow in the corner region of the tornado and not by the existence of suction vortices. Polarimetric data were combined with photographs of the tornado in order to document the lofted debris cloud and its relationship with the funnel. This analysis provided an opportunity to investigate whether recent findings describing the cross-correlation coefficient
ρ
hv
and differential reflectivity
Z
DR
signatures of the lofted debris cloud could be replicated. Regions of low
ρ
hv
at the periphery of the funnel cloud suggesting high debris loading and a column of negative
Z
DR
centered on the tornado believed to be produced by common debris alignment were noted.
Significance Statement
It is well known that some tornadoes produce smaller-scale vortices that rotate around the central axis of the main circulation. In addition, numerous aerial photographs have documented cycloidal debris marks within tornado damage tracks that traverse open fields. The prevailing theory shown in numerous textbooks is that these marks are produced by these vortices. The current study suggests that this widely accepted model for producing these marks may be incorrect. It is suggested that these cycloidal marks are produced by the main tornado circulation and not by the smaller-scale vortices in this case.
Cloud water sedimentation and drizzle in a stratocumulus-topped boundary layer are the focus of an intercomparison of large-eddy simulations. The context is an idealized case study of nocturnal ...stratocumulus under a dry inversion, with embedded pockets of heavily drizzling open cellular convection. Results from 11 groups are used. Two models resolve the size distributions of cloud particles, and the others parameterize cloud water sedimentation and drizzle. For the ensemble of simulations with drizzle and cloud water sedimentation, the mean liquid water path (LWP) is remarkably steady and consistent with the measurements, the mean entrainment rate is at the low end of the measured range, and the ensemble-average maximum vertical wind variance is roughly half that measured. On average, precipitation at the surface and at cloud base is smaller, and the rate of precipitation evaporation greater, than measured. Including drizzle in the simulations reduces convective intensity, increases boundary layer stratification, and decreases LWP for nearly all models. Including cloud water sedimentation substantially decreases entrainment, decreases convective intensity, and increases LWP for most models. In nearly all cases, LWP responds more strongly to cloud water sedimentation than to drizzle. The omission of cloud water sedimentation in simulations is strongly discouraged, regardless of whether or not precipitation is present below cloud base.
Aeolian processes are the most active processes modifying the surface of Mars under present day climatic conditions. Besides wind streak changes and dune and ripple migrations, active dust devils ...occur frequently leaving numerous tracks on the Martian surface. These dust devil tracks (DDTs) are characterized by albedo changes with respect to their surroundings and are suggested to be caused by erosion of dust exposing coarser grained material. Here we show that DDTs with a cycloidal pattern analyzed in situ in southern Peru are formed by erosion of very coarse sands at the outer margins and its subsequent annular deposition in the central parts of dust devils. Field observations are supported by large-eddy simulations using typical dust devil parameters resembling the cycloidal morphology of the DDTs. Cycloidal DDTs observed on Mars resembling the Peruvian DDTs suggest an equivalent formation mechanism. Our results imply that the formation of DDTs on Mars are not solely due to dust erosion but also depositional processes and dust devils are strong enough to redistribute coarser grained material such as sands; hence they might contribute to the modification of the present day Martian landscape.
•In situ analysis of cycloidal dust devil tracks (DDTs) in southern Peru.•Formation by erosion of coarse sands at the outer margins and its subsequent deposition in the central parts.•Field observations are supported by large-eddy simulations.•Cycloidal DDTs observed on Mars resembling the Peruvian DDTs suggest an equivalent formation mechanism.•Dust devils on Mars might be strong enough to redistribute coarse grained material.
Abstract
Past numerical simulation studies found that debris loading from sand-sized particles may substantially affect tornado dynamics, causing reductions in near-surface wind speeds up to 50%. To ...further examine debris loading effects, simulations are performed using a large-eddy simulation model with a two-way drag force coupling between air and sand. Simulations encompass a large range of surface debris fluxes that cause negligible to substantial impact on tornado dynamics for a high-swirl tornado vortex simulation.
Simulations are considered for a specific case with a single vortex flow type (swirl ratio, intensity, and translation velocity) and a fixed set of debris and aerodynamic parameters. Thus, it is stressed that these findings apply to the specific flow and debris parameters herein and would likely vary for different flows or debris parameters. For this specific case, initial surface debris fluxes are varied over a factor of 16 384, and debris cloud mass varies by only 42% of this range because a negative feedback reduces near-surface horizontal velocities. Debris loading effects on the axisymmetric mean flow are evident when maximum debris loading exceeds 0.1 kg kg−1, but instantaneous maximum wind speed and TKE exhibit small changes at smaller debris loadings (greater than 0.01 kg kg−1). Initially, wind speeds are reduced in a shallow, near-surface layer, but the magnitude and depth of these changes increases with higher debris loading. At high debris loading, near-surface horizontal wind speeds are reduced by 30%–60% in the lowest 10 m AGL. In moderate and high debris loading scenarios, the number and intensity of subvortices also decrease close to the surface.
The fifth intercomparison of the Global Water and Energy Experiment Cloud System Studies Working Group 1 is used as a vehicle for better understanding the dynamics of trade wind cumuli capped by a ...strong inversion. The basis of the intercomparison is 10 simulations by 7 groups. These simulations are supplemented by many further sensitivity studies, including some with very refined grid meshes. The simulations help illustrate the turbulent dynamics of trade cumuli in such a regime. In many respects the dynamics are similar to those found in many previous simulations of trade cumuli capped by weaker inversions. The principal differences are the extent to which the cloud layer is quasi-steady in the current simulations, evidence of weak countergradient momentum transport within the cloud layer, and the development and influence of an incipient stratiform cloud layer at the top of the cloud layer. Although many elements of the turbulent structure (including the wind profiles, the evolution of cloud-base height, the statistics of the subcloud layer, and the nature of mixing in the lower and middle parts of the cloud layer) are robustly predicted, the representation of the stratiform cloud amount by the different simulations is remarkably sensitive to a number of factors. Chief among these are differences between numerical algorithms. These sensitivities persist even among simulations on relatively refined grid meshes. Part of this sensitivity is attributed to a physically realistic positive radiative feedback, whereby a propensity toward higher cloud fractions in any given simulation is amplified by longwave radiative cooling. The simulations also provide new insight into the dynamics of the transition layer at cloud base. In accord with observations, the simulations predict that this layer is most identifiable in terms of moisture variances and gradients. The simulations help illustrate the highly variable (in both height and thickness) nature of the transition layer, and we speculate that this variability helps regulate convection. Lastly the simulations are used to help evaluate simple models of trade wind boundary layers. In accord with previous studies, mass-flux models well represent the dynamics of the cloud layer, while mixing-length models well represent the subcloud layer. The development of the stratiform cloud layer is not, however, captured by the mass-flux models. The simulations indicate that future theoretical research needs to focus on interface rules, whereby the cloud layer is coupled to the subcloud layer below and the free atmosphere above. Future observational studies of this regime would be of most benefit if they could provide robust cloud statistics as a function of mean environmental conditions.
The results of high-resolution, fully three-dimensional, unsteady simulations of the interaction of a tornado with the surface are presented. The goal is to explore some of the range of structures ...that should be expected to occur in nature within the tornadic ``corner flow''--that region where the central vortex meets the surface. The most important physical variables considered are the tornado-scale circulation and horizontal convergence, the effective surface roughness, the tornado translation speed, the low-level inflow structure, and the upper-core structure. A key ingredient of the corner flow dynamics is the radial influx of fluid in the surface layer with low angular momentum relative to that of the fluid in the main vortex above it. This low swirl fluid arises initially from outside or below the larger-scale vortex or through frictional loss of angular momentum to the surface and forms much of the vortex core flow after it exits the corner flow region. Changes in the surface layer inflow or upper-core structure can dramatically affect the level of intensification and turbulent structure in the corner flow even when the swirl ratio of the tornado vortex as a whole is unchanged. The authors define a local corner flow swirl ratio, S sub(c) , based on the total flux of low angular momentum fluid through the corner flow and show that it parameterizes the leading effects on the corner flow of changes to the flow conditions immediately outside of the corner flow. As S sub(c) decreases, the low-level vortex intensity rises to a maximal level where mean swirl velocities near the surface reach 2.5 times the maximum mean swirl velocity aloft; further decreases force a transition to a much weaker low-level tornado vortex. This sensitivity suggests that differences in the near-surface inflow layer may be a critical factor in determining whether an existing supercell low-level mesocyclone spawns a tornado or not.
This paper reports an intercomparison study on undisturbed trade wind cumulus convection under steady-state conditions as observed during the Barbados Oceanographic and Meteorological Experiment ...(BOMEX) with 10 large eddy simulation (LES) models. A main objective of this study is to obtain a quantitative assessment of the quality of the turbulent dynamics for this type of boundary layer clouds as produced by the different LES codes. A 6-h simulation shows excellent model-to-model agreement of the observed vertical thermodynamical structure, reasonable agreement of variances and turbulent fluxes, and good agreement of quantities conditionally sampled within the model clouds, such as cloud cover, liquid water, and cloud updraft strength. In the second part of this paper the LES dataset is used to evaluate simple models that are used in parameterizations of current general circulation models (GCMs). Finally, the relation of this work to subsequent LES studies of more complicated regimes is discussed, and guidance is given for the design of future observational studies of shallow cumulus boundary layers.