Search NASA⌕ Search

SEARCH · Search NASA

Results for “association mapping”

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 91 records · Page 5

Evaluation of coal-associated sediments, wastes, and AMD sludge in the Southern Appalachian Basin as feedstock materials for REE and Li recovery

Critical minerals (CM) such as rare earth elements (REE+) and Lithium (Li) are essential to technological innovation, energy transitions, global economic and defense security, necessitating the search for unconventional resources and efficient recovery methods to avert supply chain disruptions. Here, this study evaluates coal-associated sediments (underclay and roof rock) and wastes from the Pennsylvanian Pottsville Formation of the Southern Appalachian Basin (SAB) as potential feedstocks for CM recovery. A total of 34 samples (15 underclays, 12 roof rocks, 5 Acid Mine Drainage (AMD) sludges, and 2 coal mining wastes) were characterized using XRD, XRF, ICP-MS, and μ-XRF analytical methods. The REE+ and Li concentrations of these materials ranged from 46.8 to 334.4 ppm and from 11.1 to 519 ppm, respectively, with one underclay sample (Hendrix 3456) yielding the highest values for both. Bulk mineralogy for all samples was dominated by aluminosilicate clay phases, particularly illite and kaolinite. All samples exhibited REY def, rel% values >26% and C outl indices that ranged from 0.69 to 0.94, classifying their REE ore potential as Category II (Promising) as defined by Seredin and Dai (2012). Extractability tests (EPA method 3051 A) yielded low REE+ and Li recoveries, with maximum values of 3.3% and 3.6%, respectively, suggesting associations with resistant minerals like clay and phosphates. Elemental mapping indicates that REE+ is associated with phosphate, whereas statistical analysis suggests that REE+ are associated with aluminosilicates, suggesting heterogeneous associations or minimal phosphate contribution. Li also correlated positively with Al 2 O 3 , indicating an aluminosilicate host. This study highlights the potential of coal-associated sediments in the SAB.

Clay minerals↗

NLR diversity and candidate fusiform rust resistance genes in loblolly pine

Resistance to fusiform rust disease in loblolly pine (Pinus taeda) is a classic gene-for-gene system. Early resistance gene mapping in the P. taeda family 10-5 identified RAPD markers for a major fusiform rust resistance gene, Fr1. More recently, single nucleotide polymorphism (SNP) markers associated with resistance were mapped to a full-length gene model in the loblolly pine genome encoding for a nucleotide-binding site leucine-rich repeat (NLR) protein. NLR genes are one of the most abundant gene families in plant genomes and are involved in effector-triggered immunity. Inter- and intraspecies studies of NLR gene diversity and expression have resulted in improved disease resistance. To characterize NLR gene diversity and discover potential resistance genes, we assembled de novo transcriptomes from 92 loblolly genotypes from across the natural range of the species. In these transcriptomes, we identified novel NLR transcripts that are not present in the loblolly pine reference genome and found significant geographic diversity of NLR genes providing evidence of gene family evolution. We designed capture probes for these NLRs to identify and map SNPs that stably cosegregate with resistance to the SC20-21 isolate of Cronartium quercuum f.sp. fusiforme (Cqf) in half-sib progeny of the 10-5 family. We identified 10 SNPs and 2 quantitative trait loci associated with resistance to SC20-21 Cqf. The geographic diversity of NLR genes provides evidence of NLR gene family evolution in loblolly pine. The SNPs associated with rust resistance provide a resource to enhance breeding and deployment of resistant pine seedlings.

59 BASIC BIOLOGICAL SCIENCES↗

How High to Fly? Mapping Evapotranspiration from Remotely Piloted Aircrafts at Different Elevations

Recent advancements in remotely piloted aircrafts (RPAs) have made frequent, low-flying imagery collection more economical and feasible than ever before. The goal of this work was to create, compare, and quantify uncertainty associated with evapotranspiration (ET) maps generated from different conditions and image capture elevations. We collected optical and thermal data from a commercially irrigated potato (Solanum tuberosum) field in the Wisconsin Central Sands using a quadcopter RPA system and combined multispectral/thermal camera. We conducted eight mission sets (24 total missions) during the 2019 growing season. Each mission set included flights at 90, 60, and 30 m above ground level. Ground reference measurements of surface temperature and soil moisture were collected throughout the domain within 15 min of each RPA mission set. Evapotranspiration values were modeled from the flight data using the High-Resolution Mapping of Evapotranspiration (HRMET) model. We compared HRMET-derived ET estimates to an Eddy Covariance system within the flight domain. Additionally, we assessed uncertainty for each flight using a Monte Carlo approach. Results indicate that the primary source of uncertainty in ET estimates was the optical and thermal data. Despite some additional detectable features at low elevation, we conclude that the tradeoff in resources and computation does not currently justify low elevation flights for annual vegetable crop management in the Midwest USA.

54 ENVIRONMENTAL SCIENCES↗

Mapping the Planetary Wake in HD 163296 with Kinematics

Abstract We map the planetary wake associated with the embedded protoplanet creating the CO kink in the disk of HD 163296. We show that the wake can be traced by a series of correlated perturbations in the peak velocity map. The sign change of the perturbations across the disk’s major axis confirms that the wake induces predominantly radial motion, as predicted by models of planet–disk interaction. These results provide the first direct confirmation of planet wakes generated by Lindblad resonances. Mapping the wake provides a constraint on the disk aspect ratio, which is required to measure the mass of the planet.

79 ASTRONOMY AND ASTROPHYSICS↗

Geothermal Play Fairway Analysis, Part 2: GIS methodology

Play Fairway Analysis (PFA) in geothermal exploration originates from a systematic methodology developed within the petroleum industry and is based on a geologic, geophysical, and hydrologic framework of identified geothermal systems. We tailored this methodology to study the geothermal resource potential of the Snake River Plain and surrounding region, but it can be adapted to other geothermal resource settings. We adapted the PFA approach to geothermal resource exploration by cataloging the critical elements controlling exploitable hydrothermal systems, establishing risk matrices that evaluate these elements in terms of both probability of success and level of knowledge, and building a code-based ‘processing model’ to process results. A geographic information system was used to compile a range of different data types, which we refer to as elements (e.g., faults, vents, heat flow, etc.), with distinct characteristics and measures of confidence. Discontinuous discrete data (points, lines, or polygons) for each element were transformed into continuous interpretive 2D grid surfaces called evidence layers. Because different data types have varying uncertainties, most evidence layers have an accompanying confidence layer which reflects spatial variations in these uncertainties. Risk layers, as defined here, are the product of evidence and confidence layers, and are the building blocks used to construct Common Risk Segment (CRS) maps for heat, permeability, and seal, using a weighted sum for permeability and heat, but a different approach with seal. CRS maps quantify the variable risk associated with each of these critical components. In a final step, the three CRS maps were combined into a Composite Common Risk Segment (CCRS) map, using a modified weighted sum, for results that reveal favorable areas for geothermal exploration. Additional maps are also presented that do not mix contributions from evidence and confidence (to allow an isolated view of evidence and confidence), as well as maps that calculate favorability using the product of components instead of a weighted sum (to highlight where all components are present). Our approach helped to identify areas of high geothermal favorability in the western and central Snake River Plain during the first phase of study and helped identify more precise local drilling targets during the second phase of work. By identifying favorable areas, this methodology can help to reduce uncertainty in geothermal energy exploration and development.

15 GEOTHERMAL ENERGY↗

Grassmannian diffusion maps based surrogate modeling via geometric harmonics

Abstract A novel surrogate model based on the Grassmannian diffusion maps (GDMaps) and utilizing geometric harmonics (GH) is developed for predicting the response of complex physical phenomena. The method utilizes GDMaps to obtain a low‐dimensional representation of the underlying behavior of physical/mathematical systems with respect to uncertain input parameters. Using this representation, GH, an out‐of‐sample extension technique, is employed to create a global map from the input parameter space to a Grassmannian diffusion manifold. GH is further employed to locally map points on the diffusion manifold onto the tangent space of a Grassmann manifold. The exponential map is then used to project the points in the tangent space onto the Grassmann manifold, where reconstruction of the full solution is performed. The performance of the proposed surrogate model is verified with three examples. The first problem is a toy example used to illustrate the technique. In the second example, errors associated with the various mappings are assessed by studying response predictions of the electric potential of a dielectric cylinder in a homogeneous electric field. The last example applies the method for uncertainty prediction in the strain field evolution in a model amorphous material using the shear transformation zone theory of plasticity.

42 ENGINEERING↗

Double Kerr-Schild spacetimes and the Newman-Penrose map

The Newman-Penrose map, which is closely related to the classical double copy, associates certain exact solutions of Einstein’s equations with self-dual solutions of the vacuum Maxwell equations. Here we initiate an extension of the Newman-Penrose map to a broader class of spacetimes. As an example, we apply principles from the Newman-Penrose map to associate a self-dual gauge field to the Kerr-Taub-NUT-(A)dS spacetime and we show that the result agrees with previously studied examples of classical double copies. The corresponding field strength exhibits a discrete electric-magnetic duality that is distinct from its (Hodge star) self-dual property.

scattering amplitudes↗

Multiscale hierarchical and heterogeneous mechanical response of additively manufactured novel Al alloy investigated by high-resolution nanoindentation mapping

Smart alloying and microstructural engineering mitigate challenges associated with laser-powder bed fusion additive manufacturing (L-PBFAM). A novel Al–Ni–Ti–Zr alloy utilized grain refinement by heterogeneous nucleation and eutectic solidification to achieve superior performance-printability synergy. Conventional mechanical testing cannot delineate complex micromechanics of such alloys. This study combined multiscale nanomechanical and microstructural mapping to illustrate mechanical signatures associated with hierarchical heat distribution and rapid solidification of L-PBFAM. The disproportionate hardening effect imparted by Al 3 (Ti,Zr) precipitates in the pool boundaries and the semi-solid zone was successfully demonstrated. Nanomechanical response associated with heterogeneity in particle volume fraction and coherency across melt pool was interpreted from nanoindentation force–displacement curves. The hardness map effectively delineated the weakest and strongest sections in the pool with microscopic accuracy. The presented approach serves as a high throughput methodology to establish the chemistry-processing-microstructure-properties correlation of newly designed alloys for L-PBFAM.

36 MATERIALS SCIENCE↗

The intrinsic reddening of the Magellanic Clouds as traced by background galaxies – III. The Large Magellanic Cloud

ABSTRACT We present a map of the total intrinsic reddening across ≃ 90 deg2 of the Large Magellanic Cloud (LMC) derived using optical (ugriz) and near-infrared (IR; YJKs) spectral energy distributions (SEDs) of background galaxies. The reddening map is created from a sample of 222 752 early-type galaxies based on the lephareχ2 minimization SED-fitting routine. We find excellent agreement between the regions of enhanced intrinsic reddening across the central (4 × 4 deg2) region of the LMC and the morphology of the low-level pervasive dust emission as traced by far-IR emission. In addition, we are able to distinguish smaller, isolated enhancements that are coincident with known star-forming regions and the clustering of young stars observed in morphology maps. The level of reddening associated with the molecular ridge south of 30 Doradus is, however, smaller than in the literature reddening maps. The reduced number of galaxies detected in this region, due to high extinction and crowding, may bias our results towards lower reddening values. Our map is consistent with maps derived from red clump stars and from the analysis of the star formation history across the LMC. This study represents one of the first large-scale categorizations of extragalactic sources behind the LMC and as such we provide the lephare outputs for our full sample of ∼ 2.5 million sources.

79 ASTRONOMY AND ASTROPHYSICS↗

RatXcan: A framework for cross-species integration of genome-wide association and gene expression data

Genome-wide association studies (GWAS) have implicated specific alleles and genes as risk factors for numerous complex traits. However, translating GWAS results into biologically and therapeutically meaningful discoveries remains extremely challenging. Most GWAS results identify noncoding regions of the genome, suggesting that differences in gene regulation are the major driver of trait variability. To better integrate GWAS results with gene regulatory polymorphisms, we previously developed PrediXcan (also known as “transcriptome-wide association studies” orTWAS), which maps SNPs to predicted gene expression using GWAS data. In this study, we developed RatXcan, a framework that extends this methodology to outbred heterogeneous stock (HS) rats. RatXcan accounts for the close familial relationships among HS rats by modeling the relatedness with a random effect that encodes the genetic relatedness. RatXcan also corrects for polygenic-driven inflation because of the equivalence between a relatedness random effect and the infinitesimal polygenic model. To develop RatXcan, we trained transcript predictors for 8,934 genes using reference genotype and expression data from five rat brain regions. We found that the cis genetic architecture of gene expression in both rats and humans was sparse and similar across brain tissues. We tested the association between predicted expression in rats and two example traits (body length and BMI) using phenotype and genotype data from 5,401 densely genotyped HS rats and identified a significant enrichment between the genes associated with rat and human body length and BMI. Thus, RatXcan represents a valuable tool for identifying the relationship between gene expression and phenotypes across species and paves the way to explore shared biological mechanisms of complex traits.

Genetics & Heredity↗

Phase diagram of the three-dimensional subsystem toric code

Subsystem quantum error-correcting codes typically involve measuring a sequence of noncommuting parity check operators. They can sometimes exhibit greater fault tolerance than conventional codes, which use commuting checks. However, unlike subspace codes, it is unclear if subsystem codes—in particular their advantages—can be understood in terms of ground-state properties of a physical Hamiltonian. In this paper, we address this question for the three-dimensional subsystem toric code (3D STC), as recently constructed by Kubica and Vasmer [], which exhibits single-shot error correction. Motivated by a conjectured relation between single-shot properties and thermal stability, we study the zero- and finite-temperature phases of an associated noncommuting Hamiltonian. By mapping the Hamiltonian model to a pair of 3D Z 2 gauge theories coupled by a kinetic constraint, we find various phases at zero temperature, all separated by first-order transitions: There are 3D toric code-like phases with deconfined point-like excitations in the bulk, and there are phases with a confined bulk supporting a 2D toric code on the surface when appropriate boundary conditions are chosen. The latter is similar to the surface topological order present in 3D STC. However, the similarities between the single-shot correction in 3D STC and the confined phases are only partial: they share the same sets of degrees of freedom, but they are governed by different dynamical rules. Instead, we argue that the process of single-shot error correction can more suitably be associated with a path (rather than a point) in the zero-temperature phase diagram, a perspective, which inspires alternative measurement sequences enabling single-shot error correction. Moreover, since none of the above-mentioned phases survives at nonzero temperature, the single-shot error-correction property of the code does not imply thermal stability of the associated Hamiltonian phase. Published by the American Physical Society 2024

Li, Yaodong (ORCID:0000000337421944)↗

Regional Assessment for the CO 2 Storage Potential in the Northern Niagaran Pinnacle Reef Trend

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. The report describes the regional analysis of the Northern Niagaran Pinnacle Reef Trend (NNPRT) in Northern Michigan as part of the MRCSP programs Regional Assessment goals. The goal of this task was to perform an initial assessment of the geologic storage capacity and injectivity of the Northern Fairway of Michigan’s Niagaran Reef Trend. This was achieved through three main objectives: 1) development of a reef atlas, 2) understanding the regional trends in geology, and 3) estimating CO 2 and CO 2 -EOR resources. Data collected across multiple subtasks were integrated into a mappable database that could be accessed from many of the leading software. The reef atlas was used to map the reefs and associated characteristics across the entire NNPRT. Several methodologies were explored to estimate the CO 2 and CO 2 -EOR resources across the trend which included three scenarios: 1) storage only using fluid substitution and volumetric estimations, 2) CO 2 -EOR which applied proximity analysis concepts and measured performance metrics from the Core Energy reefs to predict performance at all oil reefs, and 3) enhanced storage scenario which combined CO 2 -EOR with fluid substitution to represent maximized storage after completion of CO 2 -EOR.

01 COAL, LIGNITE, AND PEAT↗

Physics-based adaptivity of a spectral method for the Vlasov–Poisson equations based on the asymmetrically-weighted Hermite expansion in velocity space

We propose a spectral method for the 1D-1V Vlasov–Poisson system where the discretization in velocity space is based on asymmetrically-weighted Hermite functions, dynamically adapted via a scaling α and shifting u of the velocity variable. Specifically, at each time instant an adaptivity criterion selects new values of α and u based on the numerical solution of the discrete Vlasov–Poisson system obtained at that time step. Once the new values of the Hermite parameters α and u are fixed, the Hermite expansion is updated and the discrete system is further evolved for the next time step. The procedure is applied iteratively over the desired temporal interval. The key aspects of the adaptive algorithm are: the map between approximation spaces associated with different values of the Hermite parameters that preserves total mass, momentum and energy; and the adaptivity criterion to update α and u based on physics considerations relating the Hermite parameters to the average velocity and temperature of each plasma species. For the discretization of the spatial coordinate, we rely on Fourier functions and use the implicit midpoint rule for time stepping. The resulting numerical method possesses intrinsically the property of fluid-kinetic coupling, where the low-order terms of the expansion are akin to the fluid moments of a macroscopic description of the plasma, while kinetic physics is retained by adding more spectral terms. Moreover, the scheme features conservation of total mass, momentum and energy associated in the discrete, for periodic boundary conditions. A set of numerical experiments confirms that the adaptive method outperforms the non-adaptive one in terms of accuracy and stability of the numerical solution.

97 MATHEMATICS AND COMPUTING↗

TROPOMI observations allow for robust exploration of the relationship between solar-induced chlorophyll fluorescence and terrestrial gross primary production

Solar-induced chlorophyll fluorescence (SIF) observed by satellites has advanced the monitoring of terrestrial photosynthesis regionally and globally. The relationship between SIF and gross primary production (GPP) at leaf, canopy, and ecosystem scales has received tremendous attention in recent years. It remains controversial whether the SIF-GPP relationship at the ecosystem scale is universal or dependent upon vegetation type. New SIF observations from the TROPOspheric Monitoring Instrument (TROPOMI) with unprecedented high spatial and temporal resolution provide a new opportunity to elucidate the SIF-GPP relationship. Here, we examine the SIF-GPP relationship for seven major vegetation types across the U.S. with TROPOMI SIF and in-situ GPP data for 83 eddy covariance flux sites. We find that TROPOMI SIF shows a strong and consistent relationship with tower-based GPP at both satellite footprint and grid-cell levels. The slope of the SIF-GPP relationship is similar among all the vegetation types except croplands, demonstrating a nearly universal (converging to ~13.5 g C m -2 d -1 /W m -2 μm -1 sr -1) rather than vegetation type-specific SIF-GPP relationship. This confirms that TROPOMI SIF can be used as a proxy for GPP across a wide variety of vegetation types and can also be used to quantify GPP by avoiding uncertainty associated with land cover maps. The C 4 crops have a much higher slope than the C 3 crops, and therefore croplands tend to have a higher slope than C 3 -dominated vegetation types (e.g., forests, shrublands, savannas). We also find that the TROPOMI SIF is well correlated with GPP under normal or wetter conditions, while their relationship becomes weaker under water stress. In conclusion, our TROPOMI-based study could improve our understanding of the SIF-GPP relationship at the ecosystem scale and advance the mapping of GPP globally with SIF observations from space.

54 ENVIRONMENTAL SCIENCES↗

The impact of urban configuration types on urban heat islands, air pollution, CO 2 emissions, and mortality in Europe: a data science approach

The world is becoming increasingly urbanized. As cities around the world continue to grow, it is important for urban planners and policymakers to understand how different urban configuration patterns affect the environment and human health. We aimed at identifying European urban configuration types, based on the Local Climate Zones categories and street design variables from Open Street Map, and evaluating their association with motorized traffic flows, Surface Urban Heat Island (SUHI) intensities, tropospheric nitrogen dioxide (NO 2 ), CO 2 per capita emissions and age-standardized mortality. We considered 946 European cities from 31 countries for the analysis defined in the 2018 Urban Audit database, of which 919 European cities were analysed. Data were collected at a 250 m × 250 m grid cell resolution. We divided all cities into five concentric rings based on the Burgess concentric urban planning model and calculated the mean values of all variables for each ring. First, to identify distinct urban configuration types, we applied the Uniform Manifold Approximation and Projection for Dimension Reduction method, followed by the k-means clustering algorithm. Next, statistical differences in exposures (including SUHI) and mortality between the resulting urban configuration types were evaluated using a Kruskal–Wallis test followed by a post-hoc Dunn's test. We identified four distinct urban configuration types characterising European cities: compact high density (n=246), open low-rise medium density (n=245), open low-rise low density (n=261), and green low density (n=167). Compact high density cities were a small size, had high population densities, and a low availability of natural areas. In contrast, green low-density cities were a large size, had low population densities, and a high availability of natural areas and cycleways. The open low-rise medium and low-density cities were a small to medium size with medium to low population densities and low to moderate availability of green areas. Motorised traffic flows and NO 2 exposure were significantly higher in compact high density and open low rise medium density cities when compared with green low density and open low-rise low density cities. Additionally, green low-density cities had a significantly lower SUHI effect compared with all other urban configuration types. Per person CO 2 emissions were significantly lower in compact high density cities compared with green low density cities. Lastly, green low density cities had significantly lower mortality rates when compared with all other urban configuration types. Our findings indicate that, although the compact city model is more sustainable, European compact cities still face challenges related to poor environmental quality and health. Our results have notable implications for urban and transport planning policies in Europe and contribute to the ongoing discussion on which city models can bring the greatest benefits for the environment, climate, and health.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

First-Principles Elucidation of Initial Dehydrogenation Pathways in Mg(BH 4 ) 2

Complex borohydrides such as Mg(BH 4 ) 2 offer one of highest capacities to chemically store hydrogen for onboard applications; however, it suffers greatly from kinetic constraints that prevent realization of full capacity and reversibility. Understanding these kinetic limitations solely from experiments is extremely challenging due to the unusual complexity of various competing elemental reaction steps involved during the de/rehydrogenation reaction. This work aims to map out the energetics associated with initial dehydrogenation of Mg(BH 4 ) 2 from first-principles simulations and to identify the preferred reaction pathways. Our calculations suggest the rate-limiting step during BH 4 – –B 3 H 8 – conversion is the formation of the B 2 H 7 – intermediate. We further emphasize and clarify that the B 3 H 8 – and H – intermediates, formed during initial Mg(BH 4 ) 2 decomposition, appear as molecular species that are embedded in the Mg–BH 4 –Mg matrix as evidenced in the nuclear magnetic resonance measurements and not as bulk MgH 2 and Mg(B 3 H 8 ) 2 as previously assumed in theoretical predictions of the thermodynamics.

08 HYDROGEN↗

First Analysis of Earth's Stratospheric γ ‐Ray Emission in Geographical Coordinates With Fermi LAT

Abstract From low‐Earth orbit, the Earth is an extremely bright γ ‐ray source compared to distant astrophysical objects. The high‐energy γ ‐ray emission, from the Earth's upper atmosphere, is produced by the interactions between cosmic rays (CRs) and air molecules. As CRs are predominantly charged particles, the geomagnetic field acts as a spectrometer, forbidding CRs below a certain energy to arrive and interact with the atmosphere at a given geomagnetic location from a given direction. We use photon data between 0.2 and 20 GeV as observed by the Fermi Large Area Telescope (LAT) to create, for the first time, stratospheric γ ‐ray emission maps of the Earth in the geographical coordinate system, for different viewing directions from LAT altitude. These maps confirm a strong association between measured γ ‐ray fluxes and geomagnetic cutoffs for CRs. We also report the spectrum of the Earth's stratospheric γ ‐ray emission in cardinal directions and use the geomagnetic spectrometer effect to infer the atmosphere's yield function of γ rays per incoming CR proton as a function of CR rigidity at the altitude of the Fermi spacecraft. This work is relevant to the study of near‐Earth radiation exposure for astronauts and spacecraft, modeling the geomagnetic field, and CR‐air interactions.

Madlee, S.↗