SEARCH · Search NASA

Results for “surface equivalent surface temperature”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

141 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Multi-metallic conduction cooled superconducting radio-frequency cavity with high thermal stability

Superconducting radio-frequency cavities are commonly used in modern particle accelerators for applied and fundamental research. Such cavities are typically made of high-purity, bulk Nb and with cooling by a liquid helium bath at a temperature of ~2 K. The size, cost and complexity of operating a particle accelerator with a liquid helium refrigerator make the current cavity technology not favorable for use in industrial-type accelerators. We have developed a multi-metallic 1.495 GHz elliptical cavity conductively cooled by a cryocooler. The cavity has a ~2 µm thick layer of Nb 3 Sn on the inner surface, exposed to the rf field, deposited on a ~3 mm thick bulk Nb shell and a bulk Cu shell, of thickness $\geqslant\! 5$ mm deposited on the outer surface by electroplating. A bolt-on Cu plate 1.27 cm thick was used to thermally connect the cavity equator to the second stage of a Gifford-McMahon cryocooler with a nominal capacity of 2 W at 4.2 K. The cavity was tested initially in liquid helium at 4.3 K and reached a peak surface magnetic field of ~36 mT with a quality factor of 2 × 10 9 . The cavity cooled by the cryocooler achieved a peak surface magnetic field of ~29 mT, equivalent to an accelerating gradient of 6.5 MV m –1 . The conduction-cooled cavity could be operated in continuous-wave with as high as 5 W dissipation in the cavity for 1 h without any thermal breakdown, because of the Cu outer layer with high thermal conductivity. This result represents a paradigm shift in the technology of superconducting accelerator cavities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Structure, dynamics, and redox reactivity of an all-purpose flavodoxin

The flavodoxin of Rhodopseudomonas palustris CGA009 (Rp9Fld) supplies highly reducing equivalents to crucial en- zymes such as hydrogenase, especially when the organism is iron-restricted. By acquiring those electrons from photodriven electron flow via the bifurcating electron transfer flavoprotein, Rp9Fld provides solar power to vital metabolic processes. To understand Rp9Fld’s ability to work with diverse partners, we solved its crystal structure. We observed the canonical fla- vodoxin (Fld) fold and features common to other long-chain Flds but not all the surface loops thought to recognize part- ner proteins. Moreover, some of the loops display alternative structures and dynamics. To advance studies of protein– protein associations and conformational consequences, we assigned the 19F NMR signals of all five tyrosines (Tyrs). Our electrochemical measurements show that incorporation of 3- 19 F-Tyr in place of Tyr has only a modest effect on Rp9Fld’s redox properties even though Tyrs flank the flavin on both sides. Meanwhile, the 19 F probes demonstrate the expected paramagnetic effect, with signals from nearby Tyrs becoming broadened beyond detection when the flavin semiquinone is formed. However, the temperature dependencies of chemical shifts and linewidths reveal dynamics affecting loops close to the flavin and regions that bind to partners in a variety of systems. These coincide with patterns of amino acid type conservation but not retention of specific residues, arguing against detailed specificity with respect to partners. We pro- pose that the loops surrounding the flavin adopt altered con- formations upon binding to partners and may even participate actively in electron transfer.

19F NMR↗

DOE OSTI2020

Material Discovery and High Throughput Exploration of Ru Based Catalysts for Low Temperature Ammonia Decomposition

High throughput experimentation has the capability to generate massive, multidimensional datasets, allowing for the discovery of novel catalytic materials. Here, we show the synthesis and catalytic screening of over 100 unique Ru-Metal-K based bimetallic catalysts for low temperature ammonia decomposition, with a Ru loading between 1–3 wt% Ru and a fixed K loading of 12 wt% K, supported on γ-Al2O3. Bimetallic catalysts containing Sc, Sr, Hf, Y, Mg, Zr, Ta, or Ca in addition to Ru were found to have excellent ammonia decomposition activity when compared to state-of-the-art catalysts in literature. Furthermore, the Ru content could be reduced to 1 wt% Ru, a factor of four decrease, with the addition of Sr, Y, Zr, or Hf, where these secondary metals have not been previously explored for ammonia decomposition. The bimetallic interactions between Ru and the secondary metal, specifically RuSrK and RuFeK, were investigated in detail to elucidate the reaction kinetics and surface properties of both high and low performing catalysts. The RuSrK catalyst had a turnover frequency of 1.78 s−1, while RuFeK had a turnover frequency of only 0.28 s−1 under identical operating conditions. Based on their apparent activation energies and number of surface sites, the RuSrK had a factor of two lower activation energy than the RuFeK, while also possessing an equivalent number of surface sites, which suggests that the Sr promotes ammonia decomposition in the presence of Ru by modifying the active sites of Ru.

McCullough, Katherine (ORCID:0000000197525740)↗

DOE OSTI2022

Dimensionless parameters for cloudy Rayleigh-Bénard convection: Supersaturation, Damköhler, and Nusselt numbers

In steady-state Rayleigh-Bénard convection, heat is transported by turbulent thermal convection from the bottom, hot surface to the top, cold surface, leading to a height-independent sensible heat flux. When water vapor is present and cloud formation occurs, there is also an additional latent heat flux. Heat transport in cloudy Rayleigh-Bénard convection depends on turbulent flow as well as the microphysical state of the clouds: specifically, whether substantial supersaturations exist and whether cloud liquid water is removed through sedimentation/precipitation. In this article we bridge between the Rayleigh-Bénard convection literature and the atmospheric literature. We express the governing equations for cloudy convection in dimensionless form, thereby explicitly identifying the governing parameters relevant to the cloudy case, including Schmidt, Damköhler, supersaturation, and sedimentation numbers. We further connect to the atmospheric literature by obtaining a Nusselt number (dimensionless heat flux) for a cloud-convection system, directly from the conservation equations for temperature and water vapor. This flux has the same form as that identified by Zhang et al. [L. Zhang, K. L. Chong, and K.-Q. Xia, J. Fluid Mech. 874, 1041 (2019)] for convection with water vapor, but is extended to the cloudy case. For equal thermal and water vapor diffusivities, the flux corresponds to the widely used atmospheric quantities equivalent temperature and moist static energy. Using large eddy simulation (LES) of an idealized cloudy Rayleigh-Bénard convection system with fixed boundary conditions, we find that the equivalent heat flux (Nusselt number) is only weakly dependent on the microphysical details of the system, such as liquid water mixing ratio and cloud droplet number concentration. Finally, from the results, we show the vertical profiles of sensible and latent heat fluxes depend on the liquid water content, whereas the equivalent heat flux remains a constant throughout the height of the chamber.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Wrought FeCrAl alloy (C26M) cladding behavior and burst under simulated loss-of-coolant accident conditions

Cladding burst experiments for FeCrAl cladding were performed in the Severe Accident Test Station facility at Oak Ridge National Laboratory. These experiments were simulated using the BISON fuel performance code to better understand the cladding plastic behavior and failure under simulated loss-of-coolant accident conditions. 3D cladding surface boundary conditions were generated using composite axial and azimuthal profiles from experiment thermocouple data. To improve the simulation analysis capabilities in BISON for cladding burst behavior, new thermal creep, plasticity, and failure stress models specific to C26M, a wrought FeCrAl alloy, were developed and implemented. Initial cladding burst results indicated a general underprediction in the failure temperature of the six cladding burst simulations versus the observed failure temperatures. Close investigation of the experiment timing versus the underlying tensile test data revealed that, compared with the tensile specimens, the cladding tubes did not experience the same long holding time at high temperatures. New tensile tests were performed at high temperatures using a temperature ramp similar to the simulated loss-of-coolant accident experiments. These new tensile curves showed an approximately 80% increase in the ultimate tensile strength of the C26M alloy, indicating that a holding time of 10 min at 700 °C and 800 °C allows annealing to change the material microstructure. Using the updated tensile properties, the burst temperatures and stresses from the simulations showed remarkable agreement with the experimental results. This study was then extended by varying the initial pressure to highlight the burst temperature difference between standard Zircaloy-4 and C26M cladding under equivalent conditions. The results show that C26M has a burst temperature that is approximately 70–130 K greater than that of Zircaloy-4. In conclusion, these modeling predictions can be further improved by collecting high-temperature tensile data for C26M beyond the temperature ranges used in this work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Mountain Basin Controls on the Snow-to-Streamflow Signal: An AIC-Weighted Multiple Linear Regression Framework

A regression-based analysis quantifies how basin characteristics modulate the snow-to-streamflow signal. First, we use the ERA5-Land reanalysis gridded product (European Centre for Medium Range Weather Forecasts reanalysis 5 -Land component) for 4,655 hydrologic unit code - 10 (HUC10) mountain basins across the western United States (US) for water years 1987–2024. Linear regressions are performed for peak snow water equivalent (SWE) and annual streamflow for each mountain basin. Models use ordinary least squares in Python’s statsmodels package. After which, an Akaike Information Criterion (AIC)–weighted ensemble multiple linear regression (MLR) framework with 47 watershed traits is used to predict the linear regression coefficient of determination (r-squared) defining the ability of peak SWE to predict annual streamflow across all mountain basin. Predictor sets are constrained to avoid multicollinearity by excluding models with variance inflation factors (VIF) greater than 5. Mountain basin traits included in the MLR include seasonal climate, topography, vegetation type and structure, and bedrock geology. Accepted models are considered if their AIC is within 2.0 of the model with the minimum AIC, or best model. To compare predictor influence across acceptable models, we computed standardized regression coefficients. To evaluate structural redundancy among models, we constructed binary inclusion vectors for each acceptable model, denoting whether a predictor was present (1) or absent (0). Core predictor variables are defined as occurring in at least 67% of the acceptable models. For this regional analysis, only one model was found acceptable, with higher snow-to-streamflow translation (higher r-squared) occurring in colder mountain basins with higher relative winter precipitation, more snow accumulation and a lower fraction of annual precipitation that falls in the spring and summer. The second component of the data package uses previously published, high-resolution output from an integrated hydrological model of the East River watershed using the U.S. Geological Survey Groundwater and Surface water Flow model (GSFLOW, doi:10.15485/1998576). East River MLR expands upon the approach described above to explore the response of five streamflow metrics—annual streamflow, runoff efficiency, 7-day minimum flow, low-flow duration, and non-perennial stream fraction to snow system indicators including peak SWE, snow-covered area, snow disappearance date, and the fraction of basin area characterized by low-to-no snow, as well as seasonal precipitation and temperature, and annual hydrologic variables representing soil moisture, evapotranspiration (ET), the partitioning of incoming precipitation to evapotranspiration (ET/P), groundwater storage, and groundwater inflow to streams. MLR was done on all water years (P0: 1987-2024) and for each period as determined in the split analysis using pooled regression techniques (P1: 1987-2011 and P2: 2012-2024) to evaluate shifting predictor variable emphasis on streamflow generation. Results indicate that since 2012, peak SWE has lost statistical strength in its prediction of annual streamflow and runoff efficiency, and the indirect influence of spring temperature has emerged as critically important. Low-flow metrics remain largely influenced by soil moisture, vegetation water use and groundwater inflows with summer precipitation becoming a direct influence on minimum summer flow. Together, these data and Python-based analysis tools provide a framework for identifying the key watershed characteristics that control how streamflow responds to snow from year to year. The package also helps quantify uncertainty in statistical models and assess how snow–streamflow relationships vary across regions and over time. This dataset contains comma-separated values files (.csv), text files (.txt), python code files (.py), figure files (.png), and shapefiles (.cpg, .dbf, .prj, .sbn, .sbx, .shp, .xml). Further details on file contents and MLR execution can be found in the readme file and the FLMD files. Work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Effects of elastic softening and helium accumulation kinetics on surface morphological evolution of plasma-facing tungsten

Based on a continuous-domain model, capable of accessing the spatiotemporal scales relevant to fuzz formation on the surface of plasma-facing component (PFC) tungsten, we report self-consistent simulation results that elucidate the effects of elastic softening and helium (He) accumulation kinetics on the surface morphological response of PFC tungsten. The model accounts for the softening of the elastic moduli in the near-surface region of PFC tungsten, including both thermal softening at high temperature and softening due to He accumulation upon He implantation. The dependence of the elastic moduli on the He content follows an exponential scaling relation predicted by molecular-dynamics simulations, while the He content in the near-surface region of PFC tungsten evolves according to a first-order saturation kinetics, consistent with experimental and simulation results reported in the literature. We establish that this elastic softening accelerates both nanotendril growth on the PFC surface and the onset of fuzz formation. We also explore the role of the rate of He accumulation to a saturation level in the near-surface region of irradiated tungsten in the onset of fuzz formation. For PFC tungsten surfaces such as W(110) where, under typical irradiation conditions, the characteristic time scale for stress-driven surface diffusion is comparable to the characteristic time scale for He accumulation, we find that accelerating the rate of He accumulation accelerates the growth rate of nanotendrils emanating from the surface. Additionally, we present a systematic parametric study of the PFC surface morphological response to explore its dependence on the He accumulation kinetics that is controlled by the irradiation conditions for low-energy implantation. Finally, we introduce an incubation time for nanotendril growth on the PFC surface, a concept equivalent to that of incubation fluence discussed in the literature, to predict and explain the minimum exposure time required to observe fuzz formation on PFC tungsten surfaces.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Evaluation of a land-atmosphere coupling metric computed from a ground-based infrared interferometer

Land-atmosphere feedbacks are a critical component of the hydrologic cycle. Vertical profiles of boundary layer temperature and moisture, together with information about the land surface, are used to compute land-atmosphere coupling metrics. Ground based remote sensing platforms, such as the Atmospheric Emitted Radiance Interferometer (AERI), can provide high temporal resolution vertical profiles of temperature and moisture. When co-located with soil moisture, surface flux, and surface meteorological observations, such as at the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site, it is possible to observe both the terrestrial and atmospheric legs of land-atmosphere feedbacks. In this study, we compare a commonly used coupling metric computed from radiosonde-based data to that obtained from the AERI to characterize the accuracy and uncertainty in the metric derived from the two distinct platforms. This approach demonstrates the AERI’s utility where radiosonde observations are absent in time and/or space. Radiosonde and AERI based observations of the Convective Triggering Potential and Low-Level Humidity Index (CTP-HI low ) were computed during the 1200 UTC observation time and displayed good agreement during both 2017 and 2019 warm seasons. Furthermore, radiosonde and AERI derived metrics diagnosed the same atmospheric preconditioning based upon the CTP-HI low framework a majority of the time. When retrieval uncertainty was considered, even greater agreement was found between radiosonde and AERI derived classification. The AERI’s ability to represent this coupling metric well enabled novel exploration of temporal variability within the overnight period in CTP and HI low . Observations of CTP-HI low computed within a few hours of 1200 UTC were essentially equivalent, however with greater differences in time arose greater differences in CTP and HI low .

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Sea ice evolution along the Northern Sea Route and implications for trans-Arctic shipping from 2021 through 2060

Arctic surface temperatures warmed at twice the global average in the second half of the 20 th century due to Arctic amplification (AA), a phenomenon predominantly caused by regional polar changes, like the melting of perennial sea ice and reduced sea ice extent (leading to more solar radiation being absorbed by the ocean surface as opposed to being reflected back to space by the ice surface). AA is projected to reach a factor of three even if the climate is stabilized by the mid 21 st century by reduced greenhouse gas emissions. In all emissions scenarios, AA is projected to lead to temperature changes at least 2.4 times larger than global mean surface temperature changes occurring between 2070 and 2100. The ice-albedo feedback, which occurs when the polar-marine surface absorbs more radiation as highly reflective sea ice melts, is reversible such that the premise of a runaway process is no longer accepted as a realistic possibility. No matter what actions are taken to reduce CO 2 concentrations in the atmosphere from now on, two different methods of predicting an ice-free Arctic suggest that perennial sea ice will mostly disappear in September by the year 2050. If there is no reduction in anthropogenic CO 2 and methane emissions, that scenario could occur sooner than 2030. Defining Arctic navigability as safe and economic passage of Polar Class 7 cargo ships without need of an escorting icebreaker, no single trans Arctic ship route will be navigable year-round in the first half of the 21 st century, including in the strong emission scenarios. However, seasonal trans-Arctic navigability will increase this century. An estimate on the number of days per year that the Northern See Route (NSR) will be navigable in the future is beyond the scope of this report. Along Northern Sea Routes 5 and 6, which run close to the Russian coast and Yamal LNG plant, an ARC 7 ice class LNG tanker, the equivalent of a Polar Class 3 (PC3) vessel, will be at low risk in December through April at some point during the current decade. May will continue to entail some risk (more than April) along relatively short segments of these routes through the end of the next decade (2030-2039). Come June, snow rapidly melts away, and thereafter the underlying sea ice begins becomes thinner and less concentrated, greatly reducing risk. However, neither path is desirable for ARC 7 tankers due to shallow bathymetry in Sannikov Strait, and a more desirable path for these ships passes to the north of the New Siberian Islands (route 20, discussed further below). Conventional LNG tankers (i.e., non-ice-strengthened vessels according to the IMO classification), without icebreaker escorts, will continue to encounter dangerous or impassable conditions along many sections of the NSR through the end of this decade for many months of the year. Through 2039, these conventional tankers will be able to operate safely along NSR 6 from August through October. By 2040-2049, the span of safe operation increases to August through November, and by 2050-2059 it increases to July through November, assuming a northward deviation to avoid the Sannikov Strait.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Data-model files associated with the manuscript "The Effects of Spatial and Temporal Resolution of Gridded Meteorological Forcing on Watershed Hydrological Responses" (Shuai et al., 2022 HESS)

This data package contains the model inputs and outputs used in "The Effects of Spatial and Temporal Resolution of Gridded Meteorological Forcing on Watershed Hydrological Responses" (Shuai et al., 2022 HESS). The data.zip file contains the data used to drive the model simulations. The model.zip file contains the XML input file for ATS. The notebook.zip file contains the Jupyter notebooks for pre- and post- processing model results. The figures.zip file contains the raw figures associated with the manuscript.Meteorological forcing plays a critical role in accurately simulating the watershed hydrological cycle. With the advancement of high-performance computing and the development of integrated watershed models, simulating the watershed hydrological cycle at high temporal (hourly to daily) and spatial resolution (10s of meters) has become efficient and computationally affordable. These hyperresolution watershed models require high resolution of meteorological forcing as model input to ensure the fidelity and accuracy of simulated responses. In this study, we utilized the Advanced Terrestrial Simulator (ATS), an integrated watershed model, to simulate surface and subsurface flow and land surface processes using unstructured meshes at the Coal Creek Watershed near Crested Butte (Colorado). We compared simulated watershed hydrologic responses including streamflow, and distributed variables such as evapotranspiration, snow water equivalent (SWE), and groundwater table driven by three publicly available, gridded meteorological forcings (GMFs) -- Daily Surface Weather and Climatological Summaries (Daymet), Parameter-elevation Regressions on Independent Slopes Model (PRISM), and North American Land Data Assimilation System (NLDAS). By comparing various spatial resolutions (ranging from 400 m to 4 km) of PRISM, the simulated streamflow only becomes marginally worse when spatial resolution of meteorological forcing is coarsened to 4 km (or 30% of the watershed area). However, the 4 km resolution has much worse performance than finer resolution in spatially distributed variables such as SWE. Using temporally disaggregated PRISM, we compared models forced by different temporal resolutions (hourly to daily), sub-daily resolution preserves the dynamic watershed responses (e.g., diurnal fluctuation of streamflow) that are absent in results forced by daily resolution. Conversely, the simulated streamflow shows better performance using daily resolution compared to that using sub-daily resolution. Our findings suggest that the choice of GMF and its spatiotemporal resolution depends on the quantity of interest and its spatial and temporal scale, which may have important implications on model calibration and watershed management decisions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Effects of 9.5 years warming on SOC concentration and composition in bulk soil and density fractions

Original data of whole-soil warming experiment after 9.5 years at Blodgett Forest Research Station. The Blodgett Forest is a mixed coniferous temperate forest with Mediterranean climate. The annual air temperature is 12.5℃ and the annual precipitation is 1774 mm yr-1- The soil is mesic ultic Alfisol of granitic origin, equivalent to Dystric Cambisol according to The World Reference Base for Soil Resources (WRB) system. The soil is warmed down to 1 m at + 4℃ by vertically installed heating cables. At the time of soil sampling on 1 May 2023, the whole-soil warming experiment had been running for approximately 9.5 years, from January 2014 to May 2023. The dataset includes: - Bulk_EA: C, N content, δ13C, and CN ratio of bulk soil; - Density_fractionation: organic carbon concentration, δ13C, and C/N ratio of free light fraction (fLF), occluded ligh fraction (oLF), and heavy fraction (HF); - PCA_DRIFT_AUC: original data of area under the curve (AUC) values of eight carbon bond types integrated on diffuse reflectance infrared fourier transform spectroscopy for each soil sample and soil fraction, which are consequently used for principal component analysis (PCA); - DRIFTS_stability_index: the calculation of aliphatic C–H (3000–2800 cm-1) to aromatic C=C (1670–1600 cm-1) ratios for each bulk soil sample and soil fraction. All data are provided in CSV format and can be viewed using Microsoft Excel.

Climate change↗

DOE OSTI2021

Improving energy yield in photovoltaic modules with photonic structures

When installed outdoors, Si solar cells typically operate 20 – 30 K above ambient conditions. These elevated operating temperatures lead to diminished performance, reducing efficiency by ~0.4%/K for crystalline Si cells and decreasing module lifetime. Since much of the elevated temperature of the modules derives from parasitic absorption of sub-bandgap sunlight by the contacts, encapsulants, or other materials, reflectors that remove the sub-bandgap radiation from the module before it is absorbed promise to reduce the operating temperature. However, there are many competing strategies to reject this light: selective reflectors could be integrated into the outer module glass, at the textured Si-encapsulant interface, or on the rear. Furthermore, many of these approaches are complex and expensive to fabricate. The issue of heating can be even more critical in bifacial modules: high fractions of rear irradiance can lead to higher operating temperatures. Prior to this project, relatively little was known about the optical strategies for temperature reduction in bifacial modules. The goal of this project was to quantitatively assess optical strategies for improving the energy yield via a combination of improved anti-reflection and sub-bandgap light rejection. We studied two complementary concepts: the performance of ideal structures to determine the upper limits to performance, and low-complexity structures made from standard materials as cost-effective implementations. We also studied the performance of these reflectors in different types of modules, including Al BSF, PERC, and bifacial. Here we build from our prior accomplishments in developing mirror optimization techniques, and make use of simulation methods that accurately capture the optical properties of the modules throughout the visible and near-infrared spectrum. We connect these optical properties to energy yield under realistic, outdoor conditions. We also experimentally fabricated mirrors, incorporated them into mini-modules, and performed outdoor testing of performance at NREL. We found definitively that spectrally-selective mirrors on the outside of the module cover glass offer the best performance, both under ideal conditions and for realistic, low- complexity mirrors. Mirrors at the Si cell/encapsulant interface are limited by the multiple reflections needed to remove sub-bandgap light. Mirrors at the cell surface and mirrors on the rear also cannot remove all the sub-bandgap light, since parasitic absorption also occurs in the cover glass, encapsulant, and front contacts, depending on the module configuration. Also, cell texturing places higher demands on any cell-surface mirror because multiple interactions with the mirror are required to reject sub bandgap light which compounds losses. For bifacial modules, mirrors on the back do provide additional cooling and energy yield improvements, and the mirrors can have the same structure. In the monofacial case up to 3.3K of cooling is possible for mirrors on the front, compared to 2.2K at the cell surface and 1.2K at the rear. For bifacial modules, up to 2.4K of cooling is possible when mirrors are included at both the front and back. Our realistic designs utilize 4 – 6 layers of common photovoltaic materials such as SiO 2 , TiO 2 , and SiN x , making them more cost effective. Finally, we used our models to study thermal management in bifacial modules with different back lamination materials, finding that glass-glass laminated modules operated hotter than equivalent glass-polymer backsheet modules, due solely to the increased absorption of light incident from the rear of the module. Improving the energy yield by 2%, which is achievable with these realistic designs, would lead to LCOE reductions of $0.03/kWh by 2030. These benefits make photovoltaic panels more accessible to the public and may reduce associated maintenance costs.

14 SOLAR ENERGY↗

DOE OSTI2023

Throughfall-Reduction Drying Effects on Forest Floor Plant Tissue Biomass (Oi+Oe soil horizons) and Moisture Content in Four Lowland Panamanian Forests

Objectives: Climatic drying is predicted for many tropical forests, yet effects on soil properties across moisture and soil gradients within tropical forests remain poorly characterized, hampering predictions of forest-climate feedbacks. We hypothesized that drying would suppress soil CO2 fluxes (i.e., respiration) in already-drier tropical forests by further reductions in soil moisture, but increases CO2 fluxes in wetter tropical forests by alleviating anaerobiosis and soil saturation. We measured soil CO2 fluxes, soil moisture, soil temperature, and forest floor biomass during wet-dry cycles (2015 – 2022) in four Panamanian forests that vary in rainfall and soil fertility. We also surveyed all tree species and identified to species in 2018 and 2019.Results: We found that soil moisture peaked in the wet season and declined in the dry season. Measured soil CO2 fluxes declined in the dry season and peaked in the early wet season ahead of peak soil moisture, resulting in a lower soil moisture optimum for respiration than previously modeled. Forest floor biomass peaked in the dry season, in contrast. Chronic throughfall exclusion also suppressed soil moisture across the four forests to 20cm depths, and also initially suppressed soil CO2 fluxes across forests. There was sustained suppression of soil CO2 fluxes after four years in the wettest forest only (-28 ± 4% during the dry season), but elevated soil CO2 fluxes in a fertile forest after four years (+75 ± 28% during the late wet season). The unexpected negative drying effect in the wettest, most infertile forest could have resulted from reduced vertical flushing of nutrients into soils, as the drying effect increased with time. Including hydro-nutrient interactions in ecosystem models could improve predictions of tropical forest-climate feedbacks (results presented in Cusack et al. 2023). Datasets included: Datasets included here include .csv and .xls files for forest floor biomass (equivalent to the Oi+Oe soil horizons), and a subset include forest floor biomass moisture content (weight/weight) in the study plots. There is also a .kml file that includes coordinates for all 32 plots included in the study of four forests (n = 4 throughfall reduction and n = 4 control plots per site). No special software is needed to open these files.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Combined Mesonet and Tracker

Title: Combined Mesonet and Trackers (UNL Mobile Mesonets) Authors University of Nebraska PI: Adam Houston, UNL Professor (ahouston2@unl.edu) Mailing Address: 126 Bessey Hall P.O. Box 880340 Lincoln, NE 68588-0340 CoMeT Overview The University of Nebraska-Lincoln operates three Combined Mesonet and Tracker (CoMeTs). CoMeTs are Ford Explorers (model years 2015, 2017, and 2019) with forward-mounted suites of meteorological sensors and dual moonroofs, combining the capability of a mobile mesonet to collect near-surface observations with the capability of an unmanned aircraft systems (UAS) tracker vehicle, which enables an observer in the second row of seats to see the aircraft and maintain compliance with Federal Aviation Administration policies on UAS operation. The CoMeTs collect observations of slow temperature and humidity at ~2 m above ground level (AGL) using a Vaisala HMP155A, fast temperature at ~2 m AGL using a Campbell Scientific 109SS-L thermistor, pressure at ~2.5 m AGL using a Vaisala PTB210, wind speed and direction at ~3.25 m AGL using an R.M. Young 05103 propeller anemometer, and vehicle heading using a KVH Industries C-100 fluxgate compass (Barbieri et al. 2019). The HMP155A and 109SS-L thermistor are shielded and aspirated within a U-tube (Waugh and Frederickson 2010; Houston et al. 2016). This list of sensors is also included in the CoMeT data file metadata. Manufacturer specifications for these instruments are given in Table 1 of Hanft and Houston (2018). The reported measured quantities are summarized below.CoMeT-3 was funded through an equipment allocation included in the NSF TORUS award (AGS-1824649). Instrument Description The specific sensors included on each CoMeT are summarized in the table at the end of this section. In general each CoMeT collects observations of slow temperature and humidity at ~2 m above ground level (AGL) using a Vaisala HMP155, fast temperature at ~2 m AGL using a Campbell Scientific 109SS-L thermistor, pressure at ~2.5 m AGL using a Vaisala PTB210 barometer with a Gill pressure port, wind speed and direction at ~3.25 m AGL using an R.M. Young 05103 propeller anemometer, position using a Garmin 19x HVS receiver, and vehicle heading using a KVH Industries C-100 fluxgate compass. The HMP155 and 109SS are shielded and aspirated within a U-tube (Waugh and Frederickson 2010; Houston et al. 2016). Fast temperature and corrected RH measurements (using sensors housed within the U-tube) have a time constant of 10-12 s based on data collected across a temperature and RH shock during the CLOUD-MAP 2017 calibration/validation tests on June 26, 2017. Vehicle speed was < 10 kts for this test. CoMeT-1 CoMeT-2 CoMeT-3 Slow Temperature Slow RH Vaisala HMP155A-L20-PT Part #: 22280-7 Vaisala HMP155E Part #: E1AA11A0B1A1A0A Vaisala HMP155E Part #: E1AA11A0B1A1A0A Fast temperature Campbell Scientific 109SS-L20-PT Part #: 21448-3 Campbell Scientific 109SS-L12-PW Part #: 21448-109 Campbell Scientific 109SS-L12-PW Part #: 21448-150 Pressure Vaisala PTB-210 Part #: A1A1B Gill Pressure Port Part #: 61002 Vaisala PTB-210 Part #: A1A1B Gill Pressure Port Part #: 61002 Vaisala PTB-210 Part #: A1A1B Gill Pressure Port Part #: 61002 Wind RM Young 05103-L20-PT Part #: 18435-310 RM Young 05103-L20-PW Part #: 18435-244 RM Young 05103-L20-PW Part #: 18435-244 GPS Garmin GPS 19x HVS (NMEA 0183) Part #: 010-01010-00 Garmin GPS 19x HVS (NMEA 0183) Part #: 010-01010-00 Garmin GPS 19x HVS (NMEA 0183) Part #: 010-01010-00 Compass KVH C-100 Part #: 01-0177-15 KVH C-100 Part #: 01-0177-15 KVH C-100 Part #: 01-0177-15 Logger Campbell Scientific CR6-NA-XT-SW Part #: 28385-9 Campbell Scientific CR6-WIFI-XT-SW Part #: 28385-6 Campbell Scientific CR6-WIFI-XT-SW Part #: 28385-6 Data Collection and Real-Time Processing The reported measured quantities are summarized in the table below. Quantity Units Source Epoch time Seconds GPS Latitude and longitude Degrees GPS Altitude m GPS Pressure hPa PTB210 Temperature (fast) deg C 109SS-L Temperature (slow) deg C HMP155 RH (slow) % HMP155 Vehicle speed m/s GPS Vehicle heading deg C-100 and GPS In addition to the measured variables, several derived variables are calculated. Corrected/fast relative humidity (%) Relative humidity is adjusted to the fast temperature following Richardson et al. (1998) and Houston et al. (2016). Water vapor mixing ratio (g/kg) Dew point temperature (°C) Potential temperature (Kelvin) Virtual potential temperature (Kelvin) Equivalent potential temperature (Kelvin) Regular intercomparisons between all three CoMeTs were performed during TORUS 2019. Comparisons were also conducted between CoMeT-1 and CoMeT-2 during LAPSE-RATE (2018) on 14 July. In these intercomparisons, the vehicles were parked adjacent to each other aligned perpendicular to (and facing into) the wind. To minimize engine heating effects, intercomparisons were only conducted when the wind speed was >10 kts. Data Format Original data files for each deployment are saved as text files and then converted to NetCDF. NetCDF versions have units that are CF compliant and may not match the original units in the txt files. The naming convention for the NetCDF files is as follows: UNL.CoMeT3.{deployment date YYYYMMDD}.{start time of observation collection in UTC HHMM}.L2.{post-processing codes}.cdf example: UNL.CoMeT3.20190627.1931.L2.g1.f1.cdf Post-processing codes are included to track modifications to the raw data. These codes are closely connected to error flags associated with each record. Each letter corresponds to a particular instrument: g: GPS p: Barometer tf: Fast temperature ts: Slow temperature rh: Relative humidity f: Compass w: Wind monitor a: All instruments Each number corresponds to a particular post-processing action described more below. Measured and derived variables are included in the following table. Variable Heading Standard Name Units time Time seconds since 00:00:00, 01-01-1970 Alt Altitude meters lat Latitude degrees north lon Longitude degrees east fast_temp Air Temperature kelvin slow_temp Air Temperature kelvin pressure Air Pressure pascals logger_RH Relative Humidity percent calc_corr_RH Relative Humidity percent wind_speed Wind Speed meters per second wind_dir Wind From Direction degrees vehicle_dir Vehicle Direction degrees dewpoint Dew Point Temperature kelvin mixing_ratio Humidity Mixing Ratio g/g theta Air Potential Temperature kelvin theta_v Virtual Potential Temperature kelvin theta_e Equivalent Potential Temperature Kelvin error_flag The error_flag variable is a string that matches the post-processing codes listed above. All instruments will have an associated code, but will have a “0” if the datum is unchanged from the initial processed value. Error Codes The following table summarizes the error codes for data collected before 2020: Error Code Relevant CoMeT Description g1 1,2,3 Exact correction. GPS position and time reprocessed from raw data g2 1 As far as we can tell this is an exact correction to an error in the GPS time. During the correct time periods the time suddenly went backwards ~250s and stayed at this offset for 750s when it corrected itself. The offset was applied to the “time warp” period. p1 2 Approximate correction. Hole in the pressure tube connecting the pressure port to the barometer. Resulted in erroneously low air pressure measurements when the vehicle was in motion. Derived variables recalculated (dew point temperature [e depends on qv and p], water vapor mixing ratio, potential temperature, virtual potential temperature, equivalent potential temperature) a1 3 Exact correction. Missing data reprocessed from raw data a2 1 Bug fix to bias correction for ts1, ts2, and rh1: water vapor mixing ratio was off by a factor of 10 and virtual potential temperature was wrong because of this. f1 3 No correction, missing data. Fluxgate compass inoperable. Wind speed and direction calculated using GPS-derived vehicle heading instead. rh1 1 Approximate correction. Constant bias of +1.7% removed from relative humidity. Derived variables recalculated (corrected/fast relative humidity, dew point temperature, water vapor mixing ratio, virtual potential temperature, equivalent potential temperature) ts1 1 Approximate correction. Constant bias of +0.6 K removed from slow temperature. Derived variables recalculated (corrected/fast relative humidity, dew point temperature, water vapor mixing ratio, virtual potential temperature, equivalent potential temperature) ts2 1 Approximate correction. Constant bias of +1.0 K removed. Derived variables recalculated (corrected/fast relative humidity, dew point temperature, water vapor mixing ratio, virtual potential temperature, equivalent potential temperature) References Bolton, D., 1980: The Computation of Equivalent Potential Temperature. Mon. Wea. Rev., 108, 1046–1053, https://doi.org/10.1175/1520-0493(1980)108<1046:TCOEPT>2.0.CO;2. Hanft, W., and A. L. Houston, 2018: An Observational and Modeling Study of Mesoscale Air Masses with High Theta-E. Mon. Wea. Rev., 146, 2503–2524, https://doi.org/10.1175/MWR-D-17-0389.1.Wexler Houston, A. L., R. J. Laurence III, T. W. Nichols, S. Waugh, B. Argrow, and C. L. Ziegler, 2016: Intercomparison of unmanned aircraft-borne and mobile mesonet atmospheric sensors. Journal of Atmospheric and Oceanic Technology. 33, 1569-1582, doi: 10.1175/JTECH-D-15-0178.1. Lowe, P. R., 1977: An Approximating Polynomial for the Computation of Saturation Vapor Pressure. J. Applied Meteorology, 16, 100–103. Richardson, S. J., S. E. Frederickson, F. V. Brock, and J. A. Brotzge, 1998: Combination temperature and relative humidity probes: Avoiding large air temperature errors and associated relative humidity errors. Preprints, 10th Symp. On Meteorological Observations and Instrumentation, Phoenix, AZ, Amer. Meteor. Soc., 278–283. Waugh, S., and S. E. Frederickson, 2010: An improved aspirated temperature system for mobile meteorological observations, especially in severe weather. 25th Conf. on Severe Local Storms, Denver, CO, Amer. Meteor. Soc., P5.2. [Available online at https://ams.confex.com/ams/25SLS/techprogram/paper_176205.htm.]

54 ENVIRONMENTAL SCIENCES↗