Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data depth”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

SU-F-T-521: Flattening-Filter-Free Beam Parameters Comparison From Different Linac Machine Types

Purpose: Novel linac machines, TrueBeam (TB) and Elekta Versa have updated head designing and software control system, include flattening-filter-free (FFF) photon and electron beams. Later on FFF beams were also introduced on C-Series machines. In this work FFF beams for same energy 6MV but from different machine versions were studied with reference to beam data parameters. Methods: The 6MV-FFF percent depth doses, profile symmetry and flatness, dose rate tables, and multi-leaf collimator (MLC) transmission factors were measured during commissioning process of both C-series and Truebeam machines. The scanning and dosimetric data for 6MV-FFF beam from Truebeam and C-Series linacs was compared. A correlation of 6MV-FFF beam from Elekta Versa with that of Varian linacs was also found. Results: The scanning files were plotted for both qualitative and quantitative analysis. The dosimetric leaf gap (DLG) for C-Series 6MV-FFF beam is 1.1 mm. Published values for Truebeam dosimetric leaf gap is 1.16 mm. 6MV MLC transmission factor varies between 1.3 % and 1.4 % in two separate measurements and measured DLG values vary between 1.32 mm and 1.33 mm on C-Series machine. MLC transmission factor from C-Series machine varies between 1.5 % and 1.6 %. Some of the measured data values from C-Series FFF beam are compared with Truebeam representative data. 6MV-FFF beam parameter values like dmax, OP factors, beam symmetry and flatness and additional parameters for C-Series and Truebeam liancs will be presented and compared in graphical form and tabular data form if selected. Conclusion: The 6MV flattening filter (FF) beam data from C-Series & Truebeam and 6MV-FFF beam data from Truebeam has already presented. This particular analysis to compare 6MV-FFF beam from C-Series and Truebeam provides opportunity to better elaborate FFF mode on novel machines. It was found that C-Series and Truebeam 6MV-FFF dosimetric and beam data was quite similar.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

XO-2b: A HOT JUPITER WITH A VARIABLE HOST STAR THAT POTENTIALLY AFFECTS ITS MEASURED TRANSIT DEPTH

The transiting hot Jupiter XO-2b is an ideal target for multi-object photometry and spectroscopy as it has a relatively bright (V-mag = 11.25) K0V host star (XO-2N) and a large planet-to-star contrast ratio (R{sub p}/R{sub s} ≈ 0.015). It also has a nearby (31.″21) binary stellar companion (XO-2S) of nearly the same brightness (V-mag = 11.20) and spectral type (G9V), allowing for the characterization and removal of shared systematic errors (e.g., airmass brightness variations). We have therefore conducted a multiyear (2012–2015) study of XO-2b with the University of Arizona’s 61″ (1.55 m) Kuiper Telescope and Mont4k CCD in the Bessel U and Harris B photometric passbands to measure its Rayleigh scattering slope to place upper limits on the pressure-dependent radius at, e.g., 10 bar. Such measurements are needed to constrain its derived molecular abundances from primary transit observations. We have also been monitoring XO-2N since the 2013–2014 winter season with Tennessee State University’s Celestron-14 (0.36 m) automated imaging telescope to investigate stellar variability, which could affect XO-2b’s transit depth. Our observations indicate that XO-2N is variable, potentially due to cool star spots, with a peak-to-peak amplitude of 0.0049 ± 0.0007 R-mag and a period of 29.89 ± 0.16 days for the 2013–2014 observing season and a peak-to-peak amplitude of 0.0035 ± 0.0007 R-mag and 27.34 ± 0.21 day period for the 2014–2015 observing season. Because of the likely influence of XO-2N’s variability on the derivation of XO-2b’s transit depth, we cannot bin multiple nights of data to decrease our uncertainties, preventing us from constraining its gas abundances. This study demonstrates that long-term monitoring programs of exoplanet host stars are crucial for understanding host star variability.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Prestack reverse time migration for 3D marine reflection seismic data

Prestack reverse time migration (RTM) is a method for imaging the subsurface using the inner product of wavefield extrapolation in shot domain and in receiver domain. It is well known that RTM is better for preserving amplitudes and phases than other prestack migrations. Since 3D seismic data is huge data volume and it needs heavy computing works, it requires parallel computing in order to have a meaningful depth image of the 3D subsurface. We implemented a parallelized version of 3D RTM for prestack depth migration. The results of numerical example for 3D SEG/EAGE salt model showed good agreement with the original geological model. We applied RTM to offshore 3D seismic reflection data. The study area is 12 × 25 km with 120 survey lines. Shot and receiver spacing is 25 m and 12.5 m. The line spacing is 100 m. Shot gathers were preprocessed to enhance signal to noise ratio and velocity model was calculated from conventional stack velocity. Both of them were used to obtain 3D image using RTM. The results show reasonable subsurface image.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The Effects of Assumed Source Depth and Shear-Wave Velocity on Moment Tensors Estimated for Small, Contained Chemical Explosions in Granite

ABSTRACT The Source Phenomenology Experiment (SPE-Arizona) included of a series of chemical explosions detonated within a copper mine in Arizona. This study focuses on ground motions from detonations in the copper mine, which are analyzed to assess the uniqueness of the resulting source representation when the source region propagation characteristics have a range of possible models. P-wave velocities are well constrained by refraction data with less constraint of the S-wave velocities. The effects of explosion source depth and VS are assessed with Green’s functions for a range of models in which VP is held constant. Propagation models with a Poisson’s value of 0.25 and a source depth 30–60 m most accurately replicate the data. The explosion was detonated at a centroid depth of 30 m, so trade-offs in depth are demonstrated. The compensated linear vector dipole and explosion components of the Green’s functions convolved with a Mueller–Murphy source function are compared. Both produce significant energy in the 2–12 Hz band, due to surface-wave contributions with no clear depth dependencies above 20 Hz. The range of propagation models is used with the observational data to invert for the frequency-domain moment tensor. Fits to the data from these inversions have cross-correlation values of 0.64, demonstrating effectiveness in replicating the observations with the assumed propagation path effects and resulting source function. Inversions produce horizontal dipoles (Mxx and Myy), roughly half the maximum amplitude of Mzz, consistent with a compensated linear vector dipole source, which is frequency dependent. Denny and Johnson, Mueller–Murphy, Walter and Ford, and the revised Mueller–Murphy source models, parameterized for granite, are compared to the moment tensors. Despite a nonisotropic moment tensor source, the revised Mueller–Murphy isotropic source model best replicates the long-period moments, overshoot, and corner frequency.

Geochemistry & Geophysics↗
DOE OSTI2024

Utah FORGE: Fiber Optic Cumulative Strain Change and Strain Change Rate Data From Well 16A Stimulation at Well 16B

This dataset includes Rayleigh Frequency Shift (RFS) Distributed Strain Sensing (DSS) cumulative strain change and change rate data. The data was acquired during the stimulation of Utah FORGE Well 16A(78)-32 in April 2024 via fiber installed in Well 16B(78)-32. The fiber optic data was acquired using Neubrex SR7000 RFS DSS Distributed Strain sensing instruments and is saved here in the format of HDF5 files (.h5 extension). The spatial sampling on the full wellbore profiles is 0.20 centimeters. The data is the far field strain change response from a baseline profile made down the 16B well on April 3, 2024, so each strain value represents the strain change or strain change rate at each depth relative to the baseline reference profile. The data arrays for each type share the same dimensions (number of channels and time stamps).

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Thermal modeling of three lakes within the continuous permafrost zone in Alaska using the LAKE 2.0 model

Abstract. Lakes in the Arctic are important reservoirs of heat with much lower albedo in summer and greater absorption of solar radiation than surrounding tundra vegetation. In the winter, lakes that do not freeze to their bed have a mean annual bed temperature >0 ∘C in an otherwise frozen landscape. Under climate warming scenarios, we expect Arctic lakes to accelerate thawing of underlying permafrost due to warming water temperatures in the summer and winter. Previous studies of Arctic lakes have focused on ice cover and thickness, the ice decay process, catchment hydrology, lake water balance, and eddy covariance measurements, but little work has been done in the Arctic to model lake heat balance. We applied the LAKE 2.0 model to simulate water temperatures in three Arctic lakes in northern Alaska over several years and tested the sensitivity of the model to several perturbations of input meteorological variables (precipitation, shortwave radiation, and air temperature) and several model parameters (water vertical resolution, sediment vertical resolution, depth of soil column, and temporal resolution). The LAKE 2.0 model is a one-dimensional model that explicitly solves vertical profiles of water state variables on a grid. We used a combination of meteorological data from local and remote weather stations, as well as data derived from remote sensing, to drive the model. We validated modeled water temperatures with data of observed lake water temperatures at several depths over several years for each lake. Our validation of the LAKE 2.0 model is a necessary step toward modeling changes in Arctic lake ice regimes, lake heat balance, and thermal interactions with permafrost. The sensitivity analysis shows us that lake water temperature is not highly sensitive to small changes in air temperature or precipitation, while changes in shortwave radiation and large changes in precipitation produced larger effects. Snow depth and lake ice strongly affect water temperatures during the frozen season, which dominates the annual thermal regime of Arctic lakes. These findings suggest that reductions in lake ice thickness and duration could lead to more heat storage by lakes and enhanced permafrost degradation.

Clark, Jason A. (ORCID:0000000248105261)↗
DOE OSTI2017

Monte Carlo simulations of time-of-flight PET with double-ended readout: calibration, coincidence resolving times and statistical lower bounds

Here, this paper demonstrates through Monte Carlo simulations that a practical positron emission tomograph with (1) deep scintillators for efficient detection, (2) double-ended readout for depth-of-interaction information, (3) fixed-level analog triggering, and (4) accurate calibration and timing data corrections can achieve a coincidence resolving time (CRT) that is not far above the statistical lower bound. One Monte Carlo algorithm simulates a calibration procedure that uses data from a positron point source. Annihilation events with an interaction near the entrance surface of one scintillator are selected, and data from the two photodetectors on the other scintillator provide depth-dependent timing corrections. Another Monte Carlo algorithm simulates normal operation using these corrections and determines the CRT. A third Monte Carlo algorithm determines the CRT statistical lower bound by generating a series of random interaction depths, and for each interaction a set of random photoelectron times for each of the two photodetectors. The most likely interaction times are determined by shifting the depth-dependent probability density function to maximize the joint likelihood for all the photoelectron times in each set. Example calculations are tabulated for different numbers of photoelectrons and photodetector time jitters for three 3 × 3 × 30 mm 3 scintillators: Lu 2 SiO 5 :Ce,Ca (LSO), LaBr 3 :Ce, and a hypothetical ultra-fast scintillator. To isolate the factors that depend on the scintillator length and the ability to estimate the DOI, CRT values are tabulated for perfect scintillator-photodetectors. For LSO with 4000 photoelectrons and single photoelectron time jitter of the photodetector J = 0.2 ns (FWHM), the CRT value using the statistically weighted average of corrected trigger times is 0.098 ns FWHM and the statistical lower bound is 0.091 ns FWHM. For LaBr 3 :Ce with 8000 photoelectrons and J = 0.2 ns FWHM, the CRT values are 0.070 and 0.063 ns FWHM, respectively. For the ultra-fast scintillator with 1 ns decay time, 4000 photoelectrons, and J = 0.2 ns FWHM, the CRT values are 0.021 and 0.017 ns FWHM, respectively. Lastly, the examples also show that calibration and correction for depth-dependent variations in pulse height and in annihilation and optical photon transit times are necessary to achieve these CRT values.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

SURVIVAL DEPTH OF ORGANICS IN ICES UNDER LOW-ENERGY ELECTRON RADIATION ({<=}2 keV)

Icy surfaces in our solar system are continually modified and sputtered with electrons, ions, and photons from solar wind, cosmic rays, and local magnetospheres in the cases of Jovian and Saturnian satellites. In addition to their prevalence, electrons specifically are expected to be a principal radiolytic agent on these satellites. Among energetic particles (electrons and ions), electrons penetrate by far the deepest into the ice and could cause damage to organic material of possible prebiotic and even biological importance. To determine if organic matter could survive and be detected through remote sensing or in situ explorations on these surfaces, such as water ice-rich Europa, it is important to obtain accurate data quantifying electron-induced chemistry and damage depths of organics at varying incident electron energies. Experiments reported here address the quantification issue at lower electron energies (100 eV-2 keV) through rigorous laboratory data analysis obtained using a novel methodology. A polycyclic aromatic hydrocarbon molecule, pyrene, embedded in amorphous water ice films of controlled thicknesses served as an organic probe. UV-VIS spectroscopic measurements enabled quantitative assessment of organic matter survival depths in water ice. Eight ices of various thicknesses were studied to determine damage depths more accurately. The electron damage depths were found to be linear, approximately 110 nm keV{sup -1}, in the tested range which is noticeably higher than predictions by Monte Carlo simulations by up to 100%. We conclude that computational simulations underestimate electron damage depths in the energy region {<=}2 keV. If this trend holds at higher electron energies as well, present models utilizing radiation-induced organic chemistry in icy solar system bodies need to be revisited. For interstellar ices of a few micron thicknesses, we conclude that low-energy electrons generated through photoionization processes in the interstellar medium could penetrate through ice grains significantly and trigger organic reactions several hundred nanometers deep-bulk chemistry thus competing with surface chemistry of astrophysical ice grains.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Synthetic Absorption Lines from Simulations of Multiphase Gas in Galactic Winds

Supernova-driven galactic winds are multiphase streams of gas that are often observed flowing at a range of velocities out of star-forming regions in galaxies. In this study, we use high-resolution 3D simulations of multiphase galactic winds modeled with the hydrodynamics code Cholla to investigate the connection between numerical studies and observations. Using a simulated interaction between a hot $T$ ~ 10 4 K supernova-driven wind and a cool $T$ ~ 10 7 K cloud of interstellar material, we create mock observables, including the optical depth ($\tau$) and covering fraction ($C$ $f$ ) of six commonly observed ions (Si II , C II , Si IV , C IV , N V , and O VI ) as a function of gas velocity. We compare our mock observables to surveys of galactic winds in the literature, finding good agreement with velocities and profiles of the low ions. We then compute "empirical" values for the optical depth and covering fraction following observational techniques and compare them to the values calculated directly from the simulation data. We find that the empirically computed values tend to underestimate the "true" value of $\tau$ for ions with high optical depth and overestimate the "true" value of $\tau$ for ions with low optical depth relative to the simulated data. In conclusion, the empirically estimated covering fractions match our direct calculations very well for the low-ionization ions; for the high-ionization ions, the empirical covering fractions underestimate the directly computed values by up to ~40%.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Improved modeling of in-ice particle showers for IceCube event reconstruction

The IceCube Neutrino Observatory relies on an array of photomultiplier tubes to detect Cherenkov light produced by charged particles in the South Pole ice. IceCube data analyses depend on an in-depth characterization of the glacial ice, and on novel approaches in event reconstruction that utilize fast approximations of photoelectron yields. Here, a more accurate model is derived for event reconstruction that better captures our current knowledge of ice optical properties. When evaluated on a Monte Carlo simulation set, the median angular resolution for in-ice particle showers improves by over a factor of three compared to a reconstruction based on a simplified model of the ice. The most substantial improvement is obtained when including effects of birefringence due to the polycrystalline structure of the ice. When evaluated on data classified as particle showers in the high-energy starting events sample, a significantly improved description of the events is observed.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Citizen Green #1 Well: MICP and Pycnometry Data, Sidewall Cores

This submission includes information on helium pycnometry & MICP measurements conducted on sidewall sub-samples from the Citizen Green #1 well drilled on King Island, CA (API 07720688) as part of the WESTCARB partnership basin characterization effort. Cores were acquired with a rotary sidewall tool at a range of depths. See the sidewallCore submission for additional formation data.

Citizen Green↗
DOE OSTI2017

First Results from the Hertzsprung SONG Telescope: Asteroseismology of the G5 Subgiant Star μ Herculis

We report the first asteroseismic results obtained with the Hertzsprung Stellar Observations Network Group Telescope from an extensive high-precision radial-velocity observing campaign of the subgiant μ Herculis. The data set was collected during 215 nights in 2014 and 2015. We detected a total of 49 oscillation modes with l values from zero to three, including some l = 1 mixed modes. Based on the rotational splitting observed in l = 1 modes, we determine a rotational period of 52 days and a stellar inclination angle of 63°. The parameters obtained through modeling of the observed oscillation frequencies agree very well with independent observations and imply a stellar mass between 1.11 and 1.15 M {sub ⊙} and an age of 7.8{sub −0.4}{sup +0.3} Gyr. Furthermore, the high-quality data allowed us to determine the acoustic depths of the He ii ionization layer and the base of the convection zone.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Identification and Photometric Classification of Extragalactic Transients in the Vera C. Rubin Observatory’s Data Preview 1

The Vera C. Rubin Observatory will soon survey the southern sky, delivering a depth and sky coverage that is unprecedented in time-domain astronomy. As part of commissioning, Data Preview 1 (DP1) has been released. It comprises a Legacy Survey of Space and Time (LSST) Commissioning Camera observing campaign between 2024 November and December with multiband imaging of seven fields, covering roughly 0.4 deg 2 each, providing a first glimpse into the data products that will become available once the LSST begins. In this work, we search three fields for extragalactic transients. We identify eight new likely supernovae (SNe), and three known ones from a sample of 369,644 difference image analysis objects. Photometric classification using Superphot+ assigns subclasses with >95% confidence to only one SN Ia and one SN II in this sample. Our findings are in agreement with SN detection rate predictions of 15 ± 4 SNe from simulations using simsurvey. The SN detection rate in the data is possibly affected by the lack of suitable templates. Nevertheless, this work demonstrates the quality of the data products delivered in DP1 and indicates that the Rubin Observatory’s LSST is well placed to fulfill its discovery potential in time-domain astronomy.

Freeburn, James [University of North Carolina, Cha↗
DOE OSTI2023

NETL Plastic Pipes Project (Final Report)

Plastic or composite pipelines have been the bane of the utility locating industry because they are neither conductive nor magnetic which are the properties traditionally used to locate buried utilities. Ground penetrating radar (GPR) is an effective geophysical tool for locating plastic/composite pipelines where resistive cover allows for adequate penetration of radar energy. However, GPR has limited applicability in areas where the soil cover is conductive due to significant clay and/or salt content. This study examines complementary near-surface geophysical methods that are potentially useful for locating buried plastic/composite pipelines, either singly or in combination. Specifically, this modeling study used computational numerical methods to forward model the response of GPR, resistivity, seismic, gravity gradiometry, and photoacoustic/thermoacoustic imaging methods to plastic/composite pipelines for various scenarios including: (1) pipe diameters ranging between 2 in. to 12 in.; (2) burial depths ranging between 3 ft. to 4 ft.; (3) various degrees in contrast in physical properties (i.e., electrical permittivity, elasticity, resistivity, density); and (4) various experimental acquisition choices (e.g., GPR radar and seismic source frequencies, electrode spacing). Numerical modeling performed herein reconfirmed that GPR is the preferred method for detecting/locating plastic pipelines. A caveat for GPR detection is that the material covering the plastic pipe (trench fill material and adjacent soil) must be sufficiently resistive to allow the two-way propagation to the required depth of investigation and back to the surface. GPR was the only method modeled in this study that can be used to directly detect plastic pipelines of 2-in.-diameter and larger when buried 3-ft-deep. GPR data processing and imaging also can determine pipe depth, pipe diameter, trench dimensions, and moisture conditions. Seismic modeling results suggest that direct detection of a 12-in.-diameter plastic pipe at 3-ft.-depth may be possible under favorable conditions; however, the associated signature would be weak (e.g., surface- to S-wave, backscattered surface-waves, and/or forward scattered surface-waves to S-wave). Direct pipe detection under field conditions with noise and strong lateral geologic heterogeneity is doubtful. Numerical modeling also suggests that plastic pipelines can be indirectly located by detecting the trench in which they are buried. GPR, direct current (DC) resistivity, and seismic methods have the potential to locate the pipeline trench if there is sufficient contrast between the trench-wall and trench-fill materials for the physical property being measured by each method (i.e., electrical permittivity for GPR; resistivity for DC resistivity; or density, compressional velocity, or shear velocity for seismic). Modeling also indicated that currently available (commercial) gravity gradiometers would be unable to directly detect/locate plastic pipelines ≤ 8-in.-diameter when buried 3-ft.-deep given the typical instrument noise floor for field surveying as well as the expected density variations due to geologic heterogeneity. The numerical modeling performed in this project did not identify a universal geophysical technology that can locate buried plastic pipelines in all parts of the United States (although GPR is suggested for all areas with resistive cover). However, the project results suggest that a towed land streamer simultaneously acquiring multiple geophysical data types including multi-offset GPR, multi-channel DC resistivity, seismic geophone- and/or distributed acoustic sensing (DAS), and potentially photoacoustic/thermoacoustic data would be an appropriate platform for locating buried plastic pipeline. Moreover, the complementary multiphysics data acquired by a towed land streamer would permit the use of joint and/or cooperative inversion frameworks for a more rigorous and consistent data interpretation.

42 ENGINEERING↗
DOE OSTI2015

A patient-specific aperture system with an energy absorber for spot scanning proton beams: Verification for clinical application

Purpose: In the authors’ proton therapy system, the patient-specific aperture can be attached to the nozzle of spot scanning beams to shape an irradiation field and reduce lateral fall-off. The authors herein verified this system for clinical application. Methods: The authors prepared four types of patient-specific aperture systems equipped with an energy absorber to irradiate shallow regions less than 4 g/cm{sup 2}. The aperture was made of 3-cm-thick brass and the maximum water equivalent penetration to be used with this system was estimated to be 15 g/cm{sup 2}. The authors measured in-air lateral profiles at the isocenter plane and integral depth doses with the energy absorber. All input data were obtained by the Monte Carlo calculation, and its parameters were tuned to reproduce measurements. The fluence of single spots in water was modeled as a triple Gaussian function and the dose distribution was calculated using a fluence dose model. The authors compared in-air and in-water lateral profiles and depth doses between calculations and measurements for various apertures of square, half, and U-shaped fields. The absolute doses and dose distributions with the aperture were then validated by patient-specific quality assurance. Measured data were obtained by various chambers and a 2D ion chamber detector array. Results: The patient-specific aperture reduced the penumbra from 30% to 70%, for example, from 34.0 to 23.6 mm and 18.8 to 5.6 mm. The calculated field width for square-shaped apertures agreed with measurements within 1 mm. Regarding patient-specific aperture plans, calculated and measured doses agreed within −0.06% ± 0.63% (mean ± SD) and 97.1% points passed the 2%-dose/2 mm-distance criteria of the γ-index on average. Conclusions: The patient-specific aperture system improved dose distributions, particularly in shallow-region plans.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

The effects of burial diagenesis on multiscale porosity in the St. Peter Sandstone: An imaging, small-angle, and ultra-small-angle neutron scattering analysis

To examine the effects of burial diagenesis on heirarchical pore structures in sandstone and compare those with the effects of overgrowth formation, we obtained samples of St. Peter Sandstone from drill cores obtained in the Illinois and Michigan Basins. The multiscale pore structure of rocks in sedimentary reservoirs and the mineralogy associated with those pores are critical factors for estimating reservoir properties, including fluid mass in place, permeability, and capillary pressures, as well as geochemical interactions between the rock and the fluid. The combination of small- and ultra-small-angle neutron scattering with backscattered electron or X ray-computed tomographic imaging, or both, provided a means by which pore structures were quantified at scales ranging from aproximately 1 nm to 1 cm—seven orders of magnitude. Larger scale (>10 µm) porosity showed the expected logarithmic decrease in porosity with depth, although there was significant variation in each sample group. However, small- and ultra-small-angle neutron scattering data showed that the proportion of small-scale porosity increased with depth. Porosity distributions were not continuous, but consisted of a series of log normal-like distributions at several distinct scales within these rocks. Fractal dimensions at larger scales decreased (surfaces smoothed) with increasing depth, and those at smaller scales increased (surfaces roughened) and pores become more isolated (higher lacunarity). Furthermore, data suggest that changes in pore-size distributions are controlled by both physical (compaction) and chemical effects (precipitation, cementation, dissolution).

58 GEOSCIENCES↗
DOE OSTI2026

Integration of ultra-low coverage whole-genome sequences for reconstructing the evolutionary history of Galapagos giant tortoises

Genomic data from contemporary and historical samples often need to be coupled for evolutionary reconstructions of multitaxon complexes. However, the genetic data recovered from historical samples may result only in ultra-low coverage whole-genome sequences (ulcWGS; <0.15× depth), leading to inaccurate evolutionary inferences given a preponderance of missing data. Using the Galapagos giant tortoise radiation as a study system (Chelonoidis spp., composed of 13 extant and four extinct lineages), we assembled a novel methodological pipeline that removes potential noise introduced by the missing data and enhances the evolutionary signal from ulcWGS samples. We leveraged existing tools for phylogenomic placement (EPA-ng), population genomic structure (smartsnp) and admixture (Admixfrog, NGSadmix) to demonstrate that the evolutionary history of samples can be uncovered with sequencing depths as low as 0.008–0.139×. Importantly, these approaches do not use genotype imputation of the ulcWGS samples, which would require extensive reference datasets. Our application to two cases of extinct lineages of Galapagos giant tortoises, with and without references from the same lineage, demonstrates the general value of the approach. We confirm where the extinct lineages from San Cristóbal and Santa Fe islands fit into the Galapagos giant tortoise radiation, and that these lineages were evolutionarily distinct entities.

ancient DNA↗
DOE OSTI2014

TU-A-9A-10: Verification of Photoacoustic Computed Tomography Perfusion Imaging Using DCE-CT

Purpose: We propose to verify quantifiable perfusion information generated by a Photoacoustic Computed Tomography (PCT) scanner using Dynamic Contrast-Enhanced CT (DCE-CT), and to investigate physicsbased models of acoustic properties of tissue and photon transport to improve quantification. These corrections first necessitate a skin identifying algorithms to reduce speed-of-sound blurring and empirical photon correction methods. Methods: Xenograft mice (n=5) of breast cancer was imaged using DCE-CT which was followed by DCE-PCT. To obtain CT perfusion data, each mouse was i.v. injected (0.2mL Isovue @0.5mL/min) and subsequent radio-opaque time curves fit to a 2-compartmental model on a voxel-wise basis. For DCE-PCT, different concentrations of ICG (250, 125, and 62.5 micro-Molar) were injected at the same rate, but also acquired at different sampling rates (3, 6, and 12 seconds). The time intensity curves from PCT were fit to a 1-compartmental model on a voxel by voxel basis. The images were coregistered (Oncentra) based on the structural similarities of the tumor vasculature after which we compared both the contrastenhanced dynamics and the vascular physiology. Results: Moderate to high doses of ICG impact the washin phase of the PCT contrast due to photon losses as a function of depth. A semi-automatic algorithm has been developed to identify the skin margin, and subsequent MC and empirical models of photon transport and variations in speed-of-sound are being evaluated. Conclusion: From our results we find that there is a need to apply photon and speed-of-sound corrections to our PCT data to improve the quantifiable image data at depth in the tumor for PCT. The dose and injection rate may help in reducing large systematic effects. Our project is partially funded by a NIH SBIR grant.

60 APPLIED LIFE SCIENCES↗