Utilizing Linear Combination Fitting (LCF) to Model Highly Reduced, Multi-Component Mixtures of Cation-Sulfide Species Relevant to Mercurian Melts
Explore the source record for details and available documents.
B A Konecke↗
SEARCH · Search NASA
Results for “linear combination fitting”
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.
Non-Targeted Effects and the Dose Response for Heavy Ion Tumorigenesis
BACKGROUND: There is no human epidemiology data available to estimate the heavy ion cancer risks experienced by astronauts in space. Studies of tumor induction in mice are a necessary step to estimate risks to astronauts. Previous experimental data can be better utilized to model dose response for heavy ion tumorigenesis and plan future low dose studies. DOSE RESPONSE MODELS: The Harderian Gland data of Alpen et al.[1-3] was re-analyzed [4] using non-linear least square regression. The data set measured the induction of Harderian gland tumors in mice by high-energy protons, helium, neon, iron, niobium and lanthanum with LET s ranging from 0.4 to 950 keV/micron. We were able to strengthen the individual ion models by combining data for all ions into a model that relates both radiation dose and LET for the ion to tumor prevalence. We compared models based on Targeted Effects (TE) to one motivated by Non-targeted Effects (NTE) that included a bystander term that increased tumor induction at low doses non-linearly. When comparing fitted models to the experimental data, we considered the adjusted R2, the Akaike Information Criteria (AIC), and the Bayesian Information Criteria (BIC) to test for Goodness of fit.In the adjusted R2test, the model with the highest R2values provides a better fit to the available data. In the AIC and BIC tests, the model with the smaller values of the summary value provides the better fit. The non-linear NTE models fit the combined data better than the TE models that are linear at low doses. We evaluated the differences in the relative biological effectiveness (RBE) and found the NTE model provides a higher RBE at low dose compared to the TE model. POWER ANALYSIS: The final NTE model estimates were used to simulate example data to consider the design of new experiments to detect NTE at low dose for validation. Power and sample sizes were calculated for a variety of radiation qualities including some not considered in the Harderian Gland data set and with different background tumor incidences. We considered different experimental designs with varying number of doses and varying low doses dependant on the LET of the radiation. The optimal design to detect a NTE for an individual ion had 4 doses equally spaced below a maximal dose where bending due to cell sterilization was < 2%. For example at 100 keV/micron we would irradiate at 0.03 Gy, 0.065 Gy, 0.13 Gy, and 0.26 Gy and require 850 mice including a control dose for a sensitivity to detect NTE with 80% power. Sample sizes could be improved by combining ions similar to the methods used with the Harderian Gland data.
Chappelli, Lori J.↗
Passive microwave remote sensing of rainfall with SSM/I: Algorithm development and implementation
A physically based algorithm sensitive to emission and scattering is used to estimate rainfall using the Special Sensor Microwave/Imager (SSM/I). The algorithm is derived from radiative transfer calculations through an atmospheric cloud model specifying vertical distributions of ice and liquid hydrometeors as a function of rain rate. The algorithm is structured in two parts: SSM/I brightness temperatures are screened to detect rainfall and are then used in rain-rate calculation. The screening process distinguishes between nonraining background conditions and emission and scattering associated with hydrometeors. Thermometric temperature and polarization thresholds determined from the radiative transfer calculations are used to detect rain, whereas the rain-rate calculation is based on a linear function fit to a linear combination of channels. Separate calculations for ocean and land account for different background conditions. The rain-rate calculation is constructed to respond to both emission and scattering, reduce extraneous atmospheric and surface effects, and to correct for beam filling. The resulting SSM/I rain-rate estimates are compared to three precipitation radars as well as to a dynamically simulated rainfall event. Global estimates from the SSM/I algorithm are also compared to continental and shipboard measurements over a 4-month period. The algorithm is found to accurately describe both localized instantaneous rainfall events and global monthly patterns over both land and ovean. Over land the 4-month mean difference between SSM/I and the Global Precipitation Climatology Center continental rain gauge database is less than 10%. Over the ocean, the mean difference between SSM/I and the Legates and Willmott global shipboard rain gauge climatology is less than 20%.
Ferriday, James G.↗
Constraining the Metallicities, Ages, Star Formation Histories, and Ionizing Continua of Extragalactic Massive Star Populations
We infer the properties of massive star populations using the far-ultraviolet stellar continua of 61 star-forming galaxies: 42 at low redshift observed with the Hubble Space Telescope and 19 at z ~ 2 from the MegaSaura sample. We fit each stellar continuum with a linear combination of up to 50 single-age and single-metallicity starburst99 models. From these fits, we derive light-weighted ages and metallicities, which agree with stellar wind and photospheric spectral features, and infer the spectral shapes and strengths of the ionizing continua. Inferred light-weighted stellar metallicities span 0.05–1.5 Z(sub ⊙) and are similar to the measured nebular metallicities. We quantify the ionizing continua using the ratio of the ionizing flux at 900 Å to the non-ionizing flux at 1500 Å and demonstrate the evolution of this ratio with stellar age and metallicity using theoretical single-burst models. These single-burst models only match the inferred ionizing continua of half of the sample, while the other half are described by a mixture of stellar ages. Mixed-age populations produce stronger and harder ionizing spectra than continuous star formation histories, but, contrary to previous studies that assume constant star formation, have similar stellar and nebular metallicities. Stellar population age and metallicity affect the far-UV continua in different and distinguishable ways; assuming a constant star formation history diminishes the diagnostic power. Finally, we provide simple prescriptions to determine the ionizing photon production efficiency (ξ(sub ion)) from the stellar population properties. The ξ(sub ion) inferred from the observed star-forming galaxies has a range of log(ξ(sub ion)) = 24.4–25.7 Hz erg(exp −1) that depends on the stellar population age, metallicity, star formation history, and contributions from binary star evolution. These stellar population properties must be observationally determined to accurately determine the number of ionizing photons generated by massive stars.
Chisholm, J.↗
Multipath calibration in GPS pseudorange measurements
Novel techniques are disclosed for eliminating multipath errors, including mean bias errors, in pseudorange measurements made by conventional global positioning system receivers. By correlating the multipath signals of different satellites at their cross-over points in the sky, multipath mean bias errors are effectively eliminated. By then taking advantage of the geometrical dependence of multipath, a linear combination of spherical harmonics are fit to the satellite multipath data to create a hemispherical model of the multipath. This calibration model can then be used to compensate for multipath in subsequent measurements and thereby obtain GPS positioning to centimeter accuracy.
Kee, Changdon↗
Improved absorption cross-sections of oxygen in the wavelength region 205-240 nm of the Herzberg continuum
The laboratory values of the Herzberg continuum absorption cross-section of oxygen at room temperature from Cheung et al. (1986) and Jenouvrier et al. (1986) are compared and analyzed. It is found that there is no discrepancy between the absolute values of these two sets of independent measurements. The values are combined in a linear least-squares fit to obtain improved values of the Herzberg continuum cross-section of oxygen at room temperature throughout the wavelength region 205-240 nm. The results are compared with in situ and other laboratory measurements.
Yoshino, K.↗
Improving Photometry and Stellar Signal Preservation with Pixel-Level Systematic Error Correction
The Kepler Mission has demonstrated that excellent stellar photometric performance can be achieved using apertures constructed from optimally selected CCD pixels. The clever methods used to correct for systematic errors, while very successful, still have some limitations in their ability to extract long-term trends in stellar flux. They also leave poorly correlated bias sources, such as drifting moiré pattern, uncorrected. We will illustrate several approaches where applying systematic error correction algorithms to the pixel time series, rather than the co-added raw flux time series, provide significant advantages. Examples include, spatially localized determination of time varying moiré pattern biases, greater sensitivity to radiation-induced pixel sensitivity drops (SPSDs), improved precision of co-trending basis vectors (CBV), and a means of distinguishing the stellar variability from co-trending terms even when they are correlated. For the last item, the approach enables physical interpretation of appropriately scaled coefficients derived in the fit of pixel time series to the CBV as linear combinations of various spatial derivatives of the pixel response function (PRF). We demonstrate that the residuals of a fit of soderived pixel coefficients to various PRF-related components can be deterministically interpreted in terms of physically meaningful quantities, such as the component of the stellar flux time series which is correlated with the CBV, as well as, relative pixel gain, proper motion and parallax. The approach also enables us to parameterize and assess the limiting factors in the uncertainties in these quantities.
Kolodzijczak, Jeffrey J.↗
Spline approximation of quantile functions
The study reported here explored the development and utility of a spline representation of the sample quantile function of a continuous probability distribution in providing a functional description of a random sample and a method of generating random variables. With a spline representation, the random samples are generated by transforming a sample of uniform random variables to the interval of interest. This is useful, for example, in simulation studies in which a random sample represents the only known information about the distribution. The spline formulation considered here consists of a linear combination of cubic basis splines (B-splines) fit in a least squares sense to the sample quantile function using equally spaced knots. The following discussion is presented in five parts. The first section highlights major results realized from the study. The second section further details the results obtained. The methodology used is described in the third section, followed by a brief discussion of previous research on quantile functions. Finally, the results of the study are evaluated.
Schiess, J. R.↗
Elemental abundances of flaring solar plasma - Enhanced neon and sulfur
Elemental abundances of two flares observed with the SMM Flat Crystal Spectrometer are compared and contrasted. The first had a gradual rise and a slow decay, while the second was much more impulsive. Simultaneous spectra of seven bright soft X-ray resonance lines provide information over a broad temperature range and are available throughout both flares, making these events unique in the SMM data base. For the first flare, the plasma seemed to be characterized by coronal abundances but, for the second, the plasma composition could not be coronal, photospheric, or a linear combination of both. A good differential emission measure fit required enhanced neon such that Ne/O = 0.32 +/- 0.02, a value which is inconsistent with the current models of coronal abundances based on the elemental first-ionization potential. Similar values of enhanced neon are found for flaring plasma observed by the SMM gamma-ray spectrometer, in (He-3)-rich solar energetic particle events, and in the decay phase of several long duration soft X-ray events. Sulfur is also enhanced in the impulsive flare, but not as dramatically as neon. These events are compared with two models which attempt to explain the enhanced values of neon and sulfur.
Schmelz, J. T.↗
An analysis of fiber-matrix interface failure stresses for a range of ply stress states
A graphite/bismaleimide laminate was prepared without the usual fiber treatment and was tested over a wide range of stress states to measure its ply cracking strength. These tests were conducted using off-axis flexure specimens and produced fiber-matrix interface failure data over a correspondingly wide range of interface stress states. The absence of fiber treatment, weakened the fiber-matrix interfaces and allowed these tests to be conducted at load levels that did not yield the matrix. An elastic micromechanics computer code was used to calculate the fiber-matrix interface stresses at failure. Two different fiber-array models (square and diamond) were used in these calculations to analyze the effects of fiber arrangement as well as stress state on the critical interface stresses at failure. This study showed that both fiber-array models were needed to analyze interface stresses over the range of stress states. A linear equation provided a close fit to these critical stress combinations and, thereby, provided a fiber-matrix interface failure criterion. These results suggest that prediction procedures for laminate ply cracking can be based on micromechanics stress analyses and appropriate fiber-matrix interface failure criteria. However, typical structural laminates may require elastoplastic stress analysis procedures that account for matrix yielding, especially for shear-dominated ply stress states.
Crews, J. H.↗
Highly Scalable Matching Pursuit Signal Decomposition Algorithm
Matching Pursuit Decomposition (MPD) is a powerful iterative algorithm for signal decomposition and feature extraction. MPD decomposes any signal into linear combinations of its dictionary elements or atoms . A best fit atom from an arbitrarily defined dictionary is determined through cross-correlation. The selected atom is subtracted from the signal and this procedure is repeated on the residual in the subsequent iterations until a stopping criterion is met. The reconstructed signal reveals the waveform structure of the original signal. However, a sufficiently large dictionary is required for an accurate reconstruction; this in return increases the computational burden of the algorithm, thus limiting its applicability and level of adoption. The purpose of this research is to improve the scalability and performance of the classical MPD algorithm. Correlation thresholds were defined to prune insignificant atoms from the dictionary. The Coarse-Fine Grids and Multiple Atom Extraction techniques were proposed to decrease the computational burden of the algorithm. The Coarse-Fine Grids method enabled the approximation and refinement of the parameters for the best fit atom. The ability to extract multiple atoms within a single iteration enhanced the effectiveness and efficiency of each iteration. These improvements were implemented to produce an improved Matching Pursuit Decomposition algorithm entitled MPD++. Disparate signal decomposition applications may require a particular emphasis of accuracy or computational efficiency. The prominence of the key signal features required for the proper signal classification dictates the level of accuracy necessary in the decomposition. The MPD++ algorithm may be easily adapted to accommodate the imposed requirements. Certain feature extraction applications may require rapid signal decomposition. The full potential of MPD++ may be utilized to produce incredible performance gains while extracting only slightly less energy than the standard algorithm. When the utmost accuracy must be achieved, the modified algorithm extracts atoms more conservatively but still exhibits computational gains over classical MPD. The MPD++ algorithm was demonstrated using an over-complete dictionary on real life data. Computational times were reduced by factors of 1.9 and 44 for the emphases of accuracy and performance, respectively. The modified algorithm extracted similar amounts of energy compared to classical MPD. The degree of the improvement in computational time depends on the complexity of the data, the initialization parameters, and the breadth of the dictionary. The results of the research confirm that the three modifications successfully improved the scalability and computational efficiency of the MPD algorithm. Correlation Thresholding decreased the time complexity by reducing the dictionary size. Multiple Atom Extraction also reduced the time complexity by decreasing the number of iterations required for a stopping criterion to be reached. The Course-Fine Grids technique enabled complicated atoms with numerous variable parameters to be effectively represented in the dictionary. Due to the nature of the three proposed modifications, they are capable of being stacked and have cumulative effects on the reduction of the time complexity.
Christensen, Daniel↗
Grain size of the surface regolith of asteroid 4 Vesta estimated from its reflectance spectrum in comparison with HED meteorites
The grain-size distribution of the regolith of asteroid 4 Vesta has been estimated by comparing its reflectance spectra (0.3-2.6 microns) with those of HED meteorites. The finest grain-size separate (less than 25 micrometers) of a particular howardite has a reflectance spectrum most similar to Vesta's. In order to better simulate Vesta's surface mineralogy, reflectance spectra of those finest HED meteorite powders were linearly combined, and Vesta's spectrum was scaled for the best fit between them. Both the albedo and the shape of reflectance spectrum of Vesta were well reproduced by regional mixtures of the finest (less than 25 micrometers) powders of HED meteorites. The result suggests the heterogeneity of Vesta's surface and provides an estimate of the visible reflectance of Vesta that is close to its Infrared Astronomical Satellite (IRAS) albedo. Thus, this suggests that fine grains can be generated and retained by relatively small bodies (Vesta is approximately 500 km in diameter).
Hiroi, Takahiro↗
Quantifying radiation quality for space relevant radiation types: Fitting excess risk models to three combined HZE-irradiated mouse datasets
Radiation health risks are predominantly derived from low linear energy transfer (LET) terrestrial exposures; however, space radiation includes exposure to high-LET and high-charge, high-energy (HZE) particles. Accurately quantifying the differences in radiation quality between the space and terrestrial radiation environments is important for assessing and predicting health risks for astronauts. Weil et al. 2009 and 2014 used two different inbred mouse strains to study differences in hepatocellular carcinoma (HCC) tumorigenesis after exposures to low- and high- LET radiation. More recently, Edmundson et al. 2020 provided valuable new tumor data in outbred mice that were exposed to low- and high-LET radiation. The present study aims to rigorously investigate a relative biological effectiveness (RBE) factor by leveraging the HCC tumor data from Weil et al. 2009, Weil et al. 2014, and Edmundson et al. 2020. The three experiments were similarly designed, allowing the raw data to be combined into a pooled dataset to estimate excess relative risk (ERR) and excess absolute risk (EAR) models using Bayesian Poisson regression. These effect estimates from the pooled data provide greater power to calculate a data driven RBE. Extensive sensitivity analyses test the robustness of RBE estimates to various model assumptions. The following questions will be explored through the sensitivity analyses: • Is the shape of the dose response different for low-LET radiation and HZE radiation, indicating that RBE is a function of dose? • Does attained age modify the effect estimates differently for low-LET radiation and HZE radiation, indicating RBE is a function of attained age? • Are the effect estimates and RBE estimates different for inbred mouse strains and outbred mouse strains? • Do assumptions about differences in ERR models and EAR models change the estimated RBE? Additional studies would be needed to validate the findings from these exploratory analyses.
Lori J. Chappell↗
Quantifying radiation quality for space relevant radiation types: Fitting excess risk models to three combined HZE-irradiated mouse datasets
Radiation health risks are predominantly derived from low linear energy transfer (LET) terrestrial exposures; however, space radiation includes exposure to high-LET and high-charge, high-energy (HZE) particles. Accurately quantifying the differences in radiation quality between the space and terrestrial radiation environments is important for assessing and predicting health risks for astronauts. Weil et al. 2009 and 2014 used two different inbred mouse strains to study differences in hepatocellular carcinoma (HCC) tumorigenesis after exposures to low- and high- LET radiation[1-2]. More recently, Edmondson et al. 2020 provided valuable new tumor data in outbred mice that were exposed to low- and high-LET radiation[3]. The present study aims to rigorously investigate a relative biological effectiveness (RBE) factor by leveraging the HCC tumor data from the three datasets[1-3]. The three experiments were similarly designed, allowing the raw data to be combined into a pooled dataset to estimate excess relative risk (ERR) and excess absolute risk (EAR) models using Bayesian Poisson regression.
Lori J. Chappell↗
Integrating Solar Induced Fluorescence and the Photochemical Reflectance Index for Estimating Gross Primary Production in a Cornfield
The utilization of remotely sensed observations for light use efficiency (LUE) and tower-based gross primary production (GPP) estimates was studied in a USDA cornfield. Nadir hyperspectral reflectance measurements were acquired at canopy level during a collaborative field campaign conducted in four growing seasons. The Photochemical Reflectance Index (PRI) and solar induced chlorophyll fluorescence (SIF), were derived. SIF retrievals were accomplished in the two telluric atmospheric oxygen absorption features centered at 688 nm (O2-B) and 760 nm (O2-A). The PRI and SIF were examined in conjunction with GPP and LUE determined by flux tower-based measurements. All of these fluxes, environmental variables, and the PRI and SIF exhibited diurnal as well as day-to-day dynamics across the four growing seasons. Consistent with previous studies, the PRI was shown to be related to LUE (r2 = 0.54 with a logarithm fit), but the relationship varied each year. By combining the PRI and SIF in a linear regression model, stronger performances for GPP estimation were obtained. The strongest relationship (r2 = 0.80, RMSE = 0.186 mg CO2/m2/s) was achieved when using the PRI and SIF retrievals at 688 nm. Cross-validation approaches were utilized to demonstrate the robustness and consistency of the performance. This study highlights a GPP retrieval method based entirely on hyperspectral remote sensing observations.
cornfield↗
North America-Pacific plate boundary, an elastic-plastic megashear - Evidence from very long baseline interferometry
Data obtained by Mark III VLBI measurements of radio signals from permanent and mobile VLBI sites for 5.5 years of observations, starting in October 1982, were used to derive a picture of the earth crust deformation near the North America-Pacific plate boundary. The data, which included the vector positions of the VLBI sites and their rate of change, were used for comparison with a number of lithospheric deformation models based upon the concept that the motions of points near the North America-Pacific plate boundary are a linear combination of North America and Pacific velocities. The best of these models were found to fit 95 percent of the variance in 139 VLBI length and transverse velocity observations. Instantaneous shear deformation associated with plate tectonics is apparently developing in a zone 450 km wide paralleling the San Andreas Fault; some of this deformation will be recovered through elastic rebound, while the rest will be permanently set through plastic processes. Because the VLBI data have not been collected for a significant fraction of the earthquake cycle, they cannot discriminate between elastic and plastic behaviors.
Ward, Steven N.↗
ASpec: A New Spectrum and Line Analysis Package for IRAF
ASpec is a new spectrum and line analysis package developed at STScI. ASpec is designed as an add-on package for IRAF and incorporates a variety of analysis techniques for astronomical spectra. ASpec operates on spectra from a wide variety of ground-based and space-based instruments. ASpec allows simultaneous handling of spectra from different wavelength regimes. ASpec accommodates non-linear dispersion relations. ASpec provides a variety of functions, individually or in combination, with which to fit spectral features and the continuum. ASpec permits the masking of known bad data. ASpec provides a powerful, intuitive graphical user interface implemented using the new IRAF Object Manager and customized to handle: data input/output (I/O); on-line 'help'; selection of relevant features for analysis; plotting and graphical interaction; and data base management.
Hulbert, Stephen J.↗
Analysis of a Hybrid Wing Body Center Section Test Article
The hybrid wing body center section test article is an all-composite structure made of crown, floor, keel, bulkhead, and rib panels utilizing the Pultruded Rod Stitched Efficient Unitized Structure (PRSEUS) design concept. The primary goal of this test article is to prove that PRSEUS components are capable of carrying combined loads that are representative of a hybrid wing body pressure cabin design regime. This paper summarizes the analytical approach, analysis results, and failure predictions of the test article. A global finite element model of composite panels, metallic fittings, mechanical fasteners, and the Combined Loads Test System (COLTS) test fixture was used to conduct linear structural strength and stability analyses to validate the specimen under the most critical combination of bending and pressure loading conditions found in the hybrid wing body pressure cabin. Local detail analyses were also performed at locations with high stress concentrations, at Tee-cap noodle interfaces with surrounding laminates, and at fastener locations with high bearing/bypass loads. Failure predictions for different composite and metallic failure modes were made, and nonlinear analyses were also performed to study the structural response of the test article under combined bending and pressure loading. This large-scale specimen test will be conducted at the COLTS facility at the NASA Langley Research Center.
Wu, Hsi-Yung T.↗