Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface complexation modeling”

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,009 records · Page 56

Three and Two Phase Rotating Field Inductive Couplers for Wireless Power Transfer with One Phase per Layer Windings

Multiphase inductive wireless charging coils have been proposed recently to improve coupler surface power density, reduce component stress and size, and provide near-constant power delivery to charge mobile electric systems. Several aspects for the fundamental characterization of multiphase coils are explored up to six phases including approximate mutual inductance with size and turn variation, induced voltage, and output power estimation. The relative component stress and size of passive components for resonant operation are compared between the multiphase variants. A combination of an experimentally validated 3D electromagnetic finite element analysis (FEA) and power electronic co-simulations are used to validate the estimated quantities approximated with a mixture of analytical equations and an artificial neural network model for mutual inductance. A novel three-phase transmitter, two-phase receiver coil pair is also proposed for electric vehicle charging to reduce the number of connections and compensation complexity on the vehicle-side with improved power output compared to a two-phase configuration.

Lewis, Donovin D. [University of Kentucky]↗

Improved Infrared Optical Constants from Pressed Pellets: II. Ellipsometric n and k Values for Ammonium Sulfate with Variability Analysis

Infrared spectral imaging is facilitated via the comparison of spectra recorded in situ to a databank of actual or synthetic infrared reflectance spectra. It has recently been shown that reference spectra corresponding to the many different morphological forms of the same chemical can be generated synthetically using the imaginary, k, and real, n, components of the complex refractive index, n~ = n + ik. One method to obtain the n and k vectors is infrared ellipsometry, which measures the changes in amplitude, tan ?, and phase, ?, of polarized light reflected from the sample both as a function of wavenumber and angle of incidence. The method requires specularly reflected light so best results are usually obtained with polished planar samples of large surface area. Due to the difficulties of obtaining such samples, however, we investigate the possibility of pressing powders of neat materials and obtaining the corresponding optical constants from the pellets. In this paper, variability in the sample pellet and preparation method is investigated, as is variability in the fitting procedure for the derived optical constants. The n/k vectors are derived from the measured ellipsometric parameters, tan ? and ?, as they are fit by an oscillator model which yield n(?) and k(?) vectors as a function of wavenumber, ?. Construction of the oscillator model is not automatic and depends on significant input from the analyst as well as the sample’s physical characteristics. For the pellet pressing, the experimental variability was found to be minimized for size-selected powdered samples as gauged by the minimal variance in ? and ? for three different pellets; similarly, the analytical precision for multiple measurements of the same pellet was also quite good, suggesting that a pressed pellet is a viable sample preparation method. The greatest variability came in the analytic fitting procedure with differences up to 20% arising from four different fits to the same data set. The final ellipsometric n/k data are compared to literature values obtained from crystalline (NH4)2SO4 samples as well as single-angle reflectance measurements that also used pressed pellets. Comparison with the previous literature values shows generally good agreement although higher k-values are observed for the independent sets of data derived from pressed pellets. These data are suggested as an improved set of optical constants for (NH4)2SO4.

optical constants, ellipsometry, ammonium sulfate,↗

Advanced Instrumentation for Metal Additive Manufacturing

Laser powder bed fusion (LPBF) is the most widely used process for metal additive manufacturing (AM), particularly where complex geometries provide performance advantages unattainable with traditional manufacturing techniques. However, LPBF is highly sensitive to innate variability in both the powder spreading and fusion steps, often leading to defects such as pores that are difficult to detect yet significantly impair component mechanical properties and fatigue life. This thesis presents a range of novel instruments enabling both precise assessment of powder layer characteristics and in-situ thermal metrology of metal AM to advance the quality control of LPBF. First, leveraging a custom X-ray microscope and a radiation-transport model developed through this work, transmission X-ray imaging is used to study spreading of thin metal powder layers. Effective layer depth is directly mapped at a process-relevant size scale, surpassing optical techniques that can only estimate local deposition from layer surface topography. Layer packing density and quality are shown to be influenced by powder flowability and particle size relative to nominal powder layer thickness. Layer quality is additionally connected to the geometry of the spreading implement and its velocity. This technique and its presented findings enable pairing feedstocks with spreading strategies that create layers with consistent packing density and uniformity. Second, a twofold approach is employed to optically interrogate the laser fusion step of LPBF for observing signatures of defect formation. Aperture division multiplexing is conceptualized, providing for simultaneous laser delivery and high-fidelity infrared (IR) process monitoring through a common optic. In-situ microscopy at 50 μm spatial resolution and at mid-wave IR wavelengths is proven readily achievable with the first purpose-built optic of this type. Next, a bespoke imaging spectrometer, along with a temperature-emissivity separation technique, is used to retrieve accurate process temperatures over a 1000 K range. Data from these instruments are correlated to porosity as fine as 4.3 μm in two LPBF test artifacts, as verified using computed tomography (CT), establishing the viability of robust optically-based component qualification.

Penny, Ryan↗

Role of Capping Agents in the Synthesis of Salicylate-Capped Zinc Oxide Nanoparticles

Capping agents are often used for controlling the size, aggregation, and properties of nanoparticles. To guide the design of improved nanomaterials for targeted performance, one can use mechanistic insights into the interactions between capping agents and nanoparticles. Here, we employ density functional theory (DFT), reactive force-field molecular dynamics (ReaxFF MD) simulations, and optical spectroscopy to study the interactions between salicylate, as a model capping agent, and zinc oxide (ZnO) nanoparticles. In this work, we find that salicylate strongly interacts with the nanoparticle via the formation of O–Zn bonds in a distorted six-membered coordination ring structure. We describe the mechanisms of capping of ZnO nanoparticles by salicylate via three different binding modes. Simulations indicate that salicylate undergoes dissociative adsorption at the highly active surface Zn sites via a hydrogen-transfer process, thereby forming a tridentate configuration. The water-mediated interaction also facilitates the dissociative adsorption, leading to two salicylate O atoms coordinating with a surface Zn atom, while the other salicylate O atom bonds with another surface Zn atom. For molecular adsorption, binding free energies indicate that salicylate binds more strongly to ZnO, often in a bidentate configuration, than water does. The formation of the salicylate–ZnO complex is substantiated by UV–visible and Fourier transform infrared spectra. We find that the C=O stretching mode of salicylate becomes softened when it interacts with the nanoparticle, suggesting chemisorption of salicylate on ZnO. Although DFT predicts strong interaction between salicylate and ZnO, ReaxFF MD simulation indicates the moderate interaction between these two components in aqueous solution. Water molecules in close contact with the nanoparticle surface undergo dissociation, thus resulting in a surface hydroxyl, a reactive oxygen species that may influence the nanoparticle’s catalytic properties. Overall, the atomic-level information provided here can guide the selection process of salicylate as an appropriate agent for ZnO nanoparticle synthesis when strong interactions of particles with capping agents are required.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deconvolution of Water-Splitting on the Triple-Conducting Ruddlesden–Popper-Phase Anode for Protonic Ceramic Electrolysis Cells

Triple-conducting materials have been proved to improve the performance of popular protonic ceramic electrolysis cells. However, partially because of the complexity of the water splitting reaction involving three charge carriers, that is, oxygen (O 2– ), proton (H + ), and electron (e – ), the triple-conducting reaction mechanism was not clear, and the reaction conducting pathways have seldom been addressed. In this study, the triple conducting Ruddlesden–Popper phase Pr 1.75 Ba 0.25 NiO 4+δ as an anode on the BaCe 0.7 Zr 0.1 Y 0.1 Yb 0.1 O 3–δ electrolyte was fabricated and its electroresponses were characterized by electrochemical impedance spectroscopy with various atmospheres and temperatures. The impedance spectra are deconvoluted by means of the distribution of the relaxation time method. The surface exchange rate and chemical diffusivity of H + and O 2– are characterized by electrical conductivity relaxation. The physical locations of electrochemical processes are also identified by atomic layer deposition with a surface inhibitor. A microkinetics model is proposed toward conductivities, triple-conducting pathways, reactant dependency, surface exchange and bulk diffusion capabilities, and other relevant properties. Lastly, the rate-limiting steps and suggestions for further improvement of electrode performance are presented.

25 ENERGY STORAGE↗

A Theoretical Study of Doping Evolution of Phonons in High-Temperature Cuprate Superconductors

Hole-doped high-temperature copper oxide-based superconductors (cuprates) exhibit complex phase diagrams where electronic orders like a charge density wave (CDW) and superconductivity (SC) appear at low temperatures. The origins of these electronic orders are still open questions due to their complex interplay and correlated nature. These electronic orders can modify the phonons in the system, which has also been experimentally found in several cuprates as a softening in the phonon frequency at the CDW vector. Recent experiments have revealed that the softening in phonons in cuprates due to CDW shows intriguing behavior with increasing hole doping. Hole doping can also change the underlying Fermi surface. Therefore, it is an interesting question whether the doping-induced change in the Fermi surface can affect the softening of phonons, which in turn can reveal the nature of the electronic orders present in the system. In this work, we investigate this question by studying the softening of phonons in the presence of CDW and SC within a perturbative approach developed in an earlier work. We compare the results obtained within the working model to some experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Climate response to off-equatorial stratospheric sulfur injections in three Earth system models – Part 2: Stratospheric and free-tropospheric response

The paper constitutes Part 2 of a study performing a first systematic inter-model comparison of the atmospheric responses to stratospheric aerosol injection (SAI) at various single latitudes in the tropics, as simulated by three state-of-the-art Earth system models – CESM2-WACCM6, UKESM1.0, and GISS-E2.1-G. Building on Part 1 (Visioni et al., 2023) we demonstrate the role of biases in the climatological circulation and specific aspects of the model microphysics in driving the inter-model differences in the simulated sulfate distributions. We then characterize the simulated changes in stratospheric and free-tropospheric temperatures, ozone, water vapor, and large-scale circulation, elucidating the role of the above aspects in the surface SAI responses discussed in Part 1. We show that the differences in the aerosol spatial distribution can be explained by the significantly faster shallow branches of the Brewer–Dobson circulation in CESM2, a relatively isolated tropical pipe and older tropical age of air in UKESM, and smaller aerosol sizes and relatively stronger horizontal mixing (thus very young stratospheric age of air) in the two GISS versions used. We also find a large spread in the magnitudes of the tropical lower-stratospheric warming amongst the models, driven by microphysical, chemical, and dynamical differences. These lead to large differences in stratospheric water vapor responses, with significant increases in stratospheric water vapor under SAI in CESM2 and GISS that were largely not reproduced in UKESM. For ozone, good agreement was found in the tropical stratosphere amongst the models with more complex microphysics, with lower stratospheric ozone changes consistent with the SAI-induced modulation of the large-scale circulation and the resulting changes in transport. In contrast, we find a large inter-model spread in the Antarctic ozone responses that can largely be explained by the differences in the simulated latitudinal distributions of aerosols as well as the degree of implementation of heterogeneous halogen chemistry on sulfate in the models. The use of GISS runs with bulk microphysics demonstrates the importance of more detailed treatment of aerosol processes, with contrastingly different stratospheric SAI responses to the models using the two-moment aerosol treatment; however, some problems in halogen chemistry in GISS are also identified that require further attention. Overall, our results contribute to an increased understanding of the underlying physical mechanisms as well as identifying and narrowing the uncertainty in model projections of climate impacts from SAI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Correlations between thermal history and keyhole porosity in laser powder bed fusion

Additive manufacturing has the potential to revolutionize the production of metallic components as it yields near net shape parts with complex geometries and minimizes waste. At the present day, additively manufactured components face qualification and certification challenges due to the difficulty in controlling defects. This has driven a significant research effort aimed at better understanding and improving processing controls – yielding a plethora of in-situ measurements aimed at correlating defects with material quality metrics of interest. In this work, we develop machine-learning methods to learn correlations between thermal history and subsurface porosity for a variety of print conditions in laser powder bed fusion. Un-normalized surface temperatures (in the form of black-body radiances) are obtained using high-speed infrared imaging and porosity formation is observed in the sample cross-section through synchrotron x-ray imaging. To demonstrate the predictive power of these features, we present four statistical machine-learning models that correlate temperature histories to subsurface porosity formation in laser fused Ti-6Al4V powder.

36 MATERIALS SCIENCE↗

Fresnel diffraction imaging of surface nanostructure using coherent resonant x-ray scattering

We investigated surface nanostructures on an antiferromagnet MnBi 2 Te 4 using a novel imaging technique, direct (real)-space and real time coherent x-ray imaging (direct-CXI). This technique has provided new insights into antiferromagnetic textures, including the formation of anti-phase antiferromagnetic (AFM) domains and thermal dynamics of AFM domains and domain walls. While this method produces real-space images of AFM textures without requiring a complex imaging retrieval process, its underlying imaging mechanism has not been fully understood, limiting a deep understanding of AFM textures and the information they contain. By investigating the well-defined structural characteristics of the nanostructures fabricated on MnBi 2 Te 4 ⁠, we elucidate the imaging principle of this novel technique. We find that the observed images can be well explained by the Fresnel diffraction integral. Using a simple model from classical optics, our calculations successfully reproduce the experimentally observed images of the nanostructures. This demonstrates that direct-CXI not only provides straightforward real-space imaging but also contains phase information through its Fresnel diffraction integral.

36 MATERIALS SCIENCE↗

A Micromorphic Length-Scale Coupling Framework for the Determination of Higher-Order Constitutive Models and the Multi-Scale Simulation of Heterogeneous Materials [Thesis]

Heterogeneous materials and materials with complex microstructures pose a unique challenge in the development of accurate models of their response to external stimuli. These difficulties principally arise due to the difficulty in characterizing and modeling the constituents and their interactions. It is usually possible, though non-trivial, to construct an explicit representation of the microstructure (a direct numeric simulation or DNS) but the method by which the complex modes of deformation and other processes can be homogenized to a reduced-order approximation is frequently unclear. Many homogenization approaches are ad-hoc and lack a strong justification beyond ease of computation. Furthermore, the homogenization approach can, in some cases, not utilize the full breadth of information available in the computation of the macro-scale stresses and deformation measures. It is also noteworthy that homogenization, by its very nature, will tend to obfuscate details of processes occurring at the lower length-scale. It is therefore of interest to include as much information as possible in the construction of the reduced-order model so as to be predictive in a variety of loading environments. We here present a length-scale bridging technique based upon the micromorphic continuum mechanics of Eringen which incorporates volume and surface area averages as a part of its construction. This approach enforces the balance equations at the micro-scale and then studies the effect of the spatially varying nature on the macro-scale. This leads, naturally, to further balance equations which are solved at the macro-scale. This work details a homogenization framework which arises naturally from the micromorphic construction of Eringen attempting to introduce no definitions beyond which are justifiable from micro-structural considerations. One of the results of this effort is the so-called “micromorphic filter” which has been developed and applied to several DNS to demonstrate its effectively. In order to determine the macroscopic degrees of freedom we utilize the special case of an overlap coupling technique where the macro-scale is fully constrained to the micro-scale. This enables us to study the resulting material properties as expected but also allows us to further study the boundary conditions on the additional degrees of freedom at the macro-scale.

42 ENGINEERING↗

Fundamental limitations on photoisomerization from thermodynamic resource theories

Small, out-of-equilibrium, and quantum systems defy simple thermodynamic expressions. Such systems are exemplified by molecular switches, which exchange heat with a bath. These molecules can photoisomerize, or change conformation, or switch, on absorbing light. The photoisomerization probability depends on kinetic details that couple the molecule’s energetics to its dissipation. Therefore, a simple, general, thermodynamic-style bound on the photoisomerization probability seems out of reach. We derive such a bound using a resource theory. The resource-theory framework is a set of mathematical tools, developed in quantum information theory, used to generalize thermodynamics to small and quantum settings. From this toolkit has been derived a generalization of the second law, the thermomajorization preorder. We use thermomajorization to upper-bound the photoisomerization probability. Then, we compare the bound with an equilibrium prediction and with a Lindbladian model. We identify a realistic parameter regime in which the Lindbladian evolution saturates the thermomajorization bound. We also quantify the energy coherence in the electronic degree of freedom, and we argue that this coherence cannot promote photoisomerization. Overall, this work illustrates how quantum-information-theoretic thermodynamics can elucidate complex quantum processes in nature, experiments, and synthetics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Roles of Carburization and Oxidation on Transition-Metal Catalyzed Surface Reactions

Many industrial chemical transformation processes are based on reactions catalyzed by transition metals. In some cases the nano-size particles used as catalysts are synthesized simultaneously or right before the initiation of the main catalytic reaction. Synthesis treatments involve complex processes such as reductions followed by calcinations in oxygenated environments coupled to the presence of gases used for the main reaction. Synthesis details are crucial because they have important consequences on the catalyst shape, size, and chemical composition, and may also affect the structure and composition of the catalyst support. All these effects may alter the main catalyzed reaction in unexpected ways. This work focuses on the synthesis of the overall catalytic system that includes catalyst and support taking place in parallel with the initial stages of the main catalyzed reaction. Two important processes are characterized: carburization and oxidation of the transition metal catalyst during synthesis or initial stages of the main reaction. The objective is to understand how the coking mechanism coupled to surface oxidation and changes in the substrate influence the dry reforming of methane taken as a model reaction. The approach uses first principles computational methods to investigate the catalyst synthesis and characterization and catalytic reaction mechanisms. These studies are complemented with experimental characterization done by the group of Dr. Renu Sharma at the National Institute of Standards and Technology (unfunded collaborator) using high resolution transmission electron microscopy. Specific goals are to identify the role of the synthesis process on the initial nanocatalyst structure and composition, and to determine how that structure and composition may evolve after being exposed to the reactants atmosphere. These studies facilitate a systematic analysis of conditions where the nanoparticle composition can be optimized regarding activity and stability. Moreover, the fundamental study of catalyst/support interactions yields valuable insights for other catalytic processes and is a first step toward controlled nanocatalyst synthesis.

36 MATERIALS SCIENCE↗

Model Intercomparison of the Impacts of Varying Cloud Droplet–Nucleating Aerosols on the Life Cycle and Microphysics of Isolated Deep Convection

The microphysical impacts of aerosol particles on scattered isolated deep convective cells near Houston, Texas, on 19 June 2013, are examined using multiple cloud-resolving model (CRM) simulations initialized with vertical profiles of low and high concentrations of cloud droplet–nucleating aerosols. These simulations formed part of the Model Intercomparison Project (MIP) conducted by the Deep Convective Working Group of the Aerosol, Cloud, Precipitation and Climate (ACPC) initiative. Each CRM generated a field of convective cells representing those observed during the case study with varying degrees of accuracy. The Tracking and Object-Based Analysis of Clouds (tobac) cell-tracking algorithm was applied to each MIP CRM simulation to track relatively long-lived convective cells (20–60 min). Most of the CRMs produced similar aerosol loading impacts on the warm phase of tracked cell properties with reduced autoconversion and accretion growth of rain, increased cloud water, reduced rainfall, and reduced near-surface evaporation of rain. The sign of aerosol impacts on the warm-phase properties of the convective cells was also quite consistent over cell lifetimes with the greatest magnitude of influence in the first half of the life cycle in most CRMs. In contrast, the ice-phase response to aerosol loading was highly variable among CRMs and included increases or decreases in ice amounts at inconsistent stages of the cell life cycle and midlevel versus upper-level changes in ice. This intermodel variability in ice is indicative both of the complex indirect interactions between aerosols and ice-phase processes in deep convection and their associated parameterizations.

Aerosol-cloud interaction↗

Charges on a suspended silicon nitride membrane under a high-energy electron beam

Thin silicon nitride (SiN x ) membranes are widely used in gas and liquid phase transmission electron microscopy (TEM) and as phase plates to enhance imaging contrast. SiN x contains trap sites for both positive and negative charges, which can be manipulated by high-energy electron irradiation, external potential biasing, or light exposure. Charge accumulation on the membrane can significantly affect in situ TEM processes, including chemical and electrochemical reactions, nanoparticle dynamics, and catalytic activity, or introduce unwanted phase shifts when used as a phase plate. Here, in this study, charge accumulation on suspended SiN x membranes was investigated using off-axis electron holography combined with model-free charge analysis, supported by custom finite element analysis (FEA) simulations. An average residual positive charge density of approximately 2.8 × 10 −4 C m −2 was measured. Localized and stable regions of both positive and negative charges were identified on the membrane. The global positive and localized positive/negative charges give rise to strong electric fields and electroosmotic slip velocities at the membrane surface, which are sufficient to induce non-Brownian particle behavior and directional fluid flow, offering a physical explanation for previously observed anomalies in particle dynamics, nucleation, and growth during gas and liquid phase TEM experiments. These results provide a benchmark for understanding charge behavior at SiNx interfaces in gas and liquid phase TEM. Furthermore, the FEA simulations establish a framework for future investigations into charge distribution, electrostatic potentials, and electrical double layers at solid–liquid interfaces, particularly in complex geometries and chemically dynamic environments.

Suspended silicon↗

LLNL FY23 Aging and Lifetimes Exit Criteria: Milestone 8650, GC #4

Silicone rubbers are important materials within the complex and are used in various applications due to wide operating temperatures, favorable aging characteristics, and chemical inertness. However, silicone rubbers are susceptible to degradation in the presence of acid or base. The goal of the project is to design outgassing vessels to probe material degradation byproducts from additively manufactured (AM) silicones compressed against Arrhenius base under dry (sub 100 ppm moisture) environments. The overarching goal of this task is to develop a mechanistic understanding of degradation from experimental data that will be used to inform predictive lifetime and aging models. An outgassing chamber for multi-material compressive aging was designed. The chamber consists of a stainless-steel enclosure, which is hermetically sealed once AM silicone samples are compressed against Arrhenius base pellets in compression rigs (Fig. 1). Outgassing chambers, compression rigs, and accessories were procured, arrived onsite, and received surface treatments. Outgassing of degradation byproducts will be monitored over 12 weeks at 35 °C and 70 °C.

36 MATERIALS SCIENCE↗

Price equation captures the role of drug interactions and collateral effects in the evolution of multidrug resistance

Bacterial adaptation to antibiotic combinations depends on the joint inhibitory effects of the two drugs (drug interaction [DI]) and how resistance to one drug impacts resistance to the other (collateral effects [CE]). Here we model these evolutionary dynamics on two-dimensional phenotype spaces that leverage scaling relations between the drug-response surfaces of drug-sensitive (ancestral) and drug-resistant (mutant) populations. We show that evolved resistance to the component drugs – and in turn, the adaptation of growth rate – is governed by a Price equation whose covariance terms encode geometric features of both the two-drug-response surface (DI) in ancestral cells and the correlations between resistance levels to those drugs (CE). Within this framework, mean evolutionary trajectories reduce to a type of weighted gradient dynamics, with the drug interaction dictating the shape of the underlying landscape and the collateral effects constraining the motion on those landscapes. We also demonstrate how constraints on available mutational pathways can be incorporated into the framework, adding a third key driver of evolution. Our results clarify the complex relationship between drug interactions and collateral effects in multidrug environments and illustrate how specific dosage combinations can shift the weighting of these two effects, leading to different and temporally explicit selective outcomes.

Gjini, Erida↗

Street-level temperature estimation using graph neural networks: Performance, feature embedding and interpretability

Estimating street-level air temperature is a challenging task due to the highly heterogeneous urban surfaces, canyon-like street morphology, and the diverse physical processes in the built environment. Though pioneering studies have embarked on investigations via data-driven approaches, many questions remain to be answered. Here, in this study, we leveraged an innovative framework and redefined the street-level temperature estimation problem using Graph Neural Networks (GNN) with spatial embedding techniques. The results showed that GNN models are more capable and consistent of estimating street-level temperature among tested locations, benefiting from its unique strength in handling extensive data over unstructured graph topology. In addition, we conducted in-depth analysis of feature importance to enhance the model interpretability. Among the urban features analyzed in this study, the time-variant canopy density and meter-level land use data emerge as crucial factors. Our findings highlight GNN 's high potential in capturing the complex dynamics between urban elements and their impacts on microclimate, thus offering valuable insights for comprehensive urban data collection and urban climate modeling in general. Collectively, this study also contributes to urban planning and policy by providing avenues to enhance city resilience against climate change, thereby advancing the agenda for environmental stewardship and urban sustainability.

54 ENVIRONMENTAL SCIENCES↗

Antarctic permafrost processes and antiphase dynamics of cold-based glaciers in the McMurdo Dry Valleys inferred from 10 Be and 26 Al cosmogenic nuclides

Abstract. Soil and sediment mixing and associated permafrost processes are not widely studied or understood in the McMurdo Dry Valleys of Antarctica. In this study, we investigate the stability and depositional history of near-surface permafrost sediments to ∼ 3 m depth in the Pearse and lower Wright valleys using measured cosmogenic 10Be and 26Al depth profiles. In Pearse Valley, we estimate a minimum depositional age of ∼ 74 ka for the active layer and paleoactive-layer sediments (< 0.65 m). Combined depth profile modelling of 10Be and 26Al gives a depositional age for near-surface (< 1.65 m) permafrost in Pearse Valley of 180 +20/-40 ka, implying that the deposition of permafrost sediments predates MIS 5 advances of Taylor Glacier. Deeper permafrost sediments (> 2.09 m) in Pearse Valley are thus inferred to have a depositional age of > 180 ka. At a coastal, lower-elevation site in neighbouring lower Wright Valley, 10Be and 26Al depth profiles from a second permafrost core exhibit near-constant concentrations with depth and indicate the sediments are either vertically mixed after deposition or sufficiently young so that post-depositional nuclide production is negligible relative to inheritance. 26Al/10Be concentration ratios for both depth profiles range between 4.0 and 5.2 and are all lower than the nominal surface production rate ratio of 6.75, indicating that prior to deposition, these sediments experienced complex, yet similar, exposure–burial histories. Assuming a single-cycle exposure–burial scenario, the observed 26Al/10Be ratios are equivalent to a total minimum exposure–burial history of ∼ 1.2 Myr. In proximity to the depth profile core site, we measured cosmogenic 10Be and 26Al in three granite cobbles from thin, patchy drift (Taylor 2 Drift) in Pearse Valley to constrain the timing of retreat of Taylor Glacier. Assuming simple continuous exposure, our minimum, zero-erosion exposure ages suggest Taylor Glacier partially retreated from Pearse Valley no later than 65–74 ka. The timing of retreat after 65 ka and until the Last Glacial Maximum (LGM) when Taylor Glacier was at a minimum position remains unresolved. The surface cobble ages and permafrost processes reveal Taylor Glacier advances during MIS 5 were non-erosive or mildly erosive, preserving the underlying permafrost sediments and peppering boulders and cobbles upon an older, relict surface. Our results are consistent with U/Th ages from central Taylor Valley and suggest changes in moisture delivery over Taylor Dome during MIS 5e, 5c, and 5a appear to be associated with the extent of the Ross Ice Shelf and sea ice in the Ross Sea. These data provide further evidence of antiphase behaviour through retreat of a peripheral lobe of Taylor Glacier in Pearse Valley, a region that was glaciated during MIS 5. We suggest a causal relationship of cold-based glacier advance and retreat that is controlled by an increase in moisture availability during retreat of sea ice and perhaps the Ross Ice Shelf, as well as, conversely, a decrease during times of sea ice and Ross Ice Shelf expansion in the Ross Sea.

58 GEOSCIENCES↗