Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “process coupling”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Isotopic fractionation of hydrogen in planetary exospheres due to ionosphere-exosphere coupling - Implications for Venus

The paper considers isotopic hydrogen fractionation processes in the Venusian exosphere due to ionosphere-exosphere coupling by addressing two deficiencies in the present theory of differential escape. First, a set of D/H isotopic fractionation curves is derived for the ion-neutral interactions of charge and collisional momentum transfer, and these are compared with the results of Gurwell and Yung (1993) for hot O collisional ejection. Then, the question of the relative importance of collisional ejection in atmospheric escape is reexamined using two simple exosphere models. It is shown that O-O collisions suppress the high energy component of the hot O distribution by more than a factor of 10. Moreover, the ballistic trajectories of fast O atoms that reach the nighttime reservoir of exospheric hydrogen favor downward scatter of D and H rather than their escape. It is concluded that, due to severe limits placed on the effectiveness of collisional ejection, the differential escape of D and H from Venus is determined by charge exchange interactions rather than the collisional ejection.

Hodges, R. R., Jr.↗
DOE OSTI2026

Observed Land Surface Influence on Atmospheric Heat and Moisture Profiles During Interstorms

Land-atmospheric (L-A) feedbacks have historically been studied using models whose structure and parameterizations influence outcomes and insights. The representation of L-A feedbacks based on observations alone remains an ongoing challenge for understanding boundary layer development and precipitation. To address this gap, we use ground-based passive remote sensing and in-situ observations to present an analysis of the atmosphere during 103 interstorm soil moisture drydown events spanning nine warm seasons (2016–2024) in the U.S. Southern Great Plains region. By separating events based on local L-A coupling signals and characterizing the profiles of atmospheric heat and moisture to surface energy flux behavior, we investigate the physical mechanisms linking land surface processes to boundary layer development. We find that during interstorm drydowns, the atmospheric column follows a consistent pattern: moisture increases within the boundary layer, peaks near its top, and declines rapidly above, while warming occurs through the depth. Drydowns that shift toward evaporation produce stronger and deeper thermodynamic responses than cases dominated by sensible heating, which are weaker and shallower. Additionally, moisture is accumulated faster within the boundary layer during shorter drydowns, with longer drydowns representing slower, moisture-limited growth. Drydowns with wetter initial soil moisture will sustain stronger moistening within and above the boundary layer, accelerating buoyancy growth and convective potential toward the next storm. These results provide observational evidence linking surface flux evolution to boundary layer thermodynamics and offer a process-level benchmark for evaluating coupled L-A representations in models and demonstrating the influence of soil moisture on short-term weather forecasting skill.

Zhang, M. S. [Massachusetts Inst. of Technology (M↗
NASA NTRS1984

Dynamics of coupled ice-ocean system in the marginal ice zone: Study of the mesoscale processes and of constitutive equations for sea ice

This study is aimed at the modelling of mesoscale processed such as up/downwelling and ice edge eddies in the marginal ice zones. A 2-dimensional coupled ice-ocean model is used for the study. The ice model is coupled to the reduced gravity ocean model (f-plane) through interfacial stresses. The constitutive equations of the sea ice are formulated on the basis of the Reiner-Rivlin theory. The internal ice stresses are important only at high ice concentrations (90-100%), otherwise the ice motion is essentially free drift, where the air-ice stress is balanced by the ice-water stress. The model was tested by studying the upwelling dynamics. Winds parallel to the ice edge with the ice on the right produce upwilling because the air-ice momentum flux is much greater that air-ocean momentum flux, and thus the Ekman transport is bigger under the ice than in the open water. The upwelling simulation was extended to include temporally varying forcing, which was chosen to vary sinusoidally with a 4 day period. This forcing resembles successive cyclone passings. In the model with a thin oceanic upper layer, ice bands were formed.

Hakkinen, S.↗
DOE OSTI2026

Highly Anisotropic Quasi‐Direct Organic Metal Halide Hybrids: A Platform for Polarization‐Sensitive Optoelectronics

Low-dimensional organic–inorganic metal halide hybrids (OMHHs) exhibit remarkable optical properties and enhanced environmental stability. We investigate a 1D OMHH with formula C 4 N 2 H 14 PbBr 4 , consisting of Pb–Br chains separated by organic cations, which shows a large Stokes shift (0.83 eV) and broadband emission. Through first-principles calculations and polarized Raman spectroscopy, we characterize the material's vibrational properties and identify the specific phonon modes that drive exciton self-trapping. Our novel GW/Bethe-Salpeter equation force formalism reveals that low-frequency phonons (∼ 100 cm −1 , primarily involving Pb–Br motions) couple strongly with excitons, with a remarkably high Huang-Rhys factor of 137 ± 4, and gives a pathway for ultrafast structural analysis during the absorption process. This phonon-exciton coupling mechanism explains the material's broadband emission and provides a pathway for controlling optical properties through vibrations and for tuning vibrations through optical excitations. The material also exhibits highly anisotropic optical properties and electronic transport, with bands that are dispersive along the Pb–Br chains but nearly flat in perpendicular directions, resulting in direction-dependent electrical conductivity that is calculated to be an order of magnitude higher along the chain direction and consistent with measurements. These combined properties make this system an excellent platform for polarization-sensitive optoelectronic devices.

36 MATERIALS SCIENCE↗
NASA NTRS1984

Review of solar fuel-producing quantum conversion processes

The status and potential of fuel-producing solar photochemical processes are discussed. Research focused on splitting water to produce dihydrogen and is at a relatively early stage of development. Current emphasis is primarily directed toward understanding the basic chemistry underlying such quantum conversion processes. Theoretical analyses by various investigators predict a limiting thermodynamic efficiency of 31% for devices with a single photosystem operating with unfocused sunlight at 300 K. When non-idealities are included, it appears unlikely that actual devices will have efficiencies greater than 12 to 15%. Observed efficiencies are well below theoretical limits. Cyclic homogeneous photochemical processes for splitting water have efficiencies considerably less than 1%. Efficiency can be significantly increased by addition of a sacrificial reagent; however, such systems are no longer cyclic and it is doubtful that they would be economical on a commercial scale. The observed efficiencies for photoelectrochemical processes are also low but such systems appear more promising than homogeneous photochemical systems. Operating and systems options, including operation at elevated temperature and hybrid and coupled quantum-thermal conversion processes, are also considered.

Peterson, D. B.↗
NASA NTRS1990

Recent trends in stratospheric ozone - Implications of ozone and temperature correlation

The relation between ozone and temperature in the upper stratosphere is used to estimate the drift in the SBUV spectrometer on Nimbus-7 and the linear trend in ozone from 1979 to 1986. The ozone and temperature variabilities are coupled through chemical and dynamical processes. In the upper stratosphere they are anticorrelated through temperature-dependent loss rates, while in the lower stratosphere their variabilities are positively correlated through dynamical coupling. These relations are used to determine the ozone sensitivity to temperature and ozone trends in terms of temperature trends. It is shown that in the upper stratosphere (approximately 2 mbar), the eight-year trend in ozone is a manifestation of dynamical perturbations. It varies from about +4 percent at the equator to about -1 to -2 percent at high latitudes in both hemispheres.

Chandra, S.↗
DOE OSTI2021

Direct Optical Lithography of CsPbX 3 Nanocrystals via Photoinduced Ligand Cleavage with Postpatterning Chemical Modification and Electronic Coupling

Microscale patterning of solution-processed nanomaterials is important for integration in functional devices. Colloidal lead halide perovskite (LHP) nanocrystals (NCs) can be particularly challenging to pattern due to their incompatibility with polar solvents and lability of surface ligands. Here, we introduce a direct photopatterning approach for LHP NCs through the binding and subsequent cleavage of a photosensitive oxime sulfonate ester (-C=N-OSOO-). The photosensitizer binds to the NCs through its sulfonate group and is cleaved at the N-O bond during photoirradiation with 405 nm light. This bond cleavage decreases the solubility of the NCs, which allows patterns to emerge upon development with toluene. Postpatterning ligand exchange results in photoluminescence quantum yields of up to 79%, while anion exchange provides tunability in the emission wavelength. Finally, the patterned NC films show photoconductive behavior, demonstrating that good electrical contact between the NCs can be established.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Modular Electrochemical Synthesis Using a Redox Reservoir Paired with Independent Half-Reactions

Electrochemical synthesis of high-value chemicals using renewable electricity could lead to sustainability advantages over conventional chemical manufacturing. The tight coupling between paired anodic and cathodic reactions in conventional electrochemical processes, however, constrains the products and their efficiencies. Here, we demonstrate a strategy to achieve modular electrochemical synthesis (ModES) of different chemicals by decoupling and pairing the oxidative and reductive half-reactions with a heterogeneous redox reservoir (RR). Specifically, sodium nickel hexacyanoferrate (NaNiHCF) can serve as a robust and reversible RR that enables the modular production of several strong oxidants, hydrogen peroxide (H 2 O2) and sodium persulfate (Na 2 S 2 O 8 ) or active chlorine, in separate RR-coupled electrochemical processes. The system’s excellent performance and operational stability is demonstrated by continuous operation over 100 cycles (~36 h) to produce H 2 O 2 and Na 2 S 2 O 8 with a high voltage efficiency. The concept presented here highlights the prospects for using RRs to enable efficient, modular, and on-demand electrochemical synthesis of valuable chemicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Evaluating the Influence of Plants on Hydrologic Cycling: Quantifying and Validating the Role of Plant Processes and Stomatal Conductance (Final Report)

Plants can exert strong controls on water cycling, both locally and across the globe. The impacts are driven by how plants control the exchange of water and energy between the land and the atmosphere, which in turn are influenced by how plants photosynthesize and grow (e.g. biogeochemical cycling). Our team has shown previously that plant responses to increasing CO 2 in the atmosphere can influence the water cycle through changes in rainfall and the amount of water on land with implications for drought. Plant responses can also alter both average and extreme flow of water in rivers, with implications for freshwater availability and flood frequency. Although the impact of plant responses on water cycling has been demonstrated to exist, significant uncertainty remains in our understanding of the magnitude and form of the plant responses themselves, as well as their ultimate impact on water availability for people and ecosystems. In this project we worked to quantify the role of plant processes in regulating water cycling (e.g., precipitation, evapotranspiration, runoff, droughts, etc.). We addressed the following objectives: quantify the control of physiological vs. radiative effects on water cycling across many coupled models in idealized and realistic scenarios, as well as to quantify how our assumptions about leaf-level processes (e.g. coupling between stomatal conductance and photosynthesis), and organism-level processes (e.g. leaf area response to high CO 2 ) contributed to uncertainties and biases in water cycling. We focused on the impact of these processes on water cycling to assess the sensitivities of precipitation, evapotranspiration, and runoff to model assumptions about plant processes at the leaf, organism, and community level. As part of this project we used measurements of stable carbon isotopes to generate observationally based estimates of plant functioning during the historical period. This work is directly relevant to DOE’s primary scientific research questions in the RGMA topic area (b) on Biogeochemical processes, feedbacks, and interactions in the Earth system.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Compaction of crushed salt for safe containment – a summary of the KOMPASS projects

Abstract. For the underground disposal of high-level nuclear waste in rock salt formations, the safety concept includes the backfilling of open cavities with crushed salt. For the prognosis of the sealing function of the backfill for the safe containment of the radioactive waste, it is crucial to have a comprehensive process understanding of the crushed-salt compaction behavior. The crushed-salt compaction process is influenced by internal properties (e.g., grain size, mineralogy, and moisture content) and boundary conditions (e.g., temperature, stress state, and compaction rate) and, therefore, involves several coupled thermal–hydro–mechanical (THM) processes (Hansen et al., 2014; Kröhn et al., 2017). With the paradigm shift from the limited release of radionuclides to safe containment due to the German Repository Site Selection Act passed in 2017, the importance of crushed salt as geotechnical barrier has increased, with a focus on the evolution of its hydraulic properties. Based on the knowledge gaps in the current process understanding, the “Compaction of crushed salt for safe containment” (KOMPASS) projects were initiated to improve the scientific basis behind using crushed salt for the long-term isolation of high-level nuclear waste within rock salt repositories. The efforts to improve the prediction of crushed-salt compaction begun during the first phase of the KOMPASS projects (Czaikowski et al., 2020) and were followed up in a second phase ending in June 2023. The primary achievements of the projects are as follows (Czaikowski et al., 2020; Friedenberg et al., 2022): specification of the KOMPASS reference material, an easily available and reproducible synthetic crushed-salt material, for generic investigations; development of pre-compaction methods and successful production of samples in the short term and under in situ loading conditions; formulation of an extended laboratory program addressing the isolated investigation of known relevant factors influencing the compaction behavior of crushed salt (Düsterloh et al., 2022); execution of long-term compaction tests addressing isotropic and deviatoric load changes, temperature, and compaction state; construction of a backfill body using the KOMPASS reference material in the Sondershausen mine through collaboration with the SAVER (Entwicklung eines salzgrusbasierten Versatzkonzepts unter der Option Rückholbarkeit) project (Schaarschmidt and Friedenberg, 2022); advancement of the tools for microstructure investigation methods (Svensson and Laurich, 2022); generation (first stages) of a microphysical process list combining literature research with our own findings; benchmarking of long-term compaction test for model development and optimization of various existing models as well as the development of new models; application of a virtual demonstrator (2D model representing a backfilled drift in rock salt) for the visualization of developments and the quantification of the models (Rabbel, 2022). In summary, the KOMPASS projects contributed to the reduction of uncertainties and the strengthening of the safety case for using crushed salt within rock salt repositories.

Friedenberg, Larissa↗
NASA NTRS1996

Aeroelastic Analysis of Modern Complex Wings

A process is presented by which aeroelastic analysis is performed by using an advanced computational fluid dynamics (CFD) code coupled with an advanced computational structural dynamics (CSD) code. The process is demonstrated on an F/A-18 Stabilator using NASTD (an in-house McDonnell Douglas Aerospace East CFD code) coupled with NASTRAN. The process is also demonstrated on an aeroelastic research wing (ARW-2) using ENSAERO (an in-house NASA Ames Research Center CFD code) coupled with a finite element wing-box structures code. Good results have been obtained for the F/A-18 Stabilator while results for the ARW-2 supercritical wing are still being obtained.

Kapania, Rakesh K.↗
DOE OSTI2020

Mass production of high-performance single atomic FeNC electrocatalysts via sequenced ultrasonic atomization and pyrolysis process

Mass production of highly efficient, durable, and inexpensive single atomic catalysts is currently the major challenge associated with the oxygen reduction reaction (ORR) for fuel cells. In this study, we develop a general strategy that uses a simple ultrasonic atomization coupling with pyrolysis and calcination process to synthesize single atomic FeNC catalysts (FeNC SACs) at large scale. The microstructure characterizations confirm that the active centers root in the single atomic Fe sites chelating to the four-fold pyridinic N atoms. The identified specific Fe active sites with the variable valence states facilitate the transfer of electrons, endowing the FeNC SACs with excellent electrochemical ORR activity. The FeNC SACs were used as cathode catalysts in a homemade Zn-air battery, giving an open-circuit voltage (OCV) of 1.43 V, which is substantially higher than that of commercial Pt/C catalysts. This study provides a simple approach to the synthesis of single atomic catalysts at large scale.

36 MATERIALS SCIENCE↗
DOE OSTI2022

TDycores-Project/TDycore

The terrestrial dycores will solve tightly coupled conservation of water and energy equations. While the land dycore will simulate processes for a three-dimensional subsurface coupled to a one-dimensional plant’s vascular system, the river dycore will simulate two-dimensional processes in the river.

Bisht, Gautam↗
NASA NTRS2012

Validation of Atmosphere/Ionosphere Signals Associated with Major Earthquakes by Multi-Instrument Space-Borne and Ground Observations

The latest catastrophic earthquake in Japan (March 2011) has renewed interest in the important question of the existence of pre-earthquake anomalous signals related to strong earthquakes. Recent studies have shown that there were precursory atmospheric/ionospheric signals observed in space associated with major earthquakes. The critical question, still widely debated in the scientific community, is whether such ionospheric/atmospheric signals systematically precede large earthquakes. To address this problem we have started to investigate anomalous ionospheric / atmospheric signals occurring prior to large earthquakes. We are studying the Earth's atmospheric electromagnetic environment by developing a multisensor model for monitoring the signals related to active tectonic faulting and earthquake processes. The integrated satellite and terrestrial framework (ISTF) is our method for validation and is based on a joint analysis of several physical and environmental parameters (thermal infrared radiation, electron concentration in the ionosphere, lineament analysis, radon/ion activities, air temperature and seismicity) that were found to be associated with earthquakes. A physical link between these parameters and earthquake processes has been provided by the recent version of Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) model. Our experimental measurements have supported the new theoretical estimates of LAIC hypothesis for an increase in the surface latent heat flux, integrated variability of outgoing long wave radiation (OLR) and anomalous variations of the total electron content (TEC) registered over the epicenters. Some of the major earthquakes are accompanied by an intensification of gas migration to the surface, thermodynamic and hydrodynamic processes of transformation of latent heat into thermal energy and with vertical transport of charged aerosols in the lower atmosphere. These processes lead to the generation of external electric currents in specific regions of the atmosphere and the modifications, by dc electric fields, in the ionosphere-atmosphere electric circuit. We retrospectively analyzed temporal and spatial variations of four different physical parameters (gas/radon counting rate, lineaments change, long-wave radiation transitions and ionospheric electron density/plasma variations) characterizing the state of the lithosphere/atmosphere coupling several days before the onset of the earthquakes. Validation processes consist in two phases: A. Case studies for seven recent major earthquakes: Japan (M9.0, 2011), China (M7.9, 2008), Italy (M6.3, 2009), Samoa (M7, 2009), Haiti (M7.0, 2010) and, Chile (M8.8, 2010) and B. A continuous retrospective analysis was preformed over two different regions with high seismicity- Taiwan and Japan for 2003-2009. Satellite, ground surface, and troposphere data were obtained from Terra/ASTER, Aqua/AIRS, POES and ionospheric variations from DEMETER and COSMIC-I data. Radon and GPS/TEC were obtaining from monitoring sites in Taiwan, Japan and Italy and from global ionosphere maps (GIM) respectively. Our analysis of ground and satellite data during the occurrence of 7 global earthquakes has shown the presence of anomalies in the atmosphere. Our results for Tohoku M9.0 earthquake show that on March 7th, 2011 (4 days before the main shock and 1 day before the M7.2 foreshock of March 8, 2011) a rapid increase of emitted infrared radiation was observed by the satellite data and an anomaly was developed near the epicenter. The GPS/TEC data indicate an increase and variation in electron density reaching a maximum value on March 8. From March 3 to 11 a large increase in electron concentration was recorded at all four Japanese ground-based ionosondes, which returned to normal after the main earthquake. Similar approach for analyzing atmospheric and ionospheric parameters has been applied for China (M7.9, 2008), Italy (M6.3, 2009), Samoa (M7, 2009), Haiti (M7.0, 2010) and Chile (M8.8, 2010) eahquakes. Results have revealed the presence of related variations of these parameters implying their connection with the earthquake process. The second phase (B) of this validation included 102 major earthquakes (M>5.9) in Taiwan and Japan. We have found anomalous behavior before all of these events with no false negatives. False alarm ratio for false positives is less then 10% and has been calculated for the same month of the earthquake occurrence for the entire period of analysis (2003-2009). The commonalities for detecting atmospheric/ionospheric anomalies are: i.) Regularly appearance over regions of maximum stress (i.e., along plate boundaries); ii.) Anomaly existence over land and sea; and iii) association with M>5.9 earthquakes not deeper than 100km. Due to their long duration over the same region these anomalies are not consistent with a meteorological origin. Our initial results from the ISTF validation of multi-instrument space-borne and ground observations show a systematic appearance of atmospheric anomalies near the epicentral area, one to seven (average) days prior to the largest earthquakes, and suggest that it could be explained by a coupling process between the observed physical parameters and the pre-earthquake preparation processes.

Ouzounov, Dimitar↗
DOE OSTI2024

Proton-Coupled Electron Transfer Mechanisms for CO 2 Reduction to Methanol Catalyzed by Surface-Immobilized Cobalt Phthalocyanine

Immobilized cobalt phthalocyanine (CoPc) is a highly promising architecture for the six-proton, six-electron reduction of CO 2 to methanol. This electroreduction process relies on proton-coupled electron transfer (PCET) reactions that can occur by sequential or concerted mechanisms. Immobilization on a conductive support such as carbon nanotubes or graphitic flakes can fundamentally alter the PCET mechanisms. We use density functional theory (DFT) calculations of CoPc adsorbed on an explicit graphitic surface model to investigate intermediates in the electroreduction of CO 2 to methanol. Our calculations show that the alignment of the CoPc and graphitic electronic states influences the reductive chemistry. These calculations also distinguish between charging the graphitic surface and reducing the CoPc and adsorbed intermediates as electrons are added to the system. This analysis allows us to identify the chemical transformations that are likely to be concerted PCET, defined for these systems as the mechanism in which protonation of a CO 2 reduction intermediate is accompanied by electron abstraction from the graphitic surface to the adsorbate without thermodynamically stable intermediates. Furthermore, this work establishes a mechanistic pathway for methanol production that is consistent with experimental observations and provides fundamental insight into how immobilization of the CoPc impacts its CO 2 reduction chemistry.

Adsorption↗
DOE OSTI2023

Pulsed-Potential Electrolysis Enhances Electrochemical C–N Coupling by Reorienting Interfacial Ions

The electrochemical processing of anthropogenic CO 2 is an emerging technology aimed at utilizing renewable energies to synthesize valuable chemicals. Recently, developments in broadening the scope of the CO 2 reduction reaction (CO 2 RR) by enabling heteroatom coupling have surged with a focus on C–N bond formation. Herein, we investigate the factors that govern the selectivity and activity in synthesizing urea from environmentally malignant chemical feedstocks (CO 2 and NO 3 – ). Through a combination of electrolyte optimization and pulsed potential electrolysis, electrochemical urea production was optimized to a Faradaic efficiency of 60.4% with current densities reaching as high as 310 μA cm –2 . This work was further supported by in situ surface enhanced infrared absorbance spectroscopy that reveals the formation of C–N-related species at low overpotentials. Density functional theory calculations revealed that the reaction progresses between early reduction intermediates for the CO 2 RR and NO 3 RR and offered insights into the impacts of pulsed-potential on substrate transport.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Spin selective charge recombination in chiral donor-bridge-acceptor triads

In this paper, we outline a physically motivated framework for describing spin-selective recombination processes in chiral systems, from which we derive spin-selective reaction operators for recombination reactions of donor–bridge–acceptor molecules, where the electron transfer is mediated by chirality and spin–orbit coupling. In general, the recombination process is selective only for spin-coherence between singlet and triplet states, and it is not, in general, selective for spin polarization. We find that spin polarization selectivity only arises in hopping-mediated electron transfer. We describe how this effective spin-polarization selectivity is a consequence of spin-polarization generated transiently in the intermediate state. Furthermore, the recombination process also augments the coherent spin dynamics of the charge separated state, which is found to have a significant effect on the recombination dynamics and to destroy any long-lived spin polarization. Although we only consider a simple donor–bridge–acceptor system, the framework we present here can be straightforwardly extended to describe spin-selective recombination processes in more complex systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗