Evaluation of deep convective transport in storms from different convective regimes during the DC3 field campaign using WRF‐Chem with lightning data assimilation
Explore the source record for details and available documents.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Convective updrafts are one of the main characteristics of convective clouds, responsible for the convective mass flux and the redistribution of energy and condensate in the atmosphere. During the early stages of their lifecycle, convective clouds experience rapid cloud-top ascent manifested by a decrease in the geostationary IR brightness temperature (𝑇𝐵 𝐼𝑅 ). Under the assumption that the convective cloud top behaves like a black body, the ascent rate of the convective cloud top can be estimated as ($\frac{∂𝑇𝐵_{𝐼𝑅}}{∂𝑡}$), and it can be used to infer the near cloud-top convective updraft. The temporal resolution of the geostationary IR measurements and non-uniform beam-filling effects can influence the convective updraft estimation. However, the main shortcoming until today was the lack of independent verification of the strength of the convective updraft. Here, Doppler radar observations from the ESCAPE and TRACER field experiments provide independent estimates of the convective updraft velocity at higher spatiotemporal resolution throughout the convective core column and can be used to evaluate the updraft velocity estimates from the IR cooling rate for limited samples. Isolated convective cells were tracked with dedicated radar (RHIs and PPIs) scans throughout their lifecycle. Radial Doppler velocity measurements near the convective cloud top are used to provide estimates of convective updrafts. These data are compared with the geostationary IR and VIS channels (from the GOES satellite) to characterize the convection evolution and lifecycle based on cloud-top cooling rates.
Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.
It is well known that convection is a dominant mass transport mechanism when materials are solidified on Earth's surface. This convection is caused by gradients in density (and therefore gravitational force) that are brought about by gradients in temperature, composition or both. Diffusion of solute is therefore dwarfed by convection and the study of fundamental parameters, such as dendrite tip shape and growth velocity in the absence of convection is nearly impossible. Significant experimental work has therefore been carried out in orbiting laboratories with the intent of minimizing convection by minimizing gravity. One of the best known experiments of this kind is the Isothermal Dendritic Growth Experiment (IDGE), supported by NASA. Naturally such experiments are costly and one objective of the present investigation is to develop an experimental method whereby convection can be- halted, in solidification and other experiments, on the surface. A second objective is to use the method to minimize convection resulting from the residual accelerations suffered by experiments in microgravity. The method to be used to minimize convection relies on the dependence of the magnetic susceptibility of a fluid on temperature or composition (whichever is driving convection). All materials experience a force when placed in a magnetic field gradient. The direction and magnitude of that force depend on the magnetic susceptibility of the material. Consequently the force will vary if the susceptibility varies with temperature or composition. With a magnetic field gradient in the right direction (typically upward) and of the right magnitude, this variation in the magnetic force can be made to exactly cancel the variation in the gravitational force. Expressed another way, normal buoyancy is exactly countered by a "magnetic buoyancy". To demonstrate the principle, a solution of MnC12 in water has been used. First the variation of the susceptibility of this paramagnetic solution with temperature and concentration was measured. Then a "cell", containing this solution and 50mm long by 15mm high by 155mm wide, was placed in a superconducting magnet at Marshall Space Flight Center. The magnetic field was measured at various positions within the bore of the magnet using a Hall effect probe. In this way, a position was found where the magnetic field gradient was predominantly upward; the magnitude of the gradient could then be adjusted by adjusting the current of the magnet. The ends of the cell consisted of machined copper blocks maintained at controlled temperatures by circulating water from constant temperature baths. The walls of the cell were of rectangular section glass tubing so that the cell contents could be seen. Velocities arising from thermal gradients within the cell were measured by particle image velocimetry (PIV). Particles used for this purpose were silver-coated hollow glass spheres of micrometers diameter and nearly the same density as the solution. A central vertical plane of the cell was illuminated by a laser beam passing through a cylindrical lens. Digital images of the particles were captured on a CCD camera and fed to a computer so that frame-to-frame movements of particles traveling with the fluid were captured. These images were employed to compute velocity maps using commercial PIV software. In a typical experiment the cold end of the cell was maintained at 10C and the warm end at 30 C. With no current in the magnet, i.e.- with natural -convection allowed to occur, the fluid was observed to circulate with an average speed of approximately 0.3 millimeters per second. It was visually apparent that this circulation was diminished as the current was increased. At currents of approximately 20A the flow was halted, to within the precision of the PIV measurements. At yet higher currents the convection was reversed with the hotter solution sinking and the cooler solution rising. At 40A this reversed convection had speeds averaging 0.43 millimeters per second. The measurements of susceptibility and density allow an estimate of the field gradient necessary to halt convection in the experiment. That estimate was 7.8T (squared) per meter and the convection was observed to halt in the magnet at a current giving 7.21T (squared) per meter from the magnetic field measurements. Calculations of the flow have been carried out using the computational fluid dynamics software FLUENT and show good agreement with the measurements.
The authors analyze the influence of Sea Surface Temperature (SST) and surface wind divergence on atmospheric thermodynamic structure and the resulting effects on the occurrence of deep convection using National Meteorological Center radiosonde data and International Satellite Cloud Climatology Program data for July 1983-July 1985. The onset of deep convection requires not only the existence of positive convective available potential energy (CAPE), but also an unstable planetary boundary layer (PBL). A stable PBL is observed to suppress deep convection even when CAPE is positive. Variations of SST have a major effect on CAPE, but surface wind divergence can also affect deep convection by changing the lapse rate in the lower troposphere and humidity in the PBL. Specifically, when SST is greater than or equal to 28 C, CAPE is always positive, and surface wind divergence does not qualitatively change the buoyancy profile above the PBL. Strong surface wind divergence, however, stabilizes the PBL so as to suppress the initiation of deep convection. In warm SST regions, CAPE is greater than 0 regardless of assumptions about condensate loading, although the pseudoadiabatic limit is more consistent with the observed deep convection than the reversible moist-adiabatic limit under these circumstances. When SST is less than 27 C, CAPE is usually negative and inhibits convection, but strong surface wind convergence can destabilize the inversion layer and moisten the PBL enough to make the atmosphere neutrally stable in the mean. As a result, deep convection is generally enhanced either when SST is greater than or equal to 28 C in the absence of strong surface wind divergence or when strong surface wind convergence occurs even if SST is less than 27 C. The anomalous suppression of deep convection in the warm area of the equatorial west Pacific lying between the intertropical convergence zone (ITCZ) and south Pacific convergence zone (SPCZ) is probably caused by dryness in the PBL and an inversion in that area. The seasonal cycles of deep convection and surface wind divergence are in phase with the maximum solar radiation and lead SST for one to three months in the central Pacific. The change of PBL relative humidity plays a critical role in the changeover to convective instability in this case. The seasonal change of deep convection and associated clouds seems not to have important effects on the seasonal change of local SST in the central Pacific.
Context. Theoretical descriptions of convective overshooting in stellar interiors often rely on a basic one-dimensional parameterization of the flow called the filling factor for convection. Several different definitions of the filling factor have been developed for this purpose, based on: (1) the percentage of the volume, (2) the mass flux, and (3) the convective flux that moves through the boundary. Aims. We examine these definitions of the filling factor with the goal of establishing their ability to explain differences between 2D and 3D global simulations of stellar interiors that include fully compressible hydrodynamics and realistic microphysics for stars. Methods. We study convection and overshooting in pairs of identical two-dimensional (2D) and three-dimensional (3D) global simulations of stars produced with MUSIC , a fully compressible, time-implicit hydrodynamics code. We examine pairs of simulations for (1) a 3 M ⊙ red giant star near the first dredge-up point, (2) a 1 M ⊙ pre-main-sequence star with a large convection zone, (3) the current sun, and (4) a 20 M ⊙ main-sequence star with a large convective core. Results. Our calculations of the filling factor based on the volume percentage and the mass flux indicate asymmetrical convection near the surface for each star with an outer convection zone. However, near the convective boundary, convective flows achieve inward-outward symmetry for each star that we study; for 2D and 3D simulations, these filling factors are indistinguishable. A filling factor based on the convective flux is contaminated by boundary-layer-like flows, making a theoretical interpretation difficult. We present two possible new alternatives to these frequently used definitions of a filling factor, which instead compare flows at two different radial points. The first alternative is the penetration parameter of Anders et al. (2022, ApJ, 926, 169). The second alternative is a new statistic that we call the plume interaction parameter. We demonstrate that both of these parameters captures systematic differences between 2D and 3D simulations around the convective boundary.
The impact of convection on the humidity and clouds in the tropical tropopause layer (TTL) during boreal summer 2007 is investigated in simulations of detailed cloud microphysical processes and their effects on the water vapor (H2O) profile along backward trajectories from the 379 K potential temperature (100‐hPa pressure) surface. Convective influence is determined by tracing the trajectories through time‐dependent fields of satellite‐based convective cloud top height. The simulated H2O mixing ratios at the 100‐hPa level and cloud occurrence fractions in the middle to upper (16–18 km) TTL exhibit a pronounced maximum over the Asian monsoon region as in observations; these local enhancements are virtually absent in the simulation without convection, indicating that convection is the dominant driver of the localized H2O and cloud maxima in the Asian summer monsoon region. Convection moistens the 100‐hPa level by 0.6 ppmv (~15%) averaged over the 10°S–50°N domain and increases tropical (10°S–30°N) mean cloud occurrence in the middle to upper TTL by ~170%. Nearly all of the convective enhancements in H2O and clouds are due to the effect of convective saturation; convectively detrained ice crystals have negligible impact. Parcels are most frequently hydrated by deep convection in the southern sector of the Asian monsoon anticyclone and subsequently dehydrated downstream of convection to the west, shifting the locations of final dehydration northwest of the cold temperature region in the northern Tropics. Infrequent, extreme deep convective systems (cloud tops exceeding 380 K) have a disproportionately large effect on TTL humidity and clouds.
Water vapor in the stratosphere is primarily controlled by temperatures in the tropical upper troposphere and lower stratosphere. However, the direct impact of deep convection on the global lower stratospheric water vapor budget is still an actively debated issue. Two complementary modeling approaches are used to investigate the convective impact in boreal winter and summer. Convective influence is diagnosed by tracing trajectories through convective cloud top altitude fields derived from global rainfall and brightness temperature data. Backward trajectory model simulations coupled with a detailed treatment of cloud microphysical processes indicate that convection moistens the global lower stratosphere by approximately 0.3 ppmv in boreal winter and summer 2010. The diurnal peak in convection is responsible for about half of the total convective moistening during winter and nearly all of the convective moistening during summer. Deep convective clouds overshooting the tropopause have relatively minor effect on global lower stratospheric water vapor. A forward trajectory model coupled with a simplified cloud module is used to estimate the relative magnitude of the interannual variability of the convective impact. Combining the results from the two models, we find that the convective impact on the global lower stratospheric water vapor during 2006-2016 is approximately 0.3 ppmv with year-to-year variations of up to 0.1 ppmv. An important mechanism of convective hydration of the lower stratosphere is via the detrainment of saturated air and ice into the tropical uppermost troposphere and the subsequent upward transport of some of these moist air parcels across relatively warm and subsaturated tropopause.
Abstract This study evaluates convective cell properties and their relationships with convective and stratiform rainfall within a season‐long convection‐permitting weather research and forecasting simulation over central Argentina using radar, satellite, and radiosonde measurements from the RELAMPAGO‐CACTI field campaign. The simulation slightly underestimates radar‐estimated rainfall over the ∼3.5‐month evaluation period but underestimates stratiform rainfall by 46% and overestimates convective rainfall by 43%. As convective available potential energy (CAPE) increases, the convective rainfall overestimation decreases, but the stratiform rainfall underestimation increases such that the contribution of convective to total rainfall remains constantly high biased by ∼26%. Overestimated convective rainfall arises from the simulation generating 2.6 times more precipitating convective cells (14,299) than observed by radar (5,662) despite similar observed and simulated cell growth processes, with relatively wide cells contributing mostly to excessive convective rainfall. Relatively shallow cells, typically reaching heights of 4–7 km, contribute most to the cell number bias. This cell number bias increases as CAPE decreases, potentially because cells and their updrafts become narrower and more under‐resolved as CAPE decreases. The gross overproduction of precipitating shallow cells leads to overly efficient precipitation and inadequate detrainment of ice aloft, thereby diminishing the formation of robust stratiform rainfall regions. Decreasing model horizontal grid spacing from 3 to 1 or 0.333 km for low (<300 J kg −1 ) and high CAPE (>1,000 J kg −1 ) cases results in minimal change to cell number, depth, and convective‐to‐stratiform partitioning biases. This suggests that improving prediction of these convective properties depends on factors beyond solely increasing model resolution.
This presentation will trace the important research developments of the last 20+ years in defining the roles of deep convection in tropospheric chemistry. The role of deep convection in vertically redistributing trace gases was first verified through field experiments conducted in 1985. The consequences of deep convection have been noted in many other field programs conducted in subsequent years. Modeling efforts predicted that deep convection occurring over polluted continental regions would cause downstream enhancements in photochemical ozone production in the middle and upper troposphere due to the vertical redistribution of ozone precursors. Particularly large post-convective enhancements of ozone production were estimated for convection occurring over regions of pollution from biomass burning and urban areas. These estimates were verified by measurements taken downstream of biomass burning regions of South America. Models also indicate that convective transport of pristine marine boundary layer air causes decreases in ozone production rates in the upper troposphere and that convective downdrafts bring ozone into the boundary layer where it can be destroyed more rapidly. Additional consequences of deep convection are perturbation of photolysis rates, effective wet scavenging of soluble species, nucleation of new particles in convective outflow, and the potential fix stratosphere-troposphere exchange in thunderstorm anvils. The remainder of the talk will focus on production of NO by lightning, its subsequent transport within convective clouds . and its effects on downwind ozone production. Recent applications of cloud/chemistry model simulations combined with anvil NO and lightning flash observations in estimating NO Introduction per flash will be described. These cloud-resolving case-study simulations of convective transport and lightning NO production in different environments have yielded results which are directly applicable to the design of lightning parameterizations for global chemical transport models. The range of mean values (factor of 3) of NO production per flash (or per meter of lightning channel length) that have been deduced from the model will be shown and compared with values of production in the literature that have been deduced using other methods, Results show that on a per flash basis, IC flashes are nearly as productive of NO as CG flashes. When combined with the global flash rate of 44 flashes per second from NASA's Optical Transient Detector (OTD) measurements, these estimates and the results from other techniques yield global NO production rates of 2-9 TgN/year. Vertical profiles of lightning NOx mass at the end of the 3-D storm simulations have been summarized to yield suggested profiles for use in global models. Simulations of the photochemistry over the 24 hours following a storm have been performed to determine the additional ozone production which can be attributed to lightning NO.
Entrainment of surrounding cooler and drier air into convective updrafts is one of the key processes that influence deep convection initiation and growth. Numerous studies have investigated the effect of entrainment on isolated convective cloud growth in idealized simulations, but the importance of this effect in realistic conditions with many interacting convective clouds remains uncertain. We examine the impact of entrainment on the depth reached by convective clouds in realistic large-eddy simulations (LES) over central Argentina during the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign. Cloudy updrafts and their associated properties are assigned to convective cells tracked with radar reflectivity signatures. Several thousand convective cells are tracked over two high convective available potential energy (CAPE) and two low CAPE cases that support cells of varying depths and intensities. Entrainment is calculated explicitly as the fluxes of air into the outer surface of each cloudy updraft. Single-predictor logistic regression models are used to determine the relative importance of updraft, near-updraft, and preconvective initiation atmospheric conditions in predicting whether convective cells become deep. We then build a multiple-predictor regression model pairing important updraft and meteorological metrics with fractional entrainment rate. The probability of cells transitioning to deep convection is most sensitive to ambient 600-hPa relative humidity (42% of total metric contribution to cloud depth predictability), followed by low-level CAPE (28%), cloud-base updraft width (19%), and fractional entrainment (11%). Thus, the initial width of the updraft along with potential buoyancy and its dilution through the midtroposphere collectively determine whether deep convection will result from shallower clouds.
Atmospheric aerosols affect the global energy budget by scattering and absorbing sunlight (direct effects) and by changing the microphysical structure, lifetime, and coverage of clouds (indirect effects). Globally, the free troposphere is a major source of nucleation- and Aitken-mode aerosols due to the enhanced new particle formation rates at high altitudes. Recent studies have shown deep convective systems are capable of transporting these small aerosols from the free troposphere to the boundary layer by strong convective downdrafts and weaker downward motions in the stratiform regions. These vertically transported aerosols can grow into cloud condensation nuclei (CCN) and play a significant role in the global climate. During the deep convective processes, existing accumulation-mode aerosols that act as coagulation sinks of smaller particles are also removed by wet scavenging. Compared to the vertical transport of these particles by entrainment mixing, which is slower but more prevalent, the deep convective downdraft processes may be more rapid and efficient in the vertical transport of aerosols. However, most of the current climate models do not include this mechanism as a source of CCN, mainly because the frequency of deep convective events varies significantly with geographic location and thus their contributions to CCN are unpredictable. We target this critical gap in understanding the vertical transport and removal of aerosols by deep convections. We proposed to analyze a multi-year, multi-site measurement record available from the U.S. Department of Energy (DOE) ARM program, including the observations from the 2014/15 Observations and Modeling of the Green Ocean Amazon (GoAmazon) field campaign, the 2017/18 Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) field campaign, the 2018/19 Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign, the 2021/22 Tracking Aerosol Convection Interactions Experiment (TRACER) field campaign, and the long-term measurements collected at the Southern Great Plains (SGP) atmospheric observatory, where deep convective clouds were frequently observed. This project is aimed at the following three objectives: (1) Gaining a detailed and quantitative understanding of the aerosols transported by a convective downdraft and their evolution in the atmosphere; (2) Examining the wet scavenging mechanisms and efficiencies of aerosols at altitudes of deep convective systems based on ground and aircraft measurements; (3) Evaluating the contribution of deep convective systems to CCN as both a source and a sink of atmospheric aerosols and its seasonal variabilities.
Models of varying complexity that simulate water vapor and clouds in the Tropical Tropopause Layer (TTL) show that including convection directly is essential to properly simulating the water vapor and cloud distribution. In boreal winter, for example, simulations without convection yield a water vapor distribution that is too uniform with longitude, as well as minimal cloud distributions. Two things are important for convective simulations. First, it is important to get the convective cloud top potential temperature correctly, since unrealistically high values (reaching above the cold point tropopause too frequently) will cause excessive hydration of the stratosphere. Second, one must capture the time variation as well, since hydration by convection depends on the local relative humidity (temperature), which has substantial variation on synoptic time scales in the TTL. This paper describes a method for obtaining high frequency (3-hourly) global convective cloud top distributions which can be used in trajectory models. The method uses rainfall thresholds, standard IR brightness temperatures, meteorological temperature analyses, and physically realistic and documented corrections IR brightness temperature corrections to derive cloud top altitudes and potential temperatures. The cloud top altitudes compare well with combined CLOUDSAT and CALIPSO data, both in time-averaged overall vertical and horizontal distributions and in individual cases (correlations of .65-.7). An important finding is that there is significant uncertainty (nearly .5 km) in evaluating the statistical distribution of convective cloud tops even using lidar. Deep convection whose tops are in regions of high relative humidity (such as much of the TTL), will cause clouds to form above the actual convection. It is often difficult to distinguish these clouds from the actual convective cloud due to the uncertainties of evaluating ice water content from lidar measurements. Comparison with models show that calculated cloud top altitudes are generally higher than those calculated by global analyses (e.g., MERRA). Interannual variability in the distribution of convective cloud top altitudes is also investigated.