Search NASA⌕ Search

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.

At least 37 records · Page 2

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 capabilities similar to missions such as TEMPO (Tropospheric Emission: Monitoring Pollution), currently planned to launch in 2023, as well as OMI (Ozone Monitoring Instrument), 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, for which 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 a maximum retrieval error of 0.25×1015 molec. 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.

Aaron Pearlman↗

Constraining the High-Energy Emission from Gamma-Ray Bursts with Fermi

We examine 288 GRBs detected by the Fermi Gamma-ray Space Telescope's Gamma-ray Burst Monitor (GBM) that fell within the field-of-view of Fermi's Large Area Telescope (LAT) during the first 2.5 years of observations, which showed no evidence for emission above 100 MeV. We report the photon flux upper limits in the 0.1-10 GeV range during the prompt emission phase as well as for fixed 30 s and 100 s integrations starting from the trigger time for each burst. We compare these limits with the fluxes that would be expected from extrapolations of spectral fits presented in the first GBM spectral catalog and infer that roughly half of the GBM-detected bursts either require spectral breaks between the GBM and LAT energy bands or have intrinsically steeper spectra above the peak of the nuF(sub v) spectra (E(sub pk)). In order to distinguish between these two scenarios, we perform joint GBM and LAT spectral fits to the 30 brightest GBM-detected bursts and find that a majority of these bursts are indeed softer above E(sub pk) than would be inferred from fitting the GBM data alone. Approximately 20% of this spectroscopic subsample show statistically significant evidence for a cut-off in their high-energy spectra, which if assumed to be due to gamma gamma attenuation, places limits on the maximum Lorentz factor associated with the relativistic outflow producing this emission. All of these latter bursts have maximum Lorentz factor estimates that are well below the minimum Lorentz factors calculated for LAT-detected GRBs, revealing a wide distribution in the bulk Lorentz factor of GRB outflows and indicating that LAT-detected bursts may represent the high end of this distribution.

Gehrels, Neil↗

Fe XXV temperatures in flares from the Yohkoh Bragg crystal spectrometer

Studies by Doschek et al. using P78-1 and Solar Maximum Misson (SMM) data have shown that the ratio of intensities of the Fe XXV and Ca XIX resonance lines can be expressed as a function of Fe XXV temperature. Using a more recent data set consisting of 13 flares observed by the Bragg crystal spectrometer (BCS) experiment on board Yohkoh, we find a nearly identical functional relationship between the same resonance line ratios and Fe XXV temperatures. We use this functional relationship to obtain resonance line ratio temperatures (T(sub RLR)) for each flare in our data set, and compare them with temperatures resulting from application of a simple spectral fitting method. (T(sub SSF)) to individal Fe XXV spectra. We also use a more involved free-parameter spectral fitting method to deduce temperatures (T(sub FSF)) from some of these spectra. On average, agreement between T(sub RLR) and T(sub SSF) improves as a flare progresses in time, with average agreements of 10.0% +/- 5.2%, 6.4% +/- 5.4%, and 5.0% +/- 3.9% over the rise, peak, and decay phases, respectively. Deviations between T(sub RLR) and T(sub FSF) are about the same or smaller. Thus, for most analysis purposes, all three methods yield virtually identical temperatures in flares. The somewhat poorer agreement between T(sub SSF) and T(sub RLR) during the earlier phases may be partially a result of difficulties in obtaining precise values for temperatures from spectral fits when blueshifts and large nonthermal broadenings are present in the spectra. Because of the high sensitivity of the Yohkoh BCS compared to that of BCS experiments on earlier spacecraft, we can for the first time consistently observe the heating phase of flares in Fe XXV.

Sterling, Alphonse C.↗

DETERMINATION OF LOW-ENERGY CUTOFFS AND TOTAL ENERGY OF NONTHERMAL ELECTRONS IN A SOLAR FLARE ON 2002 APRIL 15

The determination of the low-energy cutoff to the spectrum of accelerated electrons is decisive for the estimation of the total nonthermal energy in solar flares. Because thermal bremsstrahlung dominates the low-energy part of flare X-ray spectra, this cutoff energy is difficult to determine with spectral fitting alone. We have used anew method that combines spatial, spectral, and temporal analysis to determine the cutoff energy for the M1.2 flare observed with RHESSI on 2002 April 15. A low-energy cutoff of 24 +/- 2 keV is required to ensure that the assumed thermal emissions always dominate over nonthermal emissions at low energies (<20 keV) and that the spectral fitting results are consistent with the RHESSI light curves and images. With this cutoff energy, we obtain a total nonthermal energy in electrons of (1.6 +/- 1) x 10(exp 30) ergs that is comparable to the peak energy in the thermal plasma, estimated from RHESSI observations to be (6 +/- 0.6) x 10(exp 29) ergs assuming a filling factor of 1.

Sui, Linhui↗

Gamma-Ray Burst Afterglows as Probes of Environment and Blastwave Physics II: The Distribution of p and Structure of the Circumburst Medium

We constrain blastwave parameters and the circumburst media of a subsample of BeppoSAX Gamma-Ray Bursts. For this sample we derive the values of the injected electron energy distribution index, p, and the density structure index of the circumburst medium, k, from simultaneous spectral fits to their X-ray, optical and nIR afterglow data. The spectral fits have been done in count space and include the effects of metallicity, and are compared with the previously reported optical and X-ray temporal behaviour. Assuming the fireball model, we can find a mean value of p for the sample as a whole of 2.035. A statistical analysis Of the distribution demonstrates that the p values in this sample are inconsistent with a single universal value for p at the 3sigma level or greater. This approach provides us with a measured distribution of circumburst density structures rather than considering only the cases of k = 0 (homogeneous) and k = 2 (wind-like). We find five GRBs for which k can be well constrained, and in four of these cases the circumburst medium is clearly wind-like. The fifth source has a value of 0 less than or equal to k less than or equal to 1, consistent with a homogeneous circumburst medium.

Starling, R. L. C.↗

Gamma-Ray Burst Afterglows as Probes of Environment and Blast Wave Physics. II. The Distribution of rho and Structure of the Circumburst Medium

We constrain blast wave parameters and the circumburst media ofa subsample of 10 BeppoSAX gamma-ray bursts (GRBs). For this sample we derive the values of the injected electron energy distribution index, p, and the density structure index of the circumburst medium, k, from simultaneous spectral fits to their X-ray, optical, and NIR afterglow data. The spectral fits have been done in count space and include the effects ofmetallicity, and are compared with the previously reported optical and X-ray temporal behavior. Using the blast wave model and some assumptions which include on-axis viewing and standard jet structure, constant blast wave energy, and no evolution of the microphysical parameters, we find a mean value ofp for the sample as a whole of 9.... oa -0.003.0" 2 a_ statistical analysis of the distribution demonstrates that the p-values in this sample are inconsistent with a single universal value forp at the 3 _ level or greater, which has significant implications for particle acceleration models. This approach provides us with a measured distribution ofcircumburst density structures rather than considering only the cases of k ----0 (homogeneous) and k - 2 (windlike). We find five GRBs for which k can be well constrained, and in four of these cases the circumburst medium is clearly windlike. The fifth source has a value of 0 < k < 1, consistent with a homogeneous circumburst medium.

Starling, R. L. C.↗

New-Generation NASA Aura Ozone Monitoring Instrument (OMI) Volcanic SO2 Dataset: Algorithm Description, Initial Results, and Continuation with the Suomi-NPP Ozone Mapping and Profiler Suite (OMPS)

Since the fall of 2004, the Ozone Monitoring Instrument (OMI) has been providing global monitoring of volcanic SO2 emissions, helping to understand their climate impacts and to mitigate aviation hazards. Here we introduce a new-generation OMI volcanic SO2 dataset based on a principal component analysis (PCA) retrieval technique. To reduce retrieval noise and artifacts as seen in the current operational linear fit (LF) algorithm, the new algorithm, OMSO2VOLCANO, uses characteristic features extracted directly from OMI radiances in the spectral fitting, thereby helping to minimize interferences from various geophysical processes (e.g., O3 absorption) and measurement details (e.g., wavelength shift). To solve the problem of low bias for large SO2 total columns in the LF product, the OMSO2VOLCANO algorithm employs a table lookup approach to estimate SO2 Jacobians (i.e., the instrument sensitivity to a perturbation in the SO2 column amount) and iteratively adjusts the spectral fitting window to exclude shorter wavelengths where the SO2 absorption signals are saturated. To first order, the effects of clouds and aerosols are accounted for using a simple Lambertian equivalent reflectivity approach. As with the LF algorithm, OMSO2VOLCANO provides total column retrievals based on a set of predefined SO2 profiles from the lower troposphere to the lower stratosphere, including a new profile peaked at 13 km for plumes in the upper troposphere. Examples given in this study indicate that the new dataset shows significant improvement over the LF product, with at least 50% reduction in retrieval noise over the remote Pacific. For large eruptions such as Kasatochi in 2008 (approximately 1700 kt total SO2/ and Sierra Negra in 2005 (greater than 1100DU maximum SO2), OMSO2VOLCANO generally agrees well with other algorithms that also utilize the full spectral content of satellite measurements, while the LF algorithm tends to underestimate SO2. We also demonstrate that, despite the coarser spatial and spectral resolution of the Suomi National Polar-orbiting Partnership (Suomi-NPP) Ozone Mapping and Profiler Suite (OMPS) instrument, application of the new PCA algorithm to OMPS data produces highly consistent retrievals between OMI and OMPS. The new PCA algorithm is therefore capable of continuing the volcanic SO2 data record well into the future using current and future hyperspectral UV satellite instruments.

OMI↗

The structure of high-temperature flare plasma. II Comparison with observations

In Emslie, (1985) scaling law's which relate peak loop temperature to injected nonthermal electron flux in electron-heated models of the flare corona are derived. These predicted relationships are compared with recent coordinated observations in hard X-rays and soft X-rays. Satisfactory agreement is found for large events, while for smaller events the plasma electron temperature determined by soft X-ray spectral fitting is too high. A possible resolution of this apparent discrepancy through careful examination of the method used to determine the temperature of the soft X-ray emitting plasma is discussed. It is concluded that temperatures determined by spectral fitting over a series of lines are not necessarily representative of the true temperature of the plasma, since the fitting technique is plagued by the same difficulties as a more straightforward technique involving a single spectral line or portion of continuum (Craig and Brown, 1976). The differences between actual and derived temperatures are sufficiently large to remove the above discrepancy between observations and modeling of electron-heated coronae.

Emslie, A. G.↗

Disentangling AGN and Star Formation in Soft X-Rays

We have explored the interplay of star formation and active galactic nucleus (AGN) activity in soft X-rays (0.5-2 keV) in two samples of Seyfert 2 galaxies (Sy2s). Using a combination of low-resolution CCD spectra from Chandra and XMM-Newton, we modeled the soft emission of 34 Sy2s using power-law and thermal models. For the 11 sources with high signal-to-noise Chandra imaging of the diffuse host galaxy emission, we estimate the luminosity due to star formation by removing the AGN, fitting the residual emission. The AGN and star formation contributions to the soft X-ray luminosity (i.e., L(sub x,AGN) and L(sub x,SF)) for the remaining 24 Sy2s were estimated from the power-law and thermal luminosities derived from spectral fitting. These luminosities were scaled based on a template derived from XSINGS analysis of normal star-forming galaxies. To account for errors in the luminosities derived from spectral fitting and the spread in the scaling factor, we estimated L(sub x,AGN) and L(sub x,SF))from Monte Carlo simulations. These simulated luminosities agree with L(sub x,AGN) and L(sub x,SF) derived from Chandra imaging analysis within a 3sigma confidence level. Using the infrared [Ne ii]12.8 micron and [O iv]26 micron lines as a proxy of star formation and AGN activity, respectively, we independently disentangle the contributions of these two processes to the total soft X-ray emission. This decomposition generally agrees with L(sub x,SF) and L(sub x,AGN) at the 3 sigma level. In the absence of resolvable nuclear emission, our decomposition method provides a reasonable estimate of emission due to star formation in galaxies hosting type 2 AGNs.

Disentangling↗

Estimating the Properties of Hard X-Ray Solar Flares by Constraining Model Parameters

We wish to better constrain the properties of solar flares by exploring how parameterized models of solar flares interact with uncertainty estimation methods. We compare four different methods of calculating uncertainty estimates in fitting parameterized models to Ramaty High Energy Solar Spectroscopic Imager X-ray spectra, considering only statistical sources of error. Three of the four methods are based on estimating the scale-size of the minimum in a hypersurface formed by the weighted sum of the squares of the differences between the model fit and the data as a function of the fit parameters, and are implemented as commonly practiced. The fourth method is also based on the difference between the data and the model, but instead uses Bayesian data analysis and Markov chain Monte Carlo (MCMC) techniques to calculate an uncertainty estimate. Two flare spectra are modeled: one from the Geostationary Operational Environmental Satellite X1.3 class flare of 2005 January 19, and the other from the X4.8 flare of 2002 July 23.We find that the four methods give approximately the same uncertainty estimates for the 2005 January 19 spectral fit parameters, but lead to very different uncertainty estimates for the 2002 July 23 spectral fit. This is because each method implements different analyses of the hypersurface, yielding method-dependent results that can differ greatly depending on the shape of the hypersurface. The hypersurface arising from the 2005 January 19 analysis is consistent with a normal distribution; therefore, the assumptions behind the three non- Bayesian uncertainty estimation methods are satisfied and similar estimates are found. The 2002 July 23 analysis shows that the hypersurface is not consistent with a normal distribution, indicating that the assumptions behind the three non-Bayesian uncertainty estimation methods are not satisfied, leading to differing estimates of the uncertainty. We find that the shape of the hypersurface is crucial in understanding the output from each uncertainty estimation technique, and that a crucial factor determining the shape of hypersurface is the location of the low-energy cutoff relative to energies where the thermal emission dominates. The Bayesian/MCMC approach also allows us to provide detailed information on probable values of the low-energy cutoff, Ec, a crucial parameter in defining the energy content of the flare-accelerated electrons. We show that for the 2002 July 23 flare data, there is a 95% probability that Ec lies below approximately 40 keV, and a 68% probability that it lies in the range 7-36 keV. Further, the low-energy cutoff is more likely to be in the range 25-35 keV than in any other 10 keV wide energy range. The low-energy cutoff for the 2005 January 19 flare is more tightly constrained to 107 +/- 4 keV with 68% probability.

X-rays↗

X-ray Spectra and Multiwavelength Machine Learning Classification for Likely Counterparts toFermi3FGL Unassociated Sources

We conduct X-ray spectral fits on 184 likely counterparts to Fermi-LAT 3FGL unassociated sources. Characterization and classification of these sources allows for more complete population studies of the high-energy sky. Most of these X-ray spectra are well fit by an absorbed power law model, as expected for a population dominated by blazars and pulsars. A small subset of 7 X-ray sources ave spectra unlike the power law expected from a blazar or pulsar and may be linked to coincident stars or background emission. We develop a multiwavelength machine learning classifier to categorize unassociated sources into pulsars and blazars using gamma- and X-ray observations. Training a random forest procedure with known pulsars and blazars, we achieve a cross-validated classification accuracy of 98.6%. Applying the random forest routine to the unassociated sources returned 126 likely blazar candidates (defined as P(bzr) ≥ 90%) and 5 likely pulsar candidates (P(bzr) ≤ 10%). Our new X-ray spectral analysis does not drastically alter the random forest classifications of these sources compared to previous works, but it builds a more robust classification scheme and highlights the importance of X-ray spectral fitting. Our procedure can be further expanded with UV, visual, or radio spectral parameters or by measuring flux variability.

Stephen Kerby↗

Obtaining Electron Temperatures and Flow Speeds from Thomson Scattered Coronal Emission Observed during the 29 March 2006 Total Solar Eclipse in Libya

An experiment to measure the electron temperature and flow speed in the solar corona by observing the visible K-coronal spectrum was conducted during the total solar eclipse on 29 March 2006 in Libya. New corona1 models accounting for the effect of electron temperature and flow on the resulting K-corona spectrum were used to interpret the observations. Results show electron temperatures of 1.10 +/- 0.05, 0.98 +/- 0.12, and 0.70 +/- 0.08 MK, at l.l{\it R)$-{\odot)$ in the solar north, east and west, respectively, and 0.93 +/- 0.12 MK, at 1.2 R(sub sun) in the solar east. The corresponding outflow speeds obtained from the spectral fit are 103 +/- 92, 0 + 10, 0 + 10, and 0 + 10 km/s. Since the observations are taken only at 1.1 and 1.2 R(sub sun) these velocities , consistent with zero outflow, are in agreement with expectations and provide additional confirmation that the spectral fitting method is working.

Davila, Joseph M.↗

Electron Temperatures and Flow Speeds of the Low Solar Corona: MACS Results from the Total Solar Eclipse of 29 March 2006 in Libya

An experiment was conducted in conjunction with the total solar eclipse on 29 March 2006 in Libya to measure both the electron temperature and its flow speed simultaneously at multiple locations in the low solar corona by measuring the visible K-coronal spectrum. Coronal model spectra incorporating the effects of electron temperature and its flow speed were matched with the measured K-coronal spectra to interpret the observations. Results show electron temperatures of (1.10 +/- 0.05) MK, (0.70 +/- 0.08) MK, and (0.98 +/- 0.12) MK, at 1.1 Solar Radius from Sun center in the solar north, east and west, respectively, and (0.93 +/- 0.12) MK, at 1.2 Solar Radius from Sun center in the solar west. The corresponding outflow speeds obtained from the spectral fit are (103 +/- 92) km/s, (0 + 10) km/s, (0+10) km/s, and (0+10) km/s. Since the observations were taken only at 1.1 Solar Radius and 1.2 Solar Radius from Sun center, these speeds, consistent with zero outflow, are in agreement with expectations and provide additional confirmation that the spectral fitting method is working. The electron temperature at 1.1 Solar Radius from Sun center is larger at the north (polar region) than the east and west (equatorial region).

Reginald, Nelson L.↗

Spectral variations in the low-mass X-ray binary 4U 0614+09

We report the X-ray spectra of a binary 4U 0614+09 observed with EXOSAT on five different occasions in 1984 and 1986. In four observations in 1984, the source intensity was 'low' compared to its intensity in the 1986 observation. The source was steady on timescales of 10 to 1000 s in each of the 1984 observations but highly variable during the 1986 variations. The best-fit spectral models for the continuum are a blackbody spectrum plus an emission from a Comptonized region. In addition, a component due to reflection from a cold disk around the source was needed as an absorbiton edge due to cold or ionized iron was seen near 7 keV with good significance in most of the spectra, particularly in the 'low' state. The relation of the hardness ratios with intensity follow the behavior seen in the 'Atoll' sources. We have studied the behavior of the hardness ratios with the variations in the parameters of the best-fit spectral models. Most of the variability is associated with the blackbody component thought to arise from the boundary layer on the surface of the neutron star. During the 'high' state in 1986, the blackbody temperature is higher, its size smaller, and its contribution significantly greater than in the 'low' state. The other notable difference between the two states is the lower opacity in the Comptonized region, and a significantly higher absorption in the 'high' state. A weak line emission (equivalent width = 10 to 40 eV) is detected during the 'high' state. The line energy lies between 6 and 7 keV wh ich implies its origin from ionized iron.

Singh, K. P.↗

The detection of large amounts of cool, x ray absorbing gas in distant clusters of galaxies. What does this mean?

We present an x-ray spectral study of 12 distant (z = 0.17-0.54) rich clusters of galaxies observed with the Einstein Observatory Imaging Proportional Counter. These x-ray spectral data show evidence for substantial excess absorptions beyond those expected in the galaxy, indicating the presence of large amounts of x-ray absorbing cool gas in these distant clusters. The mean value of the excess absorptions corresponds to an absorbing gas column density approximately greater than 10(exp 21)/sq cm. We calculate the x-ray luminosities of the clusters with observed fluxes only in the 0.8-3.5 keV band where the fluxes are less effected by the absorptions, and use the temperature-to-luminosity correlation (known only for nearby clusters) to estimate the temperatures of the hot intracluster medium (ICM) in the distant clusters. These temperature estimates, together with the spectral fits, provide further constraints on the column densities in the individual clusters. For the cluster CL 0016+16, the lower limit on the column density is found to be 8 x 10(exp 20)/sq cm at the 99 percent confidence limit. We also show that the ratio of the temperature obtained from the spectral fit to the temperature expected from the correlation tends to decrease with increasing look-back time, indicating possible temperature evolution of the hot ICM in the recent past. The inclusion of this evolutionary effect further increases the absorptions required in fitting the spectra.

Wang, Qingde↗

Band analysis by spectral curve fitting

A method of estimating the values of the parameters in the models describing the positions, widths, and intensities of the lines in rotation-vibration bands of gases, without the need for line by line analysis, is described. To illustrate the technique, portions of the 1-0 bands of HCl and CO have been analyzed. The values of up to 27 parameters, their standard deviations, and the correlations between the parameters required to describe the spectra have been obtained.

Lin, C. L.↗

Determination of total column water vapor in the atmosphere at high spatial resolution from AVIRIS data using spectral curve fitting and band ratioing techniques

Techniques are developed for quantitative retrievals of high spatial resolution column atmospheric water vapor that is largely contained in the lower portion of the troposphere. One method consists of curve fitting observed spectra with simulated spectra in the 1.14 microns or the 0.94 micron water vapor band absorption region. The other method is a simple band ratioing technique, which requires less computer time than the curve fitting method. The advantage of the technique over humidity sounding by IR emission measurements is that the retrieved column water vapor amounts over land surfaces have significantly higher precision.

Gao, Bo-Cai↗

Deconstructing the Spectrum of the Soft X-ray Background

The soft X-ray background in the 0.1-1.0 keV band is known to be produced by at least three sources; the Local Hot Bubble (LHB), the extragalactic power law (EPL), and a seemingly galactic component that lies outside the bulk of the absorption that is due to the ISM of the galactic disk. This last component, which we call the Trans-Absorption Emission (TAE), has been modeled by a number of groups who have derived disparate measures of its temperature. The differences have arisen from differing assumptions about the structure of the emitting gas and unrecognized methodological difficulties. In particular, spectral fitting methods do not uniquely separate the TAE from the foreground emission that is due the LHB. This "degeneracy" can be resolved using the angular variation of the absorption of the TAE. We show that the TAE cannot be characterized by a single thermal component; no single-component model can be consistent with both the spectral energy distribution of the TAE emission and the angular variation due to absorption by the galactic disk. We use the angular anticorrelation of the ROSAT All-Sky Survey with the galactic absorption to separate local from distant emission components, and to fit the spectral energy distribution of the resulting distant emission. We find that the emission is best described by a two-thermal-component model with logT(sub S) = 6.06(sup +0.14, sub -0.12) and log T(sub H) = 6.42(sup +0.14, sub -0.12). This two-thermal-component TAE fits the ROSAT spectral energy distribution significantly better than single-component models, and is consistent with both angular variation and spectral constraints.

Kuntz, K. D.↗