
Product Quality of Climate Data Records

1D Monthly Profiles  Product Service Spec  2D Time Series  Scalar Time Series  Vertical Mean TS  Plot Description  Switch to Daily 
The 1D profile plots show differences between the retrieved GNSSRO variables and equivalent values calculated using fields from ERAI, ERA5 or operational ECMWF data. Comparisons with both model analysis (A) and model forecast (B) are made available.
The variables temperature, dry temperature, pressure and specific humidity are compared in terms of difference between absolute values, and refractivity, bending angle, optimized bending angle, pressure and specific humidity are compared in terms of their fractional differences.
In addition a small suite of diagnostic plots related to performance of 1DVar are made available at the end of the drop down menu. These are the mean adjustment in refractivity space, the 4 Desrozier ratios and the solution error relative to the a priori error.
Profiles are interpolated (cubic spline interpolation) to a fixed altitude grid following this (quadratic) form: [0 2 4 6 .. 284]^{2} m, i.e. [0, 4,16,36,64,100,144 .. 79524] m, (143 levels). Impact height, geometric altitude and geopotential height are is used as independent vertical variable for Level 1B, Level 2A and Level 2B respectively.
The fractional differences in bending angle, refractivity and pressure are calculated as (XY)/Y, where Y refer to either analysis or forecast from ERAI, ERA5 or operational ECMWF data. These plots are called "relative". Due to relatively large variability in specific humidity, the fractional differences for specific humidity are calculated as (XY)/mean(B), where B (background) refers to the ERAI, ERA5 or ECMWF operational forecast, i.e., the specific humidity difference is normalized with respect to the monthly or daily mean of the forecast. These plots are called "normalized". In this particular calculation the background specific humidity (B) is interpolated linearly to the common grid instead of using splines.
For all three plot types  absolute, relative and normalized  the monthly or daily statistics  mean, standard deviation, median and median absolute deviation  are calculated after interpolation and possible division with reference profile.
Solid black lines show mean values, dashed black lines show standard deviations, solid red lines show median and dashed red lines show the median absolute deviation (MAD). MAD(x) = k median(xmedian(x)), where k = 1.4826 is chosen such that the MAD equals the standard deviation for normally distributed data. For daily statistics, the blue lines show individual profile differences (spaghetti plots)
The plot on the right shows the percentage of data points valid at a given height. The number of occultations used in the statistics for that month or day, i.e. number that were accepted by the quality control is printed in the title.
The schematic graphs show with a color code whether different product specifications have been violated during a given month and altitude interval, for a given variable and satellite mission. If an area is hatched it means that the service specifications for the variable were violated on average during the corresponding period. The profiles are divided into latitudinal zones: tropics (30° S30° N), midlatitudes (30° S60° S and 30° N60° N), and polar regions (60° S90° S and 60° N90° N). Each zone is further separated into height layers, which depend on the product type. Monitoring of the tropopause height performance is found the scalar timeseries.
Compliance with specifications:
For each variable and day we calculate the mean and standard deviation as follows: (i) Bending angle and refractivity: "observation  background", (ii) 1DVar products: "solution  analysis". If monthly STD value exceeds the product requirements or service specification at one or more altitudes, the color code for that violation is used for the corresponding altitude range.
Product requirements and service specifications are calculated in many different ways, based on mean and standard deviations.These are summarized in the PRD and SeSp documents found at https://romsaf.eumetsat.int/product_documents.php. Three product requirement accuracy levels are specified: threshold, target, and optimal. The color codes seen in the color bar indicate that the given requirements are met, i.e., not violated; e.g., the color yellow means that threshold level is met but the target level has been violated.
These plots show monthly or daily standard deviation or mean of the chosen variable comparison, satellite and area. The plots include all data from the selected RO mission that passed the QC screening (nominal data). The monthly or daily STD and mean used to produce these plots are exactly the same as those used in "1D Profiles". The standard deviation plots of temperature, humidity and pressure, compared to ERAI, ERA5 or ECMWF operational analyses, are overlaid with contours enclosing areas where requirements (stippled areas enclosed by solid lines) or service specifications (hatched areas enclosed by purple lines) are violated. The requirements and specifications are summarized in the PRD and SeSp documents.
The scalar time series plots show differences between retrieved GNSSRO scalar variables and equivalent values calculated using fields from ERAI, ERA5 or operational ECMWF data. Comparisons are made available in 8 different latitudinal zones. For the tropopause an additional narrow tropical zone (15° S15° N) is also monitored.
The available variables, surface pressure, tropopause height and scaled cost function, are of different kinds and therefore treated differently:
The pressure is compared in terms of absolute difference with both model analysis (A) and model forecast (B). Both standard deviation and mean are shown in the same plots.
The tropopause height and temperature are viewed in separate mean and standard deviation plots (that is 4 plot types). Three tropopause heights are compared to ERAI, ERA5, or ECMWF operational analysis: Refractivity based, dry temperature cold point based and dry temperature lapse rate based. For tropopause temperature there is only two types presented (dry temperature cold point and dry temperature lapse rate).
Instead of evaluating the difference between all tropopause estimates and all model tropopause estimates, we show the difference between different model based tropopause estimates and ERAI, ERA5 or ECMWF operational lapse rate tropopause in the "Model Lapse Rate Tropopause" and "Model Cold Point Tropopause" plots.
The scaled cost function at the end of the left drop down menu is a diagnosis of how well the 1DVar retrieval of temperature, specific humidity and surface pressure performs, in a well tuned 1DVar its value should be close to unity. For more information, see the 1DVar ATBD which is found here https://romsaf.eumetsat.int/product_documents.php.
The vertical mean plots are designed to monitor the Level 1B and Level 2A products with respect to the service specifications, which are defined as vertical mean values for (OB) bias and standard deviation (relative to B for bending angle and refractivity; in kelvin for dry temperature). Here O is the retrieval and B is the ERAI, ERA5, or ECMWF operational forecast. Bias should be understood as the mean value of the absolute value of the bias. The bias difference between rising and setting occultations (R vs S) is the absolute value of the difference between rising and setting absolute biases. Service Specifications are marked with red dashed line (only for single missions  not for combined data sets).