Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Surface estimation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Improved Data Reduction Algorithm for the Needle Probe Method Applied to In-Situ Thermal Conductivity Measurements of Lunar and Planetary Regoliths

The needle probe method (also known as the' hot wire' or 'line heat source' method) is widely used for in-situ thermal conductivity measurements on soils and marine sediments on the earth. Variants of this method have also been used (or planned) for measuring regolith on the surfaces of extra-terrestrial bodies (e.g., the Moon, Mars, and comets). In the near-vacuum condition on the lunar and planetary surfaces, the measurement method used on the earth cannot be simply duplicated, because thermal conductivity of the regolith can be approximately 2 orders of magnitude lower. In addition, the planetary probes have much greater diameters, due to engineering requirements associated with the robotic deployment on extra-terrestrial bodies. All of these factors contribute to the planetary probes requiring much longer time of measurement, several tens of (if not over a hundred) hours, while a conventional terrestrial needle probe needs only 1 to 2 minutes. The long measurement time complicates the surface operation logistics of the lander. It also negatively affects accuracy of the thermal conductivity measurement, because the cumulative heat loss along the probe is no longer negligible. The present study improves the data reduction algorithm of the needle probe method by shortening the measurement time on planetary surfaces by an order of magnitude. The main difference between the new scheme and the conventional one is that the former uses the exact mathematical solution to the thermal model on which the needle probe measurement theory is based, while the latter uses an approximate solution that is valid only for large times. The present study demonstrates the benefit of the new data reduction technique by applying it to data from a series of needle probe experiments carried out in a vacuum chamber on JSC-1A lunar regolith stimulant. The use of the exact solution has some disadvantage, however, in requiring three additional parameters, but two of them (the diameter and the volumetric heat capacity of the probe) can be measured and the other (the volumetric heat capacity of the regolith/stimulant) may be estimated from the surface geologic observation and temperature measurements. Therefore, overall, the new data reduction scheme would make in-situ thermal conductivity measurement more practical on planetary missions.

JSC-1A↗
NASA NTRS1998

Seasonal Cycles of Meridional Overturning and Heat Transport of the Indian Ocean

A general circulation model of the Indian Ocean is fitted to monthly averaged climatological temperatures, salinities, and surface fluxes using the adjoint method. Interannual variability is minimized by penalizing the temporal drift from one seasonal cycle to another during a two-year integration. The resultant meridional overturning and heat transport display large seasonal variations, with maximum amplitudes of 18 and 22 (x 10(exp 6) cubic m/s) for the overturning and 1.8 and 1.4 (x 10(exp 15) W) for heat transport near 10 S and 10 N, respectively. A dynamical decomposition of the overturning and heat transport shows that the time-varying Ekman How plus its barotropic compensation can explain a large part of the seasonal variations in overturning and heat transport. The maximum variations at 10 deg N and 10 deg S are associated with monsoon reversal over the northern Indian Ocean and changes of the easterlies over the southern Indian Ocean. An external mode with variable topography has a moderate contribution where the Somali Current and the corresponding gyre reverse direction seasonally. Contribution front vertical shear (thermal wind and ageostrophic shear) is dominant near the southern boundary and large near the Somali Current latitudes. The dominant balance in the zonally integrated heat budget is between heat storage change and heat transport convergence except south of 15 S. Optimization with seasonal forcings improves estimates of sea surface temperatures, but the annual average overturning and heat transport are very similar to previous results with annual mean forcings. The annual average heat transport consists of roughly equal contributions from time-mean and time-varying fields of meridional velocities and temperatures in the northern Indian Ocean. indicating a significant rectification to the heat transport due to the time-varying fields. The time-mean and time-varying contributions are primarily due to the overturning and horizontal gyre, respectively. Inclusion of TOPEX data enhances the seasonal cycles of the estimated overturning and heat transport in the central Indian Ocean significantly and improves the estimated equatorial zonal flows but leads 10 unrealistic estimates of the velocity structure near the Indonesian Throughflow region, most likely owing to the deficiencies in the lateral boundary conditions.

Lee, Tong↗
NASA NTRS1978

The effect of surface anisotropy on the accuracy of total ozone estimates from satellite observations

The total amount of ozone in a vertical column of the earth's atmosphere is being derived from satellite measurements of the intensity of ultraviolet sunlight scattered by the earth-atmosphere system. The algorithm for deriving the ozone amount utilizes the assumption that the earth's surface reflects the incident light isotropically according to Lambert's law. Natural surface reflection deviates more or less from this law. Two extreme examples of anisotropic reflection from dark ocean and from bright snow are analyzed by means of models for their effects on the derived values of ozone.

Fraser, R. S.↗
NASA NTRS2005

Asteroid Shape Reconstruction From Radar Observations

I estimate near-Earth asteroid 1992 SK's physical properties from radar delay-Doppler images, Doppler-only echo spectra and optical lightcurves. The images are not very strong, but place up to 20 (40 m by 160 m) pixels on the asteroid. The radar tracks are confined to subradar latitudes between 20 and 40 degrees but have complete rotational phase coverage. The echo spectra and optical lightcurves span approx.80 degrees of sky motion, providing geometric leverage to constrain the pole direction. The optical lightcurves are essential to accurate determination of the asteroid's shape and spin state. The asteroid is approx.1.4 km in maximum extent and mildly asymmetric, with an elongation of approx.1.5 and relatively subdued topography. The radar albedo is about 0.13 and the optical albedo about 0.3. The circular polarization ratio for the object is about 0.34, implying typical cm-scale surface roughness. I estimate the asteroid's period to be 7.3182+/-0.0003 hours and its pole direction as (99deg+/-5deg,-3deg+/-5deg) in ecliptic coordinates. The radar-refined orbital solution accurately predicts planetary close approaches between the years 826 and 2690. I have used my model to predict salient characteristics of radar images and optical lightcurves obtainable during the asteroid's March 2006 approach.

Busch, Michael J.↗
DOE OSTI2023

Interannual variability of spring and summer monsoon growing season carbon exchange at a semiarid savanna over nearly two decades

Eddy covariance measurements of land-atmosphere energy, carbon, and water exchange now span multiple decades at some sites, supporting an improved understanding of flux interannual variability (IAV) and its ecophysiological and physical controls. Most eddy covariance IAV studies have focused on temperate forest ecosystems, where carbon fluxes are large and flux records are longest – but also where IAV is much lower than in dryland regions, which have been identified as an essential driver of the trend and variability in the global terrestrial carbon sink. In this study, we leveraged 19 years of continuous micrometeorological measurements at the AmeriFlux US-SRM mesquite savanna site in southern Arizona, USA to quantify the IAV, trends, and drivers of carbon fluxes during the distinct spring and summer growing seasons. We also assessed the ability of modern satellite and land surface models to capture the IAV of seasonal water and carbon fluxes. Annual net ecosystem production (NEP) was small and highly variable (23 +/- 64 gC m –2 yr –1 ). Precipitation and associated measures of water availability determined most of the variability in NEP, largely through their influence on annual and seasonal gross ecosystem productivity (GEP) as opposed to ecosystem respiration (ER). Root-zone soil moisture captured between 73% (spring) and 85% (summer) of GEP variability and between 73% (spring) and 58% (summer) of ER variability. Throughout the study period, soil moisture and greenness increased with associated increases in GEP, ER and NEP. These trends were strongly influenced by very productive and wet summer growing seasons during the last two years, which were characterized by abundant understory grass cover. Typically, less than half of the variability in growing season GEP and evapotranspiration was captured by satellite-based estimates and land surface model simulations with local site forcing and calibration, highlighting the ongoing utility of long-term datasets to support careful model testing and improvement.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2011

Compact Polarimetry in a Low Frequency Spaceborne Context

Compact polarimetry has been shown to be an interesting alternative mode to full polarimetry when global coverage and revisit time are key issues. It consists on transmitting a single polarization, while receiving on two. Several critical points have been identified, one being the Faraday rotation (FR) correction and the other the calibration. When a low frequency electromagnetic wave travels through the ionosphere, it undergoes a rotation of the polarization plane about the radar line of sight for a linearly polarized wave, and a simple phase shift for a circularly polarized wave. In a low frequency radar, the only possible choice of the transmit polarization is the circular one, in order to guaranty that the scattering element on the ground is illuminated with a constant polarization independently of the ionosphere state. This will allow meaningful time series analysis, interferometry as long as the Faraday rotation effect is corrected for the return path. In full-polarimetric (FP) mode, two techniques allow to estimate the FR: Freeman method using linearly polarized data, and Bickel and Bates theory based on the transformation of the measured scattering matrix to a circular basis. In CP mode, an alternate procedure is presented which relies on the bare surface scattering properties. These bare surfaces are selected by the conformity coefficient, invariant with FR. This coefficient is compared to other published classifications to show its potential in distinguishing three different scattering types: surface, doublebounce and volume. The performances of the bare surfaces selection and FR estimation are evaluated on PALSAR and airborne data. Once the bare surfaces are selected and Faraday angle estimated over them, the correction can be applied over the whole scene. The algorithm is compared with both FP techniques. In the last part of the paper, the calibration of a CP system from the point of view of classical matrix transformation methods in polarimetry is proposed.

radar polarization diversity↗
NASA NTRS1976

A shock surface geometry - The February 15-16, 1967, event

An estimated shape is presented for the surface of the flare-associated interplanetary shock of February 15-16, 1967, as seen in the ecliptic-plane cross section. The estimate is based on observations by Explorer 33 and Pioneers 6 and 7. The estimated shock normal at the Explorer 33 position is obtained by a least-squares shock parameter-fitting procedure for that satellite's data; the shock normal at the Pioneer 7 position is found by using the magnetic coplanarity theorem and magnetic-field data. The average shock speed from the sun to each spacecraft is determined along with the local speed at Explorer 33 and the relations between these speeds and the position of the initiating solar flare. The Explorer 33 shock normal is found to be severely inclined and not typical of interplanetary shocks. It is shown that the curvature of the shock surface in the ecliptic plane near the earth-Pioneer 7 region is consistent with a radius of not more than 0.4 AU.

Lepping, R. P.↗
NASA NTRS2009

Mars Entry Atmospheric Data System Modelling and Algorithm Development

The Mars Entry Atmospheric Data System (MEADS) is being developed as part of the Mars Science Laboratory (MSL), Entry, Descent, and Landing Instrumentation (MEDLI) project. The MEADS project involves installing an array of seven pressure transducers linked to ports on the MSL forebody to record the surface pressure distribution during atmospheric entry. These measured surface pressures are used to generate estimates of atmospheric quantities based on modeled surface pressure distributions. In particular, the quantities to be estimated from the MEADS pressure measurements include the total pressure, dynamic pressure, Mach number, angle of attack, and angle of sideslip. Secondary objectives are to estimate atmospheric winds by coupling the pressure measurements with the on-board Inertial Measurement Unit (IMU) data. This paper provides details of the algorithm development, MEADS system performance based on calibration, and uncertainty analysis for the aerodynamic and atmospheric quantities of interest. The work presented here is part of the MEDLI performance pre-flight validation and will culminate with processing flight data after Mars entry in 2012.

Karlgaard, Christopher D.↗
NASA NTRS2014

Mars Entry Atmospheric Data System Modeling, Calibration, and Error Analysis

The Mars Science Laboratory (MSL) Entry, Descent, and Landing Instrumentation (MEDLI)/Mars Entry Atmospheric Data System (MEADS) project installed seven pressure ports through the MSL Phenolic Impregnated Carbon Ablator (PICA) heatshield to measure heatshield surface pressures during entry. These measured surface pressures are used to generate estimates of atmospheric quantities based on modeled surface pressure distributions. In particular, the quantities to be estimated from the MEADS pressure measurements include the dynamic pressure, angle of attack, and angle of sideslip. This report describes the calibration of the pressure transducers utilized to reconstruct the atmospheric data and associated uncertainty models, pressure modeling and uncertainty analysis, and system performance results. The results indicate that the MEADS pressure measurement system hardware meets the project requirements.

Karlgaard, Christopher D.↗
DOE OSTI2021

Global Temperature Responses to Large Tropical Volcanic Eruptions in Paleo Data Assimilation Products and Climate Model Simulations Over the Last Millennium

Abstract Large volcanic eruptions are one of the dominant perturbations to global and regional atmospheric temperatures on timescales of years to decades. Discrepancies remain, however, in the estimated magnitude and persistence of the surface temperature cooling caused by volcanic eruptions, as characterized by paleoclimatic proxies and climate models. We investigate these discrepancies in the context of large tropical eruptions over the Last Millennium using two state‐of‐the‐art data assimilation products, the Paleo Hydrodynamics Data Assimilation product (PHYDA) and the Last Millennium Reanalysis (LMR), and simulations from the National Center for Atmospheric Research Community Earth System Model‐Last Millennium Ensemble (NCAR CESM‐LME). We find that PHYDA and LMR estimate mean global and hemispheric cooling that is similar in magnitude and persistence once effects from eruptions occurring in short succession are removed. The estimates also compare well to Northern‐Hemisphere reconstructions based solely or partially on tree‐ring density, which have been proposed as the most accurate proxy estimates of surface cooling due to volcanism. All proxy‐based estimates also agree well with the magnitude of the mean cooling simulated by the CESM‐LME. Differences remain, however, in the spatial patterns of the temperature responses in the PHYDA, LMR, and the CESM‐LME. The duration of cooling anomalies also persists for several years longer in the PHYDA and LMR relative to the CESM‐LME. Our results demonstrate progress in resolving discrepancies between proxy‐ and model‐based estimates of temperature responses to volcanism, but also indicate these estimates must be further reconciled to better characterize the risks of future volcanic eruptions.

Tejedor, E.↗
DOE OSTI2023

How does humidity data impact land surface modeling of hydrothermal regimes at a permafrost site in Utqia$\mathrm{\dot{g}}$vik, Alaska?

Humidity is a basic and crucial meteorological indicator commonly measured in several forms, including specific humidity, relative humidity, and absolute humidity. These different forms can be inter-derived based on the saturation vapor pressure (SVP). In past decades, dozens of formulae have been developed to calculate the SVP with respect to, and in equilibrium with, liquid water and solid ice surfaces, but many prior studies use a single function for all temperature ranges, without considering the distinction between over the liquid water and ice surfaces. These different approaches can result in humidity estimates that may impact our understanding of surface-subsurface thermal-hydrological dynamics in cold regions. In this study, we compared the relative humidity (RH) downloaded and calculated from four data sources in Alaska based on five commonly used SVP formulas. These RHs, along with other meteorological indicators, were then used to drive physics-rich land surface models at a permafrost-affected site. Here, we found that higher values of RH (up to 40%) were obtained if the SVP was calculated with the over-ice formulation when air temperatures were below freezing, which could lead to a 30% maximum difference in snow depths. The choice of whether to separately calculate the SVP over an ice surface in winter also produced a significant range (up to 0.2 m) in simulated annual maximum thaw depths. The sensitivity of seasonal thaw depth to the formulation of SVP increases with the rainfall rate and the height of above-ground ponded water, while it diminishes with warmer air temperatures. These results show that RH variations based on the calculation of SVP with or without over-ice calculation meaningfully impact physically-based predictions of snow depth, sublimation, soil temperature, and active layer thickness. Under particular conditions, when severe flooding (inundation) and cool air temperatures are present, care should be taken to evaluate how humidity data is estimated for land surface and earth system modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Field‐representative evaluation of PID‐polarization in TOPCon PV modules by accelerated stress testing

Abstract Potential‐induced degradation‐polarization (PID‐p) can reduce module power, but how to project the extent to which PID‐p may occur in field conditions considering the factors of system voltage, condensed moisture, temperature, and illumination has not been clarified. Using tunnel oxide passivated contact (TOPCon) modules, this work demonstrates a method to test full‐size crystalline silicon PV modules for PID‐p to provide field‐representative results. In initial screening tests with positive or negative 1000 V electrical bias applied at 60°C for 96 h using Al foil electrodes on the glass surfaces, the module type exhibited reversible PID‐p only on the front face when the cell circuit was in negative voltage potential. No PID was detected on the rear after testing in either polarity. We then evaluated the PID‐p sensitivity on the front side under different UV irradiances while maintaining the glass surface wet to estimate real‐world susceptibility to PID‐p. The magnitude of the observed behavior was fit using a previously developed charge transfer and depletion by light model. Whereas power loss with −1000 V applied to the cell circuit at 60°C for 96 h in the dark was about 30%, testing the module front under 0.051 W·m −2 nm −1 at 340 nm UVA irradiation using fluorescent tubes, the mean degradation was only 3%. When the modules were tested in the dark for PID‐p with in situ dark current–voltage (I‐V) characterization, the thermal activation energy for degradation was 0.71 eV; for recovery in the dark, it was 0.58 eV. Whereas recovery from the degraded state at 60°C in the dark without voltage bias was 5% absolute in 38 h, rapid recovery of about 5% absolute was observed with 1000 W·s/m 2 exposure at 25°C using a flash tester.

14 SOLAR ENERGY↗
DOE OSTI2023

Evolution of titanium particle combustion in potassium perchlorate and air

Understanding titanium particle combustion processes is critical not only for characterizing existing pyrotechnic systems but also for creating new igniter designs. In order to characterize titanium particle combustion processes, morphologies, and temperatures, simultaneous spatially-resolved electric field holography and imaging pyrometry techniques were used to capture post-ignition data at up to 7 kHz. Due to the phase and thermal distortions present in the combustion cloud, traditional digital in-line holography techniques fail to capture accurate data. Here in this work, electric field holography techniques are used in order to cancel distortions and capture the three-dimensional spatial locations and diameters of the particles. In order to estimate the projected surface temperatures of the titanium particles, an imaging pyrometry method that ratios emission at 750 and 850 nm is utilized. Using these diagnostics, joint statistics are collected for particle size, morphology, velocity, and temperature. Results show that, early in the combustion process, the titanium particles are primarily oxidized by potassium perchlorate inside the igniter cup, resulting in projected surface temperatures near 3000 K. Later in the process, the particles interact with ambient air, resulting in lower surface temperatures around 2400 K and the formation of flame zones. These results are consistent with adiabatic flame temperature predictions as well as particle morphology observations of a titanium core with a TiO 2 surface. Late stage particle expansion, star fragmentation, and molten droplet breakup events are also observed using the time-resolved morphology and temperature diagnostics. These results illustrate the different stages of titanium particle combustion in pyrotechnic environments, which can be used to inform improvements in next-generation igniters.

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

An arbitrarily high-order three-dimensional Cartesian-grid method for reconstructing interfaces from volume fraction fields

Here Tthis work describes a newly developed, arbitrarily high-order Cartesian-grid method for reconstructing material interfaces from a volume fraction field. The method begins by identifying all of the grid cells in the volume fraction field that are intersected by the interface and need to be approximated by the reconstruction scheme. Finite-differences are used to calculate the gradient of the volume fraction field and provide an estimate of the surface normal in all of the interfacial grid cells. Groups of connected grid cells are then identified which all have the same dominant component of the normal vector. This grouping by orientation determines the proper dependent variable to use in the surface reconstruction (e.g. for a 2D curve, this step determines if the surface will be approximated by a function of x or y). A cumulative integral over the surface is constructed and fit using b-splines for two-dimensional problems or tensor-product b-splines for three-dimensional problems. This construction allows for the interface to be recovered through application of the second fundamental theorem of calculus. Fitting the cumulative integral with $\mathscr{N}$ th-order b-splines (or tensor-product b-splines) yields an ($\mathscr{N}$-1) th-order convergence rate of the interface shape. Differentiation of the b-spline interface function(s) allows for the high-order approximation of the normal vector and curvature to be obtained directly anywhere along b-spline. Together, the proposed reconstruction technique can achieve arbitrarily high mesh convergence rates. Validation tests are presented with mesh convergence rates ranging from fourth- to tenth-order.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A critical look at the prediction of the temperature field around a laser-induced melt pool on metallic substrates

The study of microstructure evolution in additive manufacturing of metals would be aided by knowing the thermal history. Since temperature measurements beneath the surface are difficult, estimates are obtained from computational thermo-mechanical models calibrated against traces left in the sample revealed after etching, such as the trace of the melt pool boundary. Here we examine the question of how reliable thermal histories computed from a model that reproduces the melt pool trace are. To this end, we perform experiments in which one of two different laser beams moves with constant velocity and power over a substrate of 17-4PH SS or Ti-6Al-4V, with low enough power to avoid generating a keyhole. We find that thermal histories appear to be reliably computed provided that (a) the power density distribution of the laser beam over the substrate is well characterized, and (b) convective heat transport effects are accounted for. Poor control of the laser beam leads to potentially multiple three-dimensional melt pool shapes compatible with the melt pool trace, and therefore to multiple potential thermal histories. Ignoring convective effects leads to results that are inconsistent with experiments, even for the mild melt pools here.

42 ENGINEERING↗
DOE OSTI2020

H i imaging of dwarf star-forming galaxies: masses, morphologies, and gas deficiencies

ABSTRACT The Giant Meter-wave Radio Telescope observations of the H i 21 cm-line emission from 13 nearby dwarf star-forming galaxies are presented. These galaxies are selected from the catalogues of Wolf−Rayet galaxies having very young (≤10 Myr) star formation. The ranges of star formation rates and stellar masses of the sample galaxies are 0.03–1.7 M⊙ yr−1 and 0.04–22.3 × 108 M⊙, respectively. The H i line emission is detected from 12 galaxies with peak column density >1 × 1021 cm−2. The 3σ H i column density sensitivities per channel width of 7 km s−1 for low (60 arcsec × 60 arcsec) resolution images are in the range 0.8–1.9 × 1019 cm−2. The H i channel images, moment images, global profiles, and mass surface density profiles are presented here. The average value of the peak H i mass surface density is estimated to be ∼2.5 M⊙ pc−2, which is significantly less compared to that in massive spiral galaxies. The scaling relations of $(M_{stars} + M_{\rm H\, I} + M_{\rm He})$versus Mdyn, gas fraction versus MB, $M_{\rm H\, I}$versus Mstars, H i-to-stellar mass ratio versus Mstars, and $M_{\rm H\, I}$versus $D_{\rm H\, I}$for the sample galaxies are estimated. These scaling relations can be used to constraint the key parameters in the galaxy evolution models. These galaxies are residing in group environment with galaxy density up to eight galaxy Mpc−3. An H i mass deficiency (with DEFH i > 0.3) is noticed in majority of galaxies for their optical diameters as compared to galaxies in field environments. Clear signatures of tidal interactions in these galaxies could be inferred using the H i images. Isolated H i clouds without known optical counterparts are seen in the vicinity of several galaxies. H i emission envelope is found to be having an offset from the optical envelope in several galaxies. Consistent with the previous studies on galaxy evolution in group environments, tidal interactions seem to play an important role in triggering recent star formation.

Jaiswal, S.↗
DOE OSTI2020

Temperature Measurement of Laser-Irradiated Metals Using Hyperspectral Imaging

Accurate noncontact surface-temperature measurements during laser-based materials processing remain challenging due to the difficulty of establishing reliable emissivity values as a function of temperature and wavelength. Direct measurement of emissivity is difficult, as the emissivity may be changing constantly in the laser-material interaction region, where the temperature gradients are extreme and surface displacement can complicate the measurement. Here, we present a hyperspectral imaging method using a multiwavelength camera to capture the spectral radiance in eight different bands and fit the wavelength-dependent radiance to Planck’s law. Time-resolved temperature measurements during microsecond pulsed-laser irradiation of a metal plate made of the titanium alloy Ti -6 Al -4 V provide temperature information about the irradiated surface with an estimated accuracy of ± 10 % . The extracted wavelength-dependent emissivity slope compares well with reported results for macroscale titanium melts obtained using IR spectroscopy, which measure temperature without any built-in assumptions about the emissivity. The results are directly applicable to temperature monitoring during welding, additive manufacturing, and semiconductor growth.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Data From: "Warming and snow loss increase reliance on old groundwater in a Colorado River headwater"

This repository contains the data and code associated with the paper titled "Warming and snow loss increase reliance on old groundwater in a Colorado River headwater," published in Nature Geoscience, 2026. This study seeks to answer how various ages of groundwater interact with mountainous streamflow in mountainous headwaters such as the East River. It includes various model-data processing scripts, primarily for ParFlow-CLM analysis of simulated water years 2015-2021, and two numerical warming experiments (+2.5 and +4.0 degrees C), including run scripts, forcing scripts, and post-processing, as well as comparison to observation datasets, detailed below. This data requires the use of R (.r, .rmd), Python (.py), Jupyter Notebook or Jupyter Lab (.ipynb), ParFLOW-CLM, EcoSLIM. Further information on the use of all file formats mentioned below (e.g. .tff. .nc) are provided within the associated scripts and directory where the files are located. Contents & Usage ASO/: ​​Contains the bash and python scripts used to convert airborne snow observatory (ASO) data (ASO, 2023) in various data formats (georeferenced tiff file, NetCDF, UTM, and to latitude/longitude) then regrided to the ParFlow equivalent grid. Output data are in regrid_regll_data.zip and subsequently visualized and analyzed in plot_and_compare.py for Supplementary Figures A14 and A15. The wksht_ASO_comparison.xlsx spreadsheet is used to calculate the data for Supplementary Figure A16. EcoSLIM/: Contains the scripts and input files to run the EcoSLIM particle tracking simulations (/run_scripts) and the post-processing python script (/plot_scripts/eco_agedist_plots.ipynb). Jasechko et al./: Contains the jupyter notebook (Extract_Elevation.ipynb) to determine the outlet elevations of the 260 watersheds used in Jasechko et al. (2016), and the corresponding table, Table_S1_Watersheds_alt.csv. Used to create Supplementary Information Figure A2. PLM_Wells/: Contains the QA/QC-ed groundwater level time series of the PLM-1 and PLM-6 Monitoring Wells from Faybishenko et al. (2023), reformatted to water years used for Supplementary Figures A19 and and A20. ParFlow/: Contains the input files and run scripts to run ParFlow-CLM (/run_scripts), the python and tool command language (Tcl) scripts to create and distribute the ParFlow forcing simulation files (/forcing), and various scripts and intermediary files to analyze the model outputs (/post_process). SQUIRE/: Contains the processing scripts and intermediary files for the Surface QUantitatIve pRecipitation Estimation (SQUIRE) data (Grover, 2023) used to generate Supplementary Figure A18. USGS_Streamflow/: Contains the raw and gap-filled United States Geological Survey streamflow data (U.S. Geological Survey, 2026) used at the Almont station (site number 09112500). Gap-filling is performed in the R script with data from the Taylor station (site number 09110000). (/USGS_09112500_EAST_RIVER_AT_ALMONT_GAP_FILLED/code_almont_streamflow_gap_fill.Rmd). discharge/: Contains the gap-filled discharge data at the Watershed Function SFA East River pumphouse site (Newcomer et al., 2022) used to generate Supplementary Figure A13 and to compute hourly Nash-Sutcliffe model efficiency coefficients (NSE) in Table A4. snotel_and_flux_tower/: Contains the snow telemetry data (U.S. Department of Agriculture, 2024) from the Butte (site ID 380) and Schofield (site ID 737) stations, reformatted by water year, accessed with the snotelr R package. Used to create Supplementary Figure A17. Also contains the flux tower observational data (FluxTower_Pumphouse_ESS-DIVE.ET_only.h.txt) from Ryken et al. (2022) and sap flux transpiration data (MaxB_Transpiration_5Sites.daily_sums.h.txt) from Ryken (2021), used to create Supplementary Figures A22 and A23, respectively. Raw EcoSLIM model outputs are in excess of 24TB, and are stored on National Energy Research Scientific Computing Center (NERSC) and publicly available via the external link provided in the paper.

atmospheric warming↗