The Gas-cell Simulator models gas spectra. With just a few clicks you can accurately simulate the spectra of most of the important trace gas species in the Earth's atmosphere, such as water vapor, carbon dioxide, ozone, methane and dozens more. Spectra are calculated in the user-specified waveband with full resolution, fine enough to accurately capture the shape of even the narrowest molecular absorption line. The underlying software is scientifically rigorous and has been used in support of a number of NASA research projects.
The basic setup models the transmittance through a "gas cell" -- an imaginary sealed container with transparent walls. The transmittance is simply the fraction of radiation that makes it through the gas cell. This is illustrated below. The transmittance, τ, is the ratio of the intensity of the transmitted radiation, I, to the incident,
I0. You control the pressure, P, temperature, T and length, L, of the gas cell, the concentration of the absorbing
gas inside, q, and the range of wavenumbers over which to compute the spectrum. The Gas-cell Simulator then accurately calculates the full-resolution transmittance spectrum.
Subscribers have access to some advanced features. They can simulate mixtures of several gases and can "stack" up to six gas cells end-to-end, controlling the conditions of each. They can choose to compute emission spectra rather than transmittance spectra, and add an optional blackbody source behind the gas cells. Instrument lineshape functions that "apodize" or smooth the spectra can also be included. These features are illustrated below.
Click here for a comprehensive list of features available to subscribers.
Choose your waveband, set up your gas cells, and press Calculate! Many variations are possible, so we've divided the inputs into functional "tabs" to make things more manageable.
Observer Tab - choose the waveband and any other instrument parameters
Gas Cells Tab - specify the gases to simulate and gas cell conditions
Source Tab (optional) - include a blackbody background
Plot Options Tab (optional) - control plot features like log or linear scaling
My Settings Tab (optional) - save and retrieve simulation setups
All inputs have default values that are sensible, and simple messages will inform you if an input value is out of range. Click here for detailed instructions.
The Gas-cell Simulator uses a line-by-line model called LinePakTM to accurately model the spectra of molecular absorption lines. The shape, position and depth of each line are described by parameters in a line list such as HITRAN, and adjusted for the user-specified pressure, temperature, and gas concentration. We maintain an up-to-date copy of the HITRAN line parameters, and allow alternate line lists to be used as well. Details of the calculation are given here.
The Line List browser lets you quickly and easily see where gases have absorption lines, and how strongly they absorb, without the time and effort of performing a detailed simulation. This tool displays the contents of HITRAN, GEISA and other line list databases. You can also download the line list data here, selecting the complete datasets or just the molecules and wavebands you want.
Why do I only have to specify the length of the gas cell, not its volume?
The length of the gas cell is the only dimension that matters because the simulation is of a single beam through the gas cell. Only molecules along the line-of-sight will affect this radiation.
At all but very low wavenumbers, the LinePakTM software uses the Voigt profile. The Voigt profile combines the pressure-induced line broadening effects of molecular collisions, and the thermal-induced Doppler effect resulting from their randomly oriented motions. At very low wavenumbers (below 200 cm-1) the Van Vleck Weisskopf line shape is used.
LinePakTM can easily accommodate more complex line shapes, such as the Galatry or speed-dependent profiles, but the
spectroscopic parameters needed for these are currently available only for a small number of absorption lines. Note that these provide only a very small improvement, and in most situations, the Voigt profile adds no significant error.
The spectral resolution is chosen to be fine enough to fully capture the line shape of the thinnest absorption line in the spectrum. The Gas-cell Simulator software first surveys all absorption lines needed for the calculation, including lines centered outside the spectral bandpass whose wings may contribute to the measured spectral range. The algorithm then determines the smallest width of all the absorbing lines, adjusted for pressure and temperature. The spectral resolution is set to 1/20th of this smallest halfwidth.
Yes. In many situations, the spectrum in a specified bandpass will have significant contributions from absorption lines centered outside the bandpass. To ensure that this is modeled correctly, the Gas-cell Simulator software calculates the spectra over a larger waveband, extended by 120 cm-1 on either side of the specified range.
The Gas-cell Simulator reports spectra as would be measured in vacuum. (Outside the gas cell a vacuum is assumed). In some applications, spectra are measured in air. Because the speed of light is slower in air than in vacuum,
the observed wavenumber (1/c in cm-1) is slightly larger than it would have been if it had been observed in vacuum. To compare Gas-cell Simulator results to spectra measured in air (or some other media), this shift must be accounted for.
Subscribers can import their own spectral line lists using the My Line
Lists tab on the Linelist Browser tool. The line list must be a text file in HITRAN format. Up to 100,000
lines are allowed. The list can include only molecules/isotopologues that are already in the public HITRAN distribution. (If you have line
parameters for molecules/isotopologues not in HITRAN, please contact us and we will help get your
data uploaded.) It can take several minutes for the conversion, validation and upload process. Once your data is uploaded, your line list will be visible under My Line Lists and available for browsing or spectral calculations. Other subscribers will not be able to access your line lists.
There are many benefits of becoming a subscriber. Some of the more popular features include the ability to simulate large wavebands, multiple gases and multiple gas cells. Subscribers have access to publication quality graphics and can download the raw data in text format. They can store their simulation setups for recall in later sessions. Subscribers can also upload their own custom line lists for testing and validation. Emission spectra can be computed with an optional blackbody source function, and instrument lineshape functions can be customized and applied to the spectra. Click here for a comparison of the features available with and without a subscription.
The basic functionality of all the tools at SpectralCalc.com is freely available without subscribing. You can simulate transmittance spectra of gases with the same fast, accurate software used by many NASA remote sensing missions, and visually browse the HITRAN and other line list databases. The Blackbody Calculator, Atmosphere Browser and Solar Calculator are all freely available. Click here for a description of the many additional features available with a subscription.
Spectra with up to 1 million points can be calculated. The size of the waveband that this supports depends on the spectral resolution of the spectrum. In general, low pressure spectra have very narrow spectral features, and require very fine resolution, especially at low wavenumbers. (Note: the spectral resolution of an atmospheric path will be governed by its highest altitude.) Here are some typical examples:
Maybe nothing. This simply means there are no absorption lines for the selected gas in the specified waveband. You may wish to visit the Line List Browser to identify spectral regions in which each gas does have absorption lines, and adjust your waveband and/or gas accordingly.
Yes. On the Gas Cell calculator, select the isotopologue to be simulated and specify the desired abundance. You can simulate mixtures of different isotopologues and control the concentration of each. Note: If “All” isotopologues are selected for a molecule, they are automatically scaled by their natural abundance so that they sum to the specified VMR.
Atmospheric path calculations automatically account for the natural abundance of each isotopologue. You can add an individual isotopologue as one of your six species to further control its concentration. Note: all isotopologues of a molecule have the same vertical structure, with their individual concentrations scaled everywhere by their natural abundance. The vertical profiles of all species are available in the Atmosphere Browser.
Line intensity is a parameter that indicates the strength of a spectral line. These are used in our simulations to compute absorption and emission spectra. The Line List Browser shows the line intensities of each molecule in a spectral range. Viewing the line intensities gives a quick indication of the spectrum of the chosen molecule(s), without having to perform any lengthy calculations. This is very useful for determining which molecules absorb in a particular waveband, or for identifying spectral lines in measured spectra. A detailed description of how these intensities are used in the simulation of a spectrum is available here.
Because molecular species have widely differing concentrations in the atmosphere, we offer the option of plotting intensities scaled by the atmospheric concentration (in molecules per cm3). Simply choose an atmosphere model and altitude, and the line intensities will be scaled by the species’ concentrations at that altitude. The raw intensities, as supplied by HITRAN and other databases, have units of cm-1/(molecule/cm-2) at 296 K, while the scaled line intensities have units of cm-1/cm (wavenumber per unit path length).
When plotting line intensities of individual isotopologues of a molecule, you have the option of plotting the intensities with or without scaling by their natural abundance.
You can select an instrument function (Observer Tab) to model spectra as would be obtained from an instrument with a specified spectral resolution. Adjusting the width of the instrument function will control the resolution of the resulting spectrum. Please see the Instrument Function FAQ for details.
SpectralCalc models spectra at very fine resolution – fine enough to capture every significant spectral feature. We sometimes call these “monochromatic” spectra. Real instruments, however, always have some inherent resolution and produce somewhat smoothed spectra, effectively convolving the received monochromatic spectra with the instrument’s spectral response function. To simulate a spectrum with an instrument function, select the desired function type and width from the Observer Tab. The available types are shown here. Click on the image to see a larger version
There are several important things to note when applying instrument functions:
The specified width must be less than 250 cm-1 and less than 1/10th of the waveband. The width can be arbitrarily small, but widths smaller than the inherent resolution of the underlying spectrum are ignored.
Larger widths produce coarser spectra (40 points per instrument function width), but the underlying calculation is always done at high spectral resolution.
Instrument functions are computed and applied in wavenumbers. A width specified in microns is converted to the equivalent width in wavenumbers at the center of the bandpass. As a result, an instrument function with its width specified in microns will actually have different effective widths in different spectral bandpasses.
Instrument functions should only be applied to the final spectrum. For example, to model the combined spectrum of sequential absorbing paths, we should compute and save the monochromatic transmittance spectra of each. Then on the My Spectra Tab, we can compute their product, and finally apply an instrument function to the resulting monochromatic spectrum. Applying the instrument function to the individual spectra and then combining them will produce erroneous results.
To register, click on Subscribe, fill out the "New Users" information and press the "Register" button. Once registered, you can choose a subscription plan and immediately gain full access to all tools and features. Simply press the appropriate "Buy Now" button and proceed with check-out at the PayPal website. Note that you do not need a PayPal account, only a valid credit card. Your subscription will be active immediately after completion of the payment.
Yes. The underlying algorithm was first described in the following paper:
Gordley et al., LINEPAK: Algorithm for Modeling Spectral Transmittance and Radiance, J. Quant. Spectrosc. Radiat. Transfer Vol. 52, No. 5, pp.563-580, 1994
Administrators can manage the group account from "My Account". Here you will see a list of all users currently in the group. You can add or delete users, and can designate an additional administrator if desired. Users can also join the group by logging in and entering the group code in "My Account".
Notes:
Users must be registered.
Group accounts can have up to 12 users, including the administrator(s).
Administrator(s) will receive email confirmation when a user joins the group.
A group account can have a maximum of 2 administrators.
Group codes are chosen by the administrator at time of purchase.
To join a group account, go to "My Account" and click on "Join Group". Here you will be prompted to enter the group code. Contact the group administrator for this code.
Notes:
Users must be registered.
Group accounts can have up to 12 users, including the administrator(s).
Administrator(s) will receive email confirmation when a user joins the group.
Group codes are chosen by the administrator at time of purchase.
Subscribers can import their transmittance or radiance spectra by using the Upload tab on the
My Spectra tool. Enter the full name (including directory path) of the file to upload, or use
the Browse button to locate the file on your computer.
Choose a short descriptive name to identify the spectrum, and select Transmittance or
Radiance for the spectrum type. A description of the file can be typed into the text box
provided. Once uploaded, the spectrum will be available on the My Spectra tab, where it can be
viewed, downloaded, deleted, etc. (Other users do not have access to your stored spectra.)
File requirements: The uploaded spectrum file must be in plain text format with two columns,
the first column containing the wavenumber (cm-1) values, and the second containing the
transmittance or radiance values (W/m2/sr/cm-1). Files must be no bigger than 20 MB, and
contain no more than 400,000 spectral points. Header lines (lines that begin with something other than a space, a digit or a decimal point) are
allowed at the beginning of the file, but are ignored. The columns can be space-, tab- or comma-delimited. Any text following the second column is ignored.
Wavenumbers are assumed to be equally spaced and increasing. Wavenumber spacing must be in the range 1e-8 to 100 cm-1.
(Only the first and last wavenumber are actually stored. The grid is reconstructed each time the data are retrieved from our database.)
Transmittance values must be between -1 and 2, and radiance values must be between 0 and 1000 W/m2/sr/cm-1.
The atmosphere is modeled as graduated concentric spherical shells. The number of shells depends on the path length and altitude range. For example, a path from the ground to 120 km (the top of our model atmospheres) is split into 19 shells: 250 meters thick at the surface, growing to 10 km thick at high altitudes. The temperature and gas concentrations in each shell are taken from the selected model atmosphere. These model data can be viewed/downloaded with the Atmosphere Browser tool. Within each shell, a mass-weighted computation is performed, as described by Gordley and Marshall. Refraction effects are rigorously included in all paths. (Note that radiance spectra for up-looking and downlooking paths are different, whereas transmittance spectra are the same for both directions.)
Yes, continuum absorption/emission for N2, O2 and H2O are included.
These effects are automatically included as needed for simulations involving these molecules. Other effects, such as Rayleigh (molecular) and Mie (particulate) scattering, are not yet included.
Continuum absorption/emission effects for N2, O2 and H2O are modeled using the approaches described by Clough, et al. [1] and Lafferty et al. [2]. These continuum
models augment the standard "line-by-line" calculation, more accurately accounting for far-wing effects at high pressures.
For O2 and N2, the spectral ranges affected are 1345-1820 cm-1 and 2105-2620 cm-1 respectively. The H2O continuum is significant at wavenumbers
below 600 cm-1 and between 1200 and 2000 cm-1.
Click on an image to see a larger version
Gas cell simulations of N2, O2 and H2O spectra in any of the above wavebands will automatically include the appropriate continuum model(s).
No action is needed from the user.
For all three gases, the effect is proportional to pressure, and is negligible at pressures less than ~10mbar. Temperature dependence of the continua
is included in the model, derived from laboratory measurements between 70 and 350 K. Continuum absorption for spectra at higher temperatures is based on extrapolation
of these data, and therefore may contain some inaccuracies. Note, however that the continuum effects are inversely proportional to temperature, so the overall
magnitude become negligible for high temperatures.
[1] Clough, S. A., F. X. Kneizys, and R. W. Davies, "Lineshape and the water vapor continuum," Atmos. Res. 23, 229-241 (1989)
[2] Lafferty, W. J., A. M. Solodov, A. Weber, W. B. Olson, and J. Hartmann, "Infrared collision-induced absorption by N2 near 4.3 µm for atmospheric applications: measurements and empirical modeling," Appl. Opt. 35, 5911-5917 (1996)
Cross sections are normalized molecular absorption spectra. These are often used to simulate broad absorption features that are not well modeled using line-by-line techniques, or for which adequate line parameters are not available. Cross section spectra are typically given in units of cm2 per molecule at each wavenumber.
Infrared spectra of polyatomic "heavy" molecules such as CFCs, and UV absorption features of several important atmospheric species, most notably ozone, are better simulated using cross sections than line-by-line models. Cross section spectra for these molecules can be accessed in the Line List Browser, by selecting Plot Type-->Cross Sections.
We're working on integrating the capability of simulating spectra using cross sections to our online tools, but have not yet validated and tested a general-purpose interface for this.
The total number of points in a spectrum is the waveband size divided by the spectral resolution (point spacing in cm-1). There are several things you can do to reduce the total number of points:
Reduce the waveband. This will directly reduce the number of points. In general, the lower wavenumber regions require the most points.
Increase the pressure. Gases at low pressure (< 0.1 mbar) have very thin Doppler-broadened lines, and require very fine spacing in the calculation.
Limit the altitude of atmospheric paths. High altitudes have low pressures, and thus need fine resolution. Restricting your path to low altitudes should help.
Avoid low wavenumbers at low pressure. At low-pressure, absorption lines have Doppler widths proportional to wavenumber, so low wavenumber/low pressure spectral calculations are particularly difficult.
In some situations, the calculation may need to be split up into a number of spectral segments, and each calculated individually.
The spectral resolution (point spacing in cm-1) is chosen to be fine enough to fully capture the line shape
of the thinnest absorption line in the spectrum. Our software surveys all absorption lines needed and finds the
smallest width. The spectral resolution is set to 20% of the thinnest line's width. In general, at
pressures above 1 mbar, the lines are pressure broadened with widths proportional to pressure.
At low pressures (below 1 mbar), only Doppler broadening is important, and the widths are proportional to
wavenumber (position of line). A more detailed description is available
here.
The atmosphere models extend to 120 km. If the observer height or target height is above this altitude, the portion of the atmosphere above 120 km is represented by an isothermal layer with the same temperature and VMR conditions as 120 km. The pressure in this layer is hydrostatically decreasing and extends for 3 scale heights.
The tangent height of a limb path is the altitude of its closest approach to the Earth’s surface. Atmospheric refraction bends the path, and the refracted tangent height is lower than the “apparent” tangent height, as shown here:
The user supplies the apparent tangent height. SpectralCalc calculates the refracted tangent height. Refraction is more severe at low altitudes, and depends somewhat on the model atmosphere. Rays with very low apparent tangent heights may intersect the Earth. In this case, an error message is displayed, and you will need to increase the apparent tangent height until the refracted ray does not hit the Earth.
Of course not. However, SpectralCalc can be used to demonstrate some of the fundamental principles. Spectroscopy of greenhouse gases is critical in atmospheric modeling, and SpectralCalc accurately simulates high-resolution molecular spectra of atmospheric paths (currently for Earth or Mars). Still, while these simulations are useful and important, any sufficiently complete climate model has to include a host of other significant components, including surface albedo, clouds and aerosols, solar variability and ocean cycles. For more details on atmospheric climate modeling, a nice introduction is available here: http://geosci.uchicago.edu/~rtp1/papers/PhysTodayRT2011.pdf
Each isotopologue has a unique spectrum and must be simulated separately. Unless the user specifies otherwise, it is assumed that the gas sample includes all isotopologues of the specified molecule, in their natural abundances. The natural abundances are implicitly accounted for in our calculations by pre-scaling the spectral line intensities.
To simulate spectra of a single isotopologue, or a sample with an elevated concentration of a minor isotopologue, simply switch the isotopologue drop-down menu from All to the desired isotopologue. When a specific isotopologue is simulated, our software removes the pre-scaling from the HITRAN line intensities, so that the simulation reflects the concentration supplied by the user.
To make the simulation of atmospheric spectra more efficient, HITRAN and other databases have pre-scaled the spectral line intensities for each isotopologue by the natural abundance of the isotopologue. For example, the line intensities for the primary carbon dioxide isotopologue (16O12C16O) are scaled by a factor of 0.984204, whereas the line intensities of the secondary isotopologue (16O13C16O) are scaled by a much smaller factor, 0.0110574. This pre-scaling implicitly accounts for the relative concentration of the various isotopologues in atmospheric gases, making line-by-line calculations more efficient. In the event a single isotopologue is to be simulated, our software simply “unscales” the line-intensities from HITRAN (or other database) before proceeding with the calculation of the spectrum.
Note: any linelists uploaded by users are assumed to be pre-scaled in this same manner.
absorption--fractional amount of incident radiation absorbed by a gas. Absorption = 1 – transmittance. In general, absorption varies with the frequency of the radiation.
absorption line-- The decrease in a transmittance spectrum corresponding to the radiation absorbed by a particular transition
between quantum states of a molecule. Early spectra were recorded on strips of film, where molecular absorption features appeared as dark lines across the continuous spectrum, hence the name "absorption line".
abundance--The relative amount, or concentration, of a gas in a volume. This is typically specified by the gas' volume mixing ratio.
blackbody--a hypothetical object that absorbs all incident radiation while maintaining thermal equilibrium. No light is reflected
from or passes through a blackbody, but radiation is emitted, and is called blackbody radiation. The prefix "black" is used because at room temperature such an object would emit almost no visible light, appearing black to an observer. No physical object exactly fits this definition, but most behave at least in part as blackbodies. Planck's law gives the spectral distribution of radiation emitted by a blackbody.
cross sections--normalized molecular absorption spectra, given in units of cm2 per molecule at each wavenumber. These are used to model spectra where line parameters are insufficient or unavailable.
Doppler effect--the change in frequency of observed radiation resulting from motion of the source relative to the observer.
emissivity--a unitless spectral quantity equal to the ratio of energy emitted by an object to that emitted by a blackbody at the same temperature.
GEISA--a compilation of spectroscopic line parameters, similar to HITRAN (GEISA = Gestion et Etude des Informations Spectroscopiques Atmosphériques.)
HITEMP--the high-temperature analogue of HITRAN. The current compilation contains parameters for water vapor, carbon dioxide and carbon monoxide.
HITRAN--High-resolution transmission molecular absorption database. HITRAN is a compilation of spectroscopic parameters widely used for simulating the transmission and emission of radiation of gases in the atmosphere. See www.hitran.com for more details.
The HITRAN data can be downloaded to your computer (go to Line List Browser, Extract Data Tab). Subscribers can also upload their own line list files. The required format is the HITRAN format -- a fixed-field text file with the following fields
Parameter
Field width
FORTRAN format
molecule number
2
I2
isotopologue number
1
I1
transition frequency (cm-1)
12
F12.6
line Intensity
10
1PE10.3
Einstein A-coefficient
10
E10.3
air-broadened width
5
0PF5.4
self-broadened width
5
F5.4
lower-state energy
10
1PF10.4
temperature dependence (of air width)
4
0PF4.2
pressure shift
8
F8.6
upper vibrational quanta
15
A15
lower vibrational quanta
15
A15
upper local quanta
15
A15
lower local quanta
15
A15
error codes
6
6I1
reference codes
12
6I2
flag for line-mixing
1
A1
upper statistical weight
7
F7.1
lower statistical weight
7
F7.1
Isotopologue--Isotopologues of a molecule share the same molecular structure, but their individual atoms are different isotopes. For example, 16O12C16O and 16O13C16O are different isotopologues of CO2. Any sample of air will contain a mixture of CO2 isotopologues in fixed proportions. 98.4% of the CO2 molecules will be the primary isotopologue, 16O12C16O, while the next most abundant isotopologue, 16O13C16O, will account for about 1.1%. The remaining fraction is made up by even less common isotopologues. These fixed percentages are sometimes referred to as natural or terrestrial abundances, and are listed here for all the species in our database. Under certain circumstances, there are slight variations in these isotopic fractions, but for most atmospheric applications it is sufficient to assume the fractions are constant.
Note: The name isotope applies only to atoms, but is commonly (albeit incorrectly) used as a synonym for isotopologue.
Isotope--An isotope is a variant of an atom with a specific atomic number (i.e. number protons plus neutrons). For example, the most common isotope of hydrogen, 1H, contains one proton and no neutrons. A less common isotope is 2H, called deuterium, which contains one proton and one neutron.
Note: The name isotope applies only to atoms, but is commonly (albeit incorrectly) used as a synonym for isotopologue.
LinePakTM--The GATS library of radiative transfer software used by SpectralCalc. LinePak performs detailed "line-by-line" modeling of molecular absorption and emission. Efficient and flexible, LinePakTM is at the heart of analysis systems for many major atmospheric remote sensing missions including HALOE, LIMS,
SABER, CRISTA, CLAES and SOFIE. More details.
Line List--A list of parameters used in modeling molecular absorption spectra. Each molecular absorption line is described by a set of parameters including its center frequency, halfwidth, and intensity. These parameters allow the radiative transfer code (LINEPAK) to accurately model the absorption line shape at any temperature, pressure and gas concentration. See HITRAN for a more complete description of the parameters in a line list.
micron--10-6 meters. microns are a common unit for the wavelength of infrared radiation.
millibar (mbar)--a unit of pressure equal to 100 Pascals. Other common conversions are: 1 atmosphere = 1013.25 mbar, 1 lb/in2 = 68.9476 mbar
radiance--emitted power per unit surface area per solid angle per spectral interval. Common units are (W m-2 sr-1 cm-1) or (photons s-1 m-2 sr-1 m-1). Also referred to as emission.
Torr-- A unit of pressure equal to 1/760 atmosphere, or 101325/760 Pa.
transmittance--ratio of the intensity of radiation exiting a medium to the incident intensity. In general, transmittance varies with the frequency of radiation. Absorption = 1 – transmittance.
volume mixing ratio (vmr)--the fractional number of molecules of a species in a mixture. vmr = M/N where M is the number of molecules of the species in a volume and N is the total number of molecules. For trace gases, vmr is often given in parts per million or parts per billion.
wavenumber--the number of cycles of electromagnetic radiation per cm. The wavenumber, in cm-1, is equal to 10000 / (wavelength in microns).