Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spectral mixture analysis”

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 19 records

Prior-knowledge-based spectral mixture analysis for impervious surface mapping

In this study, we developed a prior-knowledge-based spectral mixture analysis (PKSMA) to map impervious surfaces by using endmembers derived separately for high- and low-density urban regions. First, an urban area was categorized into high- and low-density urban areas, using a multi-step classification method. Next, in high-density urban areas that were assumed to have only vegetation and impervious surfaces (ISs), the Vegetation-Impervious model (V-I) was used in a spectral mixture analysis (SMA) with three endmembers: vegetation, high albedo, and low albedo. In low-density urban areas, the Vegetation-Impervious-Soil model (V-I-S) was used in an SMA analysis with four endmembers: high albedo, low albedo, soil, and vegetation. The fraction of IS with high and low albedo in each pixel was combined to produce the final IS map. The root mean-square error (RMSE) of the IS map produced using PKSMA was about 11.0%, compared to 14.52% using four-endmember SMA. Particularly in high-density urban areas, PKSMA (RMSE = 6.47%) showed better performance than four-endmember (15.91%). The results indicate that PKSMA can improve IS mapping compared to traditional SMA by using appropriately selected endmembers and is particularly strong in high-density urban areas.

Impervious surface↗

Sensitivity of Optical Satellites to Estimate Windthrow Tree-Mortality in a Central Amazon Forest

Windthrow (i.e., trees broken and uprooted by wind) is a major natural disturbance in Amazon forests. Images from medium-resolution optical satellites combined with extensive field data have allowed researchers to assess patterns of windthrow tree-mortality and to monitor forest recovery over decades of succession in different regions. Although satellites with high spatial-resolution have become available in the last decade, they have not yet been employed for the quantification of windthrow tree-mortality. Here, we address how increasing the spatial resolution of satellites affects plot-to-landscape estimates of windthrow tree-mortality. We combined forest inventory data with Landsat 8 (30 m pixel), Sentinel 2 (10 m), and WorldView 2 (2 m) imagery over an old-growth forest in the Central Amazon that was disturbed by a single windthrow event in November 2015. Remote sensing estimates of windthrow tree-mortality were produced from Spectral Mixture Analysis and evaluated with forest inventory data (i.e., ground true) by using Generalized Linear Models. Field measured windthrow tree-mortality (3 transects and 30 subplots) crossing the entire disturbance gradient was 26.9 ± 11.1% (mean ± 95% CI). Although the three satellites produced reliable and statistically similar estimates (from 26.5% to 30.3%, p < 0.001), Landsat 8 had the most accurate results and efficiently captured field-observed variations in windthrow tree-mortality across the entire gradient of disturbance (Sentinel 2 and WorldView 2 produced the second and third best results, respectively). As expected, mean-associated uncertainties decreased systematically with increasing spatial resolution (i.e., from Landsat 8 to Sentinel 2 and WorldView 2). However, the overall quality of model fits showed the opposite pattern. We suggest that this reflects the influence of a relatively minor disturbance, such as defoliation and crown damage, and the fast growth of natural regeneration, which were not measured in the field nor can be captured by coarser resolution imagery. Our results validate the reliability of Landsat imagery for assessing plot-to-landscape patterns of windthrow tree-mortality in dense and heterogeneous tropical forests. Satellites with high spatial resolution can improve estimates of windthrow severity by allowing the quantification of crown damage and mortality of lower canopy and understory trees. However, this requires the validation of remote sensing metrics using field data at compatible scales.

54 ENVIRONMENTAL SCIENCES↗

Hurricane Maria Puerto Rico Landsat Analysis

Hurricane Maria made landfall as a strong Category 4 storm in southeast Puerto Rico on September 20th, 2018. The powerful storm traversed the island in a northwesterly direction causing widespread destruction. This study focused on a rapid assessment of Hurricane Marias impact to Puerto Ricos forests. Calibrated and corrected Landsat 8 image composites for the entire island were generated using Google Earth Engine for a comparable pre-Maria and post-Maria time period that accounted for phenology. Spectral mixture analysis (SMA) using image-derived end members was carried out on both composites to calculate the change in the non-photosynthetic vegetation (Delta-NPV) spectral response, a metric that quantifies the increased fraction of exposed wood and surface litter associated with tree mortality and crown damage from the storm. Hurricane simulations were also conducted using the Weather Research and Forecasting (WRF) regional climate model to estimate wind speeds associated with forest disturbance. Dramatic changes in forest structure across the entire island were evident from pre- and post-Maria composited Landsat 8 images. A Delta-NPV map for only the forested pixels illustrated significant spatial variability in disturbance, with patterns that associated with factors such as slope, aspect and elevation. An initial order-of-magnitude impact estimate based on previous work indicated that Hurricane Maria may have caused mortality and severe damage to 23-31 million trees. Additional field work and image analyses are required to further detail the impact of Hurricane Maria to Puerto Rico forests. A minor update to this dataset was posted on April 20, 2018. The previous version is being retired. If you need access to the prior version of the data, email ngee-tropics-archive@lbl.gov.

54 ENVIRONMENTAL SCIENCES↗

Spectral analysis of amazon canopy phenology during the dry season using a tower hyperspectral camera and modis observations

The association between spectral reflectance and canopy processes remains challenging for quantifying large-scale canopy phenological cycles in tropical forests. In this paper, we used a tower-mounted hyperspectral camera in an eastern Amazon forest to assess how canopy spectral signals of three species are linked with phenological processes in the 2012 dry season. We explored different approaches to disentangle the spectral components of canopy phenology processes and analyze their variations over time using 17 images acquired by the camera. The methods included linear spectral mixture analysis (SMA); principal component analysis (PCA); continuum removal (CR); and first-order derivative analysis. In addition, three vegetation indices potentially sensitive to leaf flushing, leaf loss and leaf area index (LAI) were calculated: the Enhanced Vegetation Index (EVI), Normalized Difference Vegetation Index (NDVI) and the entitled Green-Red Normalized Difference (GRND) index. We inspected also the consistency of the camera observations using Moderate Resolution Imaging Spectroradiometer (MODIS) and available phenological data on new leaf production and LAI of young, mature and old leaves simulated by a leaf demography-ontogeny model. The results showed a diversity of phenological responses during the 2012 dry season with related changes in canopy structure and greenness values. Because of the differences in timing and intensity of leaf flushing and leaf shedding, Erisma uncinatum, Manilkara huberi and Chamaecrista xinguensis presented different green vegetation (GV) and non-photosynthetic vegetation (NPV) SMA fractions; distinct PCA scores; changes in depth, width and area of the 681-nm chlorophyll absorption band; and variations over time in the EVI, GRND and NDVI. At the end of dry season, GV increased for Erisma uncinatum, while NPV increased for Chamaecrista xinguensis. For Manilkara huberi, the NPV first increased in the beginning of August and then decreased toward September with new foliage. Variations in red-edge position were not statistically significant between the species and across dates at the 95% confidence level. The camera data were affected by view-illumination effects, which reduced the SMA shade fraction over time. When MODIS data were corrected for these effects using the Multi-Angle Implementation of Atmospheric Correction Algorithm (MAIAC), we observed an EVI increase toward September that closely tracked the modeled LAI of mature leaves (3–5 months). Compared to the EVI, the GRND was a better indicator of leaf flushing because the modeled production of new leaves peaked in August and then declined in September following the GRND closely. Finally, while the EVI was more related to changes in mature leaf area, the GRND was more associated with new leaf flushing.

54 ENVIRONMENTAL SCIENCES↗

Hyperspectral remote sensing of fire: state-of-the-art and future perspectives

We report that fire is a widespread Earth system process with important carbon and climate feedbacks. Multispectral remote sensing has enabled mapping of global spatiotemporal patterns of fire and fire effects, which has significantly improved our understanding of interactions between ecosystems, climate, humans and fire. With several upcoming spaceborne hyperspectral missions like the Environmental Mapping And Analysis Program (EnMAP), the Hyperspectral Infrared Imager (HyspIRI) and the Precursore Iperspettrale Della Missione Applicativa (PRISMA), we provide a review of the state-of-the-art and perspectives of hyperspectral remote sensing of fire. Hyperspectral remote sensing leverages information in many (often more than 100) narrow (smaller than 20 nm) spectrally contiguous bands, in contrast to multispectral remote sensing of few (up to 15) non-contiguous wider (greater than 20 nm) bands. To date, hyperspectral fire applications have primarily used airborne data in the visible to short-wave infrared region (VSWIR, 0.4 to 2.5 μm). This has resulted in detailed and accurate discrimination and quantification of fuel types and condition, fire temperatures and emissions, fire severity and vegetation recovery. Many of these applications use processing techniques that take advantage of the high spectral resolution and dimensionality such as advanced spectral mixture analysis. So far, hyperspectral VSWIR fire applications are based on a limited number of airborne acquisitions, yet techniques will approach maturity for larger scale application when spaceborne imagery becomes available. Recent innovations in airborne hyperspectral thermal (8 to 12 μm) remote sensing show potential to improve retrievals of temperature and emissions from active fires, yet these applications need more investigation over more fires to verify consistency over space and time, and overcome sensor saturation issues. Furthermore, hyperspectral information and structural data from, for example, light detection and ranging (LiDAR) sensors are highly complementary. Finally, their combined use has demonstrated advantages for fuel mapping, yet its potential for post-fire severity and combustion retrievals remains largely unexplored.

54 ENVIRONMENTAL SCIENCES↗

Using high spatial resolution satellite imagery to map forest burn severity across spatial scales in a Pine Barrens ecosystem

As a primary disturbance agent, fire significantly influences local processes and services of forest ecosystems. Although a variety of remote sensing based approaches have been developed and applied to Landsat mission imagery to infer burn severity at 30 m spatial resolution, forest burn severity have still been seldom assessed at fine spatial scales (≤ 5 m) from very-high-resolution (VHR) data. Here we assessed a 432 ha forest fire that occurred in April 2012 on Long Island, New York, within the Pine Barrens region, a unique but imperiled fire-dependent ecosystem in the northeastern United States. The mapping of forest burn severity was explored here at fine spatial scales, for the first time using remotely sensed spectral indices and a set of Multiple Endmember Spectral Mixture Analysis (MESMA) fraction images from bi-temporal — pre- and post-fire event — WorldView-2 (WV-2) imagery at 2 m spatial resolution. We first evaluated our approach using 1 m by 1 m validation points at the sub-crown scale per severity class (i.e. unburned, low, moderate, and high severity) from the post-fire 0.10 m color aerial ortho-photos; then, we validated the burn severity mapping of geo-referenced dominant tree crowns (crown scale) and 15 m by 15 m fixed-area plots (inter-crown scale) with the post-fire 0.10 m aerial ortho-photos and measured crown information of twenty forest inventory plots. Our approach can accurately assess forest burn severity at the sub-crown (overall accuracy is 84% with a Kappa value of 0.77), crown (overall accuracy is 82% with a Kappa value of 0.76), and inter-crown scales (89% of the variation in estimated burn severity ratings (i.e. Geo-Composite Burn Index (CBI)). Lastly, this work highlights that forest burn severity mapping from VHR data can capture heterogeneous fire patterns at fine spatial scales over the large spatial extents. This is important since most ecological processes associated with fire effects vary at the < 30 m scale and VHR approaches could significantly advance our ability to characterize fire effects on forest ecosystems.

54 ENVIRONMENTAL SCIENCES↗

Nutrient Cycle Impacts on Forest Ecosystem Carbon Cycling: Improved Prediction of Climate Feedbacks from Coupled C–Nutrient Dynamics from Ecosystem to Regional Scales

One of the largest sources of uncertainty in the projected global mean temperature response to anthropogenic carbon (C) emissions is the role of the land surface in contributing to or mediating the rise in atmospheric CO2 [IPCC, 2007]. Development of ecosystem processes in global land surface models (LSMs) is far from complete, with major components such as the nitrogen (N) cycle only recently being implemented in some LSMs [Thornton et al., 2007; Sokolov et al., 2008; Xu-Ri and Prentice, 2008; Fisher et al., 2010]. Inclusion of the N cycle can reduce terrestrial C uptake by up to ~50 Pg C·y-1 globally [Figure 1: Fisher et al., 2010]. However, even these early representations of the N cycle in LSMs are missing critical plant–soil interactions that have been shown empirically to exert a major influence over ecosystem C fluxes. Here we propose to use a newly developed plant N model by PI Fisher as a theoretical and analytical framework for integrating novel field-based methods and long-term flux observations to quantify the impacts of belowground processes on fluxes of C, water and energy, from ecosystem to regional scales. Although the model (Fixation & Uptake of Nitrogen: FUN) is one of the more advanced global plant N models available, the model is dependent on theoretical C flux curves optimizing plant C allocation to root exudates, symbiotic bacteria and fungi, and retranslocation enzymes that have not been confirmed or validated with measurements—primarily because the measurement techniques had not been available, until recently. Parallel to the model developments by Fisher et al. were developments in in situ direct measurement techniques by Phillips et al. [2008; 2011], who were able to measure C fluxes from roots to soil. Not only do these techniques have the capacity to confirm and/or modify the modeling components in the Fisher et al. plant N model, but they also provide insight into how belowground C supply from roots couples the C–N cycles at the micro-site and ecosystem scales. Collectively, the integration of measurements and modeling provides empirical grounding for theoretical process behavior. Three critical components will enable us to successfully carry out this research: I) new methodological approaches developed by co-I Phillips for in situ measurement of root exudates, and belowground C fluxes that allow for testing and refinement of FUN [Phillips et al., 2008; Phillips et al., 2011]; II) an innovative research design that involves data integration across an ecological gradient, including a DOE-supported AmeriFlux tower operated by co-I Dragoni, which is essential to determining the ecosystem flux impacts [Dragoni et al., 2011]; and, III) new results from co-I Evans showing unique spectral signatures of trees in association with ectomycorrhizae (ECM) and arbuscular mycorrhizae (AM) fungi from Landsat data using spectral mixture analysis (SMA), which is key to scaling up belowground processes to regional scales [Andersson et al., 2009]. Each component provides feedback to and is intimately linked with all other components; no part operates in isolation. Our team and approach combines expertise in modeling, remote sensing, ecosystem ecology, plant physiology, soil biogeochemistry, and micrometeorology to address these inherently interdisciplinary problems and questions. This proposal directly responds to the focused Science Area of the FOA: “The role of belowground processes and mechanisms across scales (e.g., soil carbon transformation, root dynamics, mycorrhizal interactions, and plant mediated (e.g., root exudates) biogeochemical transformations) associated with a changing climate.

54 ENVIRONMENTAL SCIENCES↗

Landscape-scale consequences of differential tree mortality from catastrophic wind disturbance in the Amazon

Wind disturbance can create large forest blowdowns, which greatly reduces live biomass and adds uncertainty to the strength of the Amazon carbon sink. Observational studies from within the central Amazon have quantified blowdown size and estimated total mortality but have not determined which trees are most likely to die from a catastrophic wind disturbance. Also, the impact of spatial dependence upon tree mortality from wind disturbance has seldom been quantified, which is important because wind disturbance often kills clusters of trees due to large treefalls killing surrounding neighbors. We examine (1) the causes of differential mortality between adult trees from a 300-ha blowdown event in the Peruvian region of the northwestern Amazon, (2) how accounting for spatial dependence affects mortality predictions, and (3) how incorporating both differential mortality and spatial dependence affect the landscape level estimation of necromass produced from the blowdown. Standard regression and spatial regression models were used to estimate how stem diameter, wood density, elevation, and a satellite-derived disturbance metric influenced the probability of tree death from the blowdown event. The model parameters regarding tree characteristics, topography, and spatial autocorrelation of the field data were then used to determine the consequences of non-random mortality for landscape production of necromass through a simulation model. Tree mortality was highly non-random within the blowdown, where tree mortality rates were highest for trees that were large, had low wood density, and were located at high elevation. Of the differential mortality models, the non-spatial models overpredicted necromass, whereas the spatial model slightly underpredicted necromass. When parameterized from the same field data, the spatial regression model with differential mortality estimated only 7.5% more dead trees across the entire blowdown than the random mortality model, yet it estimated 51% greater necromass. As a result, we suggest that predictions of forest carbon loss from wind disturbance are sensitive to not only the underlying spatial dependence of observations, but also the biological differences between individuals that promote differential levels of mortality.

54 ENVIRONMENTAL SCIENCES↗

Interfacial Nanostructure and Hydrogen Bond Networks of Choline Chloride and Glycerol Mixtures Probed with X-ray and Vibrational Spectroscopies

The molecular distribution at the liquid-vapor interface and evolution of the hydrogen bond interactions in mixtures of glycerol and choline chloride are investigated using X-ray photoelectron spectroscopy (XPS) and Raman spectroscopy. Nanoscale depth profiles of supersaturated deep eutectic solvent (DES) mixtures up to ~2 nm measured by ambient-pressure XPS show the enhancement of choline cation (Ch + ) concentration by a factor of 2 at the liquid-vapor interface compared to the bulk. In addition, Raman spectral analysis of a wide range of DES mixtures reveals the conversion of gauche-conformer Ch + into the anti-conformer in relatively lower ChCl concentrations. Finally, the depletion of Ch + from the interface (probing depth = 0.4 nm) is demonstrated by aerosol-based velocity map imaging XPS measurements of glyceline and water mixtures. The nanostructure of liquid-vapor interfaces and structural rearrangement by hydration can provide critical insight into the molecular origin of the deep eutectic behavior and gas-capturing application of DESs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Intermolecular polarizabilities in H{sub 2}-rare-gas mixtures (H{sub 2}–He, Ne, Ar, Kr, Xe): Insight from collisional isotropic spectral properties

The report presents results of theoretical and numerical analysis of the electrical properties related to the isotropic part of the polarizability induced by interactions within compounds built up of a hydrogen H{sub 2} molecule and a set of noble gas atoms, Rg, ranging from the least massive helium up to the heaviest xenon perturber. The Cartesian components of the collisional polarizabilities of the H{sub 2}–Rg systems are found by means of the quantum chemistry methods and their dependence on the intermolecular distance is determined. On the basis of these data, the spherical, symmetry adapted components of the trace polarizability are derived in order to provide data sets that are convenient for evaluating collisional spectral profiles of the isotropic polarized part of light scattered by the H{sub 2}–Rg mixtures. Three independent methods of numerical computing of the spectral intensities are applied at room temperature (295 K). The properties of the roto-translational profiles obtained are discussed in order to determine the role played by contributions corresponding to each of the symmetry adapted parts of the trace polarizability. By spreading the analysis over the collection of the H{sub 2}–Rg systems, evolution of the spectral properties with the growing masses of the supermolecular compounds can be observed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deciphering phase evolution in complex metal oxide thin films via high-throughput materials synthesis and characterization

Discovery of structure-property relationships in thin film alloys of complex metal oxides enabled by high-throughput materials synthesis and characterization facilities is demonstrated here with a case-study. In this study, thin films of binary transition metal oxides (Ti–Zn) are prepared by pulsed laser deposition with continuously varying Ti:Zn ratio, creating combinatorial samples for exploration of the properties of this material family. The atomic structure and electronic properties are probed by spatially resolved techniques including x-ray absorption near edge structures (XANES) and x-ray fluorescence (XRF) at the Ti and Zn K-edge, x-ray diffraction, and spectroscopic ellipsometry. The observed properties as a function of Ti:Zn ratio are resolved into mixtures of five distinguishable phases by deploying multivariate curve resolution analysis on the XANES spectral series, under constraints set by results from the other characterization techniques. First-principles computations based on density function theory connect the observed properties of each distinct phase with structural and spectral characteristics of crystalline polymorphs of Ti–Zn oxide. Continuous tuning of the optical absorption edge as a function of Ti:Zn ratio, including the unusual observation of negative optical bowing, exemplifies a functional property of the film correlated to the phase evolution.

36 MATERIALS SCIENCE↗

Prediction of the Cu oxidation state from EELS and XAS spectra using supervised machine learning

Abstract Electron energy loss spectroscopy (EELS) and X-ray absorption spectroscopy (XAS) provide detailed information about bonding, distributions and locations of atoms, and their coordination numbers and oxidation states. However, analysis of XAS/EELS data often relies on matching an unknown experimental sample to a series of simulated or experimental standard samples. This limits analysis throughput and the ability to extract quantitative information from a sample. In this work, we have trained a random forest model capable of predicting the oxidation state of copper based on its L-edge spectrum. Our model attains an R 2 score of 0.85 and a root mean square error of 0.24 on simulated data. It has also successfully predicted experimental L-edge EELS spectra taken in this work and XAS spectra extracted from the literature. We further demonstrate the utility of this model by predicting simulated and experimental spectra of mixed valence samples generated by this work. This model can be integrated into a real-time EELS/XAS analysis pipeline on mixtures of copper-containing materials of unknown composition and oxidation state. By expanding the training data, this methodology can be extended to data-driven spectral analysis of a broad range of materials.

36 MATERIALS SCIENCE↗

Determining Reactor Fuel Type from Continuous Antineutrino Monitoring

In this work, we investigate the ability of an antineutrino detector to determine the fuel type of a reactor. A hypothetical 5-ton antineutrino detector is placed 25 m from the core and measures the spectral shape and rate of antineutrinos emitted by fission fragments in the core for a number of 90-d periods. Our results indicate that four major fuel types can be differentiated from the variation of fission fractions over the irradiation time with a true positive probability of detection at approximately 95%. In addition, we demonstrate that antineutrinos can identify the burnup at which weapons-grade mixed-oxide (MOX) fuel would be reduced to reactor-grade MOX, on average, providing assurance that plutonium-disposition goals are met. We also investigate removal scenarios where plutonium is purposefully diverted from a mixture of MOX and low-enriched uranium fuel. Lastly, we discuss how our analysis is impacted by a spectral distortion around 6 MeV observed in the antineutrino spectrum measured from commercial power reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

X-ray Spectral Analysis for Void Fraction Measurement

X-ray attenuation technique has been widely used in two-phase flow measurement. In general, an x-ray system measures the overall attenuation of the two-phase mixture existing between an x-ray source and a detector. One can obtain the void fraction from the intensity of the attenuated x-ray beam and known attenuation coefficients of each phase. Fig. 1 shows a schematic of the line detector system developed at Virginia Tech, and the reconstructed image of a sample slug bubble. The system is capable of obtaining high-resolution (100 {mu}m or less) x-ray images of two-phase flow at a frame rate exceeding 1000 Hz. Note that the intensity of the image can be converted into the line-averaged void fraction data measured by each detector pixel. Typical x-ray systems do not resolve the energy spectrum. Instead, they only measure the total energy deposited in the detector during the exposure period. However, the x-ray source generated from x-ray tubes present polychromatic characteristics. The attenuation coefficient and detector efficiency strongly depend on the energy of x-ray photons. In relation to this, this paper presents the spectral characteristics of a typical x-ray system used for void fraction measurement in two-phase flows. The photon-energy-dependent efficiency of the detector is determined by a calibration test. (authors)

42 ENGINEERING↗

Machine Learning-Based Classification of Lignocellulosic Biomass from Pyrolysis-Molecular Beam Mass Spectrometry Data

High-throughput analysis of biomass is necessary to ensure consistent and uniform feedstocks for agricultural and bioenergy applications and is needed to inform genomics and systems biology models. Pyrolysis followed by mass spectrometry such as molecular beam mass spectrometry (py-MBMS) analyses are becoming increasingly popular for the rapid analysis of biomass cell wall composition and typically require the use of different data analysis tools depending on the need and application. Here, the authors report the py-MBMS analysis of several types of lignocellulosic biomass to gain an understanding of spectral patterns and variation with associated biomass composition and use machine learning approaches to classify, differentiate, and predict biomass types on the basis of py-MBMS spectra. Py-MBMS spectra were also corrected for instrumental variance using generalized linear modeling (GLM) based on the use of select ions relative abundances as spike-in controls. Machine learning classification algorithms e.g., random forest, k-nearest neighbor, decision tree, Gaussian Naïve Bayes, gradient boosting, and multilayer perceptron classifiers were used. The k-nearest neighbors (k-NN) classifier generally performed the best for classifications using raw spectral data, and the decision tree classifier performed the worst. After normalization of spectra to account for instrumental variance, all the classifiers had comparable and generally acceptable performance for predicting the biomass types, although the k-NN and decision tree classifiers were not as accurate for prediction of specific sample types. Gaussian Naïve Bayes (GNB) and extreme gradient boosting (XGB) classifiers performed better than the k-NN and the decision tree classifiers for the prediction of biomass mixtures. The data analysis workflow reported here could be applied and extended for comparison of biomass samples of varying types, species, phenotypes, and/or genotypes or subjected to different treatments, environments, etc. to further elucidate the sources of spectral variance, patterns, and to infer compositional information based on spectral analysis, particularly for analysis of data without a priori knowledge of the feedstock composition or identity.

59 BASIC BIOLOGICAL SCIENCES↗

Chemical variability in mineralized veins observed by ChemCam on the lower slopes of Mount Sharp in Gale crater, Mars

ChemCam has observed a wide range of diagenetic features along the Curiosity rover traverse including pervasive Ca-sulfate veins. Observations by multiple instruments on Curiosity indicate that these veins are hydrated, formed during diagenetic fluid event(s). In this study, we delve into the chemical variability in these Ca-sulfate bearing veins and have identified two subsets in the Murray formation with enrichments in Fe and Fe + Mg. These chemical trends do not reflect a sampling mixture with the surrounding host rock but likely indicates the presence of authigenic phases formed during the emplacement of these veins. Based on passive reflectance spectral analysis and correlation with other elements, Fe 3+ oxides and/or sulfates are proposed to account for the Fe-rich observations in the vicinity of the Naukluft Plateau whereas the Fe + Mg trend is also observed in adjacent dark-toned features with elevated Mn and P near the Old Soaker outcrop. In conclusion, the specific localization of these observations in the Gale stratigraphy implies changing pH and redox conditions in the groundwater at the time of formation of these veins, from oxidizing and likely more acidic near the Naukluft Plateau to more reducing conditions in the upper part of the Murray formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

TH-CD-206-01: Expectation-Maximization Algorithm-Based Tissue Mixture Quantification for Perfusion MRI

Purpose: To investigate the feasibility of estimating the tissue mixture perfusions and quantifying cerebral blood flow change in arterial spin labeled (ASL) perfusion MR images. Methods: The proposed perfusion MR image analysis framework consists of 5 steps: (1) Inhomogeneity correction was performed on the T1- and T2-weighted images, which are available for each studied perfusion MR dataset. (2) We used the publicly available FSL toolbox to strip off the non-brain structures from the T1- and T2-weighted MR images. (3) We applied a multi-spectral tissue-mixture segmentation algorithm on both T1- and T2-structural MR images to roughly estimate the fraction of each tissue type - white matter, grey matter and cerebral spinal fluid inside each image voxel. (4) The distributions of the three tissue types or tissue mixture across the structural image array are down-sampled and mapped onto the ASL voxel array via a co-registration operation. (5) The presented 4-dimensional expectation-maximization (4D-EM) algorithm takes the down-sampled three tissue type distributions on perfusion image data to generate the perfusion mean, variance and percentage images for each tissue type of interest. Results: Experimental results on three volunteer datasets demonstrated that the multi-spectral tissue-mixture segmentation algorithm was effective to initialize tissue mixtures from T1- and T2-weighted MR images. Compared with the conventional ASL image processing toolbox, the proposed 4D-EM algorithm not only generated comparable perfusion mean images, but also produced perfusion variance and percentage images, which the ASL toolbox cannot obtain. It is observed that the perfusion contribution percentages may not be the same as the corresponding tissue mixture volume fractions estimated in the structural images. Conclusion: A specific application to brain ASL images showed that the presented perfusion image analysis method is promising for detecting subtle changes in tissue perfusions, which is valuable for the early diagnosis of certain brain diseases, e.g. multiple sclerosis.

60 APPLIED LIFE SCIENCES↗

Organic contamination detection for isotopic analysis of water by laser spectroscopy

Hydrogen and oxygen stable isotope ratios (δ 2 H, δ 17 O, and δ 18 O values) are commonly used tracers of water. These ratios can be measured by isotope ratio infrared spectroscopy (IRIS). However, IRIS approaches are prone to errors induced by organic compounds present in plant, soil, and natural water samples. A novel approach using 17 O-excess values has shown promise for flagging spectrally contaminated plant samples during IRIS analysis. A systematic assessment of this flagging system is needed to prove it useful. Errors induced by methanol and ethanol water mixtures on measured IRIS and isotope ratio mass spectrometry (IRMS) results were evaluated. For IRIS analyses both liquid- and vapour-mode (via direct vapour equilibration) methods are used. The δ 2 H, δ 17 O, and δ 18 O values were measured and compared with known reference values to determine the errors induced by methanol and ethanol contamination. In addition, the 17 O-excess contamination detection approach was tested. This is a post-processing detection tool for both liquid and vapour IRIS triple-isotope analyses, utilizing calculated 17 O-excess values to flag contaminated samples. Organic contamination induced significant errors in IRIS results, not seen in IRMS results. Methanol caused larger errors than ethanol. Results from vapour-IRIS analyses had larger errors than those from liquid-IRIS analyses. The 17 O-excess approach identified methanol driven error in liquid- and vapour-mode IRIS samples at levels where isotope results became unacceptably erroneous. For ethanol contaminated samples, a mix of erroneous and correct flagging occurred with the 17 O-excess method. Here our results indicate that methanol is the more problematic contaminant for data corruption. The 17 O-excess method was therefore useful for data quality control. Organic contamination caused significant errors in IRIS stable isotope results. These errors were larger during vapour analyses than during liquid IRIS analyses, and larger for methanol than ethanol contamination. The 17 O-excess method is highly sensitive for detecting narrowband (methanol) contamination error in vapour and liquid analysis modes in IRIS.

60 APPLIED LIFE SCIENCES↗