ATMOS data processing and science analysis methods

Py4CAtS—PYthon for Computational ATmospheric Spectroscopy

Radiation is a key process in the atmosphere. Numerous radiative transfer codes have been developed spanning a large range of wavelengths, complexities, speeds, and accuracies. In the infrared and microwave, line-by-line codes are crucial esp. for modeling and analyzing high-resolution spectroscopic observations. Here we present Py4CAtS—PYthon scripts for Computational ATmospheric Spectroscopy, a Python re-implemen-tation of the Fortran Generic Atmospheric Radiation Line-by-line Code GARLIC, where computationally-intensive code sections use the Numeric/Scientific Python modules for highly optimized array processing. The individual steps of an infrared or microwave radiative transfer computation are implemented in separate scripts (and corresponding functions) to extract lines of relevant molecules in the spectral range of interest, to compute line-by-line cross sections for given pressure(s) and temperature(s), to combine cross sections to absorption coefficients and optical depths, and to integrate along the line-of-sight to transmission and radiance/intensity. Py4CAtS can be used in three ways: in the (Unix/Windows/Mac) console/terminal, inside the (I)Python interpreter, or Jupyter notebook. The basic design of the package, numerical and computational aspects relevant for optimization, and a sketch of the typical workflow are presented. In conclusion, Py4CAtS provides a versatile environment for “interactive” (and batch) line-by-line radiative transfer modeling.

Quantitative Fourier transform infrared diagnostics of the gas-phase composition using the HITRAN database and the equivalent width of the spectral features

This paper presents a strategy for quantification of medium resolution Fourier transform infrared (FT-IR) spectra. The approach is based on a comparison of the values of the equivalent width of spectral features determined from the measured FT-IR spectra with those calculated from corresponding molecular spectra simulated using spectroscopic parameters tabulated in the HITRAN database. Although the equivalent-width method is routinely applied in many high-resolution experiments, its potential is often ignored when the spectral resolution of the measurements is bigger than the width of the studied molecular transitions. Here with an eye on application of the method to the study of the biomass thermal decomposition products, we demonstrate the capability of the method for analysis of FT-IR spectra with moderate resolution. The method is validated for a number of molecules (NH3, CO, CH4, C2H2, C2H4, and NO) that are important products of biomass thermal decomposition. Namely, known amounts of gases were placed in a sample cell and their concentrations were determined using the method for different FT-IR settings. The agreement between the concentration values determined using the present method and those found in the sample cell was generally better than 10%. The paper also shows examples of application of the developed methodology for the analysis of FT-IR spectra from biomass pyrolysis.

Atmospheric Trace Molecule Spectroscopy (ATMOS) Experiment Version 3 data retrievals

Version 3 of the Atmospheric Trace Molecule Spectroscopy (ATMOS) experiment data set for some 30 trace and minor gas profiles is available. From the IR solar-absorption spectra measured during four Space Shuttle missions (in 1985, 1992, 1993, and 1994), profiles from more than 350 occultations were retrieved from the upper troposphere to the lower mesosphere. Previous results were unreliable for tropospheric retrievals, but with a new global-fitting algorithm profiles are reliably returned down to altitudes as low as 6.5 km (clouds permitting) and include notably improved retrievals of H2O, CO, and other species. Results for stratospheric water are more consistent across the ATMOS spectral filters and do not indicate a net consumption of H2 in the upper stratosphere. A new sulfuric-acid aerosol product is described. An overview of ATMOS Version 3 processing is presented with a discussion of estimated uncertainties. Differences between these Version 3 and previously reported Version 2 ATMOS results are discussed. Retrievals are available at http://atmos.jpl.nasa.gov/atmos.

Apodization effects in the retrieval of volume mixing ratio profiles

In remote sensing applications, spectra measured by Fourier-transform spectrometers are routinely apodized. A rigorous analysis approach would explicitly account for correlations induced in the covariance matrix by apodization, but these correlations are often ignored to simplify and speed up the processing. Using spectra measured by the Atmospheric Trace Molecule Spectroscopy missions, we investigated the effect of apodization on the retrieval of volume mixing ratio profiles for the case in which these correlations are ignored. Minor discrepancies occur between results for apodized and unapodized spectra, particularly when lines with a low signal-to-noise ratio are fitted. A set of microwindows is reported for O3 in the range of 1550-3350 cm(-1).

Vertical profiles of aerosol volume from high-spectral-resolution infrared transmission measurements. I. Methodology

The wavelength-dependent aerosol extinction in the 800-1250-cm(-1) region has been derived from ATMOS (atmospheric trace molecule spectroscopy) high-spectral-resolution IR transmission measurements. Using models of aerosol and cloud extinction, we have performed weighted nonlinear least-squares fitting to determine the aerosol-volume columns and vertical profiles of stratospheric sulfate aerosol and cirrus cloud volume. Modeled extinction by use of cold-temperature aerosol optical constants for a 70-80% sulfuric-acid-water solution shows good agreement with the measurements, and the derived aerosol volumes for a 1992 occultation are consistent with data from other experiments after the eruption of Mt. Pinatubo. The retrieved sulfuric acid aerosol-volume profiles are insensitive to the aerosol-size distribution and somewhat sensitive to the set of optical constants used. Data from the nonspherical cirrus extinction model agree well with a 1994 mid-latitude measurement indicating the presence of cirrus clouds at the tropopause.

Diode-Laser Measurements of He-, Ar-, and N2-Broadened HF Lineshapes in the First Overtone Band

High-resolution lineshapes of the HF P(3) and P(6) transitions (2-0 band) have been recorded in a 14.9-cm long absorption cell (T = 296 K) using a pair of (InGaAsP)-distributed feedback diode lasers operating near 1.31 and 1.34 µm, respectively. The absorption lineshapes were least-squares fit to Voigt, Rautian, and Galatry profiles to determine the collision widths and velocity-changing collision frequencies for the perturbers He, Ar, and N2. Motional narrowing was observed for the He-broadened lineshapes over a range of total pressures from 50 to 646 Torr. The Ar- and N2-broadened spectra showed a strong narrowing effect for pressures less than 400 Torr. At higher pressures, narrowing parameters could not be reliably determined due to the relatively large Ar- and N2-pressure broadening. The collision widths inferred from Voigt-profile fits were systematically less than those determined from Rautian and Galatry fits (up to 25% for the He-broadened P(6) transition). The measured pressure-induced broadening and shift coefficients were compared with other experimental data and full-quantal close-coupling (CC) scattering cross-section calculations. The measured Ar-broadening half-widths were in excellent agreement with calculated values, while the measured shifts were 10 and 20% less than calculated values for the P(3) and P(6) transitions, respectively. The measured N2-broadening and shift coefficients agree reasonably well with values determined from previous experiments. Copyright 1999 Academic Press.

Isotopic fractionation of stratospheric nitrous oxide

We propose an isotopic fractionation mechanism, based on photolytic destruction, to explain the 15N/14N and 18O/16O fractionation of stratospheric nitrous oxide (N2O) and reconcile laboratory experiments with atmospheric observations. The theory predicts that (i) the isotopomers 15N14N16O and 14N15N16O have very different isotopic fractionations in the stratosphere, and (ii) laboratory photolysis experiments conducted at 205 nanometers should better simulate the observed isotopic fractionation of stratospheric N2O. Modeling results indicate that there is no compelling reason to invoke a significant chemical source of N2O in the middle atmosphere and that individual N2O isotopomers might be useful tracers of stratospheric air parcel motion.

Analysis of aircraft exhausts with Fourier-transform infrared emission spectroscopy

Because of the worldwide growth in air traffic and its increasing effects on the atmospheric environment, it is necessary to quantify the direct aircraft emissions at all altitudes. In this study Fourier-transform infrared emission spectroscopy as a remote-sensing multi-component-analyzing technique for aircraft exhausts was investigated at ground level with a double pendulum interferometer and a line-by-line computer algorithm that was applied to a multilayer radiative transfer problem. Initial measurements were made to specify the spectral windows for traceable compounds, to test the sensitivity of the system, and to develop calibration and continuum handling procedures. To obtain information about the radial temperature and concentration profiles, we developed an algorithm for the analysis of an axial-symmetric multilayered plume by use of the CO(2) hot band at approximately 2400 cm(-1). Measurements were made with several in-service engines. Effects that were due to engine aging were detected but have to be analyzed systematically in the near future. Validation measurements were carried out with a conventional propane gas burner to compare the results with those obtained with standard measurement equipment. These measurements showed good agreement to within +/-20% for the CO and NO(x) results. The overall accuracy of the system was found to be +/-30%. The detection limits of the system for a typical engine plume (380 degrees C, ? = 50 cm) are below 0.1% for CO(2), ~0.7% for H(2)O, ~20 ppmv (parts per million by volume) for CO, and ~90 ppmv for NO.

Pressure sounding of the middle atmosphere from ATMOS solar occultation measurements of atmospheric CO(2) absorption lines

A method for retrieving the atmospheric pressure corresponding to the tangent point of an infrared spectrum recorded in the solar occultation mode is described and applied to measurements made by the Atmospheric Trace Molecule Spectroscopy (ATMOS) Fourier-transform spectrometer. Tangent pressure values are inferred from measurements of isolated CO(2) lines with temperature-insensitive strengths by measuring the slant-column CO(2) amount and by adjusting the viewing geometry until the calculated column matches the observed column. Tangent pressures are determined with a spectroscopic precision of l%-3%, corresponding to a tangent-point height precision of 70-210 m. The total uncertainty is limited primarily by the quality of the spectra and ranges between 4% and 6% (280-420 m) for spectra with signal-to-noise ratios of 300:1 and between 4% and 10% for spectra with signal-to-noise ratios of 100:1. The retrieval of atmospheric pressure increases the accuracy of the retrieved-gas concentrations by minimizing the effect of systematic errors introduced by climatological pressure data, ephemeris parameters, and the uncertainties in instrumental pointing.

Observations of the infrared solar spectrum from space by the ATMOS experiment

The final flight of the Atmospheric Trace Molecule Spectroscopy experiment as part of the Atmospheric Laboratory for Applications and Science (ATLAS-3) Space Shuttle mission in 1994 provided a new opportunity to measure broadband (625-4800 cm(-1), 2.1-16 µm) infrared solar spectra at anunapodized resolution of 0.01 cm(-1) from space. The majority of the observations were obtained as exoatmospheric, near Sun center, absorption spectra, which were later ratioed to grazing atmospheric measurements to compute the atmospheric transmission of the Earth's atmosphere and analyzed for vertical profiles of minor and trace gases. Relative to the SPACELAB-3 mission that produced 4800 high Sun spectra (which were averaged into four grand average spectra), the ATLAS-3 mission produced some 40,000 high Sun spectra (which have been similarly averaged) with an improvement in signal-to-noise ratio of a factor of 3-4 in the spectral region between 1000 and 4800 cm(-1). A brief description of the spectral calibration and spectral quality is given as well as the location of electronic archives of these spectra.

Measurements of Lorentz air-broadening coefficients and relative intensities in the H(2)(16)O pure rotational and v(2) bands from long horizontal path atmospheric spectra

Lorentz air-broadening coefficients and relative intensities have been measured for forty-three lines in the pure rotational band and twenty lines in the v(2) band of H(2)(16)O between 800 and 1150 cm(-1). The results were derived from analysis of nine 0.017-cm(-1) resolution atmospheric absorption spectra recorded over horizontal paths of 0.5-1.5 km with the McMath Fourier transform spectrometer and main solar telescope operated on Kitt Peak by the National Solar Observatory. A nonlinear least-squares spectral fitting technique was used in the spectral analysis. The results are compared with previous measurements and calculations. In most cases, the measured pressure-broadening coefficients and intensities are significantly different from the values in the 1986 HITRAN line parameters compilation.