Search NASA⌕ Search

SEARCH · Search NASA

Results for “atmospheric tracers”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Simulating Wind Around Isolated Buildings with the System for Atmospheric Modeling

A method to represent the flow around obstacles such as buildings is implemented in the System for Atmospheric Modeling (SAM). The method, which we call a Quasi-Solid Box Method (QSBM), introduces a simple modification to the anelastic equations that forces the flow to stagnate within the obstacle's boundaries. The performance of the modified SAM is evaluated using CEDVAL (Compilation of Experimental Data for Validation of Microscale Dispersion Models) wind tunnel measurements of the wind and tracer dispersion around a single rectangular building. All major observed features of the flow are quite well reproduced. The dispersion of the tracer released at the building base is also simulated. The results of simulation of a flow around a cubic building rotated by 45 degrees relative to the flow, and around a building in the form of a cylinder of aspect ratio of one, are also presented.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Studying Aerosol, Clouds, and Air Quality in the Coastal Urban Environment of Southeastern Texas

A multi-agency succession of field campaigns was conducted in southeastern Texas during July 2021 through October 2022 to study the complex interactions of aerosols, clouds and air pollution in the coastal urban environment. As part of the Tracking Aerosol Convection interactions Experiment (TRACER), the TRACER- Air Quality (TAQ) campaign the Experiment of Sea Breeze Convection, Aerosols, Precipitation and Environment (ESCAPE) and the Convective Cloud Urban Boundary Layer Experiment (CUBE), a combination of ground-based supersites and mobile laboratories, shipborne measurements and aircraft-based instrumentation were deployed. These diverse platforms collected high-resolution data to characterize the aerosol microphysics and chemistry, cloud and precipitation micro- and macro-physical properties, environmental thermodynamics and air quality-relevant constituents that are being used in follow-on analysis and modeling activities. We present the overall deployment setups, a summary of the campaign conditions and a sampling of early research results related to: (a) aerosol precursors in the urban environment, (b) influences of local meteorology on air pollution, (c) detailed observations of the sea breeze circulation, (d) retrieved supersaturation in convective updrafts, (e) characterizing the convective updraft lifecycle, (f) variability in lightning characteristics of convective storms and (g) urban influences on surface energy fluxes. The work concludes with discussion of future research activities highlighted by the TRACER model-intercomparison project to explore the representation of aerosol-convective interactions in high-resolution simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The Amazon and La Plata River Basins as Moisture Sources of South America: Climatology and Intraseasonal Variability

Abstract Land‐atmosphere interactions are critical for precipitation (PPT) over South America where terrestrial evapotranspiration (ET) constitutes a significant fraction of moisture for rainfall over the ecologically and socio‐economically vital Amazon (AMZ) and La Plata (LPB) river basins. We quantify the contribution of ET from AMZ and LPB to PPT over the continent with a focus on the intraseasonal time scale. Using numerical water tracers embedded in the Weather Research and Forecasting model we track the moisture originating from the two basins. Our findings indicate that approximately 40% of annual rainfall over the eastern foothills of the Andes originates as AMZ ET, and nearly 30% of rainfall over northern Argentina originates as LPB ET. Analysis of moisture transport during both phases of the dominant intraseasonal oscillation pattern over South America reveals an intraseasonal “sloshing” of LPB moisture between the South Atlantic convergence zone (SACZ) and southeastern South America (SESA) regions. AMZ and LPB each supply approximately 6% of moisture for SACZ PPT during periods of intraseasonal enhancement (positive anomalies), highlighting the importance of moisture from the Atlantic Ocean. For the SESA region, LPB supplies 26% of the moisture for PPT during periods of intraseasonal enhancement while AMZ supplies 5%.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Data and script associated with “Shifts in Rain-Snow Partitioning Drive Faster Water Transit Times in the US Pacific Northwest”

This data package contains the data and code to use and run the Water Tracer enabled version of the Weather Research and Forecasting Hydrologic model (WT-WRF-Hydro) with the Sequential Precipitation Input Tagging (SPIT) framework. It is associated with the publication “Shifts in Rain-Snow Partitioning Drive Faster Water Transit Times in the US Pacific Northwest” published in Scientific Reports (Butler et al., 2026; https://doi.org/10.1038/s41598-026-46539-1). We use the Continental U.S. (CONUSII; Rasmussen et al., 2021) dataset to force the model with an historical climate (2006–2013) and a future climate (2086–2093) with a representative carbon pathway (RCP) 8.5 scenario. We use the model to calculate water transit times in five headwater catchments within the U.S. Pacific Northwest. We also show key hydrologic and environmental variables that affect water transit times and changes in the future. Finally, we use observed data to validate the model such as stream water isotopes, snowpack characteristics, and stream discharge. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. In addition to a readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions. The data package consists of 11 folders: (1) "Figures" contains the exported figures used in the manuscript; (2) "Model_Isotope_Date" contains the WT-WRF-Hydro isotope date used in model validation; (3) “Model_Outputs_Future” contains the WT-WRF-Hydro future climate outputs; (4) “Model_Outputs_Historical” contains the WT-WRF-Hydro historical climate outputs; (5) “Model_Outputs_Weights_Areas” contains the WT-WRF-Hydro weights per catchment used to calculate water transit times and isotopes in stream water; (6) “MODIS_data_scripts” contains data used to validate snow conditions in the study area; (7) “Observed_Flow_Data” contains the observed streamflow data used in model validation; (8) “Observed_Isotope_Data” contains the observed stream water isotope data used in model validation; (9) “Scripts” contains the Python scripts used to general results and the figures; (10) “Statistic_Outputs” contains the water transit time statistical outputs reported in this manuscript; (11) “Validation_SNOTEL” contains the SNOTEL data used in model validation. The files in this data package have the following file extensions: .tif, .txt, .csv, .pdf, .py, .jpg, and .png.

American River↗
DOE OSTI2022

Coupled abiotic-biotic cycling of nitrous oxide in tropical peatlands

Atmospheric nitrous oxide (N 2 O) is a potent greenhouse gas thought to be mainly derived from microbial metabolism as part of the denitrification pathway. Here we report that in unexplored peat soils of Central and South America, N 2 O production can be driven by abiotic reactions (≤98%) highly competitive to their enzymatic counterparts. Extracted soil iron positively correlated with in situ abiotic N 2 O production determined by isotopic tracers. Moreover, we found that microbial N 2 O reduction accompanied abiotic production, essentially closing a coupled abiotic-biotic N 2 O cycle. Anaerobic N 2 O consumption occurred ubiquitously (pH 6.4–3.7), with proportions of diverse clade II N 2 O reducers increasing with consumption rates. Our findings show that denitrification in tropical peat soils is not a purely biological process but rather a ‘mosaic’ of abiotic and biotic reduction reactions. Here, we predict that hydrological and temperature fluctuations differentially affect abiotic and biotic drivers and further contribute to the high N 2 O flux variation in the region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

DTS data from TRACER TBS

Distributed temperature-sensing data collected using a Silixa-XT DTS and 50/125 MM fiber. Data are calibrated using two water baths at the surface and an iMet aloft at the end of the fiber and tether. The data are then assigned an altitude for the TBS flight.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Inorganic Hydrogeochemistry in the 21st Century

Abstract Chemical and isotopic processes occur in every segment of the hydrological cycle. Hydrogeochemistry—the subdiscipline that studies these processes—has seen a transformation from “witch's brew” to credible science since 2000. Going forward, hydrogeochemical research and applications are critical to meeting urgent societal needs of climate change mitigation and clean energy, such as (1) removing CO 2 from the atmosphere and storing gigatons of CO 2 in soils and aquifers to achieve net‐zero emissions, (2) securing critical minerals in support of the transition from fossil fuels to renewable energies, and (3) protecting water resources by adapting to a warming climate. In the last two decades, we have seen extensive activity and progress in four research areas of hydrogeochemistry related to water‐rock interactions: arsenic contamination of groundwater; the use of isotopic and chemical tracers to quantify groundwater recharge and submarine groundwater discharge; the kinetics of chemical reactions and the mineral‐water interface's control of contaminant fate and transport; and the transformation of geochemical modeling from an expert‐only exercise to a widely accessible tool. In the future, embracing technological advances in machine learning, cyberinfrastructure, and isotope analytical tools will allow breakthrough research and expand the role of hydrogeochemistry in meeting society's needs for climate change mitigation and the transition from fossil fuels to renewable energies.

Zhu, Chen (ORCID:0000000153746787)↗
DOE OSTI2025

Entrainment, Detrainment, and Dilution of Dry and Moist Atmospheric Thermals

Here this study examines the entrainment, detrainment, and dilution of dry and moist (cloud) atmospheric thermals in large-eddy simulations. In a neutrally stable environment (with respect to dry dynamics), moist thermals have an increase in radius R with thermal height z t (α ≡ dR/dz t ) about 4 times smaller compared to dry thermals when density stratification is considered and ~2.4 times smaller without density stratification (i.e., applying the Boussinesq approximation). An analytic expression relating α to several dimensionless parameters is derived from the thermal impulse–circulation relation to clarify the factors impacting α. This expression shows that the difference in buoyancy structure between moist and dry thermals, with buoyancy concentrated in the central cores of moist thermals owing to latent heating, explains their smaller spreading rates. Individual contributions of entrainment and detrainment are analyzed using a direct parcel-based approach in the simulations. Moist thermals have similar fractional detrainment but much smaller fractional entrainment rates compared to dry thermals, consistent with the differences in α. Despite having smaller α, moist thermals are similarly dilute (quantified by a passive tracer) as dry thermals because of their greater mixing efficiency with the environment. Thus, moist thermals are substantially dilute but expand much less in size/volume as they rise compared to dry thermals. The α values for moist thermals in (dry) neutral and statically stable environments are similar, but fractional entrainment and especially detrainment rates are greater in the stable environment. Large detrainment rates are associated with a breakdown of the broader thermal vortex ring structure, especially with low environmental relative humidity, attributed in part to evaporation and buoyancy reversal.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

TRACER-Tethersonde Ozone data

TRACER-Tethersonde, a sub-campaign of TRACER, took place during the TRACER June-September 2022 intensive operating period (IOP) at the S3 Ancillary site (29.33°N, 95.74°W) near Guy, TX. The Tethered Balloon System (TBS) operated during the first two weeks of each month during the TRACER IOP, making multiple (~4) up and down vertical profiles each day that could reach as high as ~1 km. TRACER-Tethersonde included an En-Sci electrochemical cell (ECC) ozonesonde on the TBS. The ozonesonde was connected to an iMet-4RSB radiosonde, and the overall data collected included ozone, relative humidity, temperature, and altitude.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Optical Particle Measurements during EPCAPE Field Campaign Report

This campaign requested the deployment of the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) User Facility optical particle counter (OPC) at the first ARM Mobile Facility (AMF1) located at the Scripps Pier in La Jolla, California during the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE). The addition of the OPC was requested for two reasons. (1) Close the gap between the scanning mobility particle sizer (SMPS) and aerodynamic particle sizer (APS) size distribution from the Aerosol Observing System (AOS) measurements. (2) Principal investigator Petters has been working with Tracking Aerosol Convection Interaction Experiment (TRACER) data to compute particle fluxes from Doppler lidar (Petters et al. 2024). Briefly, backscatter flux is obtained using the eddy covariance technique using the Doppler vertical velocity and attenuated backscatter. Building upon prior studies, we were able to relate backscatter to particle number concentration by calibrating the lidar retrievals against optical particle counter-measured ground-based aerosol size distribution and radiosonde-interpolated relative humidity at lidar sample height. Performing similar analysis was of interest to EPCAPE to better understand the emissions and vertical transport of large particles into the overlying stratus clouds. However, as stated above, this analysis requires an optical size distribution that covers the 0.3-30-μm-diameter size range. The OPC was deployed between 2023-04-14 and 2024-02-14. The deployment, data quality analysis, and data archiving was handled by the DOE ARM instrument mentor team without additional involvement by the principal investigator. Data quality was marked as “routine” for the majority of the campaign.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

An Improved Convection Parameterization with Detailed Aerosol–Cloud Microphysics for a Global Model

Abstract A new microphysical treatment that includes aerosol–cloud interactions and secondary ice production (SIP) mechanisms is implemented in the convection scheme of the Community Atmosphere Model, version 6 (CAM6). The approach is to embed a 1D Lagrangian parcel model in the bulk convective plume of the existing deep convection parameterization. Aerosol activation, growth processes including collision/coalescence, and three processes of SIP mechanisms, two of which are normally overlooked in atmospheric models, are represented in this embedded parcel model. These microphysical processes are treated with a hybrid bin/bulk scheme and a high spatial and temporal resolution for the integration of the embedded parcel in 1D, allowing vertical velocity to determine the microphysical evolution following the in-cloud motion during ascent. Simulations of an observed case (Midlatitude Continental Convective Clouds Experiment) of a mesoscale convective system in Oklahoma, United States, with a single-column model (SCAM) version of CAM, are compared with aircraft in situ and ground-based observations of microphysical properties from the convection and precipitation. Results from the validation show the new microphysical scheme has a good representation of the ice initiation in the bulk convective plume, including the known and empirically quantified pathways of primary and secondary initiation, with benefits for the accuracy of properties of its supercooled cloud liquid. The sensitivity simulations and use of tagging tracers for the validated simulation confirm that the newly included SIP mechanisms are of paramount importance for convective microphysics and can be successfully treated in the global model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Understanding aerosol properties in convective outflows during TRACER

Convective outflow boundaries commonly form across the southeastern U.S., often originating from storms initiated along sea-breeze fronts. As outflows spread, merge, and collide, they frequently trigger secondary convection in air masses already modified by primary storm outflow. This resulting air mass carries a distinctive aerosol population, which can influence the development of secondary convective cells. During the U.S. DOE TRACER field campaign in summer 2022, we investigated how outflows affect aerosol size, composition, and hygroscopicity across Greater Houston, Texas. Using a novel detection algorithm, we identified and verified 72 outflow boundaries from surface meteorological observations. In the relatively simple continental environment northwest of Houston, a common relationship emerged between CCN-inferred hygroscopicity (κ) and size-resolved aerosol concentrations. The κ of aerosols typically decreased in size ranges where aerosol concentrations increased following outflow passage and increased where concentrations decreased. In contrast, at the Atmospheric Radiation Measurement’s (ARM) fixed site in La Porte, Texas, the aerosol response depended on outflow direction, shaped by concentrated industrial and shipping activity to the north and south/southeast along the Houston Ship Channel. Concurrent size distribution and composition changes suggested case-specific variability in κ, with either reinforcing or offsetting effects. Outflows passing over industrial/shipping sectors often promoted less hygroscopic aerosol populations while cleaner residential outflows often favored more hygroscopic aerosol populations. These results highlight key aerosol changes following outflow passage, providing a foundation for improved understanding of aerosol-cloud interactions related to secondary convection initiation along or in the wake of outflow boundaries.

Thompson, Seth A. [Texas A & M Univ., College Stat↗
DOE OSTI2023

Low-Temperature Gas-Phase Formation of Methanimine (CH 2 NH; X 1 A')–the Simplest Imine–under Single-Collision Conditions

The D1-methanimine molecule (CHDNH; X 1 A') – the simplest (deuterated) imine – has been prepared through the elementary reaction of the D1-methylidyne (CD; X 2 Π) with ammonia (NH 3 ; X 1 A 1 ) under single collision conditions. As a highly reactive species with a carbon-nitrogen double bond and a key building block of biomolecules like amino acids and nucleobases, methanimine is of particular significance in coupling the nitrogen and carbon chemistries in the interstellar medium and in hydrocarbon-rich atmospheres of planets and their moons. However, the underlying formation mechanisms of methanimine in these extreme environments are still elusive. Here, the directed, low-temperature gas-phase formation of D1-methanimine will deepen our fundamental understanding of low-temperature molecular growths processes via carbon-nitrogen bond coupling. Considering the recent detection of the interstellar D1-methylidyne radical, the investigation of the CD–NH 3 system also suggests a promising pathway for future astronomical observations of D1-methanimine as a molecular tracer of gas phase deuterium-enrichment in deep space.

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

Wind lidar operations and observations from Tundra Pigeon

Lawrence Livermore National Laboratory brought a ZX-300 profiling lidar to the Tundra Pigeon experiment for the purpose of collecting wind measurements in the lower atmospheric boundary layer. The ZX-300 lidar is a portable Doppler lidar which uses light to track naturally occurring aerosols across a measurement cone, thereby deriving wind speed (horizontal and vertical) and wind direction. The ZX-300 is programable between the heights of 10 m and 300 m above ground level and additionally has a 1 m onboard meteorological sensor for collecting measurements of air temperature, relative humidity, air pressure and wind at 1 m height. We programmed the ZX-300 to target altitudes of most interest to the tracer experiment. These heights were 10, 15, 20, 30, 38, 50, 75, 100, 125, 150 and 200 m. Note that the 38 m level is a fixed calibration height and cannot be changed. This measurement strategy prioritized winds close to the surface while ensuring that we’d also collect information about the winds aloft. The wind measurements are collected directly over the lidar instrument and represent a vertical profile of the winds over that location. However, given that there were minimal terrain and vegetation differences in the area, we’d expect these wind conditions to be representative of a larger volume area, discussed in more detail below.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Data and scripts associated with “Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics within Water-Tagging Enabled Hydrologic Models” (v3)

This data package is associated with the publication “Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics within Water-Tagging Enabled Hydrologic Models” submitted to Journal of Advances in Modeling Earth Systems (Butler et al. 2025). This study developed the Sequential Precipitation Input Tagging (SPIT) framework to tag input precipitation and estimate water transit times and hydrologic tracers. SPIT tags all precipitation events at regular intervals over an extended period (monthly tags over seven years) in a hydrologic model from 2016-2022. SPIT is applied at six National Ecological Observatory Network (NEON) sites across the continental United States to calculate transit time distributions (TTD) and derive from these mean transit times (MTT), fractions of young water (Fyw), and hydrologic tracer concentrations in stream water (δ18O) within a water-tagging enabled version of the Weather Research and Forecast (WT-WRF-Hydro) model with national water model (NWM) configurations. We go on to validate WT-WRF-Hydro estimates against Butler et al. (2023), who analyzed the same NEON sites using stable water isotope data to estimate water transit times. This new tracking method provides a detailed picture of water movement and helps improve predictions about water availability in the future. This data package was originally published in January 2025. It was updated May 2025 (v2; new and modified files) and October 2025 (v3; new and modified files). File and folder names were not revised to indicate changes. See the change history section in the readme for more details. This data package contains the data and scripts used to develop the SPIT framework WT-WRF-Hydro (Water Tagging Weather Research and Forecasting Hydrologic) model and is associated with the following GitHub repository: https://github.com/zbutler33/SPIT-Framework. This data package contains five parent folders: (1) “Manipulated_outputs”, (2) “Metadata”, (3) “Observed”, (4) “Outputs”, and (5) “Scripts”. Each of these parent folders contains additional subfolders and files. Please see the FLMD (“v*_Butler_2024_WT_WRF_Hydro_flmd.csv”) for a list of all the files contained in this data package and descriptions for each. See the data dictionary (“v*_Butler_2024_WT_WRF_Hydro_dd.csv”) for definitions and units of all of the tabular (files ending in “.csv” and ".tsv") column headers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Source Attribution of Methane Enhancements at the Southern Great Plains Atmospheric Radiation Measurement Site Field Campaign Report

Tower data from 2017 to 2020 revealed that the U.S. Department of Energy’s Southern Great Plains (SGP) Atmospheric Radiation Measurements (ARM) observatory in Oklahoma exhibited significant CH 4 enhancements compared to other U.S. tower sites. On average, near-surface CH 4 at the ARM site was higher, had larger diurnal and seasonal variations, and had sharper near-surface increases at night compared to other tower sites in the U.S. A field campaign was conducted in June 2024 to identify and quantify the surrounding CH 4 emission sources by two researchers from the University of Oklahoma. During the campaign, a mobile measurement was conducted using a LI-COR 7810 tracer gas analyzer mounted on a vehicle during June 23-26, 2024 near the SGP site (within 7 km distance). We sampled CH 4 along predefined routes surrounding a suspected source, i.e., a large animal feeding operation (AFO) farm 5.8 km northwest of the SGP ARM site. Measurements were collected between 04:30 and 09:00 local time when the planetary boundary layer (PBL) was shallow, increasing sensitivity to surface emissions. Driving speed was kept below 5 mph to minimize pressure-related artifacts. CH 4 concentrations downwind of the AFO exceeded 6000 ppb on the early morning of June 26, corresponding to an emission rate of ~95 kg⋅hr -1 estimated using the mass balance method. Additional sources identified during the campaign include nearby open-range cattle and leaks from a nearby natural gas pipeline.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Multimodel simulation of vertical gas transfer in a temperate lake

In recent decades, several lake models of varying complexity have been developed and incorporated in numerical weather prediction systems and climate models. To foster enhanced forecasting ability and verification, improvement of these lake models remains essential. This especially applies to the limited simulation capabilities of biogeochemical processes in lakes and greenhouse gas exchanges with the atmosphere. Here we present multi-model simulations of physical variables and dissolved gas dynamics in a temperate lake (Harp Lake, Canada). The five models (ALBM, FLake, LAKE, LAKEoneD, MTCR-1) considered within this most recent round of the Lake Model Intercomparison Project (LakeMIP) all captured the seasonal temperature variability well. In contrast, none of the models is able to reproduce the exact dates of ice-cover and ice-off, leading to considerable errors in the simulation of eddy diffusivity around those dates.We then conducted an additional modeling experiment with a diffusing passive tracer to isolate the effect of the eddy diffusivity on gas concentration. Remarkably, sophisticated k -" models do not demonstrate a significant difference in the vertical diffusion of a passive tracer compared to models with much simpler turbulence closures. All models simulate less intensive spring overturn compared to autumn. Reduced mixing in the models consequently leads to the accumulation of the passive tracer distribution in the water column. The lake models with a comprehensive biogeochemical module, such as ALBM and LAKE, predict dissolved oxygen dynamics adequate to the observed data. However, for the surface carbon dioxide concentration the correlation between modeled (ALBM, LAKE) and observed data is weak (~0.3). Overall our results indicate the need to improve the representation of physical and biogeochemical processes in lake models, thereby contributing to enhanced weather prediction and climate projection capabilities.

Guseva, Sofya↗
DOE OSTI2022

Variability in observed stable water isotopes in snowpack across a mountainous watershed in Colorado

In this study, isotopic information from 81 snowpits was collected over a 5-year period in a large, Colorado watershed. Data spans gradients in elevation, aspect, vegetation, and seasonal climate. They are combined with overlapping campaigns for water isotopes in precipitation and snowmelt, and a land-surface model for detailed estimates of snowfall and climate at sample locations. Snowfall isotopic inputs, describe the majority of δ18O snowpack variability. Aspect is a secondary control, with slightly more enriched conditions on east and north facing slopes. This is attributed to preservation of seasonally enriched snowfall and vapour loss in the early winter. Sublimation, expressed by decreases in snowpack d-excess in comparison to snowfall contributions, increases at low elevation and when seasonal temperature and solar radiation are high. At peak snow accumulation, post-depositional fractionation appears to occur in the top 25 ± 14% of the snowpack due to melt-freeze redistribution of lighter isotopes deeper into the snowpack and vapour loss to the atmosphere during intermittent periods of low relative humidity and high windspeed. Relative depth of fractionation increases when winter daytime temperatures are high and winter precipitation is low. Once isothermal, snowpack isotopic homogenization and enrichment was observed with initial snowmelt isotopically depleted in comparison to snowpack and enriching over time. The rate of δ18O increase (d-excess decrease) in snowmelt was 0.02‰ per day per 100-m elevation loss. Isotopic data suggests elevation dictates snowpack and snowmelt evolution by controlling early snow persistence (or absence), isotopic lapse rates in precipitation and the ratio of energy to snow availability. Hydrologic tracer studies using stable water isotopes in basins of large topographic relief will require adjustment for these elevational controls to properly constrain stream water sourcing from snowmelt.

54 ENVIRONMENTAL SCIENCES↗