Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atmospheric Science”

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.

At least 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Multiple Environmental Influences on the Lightning of Cold-Based Continental Convection. Part II: Sensitivity Tests for Its Charge Structure and Land–Ocean Contrast

In Part I, an electrification scheme was described and a simulation of an observed cold-based storm from the U.S. Great Plains was validated with electrical observations. Most charge in the storm was separated by rebounding collisions of secondary ice originating from prior graupel–snow collisions. In this Part II, sensitivity tests are performed with the control simulation (Part I) and influences from environmental factors (aerosols, temperature, moisture, and shear) on lightning are elucidated. Environmental factors [e.g., convective available potential energy (CAPE)] controlling updraft speed are salient. When ascent is reduced by 30% and 70%, flashes become 70% fewer and disappear, respectively; faster ascent promotes positive cloud-to-ground (+CGs) flashes. Since cloud base is too cold (1°C) for coalescence, cloud condensation nucleus aerosol concentrations do not influence the lightning appreciably. The electrical response to varying concentrations of active ice nuclei is limited by most ice particles being secondary and less sensitive—a natural “buffer.” Imposing a maritime sounding suggests that the land–sea contrast in lightning for such storms is due to the vertical structure of environmental temperature and humidity. Weak CAPE, and both entrainment and condensate weight from a low cloud base, suppress ascent and charging. Maritime thermodynamic conditions reduce simulated flash rates by two orders of magnitude. Reducing aerosol loadings from continental to maritime only slightly reinforces this suppression. Last, a conceptual model is provided for how any simulated storm is either normal because graupel/hail is mostly positively charged or else is inverted/anomalous because graupel/hail is mostly negatively charged, with environmental factors controlling the charging. Impacts from microphysical processes, including three processes of secondary ice production, on lightning are analyzed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Impact of Cloud-Base Turbulence on CCN Activation: Single-Size CCN

Abstract This paper examines the impact of cloud-base turbulence on activation of cloud condensation nuclei (CCN). Following our previous studies, we contrast activation within a nonturbulent adiabatic parcel and an adiabatic parcel filled with turbulence. The latter is simulated by applying a forced implicit large-eddy simulation within a triply periodic computational domain of 64 3 m 3 . We consider two monodisperse CCN. Small CCN have a dry radius of 0.01 μ m and a corresponding activation (critical) radius and critical supersaturation of 0.6 μ m and 1.3%, respectively. Large CCN have a dry radius of 0.2 μ m and feature activation radius of 5.4 μ m and critical supersaturation 0.15%. CCN are assumed in 200-cm −3 concentration in all cases. Mean cloud-base updraft velocities of 0.33, 1, and 3 m s −1 are considered. In the nonturbulent parcel, all CCN are activated and lead to a monodisperse droplet size distribution above the cloud base, with practically the same droplet size in all simulations. In contrast, turbulence can lead to activation of only a fraction of all CCN with a nonzero spectral width above the cloud base, of the order of 1 μ m, especially in the case of small CCN and weak mean cloud-base ascent. We compare our results to studies of the turbulent single-size CCN activation in the Pi chamber. Sensitivity simulations that apply a smaller turbulence intensity, smaller computational domain, and modified initial conditions document the impact of specific modeling assumptions. The simulations call for a more realistic high-resolution modeling of turbulent cloud-base activation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Disentangling the Influences of Storm-Relative Flow and Horizontal Streamwise Vorticity on Low-Level Mesocyclones in Supercells

Sufficient low-level storm-relative flow is a necessary ingredient for sustained supercell thunderstorms and is connected to supercell updraft width. Assuming a supercell exists, the role of low-level storm-relative flow in regulating supercells’ low-level mesocyclone intensity is less clear. One possibility considered in this article is that storm-relative flow controls mesocyclone and tornado width via its modulation of overall updraft extent. This hypothesis relies on a previously postulated positive correspondence between updraft width, mesocyclone width, and tornado width. An alternative hypothesis is that mesocyclone characteristics are primarily regulated by horizontal streamwise vorticity irrespective of storm-relative flow. A matrix of supercell simulations was analyzed to address the aforementioned hypotheses, wherein horizontal streamwise vorticity and storm-relative flow were independently varied. Among these simulations, mesocyclone width and intensity were strongly correlated with horizontal streamwise vorticity, and comparatively weakly correlated with storm-relative flow, supporting the second hypothesis. Accompanying theory and trajectory analysis offers the physical explanation that, when storm-relative flow is large and updrafts are wide, vertically tilted streamwise vorticity is projected over a wider area but with a lesser average magnitude than when these parameters are small. These factors partially offset one another, degrading the correspondence of storm-relative flow with updraft circulation and rotational velocity, which are the mesocyclone attributes most closely tied to tornadoes. In conclusion, these results refute the previously purported connections between updraft width, mesocyclone width, and tornado width, and emphasize horizontal streamwise vorticity as the primary control on low-level mesocyclones in sustained supercells.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Removing the Effects of Tropical Dynamics from North Pacific Climate Variability

Teleconnections from the Tropics energize variations of the North Pacific climate, but detailed diagnosis of this relationship has proven difficult. Simple univariate methods, such as regression on El Niño-Southern Oscillation (ENSO) indices, may be inadequate since the key dynamical processes involved -- including ENSO diversity in the Tropics, re-emergence of mixed layer thermal anomalies, and oceanic Rossby wave propagation in the North Pacific -- have a variety of overlapping spatial and temporal scales. Here we use a multivariate Linear Inverse Model to quantify tropical and extra-tropical multi-scale dynamical contributions to North Pacific variability, in both observations and CMIP6 models. In observations, we find that the Tropics are responsible for almost half of the seasonal variance, and almost three quarters of the decadal variance, along the North American coast and within the subtropical front region northwest of Hawaii. SST anomalies that are generated by local dynamics within the Northeast Pacific have much shorter time scales, consistent with transient weather forcing by Aleutian low anomalies. Variability within the Kuroshio-Oyashio Extension (KOE) region is considerably less impacted by the Tropics, on all time scales. Consequently, without tropical forcing the dominant pattern of North Pacific variability would be a KOE pattern, rather than the Pacific Decadal Oscillation (PDO). In contrast to observations, most CMIP6 historical simulations produce North Pacific variability that maximizes in the KOE region, with amplitude significantly higher than observed. Correspondingly, the simulated North Pacific in all CMIP6 models is shown to be relatively insensitive to the Tropics, with a dominant spatial pattern generally resembling the KOE pattern, not the PDO.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Polarimetric Radar Variables in Snowfall at Ka- and W-Band Frequency Bands: A Comparative Analysis

Dual-frequency millimeter-wavelength radar observations in snowfall are analyzed in order to evaluate differences in conventional polarimetric radar variables such as differential reflectivity (ZDR) specific differential phase shift (K DP ) and linear depolarization ratio (LDR) at traditional cloud radar frequencies at Ka and W bands (~35 and ~94 GHz, correspondingly). Low radar beam elevation (~5°) measurements were performed at Oliktok Point, Alaska, with a scanning fully polarimetric radar operating in the horizontal–vertical polarization basis. This radar has the same gate spacing and very close beam widths at both frequencies, which largely alleviates uncertainties associated with spatial and temporal data matching. It is shown that observed Ka- and W-band Z DR differences are, on average, less than about 0.5 dB and do not have a pronounced trend as a function of snowfall reflectivity. The observed Z DR differences agree well with modeling results obtained using integration over nonspherical ice particle size distributions. For higher signal-to-noise ratios, KDP data derived from differential phase measurements are approximately scaled as reciprocals of corresponding radar frequencies indicating that the influence of non-Rayleigh scattering effects on this variable is rather limited. This result is also in satisfactory agreement with data obtained by modeling using realistic particle size distributions. Observed Ka- and W-band LDR differences are strongly affected by the radar hardware system polarization “leak” and are generally less than 4 dB. Smaller differences are observed for higher depolarizations, where the polarization “leak” is less pronounced. Realistic assumptions about particle canting and the system polarization isolation lead to modeling results that satisfactorily agree with observational dual-frequency LDR data.

42 ENGINEERING↗
DOE OSTI2020

A Computationally Efficient Linear Semi-Lagrangian Scheme for the Advection of Microphysical Variables in Cloud-Resolving Models

An efficient semi-Lagrangian (SL) numerical scheme for the advection of positive-definite microphysical variables in cloud-resolving models is proposed. The scheme is a linear combination of SL schemes of the first and second order. The proposed scheme is compared with two high-order, monotonic, nonlinear Eulerian schemes, in two 2D tests. In the first test, advection of a positive scalar field with large gradients is performed within complicated idealized steering flows. In the second test, a hail storm is simulated using the 2D Hebrew University Cloud Model, describing the interactions between 15 size distributions of different microphysical quantities (each distribution is described by 43 mass bins). The proposed scheme is computationally efficient, has a low numerical diffusion and a high level of mass conservation accuracy, and preserves the sum of multiple advected variables as well as the linear relationships between the advected scalar variables. Although the proposed SL scheme is not exactly positive definite, the negative values and their number, which may appear only in the case of very strong gradients, are negligible. The scheme produces results similar to those of nonlinear Eulerian advection schemes while reducing computation time by approximately a factor of 10.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Microbial extracellular enzyme activity with simulated climate change

It is critical to understand the consequences of environmental change for the microbial regulation of carbon and nutrient cycling. Specifically, understanding microbial community traits, such as extracellular enzyme activity, can help inform nutrient cycling models and address knowledge gaps. We analyzed data on extracellular enzyme activities and litter decomposition from an 18-month experiment in which microbial communities were reciprocally transplanted along a climate gradient in Southern California. Communities were from desert, scrubland, grassland, pine–oak, and subalpine ecosystems. We aimed to test how enzyme activities responded to climate change following transplantation and how those responses related to decomposition rates. We hypothesized that microbial communities would specialize on their native climate conditions, resulting in higher enzyme activities when transplanted back into their native climate. We investigated the relationship between extracellular enzyme Vmax values, substrate mass loss, and microbial biomass as well as variation in these variables across the climate gradient. We found little evidence for climate specialization, and there was rarely a reduction in enzyme functioning after microbial communities were transplanted into new climate conditions. Moreover, observed differences in decomposition were not related to changes in extracellular enzyme potential, although there were significant differences in enzyme activities and decomposition rates across sites. These results suggest that direct, physiological impacts of climate are likely to be important for enzyme-mediated decomposition, but climate specialization will not constrain the microbial response to climate change in our system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Data and scripts associated with “Allometric scaling of hyporheic respiration across basins in the Pacific Northwest USA"

This data package is associated with the publication “Allometric scaling of hyporheic respiration across basins in the Pacific Northwest USA” submitted to JGR-Biogeosciences (Regier et al. 2025).This study used reach-scale modeled estimates of hyporheic aerobic respiration made by the River Corridor Model (Fang et al. 2020) and watershed characteristics across the Willamette and Yakima River basins to explore potential allometric scaling (i.e., power-law relationships between size and function) of cumulative hyporheic respiration across catchment-to-basin scales. Scaling was explored quantitatively via the R2, slope, and y-intercept of relationships between cumulative hyporheic respiration and watershed area, divided into hyporheic exchange flux (HEF) quantiles. We also explored relationships between allometric scaling and other watershed characteristics through linear regression, spatial patterns, and mutual information analyses. Our results also suggest variability of hyporheic respiration allometry for middle exchange flux quantiles, and in relation to land-cover. Our findings provide initial evidence that allometric scaling may be useful for predicting hyporheic biogeochemical dynamics across watersheds from reach to basin scales. This data package is associated with the GitHub repository found at https://github.com/peterregier/rc_wrb_yrb_scaling. The data package is organized into several key directories. The “data” folder contains multiple CSV files, including landscape heterogeneity, scaling analysis, and watershed boundary data. The “figures” folder has all figure files in both PDF and PNG formats. Core analysis scripts and figure generation scripts are in the “scripts” directory, systematically numbered for sequential execution. The root directory includes essential project files; please see the file ending in “flmd.csv” for a list and description of all files contained in this data package and the file ending in “dd.csv” for data dictionaries used to describe tabular column headers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Hot Droughts and Forest Tree Dynamics in the Amazon - Statistical Models, Scripts, Data, and Outputs

This package contains data, outputs, equations, and R scripts for analyses for manuscript entitled "Hot droughts in the Amazon: A window to a future hypertropical climate" by J. Chambers et al., in particular it contains statistical models and analyses for the INPA BIONTE tree mortality study. The Models folder contains details for all statistical models in PDF files. The Scripts folder contains the R scripts for Bayesian Hierarchical Models (two text files) and SEMs (one text file) are separate and reasonably annotated. All data associated with these scripts are in the data folder. The Data folder contains two of the three CSV files used for the analyses and are called by the R scripts. Two of them are part of published datasets (`BIONTE_mortality-rates.csv` from Lima et al. 2024, DOI:10.15486/ngt/1898910 and `SPEI.csv` from Pastorello et al. 2023 DOI:10.15486/ngt/1958257) and also provided in this package for convenience (please see the corresponding datasets for usage and citation terms). The third dataset (`BIONTE_gapfilled_wd.csv`) contains sensitive information and can be obtained by contacting the manuscript lead author. The Outputs folder contains the two output files that provide extra information about the analyses. The file `figuresFeb2025d.pdf` contains all the figures from the manuscript - captions are in the manuscript. The file `ChambersMS.pdf` contains primary results from Bayesian statistical models, regression analyses, and validation steps applied to the tree mortality data from the INPA experiments. The document includes visual summaries, model diagnostics, and leave-one-out (LOO) validation results. A breakdown of file contents can be found in the README file that is part of this package.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Data and scripts associated with “Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics within Water-Tagging Enabled Hydrologic Models” (v3)

This data package is associated with the publication “Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics within Water-Tagging Enabled Hydrologic Models” submitted to Journal of Advances in Modeling Earth Systems (Butler et al. 2025). This study developed the Sequential Precipitation Input Tagging (SPIT) framework to tag input precipitation and estimate water transit times and hydrologic tracers. SPIT tags all precipitation events at regular intervals over an extended period (monthly tags over seven years) in a hydrologic model from 2016-2022. SPIT is applied at six National Ecological Observatory Network (NEON) sites across the continental United States to calculate transit time distributions (TTD) and derive from these mean transit times (MTT), fractions of young water (Fyw), and hydrologic tracer concentrations in stream water (δ18O) within a water-tagging enabled version of the Weather Research and Forecast (WT-WRF-Hydro) model with national water model (NWM) configurations. We go on to validate WT-WRF-Hydro estimates against Butler et al. (2023), who analyzed the same NEON sites using stable water isotope data to estimate water transit times. This new tracking method provides a detailed picture of water movement and helps improve predictions about water availability in the future. This data package was originally published in January 2025. It was updated May 2025 (v2; new and modified files) and October 2025 (v3; new and modified files). File and folder names were not revised to indicate changes. See the change history section in the readme for more details. This data package contains the data and scripts used to develop the SPIT framework WT-WRF-Hydro (Water Tagging Weather Research and Forecasting Hydrologic) model and is associated with the following GitHub repository: https://github.com/zbutler33/SPIT-Framework. This data package contains five parent folders: (1) “Manipulated_outputs”, (2) “Metadata”, (3) “Observed”, (4) “Outputs”, and (5) “Scripts”. Each of these parent folders contains additional subfolders and files. Please see the FLMD (“v*_Butler_2024_WT_WRF_Hydro_flmd.csv”) for a list of all the files contained in this data package and descriptions for each. See the data dictionary (“v*_Butler_2024_WT_WRF_Hydro_dd.csv”) for definitions and units of all of the tabular (files ending in “.csv” and ".tsv") column headers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Probing eddy size and its effective mixing length in stably stratified roughness sublayer flows

Stably stratified roughness sublayer flows are ubiquitous yet remain difficult to represent in models and to interpret using field experiments. Here, continuous high-frequency potential temperature profiles from the forest floor up to 6.5 times the canopy height observed with distributed temperature sensing (DTS) are used to link eddy topology to roughness sublayer stability correction functions and coupling between air layers within and above the canopy. The experiments are conducted at two forest stands classified as hydrodynamically sparse and dense. Near-continuous profiles of eddy sizes (length scales) and effective mixing lengths for heat are derived from the observed profiles using a novel conditional sampling approach. The approach utilizes potential temperature isoline fluctuations from a statically stable background state. The transport of potential temperature by an observed eddy is assumed to be conserved (adiabatic movement) and we assume that irreversible heat exchange between the eddy and the surrounding background occurs along the (vertical) periphery of the eddy. This assumption is analogous to Prandtl's mixing-length concept, where momentum is transported rapidly vertically and then equilibrated with the local mean velocity gradient. A distinct dependence of the derived length scales on background stratification, height above ground, and canopy characteristics emerges from the observed profiles. Implications of these findings for (1) the failure of Monin–Obukhov similarity in the roughness sublayer and (2) above-canopy flow coupling to the forest floor are examined. The findings have practical applications in terms of analysing similar DTS data sets with the proposed approach, modelling roughness sublayer flows, and interpreting nocturnal eddy covariance measurements above tall forested canopies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Grassland carbon-water cycling is minimally impacted by a photovoltaic array

Agrivoltaic systems, whereby photovoltaic arrays are co-located with crop or forage production, can alleviate the tension between expanding solar development and loss of agricultural land. However, the ecological ramifications of these arrays are poorly known. We used field measurements and a plant hydraulic model to quantify carbon-water cycling in a semi-arid C 3 perennial grassland growing beneath a single-axis tracking solar array in Colorado, USA. Although the agrivoltaic array reduced light availability by 38%, net photosynthesis and aboveground net primary productivity were reduced by only 6–7% while evapotranspiration decreased by 1.3%. The minimal changes in carbon-water cycling occurred largely because plant photosynthetic traits underneath the panels changed to take advantage of the dynamic shading environment. Our results indicate that agrivoltaic systems can serve as a scalable way to expand solar energy production while maintaining ecosystem function in managed grasslands, especially in climates where water is more limiting than light.

14 SOLAR ENERGY↗
DOE OSTI2023

Linking repeat lidar with Landsat products for large scale quantification of fire-induced permafrost thaw settlement in interior Alaska

The permafrost–fire–climate system has been a hotspot in research for decades under a warming climate scenario. Surface vegetation plays a dominant role in protecting permafrost from summer warmth, thus, any alteration of vegetation structure, particularly following severe wildfires, can cause dramatic top–down thaw. A challenge in understanding this is to quantify fire-induced thaw settlement at large scales (>1000 km 2 ). In this study, we explored the potential of using Landsat products for a large-scale estimation of fire-induced thaw settlement across a well-studied area representative of ice-rich lowland permafrost in interior Alaska. Six large fires have affected ~1250 km 2 of the area since 2000. We first identified the linkage of fires, burn severity, and land cover response, and then developed an object-based machine learning ensemble approach to estimate fire-induced thaw settlement by relating airborne repeat lidar data to Landsat products. The model delineated thaw settlement patterns across the six fire scars and explained ~65% of the variance in lidar-detected elevation change. Our results indicate a combined application of airborne repeat lidar and Landsat products is a valuable tool for large scale quantification of fire-induced thaw settlement.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Increased salinity decreases annual gross primary productivity at a Northern California brackish tidal marsh

Tidal marshes sequester 11.4–87.0 Tg C yr –1 globally, but climate change impacts can threaten the carbon capture potential of these ecosystems. Tidal marshes occur across a wide range of salinity, with brackish marshes (0.5–18 ppt (parts per thousand)) dominating global tidal marsh extents. A diverse mix of freshwater- and saltwater-tolerant plant and microbial communities has led researchers to predict that carbon cycling in brackish wetlands may be less sensitive to changes in salinity than fresh- or saltwater wetlands. Rush Ranch, a well-monitored brackish tidal wetland of the San Francisco Bay National Estuarine Research Reserve, experiences highly variable annual salinity regimes. Within a five-year period (2014–2018), Rush Ranch experienced particularly extreme drought-induced salinization during the 2014 and 2015 growing seasons. During drought years, tidal channel salinity rose from a 15 year baseline of 4.7 ppt to growing season peaks of 10.3 ppt and 12.5 ppt. Continuous eddy covariance data from 2014 to 2018 demonstrate that during drought summers, gross primary productivity (GPP) decreased by 24%, whereas ecosystem respiration remained similar among all five years. Stepwise linear regression revealed that salinity, not air temperature or tidal height, was the dominant driver of annual GPP. A random forest model trained to predict GPP based on environmental data from low salinity years (i.e. naive to salinization) significantly over predicted GPP in drought years. When growing season salinities were doubled, annual estimates of net ecosystem exchange of CO 2 decreased by up to 30%. These results provide ecosystem-scale evidence that increased salinity influences CO 2 fluxes dominantly through reductions in GPP. This relationship provides a starting point for incorporating the effect of changes in salinity in wetland carbon models, which could improve wetland carbon forecasting and management for climate resilience.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Dynamics of Subsiding Shells in Actively Growing Clouds with Vertical Updrafts

The dynamics of a subsiding shell at the edges of actively growing shallow cumulus clouds with updrafts is analyzed using direct numerical simulation. The actively growing clouds have a fixed in-cloud buoyancy and velocity. Turbulent mixing and evaporative cooling at the cloud edges generate a subsiding shell that grows with time. A self-similar regime is observed for first- and second-order moments when normalized with respective maximum values. Internal scales derived from integral properties of the flow problem are identified. A self-similarity analysis using these scales reveals that contrary to classical self-similar flows, the turbulent kinetic energy budget terms and velocity moments scale according to the buoyancy and not with the mean velocity. The shell thickness is observed to increase linearly with time. The buoyancy scale remains time invariant and is set by the initial cloud–environment thermodynamics. The shell accelerates ballistically with a magnitude set by the saturation value of the buoyancy of the cloud–environment mixture. In this regime, the shell is buoyancy driven and independent of the in-cloud velocity. Relations are obtained for predicting the shell thickness and minimum velocities by linking the internal scales with external flow parameters. The values thus calculated are consistent with the thickness and velocities observed in typical shallow cumulus clouds. The entrainment coefficient is a function of the initial state of the cloud and the environment, and is shown to be on the same order of magnitude as fractional entrainment rates calculated for large-scale models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Analytical Investigation of the Role of Lateral Mixing in the Evolution of Nonprecipitating Cumulus. Part I: Developing Clouds

Evolution of nonprecipitating cumulus clouds (Cu) at the developing stage under the influence of lateral entrainment and mixing is studied analytically using a minimalistic analytical model. We present a model of an ascending cloud volume (a model of developing Cu) whose structure is determined by the processes of droplet diffusion growth/evaporation and entrainment mixing in the horizontal direction. Spatial and time changes of liquid water content, the adiabatic fraction, droplet concentration, and the mean volume droplet radius are calculated. It is shown that the existence of a nondiluted core in a growing cumulus cloud significantly depends on the cloud width and vertical velocity. While at the updraft velocity of 2 m s −1 the core of a 400-m-wide cloud becomes diluted at distances of a few hundred meters above cloud base, the core of a cloud of 1000-m width remains nondiluted at distances up to 1500 m above cloud base. The explanation of this result is simple: the increase in cloud width and the decrease in the updraft velocity increase the time during which the cloud is diluted due to mixing. Since lateral mixing synchronously decreases both the cloud water content and droplet concentration, the variation of the mean volume droplet radius is low inside the cloud. The approximate quantitative condition for cloud formation in updraft is derived. It is shown that a cloud can arise when its vertical velocity exceeds a critical value. To produce clouds, narrow turbulent plumes should ascend at higher velocity as compared to wider plumes. High humidity of the environment air is favorable for formation of clouds from plumes. The comparison of the obtained results with previously published observational data indicates a reasonable agreement. The results can be useful for parameterization purposes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Toward a Theory of the Evolution of Drop Morphology and Splintering by Freezing

Abstract The drop freezing process is described by a phase-field model. Two cases are considered: when the freezing is triggered by central nucleation and when nucleation occurs on the drop surface. Depending on the environmental temperature and drop size, different morphological structures develop. Detailed dendritic growth was simulated at the first stage of drop freezing. Independent of the nucleation location, a decrease in temperature within the range from ~−5° to −25°C led to an increase in the number of dendrites and a decrease in their width and the interdendritic space. At temperatures lower than about −25°C, a planar front developed following surface nucleation, while dendrites formed a granular-like structure with small interdendritic distances following bulk nucleation. An ice shell grew in at the same time (but slower) as dendrites following surface nucleation, while it started forming once the dendrites have reached the drop surface in the case of central nucleation. The formed ice morphology at the first freezing stage predefined the splintering probability. We assume that stresses needed to break the ice shell arose from freezing of the water in the interdendritic spaces. Under this assumption, the number of possible splinters/fragments was proportional to the number of dendrites, and the maximum rate of splintering/fragmentation occurred within a temperature range of about −10° to −20°C, is in agreement with available laboratory and in situ measurements. At temperatures < −25°C, freezing did not lead to the formation of significant stresses, making splintering unlikely. The number of dendrites increased with drop size, causing a corresponding increase in the number of splinters. Examples of morphology that favors drop cracking are presented, and the duration of the freezing stages is evaluated. Sensitivity of the freezing process to the surface fluxes is discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Effects of Vegetation and Topography on the Boundary Layer Structure above the Amazon Forest

Observational data from two field campaigns in the Amazon forest were used to study the vertical structure of turbulence above the forest. The analysis was performed using the reduced turbulent kinetic energy (TKE) budget and its associated two-dimensional phase space. Results revealed the existence of two regions within the roughness sublayer in which the TKE budget cannot be explained by the canonical flat-terrain TKE budgets in the canopy roughness sublayer or in the lower portion of the convective ABL. Data analysis also suggested that deviations from horizontal homogeneity have a large contribution to the TKE budget. Results from LES of a model canopy over idealized topography presented similar features, leading to the conclusion that flow distortions caused by topography are responsible for the observed features in the TKE budget. These results support the conclusion that the boundary layer above the Amazon forest is strongly impacted by the gentle topography underneath.

54 ENVIRONMENTAL SCIENCES↗