Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multivariate data”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Defining Golden Batches in Biomanufacturing Processes From Internal Metabolic Activity to Detect Process Changes That May Affect Product Quality

ABSTRACT Cellular metabolism plays a role in the observed variability of a drug substance's Critical Quality Attributes (CQAs) made by biomanufacturing processes. Therefore, here we describe a new approach for monitoring biomanufacturing processes that measures a set of metabolic reaction rates (named Critical Metabolic Parameters (CMP) in addition to the macroscopic process conditions currently being used as Critical Process Parameters (CPP) for biomanufacturing. Constraint‐based systems biology models like Flux Balance Analysis (FBA) are used to estimate metabolic reaction rates, and metabolic rates are used as inputs for multivariate Batch Evolution Models (BEM). Metabolic activity was reproducible among batches and could be monitored to detect a deliberately induced macroscopic process shift (i.e., temperature change). The CMP approach has the potential to enable “golden batches” in biomanufacturing processes to be defined from the internal metabolic activity and to aid in detecting process changes that may impact the quality of the product. Overall, the data suggested that monitoring of metabolic activity has promise for biomanufacturing process control.

Biotechnology & Applied Microbiology↗
DOE OSTI2022

Metabolome patterns identify active dechlorination in bioaugmentation consortium SDC-9™

Ultra-high performance liquid chromatography–high-resolution mass spectrometry (UPHLC–HRMS) is used to discover and monitor single or sets of biomarkers informing about metabolic processes of interest. The technique can detect 1000’s of molecules (i.e., metabolites) in a single instrument run and provide a measurement of the global metabolome, which could be a fingerprint of activity. Despite the power of this approach, technical challenges have hindered the effective use of metabolomics to interrogate microbial communities implicated in the removal of priority contaminants. Herein, our efforts to circumvent these challenges and apply this emerging systems biology technique to microbiomes relevant for contaminant biodegradation will be discussed. Chlorinated ethenes impact many contaminated sites, and detoxification can be achieved by organohalide-respiring bacteria, a process currently assessed by quantitative gene-centric tools (e.g., quantitative PCR). This laboratory study monitored the metabolome of the SDC-9™ bioaugmentation consortium during cis-1,2-dichloroethene (cDCE) conversion to vinyl chloride (VC) and nontoxic ethene. Untargeted metabolomics using an UHPLC-Orbitrap mass spectrometer and performed on SDC-9™ cultures at different stages of the reductive dechlorination process detected ~10,000 spectral features per sample arising from water-soluble molecules with both known and unknown structures. Multivariate statistical techniques including partial least squares-discriminate analysis (PLSDA) identified patterns of measurable spectral features (peak patterns) that correlated with dechlorination (in)activity, and ANOVA analyses identified 18 potential biomarkers for this process. Statistical clustering of samples with these 18 features identified dechlorination activity more reliably than clustering of samples based only on chlorinated ethene concentration and Dhc 16S rRNA gene abundance data, highlighting the potential value of metabolomic workflows as an innovative site assessment and bioremediation monitoring tool.

environmental monitoring↗
DOE OSTI2023

Kernel-based global sensitivity analysis obtained from a single data set

Results from global sensitivity analysis (GSA) often guide the understanding of complicated input–output systems. Kernel-based GSA methods have recently been proposed for their capability of treating a broad scope of complex systems. In this paper, we develop a new set of kernel GSA tools when only a single set of input–output data is available. Three key advances are made: (1) A new numerical estimator is proposed that demonstrates an empirical improvement over previous procedures. (2) A computational method for generating inner statistical functions from a single data set is presented. (3) A theoretical extension is made to define conditional sensitivity indices, which reveal the degree that the inputs carry shared information about the output when inherent input–input correlations are present. Utilizing these conditional sensitivity indices, a decomposition is derived for the output uncertainty based on what is called the optimal learning sequence of the input variables, which remains consistent when correlations exist between the input variables. Further, while these advances cover a range of GSA subjects, a common single data set numerical solution is provided by a technique known as the conditional mean embedding of distributions. The new methodology is implemented on benchmark systems to demonstrate the provided insights.

42 ENGINEERING↗
DOE OSTI2022

Development of a prediction model for radiotherapy response among patients with head and neck squamous cell carcinoma based on the tumor immune microenvironment and hypoxia signature

Abstract Introduction The immune system and hypoxia are major factors influencing radiosensitivity in patients with different cancer types. This study aimed at developing a model to predict radiotherapy response in patients with head and neck squamous cell carcinoma (HNSCC) based on the tumor immune microenvironment and hypoxia signature. Materials and Methods We first evaluated the hypoxia status and tumor immune microenvironment in the Cancer Genome Atlas (TCGA) cohort by using transcriptomic data. Differentially expressed genes (DEGs) were identified between the “high immunity and low hypoxia” and “low immunity and high hypoxia” groups and those DEGs significantly associated with disease‐specific survival in the univariate Cox regression analysis were selected as the prognostic DEGs. We selected the immune hypoxia–related genes (IHRGs) by intersecting prognostic DEGs with immune and hypoxia gene sets. We used the IHRGs to train a multivariate Cox regression model in the TCGA cohort, based on which we calculated the IHRG prognostic index (IHRGPI) for each patient and validated its efficacy in predicting radiotherapy response in the Gene Expression Omnibus cohorts. Furthermore, we explored potential mechanisms and effective combinational treatment strategies for different IHRGPI groups. Results Five IHRGs were used to construct the IHRGPI, which was used to dichotomize the cohorts. The patients with lower IHRGPI showed a better radiotherapy response across different cohorts and endpoints, including overall survival, progression‐free survival, and recurrence‐free survival ( p < 0.05). Patients with higher IHRGPI showed greater hypoxia and lesser immune cell infiltration. A lower IHRGPI indicated a better immunotherapy response, while a higher IHRGPI indicated a better chemotherapy response. Conclusions IHRGPI is promising for predicting radiotherapy response and guiding combinational treatment strategies in patients with HNSCC.

Zhu, Guang‐Li↗
DOE OSTI2023

Bullying among children with heart conditions, National Survey of Children’s Health, 2018–2020

Abstract Children with chronic illnesses report being bullied by peers, yet little is known about bullying among children with heart conditions. Using 2018–2020 National Survey of Children’s Health data, the prevalence and frequency of being bullied in the past year (never; annually or monthly; weekly or daily) were compared between children aged 6–17 years with and without heart conditions. Among children with heart conditions, associations between demographic and health characteristics and being bullied, and prevalence of diagnosed anxiety or depression by bullying status were examined. Differences were assessed with chi-square tests and multivariable logistic regression using predicted marginals to produce adjusted prevalence ratios and 95% confidence intervals. Weights yielded national estimates. Of 69,428 children, 2.2% had heart conditions. Children with heart conditions, compared to those without, were more likely to be bullied (56.3% and 43.3% respectively; adjusted prevalence ratio [95% confidence interval] = 1.3 [1.2, 1.4]) and bullied more frequently (weekly or daily = 11.2% and 5.3%; p < 0.001). Among children with heart conditions, characteristics associated with greater odds of weekly or daily bullying included ages 9–11 years compared to 15–17 years (3.4 [2.0, 5.7]), other genetic or inherited condition (1.7 [1.0, 3.0]), ever overweight (1.7 [1.0, 2.8]), and a functional limitation (4.8 [2.7, 8.5]). Children with heart conditions who were bullied, compared to never, more commonly had anxiety (40.1%, 25.9%, and 12.8%, respectively) and depression (18.0%, 9.3%, and 4.7%; p < 0.01 for both). Findings highlight the social and psychological needs of children with heart conditions.

Cardiovascular System & Cardiology↗
NASA NTRS1975

The saturated activities of Na-22, Mn-54, and Co-56 and the depth of sampling of soils

A multivariate linear regression study shows that the saturated activities of Na-22, Co-56, and Mn-54 are functions of the chemical composition and the reciprocal of the cube root of the sample weight. A procedure for taking the temporal variation in these short-lived radionuclides into account is described. The regression results are used to estimate the average depths of sampling of soil material. The data indicate that all 28 samples were in the same conditions of exposure when recovered as they had been for the last approximately 10 years. The saturation activity of Mn-54 is a very weak function of time.

Keith, J. E.↗
NASA NTRS2019

Measurement of CH3D on Titan at Submillimeter Wavelengths

We present the first radio/submillimeter detection of monodeuterated methane (CH3D) in Titan's atmosphere, using archival data from of the Atacama Large Millimeter/submillimeter Array (ALMA). The J (sub K) equals 2 (sub 1) minus 1 (sub 1) and J (sub K) equals 2 (sub 0) minus 1 (sub 0) transitions at 465.235 and 465.250 gigahertz (approximately 0.644 millimeters) were measured at significance levels of 4.6 sigma and 5.7 sigma, respectively. These two lines were modeled using the Non-linear optimal Estimator for MultivariatE spectral analySIS (NEMESIS) radiative transfer code to determine the disk-averaged CH3D volume mixing ratio equals 6.157 times 10 (sup minus 6) in Titan's stratosphere (at altitudes greater than 130 kilometers). By comparison with the CH4 vertical abundance profile measured by Cassini-Huygens mass spectrometry, the resulting value for D/H in CH4 is (1.033 plus or minus 0.081) times 10 (sup minus 4). This is consistent with previous ground-based and in situ measurements from the Cassini-Huygens mission, though slightly lower than the average of the previous values. Additional CH3D observations at higher spatial resolution will be required to determine a value truly comparable with the Cassini-Huygens CH4 measurements, by measuring CH3D with ALMA close to Titan's equator. In the post-Cassini era, spatially resolved observations of CH3D with ALMA will enable the latitudinal distribution of methane to be determined, making this an important molecule for further studies.

Thelen, Alexander E.↗
NASA NTRS1979

Multivariate approximation methods and applications to geophysics and geodesy

The first report in a series is presented which is intended to be written by the author with the purpose of treating a class of approximation methods of functions in one and several variables and ways of applying them to geophysics and geodesy. The first report is divided in three parts and is devoted to the presentation of the mathematical theory and formulas. Various optimal ways of representing functions in one and several variables and the associated error when information is had about the function such as satellite data of different kinds are discussed. The framework chosen is Hilbert spaces. Experiments were performed on satellite altimeter data and on satellite to satellite tracking data.

Munteanu, M. J.↗
DOE OSTI2023

Health outcomes among offspring of US Coast Guard responders to the Deepwater Horizon oil spill, 2010–2011

To evaluate the potential for adverse health outcomes among infants born to US Coast Guard (USCG) responders to the Deepwater Horizon (DWH) oil spill disaster. Department of Defense Birth and Infant Health Research programme data identified a cohort of singleton infants born 2010–2011 to USCG personnel in the DWH Oil Spill Coast Guard Cohort study. Infants were included if their military parent (‘sponsor’) responded to the oil spill during a selected reproductive exposure window (ie, 3 months preconception for male sponsors and periconception through pregnancy for female sponsors), or if their sponsor was a non-responder. χ 2 tests and multivariable log-binomial regression were used to compare the demographic and health characteristics of infants born to spill responders and non-responders. Altogether, 1974 infants with a male sponsor (n=182 responder, n=1792 non-responder) and 628 infants with a female sponsor (n=35 responder, n=593 non-responder) in the DWH Oil Spill Coast Guard Cohort were identified. Health outcomes were similar among the offspring of male responders and non-responders. The frequency of any poor live birth outcome (ie, low birth weight, preterm birth or birth defect) was higher among infants born to female responders (17.1%, n=6) than non-responders (8.9%, n=53); the maternal age-adjusted association was suggestively elevated (risk ratio 1.93, 95% CI 0.89 to 4.16). Infant health outcomes were comparable between the offspring of male USCG oil spill responders and non-responders. Findings were limited by the small number of infants identified, particularly among female responders, and should be interpreted with caution.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1995

A Comparison of Multivariable Control Design Techniques for a Turbofan Engine Control

This paper compares two previously published design procedures for two different multivariable control design techniques for application to a linear engine model of a jet engine. The two multivariable control design techniques compared were the Linear Quadratic Gaussian with Loop Transfer Recovery (LQG/LTR) and the H-Infinity synthesis. The two control design techniques were used with specific previously published design procedures to synthesize controls which would provide equivalent closed loop frequency response for the primary control loops while assuring adequate loop decoupling. The resulting controllers were then reduced in order to minimize the programming and data storage requirements for a typical implementation. The reduced order linear controllers designed by each method were combined with the linear model of an advanced turbofan engine and the system performance was evaluated for the continuous linear system. Included in the performance analysis are the resulting frequency and transient responses as well as actuator usage and rate capability for each design method. The controls were also analyzed for robustness with respect to structured uncertainties in the unmodeled system dynamics. The two controls were then compared for performance capability and hardware implementation issues.

Garg, Sanjay↗
DOE OSTI2025

Uncertainty quantification of a physics-informed model based on sparse identification of a Thermal Energy Distribution System

Integrated energy systems (IES)s are crucial for enhancing the economy and efficiency of power generation sources (e.g., nuclear energy) necessary to unleash American energy dominance. These systems can be integrated with thermal energy storage (TES) and intermittent renewable energies to optimize overall energy use, peak-load regulation, and demand-side responses. However, the stabilization of energy generation, transport, and utilization introduces operational complexities that exceed the challenges of managing each sub-component individually. Currently, though IESs rely on human operators for efficiency and stability, reducing human error risk and enhancing performance through automation is highly desirable. Recent advances at Idaho National Laboratory have demonstrated successful control of the Thermal Energy Distributed System (TEDS). However, the automatic control system depends on a deterministic Sparse Identification of Nonlinear Dynamics with Control (SINDyC) model, which are trained based on simulation data from physics-based simulations. Because of uncertainties in physics-based simulation, SINDyC model results in large discrepancies against experimental data and cannot be reliably used in automatic control. In this paper, we present an innovative approach to address these discrepancies by quantifying uncertainties and developing a more robust model. We first generated trajectories by using first-principles physics codes to encapsulate the experiment. Next, we trained thousands of models by randomly sampling these trajectories. We then collapsed all those models into one probabilistic SINDyC by fitting a multivariate Gaussian distribution onto the resulting coefficient’s distribution. Despite its simplicity, our approach successfully produced 95% confidence intervals that captured the experimental trajectories. It even did so with a higher probability and better U-pooling score across six of the seven relevant quantities of interest (QoIs), as compared to other classical approaches. In conclusion, ongoing research is focusing on generating new experimental trajectories to validate this approach, and on employing Bayesian calibration to refine parametric uncertainties and guide future model development efforts.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Microbial community diversity changes during voltage reversal repair in a 12-unit microbial fuel cell

Microbial fuel cell stacks (MFC-Stack) are often confronted with voltage reversals, likely due to an interplay between microbial community dynamics and insufficient electric circuit balancing. Herein, we provide new insight into voltage reversals by examining the microbiomes of twelve MFC units of a 12-liter Pilot-MFC-Stack during repair. Different biofilm repair methods (self-healing, electrostimulation, and re-acclimatization upon cross-inoculation) were used to evaluate the microbial community response. In addition, MFC-Stack simulation was performed based on Kirchhoff’s Second Law to predict values for source potentials and post-evaluate internal resistances. Analysis of the 16S rRNA amplicon sequencing data suggests that the biofilm repair methods could slowly heal damaged biofilms. Notably, severely voltage reversed MFC units had low electrogen relative abundances (18%) and positive anode potentials, while strong bioanodes and contained more than 50% electrogens and had negative anode potentials. Between-community analyses (beta diversity ordination and multinomial regression) of the voltage reversed MFC units revealed differences among biofilms in contrast to healthy/strong MFC units. Permutational multivariate analysis of variance (PERMANOVA) confirmed that reversed biofilms were, indeed, significantly (p < 0.05) different from stronger ones. Overall, these analyses demonstrated the utility of combining electrotechnical and microbial community analyses, especially beta diversity ordination and multinomial regression, to understand problematic MFC units and the potential success of a biofilm repair method. Finally, thicker biofilms were usually healthier and stronger, although thickness was no guarantee for proper structure and power function as all factors were interdependent. There was an evolutionary trend that strong anodes became stronger/healthier and others weaker. This spontaneous trend has to be considered to avoid irreversible voltage reversals and to repair electrogenic biofilms in an MFC-Stack.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Post-landing major element quantification using SuperCam laser induced breakdown spectroscopy

The SuperCam instrument on the Perseverance Mars 2020 rover uses a pulsed 1064 nm laser to ablate targets at a distance and conduct laser induced breakdown spectroscopy (LIBS) by analyzing the light from the resulting plasma. SuperCam LIBS spectra are preprocessed to remove ambient light, noise, and the continuum signal present in LIBS observations. Prior to quantification, spectra are masked to remove noisier spectrometer regions and spectra are normalized to minimize signal fluctuations and effects of target distance. In some cases, the spectra are also standardized or binned prior to quantification. To determine quantitative elemental compositions of diverse geologic materials at Jezero crater, Mars, we use a suite of 1198 laboratory spectra of 334 well-characterized reference samples. The samples were selected to span a wide range of compositions and include typical silicate rocks, pure minerals (e.g., silicates, sulfates, carbonates, oxides), more unusual compositions (e.g., Mn ore and sodalite), and replicates of the sintered SuperCam calibration targets (SCCTs) onboard the rover. For each major element (SiO 2 , TiO 2 , Al 2 O 3 , FeO T , MgO, CaO, Na 2 O, K 2 O), the database was subdivided into five “folds” with similar distributions of the element of interest. One fold was held out as an independent test set, and the remaining four folds were used to optimize multivariate regression models relating the spectrum to the composition. We considered a variety of models, and selected several for further investigation for each element, based primarily on the root mean squared error of prediction (RMSEP) on the test set, when analyzed at 3 m. In cases with several models of comparable performance at 3 m, we incorporated the SCCT performance at different distances to choose the preferred model. Shortly after landing on Mars and collecting initial spectra of geologic targets, we selected one model per element. Subsequently, with additional data from geologic targets, some models were revised to ensure results that are more consistent with geochemical constraints. The calibration discussed here is a snapshot of an ongoing effort to deliver the most accurate chemical compositions with SuperCam LIBS.

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

Fiducial and differential cross-section measurements of electroweak $W\gamma jj$ production in $pp$ collisions at $\sqrt{s} = 13$ TeV with the ATLAS detector

The observation of the electroweak production of a W boson and a photon in association with two jets, using pp collision data at the Large Hadron Collider at a centre of mass energy of $\sqrt{s} = 13$ TeV, is reported. The data were recorded by the ATLAS experiment from 2015 to 2018 and correspond to an integrated luminosity of 140 fb -1 . This process is sensitive to the quartic gauge boson couplings via the vector boson scattering mechanism and provides a stringent test of the electroweak sector of the Standard Model. Events are selected if they contain one electron or muon, missing transverse momentum, at least one photon, and two jets. Multivariate techniques are used to distinguish the electroweak $W\gamma jj$ process from irreducible background processes. The observed significance of the electroweak $W\gamma jj$ process is well above six standard deviations, compared to an expected significance of 6.3 standard deviations. Fiducial and differential cross sections are measured in a fiducial phase space close to the detector acceptance, which are in reasonable agreement with leading order Standard Model predictions from MADGRAPH5+PYTHIA8 and SHERPA. The results are used to constrain new physics effects in the context of an effective field theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2015

Assessing Wetland Hydroperiod and Soil Moisture With Remote Sensing: A Demonstration for the NASA Plum Brook Station Year 2

Primary Goal: Assist with the evaluation and measuring of wetlands hydroperiod at the PlumBrook Station using multi-source remote sensing data as part of a larger effort on projecting climate change-related impacts on the station's wetland ecosystems. MTRI expanded on the multi-source remote sensing capabilities to help estimate and measure hydroperiod and the relative soil moisture of wetlands at NASA's Plum Brook Station. Multi-source remote sensing capabilities are useful in estimating and measuring hydroperiod and relative soil moisture of wetlands. This is important as a changing regional climate has several potential risks for wetland ecosystem function. The year two analysis built on the first year of the project by acquiring and analyzing remote sensing data for additional dates and types of imagery, combined with focused field work. Five deliverables were planned and completed: 1) Show the relative length of hydroperiod using available remote sensing datasets 2) Date linked table of wetlands extent over time for all feasible non-forested wetlands 3) Utilize LIDAR data to measure topographic height above sea level of all wetlands, wetland to catchment area radio, slope of wetlands, and other useful variables 4) A demonstration of how analyzed results from multiple remote sensing data sources can help with wetlands vulnerability assessment 5) A MTRI style report summarizing year 2 results. This report serves as a descriptive summary of our completion of these our deliverables. Additionally, two formal meetings were held with Larry Liou and Amanda Sprinzl to provide project updates and receive direction on outputs. These were held on 2/26/15 and 9/17/15 at the Plum Brook Station. Principal Component Analysis (PCA) is a multivariate statistical technique used to identify dominant spatial and temporal backscatter signatures. PCA reduces the information contained in the temporal dataset to the first few new Principal Component (PC) images. Some advantages of PCA include the ability to filter out temporal autocorrelation and reduce speckle to the higher order PC images. A PCA was performed using ERDAS Imagine on a time series of PALSAR dates. Hydroperiod maps were created by separating the PALSAR dates into two date ranges, 2006-2008 and 2010, and performing an unsupervised classification on the PCAs.

Brooks, Colin↗
DOE OSTI2024

Search for non-resonant Higgs boson pair production in the $2b+2\ell +{E}_{\textrm{T}}^{\textrm{miss}}$ final state in pp collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

A search for non-resonant Higgs boson pair (HH) production is presented, in which one of the Higgs bosons decays to a b-quark pair ($b\bar{b}$) and the other decays to WW * , ZZ * , or τ + τ – , with in each case a final state with ℓ + ℓ – + neutrinos (ℓ = e, μ). The analysis targets separately the gluon-gluon fusion and vector boson fusion production modes. Data recorded by the ATLAS detector in proton-proton collisions at a centre-of-mass energy of 13 TeV at the Large Hadron Collider, corresponding to an integrated luminosity of 140 fb –1 , are used in this analysis. Events are selected to have exactly two b-tagged jets and two leptons with opposite electric charge and missing transverse momentum in the final state. These events are classified using multivariate analysis algorithms to separate the HH events from other Standard Model processes. No evidence of the signal is found. The observed (expected) upper limit on the cross-section for non-resonant Higgs boson pair production is determined to be 9.7 (16.2) times the Standard Model prediction at 95% confidence level. The Higgs boson self-interaction coupling parameter κ λ and the quadrilinear coupling parameter κ 2V are each separately constrained by this analysis to be within the ranges [–6.2, 13.3] and [–0.17, 2.4], respectively, at 95% confidence level, when all other parameters are fixed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Measurements of $\textrm{t}\overline{\textrm{t}}\textrm{W}$ differential cross sections and the leptonic charge asymmetry at $\sqrt{s}=13$ TeV

Measurements of properties of top quark-antiquark pair production in association with a W boson in proton-proton collisions at a center-of-mass energy of 13 TeV are presented, using a data sample corresponding to an integrated luminosity of 138 fb −1 , recorded by the CMS experiment at the CERN LHC. Events are selected based on the presence of either two leptons with the same electric charge or three leptons, and multiple jets and b-tagged jets. We present measurements of differential production cross sections as a function of kinematic variables sensitive to different aspects of the process modeling, using a multivariate discriminator in the two-lepton selection region and a simple selection-based method in the three-lepton region. The normalized cross section measurements are generally consistent with the standard model expectations, while we observe larger values compared to the expectations in the absolute cross section measurements, consistent with previous inclusive cross section measurements. In addition, we measure the leptonic charge asymmetry of this process, obtaining an observed value of ${A}_c^{\ell }=-{0.19}_{-0.18}^{+0.16}$, consistent with the expectation of −0.085 ± 0.006 predicted by next-to-leading order simulations.

Hadron-Hadron Scattering↗
DOE OSTI2021

Improving ChemCam LIBS long-distance elemental compositions using empirical abundance trends

The ChemCam instrument on the Curiosity rover provides chemical compositions of Martian rocks and soils using remote laser-induced breakdown spectroscopy (LIBS). The elemental calibration is stable as a function of distance for Ti, Fe, Mg, and Ca. The calibration shows small, systematically increasing abundance trends as a function of distance for Al, Na, K, and to some extent, Si. The distance effect is known to be due to a dependence with distance on the relative strengths of atomic transition lines. Emission lines representing transitions from relatively low energy levels remain intense at longer distances while emission lines representing transitions from higher energy levels decrease in intensity more rapidly as a function of distance. The multivariate algorithms used to determine elemental compositions rely on a large number of emission lines in many cases, so rather than trying to correct all emission lines, a study was made of the predicted compositions as a function of distance, in order to determine an empirical correction. Abundance trends can be well approximated by a linear trend with distance within the ranges of abundances and distances observed up to ~6 m. Data from 11 distinct geological members and data groups of the Murray formation in Gale crater, Mars, were used to form the model, selecting the members and data groups yielding the best statistics. The model was tested using data from several targets observed from two different distances, and using data from the Kimberley formation, the composition of which is significantly different from the Murray formation, showing that the model works on other compositions beyond those used to build the model. For long-distance observations up to ~6 m, corrections can be made back to an equivalent composition at the median distance of ChemCam observations (2.6 m). Finally, the model has been validated up to 6.2 m, although ChemCam is able to observe bedrock targets to >7 m, and iron meteorites to distances of >9 m.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗