Search NASA⌕ Search

SEARCH · Search NASA

Results for “difference map”

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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Indicator Patterns of Forced Change Learned by an Artificial Neural Network

Abstract Many problems in climate science require the identification of signals obscured by both the “noise” of internal climate variability and differences across models. Following previous work, we train an artificial neural network (ANN) to predict the year of a given map of annual‐mean temperature (or precipitation) from forced climate model simulations. This prediction task requires the ANN to learn forced patterns of change amidst a background of climate noise and model differences. We then apply a neural network visualization technique (layerwise relevance propagation) to visualize the spatial patterns that lead the ANN to successfully predict the year. These spatial patterns thus serve as “reliable indicators” of the forced change. The architecture of the ANN is chosen such that these indicators vary in time, thus capturing the evolving nature of regional signals of change. Results are compared to those of more standard approaches like signal‐to‐noise ratios and multilinear regression in order to gain intuition about the reliable indicators identified by the ANN. We then apply an additional visualization tool (backward optimization) to highlight where disagreements in simulated and observed patterns of change are most important for the prediction of the year. This work demonstrates that ANNs and their visualization tools make a powerful pair for extracting climate patterns of forced change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

R and DE Robotic Sensor Intern Team: Remote 3D Mapping Overview

The Department of Energy (DOE) nuclear complex includes many hazardous and challenging environments that must be inspected remotely and routinely for structural integrity evaluation and long-term planning purposes. The H-Canyon Air Exhaust (HCAEX) Tunnel biannual inspection is an example of such an environment. Presently inspections are performed by teleoperated robotic crawlers provide low resolution video feedback to the structural engineers. Advanced sensors such as Lidar and panospheric hi-res cameras have the potential to provide additional valuable data in future inspections. UT-Austin Lidar work in support of an H-Canyon deployment began with the DOE EM sponsored H-Canyon Integrated Research Project. The R and DE Robotic Sensor intern team includes two of the UT-Austin Masters students working on Lidar development to facilitate the pulling of university developed solutions to the field. The goal of this project is to develop a suite of sensors that could be used to gather advanced data in future HCAEX tunnel inspections. Of particular interest is the ability to map the tunnel walls at an accuracy that will allow for the comparison of two scans performed at different time periods to estimate the concrete erosion rate. Additionally, the development of the 3D point cloud in the GPS denied environment is the first step in being able to collect and overlay other 3D data such as radiation mapping. Project specific goals include: - Design and develop 3D mapping sensor suite. Merge Lidar point cloud with high res camera images Incorporate Inertial Measurement Unit (IMU) to assist with positioning. - Develop sensor software and user GUI Sensor data collection and integration software User GUI for reach back system. - Test solution on crawler in tunnel test bed Use crawler base similar to H-Canyon planned crawler Design and build test bed simulating tunnel conditions. To facilitate rapid sensor integration, software coding and user GUI development, the Robot Operating System (ROS) was employed as the development platform. ROS is a set of open source software libraries and tools that assist with the building of robot applications. ROS includes drivers for many sensors, state-of-the-art algorithms, and powerful developer tools. Sensors selected for the summer intern solution included assessment as to whether they possessed ROS drivers and associated sample code. Additionally, the ROS 3D Robot Visualizer (rviz) tool has been employed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1984

Selection of a map grid for data analysis and archival

Arguments for selection of a map grid are reiterated by illustrating the quantitative effects on data quality caused by using different grids. It is shown that the use of the rectangular latitude-longitude grid actually degrades data quality, increases stored data volume, and increases the complexity of data manipulation. Use of an equal-area grid is shown to be a proper way to aggregate data, but such grids are thought to be inconvenient. A solution to this dilemma is proposed by showing that the analysis and archival map grids need not be the same. The results of a proper analysis of the data on an equal-area grid can be remapped to the more 'convenient' rectangular latitude-longitude grid without loss of quality.

Rossow, W. B.↗
NASA NTRS1980

Application of satellite data for snow mapping in Norway

The author has identified the following significant results. A close quantitative relationship was found between snow covered areas and subsequent runoff for different parts of the country despite climate differences. Digital LANDSAT data can be used for areas down to approximately 10 sq km to 20 sq km for accurate measurement of snow cover extent. On large watersheds (more than 500 sq km), digital NOAA/TIROS imagery can be used for snow mapping if the area/runoff relationship is determined by using observations from previous years.

Odegaard, H. A.↗
DOE OSTI2020

A Multi-Sensor Unoccupied Aerial System Improves Characterization of Vegetation Composition and Canopy Properties in the Arctic Tundra

Changes in vegetation distribution, structure, and function can modify the canopy properties of terrestrial ecosystems, with potential consequences for regional and global climate feedbacks. In the Arctic, climate is warming twice as fast as compared to the global average (known as ‘Arctic amplification’), likely having stronger impacts on arctic tundra vegetation. In order to quantify these changes and assess their impacts on ecosystem structure and function, methods are needed to accurately characterize the canopy properties of tundra vegetation types. However, commonly used ground-based measurements are limited in spatial and temporal coverage, and differentiating low-lying tundra plant species is challenging with coarse-resolution satellite remote sensing. The collection and processing of multi-sensor data from unoccupied aerial systems (UASs) has the potential to fill the gap between ground-based and satellite observations. To address the critical need for such data in the Arctic, we developed a cost-effective multi-sensor UAS (the ‘Osprey’) using off-the-shelf instrumentation. The Osprey simultaneously produces high-resolution optical, thermal, and structural images, as well as collecting point-based hyperspectral measurements, over vegetation canopies. In this paper, we describe the setup and deployment of the Osprey system in the Arctic to a tundra study site located in the Seward Peninsula, Alaska. We present a case study demonstrating the processing and application of Osprey data products for characterizing the key biophysical properties of tundra vegetation canopies. In this study, plant functional types (PFTs) representative of arctic tundra ecosystems were mapped with an overall accuracy of 87.4%. The Osprey image products identified significant differences in canopy-scale greenness, canopy height, and surface temperature among PFTs, with deciduous low to tall shrubs having the lowest canopy temperatures while non-vascular lichens had the warmest. The analysis of our hyperspectral data showed that variation in the fractional cover of deciduous low to tall shrubs was effectively characterized by Osprey reflectance measurements across the range of visible to near-infrared wavelengths. Therefore, the development and deployment of the Osprey UAS, as a state-of-the-art methodology, has the potential to be widely used for characterizing tundra vegetation composition and canopy properties to improve our understanding of ecosystem dynamics in the Arctic, and to address scale issues between ground-based and airborne/satellite observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Heat Generation versus Heat Dissipation Maps for Batteries (Battery Safety Maps) [SWR-21-27]

NREL's Battery Safety Map (Heat Generation versus Heat Dissipation Maps for Batteries) generates a response surface that is grouped into regions where the battery is deemed safe from the propagation of thermal runaway, those where an engineering solution would render the battery safe, and finally instances where thermal runaway is inevitable. The demarcation for the different regions is provided by specifications or constraints imposed by the end user. NREL's Battery Safety Map is a user-defined subroutine that sweeps through a set of parameters describing the operating conditions and short-circuit scenarios for battery safety simulations. For each run, the electrochemical-thermal model implementing reaction kinetics and gas generation are run to evaluate the amount of heat generated. A thermal model then captures the amount of heat dissipation based on boundary conditions defined (e.g., forced convection, constant temperature at the surface, etc.). Based on the rate of heat transfer and the amount of heat generated, each run is classified into "safe," "mitigation possible," or "unsafe" zones. These categories are defined based on constraints or tolerance specifications provided by the end user.

Santhanagopalan, Shriram↗
NASA NTRS2023

Predicting Patterns of Solar Energy Buildout to Identify Opportunities for Biodiversity Conservation

The construction of solar energy facilities can have positive or negative impacts on biodiversity depending on siting and associated land use transitions. We identified drivers of solar siting and quantified patterns of buildout in states surrounding the Chesapeake Bay watershed – a biodiversity hotspot with numerous ecosystem services. Using a convolutional neural network, we mapped the footprints of ground-mounted solar arrays present in satellite imagery annually from 2017 to 2021 in Delaware, Maryland, Pennsylvania, New York, Virginia, and West Virginia. As of 2021, we identified 958 solar arrays covering 52.3 km2 built primarily on previously cultivated land, while avoiding natural landcover. We fit a binomial-Weibull model to these solar timeseries data in a hierarchical, Bayesian framework to quantify the relationship between geospatial covariates and rate of solar development. Solar array construction rate increased in cultivated areas, areas of lower agricultural suitability, lower slope, lower forest cover, lower biodiversity protection, and greater distances from roads. We also estimated changes in the rate of solar construction over time and found differences among states: acceleration in Virginia and deceleration in New York. We used parameter estimates to map the relative likelihood of future solar development across the study area. This methodology can be used to anticipate where solar is likely to be built in different landscapes and how these patterns align with conservation goals. Around the Chesapeake Bay watershed, the selection of lower quality agricultural areas for solar energy minimizes removal of important habitat and provides opportunities for native plant and pollinator restoration.

Artificial Intelligence↗
NASA NTRS1975

Utilization of satellite data for inventorying prairie ponds and lakes. LANDSAT-1 data were used to discriminate ponds and lakes for waterfowl management

The author has identified the following significant results. The mapping of open water as an indicator of waterfowl habitat quality was carried out by using two different recognition techniques, a single waveband thresholding approach and a multiple waveband approach termed proportion estimation. The single waveband technique has proven simple to implement. Its computer algorithm was rapid and accurately recognized prairie lakes and large ponds. The resultant products of this processing technique were thematic maps and statistical tabulations describing open surface water conditions. The maps served to portray visually the location and frequency of surface water bodies but usually necessitated additional interpretation.

Gilmer, D. S.↗
DOE OSTI2019

Topological string entanglement

We investigate how topological entanglement of Chern-Simons theory is captured in a string theoretic realization. Our explorations are motivated by a desire to understand how quantum entanglement of low energy open string degrees of freedom is encoded in string theory (beyond the oft discussed classical gravity limit). Concretely, we realize the Chern-Simons theory as the worldvolume dynamics of topological D-branes in the topological A-model string theory on a Calabi-Yau target. Via the open/closed topological string duality one can map this theory onto a pure closed topological A-model string on a different target space, one which is related to the original Calabi-Yau geometry by a geometric/conifold transition. We demonstrate how to uplift the replica construction of Chern-Simons theory directly onto the closed string and show that it provides a meaningful definition of reduced density matrices in topological string theory. Furthermore, we argue that the replica construction commutes with the geometric transition, thereby providing an explicit closed string dual for computing reduced states, and Rényi and von Neumann entropies thereof. While most of our analysis is carried out for Chern-Simons on $\hbox{S}^3$, the emergent picture is rather general. Specifically, we argue that quantum entanglement on the open string side is mapped onto quantum entanglement on the closed string side and briefly comment on the implications of our result for physical holographic theories where entanglement has been argued to be crucial ingredient for the emergence of classical geometry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Characterization of Rare Earth Element Minerals in Coal Utilization Byproducts

The United States currently produces over 100 million tons of coal utilization byproducts (CUB) per year in the form of fly ash, bottom ash, slag, and flue gas (American Coal Ash Association (ACCA), 2015). But this “waste material” also contains potentially useful levels of rare earth elements (REE). Rare earth elements are crucial for many existing and emerging technologies, but the U.S. lacks a domestic, sustainable REE source. Our project explored the possibility of developing a supply of REEs for U.S. technologies by extracting REEs from CUBs. This work offers the potential to reduce our dependence on other countries for supply of these critical elements (NETL, REE 2016 Project Portfolio). Geologic and diagenetic history, industrial preparation methods, and the specific combustion process all play major roles in the composition of CUB. During combustion, inorganic mineral phases of coal particles are fluidized at temperatures higher than 1400oC, so inorganic mineral materials are oxidized, fused, disintegrated, or agglomerated into larger spherical and amorphous (non-crystalline) particles. The original mineralogy of the coal-containing rock and heating/cooling of the material significantly affects the composition and morphology of the particles in the combustion byproduct (Kutchko and Kim, 2006). Thus, different types of coal/refuse/ash must be characterized to better understand mineral evolution during the combustion process. Our research focused on developing a working model to address how REE minerals behave during the combustion process: this research should help determine the most effective engineering methods for extracting REEs from CUBs. We used multimodal imaging and image processing techniques to characterize six rock and ash samples from different coal power plants with respect to morphology, grain size, presence of mineral phases, and elemental composition. The results of these characterization activities provided thresholds for realizing the occurrence of REE mineral phases in CUB and allowed us to calculate structural and volumetric estimates of REE. Collectively, the rock and coal ash samples contained minerals such as quartz, kaolinite, muscovite/illite, iron oxide (as hematite or magnetite), mullite, and clinochlore. Trace minerals included pyrite, zircon, siderite, rutile, diopside, foresterite, gypsum, and barite. We identified REE phosphate minerals monazite (Ce,La,Nd,Th)(PO 4 ,SiO 4 ), xenotime (YPO 4 ,SiO 4 ), and apatite (Ca 5 (PO 4 ) 3 (F,Cl,OH) via SEM and electron microprobe analysis: these materials generally occurred as 1-10 μm-long crystals in the rock and ash samples. As has been shown in other studies, amorphous material-aluminosilicate glass or iron oxyhydroxide-are the major components of coal fly and bottom ash. Trace amounts of amorphous calcium oxide and mixed element (e.g., Al-Si-Ca-Fe) slag are also present. Quartz, mullite, hematite, and magnetite are the crystalline phases present. We found that REEs are present as monomineralic grains dispersed within the ash, as well as fused to or encapsulated by amorphous aluminosilicate glass particles. Monazite and xenotime have relatively high melting points (>1800 °C) compared to typical combustion temperatures; our observations indicate that the REE-phosphates, which presumably contribute a large percentage of REE to the bulk ash REE pool, as measured by mass spectroscopy, are largely unaltered by the combustion. Our study shows that conventional coal combustion processes sequester REE minerals into aluminosilicate glass phases, which presents a new engineering challenge for extracting REE from coal ash. The characterization work summarized in this report provides a semi-quantitative assessments of REE in coal-containing rock and CUB. The data we obtained from 2- and 3-D imaging, elemental mapping, volumetric estimates, and advanced high-resolution pixel classification successfully identified the different mineral phases present in CUB. Further, our characterization results can guide techniques for extracting REEs from CUB, or other geologic and engineered materials. Whilst, interpretations will inform future REE separation and extraction techniques and technologies practical for commercial utilization of combustion byproducts generated by power plants.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2018

Optimizing de novo genome assembly from PCR-amplified metagenomes

Background. Metagenomics has transformed our understanding of microbial diversity across ecosystems, with recent advances enabling de novo assembly of genomes from metagenomes. These metagenome-assembled genomes are critical to provide ecological, evolutionary, and metabolic context for all the microbes and viruses yet to be cultivated. Metagenomes can now be generated from nanogram to subnanogram amounts of DNA. However, these libraries require several rounds of PCR amplification before sequencing, and recent data suggest these typically yield smaller and more fragmented assemblies than regular metagenomes. Methods. Here we evaluate de novo assembly methods of 169 PCR-amplified metagenomes, including 25 for which an unamplified counterpart is available, to optimize specific assembly approaches for PCR-amplified libraries. We first evaluated coverage bias by mapping reads from PCR-amplified metagenomes onto reference contigs obtained from unamplified metagenomes of the same samples. Then, we compared different assembly pipelines in terms of assembly size (number of bp in contigs ≥ 10kb) and error rates to evaluate which are the best suited for PCR-amplified metagenomes. Results. Read mapping analyses revealed that the depth of coverage within individual genomes is significantly more uneven in PCR-amplified datasets versus unamplified metagenomes, with regions of high depth of coverage enriched in short inserts. This enrichment scales with the number of PCR cycles performed, and is presumably due to preferential amplification of short inserts. Standard assembly pipelines are confounded by this type of coverage unevenness, so we evaluated other assembly options to mitigate these issues. We found that a pipeline combining read deduplication and an assembly algorithm originally designed to recover genomes from libraries generated after whole genome amplification (single-cell SPAdes) frequently improved assembly of contigs ≥ 10kb by 10 to 100-fold for low input metagenomes. Conclusions. PCR-amplified metagenomes have enabled scientists to explore communities traditionally challenging to describe, including some with extremely low biomass or from which DNA is particularly difficult to extract. Here we show that a modified assembly pipeline can lead to an improved de novo genome assembly from PCR-amplified datasets, and enables a better genome recovery from low input metagenomes.

Roux, Simon↗
DOE OSTI2018

Optimizing de novo genome assembly from PCR-amplified metagenomes

Background: Metagenomics has transformed our understanding of microbial diversity across ecosystems, with recent advances enabling de novo assembly of genomes from metagenomes. These metagenome-assembled genomes are critical to provide ecological, evolutionary, and metabolic context for all the microbes and viruses yet to be cultivated. Metagenomes can now be generated from nanogram to subnanogram amounts of DNA. However, these libraries require several rounds of PCR amplification before sequencing, and recent data suggest these typically yield smaller and more fragmented assemblies than regular metagenomes. Methods: Here we evaluate de novo assembly methods of 169 PCR-amplified metagenomes, including 25 for which an unamplified counterpart is available, to optimize specific assembly approaches for PCR-amplified libraries. We first evaluated coverage bias by mapping reads from PCR-amplified metagenomes onto reference contigs obtained from unamplified metagenomes of the same samples. Then, we compared different assembly pipelines in terms of assembly size (number of bp in contigs ≥ 10kb) and error rates to evaluate which are the best suited for PCR-amplified metagenomes. Results: . Read mapping analyses revealed that the depth of coverage within individual genomes is significantly more uneven in PCR-amplified datasets versus unamplified metagenomes, with regions of high depth of coverage enriched in short inserts. This enrichment scales with the number of PCR cycles performed, and is presumably due to preferential amplification of short inserts. Standard assembly pipelines are confounded by this type of coverage unevenness, so we evaluated other assembly options to mitigate these issues. We found that a pipeline combining read deduplication and an assembly algorithm originally designed to recover genomes from libraries generated after whole genome amplification (single-cell SPAdes) frequently improved assembly of contigs ≥ 10kb by 10 to 100-fold for low input metagenomes. Conclusions: PCR-amplified metagenomes have enabled scientists to explore communities traditionally challenging to describe, including some with extremely low biomass or from which DNA is particularly difficult to extract. Here we show that a modified assembly pipeline can lead to an improved de novo genome assembly from PCR-amplified datasets, and enables a better genome recovery from low input metagenomes.

Roux, Simon↗
DOE OSTI2026

The Q 10 of in situ microbial soil respiration varies with mean annual temperature, precipitation, pH, and plant cover: a meta-analysis and spatial prediction of Q 10

The temperature sensitivity of soil microbial respiration, commonly quantified using the Q 10 coefficient, is a key parameter in carbon cycle models. Uncovering how environmental factors affect in situ Q 10 values can therefore provide critical insight into potential shifts in global carbon stocks under climate change. We collected data from previously published field experiments that measured soil microbial respiration across a range of temperatures. We hypothesized that the Q 10 coefficient of in situ soil microbial respiration would vary based on environmental factors including mean annual temperature (MAT), mean annual precipitation (MAP), plant cover type, pH, soil C:N, and latitude. Linear regression revealed that Q 10 correlates negatively with MAT and MAP and positively with pH and absolute latitude. Additionally, average Q 10 varied significantly across different plant cover types; it was highest in mountain grasslands and lowest in tropical moist forests. Variation in microbial Q 10 across environmental factors may arise from underlying mechanisms such as enzyme kinetics, substrate availability and complexity, and microbial adaptation. To capture patterns in Q 10 more comprehensively, we developed a multiple linear regression model of Q 10 based on the most individually significant environmental drivers and applied it to public datasets to generate a global map of predicted Q 10 . Q 10 was higher in high-latitude and high-altitude regions, where large permafrost carbon stores are vulnerable to thawing and decomposition. We also compared fits between the Q 10 equation and a model produced from macromolecular rate theory (MMRT). We found that the MMRT model had the superior fit and may be better suited to model temperature sensitivity of complex biological reactions. Overall, our results emphasize that relationships between microbial Q 10 and environmental variables should be accounted for in climate models. Incorporating these variations in the Q 10 parameter, rather than using a fixed value, will help predict whether CO 2 emissions will be buffered or exacerbated by soil microbial respiration under climate change.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2001

Seismic Study of the Dynamics of the Solar Subsurface from SoHO Observations

In collaboration with Dr. Baudin, we have developed and refined the new observational methodology for local helioseismology known as time-distance analysis. Global helioseismology study the solar oscillations as a superposition of resonant modes, whose properties (mode frequencies) reflect the global structure of the sun (sound speed stratification, rotation rate, etc). In contrast, local helioseismology look at the solar oscillations as wave packets whose propagation will be affected by perturbations of the media sampled. These local perturbations (sound speed or velocity flows) will modify the propagation time, that in turn can be used as a diagnostic tool for a given region. From a data reduction perspective, the processing of solar dopplergrams that result in time-distance maps, i.e. propagation times as a function of distance between bounces at the surface, is radically different from the methodology used for global mode analysis. We have, in a first step, further develop the programs needed to carry out such analysis. We have then applied them to NMI data set, and explore the trade-off between various averaging and filtering approaches - steps required to improve the signal-to-noise ratio of correlation maps - and the resulting stability and precision of the fitted propagation times. While excessive averaging (whether over space, propagation distance, or time) will reduce the diagnostic potential of the method, insufficient averaging lead to unstable fits, or uncertainties so large as to hide the information we seek. In a second phase, we have developed the analysis methodology required to infer local properties from perturbation in time propagation. Namely, we have developed time-distance inversion techniques, with an emphasis on inferences of velocity flows from time anomalies. Note also that during the period covered by this grant, all the investigators on this proposal (i.e., Drs. Baudin, Eff-Darwich, Korzennik, and Noyes) took part in the organization of the SOHO 6 /GONG 99 Workshop: Structure and Dynamics of the Interior of the Sun and Sunlike Stars, held on June 1-4 1999 at the Boston Park Plaza Hotel in Boston, Massachusetts, USA. it was very well attended by more than 160 participants from 26 countries from all over the world. The proceedings were published in two volumes as ESA SP-418, with Sessions I-III in Volume 1, and Sessions IV-VI in Volume 2 (1,000 pages in total). The complete contents are also included in digital form on a CD-ROM included with Volume 1. This CD-ROM also contains additional multi-media material that complements some of the contributions.

Korzennik, Sylvain G.↗
NASA NTRS1993

A unified approach to VLSI layout automation and algorithm mapping on processor arrays

Development of software tools for designing supercomputing systems is highly complex and cost ineffective. To tackle this a special purpose PAcube silicon compiler which integrates different design levels from cell to processor arrays has been proposed. As a part of this, we present in this paper a novel methodology which unifies the problems of Layout Automation and Algorithm Mapping.

Venkateswaran, N.↗
DOE OSTI2017

Deterministic direct aperture optimization using multiphase piecewise constant segmentation

Purpose Direct aperture optimization (DAO) attempts to incorporate machine constraints in the inverse optimization to eliminate the post‐processing steps in fluence map optimization (FMO) that degrade plan quality. Current commercial DAO methods utilize a stochastic or greedy approach to search a small aperture solution space. In this study, we propose a novel deterministic direct aperture optimization that integrates the segmentation of fluence map in the optimization problem using the multiphase piecewise constant Mumford‐Shah formulation. Methods The Mumford‐Shah based direct aperture optimization problem was formulated to include an L2‐norm dose fidelity term to penalize differences between the projected dose and the prescribed dose, an anisotropic total variation term to promote piecewise continuity in the fluence maps, and the multiphase piecewise constant Mumford‐Shah function to partition the fluence into pairwise discrete segments. A proximal‐class, first‐order primal‐dual solver was implemented to solve the large scale optimization problem, and an alternating module strategy was implemented to update fluence and delivery segments. Three patients of varying complexity—one glioblastoma multiforme (GBM) patient, one lung (LNG) patient, and one bilateral head and neck (H&N) patient with 3 PTVs—were selected to test the new DAO method. For each patient, 20 non‐coplanar beams were first selected using column generation, followed by the Mumford‐Shah based DAO (DAO MS ). For comparison, a popular and successful approach to DAO known as simulated annealing—a stochastic approach—was replicated. The simulated annealing DAO (DAO SA ) plans were then created using the same beam angles and maximum number of segments per beam. PTV coverage, PTV homogeneity , and OAR sparing were assessed for each plan. In addition, high dose spillage, defined as the 50% isodose volume divided by the tumor volume, as well as conformity, defined as the van't Riet conformation number, were evaluated. Results DAO MS achieved essentially the same OAR doses compared with the DAO SA plans for the GBM case. The average difference of OAR D max and D mean between the two plans were within 0.05% of the plan prescription dose. The lung case showed slightly improved critical structure sparing using the DAO MS approach, where the average OAR D max and D mean were reduced by 3.67% and 1.08%, respectively, of the prescription dose. The DAO MS plan substantially improved OAR dose sparing for the H&N patient, where the average OAR D max and D mean were reduced by over 10% of the prescription dose. The DAO MS and DAO SA plans were comparable for the GBM and LNG PTV coverage, while the DAO MS plan substantially improved the H&N PTV coverage, increasing D99 by 6.98% of the prescription dose. For the GBM and LNG patients, the DAO MS and DAO SA plans had comparable high dose spillage but slightly worse conformity with the DAO MS approach. For the H&N plan, DAO MS was considerably superior in high dose spillage and conformity to the DAO SA . The deterministic approach is able to solve the DAO problem substantially faster than the simulated annealing approach, with a 9.5‐ to 40‐fold decrease in total solve time, depending on the patient case. Conclusions A novel deterministic direct aperture optimization formulation was developed and evaluated. It combines fluence map optimization and the multiphase piecewise constant Mumford‐Shah segmentation into a unified framework, and the resulting optimization problem can be solved efficiently. Compared to the widely and commercially used simulated annealing DAO approach, it showed comparable dosimetry behavior for simple plans, and substantially improved OAR sparing, PTV coverage, PTV homogeneity, high dose spillage, and conformity for the more complex head and neck plan.

Nguyen, Dan↗
DOE OSTI2019

Optimizing de novo genome assembly from PCR-amplified metagenomes

Background. Metagenomics has transformed our understanding of microbial diversity across ecosystems, with recent advances enabling de novo assembly of genomes from metagenomes. These metagenome-assembled genomes are critical to provide ecological, evolutionary, and metabolic context for all the microbes and viruses yet to be cultivated. Metagenomes can now be generated from nanogram to subnanogram amounts of DNA. However, these libraries require several rounds of PCR amplification before sequencing, and recent data suggest these typically yield smaller and more fragmented assemblies than regular metagenomes. Methods. Here we evaluate de novo assembly methods of 169 PCR-amplified metagenomes, including 25 for which an unamplified counterpart is available, to optimize specific assembly approaches for PCR-amplified libraries. We first evaluated coverage bias by mapping reads from PCR-amplified metagenomes onto reference contigs obtained from unamplified metagenomes of the same samples. Then, we compared different assembly pipelines in terms of assembly size (number of bp in contigs ≥ 10 kb) and error rates to evaluate which are the best suited for PCR-amplified metagenomes. Results. Read mapping analyses revealed that the depth of coverage within individual genomes is significantly more uneven in PCR-amplified datasets versus unamplified metagenomes, with regions of high depth of coverage enriched in short inserts. This enrichment scales with the number of PCR cycles performed, and is presumably due to preferential amplification of short inserts. Standard assembly pipelines are confounded by this type of coverage unevenness, so we evaluated other assembly options to mitigate these issues. We found that a pipeline combining read deduplication and an assembly algorithm originally designed to recover genomes from libraries generated after whole genome amplification (single-cell SPAdes) frequently improved assembly of contigs ≥10 kb by 10 to 100-fold for low input metagenomes. Conclusions. PCR-amplified metagenomes have enabled scientists to explore communities traditionally challenging to describe, including some with extremely low biomass or from which DNA is particularly difficult to extract. Here we show that a modified assembly pipeline can lead to an improved de novo genome assembly from PCR-amplified datasets, and enables a better genome recovery from low input metagenomes.

59 BASIC BIOLOGICAL SCIENCES↗