Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analyzed”

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,765 records · Page 98

Phase II and III Solid Waste Management Unit Assessment and Confirmatory Sampling Report for Per- and Polyfluorinated Alkyl Substances, Potential Release Location 237, Kennedy Space Center, Florida

Per- and polyfluoroalkyl substances (PFAS) Phase II and III Solid Waste Management Unit (SWMU) Assessment (SA) and Confirmatory Sampling (CS) were conducted at the John F. Kennedy Space Center (KSC), Potential Release Location (PRL) 237, located in Merritt Island, Brevard County, Florida for the National Aeronautics and Space Administration (NASA) under Contract 80KSC019D0010/80KSC019F0289. A Phase I SA and CS Report was previously submitted that identified 33 locations of concern (LOCs) at KSC associated with the storage, use, or release of materials containing PFAS. The objectives of the Phase II and III SA and CS were to: - Complete site and source characterization of identified LOCs; - Identify areas of potential concern (AOPCs) where past or present activities may have resulted in a release of PFAS to the environment; - Qualitatively characterize migration potential of released PFAS to environmental media (primarily soil, groundwater, surface water, and sediment); and - Identify potential receptors for exposure to PFAS-impacted environmental media. Surface water samples were collected from LOCs/AOPCs and locations on the KSC periphery during PFAS CS Phases I, II and III. Including duplicate samples, 143 environmental samples were analyzed for 14 PFAS compounds and 104 environmental samples were analyzed for 18 PFAS compounds. Ten of the analytes on the 14 compound list were detected at least once in surface water at KSC, with a wide range in magnitude of detections. The State of Florida has surface water screening levels (SWSLs) for perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS) in surface water, both for human and ecological receptors. The human health screening levels are lower than the ecological ones and are used in this report. PFOA was not detected in surface water at concentrations exceeding the SWSL at any LOC or Area of Potential Concern (AOPC). PFOS was detected in surface water at concentrations exceeding the SWSL at 28 LOCs/AOPCs, plus 42 locations on the KSC periphery. The data suggest a wide distribution of PFAS, particularly PFOS, in surface water but the sources of surface water contamination are not yet adequately understood. Sediment samples were collected during Phases II and III but not Phase I. Including duplicate samples, 104 environmental samples were analyzed for 14 PFAS compounds and 10 environmental samples were analyzed for 18 PFAS compounds. Twelve compounds were detected in sediment at KSC, with a wide range in magnitude of detections. There are currently no screening levels or cleanup target levels (CTLs) for sediment. The need to delineate and further characterize PFAS in sediment will be evaluated for those sites that move forward into Site Assessment. Groundwater samples were collected via DPT from LOCs/AOPCs and locations on the KSC periphery during PFAS CS Phases I, II and III. Including duplicate samples, 812 environmental DPT groundwater samples were analyzed for 14 PFAS compounds and 247 environmental samples were analyzed for 18 PFAS compounds. Twelve of the 14 compounds were detected in groundwater at KSC, with a wide range in magnitude of detections. The State of Florida has provisional groundwater cleanup target levels (GCTLs) for PFOA, PFOS and the additive concentration of PFOA+PFOS. PFOA was detected in groundwater at concentrations exceeding its provisional GCTL at 16 LOCs/AOPCs, with no exceedances on the KSC periphery. PFOS was detected in groundwater at concentrations exceeding its provisional GCTL at 30 LOCs/AOPCs, plus two locations on the KSC periphery. PFOA+PFOS was detected in groundwater at concentrations exceeding the provisional GCTL at 31 LOCs/AOPCs, plus two locations on the KSC periphery. Soil samples were collected during Phases II and III but not Phase I. Including duplicate samples, 292 environmental samples were analyzed for 14 PFAS compounds and 80 environmental samples were analyzed for 18 PFAS compounds. Fourteen compounds were detected in soil at KSC, with a wide range in magnitude of detections. The State of Florida has provisional soil cleanup target levels (SCTLs) for PFOA and PFOS. There were no PFOA results above its provisional SCTL. There was one LOC with a result above the provisional SCTL for residential settings. Recommendations are presented for each LOC and AOPC.

PFAS↗

GeoThermalCloud: Cloud Fusion of Big Data and Multi-Physics Models using Machine Learning for Discovery, Exploration, and Development of Hidden Geothermal Resources

The primary goals of this project are exploring hidden geothermal resources in the U.S.A. and designing profitable enhanced geothermal systems (EGS). Many processes and parameters control geothermal exploration and energy production from geothermal fields. Diverse datasets (e.g., geology, geochemistry, geophysics, satellite, airborne geophysics) are available to help characterize subsurface geothermal conditions. Sparse and multi-scale characteristics of these datasets prohibit properly leveraging these datasets for geothermal exploration and profitable EGS design. Recent advancements in machine learning (ML) promise to resolve these issues. The tremendous challenges and risks of geothermal exploration and production bring the demand for novel ML methods and tools that can (1) analyze large field datasets, (2) assimilate model simulations (large inputs and outputs), (3) process sparse datasets, (4) perform transfer learning (between sites with different exploratory levels), (5) extract hidden geothermal signatures in the field and simulation data, (6) label geothermal resources and processes, (7) identify high-value data acquisition targets, and (8) guide geothermal exploration and production by selecting optimal exploration, production, and drilling strategies. To address these necessities, ML-based geothermal resources exploration and enhanced geothermal systems (EGS) design tools have been developed. The exploration tool is called GeoThermalCloud and EGS design tool is called GeoDT-ML. GeoThermalCloud (https://github.com/SmartTensors/GeoThermalCloud.jl) utilizes a LANL unsupervised ML platform called SmartTensors (https://tensors.lanl.gov/) to automate data analyses and interpretations by extracting hidden signatures to identify geothermal prospects. Also, it enables the identification of critical measurements needed to identify geothermal resource signatures. Alternatively, GeoDT-ML (https://github.com/SmartTensors/GeoThermalCloud.jl/tree/master/EGS) is an ML-based alternative to GeoDT (https://github.com/GeoDesignTool/GeoDT.git), a fast, simplified multi-physics solver to evaluate EGS project designs in uncertain geologic systems. GeoDT-ML leverages recent advances in deep learning and high-performance computing. It is a faster and simpler version of GeoDT. To make this project a success, we used capabilities of LANL, PNNL, Google, Stanford, and Julia Computing. We analyzed eight datasets of the U.S.A. using GeothermalCloud and demonstrated potential highly prospective geothermal resources and identified key factors defining highly prospective sites. The first data set includes 44 locations in southwest New Mexico and 18 geological, hydrogeological, geophysical, geothermal, geochemical attributes. We defined low- and medium-temperature hydrothermal systems and discovered a new highly prospective site. The second data set analyzed 18 shallow water chemistry attributes at 14,342 locations in the Great Basin. It demarcated modestly, moderately, and highly prospective sites including key attributes for each type of prospectivity. The third data set analyzed Utah FORGE data including satellite (InSAR), geophysical (gravity, seismic), geochemical, and geothermal attributes. Here, we performed prospectivity analysis to identify future drilling locations using geological, geochemical, and geophysical attributes. Maps of temperature at depth and heat flow are constructed based on the available data. Prospectivity maps were generated, and drilling locations were proposed for future geothermal field exploration. The fourth data set analyzed 21 attributes at 120 locations in Tularosa Basin, New Mexico; data comes from past play fairway analyses in this region. ML analyses identified geothermal signatures associated with modestly, moderately, and highly hydrothermal systems. We also defined dominant attributes and spatial distribution of the geothermal signatures. The fifth, sixth, seventh, and eighth datasets include Tohatchi Springs, New Mexico, Hawaii, Brady site, Nevada, and EGS Collab, respectively. Moreover, we coupled GeothermalCloud and magnetotellurics data to pinpoint drilling locations for developing geothermal projects in the Tularosa Basin, New Mexico. GeothermalCloud found potential prospective locations for geothermal resources near White Sands Missile Range and McGregor Range at Fort Bliss. Magnetotellurics data determined the potential depth (~1800m) of geothermal prospects at McGregor Range based on apparent resistivity structures/layers in the subsurface. The McGregor Range consists of three resistivity layers and two resistivity structures. Magnetotellurics data also helps identify that the western portion of the McGregor Range has thick and low-resistivity earth materials. The low resistivity to the west is most likely for a fault system. Assuming temperature is consistent with a geothermal reservoir, the west-central part of the McGregor Range has the highest geothermal potential because of the increase in porosity and associated permeability attributed to the interpreted fault system. Also, we devised a coupling strategy between a process model and GeothermalCloud to characterize hydrogeological conditions and geothermal conditions, respectively. The process model characterizes hydrogeological and geothermal conditions on highly prospective geothermal sites provided by GeothermalCloud. We developed a physics-informed neural network (PINN) version of the Burns equation that can be easily coupled with GeothermalCloud. Furthermore, we performed an optimal design decision maximizing the economic value of an EGS power plant. This study optimized the range of well spacing between injection and production wells maximizing net present value in dollars (NPV). For this task, we used the GeoDT to simulate the Utah FORGE EGS development cycle from the initial well design to the end of production. Next, we accomplished another crucial task, which is predicting permeability of geothermal reservoirs. Predicting permeability of geothermal reservoirs is a non-trivial task because of huge computational runtime of simulation and lack of measurements. To avoid these limitations, we used easy-to-measure chemical concentrations in the subsurface as measurement data and convolutional neural network based ML model of a high-fidelity model. Next, we predicted permeability using Markov chain Monte Carlo simulation. We found that Markov chain Monte Carlo simulation predicts permeability with a high certainty if the prediction zone in the simulation area has chemical concentration data. Finally, we analyzed the DOE funded INGENIOUS and GeoDAWN projects data. For discovering hidden geothermal systems in the Great Basin, the INGENIOUS project accumulated old data, collected new data, and released them in 2022. The dataset includes a total of 24 geological, geophysical, and geochemical attributes. Data resolution and scale significantly vary prohibiting an appropriate usage. To avoid such limitations, we brought all data in the same resolution and scale by applying the inverse distance weighting interpolation technique for predicting data in unsampled locations. Subsequently, we analyzed LiDAR data of the GeoDAWN project. We received data in tiles format. The DOE’s overarching goal is to use ML on LiDAR data for finding favorable geological structures (e.g., step up faults in Brady, Nevada). To serve the purpose, we need to label favorable geologic structures that correspond to LiDAR data. We wrote an algorithm to label the LiDAR data with the favorable geologic structures.

15 GEOTHERMAL ENERGY↗

Ice nucleation measurements from DRUM impactors during TRACER campaign in the Houston TX region from July to September 2022

During TRACER, three Davis Rotating Uniform size-cut Monitors (DRUM; DRUMAir 4-DRUM) were used to collect aerosols for ice nucleation measurements in the Brooks laboratory at Texas A&M University. The three instruments were located at AMF1 in La Porte, Ancillary site in Guy (ANC), and onboard the Texas A&M Rapid Onsite Atmospheric Measurements Van (ROAM-V). 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 one of the Davis Rotating Uniform size-cut Monitors (DRUM; DRUMAir 4-DRUM). Before sampling at each location, the latitude and longitude were recorded using the GPS on the phone application “My Altitude”. Each DRUM sampler was operated at a flow rate of 23 LPM. Each DRUM has four stages where the aerodynamic diameter size cuts are as follows: stage 1 larger than 3 μm, stage 2 from 3 to 1.2 μm, stage 3 from 1.2 to 0.34 μm, and stage 4 from 0.34 to 0.15 μm. Pretreated aluminum foil was used as a substrate on all stages. At AMF1 and ANC sites, the DRUMs were operated on the shared aerosol inlet, with generous support of the DOE ARM site staff. These instruments rotated 24 mm per day and only contained aluminum foil substrates. Substrates were changed weekly and transported to Texas A&M for storage in -80C freezer until analysis. The DRUM onboard ROAM-V was operated at a faster rotation rate of 150 mm per day to clearly separate the multiple deployment locations for ROAM-V. For the ROAM-V DRUM, substrates were changed every deployment and transported to Texas A&M for storage in -80C freezer until analysis. At Texas A&M, ice nucleation experiments were conducted to measure the ice nucleation temperature of the ambient aerosol samples collected from the three DRUMs using our previously established procedures (Alsante et al., 2023; Fornea et al., 2009; Matthews et al., 2023). For ice nucleation, only samples collected on stage 3 were analyzed, given that these are the most relevant size (1.2 to 0.34 μm diameter) for potential ice nucleating particles. For the AMF1 and ANC sites, we cut and analyzed 2 mm (2-hour) samples. We analyzed the time periods of the AMF1 site instrument when the ROAM-V was deployed. A 72-hour period from July 11th at 23:49 through July 15th at 1:49 was analyzed from the ANC site instrument. For the ROAM-V instrument, we separated the daily samples by site location. Between 1- and 6-hour independent samples were analyzed at each location. We also cut the ROAM-V samples in half to allow for compositional analysis of the aerosol on the other half of the substrate. All viable samples from the ROAM-V were analyzed. Analysis was done using a custom-built ice nucleation apparatus, recently updated to include an array of 16 individual samples (Matthews et al., 2023). On our experimental setup, we used 100 μL of Ultra-High-Performance Liquid Chromatography (UHPLC) water (Sigma Aldrich, >99.9% purity) to wash off the aerosol from the DRUM substrate. Then, we micropipetted 2 μL droplet samples into each of the 16 wells of the array, using hydrophobically coated microscope slides. Experiments would cycle 28 times from 10 C to -40C over a 20-hour period. During analysis of the TRACER samples, nine experiments were conducted with UHPLC water process blanks that followed an identical preparation procedure. The average freezing temperature and standard deviation for these process blanks is -27.7±2.4 C. This data was collected for ARM Field Campaign AFC07023 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. Alsante, A. N., Thornton, D. C., & Brooks, S. D. (2023). Ice nucleation catalyzed by the photosynthesis enzyme RuBisCO and other abundant biomolecules. Communications Earth & Environment, 4(1), 51. Fornea, A. P., Brooks, S. D., Dooley, J. B., & Saha, A. (2009). Heterogeneous freezing of ice on atmospheric aerosols containing ash, soot, and soil. Journal of Geophysical Research: Atmospheres, 114(D13). DOI:10.1029/2009JD011958 Matthews, B. H., Alsante, A. N., & Brooks, S. D. (2023). Pollen Emissions of Subpollen Particles and Ice Nucleating Particles. ACS Earth and Space Chemistry.

54 ENVIRONMENTAL SCIENCES↗

The comet rendezvous asteroid flyby mission

The Comet Rendezvous Asteroid Flyby Mission (CRAF) is described. After gravity assists from Venus and Earth, the spacecraft will fly by the asteroid 46 Hestia en route to a rendezvous with P/Tempel 2 in Nov. 1996, when the comet is near aphelion. The scientific experiments for the CRAF mission are: an imaging system; a visual and infrared mapping spectrometer; an infrared radiometer; a penetrator carrying a gamma-ray spectrometer, a scanning differential calorimeter and evolved gas analyzer, accelerometers, and temperature probes; a neutral gas and thermal ion mass spectrometer; a secondary ion mass spectrometer for analyzing dust, gas, and thermal ions; a scanning electron microscope and particle analyzer; an X-ray fluorescence and gas chromatographic analyzer for collected dust and ice samples; a dust counter and velocity analyzer; a retarding potential ion mass spectrometer; a supra-thermal ion mass spectrometer and electron analyzer; a magnetometer; and a coordinated radio, electron, and plasma wave analyzer.

Neugebauer, M.↗

rHealth One Demonstration Aboard ISS: A Microfluidic Bioanalyzer Based on Sheath/Hydrodynamic Focusing Flow Cytometry

The Exploration Medical Capability (ExMC) element aims to provide astronauts with the means for their own health monitoring, diagnosis, and treatment during exploration missions. As space flight ventures further from earth, the need for autonomous medical care increases under greater constraints on size, mass, and resources. One pillar for diagnosis that crew would be separated from is laboratory analysis. Even now on ISS, blood samples must be collected and returned to earth for testing. In response, ExMC is assessing how assays for hematology, bone health, radiation exposure etc. could be addressed through commercial-off-the-shelf (COTS) and Small Business Innovation Research (SBIR) funded bioanalyzers that are miniaturizing lab technology. Since missions will reach distances where there are no timely replacements, validation on the International Space Station (ISS) is a necessary part of that assessment. In partnership with NASA Johnson Space Center (JSC) Immunology Lab and the Research Operations and Integration (ROI) element, ExMC conducted a technology demonstration on ISS of the rHEALTH ONE, a flow cytometry based bioanalyzer, to assess future devices based on this design. In flow cytometry, there are predominantly three ways to focus the cells (or particles) into a single file stream: hydrodynamic focusing, microcapillary, and acoustic focusing. The rHEALTH ONE utilizes the commercial standard, sheath-based hydrodynamic focusing. rHEALTH itself represents both a company and a suite of medical tools NASA has funded through SBIR grants towards the development of a promising diagnosis instrument for exploration missions. rHEALTH ONE is the interim version of the technology, functional as a benchtop analyzer and test bed for the next generation of rHEALTH in development. Several modifications were made to the rHEALTH ONE analyzer for operation in microgravity. For function, fluid management was key. A sheath-based analyzer uses sheath fluid to flow the sample, cleaning fluid to prevent biological contamination, and a reservoir to collect the liquid waste. The rHEALTH ONE analyzer uses bottles dependent on gravity to separate the air and liquid pathways and keep the liquids contained. It uses only 1 psig of air pressure to directly push the liquid through the device, requiring little-to-no resistance at inlet and outlet. A microtubing assembly with self-sealing luer connectors was designed – featuring 0.014 mm thick durable medical balloons to hold water inside the supply bottles – to create safe containment and easy access for the crew while maintaining the analyzer’s function. For safety, copper tape was added to the interior of the plastic housing to reduce electromagnetic interference, fluid and electrical connections were secured against vibration and leaks, and gaps were further sealed to ensure containment of the optical block and lasers. Water as the sheath and cleaning fluid and TOX 0 samples were used to reduce the biohazard risk to the crew. In May 2022, European Space Agency astronaut Samantha Cristoforetti demonstrated the rHEALTH ONE aboard ISS for its sample loading, flow cytometry, and data collection capability in microgravity. The JSC Immunology Lab provided flow cytometry expertise, designed the sample test protocol, and manufactured and benchmarked the flight samples on a ‘gold-standard’ flow cytometer. The analyzer was primed with water, purged of air, and four calibration solutions of polystyrene microparticles were tested to characterize its performance. Tightly controlled procedures and excellent execution prevented air bubble interference during air/water separation, fluid transfer, sample mixing and loading. Data showed a slight increase in signal noise on 3 of 5 channels and an anomaly of fluorophores migrating in one sample (confirmed by JSC post-flight). Flight results correctly detected the change in sample, matched the analyzer’s ground performance, and were consistent with the gold-standard’s ground results. Although areas were noted for improvement, these outcomes signified complete mission success.

R. S. Miller↗

A Reversed Impedance-Based Stability Criterion for IBR Grids

The existing impedance-based stability criterion is effective for analyzing local control interactions; however, it is difficult to scale the existing criterion to analyze wide-area control interactions among numerous IBRs through a complex power system network. The scaled version of the existing criterion requires the impedance response of each IBR in the system as well as of the network looking from all the IBRs. It is quite challenging to obtain all these impedance responses because of the computational effort and the requirement of separately scanning the impedance of the network and the IBRs. We propose a reversed criterion for the impedance-based stability analysis to address these problems. In contrast to the existing criterion, the reversed criterion analyzes the stability of a power system when an IBR is disconnected from the system. The reversed criterion estimates the impact of an IBR on the frequency and damping of power system oscillation modes using the impedance scans of only the IBR and the grid at its terminal. It can be sequentially applied at different IBRs to evaluate their impact on the power system stability. In addition to scalability, the reversed criterion gives flexibility to focus only on a few selected IBRs, depending on their rating, the magnitude of oscillations observed at their terminals, and the vendor support available for implementing stabilizing control system updates. The reversed criterion is demonstrated on a 14-bus power system with 100% IBRs.

control interactions↗

Graded Buffer Bragg Reflectors with High Reflectivity and Transparency for Metamorphic Optoelectronics

A graded buffer Bragg reflector (GBBR) is a bifunctional device component that provides the reflection of a distributed Bragg reflector and the adjustable lattice constant of a compositionally graded buffer. Prior work showed that these properties can be combined with low threading dislocation density and high reflectivity. Here, we design and demonstrate complex GBBRs for specific metamorphic solar cell applications. We design buffers that provide (1) high reflectivity over a narrow bandwidth, for quantum well solar cells, (2) reflectivity over a wide bandwidth, for optically thin solar cells, and (3) low sidelobe reflection, for multijunction devices that demand low out-of-band reflection. Apodized and triple GBBRs are demonstrated, and transparency is always considered, requiring designs with carefully engineered material combinations. A GBBR with a reflection of 99% is demonstrated, as well as a triple GBBR that has over 80% reflection for 100?nm of the spectrum around 800?nm. We also analyze potential deviations in a baseline GBBR from that of a perfect DBR using transmission electron microscopy to analyze imperfections in the material and modeling to analyze the impact of imperfect refractive index data for lattice-mismatched AlGaInAs. Minimal crosshatch roughness and unintentional thickness variation occur throughout the buffer, which likely influences reflection slightly. Small deviations between the calculated AlGaInAs and utilized AlGaAs refractive index exist, giving methods for future GBBR improvement. While the GBBR designs are intended for metamorphic solar cells, the broad and high reflection may also be useful for other optoelectronic devices such as light-emitting diodes or lasers.

41 EE - Solar Energy Technologies Office (EE-4S)↗

SAIL-Isotopologue (SAIL-iso) Field Campaign Report

The SAIL-iso project took place from November 19, 2022 to June 3, 2023 as part of the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s Surface Atmosphere Integrated Field Laboratory (SAIL) project at the main project site in Gothic, Colorado. The project sought to use measurements of the stable isotopic composition of water vapor to better understand cloud-aerosol interactions in a remote mountainous watershed. Water vapor isotopologue and humidity measurements were determined using a Los Gatos Research Triple Water Vapor Isotope Analyzer (TWVIA). The instrument was deployed in the guest van and consists of three main components: a water vapor isotope analyzer, 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. Multiple sources of uncertainty are 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) Vienna Standard Mean Ocean Water – Standard Light Antarctic Precipitation (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 δ 18 O.

54 ENVIRONMENTAL SCIENCES↗

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↗

Neural Network and Regression Approximations in High Speed Civil Transport Aircraft Design Optimization

Nonlinear mathematical-programming-based design optimization can be an elegant method. However, the calculations required to generate the merit function, constraints, and their gradients, which are frequently required, can make the process computational intensive. The computational burden can be greatly reduced by using approximating analyzers derived from an original analyzer utilizing neural networks and linear regression methods. The experience gained from using both of these approximation methods in the design optimization of a high speed civil transport aircraft is the subject of this paper. The Langley Research Center's Flight Optimization System was selected for the aircraft analysis. This software was exercised to generate a set of training data with which a neural network and a regression method were trained, thereby producing the two approximating analyzers. The derived analyzers were coupled to the Lewis Research Center's CometBoards test bed to provide the optimization capability. With the combined software, both approximation methods were examined for use in aircraft design optimization, and both performed satisfactorily. The CPU time for solution of the problem, which had been measured in hours, was reduced to minutes with the neural network approximation and to seconds with the regression method. Instability encountered in the aircraft analysis software at certain design points was also eliminated. On the other hand, there were costs and difficulties associated with training the approximating analyzers. The CPU time required to generate the input-output pairs and to train the approximating analyzers was seven times that required for solution of the problem.

Patniak, Surya N.↗

Mass Spectrometer Containing Multiple Fixed Collectors

A miniature mass spectrometer that incorporates features not typically found in prior mass spectrometers is undergoing development. This mass spectrometer is designed to simultaneously measure the relative concentrations of five gases (H2, He, N2, O2, and Ar) in air, over the relative-concentration range from 10(exp -6) to 1, during a sampling time as short as 1 second. It is intended to serve as a prototype of a product line of easy-to-use, portable, lightweight, highspeed, relatively inexpensive instruments for measuring concentrations of multiple chemical species in such diverse applications as detecting explosive or toxic chemicals in air, monitoring and controlling industrial processes, measuring concentrations of deliberately introduced isotopes in medical and biological investigations, and general environmental monitoring. The heart of this mass spectrometer is an integral combination of a circular cycloidal mass analyzer, multiple fixed ion collectors, and two mass-selective ion sources. By circular cycloidal mass analyzer is meant an analyzer that includes (1) two concentric circular cylindrical electrodes for applying a radial electric field and (2) a magnet arranged to impose a magnetic flux aligned predominantly along the cylindrical axis, so that ions, once accelerated into the annulus between the electrodes, move along circular cycloidal trajectories. As in other mass analyzers, trajectory of each ion is determined by its mass-to-charge ratio, and so ions of different species can be collected simultaneously by collectors (Faraday cups) at different locations intersected by the corresponding trajectories (see figure). Unlike in other mass analyzers, the installation of additional collectors to detect additional species does not necessitate increasing the overall size of the analyzer assembly.

Moskala, Robert↗

Fast Plasma Instrument for MMS: Simulation Results

Magnetospheric Multiscale (MMS) mission will study small-scale reconnection structures and their rapid motions from closely spaced platforms using instruments capable of high angular, energy, and time resolution measurements. The Dual Electron Spectrometer (DES) of the Fast Plasma Instrument (FPI) for MMS meets these demanding requirements by acquiring the electron velocity distribution functions (VDFs) for the full sky with high-resolution angular measurements every 30 ms. This will provide unprecedented access to electron scale dynamics within the reconnection diffusion region. The DES consists of eight half-top-hat energy analyzers. Each analyzer has a 6 deg. x 11.25 deg. Full-sky coverage is achieved by electrostatically stepping the FOV of each of the eight sensors through four discrete deflection look directions. Data compression and burst memory management will provide approximately 30 minutes of high time resolution data during each orbit of the four MMS spacecraft. Each spacecraft will intelligently downlink the data sequences that contain the greatest amount of temporal structure. Here we present the results of a simulation of the DES analyzer measurements, data compression and decompression, as well as ground-based analysis using as a seed re-processed Cluster/PEACE electron measurements. The Cluster/PEACE electron measurements have been reprocessed through virtual DES analyzers with their proper geometrical, energy, and timing scale factors and re-mapped via interpolation to the DES angular and energy phase-space sampling measurements. The results of the simulated DES measurements are analyzed and the full moments of the simulated VDFs are compared with those obtained from the Cluster/PEACE spectrometer using a standard quadrature moment, a newly implemented spectral spherical harmonic method, and a singular value decomposition method. Our preliminary moment calculations show a remarkable agreement within the uncertainties of the measurements, with the results obtained by the Cluster/PEACE electron spectrometers. The data analyzed was selected because it represented a potential reconnection event as currently published.

Figueroa-Vinas, Adolfo↗

X-Ray Photoelectron Spectroscopy

The interaction of photon and the electron goes back to the early part of 19th century emanating from the photo-electric effect depicted by none other than Albert Einstein (Ref 1) described in 1905, and the redistribution of kinetic energy resulting from the interaction of x-ray and solids reported during early part of the century (Ref.2). The spectrum resolutions obtained at that time was not sufficient to observe distinct peaks in spectra for materials. Thus, these phenomena hardly attracted any attention for many years following these discoveries. The modern X-ray Photoelectron Spectroscopy (XPS) has been possible by the extensive and significant contribution from Kai Siegbahn and others (Ref.3, 4) of Uppsala University. Siegbahn developed and employed a high-resolution electron spectrometer that revealed electron peaks in a spectrum emerging from the interaction of x-rays and solids. Eventually, Kai Siegbahn received Nobel Prize in 1981 for his contributions to XPS. Around 1958, shifts in elemental peaks were realized in compounds when the same elements are bound to other but different elements. This discovery resulted in the chemical state identification in various chemicals as well as the oxidation states of atoms in compounds. Because of these useful physical effects, the Uppsala group named XPS with a synonymous name of ESCA (Electron Spectroscopy for Chemical Analysis) used widely today and will be used here alternatively. Therefore, XPS or ESCA not only identifies the element, but also the compound these elements form, from their chemical shifts. Compared to other micro-analytical techniques such as Energy Dispersive (EDS) or Wavelength Dispersive (WDS) techniques, XPS analyzes only few atomic layers present on the surface. This was discovered early in 1966 (Ref. 5). While this has awarded a merit to the analytical technique to analyze very thin layers such as films and coatings, it often analyzes the adsorbed superficial gases and contaminations on a sample introduced to its analytical chamber. This necessitates the surface is cleaned and the underlying material, material of interest, is exposed in a clean environment such that the material of interest is analyzed. The cleaning is accomplished by a scanning ion gun within the analytical chamber of the instrument. Ion gun uses an argon gas and is commonly attached in most modern machines. Reliable and efficient vacuum systems employed in modern machines does not allow adsorbed layers to rebuild after the surface is cleaned. Development of efficient and reliable vacuum pumps over these developmental years is yet another important step in the commercialization of XPS machines. Vacuum levels of better than 10-7 torr are essential to increase the mean free path of electrons released from the sample surface. Thus, modern machines are equipped with high capacity ion, turbo or cryogenic pumps in their analytical chambers. Today, XPS has advanced from an applied physics laboratory to industry for use in quality control as well as analysis of contaminants and has taken a dominant role in microanalysis. Its uniqueness arises from the fact that it is considered non-destructive compared to other common micro-analytical techniques using the electron and ion excitation sources. Polymers and plastics could be analyzed since the binding energies of saturated and unsaturated bonds in atoms could be separated. Extremely thin layers could be analyzed including materials with layered structures. The technique, though did not advance for many years, has now opened a new window for research as well as applications in industry due to its ability to separate and measure the chemical shifts in bound elements. Principles

Panda, Binayak↗

Olivine Dissolution and Formation of Secondary phases in Ultramafic Soils

Introduction: Olivine has been proposed as an indicator for the duration of water-rock interaction within Martian rocks and sediments [1-3]. The use of olivine as a mineralogical indicator for past aqueous alteration on Mars requires interpretation of a complex combination of factors including pH, temperature, and composition [5,6]. Here, we examine the persistence of natural olivine within terrestrial ultramafic soils (Fe/Mg-rich, Al-poor) developing under different climatic conditions and the incipient dissolution of emplaced forsterite (Fo) and fayalite (Fa) surfaces to investigate environmental effects on incipient olivine dissolution, olivine persistence in soils, and formation of secondary phases. Methods: Field Sites. We examined olivine weathering and secondary material formation in ultramafic soils at 6 sites in the Klamath Mountains (KM) of northern California with a mean annual temperature of ~12.8℃ and precipitation of ~55.7-95.4 cm/year [7], and soil pH of ~6.5-7.3; 4 sites in the Tablelands (TB) of Newfoundland, Canada with a mean annual temperature of <3.9℃ and precipitation of ~120.0 cm/year [8], and soil pH of ~7.7; and at 3 sites at Pickhandle Gulch (PG), Nevada with a mean annual temperature of ~14.4℃ and precipitation of ~14.1 cm/year [7], and soil pH of ~8.5. Sampling sites span an age range of ~12.1-50+ kya in the Klamath Mountains [9,10] and ~13-30 kya in the Tablelands [11]. Pickhandle Gulch sites are undated. Parent Material and Soil Analyses. Polished thin sections of bulk soil prepared by Wagner Petrographic, Inc were carbon-coated and analyzed on a JEOL 2100 SEM in back-scattered electron mode in the EMIL lab at UNLV and at the 13-ID-E synchrotron beamline at Argonne National Laboratory using µXRF, µXRD, and XAS. Soil and parent material samples were powdered in a Fritsch pulverisette and analyzed by XRD and soil by VNIR. Soil preparation is further described in [12]. Disk Preparation, Burial, and Collection. Fo disks were cut from a column prepared via hot-pressing and Fa disks by sintering synthetic fayalite powder, see [13] for detail. Disks were polished to a 0.25-micron level with diamond grit. Disks were buried in 3 KM soils, 4 TB soils, and 3 PG soils, collected after exactly 365 days, and washed gently with 100% reagent grade ethanol to remove potential adhered soil material. Weathered disks and soil samples were stored in a -20℃ freezer until analysis. Unaltered control disks prepared identically to the buried disks were stored at -20℃ for the duration of the experiment. Disk Analyses. One Fo and Fa disk from each climate zone was analyzed on a variable pressure Zeiss Supra 40VP SEM at Northern Arizona University. A separate Fo and Fa disk from each climate zone was analyzed by XPS using a Physical Electronics VersaProbe II at the Penn State Univ. Materials Characterization Lab after a Na-dodecyl sulfate wash and ozonation to remove carbon contamination as in [14]. VNIR measurements were conducted at Johnson Space Center using an ASD FieldSpec3 under ambient lab conditions on a separate Fo and Fa disk from each climate zone. One separate Fo and Fa control sample was analyzed for each technique for comparison with weathered samples. XPS uncertainty was determined from 5 repeat measurements on controls. Results: Bedrock and Soil Results Olivine is present in the parent material in the KM and TB. Olivine is found in ~12.1 ka KM soils but is absent from all older soils, while persisting into the oldest (>20 ka) TB soil (Figure 1). In both locations, olivine is found as cores surrounded by a serpentine rind (Figure 2). VNIR spectra from the analyzed soils possess strong OH-associated spectral features at ~2.33 µm indicating the presence of Mg-rich phyllosilicates as well as ferric-oxide features at ~0.92 µm in the KM (Figure 3). Primary crystalline silicate grains mostly incorporate Fe2+, while poorly crystalline weathering rinds are best fit by ferric oxide XAS standards (Figure 4). µXRF also shows that Fe and Ni concentrate in weathering rinds and Cr remains within interior silicate grains (Figure 4). Buried Sample Results All Fo surfaces exhibited formation of dissolution features including shallow pitting not observed on controls. Dissolution features were most visually widespread on the KM disk (Figure 5). Leaching of Mg from KM and TB Fo disks was evident from <1.6 Mg/Si ratios measured by XPS (Figure 6). Fe-rich precipitates in SEM (Figure 5) and Fe presence in XPS scans (Figure 6) indicate Fe deposition onto KM and TB Fo disk surfaces. The appearance of a spectral feature at 0.55 µm in the VNIR spectra from the TB Fo suggests this Fe is ferric (Figure 7). The PG Fo appears least altered, with minimal formation of dissolution features in SEM (Figure 5), a Mg/Si ratio inconsistent with leaching (~2) (Figure 6), and VNIR spectra almost identical to the control sample. Analysis of Fa surfaces is ongoing. The higher temperatures and more acidic pH in the KM soils likely drive the faster dissolution of the Fo disks described above. While the TB soils experience greater precipitation than in the KM, the cooler temperatures and more basic soil pH facilitate observable but more limited alteration. The dry climate and basic soil pH at PG lead to minimal dissolution of the PG disk surfaces.

A D Feldman↗

Evaluation of Zero-Net-Rate Pumping Tests

Accurately estimating the distribution of aquifer properties is key to understanding contaminant movement in the subsurface. The distribution of aquifer properties is typically addressed using slug or constant-rate well tests, and the pros and cons of these tests are well known. Slug tests are appealing because they avoid removing contaminated water, but their results are affected by well skin and the small volume of displaced water limits the volume of aquifer that can be evaluated. Constant-rate well tests have the disadvantage of requiring disposal of potentially contaminated water, but they can generate properties that are more representative than slug tests, and they can be used to estimate well efficiency and storativity, which are difficult to characterize using slug tests. Periodic pumping tests are appealing because they have many of the advantages and few of the disadvantages of slug and constant-rate well tests. Periodic pumping tests involve cycling the pumping rate with a regular period and measuring the resulting response in monitoring wells. Some of these tests involve imposing a periodic rate on a constant mean rate. However, other tests involve moving water out and back into the well at rates that are balanced so the net rate after each period is zero. Some zero-net-rate (ZNR) well tests use rates that follow a sinusoidal pattern, whereas others use a square wave pattern that switches between a constant rate of pumping and a constant rate of injection to achieve a zero-net rate. ZNR well tests appear to be a useful compromise between slug and constant rate tests, but methods for conducting and analyzing the results from these tests have received limited evaluation. Periodic pumping tests are typically evaluated through recorded pressure responses to a disturbance source. Recent work has demonstrated that the strain field in the vadose zone shows a response to disturbances in the underlying aquifer. Various instrument designs are available that can record vertical and horizontal strain to very precise resolutions (10-9 ε). Measuring the strain in the vadose zone can present lower costs than measuring pressure in a monitoring well, so using strain could improve the resolution of ZNR tests, as well as make these tests potentially cheaper. A ZNR periodic pumping system was constructed and used to perform tests in the Clemson, South Carolina area. The system was designed to pump and inject with a periodic square-wave, so pumping occurs at a constant rate for half the period and is followed by injection at the same rate for half the period resulting in no net gain or loss of water from the aquifer. The system is designed to generate flow rates from 1 to 3.5 gallons per minute (gpm) with a capacity of 875-gallon per half period. Pressure and strain monitoring points are located in the vicinity of the pumping well. A series of 10 ZNR periodic pumping tests were conducted at the field site, with periods ranging from 2 to 540 minutes. Pressure data were collected and recorded in four monitoring wells in various locations around the pumping well. In addition, vertical and horizontal strain and tilt data were recorded in various locations around the pumping well. Traditional aquifer tests, including slug tests and constant-rate pumping tests were conducted at the field site to provide a baseline of aquifer parameter estimates that are compared the results of ZNR tests. Pressure data measured at monitoring wells during ZNR tests were analyzed utilizing an analytical solution (Streltsova, 1988) that assumes a confined aquifer and estimates hydraulic diffusivity and transmissivity using the time lag and amplitude of the pressure signal. The storativity can be separated out from the hydraulic diffusivity using the transmissivity estimate. The time-lag and amplitude of the pressure signals were estimated using a Fourier transform. The time-lag of the primary period of the pressure increases as a linear function of distance from the pumping well, and it increases as the square root of the period of the ZNR test over a range of periods spanning two orders of magnitudes. The first few harmonics follow similar trends, but the higher frequency harmonics diverge from this trend. Strain and tilt data from various instruments in the vadose and saturated zones are periodic during ZNR tests. The time lags of the strain components increase roughly linearly with distance from the well. Many of the time lags of the strain data from the vadose zone are similar to the time lags observed in the pressure data from a similar distance and pumping period. Theoretical experiments were performed to understand how pressure and strain responds to ZNR periodic pumping. The time-lags of the simulated pressure responses match the field data trends. The results were used to validate the Streltsova (1988) solution for periodic pumping to estimate hydraulic diffusivity. Hydraulic diffusivity was estimated to be 1.8 x 10-2 m2s-1 < Dh < 9.0 x 10-2 m2s-1 using pressure data from all the wells during ZNR tests. Transmissivity was estimated to be 0.8 x 10-4 m2s-1 < T < 4.3 x 10-4 m2s-1. Assuming Storativity = T/Dh gives 1.1 x 10-3 < S < 19 x10-3. Transmissivity is used to estimate the hydraulic conductivity, K, by dividing by the assumed aquifer thickness. Two constant-rate pumping tests were conducted and analyzed using two conceptual models: The Hantush (1961) solution was used to analyze data assuming confined conditions and the Neuman (1974) solution was used to analyze data assuming unconfined conditions. The Hantush (1961) solution estimated transmissivity to be 1.6 x 10-4 m2s-1 < T < 1.8 x 10-4 m2s-1 when using data from all the wells. Storativity was estimated to be 4.4 x 10-3 < S < 8.2 x 10-3. Assuming Dh = T/S gives 0.22 x 10-2 m2s-1 < Dh < 3.6 x 10-2 m2s-1. The Neuman (1974) solution estimated transmissivity to be 0.80 x 10-4 m2s-1 < T < 0.84 x 10-4 m2s-1. Total storativity (S + Sy) was estimated to be 22 x 10-3 < S < 110 x 10-3. Storativity excluding the Sy term was estimated to be 2.2 x 10-3 < S < 2.5 x 10-3. Assuming Dh = T/S gives 0.073 x 10-2 m2s-1 < S < 0.36 x 10-2 m2s-1. Two slug tests were performed on the pumping well, PW-2. The Bouwer and Rice (1976) solution for slug tests was used to analyze the data. This solution works for both confined and unconfined conditions, so both conceptual models were covered through the analysis. Hydraulic conductivity was estimated to be 1.95 x 10-6 ms-1 < K < 2.49 x 10-6 ms-1. Statistical comparison show generally no significant difference in the estimates of Dh, K, or S made using ZNR and conventional pumping tests.. However, the total storage estimated by the Neuman (1974) solution for unconfined settings is larger than that estimate using the Streltsove analysis applied to the ZNR tests. This is likely because the Neuman solution considers delayed yield from storage that the water table, whereas the Streltsova (1988) solution assumes confined conditions. The analysis also indicates that there is a statistically significant difference between parameters estimated with conventional slug tests and those measured with either constant-rate pumping tests or ZNR tests. Hydraulic diffusivity is estimated from the time-lag and the distance from the well. The time-lags of the strain data from the vadose zone are similar to the lags of the pressure data from a similar distance. This indicates that hydraulic diffusivities estimated from the strain data measured in the vadose zone would be similar to the estimates from the pressure data measured in the aquifer. There appear to be sign reversals in some of the strain data that were corrected to estimate the time lag. The dissolved oxygen (CO) concentration was measured during several ZNR tests to evaluate the feasibility of using the procedure to increase contaminant degradation kinetics that are related to DO. The results indicate that DO can be increased during the cycling of ZNR tests in some cases.

Smith-Jones, Austin↗

A Methodology for the Analysis of Water Oxidation Electrocatalysts in the Absence of Limiting Current that Avoids the Pitfalls of Existing Methods

Water oxidation is an important reaction studied as a way to generate electrons from water, to promote water splitting and the formation of green hydrogen. When using electrodes to drive homogeneous water oxidation catalysis, cyclic voltammograms are analyzed to provide catalytic rate constants. There are two main methods, foot-of-the-wave analysis (FOWA) and limiting current analysis. FOWA relies on approximations inherent to analyzing water oxidation catalysis, such as determining the formal potential of the catalytic intermediate, E 0 cat . Limiting current methods are the optimal way to analyze catalyst performance but rely on observable limiting current, which is virtually never seen in water oxidation. To avoid those issues, a method is proposed for analyzing nonideal cyclic voltammetry waveshapes in water oxidation: by analyzing rate data across a large range of potentials, an optimal potential, E 0 cat , can be obtained, where catalytic current, i cat , is nearly independent of scan rate and has a linear dependency on buffer concentration. Here, the method is applied to four homogeneous water oxidation catalysts with prior extensive electrochemical elucidation, all of which lack an ideal, purely kinetic waveshape in cyclic voltammetry. Application of the method avoids the biases of the other methods cited for the kinetic analyses of water oxidation catalysts.

14 SOLAR ENERGY↗

Dissolved Inorganic Carbon and Dissolved Organic Carbon Data for the East River Watershed, Colorado (2015-2025)

This data package contains mean values for dissolved organic carbon (DOC) and dissolved inorganic carbon (DIC) for water samples taken from the East River Watershed in Colorado. The East River is part of the Watershed Function Scientific Focus Area (WFSFA) located in the Upper Colorado River Basin, United States. DOC and DIC concentrations in water samples were determined using a TOC-VCPH analyzer (Shimadzu Corporation, Japan). DOC was analyzed as non-purgeable organic carbon (NPOC) by purging HCl-acidified samples with carbon-free air to remove DIC prior to measurement. After the acidified sample has been sparged, it is injected into a combustion tube filled with oxidation catalyst heated to 680 oC. The DOC in samples is combusted to CO2 and measured by a non-dispersive infrared (NDIR) detector. The peak area of the analog signal produced by the NDIR detector is proportional to the DOC concentration of the sample. DIC was determined by acidifying the samples with HCl first, and then purging with carbon-free air to release CO2 for analysis by NDIR detector. Total dissolved nitrogen (TDN) was analyzed using a Shimadzu Total Nitrogen Module (TNM-L) combined with the TOC-L analyzer (Shimadzu Corporation, Japan). TNM-L is a non-specific measurement of total nitrogen (TN). All nitrogen species in samples are combusted to nitrogen monoxide and nitrogen dioxide, then reacted with ozone to form an excited state of nitrogen dioxide. Upon returning to ground state, light energy is emitted. Then, TDN is measured using a chemiluminescence detector. All data reported are the mean values upon minimum of three replicate measurements, with a relative standard deviation < 3%. All samples were analyzed under a rigorous quality assurance and quality control (QA/QC) process. This data package contains (1) a zip file (dic_npoc_data_2014-2025.zip) containing a total of 337 files: 336 data files of DIC and NPOC data from across the Lawrence Berkeley National Laboratory (LBNL) Watershed Function Scientific Focus Area (SFA) which is reported in .csv files per location and a locations.csv (1 file) with latitude and longitude for each location; (2) a file-level metadata (v6_20250901_flmd.csv) file that lists each file contained in the dataset with associated metadata; (3) a data dictionary (v6_20250901_dd.csv) file that contains terms/column_headers used throughout the files along with a definition, units, and data type; and (4) PDF and docx files for the determiniation of Method Detection Limits (MDLs) for DIC and NPOC data, which has been updated in 2026-08. Missing values within the anion data files are noted as either "-9999" or "0.0" for not detectable (N.D.) data. There are a total of 113 locations containing DIC/NPOC data. Update on 2020-10-07: Updated the data files to remove times from the timestamps, so that only dates remain. The data values have not changed. Update on 2021-04-11: Added Determination of Method Detection Limits (MDLs) for DIC, NPOC and TDN Analyses document, which can be accessed as a PDF or with Microsoft Word.Update on 6/10/2022: versioned updates to this dataset was made along with these changes: (1) updated dissolved inorganic carbon and dissolved organic carbon data for all locations up to 2021-12-31, (2) removal of units from column headers in datafiles, (3) added row underneath headers to contain units of variables, (4) restructure of units to comply with CSV reporting format requirements, (5) added -9999 for empty numerical cells, and (6) the addition of the file-level metadata (flmd.csv) and data dictionary (dd.csv) were added to comply with the File-Level Metadata Reporting Format. Update on 2022-09-09: Updates were made to reporting format specific files (file-level metadata and data dictionary) to correct swapped file names, add additional details on metadata descriptions on both files, add a header_row column to enable parsing, and add version number and date to file names (v2_20220909_flmd.csv and v2_20220909_dd.csv). Update on 2023-08-08: Updates were made to both the data files and reporting format specific files. New available anion data was added, up until 2023-01-05. The file level metadata and data dictionary files were updated to reflect the additional data added. Update on 2024-03-11: Updates were made to both the data files and reporting format specific files. New available anion data was added, up until 2023-11-21. Further, revisions to the data files were made to remove incorrect data points (from 1970 and 2001). The reporting format specific files were updated to reflect the additional data added. Revised versions of the PDF and docx files for determination of MDLs for DIC and NPOC were added to replace previous versions. Update on 2025-05-15: Updates were made to both the data files and reporting format specific files. New available DIC and NPOC data was added, up until the end of WY2024 (September 30, 2024). International Generic Sample Numbers (IGSNs), when registered, were added to the data files. The reporting format specific files were updated to reflect the additional data added. Update on 2026-09-01: Updates were made to both the data files and reporting format specific files. New available DIC and NPOC data was added, up until the end of WY2025 (September 30, 2025). Updated versions, as of 2026-08-10, of the PDF and docx files for determination of MDLs for DIC and NPOC data were added to this dataset.

54 ENVIRONMENTAL SCIENCES↗

Solar Photovoltaic (PV) Manufacturing Expansions in the United States, 2017-2019: Motives, Challenges, Opportunities, and Policy Context

This report analyzes U.S. PV manufacturing announcements from 2017 through July-2019 and explores the potential impact of tariffs on the competitiveness of U.S. PV module manufacturing. First, the fraction of announcements that resulted in established facilities in the time frame specified is analyzed, as well as firm nationality and the existence of any other U.S. facilities for the firm. The prevalence of different PV product types among the announcements was also analyzed, as well as the type of facilities selected for development (greenfield, existing PV factories, or repurposed buildings). The cost to build these facilities was also analyzed in the context of their capacity, the value of incentives received by state and local governments, and the number of jobs created. Case studies of Auxin Solar, First Solar, Hanwha Q Cells, JinkoSolar, LG, Silfab, Solaria, and SunPower were conducted and highlight five factors influencing firms’ decisions to expand their U.S. PV production capacity. These included: 1) proximity to demand, including supply contracts and demand for domestic content 2) tariffs and other incentives. The Section 201 tariffs on PV module imports were often described helpful but tariffs on cell imports (Section 201) or materials (Section 301, Section 232) were often described as counterproductive since U.S. supply chains do not exist for many necessary goods. The corporate tax decrease and local incentives were often mentioned as helpful. 3) the ability to move quickly in order to maximize the benefit of the Section 201 tariffs. This affected decisions to select localities with faster permitting or moving into existing buildings. 4) competition with scale, where competing with Asian PV imports or avoiding imports of material from China was challenging due to the larger scale of Asian manufacturers. Additionally, smaller firms expressed concerns that larger firms may not need their U.S. facilities to profitable. 5) access to capital, which was a particular challenge for smaller firms. Several U.S. tariffs affect the economics of PV cells, modules, and balance-of-module (BOM) materials, some of which are applicable to a long list of countries and others applicable only to China or Taiwan. We analyze how the tariffs interact to affect the prices of c-Si modules sourced under various assumptions about domestic and imported components and assembly. We estimate module prices under a range of baseline scenarios featuring U.S.-assembled modules with Southeast Asian cells (not from China or Taiwan). When all BOM components are sourced from the United States, baseline module prices are $0.37/W (below the 2.5-GW Section 201 cell cap) to $0.42/W (above the 2.5-GW cell cap). If all BOM components are sourced from China, the module price in this baseline U.S.-assembly scenario can reach up to $0.60/W. Imported module prices are lower than or within the range of the U.S.-assembled baseline prices under some conditions. Importing bifacial modules from Southeast Asia is the lowest-priced alternative, at $0.35/W, because bifacial modules are excluded from the Section 201 tariffs. Importing Southeast Asia modules that contain U.S. cells results in a price of $0.38/W, also owing to a Section 201 exclusion. Importing Southeast Asian modules that have Southeast Asian cells yields a price of $0.44/W, and importing Chinese modules that have U.S. cells yields a price of $0.58/W. These results suggest that U.S. tariff policy in 2019 may encourage U.S. assembly of mono-facial c-Si modules if the 2.5-GW cell import cap is not reached and Chinese BOM imports are negligible. However, some imported module options may offer the lowest prices if the 2.5-GW cell import cap is reached and/or U.S. module manufacturers must rely on significant Chinese BOM content. The analysis shows the 2.5-GW cell cap may be exceeded in 2019 and following years. It is also likely that U.S. module manufacturers will need to use at least some Chinese BOM content, since U.S. module manufacturers have stated that U.S. supply chains do not exist (or do not exist at a sufficient scale) for most BOM components, which makes them at least partially reliant on Chinese BOM and thus subject to Section 301 tariffs.

14 SOLAR ENERGY↗