Search NASA⌕ Search

SEARCH · Search NASA

Results for “linear unmixing”

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.

18 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Deep data analysis via physically constrained linear unmixing: universal framework, domain examples, and a community-wide platform

Abstract Many spectral responses in materials science, physics, and chemistry experiments can be characterized as resulting from the superposition of a number of more basic individual spectra. In this context, unmixing is defined as the problem of determining the individual spectra, given measurements of multiple spectra that are spatially resolved across samples, as well as the determination of the corresponding abundance maps indicating the local weighting of each individual spectrum. Matrix factorization is a popular linear unmixing technique that considers that the mixture model between the individual spectra and the spatial maps is linear. Here, we present a tutorial paper targeted at domain scientists to introduce linear unmixing techniques, to facilitate greater understanding of spectroscopic imaging data. We detail a matrix factorization framework that can incorporate different domain information through various parameters of the matrix factorization method. We demonstrate many domain-specific examples to explain the expressivity of the matrix factorization framework and show how the appropriate use of domain-specific constraints such as non-negativity and sum-to-one abundance result in physically meaningful spectral decompositions that are more readily interpretable. Our aim is not only to explain the off-the-shelf available tools, but to add additional constraints when ready-made algorithms are unavailable for the task. All examples use the scalable open source implementation from https://github.com/ramkikannan/nmflibrary that can run from small laptops to supercomputers, creating a user-wide platform for rapid dissemination and adoption across scientific disciplines.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Separating Physically Distinct Mechanisms in Complex Infrared Plasmonic Nanostructures via Machine Learning Enhanced Electron Energy Loss Spectroscopy

Electron energy loss spectroscopy (EELS) enables direct exploration of plasmonic phenomena at the nanometer level. To isolate individual plasmon modes, linear unmixing methods can be used to separate different physical mechanisms, but in larger and more complex systems the interpretability of the components becomes uncertain. Here, infrared plasmonic resonances in self-assembled heterogeneous monolayer films of doped-semiconductor nanoparticles are examined beyond linear unmixing techniques, and both supervised and unsupervised machine-learning-based analyses of hyperspectral EELS datasets are demonstrated. Additionally, in the supervised approach, a human operator labels a small number of pixels in the hyperspectral dataset corresponding to features of interest which are then propagated across the entire dataset. In the unsupervised approach, non-linear autoencoders are used to create a highly-reduced latent-space representation of the dataset, within which insight into the relevant physics can be gleaned from straightforward distance metrics that do not depend on operator input and bias. The advantage of these approaches is that the labeling separates physical mechanisms without altering the data, enabling robust analyses of the influence of heterogeneities in mesoscale complex systems.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Phases and interfaces from real space atomically resolved data: Physics-based deep data image analysis

Advances in electron and scanning probe microscopies have led to a wealth of atomically resolved structural and electronic data, often with ~1–10 pm precision. However, knowledge generation from such data requires the development of a physics-based robust framework to link the observed structures to macroscopic chemical and physical descriptors, including single phase regions, order parameter fields, interfaces, and structural and topological defects. Here, we develop an approach based on a synergy of sliding window Fourier transform to capture the local analog of traditional structure factors combined with blind linear unmixing of the resultant 4D data set. This deep data analysis is ideally matched to the underlying physics of the problem and allows reconstruction of the a priori unknown structure factors of individual components and their spatial localization. We demonstrate the principles of this approach using a synthetic data set and further apply it for extracting chemical and physically relevant information from electron and scanning tunneling microscopy data. Furthermore, this method promises to dramatically speed up crystallographic analysis in atomically resolved data, paving the road toward automatic local structure–property determinations in crystalline and quasi-ordered systems, as well as systems with competing structural and electronic order parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Multifrontal Non-negative Matrix Factorization

Non-negative matrix factorization (Nmf) is an important tool in high-performance large scale data analytics with applications ranging from community detection, recommender system, feature detection and linear and non-linear unmixing. While traditional Nmf works well when the data set is relatively dense, however, it may not extract sufficient structure when the data is extremely sparse. Specifically, traditional Nmf fails to exploit the structured sparsity of the large and sparse data sets resulting in dense factors. We propose a new algorithm for performing Nmf on sparse data that we call multifrontal Nmf (Mf-Nmf) since it borrows several ideas from the multifrontal method for unconstrained factorization (e.g. LU and QR). We also present an efficient shared memory parallel implementation of Mf-Nmf and discuss its performance and scalability. We conduct several experiments on synthetic and realworld datasets and demonstrate the usefulness of the algorithm by comparing it against standard baselines. We obtain a speedup of 1.2x to 19.5x on 24 cores with an average speed up of 10.3x across all the real world datasets.

Sao, Piyush↗
DOE OSTI2020

Quantifying the Dynamics of Protein Self-Organization Using Deep Learning Analysis of Atomic Force Microscopy Data

The dynamics of protein self-assembly on the inorganic surface and the resultant geometric patterns are visualized using high-speed atomic force microscopy. The time dynamics of the classical macroscopic descriptors such as 2D fast Fourier transforms, correlation, and pair distribution functions are explored using the unsupervised linear unmixing, demonstrating the presence of static ordered and dynamic disordered phases and establishing their time dynamics. Here, the deep learning (DL)-based workflow is developed to analyze detailed particle dynamics and explore the evolution of local geometries. Finally, we use a combination of DL feature extraction and mixture modeling to define particle neighborhoods free of physics constraints, allowing for a separation of possible classes of particle behavior and identification of the associated transitions. Overall, this work establishes the workflow for the analysis of the self-organization processes in complex systems from observational data and provides insight into the fundamental mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Visualizing Charge Transport and Nanoscale Electrochemistry by Hyperspectral Kelvin Probe Force Microscopy

Charge-transport and electrochemical processes are heavily influenced by the local microstructure. Kelvin probe force microscopy (KPFM) is a widely used technique to map electrochemical potentials at the nanometer scale; however, it offers little information on local charge dynamics. In this study, we implement a hyperspectral KPFM approach for spatially mapping bias-dependent charge dynamics in timescales ranging from the sub-millisecond to the second regime. As a proof of principle, we investigate the role mobile surface charges play in a three-unit-cell LaAlO 3 /SrTiO 3 oxide heterostructure. We explore machine learning approaches to assist with visualization, pattern recognition, and interpretation of the information-rich data sets. Linear unmixing methods reveal hidden bias-dependent interfacial processes, most likely water splitting, which are essentially unnoticed by functional fitting of the dynamic response alone. Hyperspectral KPFM will be beneficial for investigating nanoscale charge transport and local reactivity in systems involving a possible combination of electronic, ionic, and electrochemical phenomena.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Gaussian process analysis of electron energy loss spectroscopy data: multivariate reconstruction and kernel control

Abstract Advances in hyperspectral imaging including electron energy loss spectroscopy bring forth the challenges of exploratory and physics-based analysis of multidimensional data sets. The multivariate linear unmixing methods generally explore similarities in the energy dimension, but ignore correlations in the spatial domain. At the same time, Gaussian process (GP) explicitly incorporate spatial correlations in the form of kernel functions but is computationally intensive. Here, we implement a GP method operating on the full spatial domain and reduced representations in the energy domain. In this multivariate GP, the information between the components is shared via a common spatial kernel structure, while allowing for variability in the relative noise magnitude or image morphology. We explore the role of kernel constraints on the quality of the reconstruction, and suggest an approach for estimating them from the experimental data. We further show that spatial information contained in higher-order components can be reconstructed and spatially localized.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Local coexistence of VO 2 phases revealed by deep data analysis

We report a synergistic approach of micro-Raman spectroscopic mapping and deep data analysis to study the distribution of crystallographic phases and ferroelastic domains in a defected Al-doped VO 2 microcrystal. Bayesian linear unmixing revealed an uneven distribution of the T phase, which is stabilized by the surface defects and uneven local doping that went undetectable by other classical analysis techniques such as PCA and SIMPLISMA. This work demonstrates the impact of information recovery via statistical analysis and full mapping in spectroscopic studies of vanadium dioxide systems, which is commonly substituted by averaging or single point-probing approaches, both of which suffer from information misinterpretation due to low resolving power.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Multi-scale integration of satellite remote sensing improves characterization of dry-season green-up in an Amazon tropical evergreen forest

In tropical forests, leaf phenology-particularly the pronounced dry-season green-up-strongly regulates biogeochemical cycles of carbon and water fluxes. However, uncertainties remain in the understanding of tropical forest leaf phenology at different spatial scales. Phenocams accurately characterize leaf phenology at the crown and ecosystem scales but are limited to a few sites and time spans of a few years. Time-series satellite observations might fill this gap, but the commonly used satellites (e.g. MODIS, Landsat and Sentinel-2) have resolutions too coarse to characterize single crowns. To resolve this observational challenge, we used the PlanetScope constellation with a 3m resolution and near daily nadir-view coverage. We first developed a rigorous method to cross-calibrate PlanetScope surface reflectance using daily BRDF-adjusted MODIS as the reference. We then used linear spectral unmixing of calibrated PlanetScope to obtain dry-season change in the fractional cover of green vegetation (GV) and non-photosynthetic vegetation (NPV) at the PlanetScope pixel level. We used the Central Amazon Tapajos National Forest k67 site, as all necessary data (from field to phenocam and satellite observations) was available. For this proof of concept, we chose a set of 22 dates of PlanetScope measurements in 2018 and 16 in 2019, all from the six drier months of the year to provide the highest possible cloud-free temporal resolution. Our results show that MODIS-calibrated dry-season PlanetScope data (1) accurately assessed seasonal changes in ecosystem-scale and crown-scale spectral reflectance; (2) detected an increase in ecosystem-scale GV fraction (and a decrease in NPV fraction) from June to November of both years, consistent with local phenocam observations with R 2 around 0.8; and (3) monitored large seasonal trend variability in crown-scale NPV fraction. Finally, our results highlight the potential of integrating multi-scale satellite observations to extend fine-scale leaf phenology monitoring beyond the spatial limits of phenocams.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

An ecologically-constrained deep learning model for tropical leaf phenology monitoring using PlanetScope satellites

In tropical forests, leaf phenology signals leaf-on/off status and exhibits considerable variability across scales from a single tree-crown to the entire forest ecosystem. Such phenology signals importantly regulate large-scale biogeochemical cycles and regional climate. PlanetScope CubeSats data with a 3-m resolution and near-daily global coverage provide an unprecedented opportunity to monitor both fine- and ecosystem-scale phenology variability along large environmental gradients. However, a scalable method that accurately characterizes leaf phenology from PlanetScope with biophysically meaningful metrics remains lacking. Here we developed an index-guided, ecologically constrained autoencoder (IG-ECAE) method to automatically derive a deciduousness metric (i.e. percentage of tree canopies with leaf-off status within an image pixel) from PlanetScope. The IG-ECAE first estimated the reflectance spectra of leafy/leafless canopies based on their spectral indices characteristics, then used the derived reflectance spectra to guide an autoencoder deep learning method with additional ecological constraints to refine the reflectance spectra, and finally used linear spectral unmixing to estimate the relative abundance of leafless canopies (or deciduousness) per PlanetScope image pixel. We tested the IG-ECAE method at 16 tropical forest sites spanning multiple continents and a large precipitation gradient (1470-2819 mm year -1 ). Among these sites, we evaluated the PlanetScope-derived deciduousness with corresponding measures derived from WorldView-2 (n=9 sites) and local phenocams (n=9 sites). Our results show that PlanetScope-derived deciduousness agrees: 1) with that derived from WorldView-2 at the patch level (90m×90m) with r 2 =0.89 across all sites; and 2) with that derived from phenocams to quantify ecosystem-scale seasonality with r 2 ranging from 0.62-0.96. These results demonstrate the effectiveness and scalability of IG-ECAE in characterizing the wide variability in deciduousness across scales from pixels to forest ecosystems, and from a single date to the full annual cycle, indicating the potential for using high-resolution satellites to track the large-scale phenological patterns and response of tropical forests to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Coso Geothermal Spectral Library for Rocks and Minerals

An integrated open mineral spectral library designed to enhance the utility and precision of mineral spectral data for geothermal exploration, developed from a reliable and comprehensive digital dataset for seamless sharing by integrating field data, the USGS spectral library, and pertinent information obtained from Coso geothermal field (Coso) in California. The ASD FieldSpec portable spectrometer was utilized for collecting spectral data, which was subsequently analyzed using the THOR Material Identification tool in ENVI, The Spectral Geologist (TSG) software by CSIRO, and the Fully Constrained Linear Spectral Unmixing algorithm (FCLSU) in MATLAB. Scanning Electron Microscopy (SEM) with a mineralogy-analyzing function was employed to assess the mineral composition of samples, ensuring precise mineralogical analysis. A portable X-ray fluorescence (pXRF) spectrometer was also utilized to gather information on elemental enrichment. A framework for developing spectra data and establishing spectral libraries for various geological cases was proposed within this study. The characteristic spectra of six alteration minerals - alunite, chalcedony, epidote, hematite, kaolinite, and opal - were acquired from Coso samples. The spectral library for the Coso alteration minerals was introduced for further application in academic study or industrial exploration. To browse the Coso Geothermal Spectral data and related figures from spreadsheets: #1 Unzip and store the following items in the same folder. 'Contact Probe Data.zip', 'Sample Photos.zip', and 'Coso spectra of higher-certainity minerals.xlsx'. #2 Open 'Coso spectra of higher-certainity minerals.xlsx'. The hyperlinks in the spreadsheet lead to the folders or figures of: spectra .asd file, spectra ASC II file, spectra plots, and sample photos. The spectra data is raw data without splice correction. Spectra .asd files require particular software to open. (These cannot be opened in GIS software such as ArcGIS.) Spectra ASC II files can be opened in a text editor or spread sheet program.

15 GEOTHERMAL ENERGY↗
DOE OSTI2018

Overhead longwave infrared hyperspectral material identification using radiometric models

Material detection algorithms used in hyperspectral data processing are computationally efficient but can produce relatively high numbers of false positives. Material identification performed as a secondary processing step on detected pixels can help mitigate false positives. A material identification processing chain for longwave infrared hyperspectral data of solid materials collected from airborne platforms is presented. The algorithms utilize unwhitened radiance data and Nelder–Meade numerical optimization to estimate the temperature, humidity, and ozone levels of the atmospheric profile. Pixel unmixing is done using constrained linear regression and Bayesian information criteria for model selection. The resulting identification product includes an optimal atmospheric profile and a full radiance material model that includes material temperature, abundance values, and several fit statistics. A logistic regression method utilizing the model parameters to improve identification is also presented. Furthermore, several examples are provided using modeled data at several noise levels.

58 GEOSCIENCES↗
DOE OSTI2018

Overhead longwave infrared hyperspectral material identification using radiometric models

Material detection algorithms used in hyperspectral data processing are computationally efficient but can produce relatively high numbers of false positives. Material identification performed as a secondary processing step on detected pixels can help separate true and false positives. This paper presents a material identification processing chain for longwave infrared hyperspectral data of solid materials collected from airborne platforms. The algorithms utilize unwhitened radiance data and an iterative algorithm that determines the temperature, humidity, and ozone of the atmospheric profile. Pixel unmixing is done using constrained linear regression and Bayesian Information Criteria for model selection. The resulting product includes an optimal atmospheric profile and full radiance material model that includes material temperature, abundance values, and several fit statistics. A logistic regression method utilizing all model parameters to improve identification is also presented. This paper details the processing chain and provides justification for the algorithms used. Several examples are provided using modeled data at different noise levels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Interpreting mixing relationships in energetic melts to estimate vapor contribution and composition

Heterogeneous, glassy debris can be produced by energetic events such as near surface nuclear explosions and asteroid impacts. In these events, the resultant debris comprises a rapidly-cooled mixture of the projectile or bomb material with proximate materials such as soil and anthropogenic structures. Glasses formed in these events preserve inter and/or intra sample compositional variation reflecting complex mixtures of the precursor materials or their chemically fractionated products. Understanding formation mechanisms and chemical processes relating to energetic glass formation based on study of the end products requires an objective method for calculating the precursor components that mixed together to produce the final product. Here, we compare a set of multivariate linear-least squares approaches to address the unmixing problem of heterogeneous glasses including classical least squares (CLS), principal component analysis (PCA) and two multivariate curve resolution – alternating least squares (MCR-ALS) approaches. We generate a synthetic data set representing spatially resolved, compositional data in a set of heterogeneously mixed glasses to understand model sensitivity to precursor abundance and composition. We explore the impact of incomplete knowledge of the system prior to the formation of the melts on the resulting confidence and accuracy of the solution (i.e., the number, composition, and abundances of the precursor materials). We conclude that a closure and equality constrained MCR-ALS approach succeeds in identifying unrecognized precursor components and is a suitable alternative to CLS approaches when the investigator has incomplete knowledge of the starting precursor chemical compositions. Finally, we demonstrate application of this MCR-ALS approach to measurements of Trinitite.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2018

Using a 2-shock 1D platform at NIF to measure the effect of convergence on mix and symmetry

In this paper, we describe the use of a robust new 1-D like implosion platform at the National Ignition Facility [G. H. Miller et al., Opt. Eng. 43, 2841 (2004)] to study the effect of convergence on mix and shape. Previous experiments suggest that nuclear yields and ion temperature degrade with increased convergence [M. D. Cable et al., Phys. Rev. Lett. 73, 2316 (1994)] due to enhanced perturbation growth and mix, but little has been reported on the distortion of the shape with time. The 2-shock platform was developed [S. F. Khan et al., Phys. Plasmas 23, 042708 (2016)] to maintain a high degree of sphericity during the whole implosion phase and has a thick, uniformly doped (1% Si) plastic CH shell to minimize the effect of mixing due to hydrodynamic feed-through from the outer ablator surface. An inner layer of deuterated plastic (CD) and hydrogen-tritium (DT) gas fill allows for the measurement of DT neutrons produced by the mix between the gas and ablator. DD neutrons provide information about the hot, unmixed CD region. By changing the fill gas density while keeping the capsule diameter, ablator thickness, and Au hohlraum conditions fixed, the x-ray hot spot convergence ratio was varied from 14 to 22. We find that the atomic mix (DT yield) grows linearly as a function of convergence, but since T ion changes as well, it does not necessarily mean that the amount or extent of mix grows linearly as well. Lastly, we also find the DD yield, which is a measurement of the shell heating, saturates above a certain convergence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Impact of Advection Schemes on Tracer Interrelationships in Large-Eddy Simulations of Deep Convection

This study investigates the preservation of tracer interrelationships during advection in large-eddy simulations of an idealized deep convective cloud, which is particularly relevant to chemistry, aerosol, and cloud microphysics models. Employing the Cloud Model 1, advection is represented using third-, fifth-, and seventh-order weighted essentially non-oscillatory schemes. As a simplified analogy for cloud hydrometeors and aerosols, several inert passive tracers following linear and nonlinear relationships are initialized after the cloud reaches ~6-km depth. Numerical mixing in the simulated turbulent convective clouds leads to significant deviations from the initial nonlinear relationships between tracers. In these simulations, a considerable fraction of the grid points where the tracers’ nonlinear relationships are altered from advection are classified as unrealistic (e.g., ~13% for the environmental tracers on average), including errors from range-preserving unmixing and overshooting. Errors in the sum of three tracers are also relatively large, ranging between ~1% and 16% for 5% of the grid points in and near the cloud. The magnitude of unrealistic mixing and errors in the sum of three tracers generally increase with the order of accuracy of the advection scheme. These results are consistent across model grid spacings ranging from 50 to 200 m, and across three different flow realizations for each combination of grid spacing and advection scheme tested. Tests employing a previously proposed scalar normalization procedure show substantially reduced errors in the sum of three tracers with a relatively small negative impact on other tracer relationships. This analysis, therefore, suggests efficacy of the normalization procedure when applied to turbulent three-dimensional cloud simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Integrating very-high-resolution UAS data and airborne imaging spectroscopy to map the fractional composition of Arctic plant functional types in Western Alaska

Widespread changes in vegetation cover and composition are driving strong impacts on Arctic ecosystem functioning and global climate feedbacks. An accurate characterization of tundra vegetation composition is required to understand how the Arctic will respond to future climate change. However, quantifying tundra vegetation composition over large areas is challenging as commonly-used satellite observations are too coarse, spatially and spectrally, to differentiate low-lying tundra vegetation types. Recent airborne and spaceborne imaging spectroscopy platforms provide better data to characterize vegetation composition. Yet, our ability to characterize vegetation composition with imaging spectroscopy remains largely unexplored in the Arctic, particularly due to a lack of ground observations needed to train and test classification models. To address this problem, we collected very-high-resolution (VHR, ~5 cm) unoccupied aerial system (UAS) imagery at three low-Arctic tundra sites located on the Seward Peninsula, western Alaska. In this paper, we examine the feasibility of integrating imagery from the UAS and the hyperspectral Airborne Visible/Infrared Imaging Spectrometer, Next Generation (AVIRIS-NG) airborne instrument to map the fractional composition of 12 key Arctic plant functional types (PFTs). To this end, we first mapped the 12 PFTs from our VHR UAS imagery using random forest classification. We then used these UAS-derived PFT maps as ground truth to develop partial least squares regression (PLSR) models to predict the fractional cover (FCover) of each PFT from AVIRIS-NG imagery. Further, we evaluated the performance of our PLSR models using reserved UAS samples, as well as by mapping PFT FCover and dominant PFT for large tundra landscapes. Our results show that 1) Arctic PFTs can be effectively mapped using VHR UAS imagery, with overall accuracy between 86% and 92%, 2) when the UAS mapped PFTs were used to inform PLSR scaling models, the FCover of the 12 PFTs could be effectively estimated from AVIRIS-NG imagery with a mean absolute error (MAE) <0.13, and 3) our PLSR models outperformed traditional, fully constrained least-squares (FCLS) linear mixture analysis and produced high-quality, spatially contiguous PFT FCover and PFT maps that captured vegetation spatial patterns with similar accuracy to those developed from UAS imagery. The developed PLSR models have the potential to be broadly applied for quantifying vegetation composition with AVIRIS-NG images to help monitor tundra vegetation dynamics and improve process-based modeling of tundra ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Development Of Nutrient And Water Recycling Capabilities In Algae Biofuels Production Systems. Final Summary Report

The objective of this project was to develop and demonstrate methods of recycling of water and nutrients for algal biofuels production. Recycling was accomplished both internal to the system and, in a broader sense, through import and reuse of municipal wastewater. Such an integrated system with wastewater input had not been demonstrated previously, and the performance was unknown, particularly in terms of influence of recycling on algal productivity and the practical extent of nutrient recovery from biomass residuals. Through long-term laboratory and pilot research, the project resulted in the following: 1. Bench-scale pretreatment of algal biomass did not sufficiently increase methane yield of nutrient solubilization during anaerobic digestion to warrant incorporation of pre-treatment into the pilot plant. The trial pretreatments were high-pressure orifice homogenization, sonication, and two types of heat treatment. 2. Solubilization of biomass particulate nutrients by lab anaerobic digesters ranged from 20% to nearly 60% for N and 40-65% for P. Subsequent aerobic degradation of the anaerobically digested biomass simulated raceways receiving whole digestate and resulted in an additional 20-55% N solubilization and additional 20% P solubilization. 3. Comparisons of laboratory and pilot digesters showed that laboratory units were reasonable proxies for pilot-scale. 4. Pilot-scale anaerobic digesters were designed, installed, and operated to digest algal biomass. Nutrient re-solubilization by the digesters was monitored and whole digestate was successfully used as a fertilizer in pilot algae raceways. 5. Unheated, unmixed digesters achieved greater methane yield and nutrient solubilization than heated, mixed digesters, presumably due to longer the solids residence times in unmixed digesters. The unmixed, unheated pilot digesters yielded 0.16 L CH4 /g volatile solids (VS) introduced with 0.15 g VS/L-d organic loading and 16oC average temperature. A conventional heated mixed lab digester yielded 0.22 L CH4 /g VS with 0.25 g VS/L-d and 30oC. The highest yield (0.30 L CH4 /g VS) was achieved by the unmixed lab digesters operated at a constant 20oC. All digesters were operated with a 40-d hydraulic residence time. 6. In general, 50-75% of initial particulate N and P could be solubilized during anaerobic digestion and available for subsequent rounds of algae cultivation. 7. Bench-scale experiments showed the recovery from hydrothermal liquefaction (HTL) wastewater of carbon via anaerobic digestion and of nutrients to grow algae. To satisfy the nitrogen demand of algae cultivation, HTL wastewater would be diluted 400-fold, which was found to eliminate inhibition of algae growth by HTL wastewater. 8. Anaerobic digestion methane yield was lower for algal biomass containing coagulants such as would be used to aid harvesting or dewatering. Depending on doses, starch-based coagulant decreased yield by 10-14% and aluminum chlorohydrate decreased it by 14-26%. The lowest yield was 0.28 L CH 4 /g volatile solids introduced to the digesters. 9. Algae harvested from raceways operated on recycled water had methane yields 13% higher than algae from raceways operated on both recycled water and nutrients provided by algae digestate. The slightly lower yield was expected due to the presence of previously digested biomass from the digestate fertilizer. 10. Defined media was replenished with nutrients and recycled repeatedly in sequential batch growth of Chlorella sorokiniana (DOE 1412). This laboratory study tested for inhibition and accumulation of inhibiting compounds (allelopathic or auto-inhibitory substances), information that would help estimate the blowdown ratio needed for an integrated system. In laboratory experiments in which water was recycled a total of five times, each successive round of reuse resulted in an average 4±3% reduction in log-phase specific growth rates. However, linear-phase growth inhibition was only observed in the final fifth round of reuse. 11. No decline in productivity was detected after 15 rounds of water recycling with nutrients provided by whole digestate in lab cultivation. Lab tests allowed for steady light and temperature, increasing the ability to detect inhibition. 12. In initial pilot inhibition studies, wastewater growth media was reused once while productivity was monitored. Media reuse was accomplished with triplicate sets of 33-m2 raceways operated in series. First-round gross productivity (based on effluent biomass flow) averaged 23 g/m 2 -day annually while second-round gross productivity averaged 19 g/m 2 -day annually. In terms of net productivity (based on raceway effluent biomass minus influent biomass), the first-round productivity averaged 15 g/m2-d and second round averaged 13 g/m 2 -d during June-September operation. The higher productivity in the first-round ponds was likely due to heterotrophic/mixotrophic growth on the wastewater organic matter. 13. In a culminating pilot experiment, coagulant was used to decrease the carry-over of unsettled algae into subsequent rounds of growth. Over nearly 8 months, 93% of the media (the equivalent of 14 rounds of water reuse) was recycled without significant productivity loss compared to controls. Ponds receiving both recycled water and nutrients had net productivities of 14-24 g/m 2 -d during fall and mid-summer, respectively. 14. Techno-economic analysis of the proposed facility found minimum fuel selling price to range from $7.01/gallon gasoline equivalent without revenue other than fuel to $3.85/GGE with revenue from wastewater treatment fees and LCFS and RIN (Low Carbon Fuel Standard and Renewable Identification Numbers) credits. 15. Life cycle assessment indicated GHG emissions of 40.7 g CO 2 /MJ fuel and a net energy ratio (energy required/energy produced) of 0.37.

09 BIOMASS FUELS↗