Fast radiative transfer using monochromatic look-up tables

AbstractLine-by-line (LBL) methods of numerically solving the equations of radiative transfer can be inhibitingly slow. Operational trace gas retrieval schemes generally require much faster output than current LBL radiative transfer models can achieve. One option to speed up computation is to precalculate absorption cross sections for each absorbing gas on a fixed grid and interpolate. This work presents a general method for creating, compressing, and validating a set of individual look-up tables (LUTs) for the 11 most abundant trace gases to use the Reference Forward Model (RFM) to simulate radiances observed by the Infrared Atmospheric Sounding Interferometer (IASI) at a more operational pace. These LUTs allow the RFM to generate radiances more than 20 times faster than LBL mode and were rigorously validated for 80 different atmospheric scenarios chosen to represent variability indicative of Earth׳s atmosphere. More than 99% of all IASI simulated spectral channels had LUT interpolation errors of brightness temperature less than 0.02K, several factors below the IASI noise level. Including a reduced spectral grid for radiative transfer speed up the computation by another factor of six at the expense of approximately doubling interpolation errors, still factors below IASI noise. Furthermore, a simple spectral compression scheme based upon linear interpolation is presented, which reduced the total LUT file size from 120Gbytes to 5.6Gbytes; a compression to just 4.4% of the original. These LUTs are openly available for use by the scientific community, whether using the RFM or to be incorporated into any forward model

ISAMS observations of stratospheric aerosol

The Improved Stratospheric and Mesospheric (ISAMS) on board the Upper Atmosphere Research Satellite (UARS) incorporates a 12.1 micron window channel for the measurement of aerosol opacity. The retrieval scheme is discussed briefly and preliminary observations of the Mt. Pinatubo aerosol cloud are presented and compared with SAGE 2 observations at 1.02 microns. The effect of aerosol on other ISAMS channels and its spectral dependence is discussed

Measurements of stratospheric constituents by ISAMS

ISAMS is a limb sounding radiometer flying on the UARS, and designed to measure temperature, pressure, O3, CO, NO, NO2, N2O5, HNO3, CH4, H2O, N2O, and aerosol. Its capabilities are described, together with the present status of validation of its data products, and plans for future improvement

Effects of stratosphere-troposphere chemistry coupling on tropospheric ozone

A new, computationally efficient coupled stratosphere-troposphere chemistry-climate model (S/T-CCM) has been developed based on three well-documented components: a 64-level general circulation model from the UK Met Office Unified Model, the tropospheric chemistry transport model (STOCHEM), and the UMSLIMCAT stratospheric chemistry module. This newly developed S/T-CCM has been evaluated with various observations, and it shows good performance in simulating important chemical species and their interdependence in both the troposphere and stratosphere. The modeled total column ozone agrees well with Total Ozone Mapping Spectrometer observations. Modeled ozone profiles in the upper troposphere and lower stratosphere are significantly improved compared to runs with the stratospheric chemistry and tropospheric chemistry models alone, and they are in good agreement with Michelson Interferometer for Passive Atmospheric Sounding satellite ozone profiles. The observed CO tape recorder is also successfully captured by the new CCM, and ozone-CO correlations are in accordance with Atmospheric Chemistry Experiment observations. However, because of limitations in vertical resolution, intrusion of CO-rich air in the stratosphere from the mesosphere could not be simulated in the current version of S/T-CCM. Additionally, the simulated stratosphere-to-troposphere ozone flux, which controls upper tropospheric OH and O3 concentrations, is found to be more realistic in the new coupled model compared to STOCHEM. © 2010 by the American Geophysical Union

The ESA MIPAS/Envisat level2-v8 dataset: 10 years of measurements retrieved with ORM v8.22

The observations acquired during the full mission of the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) instrument, aboard the European Space Agency Environmental Satellite (Envisat), have been analysed with version 8.22 of the Optimised Retrieval Model (ORM), originally developed as the scientific prototype of the ESA level-2 processor for MIPAS observations. The results of the analyses have been included into the MIPAS level-2 version 8 (level2-v8) database containing atmospheric fields of pressure, temperature, and volume mixing ratio (VMR) of MIPAS main targets H$_{2}$O, O$_{3}$, HNO$_{3}$, CH$_{4}$, N$_{2}$O, and NO$_{2}$, along with the minor gases CFC-11, ClONO$_{2}$, N$_{2}$O$_{5}$, CFC-12, COF$_{2}$, CCl$_{4}$, CF$_{4}$, HCFC-22, C$_{2}$H$_{2}$, CH$_{3}$Cl, COCl$_{2}$, C$_{2}$H$_{6}$, OCS, and HDO. The database covers all the measurements acquired by MIPAS in the nominal measurement mode of the full resolution (FR) part of the mission (from July 2002 to March 2004) and all the observation modes of the optimised resolution (OR) part (from January 2005 to April 2012). The number of species included in the MIPAS level2-v8 dataset makes it of particular importance for the studies of stratospheric chemistry. The database is considered by ESA the final release of the MIPAS level-2 products. The ORM algorithm is operated at the vertical grid coincident to the tangent altitudes of the observations or to a subset of them, spanning (in the nominal mode) the altitude range from 6 to 68 km in the FR phase and from 6 to 70 km in the OR period. In the latitude domain, FR profiles are spaced by about 4.7∘, while the OR profiles are spaced by about 3.7∘. For each retrieved species, the auxiliary data and the retrieval choices are described. Each product is characterised in terms of the retrieval error, spatial resolution, and “useful” vertical range in both phases of the MIPAS mission. These depend on the characteristics of the measurements (spectral and vertical resolution of the measurements), the retrieval choices (number of spectral points included in the analyses, number of altitudes included in the vertical retrieval grid), and the information content of the measurements for each trace species. For temperature, water vapour, ozone, and nitric acid, the number of degrees of freedom is significantly larger in the OR phase than in the FR one, mainly due to the finer vertical measurement grid. In the FR phase, some trace species are characterised by a smaller retrieval error with respect to the OR phase, mainly due to the larger number of spectral points used in the analyses, along with the reduced vertical resolution. The way of handling possible caveats (negative VMR, vertical grid representation) is discussed. The quality of the retrieved profiles is assessed through four criteria, two providing information on the successful convergence of the retrieval iterations, one on the capability of the retrieval to reproduce the measurements, and one on the presence of outliers. An easy way to identify and filter the problematic profiles with the information contained in the output files is provided. MIPAS level2-v8 data are available to the scientific community through the ESA portal (https://doi.org/10.5270/EN1-c8hgqx4)

Level 2 processor and auxiliary data for ESA Version 8 final full mission analysis of MIPAS measurements on ENVISAT

High quality long-term data sets of altitude-resolved measurements of the atmospheric composition are important because they can be used both to study the evolution of the atmosphere and as a benchmark for future missions. For the final ESA reprocessing of MIPAS (Michelson Interferometer for Passive Atmospheric Sounding) on ENVISAT (ENViromental SATellite) data, numerous improvements were implemented in the Level 2 (L2) processor Optimised Retrieval Model (ORM) version 8.22 (V8) and its auxiliary data. The implemented changes involve all aspects of the processing chain, from the modelling of the measurements with the handling of the horizontal inhomogeneities along the line of sight to the use of the optimal estimation technique to retrieve the minor species, from a more sensitive approach to detecting the spectra affected by clouds to a refined method for identifying low quality products. Improvements in the modelling of the measurements were also obtained with an update of the used spectroscopic data and of the databases providing the a priori knowledge of the atmosphere. The HITRAN_mipas_pf4.45 spectroscopic database was finalised with new spectroscopic data verified with MIPAS measurements themselves, while recently measured cross-sections were used for the heavy molecules. The Level 2 Initial Guess (IG2) data set, containing the climatology used by the MIPAS L2 processor to generate the initial guess and interfering species profiles when the retrieved profiles from previous scans are not available, was improved taking into account the diurnal variation of the profiles defined using climatologies from both measurements and models. Horizontal gradients were generated using the ECMWF ERA-Interim data closest in time and space to the MIPAS data. Further improvements in the L2 V8 products derived from the use of the L1b V8 products, which were upgraded to reduce the instrumental temporal drift and to handle the abrupt changes in the calibration gain. The improvements introduced into the ORM V8 L2 processor and its upgraded auxiliary data, together with the use of the L1b V8 products, lead to the generation of the MIPAS L2 V8 products, which are characterised by an increased accuracy, better temporal stability and a greater number of retrieved species

Level 2 processor and auxiliary data for ESA Version 8 final full mission analysis of MIPAS measurements on ENVISAT

High quality long-term data sets of altitude-resolved measurements of the atmospheric composition are important because they can be used both to study the evolution of the atmosphere and as a benchmark for future missions. For the final ESA reprocessing of MIPAS (Michelson Interferometer for Passive Atmospheric Sounding) on ENVISAT (ENViromental SATellite) data, numerous improvements were implemented in the Level 2 (L2) processor Optimised Retrieval Model (ORM) version 8.22 (V8) and its auxiliary data. The implemented changes involve all aspects of the processing chain, from the modelling of the measurements with the handling of the horizontal inhomogeneities along the line of sight to the use of the optimal estimation technique to retrieve the minor species, from a more sensitive approach to detecting the spectra affected by clouds to a refined method for identifying low quality products. Improvements in the modelling of the measurements were also obtained with an update of the used spectroscopic data and of the databases providing the a priori knowledge of the atmosphere. The HITRAN_mipas_pf4.45 spectroscopic database was finalised with new spectroscopic data verified with MIPAS measurements themselves, while recently measured cross-sections were used for the heavy molecules. The Level 2 Initial Guess (IG2) data set, containing the climatology used by the MIPAS L2 processor to generate the initial guess and interfering species profiles when the retrieved profiles from previous scans are not available, was improved taking into account the diurnal variation of the profiles defined using climatologies from both measurements and models. Horizontal gradients were generated using the ECMWF ERA-Interim data closest in time and space to the MIPAS data. Further improvements in the L2 V8 products derived from the use of the L1b V8 products, which were upgraded to reduce the instrumental temporal drift and to handle the abrupt changes in the calibration gain. The improvements introduced into the ORM V8 L2 processor and its upgraded auxiliary data, together with the use of the L1b V8 products, lead to the generation of the MIPAS L2 V8 products, which are characterised by an increased accuracy, better temporal stability and a greater number of retrieved species

The SPARC water vapour assessment II: Profile-to-profile comparisons of stratospheric and lower mesospheric water vapour data sets obtained from satellites

This work is distributed under the Creative Commons Attribution 4.0 License. Within the framework of the second SPARC (Stratosphere-troposphere Processes And their Role in Climate) water vapour assessment (WAVAS-II), profile-to-profile comparisons of stratospheric and lower mesospheric water vapour were performed by considering 33 data sets derived from satellite observations of 15 different instruments. These comparisons aimed to provide a picture of the typical biases and drifts in the observational database and to identify data-set-specific problems. The observational database typically exhibits the largest biases below 70 hPa, both in absolute and relative terms. The smallest biases are often found between 50 and 5 hPa. Typically, they range from 0.25 to 0.5 ppmv (5 % to 10 %) in this altitude region, based on the 50 % percentile over the different comparison results. Higher up, the biases increase with altitude overall but this general behaviour is accompanied by considerable variations. Characteristic values vary between 0.3 and 1 ppmv (4 % to 20 %). Obvious data-set-specific bias issues are found for a number of data sets. In our work we performed a drift analysis for data sets overlapping for a period of at least 36 months. This assessment shows a wide range of drifts among the different data sets that are statistically significant at the 2σ uncertainty level. In general, the smallest drifts are found in the altitude range between about 30 and 10 hPa. Histograms considering results from all altitudes indicate the largest occurrence for drifts between 0.05 and 0.3 ppmv decade-1. Comparisons of our drift estimates to those derived from comparisons of zonal mean time series only exhibit statistically significant differences in slightly more than 3 % of the comparisons. Hence, drift estimates from profile-to-profile and zonal mean time series comparisons are largely interchangeable. As for the biases, a number of data sets exhibit prominent drift issues. In our analyses we found that the large number of MIPAS data sets included in the assessment affects our general results as well as the bias summaries we provide for the individual data sets. This is because these data sets exhibit a relative similarity with respect to the remaining data sets, despite the fact that they are based on different measurement modes and different processors implementing different retrieval choices. Because of that, we have by default considered an aggregation of the comparison results obtained from MIPAS data sets. Results without this aggregation are provided on multiple occasions to characterise the effects due to the numerous MIPAS data sets. Among other effects, they cause a reduction of the typical biases in the observational database

Analysis of new species retrieved from MIPAS

The Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) instrument which operated on the Envisat satellite from 2002-2012 is a Fourier transform spectrometer for the measurement of high-resolution gaseous emission spectra at the Earth's limb. It operates in the near- to mid-infrared, where many of the main atmospheric trace gases have important emission features. The initial operational products were profiles of Temperature, H2O, O3, CH4, N2O, HNO3, and NO2, and this list was recently extended to include N2O5, ClONO2, CFC-11 and CFC-12. Here we present preliminary results of retrievals of the third set of species under consideration for inclusion in the operational processor: HCN, CF4, HCFC-22, COF2 and CCl4