VIRGO Data Products
VIRGO is an investigation for solar irradiance variability and helio-seismology on the ESA/NASA SOHO Mission
Table of Contents
- Instrumentation and Scientific Objectives
- VIRGO Data
- VIRGO Radiometry and Determination of TSI
- Detrending of SPM
For further information see the official virgo web site.
1. Instrumentation and Scientific Objectives
The VIRGO Experiment on the ESA/NASA SOHO Mission has three types of instruments to measure total solar irradiance (TSI) with the absolute radiometers DIARAD and PMO6V, spectral solar irradiance at 402, 500 and 862 nm with a 3-channel sun photometer (SPM) and spectral radiance at 500 nm with a luminosity oscillation imager (LOI).
- continuous high-precision, high-stability and high-accuracy measurements of the solar total and spectral irradiance and spectral radiance variation;
- frequencies, amplitudes and phases of oscillation modes in the frequency range of 0.001 to 8 mHz.
These data will be utilized to achieve the main scientific objectives of VIRGO summarized in the following list:
- detect and classify low-degree g modes of solar oscillations;
- determine the sound speed, density stratification and rotation in the solar interior, specifically determine the physical and dynamical properties of the solar core;
- study the solar atmosphere through comparison of amplitudes and phases of the p modes with these from velocity observations e.g. GOLF, SOI/MDI;
- search for the long periodicities or quasiperiodicities that have been found in other solar parameters;
- utilize the solar 'noise' signal to develop models for the global signature of stellar surface parameters;
- determine properties of the solar asphericity and its variation in time;
- study the relation between p-mode frequency changes and irradiance variations;
- study the influence of solar active regions and other large-scale structures on total and spectral irradiance;
- study the solar energy budget;
- provide accurate total and spectral irradiance data for input in terrestrial climate modelling.
The official virgo web site describes in detail the experiment, its instruments and science objectives. In the following we concentrate on data products and their evaluation .
Data Products and Evaluation
All VIRGO data are sent after reception on ground to the VIRGO Data Center which is operated by and located at the Instituto Astrofisica de Canarias in La Laguna on Tenerife. After reception the raw data are separated into instrument files forming the level-0 data sets organized in daily files. During this process possible gaps in the science data due to transmission problems are filled from the 12 and 24-hour delayed data packets, which are stored on board and transmitted together with real-time packets. Level-1 data are those that have been converted to physical units; they include the calibrations and contain all the corrections known a priori for instrument-related effects, such as the influence of temperature variation. The signals are also reduced to 1 AU distance and to zero radial velocity. The processing is based on the algorithms developed and maintained by the CoI responsible for the instrument concerned (Instrument CoIs: Th.Appourchaux for LOI, Seven Dewitte (D. Crommelynck until 1998) for DIARAD, I. Ruedi (J.Romero and M. Anklin until 1996 and 1999, respectively) for PMO6-V and Ch. Wehrli for SPM). These Level-1 data are publically available at the SOHO archive and are still raw in the sense that they do not contain corrections for, e.g. degradation, which can be calculated only from a posteriori analysis.
In a next step corrections for exposure and non-exposure dependent changes have to be determined. In the case of TSI the former ones can be deduced for each radiometer individually and yields two level-1.8 time series, whereas the latter is only accessible by comparison of the two. A combination of the two yields then the final VIRGO total solar irradiance (TSI) as explained in Section VIRGO Radiometry below. For the SPM data a similar procedure for the correction of the degradation is not possible, mainly because these instruments have much larger overall changes and the results from the rarely exposed back-up SPM do not provide a reliable source for the assessment of the long-term changes. Nevertheless the time series can be detrended and allow to assess periods of up to 1 year; they are labelled as level-2. These data (level-1.8 and 2) are available from our ftp server as described below.
Description and Availability of Level-1 Data Products
The level-1 data are organized in daily fits files with the data in ascii extensions. The data and format are described in the corresponding headers. As it is a standard FITS format, the files can be read with the FITS library for the programming language you use. An example of the FITS headers of SPM, DIARAD, PMO6V and LOI files for 1 January 2003 can be retrieved as FITS_Headers.txt from our ftp server; and there is also a collection of IDL routines for reading the VIRGO daily FITS files. The SPM, PMO6V and LOI data are 1-minute samples (1440 values per day) and the DIARAD data 3-minute samples (480 values per day). As MDI and GOLF, VIRGO should have had its one-minute values centered at 00:00:00 TAI; but early in the mission the midnight pulse which is used by VIRGO for synchronzation of the timing chain was set to 00:00:00 UTC instead of TAI, which yields a delay of 32 seconds (the sum of leap seconds since start of TAI). Thus the 1-minute averages of SPM and LOI are centered at 32s after the TAI minute (1st measurment @ 00:00:32 TAI). The 20-second averages of PMO6V every minute are centered at 22s after the TAI minute (1st measurment @ 00:00:18 TAI) and the 10-second averages of DIARAD every 3 minute at 87 seconds after the start of a 3 minute period (1st measurment @ 00:00:87 TAI).
From the data files the following information is retrieved by the idl routines:
- yymmdd_1.P00: 1440 values of PMO6V-A TSI in Wm-2 and a quality flag as 8-bit unsigned integer with the following meaning of the bit values: 64: cover closed; 32: SOHO off-pointing; 16: DIARAD calibration on (increases the noise of PMO6V-A, TSI values could be outliers), 8: period between reference values is more than the nominal 8 hours, values are less reliable, 4: smoothness condition failed, values are less reliable; 2: some temperatures from the HK packets are missing, values interpolated and essentially no problem for the data quality; 1: no distance correction. Contact: Wolfgang Finsterle
- yymmdd_1.D01: 480 values of DIARAD-L TSI in Wm-2 . The quality flag is not read as it contains status information and not the quality, the quality is reflected is either good if a value of TSI is available or bad if the value returned is NaN. Contact: Steven Dewitte
- yymmdd_1.S02: 1440 values of the red, green and blue channels in Wm-2nm-1 and one quality flag as 3-bit unsigned integer with the following meaning of the bit values: 4: no depointing correction could be applied (because of lack of the corresponding ancillary S/C data, this bit is always set, and the de-pointing angle set to zero); 2: no distance correction could be applied; 1: no temperature correction could be applied. Contact: Christoph Wehrli
- yymmdd_1.L00: 1440 values of the 12 main pixels (microamp) with one quality flag, 4 guiding pixels (microamp) with one quality flag, and the EW and NS diameter in arcmin. The quality flag is given as ascii 'GOOD' or 'POOR'. The following figure identifies the numbering for the LOI pixels. Contact: Thierry Appourchaux
The level-1 data products are normally updated on a weekly basis and transferred from VDC to the SOHO archive, from where they can be downloaded.
Description and Availability of Level-1.8 and 2 Data Products
The level 2 TSI time series for VIRGO, PMO6V and DIARAD are available as hourly and daily values. For a detailed description of the evaluation see VIRGO Radiometry below. Note that since version 5_005 the date of the update is added to the version. Since version 6_000 the algorithms for the evaluation have changed for the PMO6V evaluation from level 1 to 1.8 and for the corrections of DIARAD level 1.8 to 2.0. Since version 6_004 a new algorithm has been used to produce the level 2.0 from level 1.0 data. The actual VIRGO TSI is shown in Figure 2.1 which is also part of the composite TSI.
Figure 2.1: VIRGO level-2 time series. (pdf-Figure )
The level-2 data of the SPM are also available as hourly and daily values. A section describing the evaluation will be added in future; in the mean time the paper presented at SOHO-11 may be used as reference. The actual time series for the red, green and blue channels are plotted on Figure 2.2 together with TSI detrended with the same polynomials.
Figure 2.2: Detrended SPM Level-2 data . ( pdf-Figure )
Description and Availability of Data Products for Helioseismic Studies
Detrended and cleaned 1-minute time series for helioseismic studies are available for the SPM red, green and blue channel, the PMO6V TSI and the 12 pixels of LOI as idl save files. The SPM and LOI time series are detrended with a triangular filter with a 2-day base. The PMO6A time series is corrected for degradation (with the corresponding correction table ) and gaps shorter than 8 hours are filled with ARIMA extrapolations from both sides. These time series span somewhat different periods and are updated sporadically.
3. VIRGO Radiometry and Determination of TSI
From raw measurements to a reliable record of total solar irradiance (TSI)
The VIRGO Experiment has two types of radiometers to measure total solar irradiance (TSI): DIARAD and PMO6V. The former has been developed and built by IRMB, the latter by PMOD/WRC. For the VIRGO radiometers detailed descriptions of the principle and performance can be found in Fröhlich et al. (1995), Fröhlich et al. (1997). The scientific objective of VIRGO related the total solar irradiance, as stated in the original proposal and in Fröhlich et al. (1995), is to provide a reliable time series of TSI, which takes full advantage of both types of radiometers and their backup. This final product - the data level-2 - includes corrections for changes due to exposure to the Sun and other effects. The analysis of the data up to level-2 is now at version 6.4 and described in some detail in an aide-memoire VIRGO_TSI-vers64.pdf and summarized in the following.
From Level-0 to Level-1
The raw data (level-0) as transmitted from the spacecraft to ground are transformed into physical units using the radiometric factor which is determined for each individual radiometer by characterization experiments in the laboratory. Furthermore corrections for all a-priory known effects, are applied such as temperature, distance and radial velocity to the Sun. These data are labelled level-1. The algorithms for the calculation of level-1 data were developed by IRMB for DIARAD and by PMOD/WRC for the PMO6V radiometers. The original algorithms are described in Fröhlich et al. (1997) which are still being changed to cope with changes in operation and with the better understanding of their performance. Fig.3.1 shows the result for the present version, updated to YYMM, the date in the filename. These time series obviously do not show the variations of the solar irradiance alone, but also changes of the radiometric sensitivity which are quite different for the two types radiometers. It is also obvious that DIARAD shows less sensitivity changes than PMO6V and it is interesting to note that the decrease of sensitivity, normally termed as degradation, corresponds almost exactly to the one observed during the EURECA mission. In contrast DIARAD, which flew also on EURECA shows on SOHO much less change than then.
Figure 3.1: Level 1 data of the VIRGO radiometer: DIARAD-L, DIARAD-R, PMO6V-A and PMO6V-B.
( pdf-Figure )
From Level-1 to Level-2
The early behaviour of PMO6V is characterized by a rapid increase which is the result of a rapid darkening of the precision aperture surface which reduces the amount of stray-light from the baffle. The darkening is known since the apertures of the retrieved SOVA radiometer of the EURECA mission were inspected but the early increase was thought to be due to a change of the aperture heating. This, however, can never explain the effect and stray-light and IR radiation from the heated baffle back to the cavity is of the order of more than 3000ppm which can explain the early increase and also the high absolute values of TSI relative to TIM on SORCE. This effect is observed in most of the radiometer with the classical arrangement of the precison and view-limiting apertures as HF, ERBS and ACRIMs and so also their high absolute values. DIARAD does not show an early increase, mainly because the precision aperture is highly polished, electrolytically coated with nickel and curved to provide an image of the sun in the center of the view-limiting aperture, which avoids strai-light completely. So the high absolute value of DIARAD needs another explation and the reason was found in a substantially larger area covered by the electrical heater relative to the illuminated area which also reduces the absolute value by about 3000ppm. The revised characterization of both radiometers is described in VIRGO_Char2Space.pdf.
The new evaluation is no longer based on the straightforward correction of the PMO6V and DIARAD by their backup because PMO6V-B needs to be corrected for its early increase before it can be used to correct PMO6V-A. So, we start with the correction DIARAD-L by R and then compare the corrected DIARAD with PMO6V-B by fitting an exponential correction for the non-exposure dependent change of DIARAD and a dose and temperature dependent correction for the early increase of PMO6V-B simulataneously. This yields a fully corrected DIARAD-L which corresponds to the final level-2 and a corrected PMO6V-B which can now be used to correct PMO6V-A to level-2. In the following the different steps are illustrated.
The first step is the comparison of DIARAD-L with the backup R, illustrated in Fig. 3.2. In contrast to earlier versions the interpolation of the monthly values of DIARAD L/R is now done with splines which reveals more details and shows that also the DIARAD degradation shows a temperature dependence with its annual variation on SOHO. The correction is fundmentally different from Crommelynck et al. (2004) and the IRMB/DIARAD page). DIARAD needs also a few other corrections due to switch-off of VIRGO. They have been re-determined and showed also slips and the first observed exponential recovery is not always found after later switch-off.
Figure 3.2: The correction from the individual comparisons of DIARAD-L and R is shown together with former results. Also shown are the corrections for the changes after switch-off of VIRGO, which show also slips in contrast to the earlier corrections.
Now we have a corrected DIARAD which can be compared to PMO6V-B as shown in the top panel of Fig.3.3 together with the fitted curve which is based on an exponential function for DIARAD and a dose and temperature dependent increase for PMO6V-B. Then the change over the SOHO vacations is determined by comparison with ACRIM-I during 500 days (middle panel). If the so determined difference is applied, it turns out that the full difference is due to DIARAD. A posteriori we can now assume that the PMO6V-B does not change its sensitivity over the SOHO vacations as shown in the bottom panel.
Figure 3.3: The top panel shows the difference between the corrected DIARAD and PMO6V-B and the fit with two parameters for the exponential function of DIARAD, one parameter for the step over the gap and eight parameters for PMO6V-B (Eq.11 of VIRGO_TSI-vers64.pdf). The middle panel shows the the comparison with ACRIM-I and bottom panel the final result. Note, the zero change of PMO6V-B over the SOHO vacations.
Before we can use the corrected PMO6V-B to correct PMO6V-A we need to correct PMO6V-A for the early increase. The main part of the increase is during the shuttered operation which provides a reference to the original absolute calibration and later the change in calibration due the change in operation. For this early increase correction we need TSI values to compare with. As DIARAD starts at mission day 68 only, we expand these data with the proxy model to mission day 48. The result is shown in Fig.3.4.
Figure 3.4: Top panel: Ratio of the irradiance of PMO6V-A to a combination of the proxy model and DIARAD, which starts only at 7 February 1996. Apart from the correction due to the early increase, the change of absolute values due to change from shuttered to covered operation is also determined. Bottom panel: Corrected PMO6V-A values 'with new corrections´ compared to the former ones as 'with old corrections´. Also shown are the shifted values as used in the further analysis which refer the absolute values to the first-light value.
With the level-2 PMO6V-B series the corrections for PMO6V-A can now be determined over the whole mission and interpolated by the method of splines to the hourly values as shown in Fig.3.5 (bottom panel) and the change over the SoHO interruption is determined by comparison with ACRIM (top panel). Again, the PMO6V-A has not changed during SOHO vacations similar to PMO6V-B and we could indeed do the analysis without using ACRIM-I for comparison.
Figure 3.5: The bottom panel shows the corrections for both versions, the 'old´ 6.2 and the 'new´ 6.4 of PMO-A from comparison with the corrected PMO-B and the top panel the comparison with ACRIM for the determination of the change over the gap.
VIRGO TSI is defined as the average of the PMO6V and DIARAD values and is determined by weighting. The weights are deduced from the difference of the variances of PMO6V minus the one of DIARAD, each determined from a 81-day-running period and shown by the green line on the top panel of Fig.3.6, smoothed with a 131-day boxcar. It is normalized in such a way that the absolute maximum is set to 0.5 and a positive difference means that the noise of PMO6V is higher and its weight correspondingly lower and a negative one that the noise of DIARAD is higher. So, weights from 0...1 are shown as red and blue dashed lines (with the left scale). On the other hand, DIARAD has in general less noise than PMO6V with e.g. the spikes during the keyholes due to the fact that still the original level-0 to 1 algorithm. For the moment, it seemed more important to have an internally consistent way to produce the level-2 data, than improving on the noise of PMO6V. The new evaluation of level-0 to 1 is under way and will be available in due time.
Figure 3.6: The bottom panel shows the corrections for both versions, the 'old´ 6.2 and the 'new´ 6.4 of PMO-A from comparison with the corrected PMO-B and the top panel the comparison with ACRIM for the determination of the change over the gap.
This leads to the final VIRGO TSI record as
hourly and daily TSI. These data files have a similar format as before, the first TSI column is now VIRGO TSI on the new absolute scale (see next Section) and the others as before the VIRGO, PMO6V and DIARAD TSI values on the original scale. Please note that we have not changed the version, because a scale change is only a factor and not related to the way the data are analysed.
Now we can compare the VIRGO TSI with ACRIM-3 and TIM as shown in Fig.3.7. Besides the absolute value, the major results can be summarized as follows:
(i) Relative to TIM VIRGO TSI has a downward trend of $13.1 \pm 1.9$\,ppm/a corresponding to a difference of about 150\,ppm over a solar cycle.
(ii) Relative to ACRIM-3 VIRGO TSI has during the period of TIM an upward trend of 20.2 +/- 1.9ppm/a.
(iii) VIRGO TSI has still the artefact during the keyholes of SOHO which are coming from the PMO6V data due to some inadequate temperature interpolation of missing values in the present level0 to 1 evaluation. In general the PMO6V values have less medium- to long-term deviations relative to TIM than DIARAD, but the noise of PMO6V is higher.
Figure 3.7: VIRGO TSI (original scale) is compared to ACRIM and TIM. The ACRIM values are on the scale of ACRIM\,3 (version 1311). It is interesting to note that the slope of VIRGO relative to ACRIM for the whole period and for the TIM period are similar although the time series is composed of ACRIM\,2 and 3.
If you want to use the 1-minute data of PMO6V or the 3-minute data for DIARAD, you can fetch them from the SOHO archive
and then you can use the
correction table to get the level-2 product. The file may also be taken from ftp://ftp.pmodwrc.ch/pub/data/irradiance/virgo/1-minute_Data/VIRGO_1min_6004_1410.idl, updated until October 2014 as an IDL save file.
Revised Characterization and new absolute value of TSI
Since the advent of the results of TIM on SORCE it became clear that the classical radiometers may have a problem with the absolute scale. From the beginning it was suspected that it may have to do with the relative arrngement of the precision and view limiting apertures. The classical radiometer have the precision aperture directly in front of the cavity whereas the new generation of radiometers like TIM have it at the very entrance of the radiometer and use the cavity as a view-limiting aperture. The main difference is that in the former arrangement sunlight illuminates the precision aperture and reflects the light into the baffle, which in turn produces extra radiation in the cavity as strailight and infrared radiation due to local heating. This effect was confirmed by investigating spare radiometers of ACRIM and PMO6V at the TRF and a corresponding re-evaluation of the characterization of the PMO6V radiometers has been performed which is described in detail in a aide-memoire VIRGO_Char2Space
. The situation with DIARAD is different and the higher value is not due to the illumination of the porecision aperture, but due to the fact that the electrical heater does not cover the same area as the radiation, which produces a large non-equivalence not taken into account. A detailed calculation of this effect is also presented in the aide-memoire. The final result of this re-evaluation of the characterization of both type of radiometers yields an absolute value of VIRGO TSI which is about 5.8 Wm-2 or 0.42% lower than the original scale. During the last minimum VIRGO TSI on the new scale amounts to 1359.66 ± 2.47 Wm-2 which is within less than half the uncertainty equal to the value of TIM of 1360.52 Wm-2.
Anklin, M., Fröhlich, C., Finsterle, W. Crommelynck, D. A. and Dewitte, S.; 1998, Assessment of Degradation of VIRGO Radiometers Onboard SOHO, Metrologia 35, pp. 685-688, 1999.
Crommelinck, D., Dewitte, S. and Chevalier, A.: 2004, Total Solar Irradiance from VIRGO on SOHO: 5 Years of Operation of DIARAD, submitted to JGR
Fröhlich, C., Romero, J., Roth, H., Wehrli, C., Andersen, B. N., Appourchaux, T., Domingo, V., Telljohann, U., Berthomieu, B., Delache, P., Provost, J., Toutain, T., Crommelynck, D., Chevalier, A., Fichot, A., Däppen, W., Gough, D. O., Hoeksema, T., Jiménez, Gómez, M., Herreros, J., Roca-Cortés, T., Jones, A. R., Pap, J.and Willson, R. C.: 1995, VIRGO: Experiment for Helioseismology and Solar Irradiance Monitoring, Sol. Phys. 162, 101-128.
Fröhlich, C., Crommelynck, D., Wehrli, C., Anklin, M., Dewitte, S., Fichot, A., Finsterle, W., Jiménez, A., Chevalier, A. and Roth, H. J.; 1997, In-Flight Performances of VIRGO Solar Irradiance Instruments on SOHO, Sol. Phys. 175, 267-286.
Fröhlich, C. and Finsterle, W.: 2001, VIRGO Radiometry and Total Solar Irradiance 1996-2000 Revised, in: Recent Insights into the Physics of the Sun and Heliosphere: Highlights from SOHO and other Space Missions, Eds: P. Brekke, B. Fleck, J.B. Gurman, Astr.Soc.Pacific, IAU Symp. 203, 105-110
Fröhlich, C.: 2003, Long-Term Behaviour of Space Radiometers, Metrologia
Responsible for this page is Claus Fröhlich ; last update of text (without references): October 2014
<< overview space projects