RHOD Site - PNNL Surface Flux Station at 4m / Reviewed Data
This dataset contains preliminary processed turbulent fluxes and relevant statistics; 30-minute estimates
SEARCH · Engineering Papers
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.
This dataset contains preliminary processed turbulent fluxes and relevant statistics; 30-minute estimates
This dataset provides Level 1 (L1) orthorectified at-sensor radiance derived from measurements collected by the Imaging Spectrometer-1 (NIS-1) onboard the NEON (National Ecological Observatory Network) Airborne Observation Platform (AOP) for the 2025 Colorado Headwaters Ecological Spectroscopy Study (CHESS). NIS-1 captures light reflected from the Earth’s surface in 426 discrete wavelength bands as raw digital numbers (DNs; Level 0). These data are then calibrated to physical units (uW/cm²·sr·nm) following the processing steps described in the NEON Imaging Spectrometer Level 1B Calibrated Radiance Algorithm Theoretical Basis Document (ATBD; Gallery 2022). The data delivered here are the primary inputs for the surface reflectance product in “Custom surface reflectance, shade masks, and equivalent water thickness maps for the Colorado Headwaters Ecological Spectroscopy Study” (Carroll et al. 2026). For intertemporal comparison, the radiance data here are most directly relatable to the v2 radiance data in “NEON AOP Imaging Spectroscopy Survey of Upper East River Colorado Watersheds: Raw-Space Radiance and Observational Variable Dataset” (Goulden et al. 2018), to which the same processing methodology was applied. Together, the radiance and reflectance data enable users to exploit the unique reflection signatures of different surface objects for land cover classification, foliar trait mapping, plant vigor assessment, water content estimation, trace-element identification, and other scientific applications. The data were acquired over three study domains in the Upper Gunnison river basin: the upper East River watershed (CRBU); Almont Triangle and Taylor Canyon (ALMO); and Upper Taylor River watershed (UPTA) between 2025-06-13 and 2025-07-15. Within each domain, data are delivered by flightline as orthorectified and calibrated hyperspectral rasters in Hierarchical Data Format version 5 (HDF5) format, with radiance values provided in uW/cm²·sr·nm on a fixed, uniform Universal Transverse Mercator (UTM) grid at 1 meter spatial resolution. The radiance rasters include all 426 NIS-1 spectral bands, along with associated quality-assurance (QA) and diagnostic and ancillary layers needed for atmospheric correction workflows. Orthorectified radiance is produced from pushbroom spectrometer observations by applying NEON’s radiometric calibration (including bad pixel masking, dark subtract, dark pedestal shift correction, electronic panel ghost correction, grating ghost correction, deblur correction and flat-fielding) and spectral calibration (using spectral response function band centers and full-width at half-maximum intensity), followed by geolocation and regridding to the fixed grid. CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. Funding Acknowledgement: Field and remote-sensing data acquisition was performed under a grant from the National Aeronautics and Space Administration (80NSSC24K1005). This work was also 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.
(Fe 0.7 Co 0.3 ) 2 B are potential permanent magnets material due to its large saturation magnetization and high Curie temperature. However, it has moderate magnetocrystalline anisotropy (MCA) and low coercivity. One way to improve its coercivity is to combine the contributions from magnetocrystalline- and magnetic-shape anisotropy by preparing (Fe 0.7 Co 0.3 ) 2 B nanowires. We study the effects of size, morphology, and surface defects on the hard magnetic properties of nanowires using micromagnetic simulation. The hard magnetic properties of (Fe 0.7 Co 0.3 ) 2 B nanowire-bonded magnets are estimated, including the role of inter-wire magnetostatic interaction. By considering the existence of local reductions in MCA energy of up to 30% on the surface layer of nanowires, the anisotropic bonded magnet with a 65% vol. of (Fe 0.7 Co 0.3 ) 2 B nanowires would have typical remanence, B r = 7.6–8.4 kG, coercivity, H ci = 9.6–9.9 kOe, and maximum energy product, (BH) m = 14–17.8 MGOe. Developing effective technology for synthesizing nanowires and fabricating corresponding bonded magnets is promising for manufacturing practical magnets based on the magnetic phase with a relatively low or moderate MCA, such as (Fe 0.7 Co 0.3 ) 2 B.
Focusing on conditions of subsidence when low clouds are present, ground-based observations in both the North Atlantic and the Southern Ocean reveal strong relationships between cloud boundary (base and top heights) and different measures of lower tropospheric instability. The difference in potential temperature between the surface and 800 hPa (a metric called M) provides a stronger relationship than measures of inversion strength such as the lower tropospheric stability and estimated inversion strength. This is because (1) inversion strength itself does not correlate well with cloud boundaries, and (2) M contains information that appears important for cloud boundaries. These include the surface forcing through the use of sea surface rather than near-surface air temperature and an upper level close to the real cloud top. These results expand upon previous work on the importance of M as a predictor of cloud morphology. However, important differences are found in low-cloud conditions for the North Atlantic as compared to the Southern Ocean (for a given value of M): stronger inversions, deeper boundary layers, and much larger sea level pressures. Therefore, the relationship between cloud boundaries and M differs between the two regions. A general circulation model provides similar relationships as observed between M and both cloud top height and temperature but tends to place clouds higher and at colder temperatures than observed for a given M. This might cause issues with the representation of precipitation, cloud cover and radiation in the Southern Ocean.
Not Available
Not Available
Culverts allow roads to safely traverse small streams or drainage ditches, and their proper design is critical to ensure a safe and reliable transportation network. A correct modelization of these hydraulic structures becomes crucial in the assessment of flood footprints or discharge peak estimation in a risk evaluation plan. The question of how to include culverts comes up frequently when assembling a hydraulic model that requires the presence of as many singular elements as possible. In this work, three different culvert integrations with the surface domain are studied and compared in the context of a 2D shallow water (SW) model. All of them are based on the Federal Highway Administration (FHWA) formulation for the culvert discharge estimation but differ in complexity and in the interaction with the numerical model for surface flow, some of them as internal boundary conditions. Several steady and unsteady validation test cases are presented and the numerical results are compared with the predictions from HEC-RAS 1D and HY-8 software. Furthermore, the culvert area, shape and their sensitivity to the 2D computational mesh is also analyzed.
Information on menisci dynamics and equilibrium contact angles, θ e , is needed for modeling multiphase flow of fluids in geologic formations. The wettability of fracture faces is relevant to applications such as waterflooding in enhanced oil recovery and evaluation of caprock integrity for geologic carbon sequestration. We investigated the vertical capillary rise of air-water menisci on exposed fracture faces for a sandstone, a shale, and three granites. The experimental design approximated that employed in the Wilhelmy plate method. Replicate fracture faces were prepared by fracturing cylindrical rock cores using the Brazilian method and splitting the fractured cores apart along their fracture planes. Meniscus onset was imaged from a cross-sectional perspective using dynamic neutron radiography, and quantified with change point analysis. An existing model for meniscus onset fitted the experimental data very well (median R 2 = 0.96). Capillary rise was typified by t 2/3 scaling (where t is time), followed by a constant equilibrium height, z e . The average time taken for the menisci to reach z e was 0.39 s, with no differences between rock types. The menisci achieved a greater vertical extent on the fracture faces of the sedimentary rocks as compared to those of the granites. Apparent θ e values were calculated from the best fit estimates of z e using the Wilhelmy plate equation for a smooth and planar vertical surface. Despite the fractal roughness of some of the fracture faces, there was a significant linear relationship (R 2 = 0.80) between the apparent θ e values and intrinsic θ e values measured on polished surfaces of the same rocks using the sessile drop method. The regression equation indicated the apparent θ e values were consistently greater than the intrinsic values. This overestimation may be due to averaging of pixels at the air-water-solid interface, a required step in the image analysis. Further exploration of the range of applicability and measurement accuracy of this proposed new technique would be valuable.
In many thin film materials, nuanced interplays of interfacial energies control the surface morphology and rearrangement. This work evaluates polymer−solvent interactions and solvent-driven surface reconstructions via ex situ and in situ fluid cell Atomic Force Microscopy (fc-AFM) analysis of amphiphilic block copolymer (BCP) thin films upon exposure to deionized (DI) water. We examine the differences in surface morphology, whole-film swelling, and force response in thin films of polystyrene-block-poly(ethylene oxide) (PS-b-PEO) and polystyrene- block-poly[(allyl glycidyl ether)-co-(ethylene oxide)] (PS-b- P[AGE-co-EO]) processed into standing-up cylinder morphologies perpendicular to a silicon substrate (⊥C). Using Amplitude Modulation AFM (AM-AFM) and Amplitude-Phase Distance (APD) force spectroscopy, this work probes the mechanoresponsive nature of the dynamic surface layers of these films, unveiling surface layer stratification and surface chain rearrangement via minimal tip−sample stimulation. To help rationalize the observed reconfigurations, the energetic driving forces were estimated using the harmonic mean approximations of interfacial energies. Given the nonionizable nature of the minority P(AGE-co-EO) block and the energetic driving forces for chain mobility, this work shows how the elimination of unfavorable PS−water interfaces drives chain rearrangement and coverage of the PS surface by chains of the hydrophilic block. This work highlights considerations for increasing the heterogeneity and complexity of BCP thin films via random blocks and how those changes to local interfacial energies may drive larger scale film morphology reconstructions, with broader implications for tuning interface hydrophilicity.
Microsized pore parameters, such as pore size and distance between pores in a series of model EPDM rubbers, were determined in situ under the pressure of 500 psi using 129 Xe nuclear magnetic resonance (NMR) techniques: spin–lattice (T 1 ) and spin–spin (T 2 ) relaxation measurements, pulsed-field gradient (PFG) NMR, and two-dimensional exchange spectroscopy (2D EXSY). The T 1 /T 2 (>>1) ratio for the xenon confined in the pores is larger than that for nonconfined free xenon. This suggests that almost the entire pore surface interacts with xenon atoms like a closed pore. While these pores still connect each other through very narrow diffusion/exchange channels, it is possible to observe the echo decay in PFG-NMR and cross-peaks in 2D EXSY. The results show that both diffusion (D pore ≈ 2.1 × 10 –10 m 2 /s) and exchange (exchange rate, τ exch = a few tens of milliseconds) of xenon between a pore within the material and outer surface are prolonged. The exchange distances (l), which correspond to the xenon gas penetration depth, were estimated to be 70–100 μm based on the measured diffusion coefficients and exchange rate (1/τ exch ). NMR diffraction analysis reveals that pore size (a) and pore distance (b) are on the order of magnitude of micrometers and tens of micrometers, while the diffusion coefficients of xenon gas in the diffusion channels (D eff ) are about 10 –8 m 2 /s. Overall, this study suggests that the pores with a few micrometers connected through very narrow flowing channels with the length of several tens of micrometers are developed 70 to 100 μm below the rubber surface. Furthermore, the overall steady-state diffusion of xenon is slower, approximately 2 orders of magnitudes, than the diffusion in the channel between the pores. Finally, the pore and exchange distances correlated with the composition of rubbers showed that the properties of EPDM rubber as a high-pressure gas barrier could be improved by reducing the size of cracks and the depth of gas penetration by the addition of both carbon black and silica fillers.
Observations collected during cold-season precipitation periods at Utquagvik, Alaska and at the multidisciplinary drifting observatory for the study of Arctic climate (MOSAiC) are used to statistically analyze the relations among the atmospheric water cycle parameters including the columnar supercooled liquid and ice amounts (expressed as liquid-water and ice-water paths, i.e., LWP and IWP), the integrated water vapor (IWV) and the near-surface snowfall rate. Data come from radar and radiometer-based retrievals and from optical precipitation sensors. While the correlation between snowfall rate and LWP is rather weak, correlation coefficients between radar-derived snowfall rate and IWP are high (~0.8), which is explained, in part, by the generally low LWP/IWP ratios during significant precipitation. Correlation coefficients between snowfall rate and IWV are moderate (~0.45). Correlations are generally weaker if snowfall is estimated by optical sensors, which is, in part, due to blowing snow. Correlation coefficients between near-surface temperature and snowfall rates are low (r < 0.3). The results from the Alaska and MOSAiC sites are generally similar. These results are not very sensitive to the amount of time averaging (e.g., 15 min averaging versus daily averages). Observationally based relations among the water cycle parameters are informative about atmospheric moisture conversion processes and can be used for model evaluations.
Atmospheric organic compounds may deposit to Earth's surfaces via dry deposition, driven by concentration gradients, and wet deposition, driven by the scavenging of compounds by precipitation. Their removal by deposition has downstream impacts on concentrations of secondary organic aerosol by removing potential aerosol precursors that would otherwise oxidize to form particulate matter. The impact of deposition processes can consequently be considered as competition between rates of oxidation and deposition, but timescales for deposition are not well constrained. While timescales for dry deposition have been estimated and experimentally validated in the past, understanding of wet deposition of organics is still very limited. In this work, we estimate the wet deposition timescale for gas-phase organic compounds in the atmosphere as a function of Henry's law constants, H, using real-world precipitation frequency and size distributions at five globally-distributed sites. The wet deposition timescale decreases significantly with the increase of H until reaching a stable minimum for compounds with H > 10 5 M per atm. We estimate that the median wet deposition timescale for highly soluble gases is approximately 5 hours during a continuous rain event for all sites. However, median estimated timescales ranged from 80 to 200 hours, depending on location. Timescales are found to depend primarily on the frequency and duration of precipitation events rather than their intensity or size characteristics. Based on these data, we demonstrate that timescales for wet deposition of gases can be estimated at any given location using only basic precipitation information, without detailed or high-precision measurements.
Abstract According to recent field studies, almost half of the New Particle Formation (NPF) events occur aloft, in a residual layer, near the top of the boundary layer. Therefore, measurements of the meteorological parameters, precursor gas concentrations, and aerosol loadings conducted at the ground level are often not representative of the conditions where the NPFs take place. This paper presents new measurements obtained during the Turbulent Flux Measurements of the Residual Layer Nucleation Particles, conducted at the Southern Great Plains research site. Vertical turbulent fluxes of 3–10 nm‐sized particles were measured using a sonic anemometer and two condensation particle counters with nominal cutoff diameters of 3 nm and 10 nm mounted at the top of the 10‐m telescoping tower. Aerosol number size distribution (5–300 nm) was determined through the ground‐based Scanning Mobility Particle Sizers. The size selected (15–50 nm) particle hygroscopicity was derived with the Humidified Tandem Differential Mobility Analyzer. The ground level observations were supplemented by vertically‐resolved measurements of horizontal and vertical wind speeds and aerosol backscatter. The data analysis suggests that (a) turbulent flux measurements of 3–10 nm particles can distinguish between near‐surface and residual‐layer small particle events; (b) sub‐50 nm particles had a hygroscopicity value of 0.2, suggesting that organic compounds dominate atmospheric nanoparticle chemical composition at the site; and (c) current methodologies are inadequate for estimating the dry deposition velocity of sub‐10 nm particles because it is not feasible to measure particle concentration very near the surface, in the diffusion sublayer.
The effect of granulated blast furnace slag (GBFS) addition on the threshold chloride concentration (TCC) for rebar corrosion was investigated. A modified diffusion model, coupled with a time-dependent effective diffusion coefficient and surface chloride concentration, was proposed to predict the chloride profile. The corrosion initiation time was estimated based on the model predictions and the measured TCC. The results indicate that adding GBFS decreases the TCC by lowering the pH value of the pore solution. The evolution of corrosion potential and current density is found to obey a 3-parameter Weibull distribution. MIP tests show that adding GBFS contributes to refinement of pore structure by decreasing the fraction of large capillary pores. The time-dependent model exhibits good predictive strength and helps understand how GBFS addition delays the corrosion initiation by retarding the chloride diffusion, though a lower TCC is obtained.
The success of polyacrylic acid (PAA) to suppress Li dendrite growth suggests that the mechanical properties of polymer-based coatings, including the modulus, toughness, and interfacial adhesion are important design criteria. However, the measurement of the adhesion of thin PAA, as well as other polymer coatings to the reactive Li-metal anode surface is limited experimentally and challenging computationally. In this paper, a strategy was proposed to estimate the adhesion and delamination of the PAA(polymer)/Li interface, based on the bonding nature at the simpler PAA (oligomer)/Li interfaces using density functional theory calculations. It has been shown that the carboxylic acid groups in PAA reacted strongly with metallic Li, which significantly enhances the interfacial adhesion through the Li–O bonds formation, Li ionization and its incorporation into PAA, and –H or –OH termination of Li after decomposition of the COOH functional group. During delamination, it was found that the most likely PAA delamination route involved breaking partial Li–O bonds and lifting some ionized Li atoms from the Li-metal, especially for the Li atoms that showed a charge closer to +1 or are bonded with two O atoms from PAA. Based on the average bonding energies from PAA(oligomer)/Li interface delamination calculations, the work of separation, W sep , of the PAA(polymer)/Li interface was estimated to be ~1.0 (J/m 2 ). The high W sep of PAA (polymer)/Li was comparable with the Li 2 O/Li interface and higher than Li 2 CO 3 /Li and LiF/Li interfaces. This order correlated well with the areal density of Li–O bonds, which can serve as a descriptor for the interfacial adhesion. Furthermore, this computational approach can be applied to other interfaces with polymer-based coatings.
Abstract The companion part A paper (Stangeby et al 2022) reports a number of independent estimates indicating that high-duty-cycle DT tokamaks starting with pilot plants will likely experience rates of net erosion and deposition of solid PFC, plasma facing component, material in the range of 10 3 to 10 4 kg yr −1 , regardless of the material used. The subsequent redeposition of such large quantities of material has the potential for major interference with tokamak operation. Similar levels and issues will be involved if ∼continuous low-Z powder dropping is used for surface conditioning of DT tokamaks, independent of the material used for the PFC armor. In Stangeby et al (2022) (part A) it is proposed that for high-duty-cycle DT tokamaks, non-metallic low-Z refractory materials such as ceramics (graphite, SiC, etc) used as in situ replenishable, relatively thin—of order mm—claddings on a substrate which is resistant to neutron damage could provide a potential solution for protecting the main walls, while reducing the risk of degrading the confined plasma. Assessment of whether such an approach is viable will require information, much of which is not available today. Section 6 of part A identifies a partial list of major physics questions that will need to be answered in order to make an informed assessment. This part B report describes R&D needed to be done in present tokamaks in order to answer many of these questions. Most of the required R&D is to establish better understanding of low-Z slag generation and to identify means to safely manage it. Powder droppers provide a unique opportunity to carry out controlled studies on the management of low-Z slag in current tokamaks, independent of whether their protection tiles use low-Z or high-Z material.
Ozone (O 3 ) is a secondary pollutant in the atmosphere formed by photochemical reactions that endangers human health and ecosystems. O 3 has aggravated in Asia in recent decades and will vary in the future. In this study, to quantify the impacts of future climate change on O 3 pollution, near-surface O 3 concentrations over Asia in 2020–2100 are projected using a machine learning (ML) method along with multi-source data. The ML model is trained with combined O 3 data from a global atmospheric chemical transport model and real-time observations. The ML model is then used to estimate future O 3 with meteorological fields from multi-model simulations under various climate scenarios. The near-surface O 3 concentrations are projected to increase by 5 %–20 % over South China, Southeast Asia, and South India and less than 10 % over North China and the Gangetic Plains under the high-forcing scenarios in the last decade of 21st century, compared to the first decade of 2020–2100. The O 3 increases are primarily owing to the favorable meteorological conditions for O 3 photochemical formation in most Asian regions. We also find that the summertime O 3 pollution over eastern China will expand from North China to South China and extend into the cold season in a warmer future. Our results demonstrate the important role of a climate change penalty on Asian O 3 in the future, which provides implications for environmental and climate strategies of adaptation and mitigation.
The original data were collected during the field campaign of “Turbulent layers promoting New Particle Formation” experiment (NPFTURBULENCE; https://www.arm.gov/research/campaigns/aaf2024npfturbulence) over the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) Atmospheric Observatory (https://www.arm.gov/capabilities/observatories/sgp ) in north-central Oklahoma. The ARM Aerial Facility ArcticShark uncrewed aerial system (UAS) was based at Blackwell–Tonkawa Municipal Airport (IATA: BWL, ICAO: KBKN, FAA LID: BKN, 36.74475° N, 97.34918° W, 313.9 m MSL), for the field campaign from May 5 through May 29, 2024. The ArcticShark UAS performed 11 flights, including 10 research flights over the Central Facility of the ARM SGP to measure atmospheric state, turbulence, surface IT temperature and imagery, aerosol number concentration and size distribution. The current data set presents Best Estimate Aerosol Size Distribution: a merged aerosol size distribution composed of the data from 2 sensors: miniaturized Scanning Electrical Mobility Sizer (mSEMS) and Portable Optical Particle Spectrometer (POPS). The mSEMS data were interpolated to 1 second from “native” time resolution of about 15 second to match the other probe. The POPS data were converted from equivalent optical size into geometric size using value of aerosol refractive index of 1.477 from the HISCALE field campaign (same geographical area, altitudes, and time of year; http://www.arm.gov/campaigns/aaf2016hiscale ).