Level 2 Precision Analysis


The precision of OSIRIS Level 2 data is a measure of the uncertainty in the values retrieved from OSIRIS measurements of light scattering in the atmosphere.  This uncertainty comes from measurement noise and random noise generated through the SASKTRAN forward model in the SaskMART retrieval algorithm.  Using other retrieval methods the precision information is generated as a biproduct naturally in the form of a covariance matrix but by its nature SaskMART does not. The covariance matrix can be obtained after the retrieval using the original measrements' uncertainty and the retrieved result. This method has been validated using matching pair statics.  In this method stratospheric ozone in the tropics is considered to be stable over the timescale of a day so retrieved ozone profiles in close geographical proximity and within 1.01 days were compared.  The variation in these matching pairs compared favorably with the precision values.  This is the expected result because the precision should be a measure of the random error in the retrievals.  

Precision Value Retrieval


  This is a simplified explanation of the precision retrieval, for a more detailed explanation please see the JGR-Atmospheres article by A.E. Bourassa <link to final article once published>. It is difficult to measure things like atmospheric ozone and aerosol density directly because rocket and balloon measurements are expensive and localized.  Because of this we measure these things indirectly from light that is scattered from the atmosphere.  Modelling how light is scattered based on a given set of atmospheric conditions is difficult but possible (this is done using a "forward model" like SASKTRAN) .  The inverse problem is more difficult because atmospheric radiative transfer is not a linear problem so the equations are not invertable.  

Methods like SaskMART are numeric itterative solutions that use the forward model, and the observations to retrieve the desired state parameter like ozone. SaskMART takes an initial guess at the atmospheric state, runs it through SASKTRAN to predict the observed radiance at a set of wavelengths.  The guess is then updated based on the difference between the predicted radiance, and the actual OSIRIS measurements.  During this process a weighting matrix 'W' is used to determine how much the differences in the observation at one altitude affect the state parameter at other altitudes and how important each wavelength is.  We have uncertainty values associated with the OSIRIS measurements, and from this we can make a covariance matrix for each wavelength (which is specified by the index 'i'): 

The degree to which the uncertainty in our retrieved state parameter (like ozone) affects the predicted radiance given by OSIRIS is described by a kernel matrix where an element at row 'l' and column 'm' is:  

Using these things we can calculate the intermediate covariance of our retrieved state parameter based on each wavelength 'i':

To get elements of the total covariance matrix of the retrieved state parameter we use the weighting matrix 'W' from SaskMART:

Validation of Precision - Matching Pair Statistics



The obtained precision values were validated following a method presented by Piccolo and Dudhia (2007). It is generally accepted that at low latitudes ozone concentrations are relatively stable above the tropopause. Based on this we assume, as did the study by Toohey et al. (2010), that two measurements of low latitude stratospheric ozone taken a day apart should be approximately the same.  Any difference between the profiles could be attributed to random measurement and retrieval noise.   We assemble a set of the differences between matched pairs and find the variance of these differences.  This variance is related to  the random error associated with a profile.  

Mathematically we can say that the difference between two profiles is  so we find the variance in the set of differences:

For a large set of matched pairs the mean uncertainty in the profiles is

and is related to the variance of differences by the equation   

or we can write

and this value should match the average precision value.  

They are compared in the figure below using 


For this study we considered scans to be a matched pair if they were taken within a 1.01 day window, and were within 1.5° longitude  and 0.5° latitude of each other.

Sample Results of Precision Validation



This figure demonstrates that estimated precision (represented by the "Percent Mean Uncertainty") agrees well with the deviation between what should approximately represent measurements of identical atmospheres (represented by "Percent Standard Deviation").  This supports the precision retrieval process.

At lower altitudes (at and below the tropopause) the atmosphere varies too greatly over the course of 24 hours for the Percent Standard Deviation to represent the uncertainty.  This is why the red stars no longer match the grey area at low altitudes.

Precision Retrieval in the Processing Chain


The in the OSIRIS MART retrieval software a C++ class called MART_OSIRISLevel2 acts as the main driver of events.  Given a set of initial conditions (such as a specific OSIRIS scan and the desired output species) this class sets up the system and coordinates the operations necessary for the retrieval.  To carry out the precision retrievals two new functions were added called ErrorAnalysis_HybridAerosol and ErrorAnalysis_O3.  This level of the program handles the OSIRIS-specific operations such as obtaining the necessary OSIRIS Level 1 data.

When these functions are called they obtain critical information from other components of the program (species profile, observed radiance, and uncertainty in observed radiance) and call MARTRetrieval::RetrieveError.  The MARTRetrieval class uses this information to perform the precision retrieval that is outlined in previous pages.  This level of the program is designed to operate independant of the source of the input data.  This means that the RetrieveError function could be easily applied to other similar problems regardless of the source of the data.

To completely retrieve ozone and aerosol along with their precisions the following parameters are retrieved in this order:

  1. NO2
  2. Albedo
  3. Aerosol
  4. Albedo
  5. Aerosol
  6. Aerosol Precision
  7. O3
  8. O3 Precision