next previous
Up: Active region effects

3. Calibration in velocity

The passage from the photometric ratio r to a velocity is not trivial, but requires a calibration procedure which has been argument of many papers reporting oscillation measures from resonant spectrometers (see e.g. van der Raay et al. 1985; Pallé et al. 1993).

3.1. Standard procedure

The calibration is based on a best fitting procedure of the relation r = f(V), where V represents a total line of sight velocity. The difficulties arise since the velocity V is partly unknown and partly spurious. In effect, according to Pallé et al. (1993), we can write
displaymath1305
where tex2html_wrap_inline1317 is the unknown solar oscillation velocity, tex2html_wrap_inline1319 is an unknown solar offset, and V0 is the known part of V. tex2html_wrap_inline1319 includes the integrated signal produced by the active region as well as by the solar velocity fields at different scales. Pallé et al. claim that tex2html_wrap_inline1319 is roughly constant during a 12-hour observation (tex2html_wrap_inline1329) while it changes by as much as tex2html_wrap_inline1331 in 7 days. In the case of ground-based measurements, V0 is
displaymath1306
where the first two terms are due to Earth-Sun relative motion, tex2html_wrap_inline1335 is the signal produced by the Earth atmosphere differential extinction, and tex2html_wrap_inline1337 is the gravitational redshift in velocity unit (tex2html_wrap_inline1339). In the case of the GOLF measurements, V0 is
displaymath1307
where tex2html_wrap_inline1343 is the component relative to the satellite halo orbit around the Lagrangian point tex2html_wrap_inline1345, and tex2html_wrap_inline1347 is the same as for ground-based observation with the approximation that the Lagrangian point is placed at the Earth position. In both cases, V0 can be calculated with great accuracy.

The calibration procedure developed by Pallé et al. (1993) for the IRIS data consists of two steps: i) first a table tex2html_wrap_inline1351 vs. V0 was built, where tex2html_wrap_inline1351 is the mean of r over tex2html_wrap_inline1359 bins of V0, in which the variations of tex2html_wrap_inline1317 and tex2html_wrap_inline1319 are minimized, and an average calibration function tex2html_wrap_inline1367 was determined by fitting the function tex2html_wrap_inline1369 to the table tex2html_wrap_inline1351 vs. V0; ii) the tex2html_wrap_inline1375 and tex2html_wrap_inline1377 coefficients derived in the previous step are used to calculate the linearized ratio tex2html_wrap_inline1379; and then the fit of tex2html_wrap_inline1381 vs. tex2html_wrap_inline1383 to the straight line A V + B for each observing day provides the daily varying sensitivity A and offset B.

3.2. Crossed method

Boumier et al. (1994) discussed a method of analyzing data in the case of a resonant spectrometer sampling the sodium line profiles at 4 points, 2 in the red wing and 2 in the blue wing. A 4-point spectrometer is realized by modulating the instrumental magnetic field B by a small amount tex2html_wrap_inline1393 around its working value (Isaak & Jones 1988). In this case, it is possible to define the two crossed ratios tex2html_wrap_inline1395, and tex2html_wrap_inline1397, with the opposite signs denoting the two magnetic configurations tex2html_wrap_inline1399 and tex2html_wrap_inline1401. The calibration, tex2html_wrap_inline1403, is defined in terms of the average ratio tex2html_wrap_inline1405, and is obtained from tex2html_wrap_inline1407, where tex2html_wrap_inline1409 is the velocity displacement corresponding to the magnetic modulation, and tex2html_wrap_inline1411. Rigorously, the coefficient tex2html_wrap_inline1377 depends on the observed velocity itself, and then has to be estimated to carry on the procedure. To follow the time evolution of tex2html_wrap_inline1415, Boumier et al. fit the calibration factor with a 9th order polynomial on their one day dataset.

Looking at the expression of tex2html_wrap_inline1415, we note that its main time dependence comes from the inverse dependence on tex2html_wrap_inline1419, because it is tex2html_wrap_inline1421, and tex2html_wrap_inline1423 (see Boumier et al. 1994, their Eq. (8)). In fact, in our simulation of the IRIS offset velocities we have experimented the crossed method, following the whole time evolution of tex2html_wrap_inline1415 from August 1 to November 29, 1991, and we got the result that the active region signal was essentially canceled by this calibration. Therefore, in our opinion, one should apply to the crossed method the caution that the higher the order of the polynomial representing the calibration factor, the more precise the observed velocity fit, but also the higher the filtering probability of possible signals, as, on the other hand, Boumier et al., in the same paper, have stated for the polynomial method and the 2-point measure.

Eventually, we have selected for our simulation a calibration method, which is based on the Pallé et al. calibration procedure with one main difference: on one hand, we have at disposal only one ratio per day, because there is one solar image in the CaII K line per day, and, on the other hand, we have the possibility to build the ratio tex2html_wrap_inline1427 corresponding to the Sun without active regions; therefore, we fit directly the calibration function tex2html_wrap_inline1429 to the table tex2html_wrap_inline1427 vs. V0.

In conclusion, we emphasize that it is important to understand and model the different causes of solar noise in order to remove them directly without the risks inherent in methods based essentially on data smoothing. In other words, the best approach to calibration is to reduce as much as possible the unknown part of the velocities used for calibration.


next previous
Up: Active region effects

Copyright by the European Southern Observatory (ESO)
This email address is being protected from spambots. You need JavaScript enabled to view it.