Search NASA⌕ Search

SEARCH · Search NASA

Results for “LEADING ABSTRACT”

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,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Impact of ideal MHD stability limits on high-beta hybrid operation

Here, the hybrid scenario is a candidate for stationary high-fusion gain tokamak operation in ITER and DEMO. To obtain such performance, the energy confinement and the normalized pressure ${{\beta}_{N}}$ must be maximized, which requires operating near or above ideal MHD no-wall limits. New experimental findings show how these limits can affect hybrid operation. Even if hybrids are mainly limited by tearing modes, proximity to the no-wall limit leads to 3D field amplification that affects plasma profiles, e.g. rotation braking is observed in ASDEX Upgrade throughout the plasma and peaks in the core. As a result, even the small ASDEX Upgrade error fields are amplified and their effects become visible. To quantify such effects, ASDEX Upgrade measured the response to 3D fields applied by $8\times 2$ non-axisymmetric coils as ${{\beta}_{N}}$ approaches the no-wall limit. The full n = 1 response profile and poloidal structure were measured by a suite of diagnostics and compared with linear MHD simulations, revealing a characteristic feature of hybrids: the n = 1 response is due to a global, marginally-stable n = 1 kink characterized by a large m = 1, n = 1 core harmonic due to q min being just above 1. A helical core distortion of a few cm forms and affects various core quantities, including plasma rotation, electron and ion temperature, and intrinsic W density. In similar experiments, DIII-D also measured the effect of this helical core on the internal current profile, providing information useful to understanding of the physics of magnetic flux pumping, i.e. anomalous current redistribution by MHD modes that keeps ${{q}_{\text{min}}}>1$ . Thanks to flux pumping, a broad current profile is maintained in DIII-D even with large on-axis current drive, enabling fully non-inductive operation at high ${{\beta}_{N}}$ up to 3.5–4.

3D equilibrium↗
DOE OSTI2023

Why the observed spin evolution of older-than-solar-like stars might not require a dynamo mode change

ABSTRACT The spin evolution of main-sequence stars has long been of interest for basic stellar evolution, stellar ageing, stellar activity, and consequent influence on companion planets. Observations of older-than-solar late-type main-sequence stars have been interpreted to imply that a change from a dipole-dominated magnetic field to one with more prominent higher multipoles might be necessary to account for the data. The spin-down models that lead to this inference are essentially tuned to the Sun. Here, we take a different approach that considers individual stars as fixed points rather than just the Sun. We use a time-dependent theoretical model to solve for the spin evolution of low-mass main-sequence stars that includes a Parker-type wind and a time-evolving magnetic field coupled to the spin. Because the wind is exponentially sensitive to the stellar mass over radius and the coronal base temperature, the use of each observed star as a separate fixed point is more appropriate and, in turn, produces a set of solution curves that produces a solution envelope rather than a simple line. This envelope of solution curves, unlike a single line fit, is consistent with the data and does not unambiguously require a modal transition in the magnetic field to explain it.

Astronomy & Astrophysics↗
DOE OSTI2019

Bacterial Glutathione S‐Transferases (GSTs) Involved in Breaking the β‐Aryl Ether Bond of Lignin Reveal Novel Catalytic Abilities and Reaction Mechanisms within the GST Superfamily

The glutathione S‐transferase (GST) protein superfamily contains thousands of known or predicted members, all of which contain an N‐terminal thioredoxin domain that binds a glutathione (GSH) molecule and activates it for catalysis. GSTs are important in many organisms for conjugating GSH to toxic molecules, making them more soluble and easier to excrete from cells. However, GSTs are also capable of catalyzing a wide and continually expanding range of other reactions using GSH as cofactor. Sphingomonad bacteria are capable of metabolizing low molecular weight oligomers and monomers derived from lignin, a heterogeneous polymer of aromatic subunits that is a major component of plant biomass. These bacteria and their enzymes have potential applications in converting naturally abundant lignin into valuable commodities. Sphingomonads use GSTs (called β‐etherases and GSH lyases) to catalyze two consecutive steps in their pathway for breaking the β‐aryl ether bond commonly found between aromatic subunits in the lignin polymer: β‐etherases stereospecifically break the β‐aryl ether bond by replacing it with a thioether bond involving GSH, and GSH lyases remove the GSH moiety from the resulting β‐thioether GSH conjugate. We have identified a heterodimeric β‐etherase (called BaeAB) that provides insight into the diversity of GSTs capable of acting as β‐etherases and into the evolution of the GST ability to break the β‐aryl ether bond. In addition, BaeAB has properties that are unique amongst characterized GSTs in general, including the majority of its catalytic activity apparently residing in only one of its subunits, and using an active site asparagine for catalysis (whereas most characterized GSTs use an amino acid residue with a side chain hydroxyl for catalysis). We have also identified a GSH lyase (called NaGST Nu ), which has close homologues in many different organisms, though the physiological roles of these homologues are largely unknown in other organisms besides sphingomonads. We found that the two homologues from Escherichia coli were able to catalyze the GSH lyase reaction in the pathway for breaking the β‐aryl ether bond, despite the fact that E. coli is not known to be able to metabolize lignin‐derived oligomers. This suggests that GSH lyase activity may be a common feature of these homologues (and that they likely act on GSH‐conjugated compounds naturally found within their respective organisms). In addition, the reaction mechanism of NaGST Nu involves binding two GSH molecules (whereas most characterized GSTs can only bind one GSH), as well as amino acid residues located in a channel leading from the active site to the protein exterior. Altogether, besides helping to identify novel GSTs that could potentially be used to generate value from lignin, our work is helping to uncover unique features of the large and widely studied, but relatively uncharacterized, GST superfamily. Support or Funding Information This work was supported by U.S. Department of Energy (DOE) Great Lakes Bioenergy Research Center grants (DOE Office of Science BER DE‐FC02‐07ER64494 and DE‐SC0018409). This abstract is from the Experimental Biology 2019 Meeting. There is no full text article associated with this abstract published in The FASEB Journal .

Kontur, Wayne S.↗
DOE OSTI2023

Simulation of African Easterly Waves in a Global Climate Model

Abstract African easterly waves (AEWs) exert significant influence on local and downstream high-impact weather including tropical cyclone (TC) genesis over the Atlantic. Accurate representation of AEWs in climate and weather prediction models therefore is necessary for skillful predictions. In this study, we examine simulated AEWs, including their evolution, vertical structure, and linkage to tropical cyclone genesis, in the NASA Goddard Earth Observing System Model, version 5 (GEOS-5) atmospheric global climate model. Identified by the leading empirical orthogonal function mode of time-filtered precipitation, the observed westward propagating AEWs along the southern track over the Atlantic are largely captured in GEOS-5, but with a slower phase speed and significantly weaker amplitude downstream off the West Africa coast. The weak downstream development of AEWs in GEOS-5 is accompanied by much reduced TC genesis over the main development region. Further analyses suggest that the slow westward propagation and weaker AEW amplitude downstream can be ascribed to a weak African easterly jet, while overestimated negative (positive) meridional potential vorticity (PV) gradients appear to the north (south) of 10°N in GEOS-5. The greatly overestimated positive meridional PV gradient to the south of 10°N is expected to generate strong horizontal stretching in the AEW wave pattern in the model, which hinders organization of convection and its feedback to sustain the AEW development. Persistent and vigorous AEW precipitation in the Guinea Highlands of the West Africa coast could also be responsible for reduced westward propagation of AEWs in the model.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Isochronic Evolution and the Radioactive Decay of r-process Nuclei

Abstract We report on the creation and application of a novel decay network that uses the latest data from experiment and evaluation. We use the network to simulate the late-time phase of the rapid neutron capture ( r ) process. In this epoch, the bulk of nuclear reactions, such as radiative capture, have ceased, and nuclear decays are the dominant transmutation channels. We find that the decay from short-lived to long-lived species naturally leads to an isochronic evolution in which nuclei with similar half-lives are populated at the same time. We consider random perturbations along each isobaric chain to initial solar-like r -process compositions to demonstrate the isochronic nature of the late-time phase of the r -process. Our analysis shows that detailed knowledge of the final isotopic composition allows for the prediction of late-time evolution with a high degree of confidence despite uncertainties that exist in astrophysical conditions and the nuclear physics properties of the most neutron-rich nuclei. We provide the time-dependent nuclear composition in the Appendix as supplemental material.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2000

Regolith Evolved Gas Analyzer

The Regolith Evolved Gas Analyzer (REGA) is a high-temperature furnace and mass spectrometer instrument for determining the mineralogical composition and reactivity of soil samples. REGA provides key mineralogical and reactivity data that is needed to understand the soil chemistry of an asteroid, which then aids in determining in-situ which materials should be selected for return to earth. REGA is capable of conducting a number of direct soil measurements that are unique to this instrument. These experimental measurements include: (1) Mass spectrum analysis of evolved gases from soil samples as they are heated from ambient temperature to 900 C; and (2) Identification of liberated chemicals, e.g., water, oxygen, sulfur, chlorine, and fluorine. REGA would be placed on the surface of a near earth asteroid. It is an autonomous instrument that is controlled from earth but does the analysis of regolith materials automatically. The REGA instrument consists of four primary components: (1) a flight-proven mass spectrometer, (2) a high-temperature furnace, (3) a soil handling system, and (4) a microcontroller. An external arm containing a scoop or drill gathers regolith samples. A sample is placed in the inlet orifice where the finest-grained particles are sifted into a metering volume and subsequently moved into a crucible. A movable arm then places the crucible in the furnace. The furnace is closed, thereby sealing the inner volume to collect the evolved gases for analysis. Owing to the very low g forces on an asteroid compared to Mars or the moon, the sample must be moved from inlet to crucible by mechanical means rather than by gravity. As the soil sample is heated through a programmed pattern, the gases evolved at each temperature are passed through a transfer tube to the mass spectrometer for analysis and identification. Return data from the instrument will lead to new insights and discoveries including: (1) Identification of the molecular masses of all of the gases liberated from heated soil samples; (2) Identification of the asteroid soil mineralogy to aid in the selection process for returned samples; (3) Existence of oxygen in the asteroid soil and the potential for in-situ resource utilization (ISRU); and (4) Existence of water and other volatiles in the asteroid soil. Additional information is contained in the original extended abstract.

Hoffman, John H.↗
DOE OSTI2022

Not All Aggregates Are Made the Same: Distinct Structures of Solution Aggregates Drastically Modulate Assembly Pathways, Morphology, and Electronic Properties of Conjugated Polymers

Abstract Tuning structures of solution‐state aggregation and aggregation‐mediated assembly pathways of conjugated polymers is crucial for optimizing their solid‐state morphology and charge‐transport property. However, it remains challenging to unravel and control the exact structures of solution aggregates, let alone to modulate assembly pathways in a controlled fashion. Herein, aggregate structures of an isoindigo–bithiophene‐based polymer (PII‐2T) are modulated by tuning selectivity of the solvent toward the side chain versus the backbone, which leads to three distinct assembly pathways: direct crystallization from side‐chain‐associated amorphous aggregates, chiral liquid crystal (LC)‐mediated assembly from semicrystalline aggregates with side‐chain and backbone stacking, and random agglomeration from backbone‐stacked semicrystalline aggregates. Importantly, it is demonstrated that the amorphous solution aggregates, compared with semicrystalline ones, lead to significantly improved alignment and reduced paracrystalline disorder in the solid state due to direct crystallization during the meniscus‐guided coating process. Alignment quantified by the dichroic ratio is enhanced by up to 14‐fold, and the charge‐carrier mobility increases by a maximum of 20‐fold in films printed from amorphous aggregates compared to those from semicrystalline aggregates. This work shows that by tuning the precise structure of solution aggregates, the assembly pathways and the resulting thin‐film morphology and device properties can be drastically tuned.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Density Dependence in Three Snake River Sockeye Salmon Nursery Lakes in Central Idaho

Abstract Snake River Sockeye Salmon Oncorhynchus nerka, listed as an endangered species in 1991, currently inhabit three nursery lakes (Redfish, Pettit, and Alturas lakes) in the Sawtooth Valley, Idaho. Conspecific kokanee (lacustrine Sockeye Salmon) are also present in the lakes. Snake River Sockeye Salmon recovery efforts, initially focused on genetic conservation, are now attempting to rebuild naturally spawning populations using hatchery supplementation. However, in Sockeye Salmon nursery lakes, density dependence is frequently observed when elevated O. nerka abundance leads to declines in zooplankton biomass, body size, and shifts in community composition. In turn, these changes lead to reductions in juvenile O. nerka growth rates, survival, and adult returns. We examined a long-term data set of O. nerka population metrics and associated zooplankton community metrics. We found evidence of density dependence within and among nursery lakes. We detected differences in zooplankton biomass, lengths of preferred zooplankton prey (Daphnia spp. and cyclopoid copepods), parr growth rates, and age-1 smolt size among the three lakes. We found negative relationships between O. nerka density and zooplankton biomass and size. We identified positive relationships between zooplankton biomass and two response variables: smolt size at migration and growth rates of hatchery parr. The relationships were generally similar among lakes. Variable outcomes were a result of differences in O. nerka density (or zooplankton biomass), controlled primarily by the relative proportion of spawning and rearing habitat in each lake. Understanding unique lake habitats, ecological interactions, and the role of density dependence is germane to management of Snake River Sockeye Salmon populations.

Griswold, Robert G.↗
DOE OSTI2018

Wafer integrated micro‐scale concentrating photovoltaics

Abstract A novel micro‐scale photovoltaic concept, Wafer Integrated Micro‐scale Photovoltaics (WPV), is proposed, analyzed, and experimentally demonstrated. The WPV concept seamlessly integrates multijunction micro‐cells with a multi‐functional silicon platform that simultaneously provides optical concentration, hybrid PV/CPV architecture, and mechanical alignment features. Fabrication and optical performance characterization of the Si platform are described in this paper. Over 100% improvement in the concentration‐acceptance‐angle product ( CAP ) is demonstrated using the wafer‐embedded micro‐concentrating elements, leading to significantly reduced module material and fabrication costs, sufficient angular tolerance for low‐cost trackers, and an ultra‐compact optical architecture compatible with commercial flat panel infrastructures. The development of a prototypical module with a 400× concentration ratio is described. Outdoor optical characterization of the module shows acceptance angles of ±1.7° and ±2.5° for 90% of on‐axis power and full‐width‐half‐maximum, respectively. The projected performance of the PV/CPV hybrid architecture illustrates its potential for cost‐effective collection of both direct and diffuse sunlight, thereby extending the geographic and market domains for cost‐effective PV system deployment. Leveraging low‐cost micro‐fabrication and high‐level integration techniques, the WPV approach presents a promising route to combine the high performance of multijunction solar cells and the low costs of flat‐plate Si PV systems.

Li, Duanhui↗
DOE OSTI2019

Understand the Direct Effect of CO 2 Increase on Tropical Circulation and TC Activity: Land Surface Warming Versus Direct Radiative Forcing

Abstract The direct effect of increased CO 2 involves contributions from both land warming and direct radiative forcing. Here, their relative impacts on tropical circulation and tropical cyclones (TCs) are quantified by increasing CO 2 over land and ocean separately in a high‐resolution atmosphere‐only model. It is found that land warming induces substantial vertical velocity changes over nearby oceans and such local changes are compensated by opposite motions within tropical ascending regions. The cloud mask effect leads to smaller CO 2 ‐induced radiative heating over ascending regions. Such inhomogeneity in radiative forcing dominates compensated changes in surface fluxes and gross moist stability, leading to the slowdown of the tropical overturning circulation. TC activity is tightly linked with the large‐scale ascent because of the influence on the atmospheric humidity. Regional circulation changes caused by land warming strongly suppresses TC activity over the northwest Pacific, while weakened ascent from direct radiative forcing causes an overall TC reduction.

Zhou, Wenyu↗
DOE OSTI2023

Quantifying electron cyclotron power deposition broadening in DIII-D and the potential consequences for the ITER EC system

Abstract The injection of electron cyclotron (EC) waves fulfills a number of important tasks in nuclear fusion devices for which detailed knowledge of the spatial power deposition profile is critical. This deposition profile is commonly determined using forward models such as beam or ray tracing. Recent numerical and experimental studies have shown that small-angle scattering of the EC beam as it passes through the turbulent plasma edge can cause significant broadening of the effective deposition profile, leading to considerable underestimation of the deposition width by forward methods. However, traditional inverse methods to determine the deposition profile from measurements overestimate the deposition profile width due to transport broadening. In this work, we implement three novel methods to resolve the EC power deposition profile from measurements that counteract transport broadening by simultaneously resolving transport and power deposition. We validate their assumptions and compare the results from these methods to the traditional break-in-slope method as well as to the TORAY ray-tracing code in a set of DIII-D discharges spanning five different confinement modes. We show that the four different inverse methods, novel and established, paint a consistent picture of deposition broadening. Specifically, we show that the measured power deposition profile is between 1.6 and 3.6 times wider than the TORAY profiles. Moreover, we show the considerable consequences that this level of broadening can have for ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Efficient data acquisition and training of collisional-radiative model artificial neural network surrogates through adaptive parameter space sampling

Abstract Effective plasma transport modeling of magnetically confined fusion devices relies on having an accurate understanding of the ion composition and radiative power losses of the plasma. Generally, these quantities can be obtained from solutions of a collisional-radiative (CR) model at each time step within a plasma transport simulation. However, even compact, approximate CR models can be computationally onerous to evaluate, and in-situ evaluation of these models within a larger plasma transport code can lead to a rigid bottleneck. As a way to bypass this bottleneck, we propose deploying artificial neural network (ANN) surrogates to allow rapid evaluation of the necessary plasma quantities. However, one issue with training an accurate ANN surrogate is the reliance on a sufficiently large and representative training and validation data set, which can be time-consuming to generate. In this work we explore a data-driven active learning and training routine to allow autonomous adaptive sampling of the problem parameter space to ensure a sufficiently large and meaningful set of training data is assembled for the network training. As a result, we can demonstrate approximately order-of-magnitude savings in required training data samples to produce an accurate surrogate.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Pot1 promotes telomere DNA replication via the Stn1-Ten1 complex in fission yeast

Abstract Telomeres are nucleoprotein complexes that protect the chromosome-ends from eliciting DNA repair while ensuring their complete duplication. Pot1 is a subunit of telomere capping complex that binds to the G-rich overhang and inhibits the activation of DNA damage checkpoints. In this study, we explore new functions of fission yeast Pot1 by using a pot1-1 temperature sensitive mutant. We show that pot1 inactivation impairs telomere DNA replication resulting in the accumulation of ssDNA leading to the complete loss of telomeric DNA. Recruitment of Stn1 to telomeres, an auxiliary factor of DNA lagging strand synthesis, is reduced in pot1-1 mutants and overexpression of Stn1 rescues loss of telomeres and cell viability at restrictive temperature. We propose that Pot1 plays a crucial function in telomere DNA replication by recruiting Stn1-Ten1 and Polα-primase complex to telomeres via Tpz1, thus promoting lagging-strand DNA synthesis at stalled replication forks.

Carvalho Borges, Pâmela C. (ORCID:0000000244919874↗
DOE OSTI2021

A new open‐path eddy covariance method for nitrous oxide and other trace gases that minimizes temperature corrections

Abstract Low‐power, open‐path gas sensors enable eddy covariance (EC) flux measurements in remote areas without line power. However, open‐path flux measurements are sensitive to fluctuations in air temperature, pressure, and humidity. Laser‐based, open‐path sensors with the needed sensitivity for trace gases like methane (CH 4 ) and nitrous oxide (N 2 O) are impacted by additional spectroscopic effects. Corrections for these effects, especially those related to temperature fluctuations, often exceed the flux of gases, leading to large uncertainties in the associated fluxes. For example, the density and spectroscopic corrections arising from temperature fluctuations can be one or two orders of magnitude greater than background N 2 O fluxes. Consequently, measuring background fluxes with laser‐based, open‐path sensors is extremely challenging, particularly for N 2 O and gases with similar high‐precision requirements. We demonstrate a new laser‐based, open‐path N 2 O sensor and a general approach applicable to other gases that minimizes temperature‐related corrections for EC flux measurements. The method identifies absorption lines with spectroscopic effects in the opposite direction of density effects from temperature and, thus, density and spectroscopic effects nearly cancel one another. The new open‐path N 2 O sensor was tested at a corn ( Zea mays L.) field in Southwestern Michigan, United States. The sensor had an optimal precision of 0.1 ppbv at 10 Hz and power consumption of 50 W. Field trials showed that temperature‐related corrections were 6% of density corrections, reducing EC random errors by 20‐fold compared to previously examined lines. Measured open‐path N 2 O EC fluxes showed excellent agreement with those made with static chambers ( m = 1.0 ± 0.3; r 2 = .96). More generally, we identified absorption lines for CO 2 and CH 4 flux measurements that can reduce the temperature‐related corrections by 10–100 times compared to existing open‐path sensors. The proposed method provides a new direction for future open‐path sensors, facilitating the expansion of accurate EC flux measurements.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Positioning Accuracy in a Concurrent Robot-CNC Hybrid Manufacturing System

Abstract Additive manufacturing (AM) has gained notoriety for offering advantages over traditional manufacturing methods, such as increased design complexity and flexibility. However, it has not found widespread use beyond rapid prototyping. One hindrance to the acceptance of AM processes in industry is the time and cost of fabrication per component. While metal AM by itself can be inexpensive, extra manufacturing steps in the form of subtractive manufacturing (SM) may need to be performed to reach final part tolerances, leading to hybrid additive-subtractive manufacturing (HASM) of a part, which increases time and cost. A potential area to reduce cost is through increasing the efficiency of the HASM process by conducting additive and subtractive manufacturing simultaneously. Usually, HASM is performed in a process where AM is completed in one machine or cell and transferred to another machine or cell for SM in a sequential assembly line process. This efficiency decreases part cost, but high aspect ratio parts or parts with internal geometry that require interleaved additive deposition and machining cannot be produced. One unexplored solution to simultaneous HASM that allows for interleaved operations is to operate the deposition head and machining spindle concurrently within the same machine envelope, known as concurrent HASM (CHASM). In this type of process, both AM and SM occur simultaneously on a batch of small parts or a single large part, maintaining a high efficiency without sacrificing the full range of complex geometries that AM allows for. A potential approach to the single-machine method could be to combine a robot and mill within the same envelope. A challenge to this approach, however, is control of both systems. Most machine controllers have limited external communication or, if a robot has been integrated, only offer movement of either the robot or mill at any given time. As a result, systems must pause either the AM or SM process to switch between them rather than working simultaneously. The present work investigates the positional accuracy of such a CHASM system comprised of a robotic arm and a 3-axis mill. Open-loop tests with limited communication between machines are performed on the system to verify positional error during concurrent robot-mill movements. Under certain conditions, it is demonstrated that position error can stay within 2 mm for the duration of a single layer; however, these tests show that, generally, the open-loop positioning performance of the system is inadequate for CHASM without part-specific hand-tuning of parameters. Based on these results, a set of requirements for successful robot-CNC CHASM is proposed for future integrations.

Goodwin, Jesse↗
DOE OSTI2024

Comparing matching prescriptions between pre-equilibrium and hydrodynamic models in high-energy nuclear collisions

Abstract State-of-the-art simulations of high-energy nuclear collisions rely on hybrid setups, involving in particular a pre-equilibrium stage to let the system evolve from a far-from-equilibrium initial condition towards a near-equilibrated state after which fluid dynamics can be applied meaningfully. A known issue is the mismatch between the equation of state in the fluid-dynamical evolution and the effective one in the previous stage, which leads to discontinuities at the interface between the two models. Here we introduce a new matching prescription at this interface, based on the entropy, and we compare it with the standard one relying on local energy conservation. We study the behavior of various quantities at the switching time between the models and investigate a number of final-state hadronic observables. For the latter, we show that they are not modified significantly by the choice of matching prescription, provided an appropriate normalization is chosen for the initial state. In turn, our approach reduces sizeably the ratio of bulk over thermodynamic pressure at the beginning of the fluid-dynamical stage.

Borghini, Nicolas (ORCID:0000000239063258)↗
DOE OSTI2020

Evaluating integrated surface/subsurface permafrost thermal hydrology models in ATS (v0.88) against observations from a polygonal tundra site

Numerical simulations are essential tools for understanding the complex hydrologic response of Arctic regions to a warming climate. However, strong coupling among thermal and hydrological processes on the surface and in the subsurface and the significant role that subtle variations in surface topography have in regulating flow direction and surface storage lead to significant uncertainties. Careful model evaluation against field observations is thus important to build confidence.We evaluate the integrated surface/subsurface permafrost thermal hydrology models in the Advanced Terrestrial Simulator (ATS) against field observations from polygonal tundra at the Barrow Environmental Observatory. ATS couples a multiphase, 3D representation of subsurface thermal hydrology with representations of overland nonisothermal flows, snow processes, and surface energy balance. We simulated thermal hydrology of a 3D ice-wedge polygon with geometry that is abstracted but broadly consistent with the surface microtopography at our study site. The simulations were forced by meteorological data and observed water table elevations in ice-wedge polygon troughs. With limited calibration of parameters appearing in the soil evaporation model, the 3-year simulations agreed reasonably well with snow depth, summer water table elevations in the polygon center, and high-frequency soil temperature measurements at several depths in the trough, rim, and center of the polygon. Upscaled evaporation is in good agreement with flux tower observations. The simulations were found to be sensitive to parameters in the bare soil evaporation model, snowpack, and the lateral saturated hydraulic conductivity.Timing of fall freeze-up was found to be sensitive to initial snow density, illustrating the importance of including snow aging effects. The study provides new support for an emerging class of integrated surface/subsurface permafrost simulators.

54 ENVIRONMENTAL SCIENCES↗