Results for “Spatial Characterization” — Page 15

DOE OSTI2022

Autonomous measurements of carbon dioxide uptake in a blue carbon habitat

A better understanding of carbon cycling dynamics in coastal ecosystems is critical to quantify “blue carbon” storage. However, obtaining the necessary measurements is complicated due to the dynamic conditions that characterize these ecosystems, and requires both spatially and temporally resolved measurements of carbon dioxide concentrations and complementary water quality parameters. In this study, we evaluate the use of a utility-class autonomous surface vehicle as a platform for measuring carbon dioxide dynamics in a shallow, intertidal, blue carbon ecosystem. We integrated a carbon dioxide sensor, a water quality sonde, and hydroacoustic sensors with an autonomous surface vehicle and performed measurements in and around a meadow of Zostera marina eelgrass. Results demonstrate the utility of this system to provide insight into complex carbon cycling dynamics in and around a shallow, intertidal blue carbon habitat.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

TRACER-BC2-SP2 Field Campaign Report

The focus of this project was to assess black carbon and brown carbon aerosol in Houston, Texas during the U.S. Department of Energy’s Atmospheric Radiation Measurement (ARM) user facility Tracking Aerosol Convection Interactions Experiment (TRACER) through the co-located deployment of a single-particle soot photometer (SP2) with the University of Houston/Baylor aerosol optical instrumentation in the Texas (BC)2 network in and around Houston. The deployment of the SP2 enables us to expand the spatial assessment of black carbon (BC) in Houston during TRACER by permitting a more thorough characterization of BC, including sources, transport, and potential processing. BC is well known to undergo physical and chemical changes during local and long-range transport that will impact its optical properties. By deploying multiple SP2s in the Houston area, the impacts of local and long-range transport can be assessed. The TRACER-BC2-SP2 campaign will expand spatial and temporal characterization of BC in the Houston area during TRACER and link the TRACER-carbon analysis to ongoing black and brown carbon monitoring (Texas [BC]2 network) in Houston.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Determining Optimal Magnetometer Configuration on MAGIS-100

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.] (ORC↗

DOE OSTI2026

Determining Optimal Magnetometer Configuration on MAGIS-100

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.] (ORC↗

DOE OSTI2021

Tradeoffs and Synergies in Tropical Forest Root Traits and Dynamics for Nutrient and Water Acquisition: Field and Modeling Advances

Vegetation processes are fundamentally limited by nutrient and water availability, the uptake of which is mediated by plant roots in terrestrial ecosystems. While tropical forests play a central role in global water, carbon, and nutrient cycling, we know very little about tradeoffs and synergies in root traits that respond to resource scarcity. Tropical trees face a unique set of resource limitations, with rock-derived nutrients and moisture seasonality governing many ecosystem functions, and nutrient versus water availability often separated spatially and temporally. Root traits that characterize biomass, depth distributions, production and phenology, morphology, physiology, chemistry, and symbiotic relationships can be predictive of plants’ capacities to access and acquire nutrients and water, with links to aboveground processes like transpiration, wood productivity, and leaf phenology. In this review, we identify an emerging trend in the literature that tropical fine root biomass and production in surface soils are greatest in infertile or sufficiently moist soils. We also identify interesting paradoxes in tropical forest root responses to changing resources that merit further exploration. For example, specific root length, which typically increases under resource scarcity to expand the volume of soil explored, instead can increase with greater base cation availability, both across natural tropical forest gradients and in fertilization experiments. Also, nutrient additions, rather than reducing mycorrhizal colonization of fine roots as might be expected, increased colonization rates under scenarios of water scarcity in some forests. Efforts to include fine root traits and functions in vegetation models have grown more sophisticated over time, yet there is a disconnect between the emphasis in models characterizing nutrient and water uptake rates and carbon costs versus the emphasis in field experiments on measuring root biomass, production, and morphology in response to changes in resource availability. Closer integration of field and modeling efforts could connect mechanistic investigation of fine-root dynamics to ecosystem-scale understanding of nutrient and water cycling, allowing us to better predict tropical forest-climate feedbacks.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Effects of the large-scale circulation on temperature and water vapor distributions in the Π Chamber

Abstract. Microphysical processes are important for the development of clouds and thus Earth's climate. For example, turbulent fluctuations in the water vapor mixing ratio, r, and temperature, T, cause fluctuations in the saturation ratio, S. Because S is the driving factor in the condensational growth of droplets, fluctuations may broaden the cloud droplet size distribution due to individual droplets experiencing different growth rates. The small-scale turbulent fluctuations in the atmosphere that are relevant to cloud droplets are difficult to quantify through field measurements. We investigate these processes in the laboratory using Michigan Tech's Π Chamber. The Π Chamber utilizes Rayleigh–Bénard convection (RBC) to create the turbulent conditions inherent in clouds. In RBC it is common for a large-scale circulation (LSC) to form. As a consequence of the LSC, the temperature field of the chamber is not spatially uniform. In this paper, we characterize the LSC in the Π Chamber and show how it affects the shape of the distributions of r, T, and S. The LSC was found to follow a single roll with an updraft and downdraft along opposing walls of the chamber. Near the updraft (downdraft), the distributions of T and r were positively (negatively) skewed. At each measuring position, S consistently had a negatively skewed distribution, with the downdraft being the most negative.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Technical note: Using long short-term memory models to fill data gaps in hydrological monitoring networks

Abstract. Quantifying the spatiotemporal dynamics in subsurface hydrological flows over a long time window usually employs a network of monitoring wells. However, such observations are often spatially sparse with potential temporal gaps due to poor quality or instrument failure. In this study, we explore the ability of recurrent neural networks to fill gaps in a spatially distributed time-series dataset. We use a well network that monitors the dynamic and heterogeneous hydrologic exchanges between the Columbia River and its adjacent groundwater aquifer at the U.S. Department of Energy's Hanford site. This 10-year-long dataset contains hourly temperature, specific conductance, and groundwater table elevation measurements from 42 wells with gaps of various lengths. We employ a long short-term memory (LSTM) model to capture the temporal variations in the observed system behaviors needed for gap filling. The performance of the LSTM-based gap-filling method was evaluated against a traditional autoregressive integrated moving average (ARIMA) method in terms of error statistics and accuracy in capturing the temporal patterns of river corridor wells with various dynamics signatures. Our study demonstrates that the ARIMA models yield better average error statistics, although they tend to have larger errors during time windows with abrupt changes or high-frequency (daily and subdaily) variations. The LSTM-based models excel in capturing both high-frequency and low-frequency (monthly and seasonal) dynamics. However, the inclusion of high-frequency fluctuations may also lead to overly dynamic predictions in time windows that lack such fluctuations. The LSTM can take advantage of the spatial information from neighboring wells to improve the gap-filling accuracy, especially for long gaps in system states that vary at subdaily scales. While LSTM models require substantial training data and have limited extrapolation power beyond the conditions represented in the training data, they afford great flexibility to account for the spatial correlations, temporal correlations, and nonlinearity in data without a priori assumptions. Thus, LSTMs provide effective alternatives to fill in data gaps in spatially distributed time-series observations characterized by multiple dominant frequencies of variability, which are essential for advancing our understanding of dynamic complex systems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

In-Situ Magnetic Field Reconstruction in the MAGIS-100 Experiment

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.; Ferm↗

DOE OSTI2020

Characterizing Peatland Microtopography Using Gradient and Microform-Based Approaches

Abstract Peatlands represent an important component of the global carbon cycle, storing 180–621 Gt of carbon (C). Small-scale spatial variations in elevation, frequently referred to as microtopography, influence ecological processes associated with the peatland C cycle, including Sphagnum photosynthesis and methane flux. Microtopography can be characterized with measures of topographic variability and by using conceptual classes (microforms) linked to function: most commonly hummocks and hollows. However, the criteria used to define these conceptual classes are often poorly described, if at all, and vary between studies. Such inconsistencies compel development of explicit quantitative methods to classify microforms. Furthermore, gradient-based characterizations that describe spatial variability without the use of microforms are lacking in the literature. Therefore, the objectives of this study were to (1) calculate peatland microtopographical elevation gradients and measures of spatial variability, (2) develop three microform classification methods intended for specific purposes, and (3) evaluate and contrast classification methods. Our results suggest that at spatial scales much larger than microforms, elevation distributions are unimodal and are well approximated with parametric probability density functions. Results from classifications were variable between methods and years and exhibited significant differences in mean hollow areal coverages of a raised ombrotrophic bog. Our results suggest that the conceptualization and classification of microforms can significantly influence microtopographic structural metrics. The three explicit methods for microform classification described here may be used and built upon for future applications.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Tunable chiral and nematic states in the triple-Q antiferromagnet Co 1/3 TaS 2

Complex spin configurations in magnetic materials, ranging from collinear single-Q to non-coplanar multi-Q states, exhibit rich symmetry and chiral properties. However, their detailed characterization is often hindered by the limited spatial resolution of neutron diffraction techniques. Here we employ magnetic circular dichroism and magnetic linear dichroism to investigate the triangular lattice antiferromagnet Co 1/3 TaS 2 , revealing three-state (Z 3 ) nematicity and also spin chirality across its multi-Q magnetic phases. At intermediate temperatures, the presence of linear dichroism identifies nematicity arising from a single-Q stripe phase, while at high magnetic fields and low temperatures, a phase characterized solely by circular dichroism emerges, signifying a purely chiral non-coplanar triple-Q state. Notably, at low temperatures and small fields, we discover a unique phase where both chirality and nematicity coexist. A theoretical analysis based on a continuous multi-Q manifold captures the emergence of these distinct magnetic phases, as a result of the interplay between four-spin interactions and weak magnetic anisotropy. Additionally, both circular and linear dichroism microscopy spatially resolves the chiral and nematic domains. Our findings establish Co 1/3 TaS 2 as a rare platform hosting diverse multi-Q states with distinct combinations of spin chirality and nematicity while demonstrating the effectiveness of polarized optical techniques in characterizing complex magnetic textures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Time-resolved X-ray imaging of the current filamentation instability in solid-density plasmas

The streaming of energetic charged particles can magnetize astrophysical and laboratory plasmas via the current filamentation instability. Despite its importance, the experimental characterization of this instability has remained a challenge. Here, we report an experiment combining a high-intensity optical laser with a high-brightness X-ray free electron laser that successfully images the instability in solid-density plasmas with 200 nm spatial and 50 fs temporal resolution. We characterize the development of μm-scale filamentary structures and their evolution over tens of picoseconds through a non-linear merging process. The measured plasma density modulations and long merging time reveal the critical importance of space-charge effects and ion motion on this electron-driven instability. Supporting theoretical analysis and kinetic simulations help distinguish the relative role of space-charge and resistive effects. Our findings indicate that magnetic fields on the order of 10 megagauss are produced, with important implications for transport and radiation emission of energetic particles in plasmas.

Schoenwaelder, Christopher [SLAC National Accelera↗

DOE OSTI2021

Talbot-Lau x-ray deflectometer: Refraction-based HEDP imaging diagnostic

Talbot-Lau x-ray interferometry has been implemented to map electron density gradients in High Energy Density Physics (HEDP) experiments. X-ray backlighter targets have been evaluated for Talbot-Lau X-ray Deflectometry (TXD). Cu foils, wires, and sphere targets have been irradiated by 10–150 J, 8–30 ps laser pulses, while two pulsed-power generators (∼350 kA, 350 ns and ∼200 kA, 150 ns) have driven Cu wire, hybrid, and laser-cut x-pinches. A plasma ablation front generated by the Omega EP laser was imaged for the first time through TXD for densities >1023 cm−3. Backlighter optimization in combination with x-ray CCD, image plates, and x-ray film has been assessed in terms of spatial resolution and interferometer contrast for accurate plasma characterization through TXD in pulsed-power and high-intensity laser environments. The results obtained thus far demonstrate the potential of TXD as a powerful diagnostic for HEDP.

Valdivia, M. P.↗

DOE OSTI2025

A self-supervised robotic system for autonomous contact-based spatial mapping of semiconductor properties

Integrating robotically driven contact-based material characterization techniques into self-driving laboratories can enhance measurement quality, reliability, and throughput. While deep learning models support robust autonomy, current methods lack reliable pixel-precision positioning and require extensive labeled data. To overcome these challenges, we propose an approach for building self-supervised autonomy into contact-based robotic systems that teach the robot to follow domain expert measurement principles at high throughputs. We demonstrate the performance of this approach by autonomously driving a 4-DOF robotic probe for 24 hours to characterize semiconductor photoconductivity at 3025 uniquely predicted poses across a gradient of drop-casted perovskite film compositions, achieving throughputs of more than 125 measurements per hour. Spatially mapping photoconductivity onto each drop-casted film reveals compositional trends and regions of inhomogeneity, valuable for identifying manufacturing defects. With this self-supervised neural network–driven robotic system, we enable high-precision and reliable automation of contact-based characterization techniques at high throughputs, thereby allowing measurement of previously inaccessible yet important semiconductor properties for self-driving laboratories.

Science & Technology - Other Topics↗

DOE OSTI2020

Energy storage emerging: A perspective from the Joint Center for Energy Storage Research

Energy storage is an integral part of modern society. A contemporary example is the lithium (Li)-ion battery, which enabled the launch of the personal electronics revolution in 1991 and the first commercial electric vehicles in 2010. Most recently, Li-ion batteries have expanded into the electricity grid to firm variable renewable generation, increasing the efficiency and effectiveness of transmission and distribution. Important applications continue to emerge including decarbonization of heavy-duty vehicles, rail, maritime shipping, and aviation and the growth of renewable electricity and storage on the grid. This perspective compares energy storage needs and priorities in 2010 with those now and those emerging over the next few decades. The diversity of demands for energy storage requires a diversity of purpose-built batteries designed to meet disparate applications. Advances in the frontier of battery research to achieve transformative performance spanning energy and power density, capacity, charge/discharge times, cost, lifetime, and safety are highlighted, along with strategic research refinements made by the Joint Center for Energy Storage Research (JCESR) and the broader community to accommodate the changing storage needs and priorities. Innovative experimental tools with higher spatial and temporal resolution, in situ and operando characterization, first-principles simulation, high throughput computation, machine learning, and artificial intelligence work collectively to reveal the origins of the electrochemical phenomena that enable new means of energy storage. This knowledge allows a constructionist approach to materials, chemistries, and architectures, where each atom or molecule plays a prescribed role in realizing batteries with unique performance profiles suitable for emergent demands.

Trahey, Lynn↗

DOE OSTI2023

Multiscale spatial analysis of fracture nodes in two dimensions

Spatial arrangement of fractures as a function of scale is an important component of fracture quantification for inferential and predictive modeling. Available methods that analyze fracture spatial arrangement are based on one-dimensional spacing data; therefore, they are limited to semi-parallel fractures. Such methods cannot be applied to fracture networks in higher dimensions, particularly when fractures have different orientations. Here to characterize fracture arrangements in two dimensions, we propose using Ripley’s K-function, as a method of point pattern analysis, to quantify spatial arrangement of fracture nodes. Fracture nodes, such as barycenters, intersection points, and tips, are point-based representations of fracture locations and connectivity within the fracture network. We introduce formulations for isotropic as well as directional analyses of spatial arrangement. In addition, we derive formulations for edge correction in circular and rectangular study domains. Finally, we demonstrate applications of Ripley’s K-function on two natural fracture datasets. Our proposed method supports quantification and characterization of fracture spatial arrangements that allow practitioners to build representative models of fractures in the subsurface.

02 PETROLEUM↗

DOE OSTI2020

Fill Factor Loss in Fielded Photovoltaic Modules Due to Metallization Failures, Characterized by Luminescence and Thermal Imaging

Understanding the causes of photovoltaic module failure in the field can be achieved using spatially-resolved imaging tools in combination with current-voltage (IV) characterization. We evaluate 23 field-degraded modules by comparing their IV characteristics with module images obtained by electroluminescence, photoluminescence, UV-fluorescence, and dark lock-in thermography. We identify primarily metallization failures (metallization series resistance, solder point heating, broken fingers, disrupted interconnects) in these modules that has led to a severe drop in the fill factor and up to 35% power loss over only 2 years of deployment.

14 SOLAR ENERGY↗

DOE OSTI2025

A quantum eigenvalue solver based on tensor networks

Electronic ground states are of central importance in chemical simulations, but have remained beyond the reach of efficient classical algorithms except in cases of weak electron correlation or one-dimensional spatial geometry. We introduce a hybrid quantum-classical eigenvalue solver that constructs a wavefunction ansatz from a linear combination of matrix product states in rotated orbital bases, enabling the characterization of strongly correlated ground states with arbitrary spatial geometry. The energy is converged via a gradient-free generalized sweep algorithm based on quantum subspace diagonalization, with a potentially exponential speedup in the off-diagonal matrix element contractions upon translation into compact quantum circuits of linear depth in the number of qubits. Chemical accuracy is attained in numerical experiments for both a stretched water molecule and an octahedral arrangement of hydrogen atoms, achieving substantially better correlation energies compared to a unitary coupled-cluster benchmark, with orders of magnitude reductions in quantum resource estimates and a surprisingly high tolerance to shot noise. This proof-of-concept study suggests a promising new avenue for scaling up simulations of strongly correlated chemical systems on near-term quantum hardware.

chemistry↗

DOE OSTI2020

Characterizing Ocean Acidification and Atmospheric Emission caused by Methane Released from Gas Hydrate Systems along the US Atlantic Margin (Final Report)

To assess the environmental impact of decomposing gas hydrates, this project characterized the atmospheric emission and ocean acidification associated with seafloor-released methane along the U.S. Atlantic Margin. This project leveraged the newly discovered seafloor seeps in this region, many of which are located near the landward limit of gas hydrate stability, to investigate the fate and environmental impact of methane once released from the seafloor. Several research directions were pursued during these investigations. First, high-spatial resolution surface surveys were conducted to characterize the degree to which sea-to-air emission of methane is correlated with acoustically detected seafloor bubble emissions. Second, the natural isotopic signatures of both radiocarbon and δ 13 C were measured throughout the water column to observe if seafloor-released methane was detected in surface waters or if it was replaced by a different source. Third, the extent of methane oxidation in the water column was constrained to assess the ability of this process to mitigate atmospheric methane emissions. Fourth, since aerobic methane oxidation ultimately converts methane to carbon dioxide, the extent of ocean acidification relating to this process was assessed. Ultimately, this project was conducted in four stages. Stage one involved preparations for the research expedition along the U.S. Atlantic Margin. This involved calibrating, modifying, and validating equipment needed for sample and data collection in this unique environment. Stage two involved executing the at-sea science. This research was conducted from 25 August to 5 September 2017 on the R/V Hugh R. Sharp. This expedition left from Lewes, Delaware and investigated the continental slope from approximately this latitude down to Cape Hatteras, North Carolina. Third, any samples unanalyzed on the research vessel were returned to the laboratory and measured. Fourth, the data were interpreted, and publications and presentations were prepared.

03 NATURAL GAS↗