Development of Graphical User Interface and Spectral Fitting Tool for NEQAIR
Create a GUI and spectral fitting tool for NEQAIR
spectral fitting↗
SEARCH · Search NASA
Results for “spectral 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.
XSPEC: An x ray spectral fitting package. Version 2 of the user's guide
The user guide for XSPEC, a command driven, interactive, X-ray spectral fitting program is presented. It is designed to be completely detector independent so it can be used for any X-ray spectral instrument. An overview of the program commands and a walk through of an XSPEC session is presented. Individual commands and descriptions of the spectral models are given. For the more experienced user details of some of the programs associated with XSPEC, the command parser, and the addition of models to XSPEC are included. The PLT plotting package used by XSPEC is described together with the maximum likelihood option for users possessing X-ray spectra with small numbers of counts per bin. Details on the various files used by XSPEC and on the general file structure are also given.
Shafer, R. A.↗
Use of Machine Learning and Principal Component Analysis to Retrieve Nitrogen Dioxide (NO 2 ) With Hyperspectral Imagers and Reduce Noise in Spectral Fitting
Nitrogen dioxide (NO 2 ) is an important trace-gas pollutant and climate agent whose presence also leads to spectral interference in ocean color retrievals. NO 2 column densities have been retrieved with satellite UV–Vis spectrometers such as the Ozone Monitoring Instrument (OMI) and the Tropospheric Monitoring Instrument (TROPOMI) that typically have spectral resolutions of the order of 0.5 nm or better and spatial footprints as small as 3.6 km × 5.6 km. These NO 2 observations are used to estimate emissions, monitor pollution trends, and study effects on human health. Here, we investigate whether it is possible to retrieve NO 2 amounts with lower-spectral-resolution hyperspectral imagers such as the Ocean Color Instrument (OCI) that will fly on the Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) satellite set for launch in early 2024. OCI will have a spectral resolution of 5 nm and a spatial resolution of ∼ 1 km with global coverage in 1–2 d. At this spectral resolution, small-scale spectral structure from NO 2 absorption is still present. We use real spectra from the OMI to simulate OCI spectra that are in turn used to estimate NO 2 slant column densities (SCDs) with an artificial neural network (NN) trained on target OMI retrievals. While we obtain good results with no noise added to the OCI simulated spectra, we find that the expected instrumental noise substantially degrades the OCI NO 2 retrievals. Nevertheless, the NO 2 information from OCI may be of value for ocean color retrievals. OCI retrievals can also be temporally averaged over timescales of the order of months to reduce noise and provide higher-spatial-resolution maps that may be useful for downscaling lower-spatial-resolution data provided by instruments such as OMI and TROPOMI; this downscaling could potentially enable higher-resolution emissions estimates and be useful for other applications. In addition, we show that NNs that use coefficients of leading modes of a principal component analysis of radiance spectra as inputs appear to enable noise reduction in NO 2 retrievals. Once trained, NNs can also substantially speed up NO 2 spectral fitting algorithms as applied to OMI, TROPOMI, and similar instruments that are flying or will soon fly in geostationary orbit.
NO2↗
FIGS: spectral fitting constraints on the star formation history of massive galaxies since the cosmic noon
We constrain the stellar population properties of a sample of 52 massive galaxies – with stellar mass log (M(s)/M(⊙)) ≳ 10.5 – over the redshift range 0.5 < z < 2 by use of observer-frame optical and near-infrared slitless spectra from Hubble Space Telescope’s ACS and WFC3 grisms. The deep exposures (∼100 ks) allow us to target individual spectra of massive galaxies to F160W = 22.5 AB. Our spectral fitting approach uses a set of six base models adapted to the redshift and spectral resolution of each observation, and fits the weights of the base models, including potential dust attenuation, via a Markov Chain Monte Carlo method. Our sample comprises a mixed distribution of quiescent (19) and star-forming galaxies (33). We quantify the width of the age distribution (Δt) that is found to dominate the variance of the retrieved parameters according to principal component analysis. The population parameters follow the expected trend towards older ages with increasing mass, and Δt appears to weakly anticorrelate with stellar mass, suggesting a more efficient star formation at the massive end. As expected, the redshift dependence of the relative stellar age (measured in units of the age of the Universe at the source) in the quiescent sample rejects the hypothesis of a single burst (aka monolithic collapse). Radial colour gradients within each galaxy are also explored, finding a wider scatter in the star-forming subsample, but no conclusive trend with respect to the population parameters.
Ignacio Ferreras↗
Intensity Conserving Spectral Fitting
The detailed shapes of spectral line profiles provide valuable information about the emitting plasma, especially when the plasma contains an unresolved mixture of velocities, temperatures, and densities. As a result of finite spectral resolution, the intensity measured by a spectrometer is the average intensity across a wavelength bin of non-zero size. It is assigned to the wavelength position at the center of the bin. However, the actual intensity at that discrete position will be different if the profile is curved, as it invariably is. Standard fitting routines (spline, Gaussian, etc.) do not account for this difference, and this can result in significant errors when making sensitive measurements. Detection of asymmetries in solar coronal emission lines is one example. Removal of line blends is another. We have developed an iterative procedure that corrects for this effect. It can be used with any fitting function, but we employ a cubic spline in a new analysis routine called Intensity Conserving Spline Interpolation (ICSI). As the name implies, it conserves the observed intensity within each wavelength bin, which ordinary fits do not. Given the rapid convergence, speed of computation, and ease of use, we suggest that ICSI be made a standard component of the processing pipeline for spectroscopic data.
Klimchuk, J. A.↗
A Cloud Algorithm Based on the O2-O2 477 nm Absorption Band Featuring an Advanced Spectral Fitting Method and the Use of Surface Geometry-Dependent Lambertian-Equivalent Reflectivity
We discuss a new cloud algorithm that retrieves an effective cloud pressure, also known as cloud optical centroid pressure (OCP), from oxygen dimer (O2-O2) absorption at 477nm after determining an effective cloud fraction (ECF) at 466nm, a wavelength not significantly affected by trace-gas absorption and rotational Raman scattering. The retrieved cloud products are intended for use as inputs to the operational nitrogen dioxide (NO2) retrieval algorithm for the Ozone Monitoring Instrument (OMI) flying on the Aura satellite. The cloud algorithm uses temperature-dependent O2-O2 cross sections and incorporates flexible spectral fitting techniques that account for specifics of the surface reflectivity. The fitting procedure derives O2-O2 slant column densities (SCDs) from radiances after O3, NO2, and H2O absorption features have been removed based on estimates of the amounts of these species from independent OMI algorithms. The cloud algorithm is based on the frequently used mixed Lambertian-equivalent reflectivity (MLER) concept. A geometry-dependent Lambertian-equivalent reflectivity (GLER), which is a proxy of surface bidirectional reflectance, is used for the ground reflectivity in our implementation of the MLER approach. The OCP is derived from a match of the measured O2-O2 SCD to that calculated with the MLER method. Temperature profiles needed for computation of vertical column densities are taken from the Global Modeling Initiative (GMI) model. We investigate the effect of using GLER instead of climatological LER on the retrieved ECF and OCP. For evaluation purposes, the retrieved ECFs and OCPs are compared with those from the operational OMI cloud product, which is also based on the same O2-O2 absorption band. Impacts of the application of the newly developed cloud algorithm to the OMI NO2 retrieval are discussed.
cloud fraction and pressure↗
The X-ray spectra of galaxies. I - Spectral fits of individual galaxies and X-ray colors
The X-ray spectra of normal galaxies is investigated systematically by using the Imaging Proportional Counter (IPC) data in the Einstein data bank. In addition to the standard model fitting technique, X-ray colors are introduced in order to extract spectral information from the fainter X-ray sources. Spectral parameters for 43 galaxies and X-ray colors for 127 galaxies are presented.
Kim, D.-W.↗
Spectral fitting, shock layer modeling, and production of nitrogen oxides and excited nitrogen
An analysis was made of N2 emission from 8.72 MJ/kg shock layer at 2.54, 1.91, and 1.27 cm positions and vibrational state distributions, temperatures, and relative electronic state populations was obtained from data sets. Other recorded arc jet N2 and air spectral data were reviewed and NO emission characteristics were studied. A review of operational procedures of the DSMC code was made. Information on other appropriate codes and modifications, including ionization, were made as well as a determination of the applicability of codes reviewed to task requirement. A review was also made of computational procedures used in CFD codes of Li and other codes on JSC computers. An analysis was made of problems associated with integration of specific chemical kinetics applicable to task into CFD codes.
Blackwell, H. E.↗
FIGS: spectral fitting constraints on the star formation history of massive galaxies since the cosmic noon
Explore the source record for details and available documents.
Ignacio Ferreras↗
Soft X-ray spectral fits of Geminga with model neutron star atmospheres
The spectrum of the soft X-ray pulsar Geminga consists of two components, a softer one which can be interpreted as thermal-like radiation from the surface of the neutron star, and a harder one interpreted as radiation from a polar cap heated by relativistic particles. We have fitted the soft spectrum using a detailed magnetized hydrogen atmosphere model. The fitting parameters are the hydrogen column density, the effective temperature T(sub eff), the gravitational redshift z, and the distance to radius ratio, for different values of the magnetic field B. The best fits for this model are obtained when B less than or approximately 1 x 10(exp 12) G and z lies on the upper boundary of the explored range (z = 0.45). The values of T(sub eff) approximately = (2-3) x 10(exp 5) K are a factor of 2-3 times lower than the value of T(sub eff) obtained for blackbody fits with the same z. The lower T(sub eff) increases the compatibility with some proposed schemes for fast neutrino cooling of neutron stars (NSs) by the direct Urca process or by exotic matter, but conventional cooling cannot be excluded. The hydrogen atmosphere fits also imply a smaller distance to Geminga than that inferred from a blackbody fit. An accurate evaluation of the distance would require a better knowledge of the ROSAT Position Sensitive Proportional Counter (PSPC) response to the low-energy region of the incident spectrum. Our modeling of the soft component with a cooler magnetized atmosphere also implies that the hard-component fit requires a characteristic temperature which is higher (by a factor of approximately 2-3) and a surface area which is smaller (by a factor of 10(exp 3), compared to previous blackbody fits.
Meyer, R. D.↗
CARS Spectral Fitting with Multiple Resonant Species using Sparse Libraries
The dual pump CARS technique is often used in the study of turbulent flames. Fast and accurate algorithms are needed for fitting dual-pump CARS spectra for temperature and multiple chemical species. This paper describes the development of such an algorithm. The algorithm employs sparse libraries, whose size grows much more slowly with number of species than a conventional library. The method was demonstrated by fitting synthetic "experimental" spectra containing 4 resonant species (N2, O2, H2 and CO2), both with noise and without it, and by fitting experimental spectra from a H2-air flame produced by a Hencken burner. In both studies, weighted least squares fitting of signal, as opposed to least squares fitting signal or square-root signal, was shown to produce the least random error and minimize bias error in the fitted parameters.
Cutler, Andrew D.↗
An Archival Chandra and XMM-Newton Survey of Type 2 Quasars
In order to investigate obscuration in high-luminosity type 2 active galactic nuclei (AGNs), we analyzed Chandra and XMM-Newton archival observations for 71 type 2 quasars detected at 0.05 < z < 0.73, which were selected based on their [O III] lambda5007 emission lines. For 54 objects with good spectral fits, the observed hard X-ray luminosity ranges from 2 × 10(exp 41) to 5.3 × 10(exp 44) erg s(exp −1), with a median of 1.1 × 10(exp 43) erg s(exp −1). We find that the means of the column density and photon index of our sample are log N(sub H) = 22.9 cm(exp −2) and gamma = 1.87, respectively. From simulations using a more physically realistic model, we find that the absorbing column density estimates based on simple power-law models significantly underestimate the actual absorption in approximately half of the sources. Eleven sources show a prominent Fe K alpha emission line (EW>100 eV in the rest frame) and we detect this line in the other sources through a joint fit (spectral stacking). The correlation between the Fe K alpha and [O III] fluxes and the inverse correlation of the equivalent width of the Fe Kα line with the ratio of hard X-ray and [O III] fluxes is consistent with previous results for lower luminosity Seyfert 2 galaxies. We conclude that obscuration is the cause of the weak hard X-ray emission rather than intrinsically low X-ray luminosities. We find that about half of the population of optically selected type 2 quasars are likely to be Compton thick. We also find no evidence that the amount of X-ray obscuration depends on the AGN luminosity (over a range of more than three orders of magnitude in luminosity).
quasars: general↗
Spectral Analysis in High Radiation Space Backgrounds with Robust Fitting
Spectral analysis software is tested for its ability to fit spectra from space. The approach, which emphasizes the background shape function, is uniquely suited to the identification of weak-strength nuclides in high-radiation background environments.
Lasche, G. P.↗
Correlative Analysis of GRBs Detected by Swift and Suzaku- WAM
It is now well known that a complete understanding of the energetics of the prompt phase of gamma-ray bursts (GRBs) requires full knowledge of the spectrum, extending at least as high as the peak energy (Epeak) of the vF(v) spectrum. Since most gamma-ray bursts (GRBs) have Epeak above the energy range (15-150 keV) of the Burst Alert Telescope (BAT) on Swift, a full understanding of the prompt emission from Swift GRBs requires spectral fits over as broad an energy range as possible. This can be completed for bursts which are simultaneously detected by Swift BAT and the Suzaku Wide-band All-Sky Monitor (WAM), which covers the energy range from 50-5000 keV. Between the launch of Suzaku in July 2005 and the end of 2008, there were 44 gamma-ray bursts (GRBs) which triggered both Swift and WAM and an additional 41 bursts which triggered Swift and were detected by WAM, but did not trigger. A joint BAT-WAM team has cross-calibrated the two instruments using GRBs, and we are now able to perform joint fits on these bursts to determine spectral parameters including Epeak. The results of broad spectral fits allows us to understand the distribution of Epeak for Swift bursts and to calibrate Epeak estimators when Epeak is within the BAT energy range. For those bursts with spectroscopic redshifts, we can calculate the isotropic energy and study various correlations between Epeak and other global burst parameters. Here we present the results of joint Swift/BAT-Suzaku/WAM spectral fits for 77 of the bursts jointly detected by the two instruments. We show that the distribution of spectral fit parameters is consistent with distributions from earlier missions and confirm that Swift bursts are consistent with earlier reported relationships between Epeak and isotropic energy. We show through time-resolved spectroscopy that individual burst pulses are also consistent with this relationship.
Krimm, H.A.↗
Version 2 Ozone Monitoring Instrument SO2 Product (OMSO2 V2): New Anthropogenic SO2 Vertical Column Density Dataset
The Ozone Monitoring Instrument (OMI) has been providing global observations of SO2 pollution since 2004. Here we introduce the new anthropogenic SO2 vertical column density (VCD) dataset in the version 2 OMI SO2 product (OMSO2 V2). As with the previous version (OMSO2 V1.3), the new dataset is generated with an algorithm based on principal component analysis of OMI radiances, but features several updates. The most important among those is the use of expanded lookup tables and model a priori profiles to estimate SO2 Jacobians for individual OMI pixels, in order to better characterize pixel-to-pixel variations in SO2 sensitivity, including over snow and ice. Additionally, new data screening and spectral fitting schemes have been implemented to improve the quality of the spectral fit. As compared with the planetary boundary layer SO2 dataset in OMSO2 V1.3, the new dataset has substantially better data quality, especially over areas that are relatively clean or affected by the south Atlantic anomaly. The updated retrievals over snow/ice yield more realistic seasonal changes in SO2 at high latitudes and offer enhanced sensitivity to sources during wintertime. An error analysis has been conducted to assess uncertainties in SO2 VCDs from both the spectral fit and Jacobian calculations. The uncertainties from spectral fitting are reflected in SO2 slant column densities (SCDs) and largely depend on the signal-to-noise ratio of the measured radiances, as implied by the generally smaller SCD uncertainties over clouds or for smaller solar zenith angles. The SCD uncertainties for individual pixels are estimated to be~0.15-0.3 DU (Dobson Units) between ~40°S and ~40°N and to be~0.2-0.5 DU at higher latitudes. The uncertainties from the Jacobians are approximately ~50-100% over polluted areas, and primarily attributed to errors in SO2 a priori profiles and cloud pressures, as well as the lack of explicit treatment for aerosols. Finally, the daily mean and median SCDs over the presumably SO2-free equatorial East Pacific have increased by only~0.0035 DU and ~0.003 DU respectively over the entire 15-year OMI record; while the standard deviation of SCDs has grown by only~0.02 DU or ~10%. Such remarkable long-term stability makes the new dataset particularly suitable for detecting regional changes in SO2 pollution.
OMI, SO2, Remote Sensing↗
A search for iron features in the EXOSAT spectral survey sources
The results of the search for iron features performed on 40 of the 48 EXOSAT spectral survey sources are discussed. In thrirteen of the sources the spectral fit was improved by the inclusion of an iron emission line at approximately 6.5 keV. In three of these, an ionized absorption edge improved the spectral fit. In one source, the spectral fit was improved by the addition of a cold iron absorption edge at approximately 7.1 keV. Line equivalent widths are not well determined. A statistical analysis of the line fits suggest an inverse correlation between line equivalent width and intrinsic source luminosity. Comparison of the host galaxy axial ratio with the line equivalent width suggests that larger equivalent widths occur in sources with more nearly face on host galaxies.
Leighly, K. M.↗
Testing the E(sub peak)-E(sub iso) Relation for GRBs Detected by Swift and Suzaku-WAM
One of the most prominent, yet controversial associations derived from the ensemble of prompt-phase observations of gamma-ray bursts (GRBs) is the apparent correlation in the source frame between the peak energy (E(sub peak)) of the nuF(nu) spectrum and the isotropic radiated energy, E(sub iso). Since most gamma-ray bursts (GRBs) have E(sub peak) above the energy range (15-150 keV) of the Burst Alert Telescope (BAT) on Swift, determining accurate E(sub peak) values for large numbers of Swift bursts has been difficult. However, by combining data from Swift/BAT and the Suzaku Wide-band All-Sky Monitor (WAM), which covers the energy range from 50-5000 keV, for bursts which are simultaneously detected ; one can accurately fit E(sub peak) and E(sub iso) and test the relationship between them for the Swift sample. Between the launch of Suzaku in July 2005 and the end of March 2009, there were 45 gamma-ray bursts (GRBs) which triggered both Swift/BAT and WAM and an additional 47 bursts which triggered Swift and were detected by WAM, but did not trigger. A BAT-WAM team has cross-calibrated the two instruments using GRBs, and we are now able to perform joint fits on these bursts to determine spectral parameters. For those bursts with spectroscopic redshifts.. we can also calculate the isotropic energy. Here we present the results of joint Swift/BAT-Suzaku/WAM spectral fits for 86 of the bursts detected by the two instruments. We show that the distribution of spectral fit parameters is consistent with distributions from earlier missions and confirm that Swift, bursts are consistent with earlier reported relationships between Epeak and isotropic energy. We show through time-resolved spectroscopy that individual burst pulses are also consistent with this relationship.
Krimm, H. A.↗
Polarization performance simulation for the GeoXO atmospheric composition instrument: NO2 retrieval impacts
NOAA’s Geostationary Extended Observations (GeoXO) constellation will continue and expand on the capabilities of the current generation of geostationary satellite systems to support US weather, ocean, atmosphere, and climate operations. It is planned to consist of a dedicated atmospheric composition instrument (ACX) to support air quality forecasting and monitoring by providing similar capabilities to missions such as TEMPO (Tropospheric Emission: Monitoring Pollution), currently planned to launch in 2023, and Ozone Monitoring Instrument (OMI), TROPOMI (TROPOspheric Monitoring Instrument), and GEMS (Geostationary Environment Monitoring Spectrometer) currently in operation. As the early phases of ACX development are progressing, design trade-offs are being considered to understand the relationship between instrument design choices and trace gas retrieval impacts. Some of these choices will affect the instrument polarization sensitivity (PS), which can have radiometric impacts on environmental satellite observations. We conducted a study to investigate how such radiometric impacts can affect NO2 retrievals by exploring their sensitivities to time of day, location, and scene type with an ACX instrument model that incorporates PS. The study addresses the basic steps of operational NO2 retrievals: the spectral fitting step and the conversion of slant column to vertical column via the air mass factor (AMF). The spectral fitting step was performed by generating at-sensor radiance from a clear sky scene with a known NO2 amount, the application of an instrument model including both instrument PS and noise, and a physical retrieval. The spectral fitting step was found to mitigate the impacts of instrument PS. The AMF-related step was considered for clear sky and partially cloudy scenes, where instrument PS can lead to errors in interpreting the cloud content, propagating to AMF errors and finally to NO2 retrieval errors. For this step, the NO2 retrieval impacts were small but non-negligible for high NO2 amounts; we estimated that a typical high NO2 amount can cause maximum of 0.25 x 10^15 molecules/cm^2 for a PS of 5%. These simulation capabilities were designed to aid in the development of a GeoXO atmospheric composition instrument that will improve our ability to monitor and understand the Earth’s atmosphere.
trace gases↗