Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distribution factor”

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 1,045 records · Page 58

Hard-photon-triggered jets in 𝑝−𝑝 and 𝐴−𝐴 collisions

An investigation of high-transverse-momentum (high-𝑝 𝑇 ) photon-triggered jets in proton-proton (𝑝−𝑝) and ion-ion (𝐴−𝐴) collisions at $\sqrt{s_{NN}}$=0.2 and 5.02TeV is carried out, using the multistage description of in-medium jet evolution. Monte Carlo simulations of hard scattering and energy loss in heavy-ion collisions are performed using parameters tuned in a previous study of the nuclear modification factor (𝑅 𝐴⁢𝐴 ) for inclusive jets and high-𝑝𝑇 hadrons. We obtain a good reproduction of the experimental data for photon-triggered jet 𝑅 𝐴⁢𝐴 , as measured by the ATLAS detector, the distribution of the ratio of jet to photon 𝑝 𝑇 (𝑋 𝐽⁢𝛾 ), measured by both CMS and ATLAS, and the photon-jet azimuthal correlation as measured by CMS. We obtain a moderate description of the photon-triggered jet 𝐼 𝐴⁢𝐴 , as measured by STAR. A noticeable improvement in the comparison is observed when one goes beyond prompt photons and includes bremsstrahlung and decay photons, revealing their significance in certain kinematic regions, particularly at 𝑋 𝐽⁢𝛾 >1. Moreover, azimuthal angle correlations demonstrate a notable impact of bremsstrahlung photons on the distribution, emphasizing their role in accurately describing experimental results. This work highlights the success of the multistage model of jet modification to straightforwardly predict (this set of) photon-triggered jet observables. This comparison, along with the role played by bremsstrahlung photons, has important consequences on the inclusion of such observables in a future Bayesian analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Final cooling with thick wedges for a muon collider

In the final cooling stages for a muon collider, the transverse emittances are reduced while the longitudinal emittance is allowed to increase. In previous studies, Final 4-D cooling used absorbers within very high field solenoids to cool low-momentum muons. Simulations of the systems did not reach the desired cooling design goals. In this study, we develop and optimize a different conceptual design for the final 4D cooling channel, which is based on using dense wedge absorbers. We used G4Beamline to simulate the channel and Python to generate and analyze particle distributions. We optimized the design parameters of the cooling channel and produced conceptual designs (corresponding to possible starting points for the input beam) which achieve transverse cooling in both x and y by a factor of ~3.5. These channels achieve a lower transverse and longitudinal emittance than the best design previously published.

43 PARTICLE ACCELERATORS↗

Nucleon and pion PDFs: large-x asymptotics meets functional mimicry

We review strategies to unveil the primordial large- x x structure of the nucleons as well as the pion from hard-scattering experiments. Ideas are presented for learning about the x\to 1 x → 1 limit of nonperturbative QCD dynamics at energy scales of order 1 GeV from collider experiments at much higher scales. The behavior of parton distributions at x\to 1 x → 1 predicted by the quark counting rules and other low-energy theoretical approaches is contrasted with phenomenological PDFs. Polynomial mimicry of PDF parametrizations is one of many factors that influence the apparent power of the (1-x) ( 1 − x ) falloff. We discuss implications of the mimicry for the large- x x falloff of the pion PDFs.

Courtoy, Aurore↗

Develop and Connect TRISO Failure Analysis and Uncertainty Quantification to Fission Product Release Calculation Capability

The U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel cycle systems. This program has been providing engineering-scale support for the development of BISON, a high-fidelity and high-resolution fuel performance tool. Stress-based failure probability has been developed and analyzed to assess the integrity of tri-structural isotropic (TRISO) fuel particles during fuel life cycles. While simple, stress-based approaches to failure probability leveraging the Weibull statistical distribution entails a number of drawbacks when stress concentration occurs near crack tips, including finite element mesh size dependency. In this report, we use an interaction integral approach to the computation of stress intensity factors in functionally graded materials (FGM) for axisymmetric models. The inner pyrolytic carbon (IPyC) cracking induced silicon carbide (SiC) failure is one of the dominated failure modes in TRISO failure analysis. In this study, we consider a crack in the IPyC layer perpendicular to the SiC layer. The interface between these two TRISO layers is considered to be porous, which we simulate considering a transition of mechanical properties over the porous length. These aspects are considered in the computation of stress intensity factor (SIF) from a fracture mechanics approach and compared with the known stress-based failure probability approach.

42 ENGINEERING↗

Hofmann Stability Charts Revisited for PIP-II: From Classical Theory to Assumption-Free and ML-Driven Maps

The Hofmann stability chart remains a standard for visualizing parametric resonances in space-charge–dominated linacs, but its use typically relies on non-oscillatory Vlasov dispersion relations with simplifying assumptions (continuous focusing, KV phase space, linear optics, limited transverse–longitudinal coupling). We revisit the chart for the PIP-II linac along three tracks. (1) We reproduce the conventional maps in the (νz/νx, νx/ν0x) plane for relevant εz/εx, providing a validated reference. (2) We remove key assumptions by deriving stability surfaces directly from multi-particle tracking with realistic lattice discreteness, RF defocusing, solenoid/quad optics, and bunched-beam dynamics; local tunes and early-time growth rates are estimated from envelope oscillations and projected to the same coordinates. These assumption-reduced maps recover the canonical stopbands while revealing shifts and broadenings driven by tune modulation, non-KV distributions, and transverse–longitudinal coupling at PIP-II intensities. (3) We train a compact machine-learning surrogate that emulates the growth surface from zero-current optics, tune depression, emittance ratio, bunching factor, and selected lattice descriptors, enabling rapid scans and online working-point selection. We compare the three representations on representative PIP-II sections and discuss implications for commissioning guard bands, resonance avoidance, and routine operations.

Pathak, Abhishek [Fermilab] (ORCID:000000021704208↗

From noise to information: The transfer function formalism for uncertainty quantification in reconstructing the nuclear density

The neutron distribution of neutron-rich nuclei provides critical information on the structure of finite nuclei and neutron stars. Parity violating experiments—such as PREX and CREX—provide a clean and largely model-independent determination of neutron densities. Such experiments, however, are challenging and expensive, which is why sound statistical arguments are required to maximize the information gained. We introduce a new framework, the “transfer function formalism,” aimed at uncertainty quantification, model selection, and experimental design in the context of neutron densities. The transfer functions (TFs) are built analytically by expressing the linear response of the objective function (e.g., χ2) to small perturbations of the data. Using the TF formalism, we are able to analyze the expected overall uncertainty—quantified in terms of bias and variance—of the mean square radius and interior density of 48 Ca and 208 Pb. Using relativistic mean field models as a proxy for the weak-charge density—and assuming that a total of five measurements could be performed on the weak form factor of 48 Ca and 208 Pb—we identify the optimal models and experimental locations that minimize the uncertainty in the extraction of the radius and interior density. We also explore the use of the TF formalism to understand the influence of prior distributions for the model parameters, as well as the optimization of model hyperparameters not constrained by the data. Here,we establish how the choice of experimental locations and the model that is used can have a significant impact on the final uncertainties of the extracted quantities of interest. For challenging experiments such as CREX and PREX, a proper quantification of such uncertainties is critical. We have demonstrated how the TF formalism provides several advantages for this type of analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Columnar Localization and Laminar Origin of Cortical Surface Electrical Potentials

Electrocorticography (ECoG) methodologically bridges basic neuroscience and understanding of human brains in health and disease. However, the localization of ECoG signals across the surface of the brain and the spatial distribution of their generating neuronal sources are poorly understood. To address this gap, we recorded from rat auditory cortex using customized μECoG, and simulated cortical surface electrical potentials with a full-scale, biophysically detailed cortical column model. Experimentally, μECoG-derived auditory representations were tonotopically organized and signals were anisotropically localized to less than or equal to ±200 μm, that is, a single cortical column. Biophysical simulations reproduce experimental findings and indicate that neurons in cortical layers V and VI contribute ∼85% of evoked high-gamma signal recorded at the surface. Cell number and synchrony were the primary biophysical properties determining laminar contributions to evoked μECoG signals, whereas distance was only a minimal factor. Thus, evoked μECoG signals primarily originate from neurons in the infragranular layers of a single cortical column. SIGNIFICANCE STATEMENT ECoG methodologically bridges basic neuroscience and understanding of human brains in health and disease. However, the localization of ECoG signals across the surface of the brain and the spatial distribution of their generating neuronal sources are poorly understood. We investigated the localization and origins of sensory-evoked ECoG responses. We experimentally found that ECoG responses were anisotropically localized to a cortical column. Biophysically detailed simulations revealed that neurons in layers V and VI were the primary sources of evoked ECoG responses. These results indicate that evoked ECoG high-gamma responses are primarily generated by the population spike rate of pyramidal neurons in layers V and VI of single cortical columns and highlight the possibility of understanding how microscopic sources produce mesoscale signals.

Baratham, Vyassa L.↗

Measurements of the suppression and correlations of dijets in Pb+Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

Studies of the correlations of the two highest transverse momentum (leading) jets in individual Pb+Pb collision events can provide information about the mechanism of jet quenching by the hot and dense matter created in such collisions. In Pb+Pb and $pp$ collisions at $\sqrt{s_{NN}}$ = 5.02 TeV, measurements of the leading dijet transverse momentum ($p_T$) correlations are presented. Additionally, measurements in Pb+Pb collisions of the dijet pair nuclear modification factors projected along leading and subleading jet $p_T$ are made. The measurements are performed using the ATLAS detector at the LHC with 260 pb –1 of pp data collected in 2017 and 2.2 nb –1 of Pb+Pb data collected in 2015 and 2018. An unfolding procedure is applied to the two-dimensional leading and subleading jet $p_T$ distributions to account for experimental effects in the measurement of both jets. Results are provided for dijets with leading jet $p_T$ greater than 100 GeV. Measurements of the dijet-yield-normalized $x_{J}$ distributions in Pb+Pb collisions show an increased fraction of imbalanced jets compared to pp collisions; these measurements are in agreement with previous measurements of the same quantity at 2.76 TeV in the overlapping kinematic range. Measurements of the absolutely normalized dijet rate in Pb+Pb and pp collisions are also presented, and show that balanced dijets are significantly more suppressed than imbalanced dijets in Pb+Pb collisions. It is observed in the measurements of the pair nuclear modification factors that the subleading jets are significantly suppressed relative to leading jets with $p_T$ between 100 and 316 GeV for all centralities in Pb+Pb collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Machine-Learning-Based Mapping and Modeling of Solar Energy with Ultra-High Spatiotemporal Granularity

Despite the rapid growth of solar energy, we still lack a dynamic, high-fidelity database that tracks the spatiotemporal variations of solar PVs and their associated infrastructures across different places at a spatially resolved scale. The absence of such data presents a barrier to various applications such as solar PV growth projection, solar energy integration, solar incentive design, and climate risk assessment. In this project, we aim to bridge this gap by developing AI-based algorithms to extract granular information about solar PV installations and their associated infrastructures (i.e., distribution grids) from widely available unstructured data like remote sensing images and street views. As a result, we have built the Solar Energy Atlas, a fine-grained, large-scale geospatial overlay of distributed solar PVs and distribution grids. On top of it, we have advanced the understanding of solar adoption and distribution grid vulnerability to climate-induced extremes. Our major contributions can be summarized as follow: (1) By developing new AI algorithms, we have built the most comprehensive solar PV spatiotemporal database covering the entire US. This is the first time we obtained the exact GPS locations, size, subtype, and installation year information for rooftop solar PVs across the US. This database can be used for solar PV growth projection, solar energy integration, solar energy policy analysis and design, and spatially-resolved climate risk assessment. (2) Leveraging this database, we have uncovered the socioeconomic driving factors that are correlated with earlier onset of solar adoption and higher saturated adoption levels. We have identified the heterogeneity in the effects of different types of financial incentives on solar adoption and provided implications for tailoring incentive design based on local income levels to promote equitable solar adoption. (3) We have developed a distribution grid GIS mapping algorithm which can obtain granular geospatial and topology information about distribution grids using multi-modal open data, reducing the dependency on hard-to-obtain smart meter data of conventional approaches. It shows effectiveness in both the U.S. and Sub-Saharan Africa. Using this algorithm, we have uncovered the non-uniform vulnerability of distribution grids to wildfires in California in the aspects of undergrounding protection and Distributed Energy Resources (DER) preparedness. This has provided important implications for improving the affordability and equity of grid adaptation approaches. (3) We have made our produced database publicly available and provided user-friendly interface to enable various stakeholders and the general public to interact with the data. We have also integrated the produced data into the Data Commons platform to enable the public to access the data and correlate it with other location-specific characteristics simply using natural language as queries. The impact of our project is three-fold: (1) New algorithms for mapping solar PVs and distribution grids across space and time, which are open source to facilitate researchers and industry; (2) New databases of solar PVs and distribution grids that have been made publicly available for engineering, social, and policy applications; (3) New understandings and actionable insights on the potential approaches to promoting solar adoption and reducing energy infrastructure vulnerabilities. In this report, we start by discussing the project background and motivation (section 5), followed by the overview of project objectives (section 6). Results and discussion for each task are presented in section 7. Significant accomplishments are summarized in section 8. This report will be concluded by discussing the paths forwards (section 9), products (section 10), and team roles (section 11).

14 SOLAR ENERGY↗

Connected Communities: A Multi-Building Energy Management Approach

Connected Communities (CCs) are collections of buildings and distributed energy resources (DERs) that incorporate integrated energy management strategies at the multi-building scale. Based on the research conducted for this report, RMI and NREL have identified examples of CCs that are enabling increasingly integrated solutions, sharing resources among buildings, and providing services back to the grid. As demonstrated by these examples, CCs can provide real and scalable solutions to mitigate energy costs, enable more flexibility in grid operations, and offer choices in meeting greenhouse gas (GHG) emissions goals. Nationwide examples of these solutions are included in the appendix of this report. In order to understand the most promising practices for the design, operation, and business models associated with CCs, a search was performed to create a nationwide inventory of CCs and multiple stakeholder interviews were conducted. Based on that research, this paper explores the many factors that influence the development and operation of CCs and evaluates their potential value.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Outdoor Test Bed Performance of a Power Line Sensor Using a Real-Time Event Simulator

This report summarizes the design and application of an outdoor power line sensor testbed (OPLST) with a real-time simulator and power meter to compare a potential transformer (PT) and current transformer (CT) versus an advanced line-post sensor. The OPLST was created to validate advanced medium-voltage (20/34.5 kV) outdoor power line sensors used in electrical distribution systems. Electrical utilities have installed metering and relay-protection transformers like PTs and CTs for several decades. The PTs/CTs are iron core measurement transformers based on the electromagnetic induction principle and provide reliable data in normal grid operation. However, new outdoor power line sensors (OPLS) using other technologies like voltage divider, Rogowski coil and optical principles have become available and may have favorable performance and costs compared to PTs/CTs. Therefore, the importance of testing these technologies with the PT/CT, to compare the measured phase voltage/current at different power grid scenarios is crucial to understand the performance of these new OPLS. For this study a G&W Electric Model CVS-36-O power line sensor was chosen as the OPLS to test with voltage/current signals. An OPAL-RT Technologies Model OP4510 real-time simulator and Schweitzer Engineering Laboratories Model SEL-735 power meter were installed with the 20/34.5 kV OPLST to compare the measured transient events collected from an advanced OPLS and the PT/CT. This system is installed at the Distributed Energy Communications and Control (DECC) lab, Oak Ridge National Laboratory (ORNL). The simulator generated different power grid scenarios (electrical faults, capacitor bank operation, service restoration, etc.), and its analog-output signals were connected to the voltage/current amplifiers that feed the 20/34.5 kV aerial cable loop through the PT/CT devices. Additional PT/CT devices were also wired with the medium voltage aerial cable loop to measure the phase current/voltage signals and server as references. After each test, common format for transient data exchange (COMTRADE) files were collected from the SEL-735 power meter and used to compare the performance of the OPLS with the PT/CT. The behavior of analog signals, harmonic components, total harmonic distortion and crest factors were assessed and found favorable for the G&W sensor as compared to the reference PT/CT.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Discovering hidden geothermal signatures using non-negative matrix factorization with customized k-means clustering

Discovery of hidden geothermal resources is challenging. It requires the mining of large datasets with diverse data attributes representing subsurface hydrogeological and geothermal conditions. The commonly used play fairway analysis approach typically incorporates subject-matter expertise to analyze regional data to estimate geothermal characteristics and favorability. We demonstrate an alternative approach based on machine learning (ML) to process a geothermal dataset from southwest New Mexico (SWNM). The study region includes low- and medium-temperature hydrothermal systems. Several of these systems are not well characterized because of insufficient existing data and limited past explorative work. This study discovers hidden patterns and relations in the SWNM geothermal dataset to improve our understanding of the regional hydrothermal conditions and energy-production favorability. This understanding is obtained by applying an unsupervised ML algorithm based on non-negative matrix factorization coupled with customized k-means clustering (NMFk). NMFk can automatically identify (1) hidden signatures characterizing analyzed datasets, (2) the optimal number of these signatures, (3) the dominant data attributes associated with each signature, and (4) the spatial distribution of the extracted signatures. Here, in this study, NMFk is applied to analyze 18 geological, geophysical, hydrogeological, and geothermal attributes at 44 locations in SWNM. Using NMFk, we find data patterns and identify the spatial associations of hydrothermal signatures within two physiographic provinces (Colorado Plateau and Basin and Range) and two sub-regions of these provinces (the Mogollon-Datil volcanic field and the Rio Grande rift) in SWNM. The ML algorithm extracted five hydrothermal signatures in the SWNM datasets that differentiate between low (<90°C) and medium (90-150°C)-temperature hydrothermal systems. The algorithm also suggests that the Rio Grande rift and northern Mogollon-Datil volcanic field are the most favorable regions for future geothermal resource discovery. NMFk also identified critical attributes to identify medium-temperature hydrothermal systems in the study area. The resulting NMFk model can be applied to predict geothermal conditions and their uncertainties at new SWNM locations based on limited data from unexplored regions. The code to execute the performed analyses as well as the corresponding data can be found at https://github.com/SmartTensors/GeoThermalCloud.jl.

15 GEOTHERMAL ENERGY↗

Anthropogenic Contributions to the 2021 Pacific Northwest Heatwave

Abstract Daily maximum temperatures during the 2021 heatwave in the Pacific Northwest United States and Canada shattered century old records. Multiple causal factors, including anthropogenic climate change, contributed to these high temperatures, challenging traditional methods of attributing human influence. We demonstrate that the observed 2021 daily maximum temperatures are far above the bounds of Generalized Extreme Value distributions fitted from historical data. Hence, confidence in Granger causal inference statements about the human influence on this heatwave is low. Alternatively, we present a more conditional hindcast attribution study using two regional models. We performed ensembles of simulations of the heatwave to investigate how the event would have changed if it had occurred without anthropogenic climate change and with future warming. We found that global warming caused a ∼0.8°C–1°C increase in heatwave temperatures. Future warming would lead to a ∼5°C increase in heatwave temperature by the end of the 21st century.

54 ENVIRONMENTAL SCIENCES↗

Quantum and anharmonic effects in non-adiabatic transition state theory

Quantitative descriptions of non-adiabatic transition rates at intermediate temperatures are challenging due to the simultaneous importance of quantum and anharmonic effects. Here, in this paper, the interplay between quantum effects—for motion across or along the seam of crossing—and anharmonicity in the seam potential is considered within the weak coupling limit. The well-known expression for quantized 1-D motion across the seam (i.e., tunneling) in the linear terms approximation is derived in the thermal domain using the Lagrangian formalism, which is then applied to the case when tunneling is distributed along the seam of crossing (treating motion along the seam classically). For high-frequency quantum modes, a vibrationally adiabatic (VA) approach is developed that introduces to the non-adiabatic rate constant a factor associated with high-frequency wavefunction overlap; this approach treats the high-frequency motion along the seam quantum mechanically. To test these methodologies, the reaction N 2 O ↔ N 2 + O( 3 P) was chosen. CCSD(T)-F12b/cc-pVTZ-F12 explorations of the 3 A'- 1 A' seam of N 2 O revealed that seam anharmonicity has a strong effect on the rate constant (a factor of ~20 at 2000 K). Several quantum effects were found to be significant at intermediate/lower temperatures, including the quantum N–N vibration that was coupled with seam anharmonicity using the VA approach. Finally, a 1-D approximation to non-adiabatic instanton theory is presented to estimate the validity limit of the linear terms model at low temperatures (⁠~ 250 K for N 2 O). We recommend that the assumptions built into many statistical theories for non-adiabatic reactions—harmonic behavior, classical motion, linear terms, and weak coupling—should be verified on a case-by-case basis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interpretative transport modeling of the WEST boundary plasma: main plasma and light impurities

Understanding impurity transport in tokamak plasmas is crucial to control radiative losses and material migration in future magnetic fusion reactors. In this work we deploy the SolEdge2D-EIRENE code to model the boundary plasma in a WEST discharge, satisfactorily reproducing measurements of both upstream and divertor plasma conditions. Additionally, the spatial distribution of oxygen, studied here as a representative light impurity, is compared to vacuum ultraviolet spectroscopy measurements acquired with an oscillating line of sight. The simulation captures a key feature of the experiment, namely a factor of sime2 higher oxygen brightness in the inner divertor region compared to the outer one. This spatial asymmetry in oxygen concentration is interpreted by analyzing the balance of friction forces and thermal gradient forces that the light impurity exchanges with the main plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hydrogen Blending into Natural Gas Pipeline Infrastructure: Review of the State of Technology

Hydrogen is an energy carrier that could play an important role in reducing emissions associated with difficult-to-decarbonize sectors including peaking and load-following electricity and industrial heating. Blending hydrogen into the natural gas pipelines has been proposed as an approach for achieving near-term emissions reductions and early-market access for hydrogen technologies such as electrolyzers. Numerous challenges and uncertainties complicate this approach to natural gas decarbonization, however, and this review summarizes current research on the material, economic, and operational factors that must be considered for hydrogen blending. First, this review explores previous research regarding the effects of blending hydrogen on gas mixture fluid and thermodynamic properties, pipeline materials and equipment performance within transmission and distribution networks, and supporting facilities such as underground storage and end-use hydrogen separation. We also investigate and summarize studies that developed mathematical models of natural gas pipeline networks with hydrogen blending, and the operational and techno-economic findings of these network studies. Finally, we discuss notable hydrogen blending demonstrations and their key outcomes. Many blending demonstrations internationally have proven that low hydrogen percentage blending is feasible under very specific scenarios with limited end-usage applications on both high-pressure transmission lines and low-pressure distribution lines. This report summarizes findings from literature into key areas of consensus and disagreement. Research gaps and disagreements between the literature are highlighted to provide directions for future hydrogen blending research.

03 NATURAL GAS↗

Geometry-Dependent Solidification Regimes in Metal Additive Manufacturing

Recent modeling and experimental work in additive manufacturing has suggested that cross-sectional geometry may play a significant role in the local development of the solidification structure through its influence on the heat source path. This effect has been rationalized as the transition from a quasi-static point heat source regime, to a regime dominated by the quasi-static motion of an equivalent line source. This work provides a simple analytical framework for determining the conditions under which a system transitions between these regimes. Here, a transient semi-analytical heat transfer model is used to examine a wide range of process conditions and material properties. A simple analytical expression is derived and shown to accurately predict the transition between solidification regimes over these conditions. The functional form of this expression is then used to help understand the importance of various material properties, process parameters, and geometric factors on the characteristics of the solidification conditions. This approach may be used as a simple guideline for optimizing process conditions in response to variations in cross-sectional geometry to produce more consistent microstructural distributions in additively manufactured components.

36 MATERIALS SCIENCE↗

Atmospheric particle abundance and sea salt aerosol observations in the springtime Arctic: a focus on blowing snow and leads

Abstract. Sea salt aerosols play an important role in the radiation budget and atmospheric composition over the Arctic, where the climate is rapidly changing. Previous observational studies have shown that Arctic sea ice leads are an important source of sea salt aerosols, and modeling efforts have also proposed blowing snow sublimation as a source. In this study, size-resolved atmospheric particle number concentrations and chemical composition were measured at the Arctic coastal tundra site of Utqiaġvik, Alaska, during spring (3 April–7 May 2016). Blowing snow conditions were observed during 25 % of the 5-week study period and were overpredicted by a commonly used blowing snow parameterization based solely on wind speed and temperature. Throughout the study, open leads were present locally. During periods when blowing snow was observed, significant increases in the number concentrations of 0.01–0.06 µm particles (factor of 6, on average) and 0.06–0.3 µm particles (67 %, on average) and a significant decrease (82 %, on average) in 1–4 µm particles were observed compared to low wind speed periods. These size distribution changes were likely caused by the generation of ultrafine particles from leads and/or blowing snow, with scavenging of supermicron particles by blowing snow. At elevated wind speeds, both submicron and supermicron sodium and chloride mass concentrations were enhanced, consistent with wind-dependent local sea salt aerosol production. At moderate wind speeds below the threshold for blowing snow as well as during observed blowing snow, individual sea spray aerosol particles were measured. These individual salt particles were enriched in calcium relative to sodium in seawater due to the binding of this divalent cation with organic matter in the sea surface microlayer and subsequent enrichment during seawater bubble bursting. The chemical composition of the surface snowpack also showed contributions from sea spray aerosol deposition. Overall, these results show the contribution of sea spray aerosol production from leads on both aerosols and the surface snowpack. Therefore, if blowing snow sublimation contributed to the observed sea salt aerosol, the snow being sublimated would have been impacted by sea spray aerosol deposition rather than upward brine migration through the snowpack. Sea spray aerosol production from leads is expected to increase, with thinning and fracturing of sea ice in the rapidly warming Arctic.

54 ENVIRONMENTAL SCIENCES↗