Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HYPOTHESIS”

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 505 records · Page 28

Biological mechanisms may contribute to soil carbon saturation patterns

Increasing soil organic carbon (SOC) storage is a key strategy to mitigate rising atmospheric CO 2 , yet SOC pools often appear to saturate, or increase at a declining rate, as carbon (C) inputs increase. Soil C saturation is commonly hypothesized to result from the finite amount of reactive mineral surface area available for retaining SOC, and is accordingly represented in SOC models as a physicochemically determined SOC upper limit. However, mineral-associated SOC is largely microbially generated. In this perspective, we present the hypothesis that apparent SOC saturation patterns could emerge as a result of ecological constraints on microbial biomass—for example, via competition or predation—leading to reduced C flow through microbes and a reduced rate of mineral-associated SOC formation as soil C inputs increase. Microbially explicit SOC models offer an opportunity to explore this hypothesis, yet most of these models predict linear microbial biomass increases with C inputs and insensitivity of SOC to input rates. Synthesis of 54 C addition studies revealed constraints on microbial biomass as C inputs increase. Different hypotheses limiting microbial density were embedded in a three-pool SOC model without explicit limits on mineral surface area. As inputs increased, the model demonstrated either no change, linear, or apparently saturating increases in mineral-associated and particulate SOC pools. Taken together, our results suggest that microbial constraints are common and could lead to reduced mineral-associated SOC formation as input rates increase. We conclude that SOC responses to altered C inputs—or any environmental change—are influenced by the ecological factors that limit microbial populations, allowing for a wider range of potential SOC responses to stimuli. Understanding how biotic versus abiotic factors contribute to these patterns will better enable us to predict and manage soil C dynamics.

54 ENVIRONMENTAL SCIENCES↗

Simulating evaporative wet and dry cycles in Gale crater, Mars using thermochemical modeling techniques

The aim of this work is to provide a model‐backed hypothesis for the formation of evaporites—sulfates, borates—in Gale crater using thermochemical modeling to determine constraints on their formation. We test the hypothesis that primary evaporites required multiple wet–dry cycles to form, akin to how evaporite assemblages form on Earth. Starting with a basalt‐equilibrated Mars fluid, Mars‐relevant concentrations of B and Li were added, and then equilibrated with Gale lacustrine bedrock. We simulated the cycles of evaporation followed by groundwater recharge/dilution to establish an approximate minimum number of wet–dry cycles required to form primary evaporites. We determine that a minimum of 250 wet–dry cycles may be required to start forming primary evaporites that consist of borates and Ca‐sulfates. We estimate that ~14,250 annual cycles (~25.6 k Earth years) of wet and dry periods may form primary borates and Ca‐sulfates in Gale crater. These primary evaporites could have been remobilized during secondary diagenesis to form the veins that the Curiosity rover observes in Gale crater. No Li salts form after 14,250 cycles modeled for the Gale‐relevant scenario (approximately 10 6 cycles would be needed) which implies Li may be leftover in a groundwater brine after the time of the lake. No major deposits of borates are observed to date in Gale crater which also implies that B may be leftover in the subsequent groundwater brine that formed after evaporites were remobilized into Ca‐sulfate veins.

58 GEOSCIENCES↗

Worldwide phylogeography of rough-toothed dolphins ( Steno bredanensis ) provides evidence for subspecies delimitation

Rough-toothed dolphins (Steno bredanensis) have a global tropical and subtropical distribution with oceanic, neritic, and island-associated populations. Here, to inform conservation and management for this species, we used sequences from the mtDNA control region (n = 360), mitogenomes (n = 19), and six nuclear introns (n = 35) to provide multiple lines of evidence to critically evaluate the potential taxonomic status of rough-toothed dolphins. Using samples from the Pacific, Indian, and Atlantic Oceans, we examined the null hypothesis that rough-toothed dolphins are one panmictic species and the alternate hypothesis of oceanic subspecies. Phylogenetic analyses of mitogenomes revealed a private Atlantic clade sister to a larger cosmopolitan clade including individuals from all tropical and subtropical oceans. We dated the split between the Atlantic clade and the cosmopolitan clade to 890,000 years ago. We determined that Atlantic rough-toothed dolphins could be correctly diagnosed with 98% accuracy with the mtDNA control region and calculated the net nucleotide divergence as 0.02. Population level analyses revealed significant genetic differentiation using mtDNA among most regions, while significant differentiation using nuclear markers occurred only between the Atlantic and the Indian/Pacific regions. Therefore, the oceanic divergence and diagnosability of rough-toothed dolphins in the Atlantic and the Indian/Pacific Oceans meet proposed criteria for recognition as two subspecies.

60 APPLIED LIFE SCIENCES↗

Deciphering the shifting role of intrinsic and extrinsic drivers on moss decomposition in peatlands over a 5-year period

Bryophytes are largely responsible for globally significant carbon accumulation in peatland ecosystems. This accumulation is primarily caused by the slow decomposition of these mosses, which can be attributed to a combination of intrinsic (chemical) characteristics of decaying mosses and extrinsic (environmental) influences. Here we investigated the importance of intrinsic and extrinsic drivers of moss decomposition in peatlands and hypothesized that the early stages of decay will be driven primarily by intrinsic characteristics, while extrinsic drivers will become more important in the later phases of decomposition. We tested this hypothesis by placing litterbags of three moss types (Sphagnum divinum, Sphagnum angustifolium/fallax and Polytrichum sp.) into hummock and hollow microtopographies in a bog and a lawn microtopography in a nearby poor fen in northern Minnesota, USA. Decomposition was measured over a 5-year period and confirmed our hypothesis; while intrinsic and extrinsic variables were both important in driving decay rates over time, decay in the first year was primarily predicted by litter type, while extrinsic influences (e.g. peat temperature, moisture, nutrient limitation) that varied at peatland and microtopographic scales became more important by year 5. Though the nutrient (i.e. nitrogen [N] and phosphorus [P]) content of decaying mosses did not explain the variation in decay rates among litter types, analysis of N:P in moss litter suggested that litter-decaying microorganisms became increasingly P-limited over time. The majority of litter decomposition studies in peatlands are generally of short duration (< 3 years). Our study suggests that longer-term studies are needed to reveal the roles of multiple drivers of moss decay. Improving the mechanistic understanding of decomposition processes in peatlands will be critical to developing more robust representations of peatlands and their ecosystem processes in carbon and climate models.

29 poor fen↗

Contrasting anther glucose–6–phosphate dehydrogenase activities between two bean varieties suggest an important role in reproductive heat tolerance

Common beans (Phaseolus vulgaris) are highly sensitive to elevated temperatures, and rising global temperatures threaten bean production. Plants at the reproductive stage are especially susceptible to heat stress due to damage to male (anthers) and female (ovary) reproductive tissues, with anthers being more sensitive to heat. Heat damage promotes early tapetal cell degradation, and in beans this was shown to cause male infertility. In this study, we focus on understanding how changes in leaf carbon export in response to elevated temperature stress contribute to heat-induced infertility. We hypothesize that anther glucose-6-phosphate dehydrogenase (G6PDH) activity plays an important role at elevated temperature and promotes thermotolerance. To test this hypothesis, we compared heat-tolerant and susceptible common bean genotypes using a combination of phenotypic, biochemical, and physiological approaches. Our results identified changes in leaf sucrose export, anther sugar accumulation and G6PDH activity and anther H 2 O 2 levels and antioxidant-related enzymes between genotypes at elevated temperature. Further, anther respiration rate was found to be lower at high temperature in both bean varieties. Overall, our results support the hypothesis that enhanced male reproductive heat tolerance involves changes in the anther oxidative pentose phosphate pathway, which supplies reductants to critical H 2 O 2 scavenging enzymes.

59 BASIC BIOLOGICAL SCIENCES↗

Arabidopsis stromal carbonic anhydrases exhibit non–overlapping roles in photosynthetic efficiency and development

Carbonic anhydrases (CAs) are ubiquitous enzymes that accelerate the reversible conversion of CO 2 to HCO 3 – . The Arabidopsis genome encodes members of the α-, β- and γ-CA families, and it has been hypothesized that βCA activity has a role in photosynthesis. In this work, we tested this hypothesis by characterizing the two plastidial βCAs, βCA1 and βCA5, in physiological conditions of growth. We conclusively established that both proteins are localized in the chloroplast stroma and that the loss of βCA5 induced the expression of βCA1, supporting the existence of regulatory mechanisms to control the expression of stromal βCAs. We also established that βCA1 and βCA5 have markedly different enzymatic kinetics and physiological relevance. Specifically, we found that βCA5 had a first-order rate constant ~10-fold lower than βCA1, and that the loss of βCA5 is detrimental to growth and could be rescued by high CO 2 . Furthermore, we established that, while a βCA1 mutation showed near wild-type growth and no significant impact on photosynthetic efficiency, the loss of βCA5 markedly disrupted photosynthetic efficiency and light-harvesting capacity at ambient CO 2 . Therefore, we conclude that in physiological autotrophic growth, the loss of the more highly expressed βCA1 does not compensate for the loss of a less active βCA5, which in turn is involved in growth and photosynthesis at ambient CO 2 levels. These results lend support to the hypothesis that, in Arabidopsis,βCAs have non-overlapping roles in photosynthesis and identify a critical activity of stromal βCA5 and a dispensable role for βCA1.

59 BASIC BIOLOGICAL SCIENCES↗

On the Cutting of Metals: A Mechanics Viewpoint

The mechanics of large-strain deformation in cutting of metals is discussed, primarily from viewpoint of recent developments in in situ analysis of plastic flow and microstructure characterization. It is shown that a broad range of deformation parameters can be accessed in chip formation—strains of 1–10, strain rates of 10–10 5 /s, and temperatures up to 0.7T m —and controlled. This range is far wider than achievable by any other single-stage, severe plastic deformation (SPD) process. The resulting extreme deformation conditions produce a rich variety of microstructures in the chip. Four principal types of chip formation—continuous, shear-localized, segmented, and mushroom-type—as elucidated first by Nakayama (1974, “The Formation of ‘Saw-Toothed Chip’ in Metal Cutting,” Proceedings of International Conference on Production Engineering, Tokyo, pp. 572–577) are utilized to emphasize the diverse plastic flow phenomena, especially unsteady deformation modes that prevail in cutting. These chip types are intimately connected with the underlying flow, each arising from a distinct mode and triggered by an instability phenomenon. The role of plastic flow instabilities such as shear banding, buckling, and fracture in mediating unsteady flow modes is expounded, along with consequences of the flow modes and chip types for the cutting. Sinuous flow is shown to be the reason why gummy (highly strain-hardening) metals, although relatively soft, are so difficult to cut. Synthesizing the various observations, a hypothesis is put forth that it is the stability of flow modes that determines the mechanics of cutting. Furthermore, this leads to a flow-stability phase diagram that could provide a framework for predicting chip types and process attributes.

42 ENGINEERING↗

Describing the Mechanism of Instability Suppression Using a Central Pilot Flame With Coupled Experiments and Simulations

Pilot flames are commonly used to extend combustor operability limits and suppress combustion oscillations in low-emissions gas turbines. Combustion oscillations, a coupling between heat release rate oscillations and combustor acoustics, can arise at the operability limits of low-emissions combustors where the flame is more susceptible to perturbations. While the use of pilot flames is common in land-based gas turbine combustors, the mechanism by which they suppress instability is still unclear. In this study, we consider the impact of a central jet pilot on the stability of a swirl-stabilized flame in a variable-length, single-nozzle combustor. Previously, the pilot flame was found to suppress the instability for a range of equivalence ratios and combustor lengths. We hypothesize that combustion oscillation suppression by the pilot occurs because the pilot provides hot gases to the vortex breakdown region of the flow that recirculate and improve the static, and hence dynamic, stability of the main flame. This hypothesis is based on a series of experimental results that show that pilot efficacy is a strong function of pilot equivalence ratio but not pilot flow rate, which would indicate that the temperature of the pilot products as well as the combustion intensity of the pilot flame play more of a role in oscillation stabilization than the length of the pilot flame relative to the main flame. Further, the pilot-flame efficacy increases with pilot-flame equivalence ratio until it matches the main-flame equivalence ratio; at pilot equivalence ratios greater than the main equivalence ratio, the pilot-flame efficacy does not change significantly with pilot equivalence ratio. To understand these results, we use large-eddy simulation (LES) to provide a detailed analysis of the flow in the region of the pilot flame and the transport of radical species in the region between the main flame and pilot flame. Furthermore, the simulation, using a flamelet/progress variable-based chemistry tabulation approach and standard eddy viscosity/diffusivity turbulence closure models, provides detailed information that is inaccessible through experimental measurements.

42 ENGINEERING↗

Cryptic cycling by electroactive bacterioplankton in Trout Bog Lake

The potential for extracellular electron transfer (EET) is a prevailing genomic feature of humic lake bacterioplankton. However, there has been little evidence for the substantial ecological contribution predicted by genetics. We hypothesized that anoxygenic phototrophic electrotrophs and accompanying heterotrophic electrogens cycle dissolved organic matter (DOM) between oxidized and reduced states. We predicted that such bacterioplankton would exhibit diel-scale oscillations due to the light dependency of photosynthesis. Using Trout Bog Lake in Wisconsin, USA, as our model ecosystem, we profiled the water column with depth-discrete metagenomic, physiochemical, and electrochemical analyses. We observed variation in oxidation reduction potential (ORP) in response to sunlight, initiating at depths populated by anoxygenic phototrophs with EET genes. We developed an automated buoy to measure electric current flow between many pairs of electrodes simultaneously, observing correlation in electron consumption to sunlight. Our results, combined with published metatranscriptomic analysis, indicate the occurrence of electron cycling between phototrophic oxidation (electrotrophic metabolism) by Chlorobium and anaerobic respiration (electrogenic metabolism) by Geothrix, involving DOM. We also repeatedly observed gradual seasonal increases in hypolimnion ORP throughout summer. These diel and seasonal patterns imply that electroactive DOM mediates the ecology of electroactive bacteria in lakes, controlling humic lake methane emissions.IMPORTANCEWe investigated the physical, chemical, and redox characteristics of a bog lake and electrodes hung therein to test the hypothesis that dissolved organic matter is being cycled between oxidized and reduced states by electroactive bacterioplankton powered by phototrophy. To do so, we performed field-based analyses on multiple timescales using both established and novel instrumentation. We paired these analyses with recently developed bioinformatics pipelines for metagenomics data to investigate genes that enable electroactive metabolism and accompanying metabolisms. Our results are consistent with our hypothesis and yet upend some of our other expectations. Our findings have implications for understanding greenhouse gas emissions from lakes, including electroactivity as an integral part of lake metabolism throughout more of the anoxic parts of lakes and for a longer portion of the summer than expected. Our results also give a sense of what electroactivity occurs at given depths and provide a strong basis for future studies.

carbon emissions↗

Observation of new excited ${B} ^0_{s}$ states

A structure is observed in the ${B} ^{\pm }{K} ^{\mp }$ mass spectrum in a sample of proton–proton collisions at centre-of-mass energies of 7, 8, and 13 TeV, collected with the LHCb detector and corresponding to a total integrated luminosity of 9$\,\text {fb} ^{-1}$. The structure is interpreted as the result of overlapping excited ${B} ^0_{s}$ states. With high significance, a two-peak hypothesis provides a better description of the data than a single resonance. Under this hypothesis the masses and widths of the two states, assuming they decay directly to ${B} ^{\pm }{K} ^{\mp }$, are determined to be $\begin{aligned} m_1&= 6063.5 \pm 1.2 \text { (stat)} \pm 0.8\text { (syst)}\,\text {Me}\text {V}, \\ \Gamma _1&= 26 \pm 4 \text { (stat)} \pm 4\text { (syst)}\,\text {Me}\text {V}, \\ m_2&= 6114 \pm 3 \text { (stat)} \pm 5\text { (syst)}\,\text {Me}\text {V}, \\ \Gamma _2&= 66 \pm 18 \text { (stat)} \pm 21\text { (syst)}\,\text {Me}\text {V}. \end{aligned}$ Alternative values assuming a decay through ${B} ^{*\pm }{K} ^{\mp }$, with a missing photon from the $B^{*\pm } \rightarrow B^{\pm }\gamma $ decay, which are shifted by approximately 45 $\,\text {Me}$V, are also determined. The possibility of a single state decaying in both channels is also considered. The ratio of the total production cross-section times branching fraction of the new states relative to the previously observed $B_{s2}^{*0}$ state is determined to be $0.87 \pm 0.15 \text { (stat)} \pm 0.19 \text { (syst)}$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Measurement of the CP properties of Higgs boson interactions with $\tau$-leptons with the ATLAS detector

A study of the charge conjugation and parity ($CP$) properties of the interaction between the Higgs boson and $\tau$-leptons is presented. The study is based on a measurement of $CP$-sensitive angular observables defined by the visible decay products of $\tau$-leptons produced in Higgs boson decays. The analysis uses 139 fb -1 of proton–proton collision data recorded at a centre-of-mass energy of $\sqrt{s}$ = 13 TeV with the ATLAS detector at the Large Hadron Collider. Contributions from $\tau$-violating interactions between the Higgs boson and $\tau$-leptons are described by a single mixing angle parameter Φ $\tau$ in the generalised Yukawa interaction. Without constraining the $H$ $\rightarrow$ $\tau$$\tau$ signal strength to its expected value under the Standard Model hypothesis, the mixing angle Φ $\tau$ is measured to be 9° ± 16°, with an expected value of 0° ± 28° at the 68% confidence level. The pure $CP$-odd hypothesis is disfavoured at a level of 3.4 standard deviations. The results are compatible with the predictions for the Higgs boson in the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Search for baryon junctions in e+A collisions at the electron ion collider

Constituent quarks in a nucleon are the essential elements in the standard “quark model” associated with the electric charge, spin, mass, and baryon number of a nucleon. Quantum chromodynamics (QCD) describes nucleon as a composite object containing current quarks (valence quarks and sea (anti-)quarks) and gluons. These subatomic elements and their interactions are known to contribute in complex ways to the overall nucleon spin and mass. In the early development of QCD theory in the 1970s, an alternative hypothesis postulated that the baryon number might manifest itself through a non-perturbative configuration of gluon fields forming a Y-shaped topology known as the gluon junction. In this work, we propose to test such hypothesis by measuring (i) the Regge intercept of the net-baryon distributions for e+(p)Au collisions, (ii) baryon and charge transport in the isobaric ratio between e+Ru and e+Zr collisions, and (iii) target flavor dependence of proton and antiproton yields at large rapidity, transported from the hydrogen and deuterium targets in e+p(d) collisions. Our study indicates that these measurements at the EIC can help determine what carries the baryon number.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Ensuring statistical reproducibility of ocean model simulations in the age of hybrid computing

Novel high performance computing systems that feature hybrid architectures require large scale code refactoring to unravel underlying exploitable parallelism. Such redesign can often be accompanied with machine-precision changes as the order of computation cannot always be maintained. For chaotic systems like climate models, these round-off level differences can grow rapidly. Systematic errors may also manifest initially as machine-precision differences. Isolating genuine round off level differences from such errors remains a challenge. Here, we apply two-sample equality of distribution tests to evaluate statistical reproducibility of the ocean model component of US Department of Energy's Energy Exascale Earth System Model (E3SM). A 2-year control simulation ensemble is compared to a modified ensemble as a test case - after a known non-bit-for-bit change in a model component is introduced - to evaluate the null hypothesis that the two ensembles are statistically indistinguishable. To quantify the false negative rates of these tests, we conduct a formal power analysis using a targeted suite of short simulation ensembles. The ensemble suite contains several perturbed ensembles, each with a progressively different climate than the baseline ensemble - obtained by perturbing the magnitude of a single model tuning parameter, the Gent and McWilliams κ, in a controlled manner. The null hypothesis is evaluated for each of perturbed ensembles using these tests. The power analysis informs on the detection limits of the tests for given ensemble size allowing model developers to evaluate the impact of an introduced non-bit-for-bit change to the model.

Mahajan, Salil↗

matsim-agents v1.0

matsim-agents is a multi-agent AI framework for atomistic materials simulation and discovery. It orchestrates large language models (LLMs), machine-learned interatomic potentials (MLIPs), and DFT codes into a single agentic loop running on laptops and DOE leadership-class supercomputers. MULTI-AGENT ORCHESTRATION A LangGraph state machine with three nodes: a Planner that converts a natural-language research objective into structured tasks; an Executor that dispatches atomistic tools and loops until the queue is empty; and an Analyst that summarizes results into a human-readable report. State is checkpointed after every step and human-in-the-loop gates can be inserted at any edge. HYPOTHESIS-DRIVEN DISCOVERY CHAT An interactive REPL (matsim-agents chat) that couples LLM dialogue with atomistic simulation. Chemical formulas are automatically detected in conversation turns and trigger a full crystal-phase exploration: structure generation → relaxation → stability scoring → result injection back into the conversation, creating a closed hypothesis-refinement loop. CRYSTAL PHASE ENUMERATION Given a composition, the phase explorer enumerates prototypes by stoichiometry: elemental (fcc/bcc/hcp/sc/diamond), binary 1:1 (rocksalt/CsCl/zincblende/ wurtzite/fluorite/rutile), ternary 1:1:3 (cubic perovskite), ternary 1:2:4 (perovskite + spinel), quaternary 1:1:2:6 (Fm-3m double perovskite). 2-D prototypes (graphene, h-BN, MoS2 2H/1T) and multilayer stacking are also supported via --include-2d and --num-layers. SUPERCELL GENERATION AND SITE DECORATION Auto-tiling to a minimum atom count (--min-atoms), explicit NxNxN tiling (--supercell), symmetry-distinct site decorations (--n-orderings), and isotropic lattice-scale sweeps (--lattice-scales) for volume bracketing. MLFF RELAXATION AND STABILITY SCORING HydraGNN (multi-headed GNN) drives structure relaxation via ASE with FIRE, BFGS, or BFGSLineSearch. Stability output: delta-E/atom ranking across phases and a max-residual-force dynamical-stability proxy. Other MLIPs (MACE, NequIP, Orb) can be plugged in through the same interface. DFT BACKENDS Quantum ESPRESSO pw.x and VASP 6.6 are first-class labellers. Both have validated GPU builds and SLURM/PBS launchers for three DOE platforms: Frontier (AMD MI250X, ROCm), Aurora (Intel PVC, oneAPI), Perlmutter (NVIDIA A100, CUDA). QE produces ~100 binaries (pw.x, ph.x, epw.x, ...). VASP supports scf, relax, vc-relax, and vc-relax-shape run types. ACTIVE-LEARNING LOOP matsim-agents al run CONFIG.yaml drives an iterative HydraGNN-DFT loop: MD generates candidates → ensemble/MC-dropout uncertainty selects the most informative → DFT labels them in parallel inside one allocation → dataset grows → HydraGNN retrains → repeat. DFT backend is a single YAML toggle (dft.backend: vasp | qe). LLM-generated seed structures are supported (no curated POSCAR library needed). Config uses ${VAR}, ${VAR:-default}, ${VAR:?msg} shell-style substitution for cross-user/cross-site portability. LLM BACKENDS Ollama (local, default), vLLM (HPC multi-GPU serving), OpenAI, Anthropic, HuggingFace Transformers+Accelerate. Selected at runtime via flag or env var with no code changes. HPC PORTABILITY Same Python entry points run on Frontier (ROCm 7.2), Aurora (oneAPI), and Perlmutter (CUDA 12). DFT and ML stacks are never co-loaded in the same shell; they couple through the scheduler and filesystem. Advanced multi-node launchers (serve, discovery-chat, single-relaxation, active-learning, QE warm-start) are provided for all three platforms. CODABENCH COMPETITION BUNDLE A self-contained benchmark: 159 atomistic test structures across 11 material classes, 5 tasks (formation energy, forces, ML relaxation, AI-DFT relaxation, phase stability ranking), public/private leaderboard split (30/70), and four ready-to-run baselines: MACE-MP-0, HydraGNN, UMA, AllScAIP.

Lupo Pasini, Massimiliano [Oak Ridge National Labo↗

Automated Identification of Characteristic Droplet Size Distributions in Stratocumulus Clouds Utilizing a Data Clustering Algorithm

Abstract Droplet-level interactions in clouds are often parameterized by a modified gamma fitted to a “global” droplet size distribution. Do “local” droplet size distributions of relevance to microphysical processes look like these average distributions? This paper describes an algorithm to search and classify characteristic size distributions within a cloud. The approach combines hypothesis testing, specifically, the Kolmogorov–Smirnov (KS) test, and a widely used class of machine learning algorithms for identifying clusters of samples with similar properties: density-based spatial clustering of applications with noise (DBSCAN) is used as the specific example for illustration. The two-sample KS test does not presume any specific distribution, is parameter free, and avoids biases from binning. Importantly, the number of clusters is not an input parameter of the DBSCAN-type algorithms but is independently determined in an unsupervised fashion. As implemented, it works on an abstract space from the KS test results, and hence spatial correlation is not required for a cluster. The method is explored using data obtained from the Holographic Detector for Clouds (HOLODEC) deployed during the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) field campaign. The algorithm identifies evidence of the existence of clusters of nearly identical local size distributions. It is found that cloud segments have as few as one and as many as seven characteristic size distributions. To validate the algorithm’s robustness, it is tested on a synthetic dataset and successfully identifies the predefined distributions at plausible noise levels. The algorithm is general and is expected to be useful in other applications, such as remote sensing of cloud and rain properties. Significance Statement A typical cloud can have billions of drops spread over tens or hundreds of kilometers in space. Keeping track of the sizes, positions, and interactions of all of these droplets is impractical, and, as such, information about the relative abundance of large and small drops is typically quantified with a “size distribution.” Droplets in a cloud interact locally, however, so this work is motivated by the question of whether the cloud droplet size distribution is different in different parts of a cloud. A new method, based on hypothesis testing and machine learning, determines how many different size distributions are contained in a given cloud. This is important because the size distribution describes processes such as cloud droplet growth and light transmission through clouds.

54 ENVIRONMENTAL SCIENCES↗

A Formula for the Maximum Vertical Velocity in Supercell Updrafts

In supercell environments, previous authors have shown strong connections between the vertical wind shear magnitude, updraft width, and entrainment. Based on these results, it is hypothesized that the influences of entrainment-driven dilution on buoyancy and maximum updraft vertical velocity w in supercell environments are a predictable function of the vertical wind shear profile. It is also hypothesized that the influences of pressure perturbation forces on maximum updraft w are small because of a nearly complete offset between upward dynamic pressure forces and downward buoyant pressure forces. To address these hypotheses, we derive a formula for the maximum updraft w that incorporates the effects of entrainment-driven dilution on buoyancy but neglects pressure gradient forces. Solutions to this formula are compared with output from previous numerical simulations. This formula substantially improves predictions of maximum updraft w over past CAPE-derived formulas for maximum updraft w, which supports the first hypothesis. Furthermore, integrated vertical accelerations along trajectories show substantial offsets between dynamic and buoyant pressure forces, supporting the second hypothesis. It is argued that the new formula should be used in addition to CAPE-derived measures for w in forecast and research applications when accurate diagnosis of updraft speed is required.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Influences of Environmental Relative Humidity and Horizontal Scale of Subcloud Ascent on Deep Convective Initiation

This study examines two factors impacting initiation of moist deep convection: free-tropospheric environmental relative humidity (φ E ) and horizontal scale of subcloud ascent (R sub ), the latter exerting a dominant control on cumulus cloud width. A simple theoretical model is used to formulate a “scale selection” hypothesis: that a minimum R sub is required for moist convection to go deep, and that this minimum scale decreases with increasing φ E . Specifically, the ratio of R 2 sub to saturation deficit (1 - φ E ) must exceed a certain threshold value that depends on cloud-layer environmental lapse rate. Idealized, large-eddy simulations of moist convection forced by horizontally varying surface fluxes show strong sensitivity of maximum cumulus height to both φ E and R sub consistent with the hypothesis. Increasing R sub by only 300–400 m can lead to a large increase (>5 km) in cloud height. A passive tracer analysis shows that the bulk fractional entrainment rate decreases rapidly with R sub but depends little on φ E . However, buoyancy dilution increases as either R sub or φ E decreases; buoyancy above the level of free convection is rapidly depleted in dry environments when R sub is small. While deep convective initiation occurs with an increase in relative humidity of the near environment from moistening by earlier convection, the importance of this moisture preconditioning is inconclusive as it is accompanied by an increase in R sub . Overall, it is concluded that small changes to R sub driven by external forcing or by convection itself could be a dominant regulator of deep convective initiation.

54 ENVIRONMENTAL SCIENCES↗

Monsoonal MCS Initiation, Rainfall, and Diurnal Gravity Waves over the Bay of Bengal: Observation and a Linear Model

Abstract Previous observational studies have indicated that mesoscale convective systems (MCSs) contribute the majority of precipitation over the Bay of Bengal (BoB) during the summer monsoon season, yet their initiation and propagation remain incompletely understood. To fill this knowledge gap, we conducted a comprehensive study using a combination of 20-yr satellite observations, MCS tracking, reanalysis data, and a theoretical linear model. Satellite observations reveal clear diurnal propagation signals of MCS initiation frequency and rainfall from the west coast of the BoB toward the central BoB, with the MCS rainfall propagating slightly slower than the MCS initiation frequency. Global reanalysis data indicate a strong association between the offshore-propagating MCS initiation frequency/rainfall and diurnal low-level wind perturbations, implying the potential role of gravity waves. To verify the hypothesis, we developed a 2D linear model that can be driven by realistic meteorological fields from reanalysis. The linear model realistically reproduces the characteristics of offshore-propagating diurnal wind perturbations. The wind perturbations, as well as the offshore propagation signals of MCS initiation frequency and rainfall, are associated with diurnal gravity waves emitted from the coastal regions, which in turn are caused by the diurnal land–sea thermal contrast. The ambient wind speed and vertical wind shear play crucial roles in modulating the timing, propagation, and amplitude of diurnal gravity waves. Using the linear model and satellite observations, we further show that the stronger monsoonal flows lead to faster offshore propagation of diurnal gravity waves, which subsequently control the offshore propagation signals of MCS initiation and rainfall. Significance Statement Rainfall over the Bay of Bengal (BoB) is primarily contributed by large and organized rainfall systems in the summer monsoon season. During this season, these systems are commonly observed over the east coast of India around midnight, the western BoB in the morning, and the central BoB in the afternoon. This eastward rainfall propagation is confirmed by observations, reanalysis data, and a theoretical model to have a strong association with atmospheric diurnal gravity waves. These waves are caused by land–sea thermal contrast and can trigger rainfall systems over the offshore regions. We also found that the diurnal gravity waves, as well as their triggered rainfall systems, can be greatly modulated by the large-scale monsoonal flows. All the above findings improve our understanding of diurnal rainfall cycle over the BoB and may contribute to the future improvement of rainfall forecast over the region.

Meteorology & Atmospheric Sciences↗