Results for “environmental tracers” — Page 19

DOE OSTI2023

Assessment of Baseline Monitoring Data for the Remote-Handled Low-Level Waste Disposal Facility at Idaho National Laboratory

The purpose of this document is to summarize environmental monitoring data collected at the Remote-Handled Low-Level Waste (RHLLW) Disposal Facility during the first four years of facility operations (FY 2019 through FY 2022). This summary provides a “baseline” condition for the facility, which can be used to distinguish contaminant releases from the RHLLW Disposal Facility from pre existing contamination as well as potential future releases from other sources (i.e., upgradient aquifer sources), as specified in the facility monitoring plan. Sufficient data has been collected over the first four years of operation of the RHLLW Disposal Facility to establish baseline conditions for future compliance monitoring of aquifer wells and performance monitoring of vadose zone lysimeters. Except for seven elevated tritium measurements from lysimeter HFEF-South believed to have been impacted by waste disposals, all data was deemed appropriate for establishing baseline concentrations. Monitoring of the aquifer detected indicator analytes gross alpha and gross beta and target analytes tritium and C-14. There were a few C-14 detections prior to an increase in the required detection level (RDL) after the first round of sampling. I-129 and Tc-99 were not detected above RDLs in aquifer samples. Lysimeter samples detected indicator analytes gross alpha and gross beta and target analyte tritium. Target analytes C-14, I-129, and Tc-99 were not detected above RDLs in lysimeter samples. Except for tritium, baseline concentrations of indicator and target analytes measured during the first four years of RHLLW facility operations are expected to represent conditions through the expected 20 year operating period of the facility. Elevated tritium in the aquifer, a result of past discharges of tritium at upgradient facilities, showed a decline over the first four years of operations, and levels are expected to continue to decline with time as a result of dilution and decay. Current tritium levels in lysimeter samples may be elevated due to tritium in the water applied during construction and infiltration testing of the facility. As a result, tritium concentrations in lysimeters may decline with time as a result of dilution and decay. Statistical measures (i.e., mean, standard deviation, and 99% upper confidence level [UCL]) were calculated for all detected analytes using the four years of concentration data. These measures will be used to evaluate future monitoring results and to demonstrate the facility is performing, or not performing, as established in the facility performance assessment (PA). Future monitoring results will also be combined with the baseline data in this report to further establish temporal trends and natural variability in the measurements. Based on the data collected during the baseline period, it is recommended the gross alpha action level for performance monitoring increase from 10 to 20 pCi/L due to several baseline measurements exceeding the preliminary action level of 10. An action level of 20 pCi/L is slightly greater than the 99% UCL, is protective of the aquifer, and would reduce unnecessary sampling. It is also recommended that tritium be added to the lysimeter analyte list for post-baseline period monitoring. Tritium, while not a dose concern, is a good tracer that can provide valuable information on water flow in and around the RHLLW Disposal Facility. Recommendations for a tritium action level are provided if considered necessary.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Modeling Functional Organic Chemistry in Arctic Rivers: An Idealized Siberian System

Rivers of the Arctic will become ever more important for the global climate, since they carry a majority of continental dissolved organic carbon flux into the rapidly changing polar ocean. Aqueous organics comprise a wide array of functional groups, several of which are likely to impact coastal and open water biophysical properties. Light attenuation, interfacial films, aerosol formation, gas release and momentum exchange can all be cited. We performed Lagrangian kinetic modeling for the evolution of riverine organic chemistry as the molecules in question make their way from the highlands to Arctic outlets. Classes as diverse as the proteins, sugars, lipids, re-condensates, humics, bio-tracers and small volatiles are all included. Our reduced framework constitutes an idealized northward flow driving a major hydrological discharge rate and primarily representing the Russian Lena. Mountainous, high solute and tundra sources are all simulated, and they meet up at several points between soil and delta process reactors. Turnover rates are parameterized beginning with extrapolated coastal values imposed along a limited tributary network, with connections between different terrestrial sub-ecologies. Temporal variation of our total dissolved matter most closely resembles the observations when we focus on the restricted removal and low initial carbon loads, suggesting relatively slow transformation along the water course. Thus, channel combinations and mixing must play a dominant role. Nevertheless, microbial and photochemical losses help determine the final concentrations for most species. Chemical evolution is distinct for the various functionalities, with special contributions from pre- and post-reactivity in soil and delta waters. Several functions are combined linearly to represent the collective chromophoric dissolved matter, characterized here by its absorption. Tributaries carry the signature of lignin phenols to segregate tundra versus taiga sources, and special attention is paid to the early then marine behaviors of low molecular weight volatiles. Heteropolycondensates comprise the largest percentage of reactive carbon in our simulations due to recombination/accumulation, and they tend to be preeminent at the mouth. Outlet concentrations of individual structures such as amino acids and absorbers lie above threshold values for biophysical influence, on the monolayer and light attenuation. The extent of coastal spreading is examined through targeted regional box modeling, relying on salinity and color for calibration. In some cases, plumes reach the scale of peripheral arctic seas, and amplification is expected during upcoming decades. Conclusions are mapped from the Lena to other boreal discharges, and future research questions are outlined regarding the bonding type versus mass release as permafrost degrades. Dynamic aqueous organic coupling is recommended for polar system models, from headwaters to coastal diluent.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Water vapor isotopic composition during TRACER IOP

The TRACER-Iso project took place from June 1 to September 30, 2022 as part of the intensive observation period (IOP) of the Tracking Aerosol Convection Interactions (TRACER) project at the main project site at the La Porte, Texas, municipal airport. The project sought to use measurements of the stable isotopic composition of water vapor to better understand cloud-aerosol interactions in deep convection in a polluted urban setting. Water vapor isotopologue and humidity measurements were determined using a Picarro L2130 Water Vapor Isotope Analyzer. The instrument was deployed in the guest van and consists of three main components: a water vapor isotope analyzer, and a Los Gatos Research Water Vapor Isotope Standard Source (WVISS), and a Los Gatos Research Dry Air Source. Ambient air samples were collected through an inlet located about one meter above roof of the shipping container. Samples were then delivered to the analyzer through teflon tubing by use of an external pump to minimize the transport time between the inlet and analyzer. The inlet tubing was surrounded by a Briskheat heat trace and insulating material to ensure ambient air samples did not fall below the dew point and result in condensation. The analyzer uses cavity ringdown spectroscopy to report isotopic ratios of ambient air samples at a frequency of 1 Hz. The calibration unit of the instrument was used in conjunction with the dry air source to periodically measure the δ values of standard waters. It uses a nebulizer to push small water droplets into a hot chamber that vaporizes the water without fractionation. This vapor was then transported to the analyzer using a built-in compressor and the dry air source, which allows each standard with known δ values to be measured at a wide range of humidity values for post-measurement calibration of ambient air samples. There are multiple sources of uncertainty introduced at different stages during the data collection and processing. These include (1) instrument precision, (2) uncertainty in the secondary standards, (3) humidity-correction uncertainty, and (4) VSMOW-SLAP calibration uncertainty. Uncertainty from each step is propagated in quadrature to calculate a total uncertainty of each isotopologue. Uncertainty was determined to be 2.3‰ for δD and 0.4‰ for δ18O.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Investigating Secondary Aerosol Processes in the Amazon through Molecular-level Characterization of Semi-Volatile Organics

In areas where biogenic emissions are oxidized in the presence of anthropogenic pollutants, it has become increasingly apparent that secondary organic aerosol (SOA) formation from biogenic volatile organic compounds (BVOCs) is substantially enhanced. The Amazon forest is the dominant source of BVOCs globally, and the forest is rapidly being converted to urban and agricultural uses. We participated in a collaborative field study located in the Amazon region around Manaus, Brazil, in 2014 (Green Ocean Amazon; GoAmazon). This study was designed to comprehensively examine how BVOC emissions (specifically, the most highly emitted non-methane hydrocarbon, isoprene, and lesser studied sesquiterpenes) and their interaction with anthropogenic pollution (i.e., nitrogen oxides, sulfur dioxide, and black carbon) alter the atmosphere’s oxidative capacity, influence secondary aerosol formation, atmospheric composition and, ultimately, affect the earth’s radiation balance and climate. We provide the first hourly, in-situ measurement of 30 sesquiterpenes and 4 diterpenes, roughly estimating that sesquiterpene oxidation contributes to 10-14% of ozone reactive loss by terpenes and at least 0.4–5% (median 1 %) of total submicron OA mass. However, this is likely a low-end estimate, as evidence for additional unaccounted sesquiterpenes and their oxidation products clearly exists. We also provide new perspectives on sulfate, NO x , and particle acidity influencing isoprene-derived SOA, comparing the central Amazon environment with Southeastern U.S.A. summer, representing clean to polluted conditions, respectively. We find that summed concentrations of isoprene-derived SOA tracers correlated with particulate sulfate spanning three orders of magnitude, suggesting that 1 μg m -3 reduction in sulfate corresponds with at least ~0.5 μg m -3 reduction in isoprene-derived SOA. We also find that SOA mass derived from isoprene oxidation in the presence of NO x is primarily comprised of aerosol sulfate bound with isoprene oxidation products (i.e. organosulfates), ~97% in the Amazon and ~55% in the Southeastern U.S. We infer under natural conditions in high isoprene emission regions, preindustrial aerosol sulfate was almost exclusively isoprene-derived organosulfates, which are traditionally thought as representative of anthropogenic influence. We further report the first field observations showing that particle acidity impacts the distribution of isoprene oxidation tracers in the gas and particle phases, providing physicochemical insight into isoprene SOA formation chemistry. Coupled with other co-located measurements by the DOE Atmospheric Radiation Measurement (ARM) team and collaborating PIs at our measurement site (called T3), on the G1 aircraft, and at additional field sites (T0, T1, T2) this data set has been used to understand emission, oxidation, and particle formation pathways at the molecular level, to elucidate how these pathways change when influenced by anthropogenic sources, and to evaluate their importance for the atmospheric budget of aerosols and the earth’s radiation balance on regional scales. Recommendations: Further constraint of the role of sesquiterpenes on ozone (O 3 ) reactivity and SOA formation is limited by the availability of authentic standards and synthesized compounds of sesquiterpenes and their oxidation products. Future research efforts should focus on making such standards available via custom synthesis to allow for accurate quantification and focused oxidation experiments that will fully elucidate the chemistry of these compounds in the atmosphere. Further, reductions in aerosol sulfate (via SO 2 emissions control) continues to be one of the most effective methods to reduce SOA formation from isoprene, but the concerted role of ammonia (NH 3 ) gas emissions from (agricultural) industry and other sources requires more investigation to better understand the role of aerosol acidity in SOA formation and potential controls. As GoAmazon particles can be more acidic compared to areas with greater anthropogenic NH 3 emissions, this promotes relatively greater organosulfate (OS) formation from isoprene even in the presence of NO x . Because OS and inorganic sulfate differ in water uptake properties, this implies preindustrial aerosol sulfate (albeit less abundant) may have been less reflective than current earth system models assume. Further research should investigate the radiative impacts of organic vs inorganic sulfate in aerosols to improve model representation. Finally, future work (currently underway) should include constraint of the contributions of monoterpene (C 10 H 16 ) BVOCs and biomass burning VOCs as precursors to SOA formation in the region as well as their impacts on radiative forcing in the climate system.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Findings on Subtask 3.1 - Bakken Rich Gas Enhanced Oil Recovery Project

Total in-place oil for the Bakken petroleum system (BPS) (which includes the Bakken and Three Forks Formations) has been estimated to be 600 billion barrels (bbl). However, BPS wells have decline rates as high as 85% over the first 3 years of their lives, and primary recovery factors typically range from 3% to 10% of original oil in place. Given the low initial recovery rates, even small incremental productivity improvements could dramatically increase technically recoverable oil in the BPS. One potential solution is enhanced oil recovery (EOR) using gas injection, such as carbon dioxide (CO2) or hydrocarbon (HC) gases. While commonly used in conventional reservoirs, CO2 EOR in unconventional tight oil reservoirs has been limited to pilot tests. EOR using rich gas (mixture of methane, ethane, and propane) has also been employed in numerous pilots in several unconventional plays and has recently been successfully applied in the Eagle Ford play. If successful, large-scale gas-based EOR in the BPS could dramatically increase oil productivity and recovery factors and extend the life of the play for decades. While CO2 may be a technically suitable working fluid for EOR in the BPS, supplies are limited and costs for using CO2 in EOR pilots are prohibitively high. Meanwhile, produced gas flaring has presented challenges for BPS operators in North Dakota. Analysis conducted by the North Dakota Pipeline Authority indicates that the current gas-gathering infrastructure in North Dakota is insufficient to accommodate all of the associated gas that is produced from the BPS. The geographically isolated location of North Dakota relative to large natural gas markets, combined with suppressed natural gas prices, has made it economically challenging for industry to invest capital in expanding gas-gathering infrastructure in the state. These circumstances led to a research program conducted by the Energy & Environmental Research Center (EERC) in partnership with Liberty Resources Management Company LLC (LR) to examine the potential to use rich gas injection for EOR and mitigate flaring. A rich gas EOR pilot test was designed and executed by LR at its Stomping Horse development area in Williams County, North Dakota. From July 2018 through May 2019, a total of 160 million standard cubic feet (MMscf) of rich produced gas was injected into the BPS using five different wells in a sequential injection strategy. LR’s Leon–Gohrick drill spacing unit (DSU) was used as the test site. Regulatory oversight was provided by the North Dakota Industrial Commission (NDIC). Technical support was provided by the EERC through a series of laboratory, modeling, and field-based activities, and additional post-pilot research activities incorporated learnings from the test, developed new laboratory data, improved fracture modeling methods, and developed machine learning and big data analytics. The results from the Stomping Horse rich gas EOR pilot activities indicate that developing an effective, economical EOR approach for the BPS will require more field tests. Another key lesson learned from the Stomping Horse tests is that detailed pre- and posttest data on reservoir conditions and fluids production are essential. Robust reservoir characterization provides information that is crucial to creating realistic geomodels and conducting valid dynamic simulations of potential EOR scenarios. A detailed understanding of the completions and production history of offset wells is also necessary for valid test result interpretations. This knowledge is essential to designing the operational parameters of injectivity tests and interpreting the results. A conformance control strategy is also essential to success. Laboratory-based examinations of rich gas interactions with reservoir fluids and rocks were conducted, with an emphasis on determining the ability to mobilize oil in the tight reservoir rocks and shales of the BPS. Injection fluid composition was shown to have a positive impact on reducing reservoir oil minimum miscibility pressure (MMP), reducing interfacial tension (IFT), and altering wettability. IFT and contact angle measurements demonstrated that wettability can be altered in the presence of rich gas, suggesting the potential to improve oil recovery. Iterative modeling of surface infrastructure and reservoir performance using data generated by the various project activities was conducted. A geologic model of the Stomping Horse area was built; history-matched oil, gas, and water production was used in simulations of various EOR scenarios. Early programmatic modeling results were used to support LR’s design and operation of the EOR pilot and to provide insight regarding optimization of future commercial-scale BPS EOR design and operations. Post-pilot modeling focused on alternative methods of understanding complex fracture networks and accelerating simulation time. These led to improved simulation run times and provide excellent history-matching results. Several of these iterative models were used as the bases for developing algorithms into machine learning and big data analytics. History matching in reservoir simulation is time-consuming and computer processing-intensive. Machine learning algorithms were created, and an automated history-matching tool was developed. A large set of synthetic reservoir simulations were created to generate well responses (oil, gas, and water production, well bottomhole pressure [BHP], and tracer or propane breakthrough) for a set of EOR operating parameters that included offset well status (open or closed), injectate (rich gas or propane), injection rate, and injection well BHP. A user interface was developed to provide real-time visualization. Machine learning-based models were developed to provide rapid forecasting of well performance given a set of user-defined EOR operating parameters. These predictive models allow the user to modify the offset well status, injection rate, and injection well BHP and rapidly forecast future production performance. The combination of real-time visualization tools with real-time forecasting tools provides a framework for real-time control—operational changes that the EOR site operator can enact (e.g., changing gas injection rates) to affect the observed performance and potentially improve the EOR outcome. There is great reason to be optimistic about the future of EOR in the Bakken. The results of the laboratory studies suggest significant potential for high rates of oil mobilization using produced field gas injection under the right conditions. The results of the lab studies, combined with rigorous statistical analysis of well production data and associated modeling efforts, confirm the notion that fluid mobility within the reservoir is controlled by fractures. As more knowledge is gained about the nature and distribution of fracture networks in the Bakken, the industry will be in a better position to predict and, ultimately, influence fluid mobility. New field tests are necessary to develop a more complete understanding of those conditions. Thoughtful and creatively engineered field tests within a well-characterized geologic setting will yield the fundamental knowledge needed to take Bakken oil production to the next level. This subtask was cofunded through the EERC–U.S. Department of Energy Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Industrial Commission’s Oil and Gas Research Program and Computer Modelling Group.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2024

Chemical properties and single-particle mixing state of soot aerosol in Houston during the TRACER campaign

Abstract. A high-resolution soot particle aerosol mass spectrometer (SP-AMS) was used to selectively measure refractory black carbon (rBC) and its associated coating material using both the ensemble size-resolved mass spectral mode and the event trigger single particle (ETSP) mode in Houston, Texas, in summer 2022. This study was conducted as part of the Department of Energy Atmospheric Radiation Measurement (ARM) program's TRacking Aerosol Convection interactions ExpeRiment (TRACER) field campaign. The study revealed an average (±1σ) rBC concentration of 103 ± 176 ng m−3. Additionally, the coatings on the BC particles were primarily composed of organics (59 %; 219 ± 260 ng m−3) and sulfate (26 %; 94 ± 55 ng m−3). Positive matrix factorization (PMF) analysis of the ensemble mass spectra of BC-containing particles resolved four distinct types of soot aerosol, including an oxidized organic aerosol (OOABC,PMF) factor associated with processed primary organic aerosol, an inorganic sulfate factor (SO4,BC,PMF), an oxidized rBC factor (O-BCPMF), and a mixed mineral dust–biomass burning aerosol factor with significant contribution from potassium (K-BBBC,PMF). Additionally, K-means clustering analysis of the single-particle mass spectra identified eight different clusters, including soot particles enriched in hydrocarbon-like organic aerosol (HOABC,ETSP), sulfate (SO4,BC,ETSP), two types of rBC, OOA (OOABC,ETSP), chloride (ClBC,ETSP), and nitrate (NO3,BC,ETSP). The single-particle measurements demonstrate substantial variation in BC coating thickness with coating-to-rBC mass ratios ranging from 0.1 to 100. The mixing state index (χ), which denotes the degree of homogeneity of the soot aerosol, varied from 4 % to 94 % with a median of 40 %, indicating that the aerosol population lies in between internal and external mixing but has large temporal and source type variability. In addition, a significant fraction of BC-containing particles, a majority enriched with oxidized organics and sulfate, exhibit sufficiently high κ values and diameters conducive to activation as cloud nuclei under atmospherically relevant supersaturation conditions. This finding bears significance in comprehending the aging processes of rBC-containing particles and their activation into cloud droplets. Our analysis highlights the complex nature of soot aerosol and underscores the need to comprehend its variability across different environments for accurate assessment of climate change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Characterizing Greater Houston's Aerosol by Air Mass During TRACER

During the TRacking Aerosol and Convection interaction ExpeRiment (TRACER), a suite of aerosol and cloud measurements were made using the Texas A&M Rapid Onsite Atmospheric Measurement Van (ROAM-V) mobile instrument platform. Joint ROAM-V/radiosonde deployments focused on sampling polluted marine and continental air masses to understand summertime convection along and across the sea-breeze front as it propagated through Houston. The polluted marine air mass at Seawolf Park, Galveston, TX was defined by two characteristic aerosol populations. Although the background total particle concentrations averaged 2,500 cm −3 , the air mass was also influenced by frequent but irregular periods of ship emission with significantly higher aerosol concentrations at times exceeding 34,000 cm −3 . Ship emission influenced periods, typically lasting <10 min, contained smaller, less hygroscopic particles resulting in a 69% relative decline in cloud condensation nuclei (CCN) activated fraction at 1% supersaturation compared to background periods. Measurements in continental air masses northwest of Houston revealed an average aerosol concentration of 5,208 cm −3 with a κ value near 0.1 reasonably describing the CCN population, independent of particle size. In continental air masses, substantial differences in particle size and CCN activation over small distances (<42 miles between sites) suggest considerable site-to-site variability in addition to the expected day-to-day differences. This small-scale variability makes it difficult to generalize continental air mass aerosol properties. Both coastal and inland locations had effective ice nucleating particles, but inland deployments observed the warmest nucleation temperature at −15.6°C compared to −17.8°C close to the coast. These measurements can reduce uncertainties in regional convection allowing models to improve understanding of aerosol-cloud interactions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

SMPS data from TAMU TRACER campaign in the Houston TX region from July to September 2022

During TRACER, the Texas A&M Rapid Onsite Atmospheric Measurements Van (ROAM-V) was deployed to capture airmasses behind (maritime) and ahead (continental) of the passage of the sea-breeze front through Houston. On select sampling days, ROAM-V sampled in the morning/mid-day on the coast and then transited to a second inland site for the afternoon/evening. The suite of instruments deployed on ROAM-V included a Condensation Particle Counter (CPC; GRIMM Model 5.403 CPC), Scanning Mobility Particle Sizer (SMPS; TSI 3750 detector, TSI 3082 classifier, TSI 3088 neutralizer, TSI 3081A Differential Mobility Analyzer), Cloud Condensation Nuclei counter (Droplet Measurement Technologies CCN Counter), micro pulse lidar (Droplet Measurement Technologies Micro Pulse LiDAR (miniMPL)), and a Davis Rotating Uniform size-cut Monitor (DRUM; DRUMAir 4-DRUM). Before sampling at each location, the latitude and longitude were recorded using the GPS on the phone application “My Altitude”. Onboard the ROAM-V, aerosol samples are drawn through a shared isokinetic inlet at a flow rate ranging from 3.5 to 7.0 LPM. A portion of this flow, 1.0 LPM, is directed through TSI's 0.071 cm impactor attached to the classifier of the SMPS setup. The SMPS’s DMA and CPC are connected through a 20.3 cm length of 0.48 cm diameter tubing. Measured SMPS size distributions were used to calculate size-dependent particle losses for each SMPS scan. Particle losses from diffusion (based on Kesten, 1991 and Gormley, 1949) and inertial impaction in 90-degree bends (based on Aerosol Measurement, 2011 and Crane, 1977) were included in the loss calculation. This data was collected for ARM Field Campaign AFC07055 and supported by DOE ASR grant DE-SC0021047. For any further questions, please feel free to contact the instrument PI, Sarah D. Brooks, sbrooks@tamu.edu. Kesten et. al. Calibration of a TSI Model 3025 Ultrafine Condensation Particle Counter. Aerosol Science and Technology, 15:2, 107-111, 1991. Gormley et. al. Diffusion from a Stream Flowing through a Cylindrical Tube. Proceedings of the Royal Irish Academy, Vol 52, 163-169, 1948. Aerosol Measurement: Principles, Techniques, and Applications, Third Edition. John Wiley & Sons, Inc, 2011. Crane et. al. Inertial Deposition of Particles in a Bent Pipe. Journal of Aerosol Science, Vol 8, 161-170, 1977.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Aerosol-deep convection interaction based on joint cell-thermal tracking in Large Eddy Simulations during the TRACER campaign

In cumulus clouds, aerosol concentrations control cloud droplet concentrations, modifying cloud radiative properties, precipitation processes, and cloud electrification. However, mechanisms of aerosol-deep convection interactions are not well understood due to complex cloud dynamics and microphysics. We investigate the interaction of aerosols with isolated deep convection using Large Eddy Simulations of two cases during the TRacking Aerosol Convection interactions ExpeRiment (TRACER) near Houston, Texas, using a joint cell-thermal tracking algorithm. Cumulus thermals are droplet generators, since supersaturation and droplet nucleation coincide with thermal centers, where the strongest updrafts occur. Primary ice crystal formation does not take place inside thermals, but at layers where previous thermals detrained moisture. As subsequent thermals containing supercooled droplets penetrate these layers, hail and graupel form at or near these thermals. Higher aerosol concentrations result in higher droplet concentrations that suppress drizzle, delay warm rain processes, and transport more moisture aloft. This increases snow and ice amount, as well as graupel and hail, leading to more lightning. Polluted thermals initiate at slightly higher altitudes, and are slightly larger and faster, suggesting a weak invigoration. We also find more thermals per cell, but fewer isolated cells, since convection is more aggregated and intense, especially near the end of the 24 h simulation. Non-linear mesoscale feedback likely triggered by temperature and moisture responses to aerosol-thermal interactions causes the aggregation. Time-lagged aerosol-reinitialization experiments show that the mesoscale response is the predominant forcing for the invigoration. These changes happen within one day, on a smaller scale than previously suggested.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Bayesian Approach to Estimation of Water Table Elevations Using Historical Rasters as Prior Information 2019 - 20430

In cases of complex but only partially known geology and a lack of spatial control in observation well locations, water table elevation estimation is very challenging. In some cases, auxiliary information, such as observations of the movement of tracers, operation of injection and extraction wells, and calibration of groundwater models against the historical elevation data, can be combined with expert judgement to estimate flow directions in areas of sparse data and to aid in the production of more reliable contour maps (and associated rasters) than could be produced by relying on sparse well elevation data alone. Given a historical sequence of these raster maps, the question arises how to automate, to the extent possible, the process of producing new raster maps to reflect data from previous times, the current data and the operation of expert judgement. One solution is to adopt a Bayesian point of view and to regard the historical well elevation data, auxiliary information and historical raster maps as prior information. The well elevations for water table wells, as well as those for injection/extraction wells and the data associated with other relevant variables, can be viewed as predictors for the raster surface. From this prior information, we can, conditional on the values of these predictors for a new time period, compute an expected value map and a standard deviation map for the new raster. These then can be taken to specify a prior predictive distribution for the pixels in the new raster map. Then we condition the pixels, corresponding to water level observation wells within the raster, on the observed values in those wells (which in general will differ from the regression estimate) for the new time period. Given the smoothness of the water table surface, we then smooth the surface of deviations from the mean surface, based on the variograms of the historical rasters, and add this smoothed surface to the regression mean surface. The error structure of the produced raster map is defined by the regression error structure and the error due to smoothing based on the estimated variograms. This methodology has been developed and is being further refined for groundwater monitoring and remediation at LANL. It is a very flexible method that can also be applied with a variety of other predictors applied to model the water level wells in the area of interest over the historical record. The smoothness of the spatial process and its possible evolution over time can then be estimated from the residuals from this regression. This can be augmented by expert hydrogeological opinion based on site topography and hydrogeology. (authors)

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

The sources and diurnal variations of submicron aerosols in a coastal–rural environment near Houston, US

Aerosol properties were characterized at a rural site southwest of Houston from May to September 2022 during the intensive operation periods (IOPs) of the Tracking Aerosol Convection Interactions ExpeRiment (TRACER). Backward trajectory analysis reveals three major air mass types: marine air mass from the Gulf, urban air mass influenced by urban emissions, and regional air mass. Marine aerosols typically show a bimodal size distribution and have the lowest particle number and mass concentrations of PM 1 (particulate matter with an aerodynamic diameter of less than 1 µm), while aerosols from air masses strongly influenced by urban emissions exhibit the highest concentrations. Organic aerosol (OA) accounts for more than 50 % of PM 1 for urban and regional air masses, whereas sulfate is comparable to OA in marine air masses. Positive matrix factorization (PMF) analysis of aerosol mass spectra identifies 6 OA factors: hydrocarbon-like OA (HOA), OA from the oxidation of monoterpenes (MT-SOA), OA from the reactive uptake of isoprene epoxydiols by acidic sulfate particles (isoprene-SOA), oxygenated OA arising from shipping emissions (shipping-OOA), and two oxygenated OA factors with high O : C ratios (OOA1 and OOA2). OOA2 has the highest O : C ratio and exhibits elevated mass concentration in the afternoon. Similar diurnal variation of highly oxidized OA factors was commonly observed in the Houston area during previous studies and attributed to the SOA formation by photochemistry and mixing from aloft. Here, using air mass backward trajectories and a 1-D box model, we show the diurnal trend of OOA2 mass concentration is instead driven by changes in air mass arriving at the rural site. The air mass changes are likely caused by the shift between land breezes and sea/bay breezes. Within the same air mass type (e.g., either urban or marine air mass), OOA2 mass concentration is largely independent of wind direction and shows essentially no diurnal variation, suggesting OOA2 is related to aged OA with minimal influence by local emissions. This study helps identify the major sources of OA in the Houston region and highlights the impacts of both atmospheric chemistry and meteorology on aerosol properties in the coastal–rural environment.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Surface-to-Interior Transport Timescales and Ventilation Patterns in a Time-Dependent Circulation Driven by Sustained Climate Warming

Abstract The effect of climate warming in response to rising atmospheric CO 2 on the ventilation of the ocean remains uncertain. Here we make theoretical advances in elucidating the relationship between ideal age and transit time distribution (TTD) in a time-dependent flow. Subsequently, we develop an offline tracer-transport model to characterize the ventilation patterns and time scales in the time-evolving circulation for the 1850–2300 period as simulated with the Community Earth System Model version 1 (CESMv1) under a business-as-usual warming scenario. We found that by 2300 2.1% less water originates from the high-latitude deep water formation regions (both hemispheres) compared to 1850. In compensation, there is an increase in the water originating from the subantarctic. We also found that slowing meridional overturning circulation causes a gradual increase in mean age during the 1850–2300 period, with a globally averaged mean-age increase of ∼110 years in 2300. Where and when the water will be re-exposed to the atmosphere depends on the post-2300 circulation. For example, if we assume that the circulation persists in its year-2300 state (scenario 1), the mean interior-to-surface transit time in year 1850 is ∼1140 years. In contrast, if we assume that the circulation abruptly recovers to its year-1850 state (scenario 2), the mean interior-to-surface transit time in 1850 is only ∼740 years. By 2300, these differences become even larger; in scenario 1, the mean interior-to-surface transit time increases by ∼200 years, whereas scenario 2 decreases by ∼80 years. The dependence of interior-to-surface transit time on the future ocean circulation produces an additional unavoidable uncertainty in the long-term durability of marine carbon dioxide removal strategies. Significance Statement The ocean’s circulation, when altered by climate warming, can affect its capacity to absorb heat and CO 2 , which are crucial for the global climate. In our study, we investigated how global warming, caused by rising CO 2 levels, might impact the ocean circulation—the way water moves from deep ocean to the surface and vice versa. We discovered that by 2300, if we continue on our current warming trajectory, the origins of water within the ocean will shift, with less coming from deep, cold zones near the poles and more from subantarctic regions. As a result, deep water will take longer time before it resurfaces than shallow water. How quickly this water travels from deep regions to the surface could change, depending on the state of future ocean circulation. If the circulation remains as predicted in 2300, this journey will take longer. Conversely, if it reverts to the pattern in 1850, the process will be quicker. This variability introduces added uncertainty to strategies aimed at mitigating climate change by storing CO 2 in the ocean. Our work highlights the intricate ways in which climate change can influence our oceans, potentially affecting our plans to mitigate global warming.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

A dendritic strontium river isoscape for fisheries applications in the Sacramento River basin, California, USA

Objective Understanding the origins and movements of fish is fundamental to effective conservation and fisheries management. Strontium isotope ratios ( 87 Sr/ 86 Sr) in otoliths provide a powerful tracer of natal origin and migratory pathways. However, existing 87 Sr/ 86 Sr isoscapes for the Sacramento River basin, an ecosystem that supports ecologically and economically important salmon populations, rely on discrete classification approaches that overlook unsampled habitats and do not incorporate spatial uncertainty. Our objective was to develop a continuous, network-explicit 87 Sr/ 86 Sr isoscape with quantified uncertainty to fill in data gaps and enable probabilistic assignments of fish origin and movement. Methods We used river water 87 Sr/ 86 Sr data from 106 sites (1997–2021) to develop spatial stream network models that use dendritic connectivity and watershed characteristics (lithology, bedrock age, and land cover) to predict river water 87 Sr/ 86 Sr throughout the basin. Models were fitted using maximum and restricted likelihood and were evaluated via Akaike’s information criterion and leave-one-out cross validation. We produced both historical (pre-dam) and present-day (below-dam) isoscapes, delineated uncertainty-informed isotopic ranges using k -means clustering, and applied a proof-of-concept Bayesian assignment to estimate natal origins and early rearing habitats for two endangered winter-run Chinook Salmon Oncorhynchus tshawytscha. Results Cross validation indicated strong performance of the 87 Sr/ 86 Sr model (leave-one-out cross validation: R 2 = 0.91; root mean square error = 0.0005). Uncertainty-informed clustering identified 19 isotopic “suites” (reaches with indistinguishable 87 Sr/ 86 Sr values) in present-day anadromous habitats and 25 suites in the historical network. Example natal and early rearing assignments included predictions that challenged expectations for juvenile salmon migration based on predicted river 87 Sr/ 86 Sr compositions. Conclusions This study developed a continuous, network-explicit 87 Sr/ 86 Sr isoscape that integrates existing river data to predict 87 Sr/ 86 Sr in unsampled reaches and the likely achievable range and resolution of otolith-based origin and life history inference. The resulting river isoscape provides a valuable tool to predict salmon movements and identify habitats supporting their survival and growth that otherwise might remain undetected. Coupling these predictions with complementary approaches that ground-truth juvenile presence (e.g., targeted fish surveys) represents an important step toward science-informed restoration and management of critical habitats throughout the Sacramento River basin.

Environmental sciences↗

DOE OSTI2020

A Comparison of Groundwater Sampling Technologies, Including Passive Diffusion Sampling, for Radionuclide Contamination - 20205

Using traditional high-flow purge methods for long-term water quality monitoring of deep groundwater wells can be expensive, affect contaminant migration, and produce excessive volumes of discharge water that can be difficult to manage. The use of low-flow pumping methods and depth discrete bailers (DDBs) can reduce the cost of sampling deep groundwater wells. In general, using different pumping methods to obtain reproducible and representative groundwater can be challenging. Passive diffusion samplers (PDSs) have successfully been used in long-term monitoring for volatile organic compounds, major ions, and trace-elements, but application has been limited for stable and radioactive isotopes. Beginning in 2018, the United States Geological Survey (USGS) conducted three sampling events to test the ability of PDSs and DDBs to obtain reproducible and representative groundwater samples. The sampled well is completed in a regional, permeable carbonate aquifer and is one in a set of wells that have historically been used for tritium tracer testing. All samples were obtained at 180 m (590 ft) below land surface (bls) in a 13.97 cm (5.5 inch) uncased well that has a total depth of 202 m (662 ft) bls. The first sampling event deployed a regenerated cellulose dialysis membrane (RCDM) PDS for 14 days with deionized water as the blank. The second sampling event deployed a RCDM for 27 days also with deionized water as the blank. The third sampling event deployed a Dual Membrane (DM) PDS for 65 days using a blank of tritium-dead carbonate water. The DM PDS was used to assess the effect of longer term deployment on tritium concentrations and address whether or not the PDSs reached equilibrium with ambient groundwater. The day after each of the three passive samplers were retrieved a DDB was used to obtain discrete non-integrated groundwater samples. For each DDB sampling day, the bailer was lowered into the well 10 consecutive times to determine if the water chemistry changed from the first to the last bailed sample. Quality assurance samples including blanks and duplicates were obtained during all three sampling events. All blank waters had tritium concentrations less than 21±33 pCi/L. Major ion (e.g. calcium, chloride, sodium, and sulfate) results were compared between all samples obtained with RCDM and DDB. Major ion concentrations showed a coefficient of variation of less than 6% between all RCDM and DDB samples; however, the coefficient of variation between the different deployment times and the two different methods for trace-element concentrations was much larger, particularly for manganese (82%), lead (41%), and zinc (33%). Stable isotope values were compared between the RCDM and DDB samples. The DDB sample results all fell within analytical uncertainty and were considered representative of the formation groundwater. The stable isotope values from the RCDM samples indicated that a longer deployment time was necessary to gain equilibrium and to obtain representative groundwater samples. Tritium results from groundwater samples obtained from the RCDM, DM, and DDB indicate that groundwaters obtained with PDSs produced tritium concentrations 2 to 3 times higher (between 1,426 and 3,060±87 pCi/L) than groundwaters obtained with a DDB (479 to 1,219±51 pCi/L). The longer the passive diffusion samplers were deployed, the higher the tritium concentration, suggesting that equilibrium with tritium was not reached within a 27-day deployment. DDB samples showed tritium results declining from the first to the last bailed sample for all three sampling events. This research suggests that tritium results from groundwater samples obtained from PDSs are more reproducible than samples obtained from DDBs. Also, PDSs likely do not accumulate isotopes of water but rather equilibrate with the ambient groundwater. On the other hand, both PDSs and DDBs were able to provide representative groundwater samples for major ions and have the potential to produce representative groundwater samples for stable isotopes of water (δD and δ{sup 18}O). Utilizing the most effective method for sampling groundwater is important for long-term groundwater monitoring. Quantifying the ability of passive diffusion groundwater samplers to obtain representative radioactive chemistry samples will help determine if these methods can be a useful alternative to pumping and purging wells for long-term monitoring. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

CPC data from TAMU TRACER campaign in the Houston, TX region from July to September 2022

During TRACER, the Texas A&M Rapid Onsite Atmospheric Measurements Van (ROAM-V) was deployed to capture airmasses behind (maritime) and ahead (continental) of the passage of the sea-breeze front through Houston. On select sampling days, ROAM-V sampled in the morning/mid-day on the coast and then transited to a second inland site for the afternoon/evening. The suite of instruments deployed on ROAM-V included a Condensation Particle Counter (CPC; GRIMM Model 5.403 CPC), Scanning Mobility Particle Sizer (SMPS; TSI 3750 detector, TSI 3082 classifier, TSI 3088 neutralizer, TSI 3081A Differential Mobility Analyzer), Cloud Condensation Nuclei counter (Droplet Measurement Technologies CCN Counter), micro pulse lidar (Droplet Measurement Technologies Micro Pulse LiDAR (miniMPL)), and a Davis Rotating Uniform size-cut Monitor (DRUM; DRUMAir 4-DRUM). Before sampling at each location, the latitude and longitude were recorded using the GPS on the phone application “My Altitude”. Onboard the ROAM-V, aerosol samples are drawn through a shared isokinetic inlet at a flow rate ranging from 3.5 to 7.0 LPM. A portion of this flow is directed through a cyclone impactor (Brechtel, Inc. Model SCC 0.732) and 0.3 LPM is directed to the CPC. To correct the data for particle losses, we used a two-step method. First, the measured SMPS size distributions were used to calculate total particle loss through the inlet for every SMPS scan at a single deployment location. Second, the CPC data was corrected for particle losses using the average of the total losses per scan. This calculation was done separately for each deployment location due to changes in the measured size distributions between locations. Particle loss from diffusion (based on Kesten, 1991 and Gormley, 1949), inertial impaction in 90-degree bends (based on Aerosol Measurement, 2011 and Crane, 1977), and cyclone impactor efficiency (based on Dirgo, 1985) were included in the loss calculation. When the SMPS was not sampling at a location (in the case of an instrument malfunction or operator error), the reported CPC data was corrected with an average of the total losses for the entire campaign at the specified deployment location (e.g., if we needed to correct Galveston data, then the average of all calculated losses at Galveston was taken). These flatline corrections were used for all data on 22/07/13, 22/07/20, 22/07/22, and the data from Galveston on 22/08/09. An additional flatline correction of 14% was applied to all CPC data based on laboratory calibration prior to and after the campaign. In laboratory calibrations, we identified that for the same sample of air the field CPC (GRIMM Model 5.403 CPC) undercounted the total concentration of particles by 14% compared to the ground truth laboratory CPC (TSI Model 3750). This data was collected for ARM Field Campaign AFC07055 and supported by DOE ASR grant DE-SC0021047. For any further questions, please feel free to contact the instrument PI, Sarah D. Brooks, sbrooks@tamu.edu. Kesten et. al. Calibration of a TSI Model 3025 Ultrafine Condensation Particle Counter. Aerosol Science and Technology, 15:2, 107-111, 1991. Gormley et. al. Diffusion from a Stream Flowing through a Cylindrical Tube. Proceedings of the Royal Irish Academy, Vol 52, 163-169, 1948. Aerosol Measurement: Principles, Techniques, and Applications, Third Edition. John Wiley & Sons, Inc, 2011. Crane et. al. Inertial Deposition of Particles in a Bent Pipe. Journal of Aerosol Science, Vol 8, 161-170, 1977. Dirgo et. al. Cyclone Collection Efficiency: Comparison of Experimental Results with Theoretical Predictions. Aerosol Science and Technology, 4:4, 401-415, 1985.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

CCN data from TAMU TRACER campaign in the Houston TX region from July to September 2022

During TRACER, the Texas A&M Rapid Onsite Atmospheric Measurements Van (ROAM-V) was deployed to capture airmasses behind (maritime) and ahead (continental) of the passage of the sea-breeze front through Houston. On select sampling days, ROAM-V sampled in the morning/mid-day on the coast and then transited to a second inland site for the afternoon/evening. The suite of instruments deployed on ROAM-V included a Condensation Particle Counter (CPC; GRIMM Model 5.403 CPC), Scanning Mobility Particle Sizer (SMPS; TSI 3750 detector, TSI 3082 classifier, TSI 3088 neutralizer, TSI 3081A Differential Mobility Analyzer), Cloud Condensation Nuclei counter (Droplet Measurement Technologies CCN Counter), micro pulse lidar (Droplet Measurement Technologies Micro Pulse LiDAR (miniMPL)), and a Davis Rotating Uniform size-cut Monitor (DRUM; DRUMAir 4-DRUM). Before sampling at each location, the latitude and longitude were recorded using the GPS on the phone application “My Altitude”. Onboard the ROAM-V, aerosol samples are drawn through a shared isokinetic inlet at a flow rate ranging from 3.5 to 7.0 LPM. A portion of this flow is directed through a cyclone impactor (Brechtel, Inc. Model SCC 0.732) and 0.5 LPM is directed to the CCN. To calculate particle losses, we used a two-step method. First, the measured SMPS size distributions were used to calculate particle loss through the inlet during sampling. Second, the corrected SMPS data was used to calculate the average of the total losses per scan down the CCN line. Then, the correction was applied to the CCN data. This calculation was done separately for each deployment location due to changes in the measured size distributions between locations. Particle loss from diffusion (based on Kesten, 1991 and Gormley, 1949), inertial impaction in 90-degree bends (based on Aerosol Measurement, 2011 and Crane, 1977), and cyclone impactor efficiency (based on Dirgo, 1985) were included in the loss calculation. When the SMPS was not sampling at a location (in the case of an instrument malfunction or operator error), the reported CPC data was corrected with an average of the total losses for the entire campaign at the specified deployment location (e.g., if we needed to correct Galveston data, then the average of all calculated losses at Galveston was taken). These flatline corrections were used for all data on 22/07/13, 22/07/20, 22/07/22, and the data from Galveston on 22/08/09. The supersaturation uncertainty is estimated conservatively at +/- 0.03%, where variation in the inlet temperature, pressure, and calibration technique prevents a more accurate measurement. Confidence in the reported supersaturation measurements is based on a pre-campaign calibration (following the methods from our previous work and Deng, 2014 based on Rose, 2008) in addition to inter-comparisons with the DOE for two days (22/08/18 and 22/09/01) where TAMU was co-located with AMF1. The inter-comparisons show good agreement between our instrument and the DOEs instrument on both days at all supersaturations. After the last inter-comparison on 22/09/01, there was no indication of a malfunction by our instrument through the rest of the campaign. Unfortunately, the instrument was dropped during demobilization. A post-campaign calibration was conducted, which showed a substantial departure from the pre-campaign calibration. The drop may have damaged the instrument’s ability to produce the desired supersaturations. Therefore, we do consider the data after 22/09/01 to be correct, but it should be used with caution. The CCN counter sampled for 3 minutes at each supersaturation setpoint (0.2, 0.4, 0.6, 0.8, 1.0, and 1.2%). At the end of a cycle, the instrument was set to 0.01% supersaturation for 5 minutes. The data is comprised of the last 60 seconds of each supersaturation set point (0.2, 0.4, 0.6, 0.8, 1.0, and 1.2%) to ensure the instrument stabilized and was able to reach thermal equilibrium. We removed the data during the periods where there were operational difficulties, setup, or maintenance. This data was collected for ARM Field Campaign AFC07055 and supported by DOE ASR grant DE-SC0021047. For any further questions, please feel free to contact the instrument PI, Sarah D. Brooks, sbrooks@tamu.edu. Rose et. al. Calibration and Measurement Uncertainties of a Continuous-Flow Cloud Condensation Nuclei Counter (DMT-CCNC): CCN Activation of Ammonium Sulfate and Sodium Chloride Aerosol Particles in Theory and Experiment. Atmos. Chem. Phys., 8, 1153-1179, 2008. Deng et. al. Using Raman Microspectroscopy to Determine Chemical Composition and Mixing State of Airborne Marine Aerosols over the Pacific Ocean. Aerosol Science and Technology, Vol 48, Issue 2, 2014. Kesten et. al. Calibration of a TSI Model 3025 Ultrafine Condensation Particle Counter. Aerosol Science and Technology, 15:2, 107-111, 1991. Gormley et. al. Diffusion from a Stream Flowing through a Cylindrical Tube. Proceedings of the Royal Irish Academy, Vol 52, 163-169, 1948. Aerosol Measurement: Principles, Techniques, and Applications, Third Edition. John Wiley & Sons, Inc, 2011. Crane et. al. Inertial Deposition of Particles in a Bent Pipe. Journal of Aerosol Science, Vol 8, 161-170, 1977. Dirgo et. al. Cyclone Collection Efficiency: Comparison of Experimental Results with Theoretical Predictions. Aerosol Science and Technology, 4:4, 401-415, 1985.

54 ENVIRONMENTAL SCIENCES↗