Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model Parallelism”

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.

956 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS2026

Hydrology Copilot: A Cloud-Native Ai System for Hydrological Data Analysis

The emergence of AI-driven Earth observation systems promises to broaden access to petabyte-scale geospatial data beyond domain specialists. However, translating this vision into operational scientific infrastructure requires addressing fundamental challenges in data virtualization, code transparency, and domain-specific reasoning. We present Hydrology Copilot, a cloud-native AI framework for natural-language-driven analysis of Earth observation data. To demonstrate operational capabilities at scale, we implement the system using NASA's North American Land Data Assimilation System version 3 (NLDAS-3), which provides surface meteorological forcing and land-surface model output across North and Central America at 1-km resolution, from which drought diagnostics are derived. The system integrates five core contributions: (1) scalable data virtualization using Kerchunk-based cloud optimized access, achieving a 1.5 to 4.6 times improvement in I/O latency across benchmark queries spanning regional single-day extractions (4.6 times speedup) to continental monthly aggregations (1.5 times speedup); (2) transparent code generation through Microsoft Azure AI Foundry agents that expose executable Python workflows for scientific verification; (3) persistent conversational memory enabling multi-turn analytical discourse across sessions; (4) intelligent query validation that enforces dataset boundaries and resolves ambiguous requests before execution; and (5) a multi-agent architecture coordinating query parsing, code generation, and visualization. We evaluate the system through drought-monitoring workflows, demonstrating reliable code generation, accurate results validated against reference computations and the operational U.S. Drought Monitor, and efficient operation across increasingly complex tasks. By bridging natural-language interfaces with rigorous hydrological analysis, Hydrology Copilot advances beyond proof-of-concept demonstrations to provide a deployable framework for operational Earth science applications.

Data virtualization↗
NASA NTRSDate not supplied

In-Situ Scanning Electron Microscope Experiments for Microscale Mechanical Testing and Validated Modeling of Fiber Reinforced Thermoplastics

A novel, in-situ, scanning electron microscope (SEM) mechanical testing capability for materials at the microscale which provides experimental validation to a machine learning (ML) toolset for full-field validation of physics-based micromechanics models is being developed by researchers at NASA Glenn Research Center. These are enabling technologies for the integration of multiscale digital twins for materials into system level models which will result in the improved performance, material discovery, reduced production cost and time, rapid characterization, and prognostic structural health monitoring (SHM) for materials and structures for extreme environments in support of NASA space exploration missions. In order to bridge the material structure-to-system gap for digital twins, physics-based models must be experimentally validated at multiple length scales. Seminal microscale experiments, conducted at the Air Force Research Laboratory (AFRL), were limited to transverse compression of single-layer, unidirectional thermoset polymer matrix composite (PMC) micropillar specimens [1]. The early phases of the current project followed those initial results and setup to reproduce the compression testing of PMC material on the custom-built piezoelectric actuated micromechanical testing rig built by MicroTesting Solutions LLC. In this work, samples of thermoplastic PMC material were first machined into 3 mm cubes, and then further machining and final milling was done using a Focused Ion Beam (FIB). The initial experiment was done on a pillar roughly 20 µm x 20 µm x 40 µm tall. Additional pillars were milled with final sizes ranging from 20 µm x 20 µm x 40 µm tall to 40 µm x 40 µm x 65 µm tall. A speckle pattern for in-situ full-field measurements using Digital Image Correlation (DIC) was applied with platinum, which was coated on the surface, and then the FIB was used to mill away some of the coating to produce an irregular pattern of Pt on the pillar surface. The samples were loaded into the custom testing rig and placed into the SEM and loaded under compression until failure. Images were collected in the SEM during testing. Post-processing of the images was conducted using DIC to obtain full-field displacement and strain measurements elucidating the role of the matrix as well as fiber-fiber interaction at the microscale within the composite subjected to compression loading well into the non-linear regime of the material. Moreover, the evolution of fiber-matrix debonding and matrix cracking is observed in-situ at the microscale. This data, along with images segmented with a newly developed ML toolset [2], was used to create and validate physics-based micromechanics models. An image of the failed micropillar is shown in Figure 1. The techniques developed in the initial compression experiment was tailored to the validation needs of the models and expanded to include different sized samples as well as possibly tension and fatigue.

Laura Wilson↗
DOE OSTI2025

Egg yolk as a model for gelation: From rheometry to flow physics

Egg yolks are an excellent model for studying sol-gel transitions, particularly the power law viscoelasticity that defines the critical point of gelation. However, prior studies lack comprehensive datasets and fail to visualize flow behavior linked to temperature and time-dependent linear and nonlinear rheology. Here, we present a detailed dataset characterizing egg yolk viscoelasticity across temperature, time, and forcing amplitude using oscillatory shear, step strain, step stress, and constant high strain rate. Novel protorheology visualizations link rheological properties with observable flow behavior. Our findings highlight the nuanced determination of the critical gel point, emphasizing observation timescale dependencies. We compare methods to identify critical temperatures for gelation, including power law viscoelasticity, moduli crossover, diverging zero-shear viscosity, and emerging equilibrium elastic modulus, while visualizing flow consequences near these transitions. Egg yolk is an accessible, realistic, and nontoxic material relevant to the physicist and the chef alike, making it ideal for understanding the rheology of critical gels. By integrating protorheology photos and videos with rigorous rheometric data, we deepen the understanding of critical gels, with broader impacts for studying other materials with sol-gel transitions.

Marsh, Maxwell C. [Department of Mechanical Scienc↗
DOE OSTI2026

Modeling a generic TRISO-fueled heat pipe microreactor using SCALE: Depletion, transportation criticality, and shielding

This paper demonstrates the applicability of the SCALE code system to tristructural-isotropic (TRISO)-fueled heat pipe microreactors through depletion, transportation criticality, and shielding analyses of a generic reference design. The study conducted supports US Nuclear Regulatory Commission code readiness efforts for advanced non–light-water reactor technologies and is intended as a code capability demonstration rather than as an optimization of a specific microreactor design. The modeled reactor employs high-assay low-enriched uranium (HALEU) uranium oxycarbide (UCO) TRISO fuel and beryllium oxide (BeO) reflectors and operates at 7.5 MWth with a nominal lifetime of about 3 effective full power years. Representative cases for fresh and irradiated cores were selected to exercise SCALE methods relevant to reactor operation and post-irradiation transport. The discharged-core decay heat is approximately 6% of operating power immediately after shutdown. Transportation criticality calculations show that internal water ingress is the dominant reactivity effect, with fully flooded fresh core and irradiated core configurations remain above the subcriticality criterion, even with the available control mechanisms. Shielding calculations for a simplified transportation package indicate that normal-condition dose rates are governed mainly by shielding thickness and cooling time, whereas the breached hypothetical accident case is governed primarily by cooling time. Overall, the study shows that SCALE supports depletion, transportation criticality, and shielding evaluations efficiently for TRISO-fueled heat pipe microreactors within a single code system.

Criticality↗
DOE OSTI2024

WRF-Comfort: simulating microscale variability in outdoor heat stress at the city scale with a mesoscale model

Abstract. Urban overheating and its ongoing exacerbation due to global warming and urban development lead to increased exposure to urban heat and increased thermal discomfort and heat stress. To quantify thermal stress, specific indices have been proposed that depend on air temperature, mean radiant temperature (MRT), wind speed, and relative humidity. While temperature and humidity vary on scales of hundreds of meters, MRT and wind speed are strongly affected by individual buildings and trees and vary on the meter scale. Therefore, most numerical thermal comfort studies apply microscale models to limited spatial domains (commonly representing urban neighborhoods with building blocks) with resolutions on the order of 1 m and a few hours of simulation. This prevents the analysis of the impact of city-scale adaptation and/or mitigation strategies on thermal stress and comfort. To solve this problem, we develop a methodology to estimate thermal stress indicators and their subgrid variability in mesoscale models – here applied to the multilayer urban canopy parameterization BEP-BEM within the Weather Research and Forecasting (WRF) model. The new scheme (consisting of three main steps) can readily assess intra-neighborhood-scale heat stress distributions across whole cities and for timescales of minutes to years. The first key component of the approach is the estimation of MRT in several locations within streets for different street orientations. Second, mean wind speed and its subgrid variability are downscaled as a function of the local urban morphology based on relations derived from a set of microscale LES and RANS simulations across a wide range of realistic and idealized urban morphologies. Lastly, we compute the distributions of two thermal stress indices for each grid square, combining all the subgrid values of MRT, wind speed, air temperature, and absolute humidity. From these distributions, we quantify the high and low tails of the heat stress distribution in each grid square across the city, representing the thermal diversity experienced in street canyons. In this contribution, we present the core methodology as well as simulation results for Madrid (Spain), which illustrate strong differences between heat stress indices and common heat metrics like air or surface temperature both across the city and over the diurnal cycle.

Geology↗
NASA NTRSDate not supplied

The Future of Oaks in the Santa Monica Mountains: A Case Study in Using Remote Sensing Data for Species Distributions Models

The Woolsey Fire began on November 8, 2018, and lasted for almost two weeks, during which it burned almost 100,000 acres of valuable landscape and habitat, including a vast area of woodland. The persistence of key woodland species provides aesthetic, monetary, and ecological value to the landscape through carbon sequestration, air temperature moderation, and erosion mitigation, among other ecosystem services. This study investigated the impact of the Woolsey Fire on native woodland species distributions and identified areas suitable for restoration within the Santa Monica Mountains National Recreation Area. The team partnered with the Resource Conservation District of the Santa Monica Mountains; National Park Service, Santa Monica Mountains National Recreation Area; California Department of Parks and Recreation, Los Angeles County Division; County of Los Angeles Fire Department, Prevention Services Bureau, Forestry Division; County of Los Angeles Department of Regional Planning; and the University of Montana. The Earth observations used include data from Landsat 8 Operational Land Imager, NASA ER-2 Jet Airborne Visible InfraRed Imaging Spectrometer, Shuttle Radar Topography Mission, and RapidEye. The team produced maps of burn severity from the Woolsey Fire, its impact on plant species distributions, and habitat suitability projections for 2050 and 2099 to assist partners in prioritizing areas for restoration. A plant community classification was successfully created using Multiple Endmember Spectral Mixture Analysis (MESMA). Overall accuracy was assessed at 90.54% by comparing the classification to validation pixels derived from ground truth information provided by our partners.

Roger Ly↗
NASA NTRSDate not supplied

Predicting Team Functioning in Long Term Space Missions Using Acoustic and Linguistic Measures

Maintaining optimal team functioning is critical for long-duration space exploration missions, yet traditional monitoring methods, such as self-reports and wearable sensors, often impose operational burdens or suffer from bias. This paper investigates a non-intrusive speech-based artificial intelligence (AI) framework to predict degradations in team functioning using data from the Human Exploration Research Analog (HERA) of the U.S. National Aeronautics and Space Administration (NASA). Using acoustic features, linguistic descriptors, and semantic embeddings, we evaluate static non-linear and temporal machine learning models to predict both objective (task accuracy) and subjective (self-reported efficacy and cohesion) team functioning outcomes. Results indicate that temporal models outperform static approaches, with prediction of objective task accuracy in Team Interaction Battery (TIB) improving from near chance to 71%. Self-reported outcomes, including team efficacy and cohesion, are predicted more reliably than task performance, achieving balanced accuracies of up to 85.56% and 78.12%, respectively, and are found to be most strongly associated with acoustic features. In a second interdependent task, the MMSEV–EVA, accuracies of up to 78% are achieved using temporal models with acoustic features. Furthermore, incorporating just 1–2 days of team-specific historical data systematically improved performance, and acoustic markers from informal pre-task interactions provided modest predictive gains. Finally, while automated preprocessing yielded viable accuracy, humancorrected data provided moderate performance gains, though transcription error rates did not significantly correlate with model performance. These findings highlight the potential of speech as a passive, high-fidelity monitoring tool for autonomous habitats.

Temporal modeling↗
DOE OSTI2025

Role of Matter Interactions in Superradiant Phenomena

The superradiant phenomenon, usually described by the Dicke model, is a hallmark of strong light-matter interaction. Here, we explore how matter-matter interactions influence this phenomenon by performing ground-state simulations of Dicke-like models with both isotropic and anisotropic interactions. We find that Ising-type interactions produce two qualitatively distinct phase boundaries, one of which gives rise to an antiferromagnetic-normal phase connected to the superradiant regime via a first-order phase transition. Under anisotropic couplings, we uncover a strongly correlated phase where in-plane spin order coexists with superradiance, exhibiting sublinear scaling of the photon occupation per site and power-law decay of spin correlations. Furthermore, superradiance is strengthened by tuning either isotropic or anisotropic interactions, highlighting the role of intrinsic many-body correlations in shaping light-matter quantum phases.

Dicke model↗
DOE OSTI2026

First-principles DFT modeling of nitrobenzene adsorption on the Ag(111) surface at varying monolayer coverages

Here, the adsorption behavior of nitrobenzene on the Ag(111) surface as a function of coverage is investigated using density functional theory. Adsorption energies and optimized geometries are analyzed together with isolated intermolecular interaction calculations and electronic structure analysis, including Bader charge partitioning and projected density of states, to disentangle the roles of adsorbate–surface bonding and through-space adsorbate–adsorbate interactions. At low and intermediate coverages (θ = 1/9 and θ = 2/9), bidentate adsorption at top sites is favored due to strong adsorbate–surface interactions, with additional stabilization at θ = 2/9 arising from favorable intermolecular separations. At higher coverage (θ = 1/3), bidentate adsorption is destabilized by strong intermolecular repulsion, and monodentate adsorption becomes energetically preferred as rotational freedom allows more favorable intermolecular spacing. Charge density differences, Bader charge, and density of states analyses show that charge transfer from the Ag surface is localized primarily on the nitro group and increases with coverage and adsorption denticity, although this increase does not directly correlate with adsorption strength at high coverage due to competing intermolecular interactions. These results demonstrate that surface coverage can induce a transition in preferred adsorption denticity driven by intermolecular interactions, highlighting the importance of adsorbate packing in organic molecule adsorption at metal interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Generative learning of densities on manifolds

A generative modeling framework is proposed that combines diffusion models and manifold learning to efficiently sample data densities on manifolds. The approach utilizes Diffusion Maps to uncover possible low-dimensional underlying (latent) spaces in the high-dimensional data (ambient) space. Two approaches for sampling from the latent data density are described. The first is a score-based diffusion model, which is trained to map a standard normal distribution to the latent data distribution using a neural network. The second one involves solving an Itô stochastic differential equation in the latent space. Additional realizations of the data are generated by lifting the samples back to the ambient space using Double Diffusion Maps , a recently introduced technique typically employed in studying dynamical system reduction; here the focus lies in sampling densities rather than system dynamics. The proposed approaches enable sampling high dimensional data densities restricted to low-dimensional, a priori unknown manifolds. The efficacy of the proposed framework is demonstrated through a benchmark problem and a material with multiscale structure.

Double diffusion maps↗
DOE OSTI2026

Coalition for Community-Supported Affordable Geothermal Energy Systems (C2SAGES)

The C2SAGES project evaluated the feasibility of a community geothermal system for the planned Windy Ridge affordable housing development in Hinesburg, Vermont. Led by GTI Energy with Vermont Gas Systems, LN Consulting, NREL, and Frontier Energy, the work assessed technical design, energy performance, costs, business models, community engagement, maintenance, workforce development, and permitting. The proposed system was designed to serve 100% of the development’s heating, cooling, and domestic hot water loads. Compared with a baseline using air-source heat pumps and natural gas water heating, the geothermal system was estimated to reduce HVAC and domestic hot water energy use by about 45% to 48%, lower operating and maintenance costs, and reduce 30-year life-cycle costs by 37% for Phase 1 and 10% for Phase 2. Technical testing and modeling indicated that the Windy Ridge site is suitable for a community-scale geothermal system. The project also developed borehole field layouts, piping concepts, pump house designs, controls, maintenance plans, and supporting engineering drawings. The business model analysis found that first cost, ownership structure, and customer affordability remain major deployment challenges. Utility-led maintenance and operation were viewed favorably, but traditional utility cost-recovery models may require subsidy or revised financing structures to be practical for affordable housing. Community engagement highlighted the need for clear public education, transparent financing, reliable long-term maintenance, trained technicians, and the potential to pair geothermal systems with weatherization. Overall, the report concludes that community geothermal is technically feasible and offers meaningful energy, emissions, and life-cycle cost benefits, but broader deployment will depend on workable financing models and workforce readiness.

15 GEOTHERMAL ENERGY↗
NASA NTRS2026

Machine Learning for Predicting Team Functioning in HERA Missions

Team functioning is integral to success in future long term space exploration missions. Proactively detecting declines in team functioning can mitigate conflict and ensure mission success. This project developed a speech-based artificial intelligence (AI) system that unobtrusively predicts degradation in team functioning, including performance and cohesion, in the Human Exploration Research Analog (HERA) Campaigns 4 and 5. The AI system conducted automated analysis of the prosodic (tone of voice) and linguistic (language content) components of speech, modeling interpersonal dynamics at both the turn-taking and day-wide levels. We investigated team functioning via observing structured interactions (i.e., multi-mission space exploration vehicle-extra vehicular activity [MMSEV-EVA], team interaction battery [TIB]) and unstructured interactions before the MMSEV-EVA task. We developed machine learning models to predict team functioning (objective task accuracy, self reported team efficacy and self reported team cohesion) by analyzing OpenSmile acoustic features, linguistic descriptors extracted via the linguistic inquiry and word count (LIWC) dictionary, and semantic embeddings. In the TIB, static models using logistic regression and random forests were not able to predict task accuracy, but predicted team efficacy and cohesion during both the decision making and relational tasks to a moderate level (60-70%). Majority voting on the individual turns to predict day long team efficacy further increased accuracies (70-80%). Finally, long short-term memory (LSTM) models showed the best performance across all variables (80-91%), including task performance. In the MMSEV-EVA, static models achieved an accuracy of 60% with majority voting, which increased to 80% through the incorporation of mission day as a variable, accounting for the learning effect. A key finding across both tasks was the "team-dependent" nature of these interactions; models achieved much higher accuracy when trained on prior days of the same team's data rather than attempting to generalize across entirely different teams, with even 1-2 days of prior data per team achieving 5-15% improvement over team-independent models. In addition, the incorporation of pre-task data from the same team also improves model performance, e.g., incorporating data from the decision-making task of the TIB, which preceded the relational task, improved the prediction of team efficacy and cohesion during the latter. We compared model performance when trained on machine-generated data compared to data that had been further corrected by human annotators. Overall, models trained on human-corrected data exhibited a modest improvement in performance, particularly when acoustic features were used. We found no significant correlation between word error rate (WER) and model accuracy (r(55) = -0.08, p = 0.51), but model’s accuracy was significantly higher for medium/high quality transcription (0.74 (SD = 0.48)) compared to the low-quality group (0.64 (SD = 0.36)) (t(63)=2.82, p = 0.006). Based on these, several design recommendation emerge, that could inform Standards at NASA. Models predicting team functioning should incorporate at least one to two days of historical interaction data, include brief pre-task discussions, and explicitly model temporal learning effects, especially for longer operational tasks. Minimum quality standards for automated speech-processing pipelines are needed, given the performance gains observed with manually corrected acoustic data. Finally, systems should leverage both acoustic features and language embeddings in complementary ways, with modality choices and fusion strategies tailored to mission context, task demands, and data quality requirements.

Shrivatsa Mishra↗
DOE OSTI2026

Paleotribological models using preserved fossil tissue properties reveal functional significance of Eurasian mammoth dental evolution

Eurasian mammoths (Mammuthus) underwent substantial modifications in molar morphology as later-diverging species evolved progressively thinner enamel and increased enamel crest complexity. These features have been hypothesized to reduce whole-tooth wear and extend dental longevity as increasingly graze-dominated diets evolved within the lineage. This hypothesis has yet to be directly tested. Here, in this study, we developed an in-silico wear model using experimentally derived wear rates from fossil and extant proboscidean dental tissues. The models revealed that shifts in tissue topology do not affect whole-tooth wear rate, as inverse trends in lamellar frequency and enamel thickness preserve a consistent surface enamel area fraction; the determining factor of wear. Rather, topological shifts produce a wear-emergent secondary occlusal surface with greater numbers of triturating crests that create a regular, low-relief, file-like shearing pavement. These changes in occlusal architecture likely directly impacted the mastication capacity of Mammuthus dentitions, facilitating their dietary expansion to incorporate fibrous, lower-nutrient graze.

Dental wear↗
DOE OSTI2026

Physical Modeling and Design of a Nonvolatile Optically Gated High‐Power Diamond Transistor

In this work, we present the theory and modeling framework of a diamond optically gated junction field‐effect transistor (DOGFET). The device utilizes nitrogen substitutional centers in type‐1b diamond to optically modulate a p‐ boron doped diamond channel. Using sub‐gap lasers with intensities as low as 100 W/cm 2 , electrons are optically excited from substitutional nitrogen sites to the conduction band of the diamond substrate, thus enabling the optical gate to exercise control on modulating the space‐charge region at the junction and therefore the channel conductivity. We show that the device can deliver a current of 7 μA/μm, or equivalently 1750 A/cm 2 , while switching at a frequency greater than 100 kHz, in a form factor of 5 μm 2 . The breakdown voltage is found to be greater than 1850 V, with a breakdown field strength of ~13 MV/cm. Moreover, the device supports nonvolatile operation with a “memory effect” enabling single transistor state retention. The presented simulation framework provides a physically grounded insight into the limits and opportunities of optoelectronic diamond systems.

Engineering - Electronic and electrical engineerin↗
NASA NTRSDate not supplied

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

The environments experienced by vehicles upon entry into planetary atmospheres generally result in severe aerothermodynamic loading. These flows are characterized by the formation of strong shock waves, behind which high temperatures and non-equilibrium environments are generated, and the accurate prediction of internal energy relaxation and chemical kinetics becomes important. This manuscript focuses on the internal relaxation of the rotational energy of polyatomic molecules. Previous work has described the implementation of a discrete rotational energy model for diatomic molecules (linear rotators). Later, Gimelshein et. al extended this model to include polyatomic molecular internal energies. However, they did not include details on how to generate the particles at equilibrium conditions and centrifugal distortion was not included in their model. The purpose of the present study is to extend the afore-mentioned models to the treatment of discrete rotational energy in polyatomic molecules with the addition of centrifugal distortion. An overview of the background theory needed for the implementation of the model will be given followed by the model itself and verification of the model. For example, sampled versus Boltzmann distributions are compared in Figure 1 for CH3 (oblate rotator – two quantum numbers) and examples of thermal relaxation behavior for linear and non-linear rotators are presented in Figure 2. In addition, comparisons of thermophysical properties will be presented, as well as examples of adiabatic relaxation and application to hypersonic flow.

Rarefied Gas Dynamics↗
DOE OSTI2026

Predicting Atomistic Transitions with Transformers

Accurate knowledge of the atomistic transition pathways in materials and material surfaces is crucial for many material science problems. However, conventional simulation techniques used to find these transitions are extremely computationally intensive. Even with large-scale, accelerated material simulations, the computational cost constrains the applicable domain in practice. Machine learning models, with the potential to learn the complex emergent behaviors governing atomistic transitions as a fast surrogate model, have great promise to predict transitions with a vastly reduced computational cost. Here, we demonstrate how transformers can be trained to predict atomistic transitions in nano-clusters. We show how we evaluate physical validity of the predictions and how a multitude of additional, different microstates can be generated by slightly varying the data provided to the model.

36 MATERIALS SCIENCE↗
NASA NTRS2026

Global Model Estimates of Atmospheric Al, Ca, Fe, Si, and Ti from Dust and Non-Dust Aerosols Informed by EMIT Surface Mineralogy and Evaluated Against Observations

Atmospheric deposition of micro-nutrients like Fe has been shown to be important for ocean biogeochemistry. The largest source of atmospheric Fe and other elements (e.g., Ca, Al, Si, and Ti) is desert dust, although there are significant non-dust sources in some regions. However, past estimates of these elements have been substantially uncertain due to limited information about the composition of the desert source regions. Here we use elemental distributions estimated from new Earth Surface Mineral Dust Source Investigation (EMIT) observations, which provide mineralogical composition at the surface of the Earth based on imaging spectroscopy measurements from the International Space Station. We add in other sources of these elements (anthropogenic and natural) and compare to a compilation of available surface concentration data from stations over land and from shipborne observations. Our results suggest that the modeled distribution is similar to available observations, but discrepancies still exist in both natural desert dust regions as well as regions dominated by anthropogenic sources. Global budgets for the elements Ca, Al, Fe, Si, and Ti suggest that desert dust remains the dominant source for these elements but anthropogenic or volcanic sources are also important for these elements. Changes in elemental distributions since preindustrial times were also estimated.

aerosols↗