Statistical comparison of the observed and synthetic hydrogen Lyman line profiles in solar prominences
S. Gunár^{1}  P. Schwartz^{1}  B. Schmieder^{2}  P. Heinzel^{1}  U. Anzer^{3}
1  Astronomical Institute, Academy of Sciences of the Czech Republic,
25165 Ondrejov, Czech Republic
2  Observatoire de Paris, Section de Meudon, 92195 Meudon Principal
Cedex, France
3  MaxPlanckInstitut für Astrophysik, KarlSchwarzschildStrasse 1,
85740 Garching, Germany
Received 6 October 2009 / Accepted 23 January 2010
Abstract
Aims. We analyse a unique set of prominence
SOHO/SUMER Lyman spectra by comparing it with synthetic spectra
obtained by 2D multithread prominence finestructure models.
Methods. We employed a novel statistical approach to
the analysis of the observed and synthetic Lyman spectra. We compared
the statistical distributions of the line properties of the observed
and synthetic Lyman spectra using a set of four statistical
criteria.
Results. We demonstrate the very good agreement
between the observed Lyman spectra and synthetic spectra obtained by
modelling.
Conclusions. Our set of statistical criteria is
wellsuited to analyses of the prominence Lyman spectra because of its
sensitivity to a number of different parameters governing the
conditions in the prominence fine structures.
Key words: Sun: filaments, prominences  radiative transfer  line: profiles  techniques: spectroscopic
1 Introduction
With progress in solar observations from space, reliable UV and EUV data on prominences have become available. This offered new possibilities for spectroscopic investigation of the prominence plasma and magnetic fields. Hydrogen Lyman lines and continuum observations, in particular, represent a unique probe to the structure and physical properties of quiescent prominences. A substantial amount of prominence spectral data in Lyman lines was provided by SUMER (Solar Ultraviolet Measurements of Emitted Radiation) UVspectrograph (Wilhelm et al. 1995) onboard SOHO (Solar and Heliospheric Observatory). The SOHO prominence observations were reviewed by Patsourakos & Vial (2002); see also a review of the prominence fine structures by Heinzel (2007).
The interpretation of the prominence observations in the Lyman lines requires sophisticated models of the prominence plasma in the magnetic field and complex nonLTE radiative transfer computations. The first nonLTE 1D model in magnetohydrostatic (MHS) equilibrium based on the KippenhahnSchlüter magnetic structure (Kippenhahn & Schlüter 1957) with multilevel radiative transfer computations was developed by Heasley & Mihalas (1976). Subsequent development of similar isothermalisobaric 1D models was led by a number of authors (Gouttebroze, Fontenla, Heinzel, Vial). The substantial influence of the prominencecorona transition region (PCTR) on the Lymanline intensities was demonstrated by Anzer & Heinzel (1999), by showing that a departure from isothermal models was needed. Thus Heinzel & Anzer (2001) constructed the 2D nonLTE model of an individual prominence finestructure thread in the MHS equilibrium with the PCTR, which exhibits different shapes along and across the magnetic field. This considerably affects the resulting Lyman line profiles. Line profiles obtained across the magnetic field (narrow PCTR with a steep temperature gradient) tend to be significantly reversed, while profiles obtained along the magnetic field (wide PCTR with gradually increasing temperature) are usually purely emission. This corresponds to the assumptions of Heinzel et al. (2001) and is further corroborated by Heinzel et al. (2005) and Schmieder et al. (2007). The importance of the multithread prominence finestructure models is shown by Gunár et al. (2007), who used a trialanderror method to reproduce the observed Lyman spectrum. Furthermore, Gunár et al. (2008) added random lineofsight (LOS) macroscopic velocities to each thread of the multithread model, thus showing that even relatively low LOS velocities (of the order of 10 km s^{1}) have a significant effect on the asymmetries of the Lyman line profiles.
In the present work we compare the observed Lyman spectra with the synthetic spectra resulting from multithread modelling of the prominence fine structures based on Gunár et al. (2008). However, instead of comparing individual Lyman line profiles we analyse statistical distributions of line properties of the observed and synthetic spectra characterized by following statistical criteria: the integrated intensities of the Lyman lines, the Lyman decrement (ratios of the integrated intensities of Lyman lines to Lyman), the redtoblue peak intensity ratio that characterizes the asymmetries of the line profiles, and finally the centre to mean peak intensity ratio characterizing the depth of the profile reversals. For statistical purposes, we use large sets of observed and synthetic spectra.
The paper is organized as follows. Section 2 gives details regarding the observed spectra. Section 3 briefly presents the multithread models and provides details on the input parameters of the models used in the present study. Section 4 shows results of the comparison of the observed and synthetic spectra for four statistical criteria, and Sect. 5 contains discussion and our conclusions.
2 Observed spectra
The observed spectra was obtained by the SUMER on May 25 and 26, 2005. The observations were made in the sitandstare mode with the spectrograph slit crossing the southern part of the prominence, as shown in Fig. 1. The observed data set consists of the whole hydrogen Lyman line series obtained using four wavelength windows taken with the detector A. The exposure time of 115 s was used for each observation. The Lyman line was obtained outside the attenuator in the first spectral window. The Lyman line was obtained in the second, the Lyman with the Lyman in the third, and higher Lyman lines up to the Lyman continuum were obtained in the fourth spectral window.
The observations were made in three time blocks with the correction for solar rotation applied between them. The first block was made on May 25, 2005 between 16:03:15 UT and 20:40:59 UT, the second block of observations was made between 20:46:27 UT and 23:01:26 UT, and the third block was made between 23:04:54 UT on May 25, 2005 and 01:21:51 UT on May 26, 2005. The number of observations in each spectral window in a given observational block is indicated in Table 1.
Figure 1: Raster obtained by SOHO/CDS (Coronal Diagnostic Spectrometer) (Harrison et al. 1995) in the He I 584.33 Å line on May 25, 2005 between 15:59:33 and 16:52:26 UT. The position of the SUMER slit is marked by the vertical bar crossing the southern leg of the prominence. Fullline part of the bar marks the working part of detector A. 

Open with DEXTER 
Three different values of the detector voltage were applied in different times of Lyman observations. However, we did not detect any significant influence of the different voltages on our statistical analyses of the observed spectra, so we took all Lyman observations into account regardless of the different voltages. An example of the observed spectra is shown in Fig. 2.
The observed spectra were reduced and calibrated using standard SolarSoft procedures for the SUMER data. The following procedures were applied in this order: decompression of binary data saved in IDLsave files, deadtime correction, flatfielding, localgain correction, and correction for geometrical distortion of the detector. Finally, the data were calibrated to units with the radiometry procedure (Schühle 2003, private communication). For detailed technical informations about the instrument, corrections, and procedures see http://www.mps.mpg.de/projects/soho/sumer/text/webluca/ch_iust.html and references therein. The lines were identified with the SUMER spectral catalogue (Curdt et al. 2001). The spectral image on the detector is inclined with respect to the detector's horizontal axis due to different orientations of the grating and the detector. Thus the spectra of different lines are vertically shifted with respect to each other. In addition there is a vertical shift due to the displacement of the slit image on the detector caused by the nonlinearity of the grating focus mechanism (Schühle 2003). The shifts caused by both effects were computed for all the lines using the SolarSoft procedure delta_pixel and these corrections were taken into account in our modelling.
For the purpose of this statistical analysis we used only the observed spectra between pixels 5 and 35 for Lyman to Lyman (Fig. 2). In this part all the lines have high intensities, so the region is certainly part of the observed prominence. We consider only the Lyman to Lyman lines because there are fewer observations of the higher members of the hydrogen Lyman line series, and they exhibit higher levels of noise.
Figure 2: Example of the observed spectra obtained on May 25, 2005 by SOHO/SUMER. Lyman line was observed at 16:22:45, Lyman at 16:20:27, and Lyman and Lyman at 16:18:18 UT. Dashed lines show the portion of the spectra we used in our analyses. 

Open with DEXTER 
Table 1: Number of observations in each spectral window within each time block.
3 Synthetic spectra
To obtain the synthetic Lyman lines, we employ the same multithread prominence finestructure models as in Gunár et al. (2008). These consist of sets of identical 2D singlethread models (developed by Heinzel & Anzer 2001) without any mutual radiative interaction. Single 2D threads are vertically infinite and are in MHS equilibrium. The empirical temperature structure is characterized by a steep gradient of the temperature across the magnetic field lines and by a gradual rise of the temperature along the magnetic field from the central minimum temperature to the PCTR temperatures at the boundaries. For the detailed description of these 2D thread models, we refer the reader to Heinzel & Anzer (2001) and Heinzel et al. (2005).
The 2D finestructure threads are arranged along individual
magnetic field lines perpendicular to the LOS. Each thread is randomly
shifted in its position along given magnetic field line and has
randomly assigned LOS velocity
(Fig. 3;
for example of one particular realization see also Fig. 2 of Gunár et al. 2008).
The emerging intensity from such multithread model at a given
wavelength and position along the length of the foremost thread
(thread N) in the observers frame is
described by a formula
where and are defined as
Here is the Doppler shift for the thread n corresponding to its LOS velocity given by
where represents the linecentre wavelength of a given spectral line, and is the LOS velocity of thread n. Positive values of represent velocities oriented towards the observer. The intensity represents the intensity of radiation emerging from a single thread n obtained by solving the nonLTE radiative transfer problem within a 2D finestructure thread illuminated by the solar surface (see Heinzel & Anzer 2001). The mutual radiative interaction of the individual threads is not taken into account.
Figure 3: Scheme of the multithread model with N randomly shifted 2D threads and randomly assigned LOS velocities (). B is the horizontal component of the magnetic field. 

Open with DEXTER 
For the purpose of this statistical study, we have also modified the formal solver to compute the emerging radiation for angles different from the directions across or along the magneticfield lines. For any given realization of the multithread model (similar to the one shown in Fig. 3 but with the LOS inclined to the direction perpendicular to the magnetic field), the positions of intersections of the LOS with the surface of each 2D thread were calculated. Then the distribution of the source function and the opacity along the inclined LOS for each thread was obtained using a linear interpolation from the source function and opacity distributions obtained by solving the nonLTE radiative transfer problem within a single 2D finestructure thread illuminated by the solar surface. The optical thickness along the inclined LOS was computed for each thread by integrating the interpolated opacity along the geometrical intersection of the LOS and the 2D thread. The emerging intensity from each 2D thread in the direction of the inclined LOS ( from Eq. (1)) was computed by formal solution of the radiative transfer using the quadrature formula of Mihalas et al. (1978). The total emerging intensity in the direction of the inclined LOS was then computed using Eq. (1). Again, no mutual radiative interaction of the threads is taken into account. Note that in the case of the inclined LOS, the actual values of the LOS velocities are those of the velocities with LOS perpendicular to the magnetic field reduced by appropriate inclination factor.
3.1 Input parameters
Table 2: List of input parameters of the 2D singlethread models.
It is not the aim of this paper to look for the prominence model with the best agreement between particular synthetic and observed Lyman line profiles. Therefore, we focus our analyses on variation of several input parameters of our multithread models centered on the model labeled hereafter as M ODEL1. This model was found to have the best agreement between observed and synthetic Lyman line profiles for the same prominence by Gunár et al. (2007) and it was used for the analysis of the Lyman line asymmetries by Gunár et al. (2008). The other two models presented in this study are minor variations of the reference M ODEL1 and are thus named V AR1 and V AR2.
The input parameters of a given multithread model fall into two distinctive categories. First category involves the physical parameters of the individual 2D fine structure thread, while the second category represents the layout of the threads within the multithread model. The input parameters of the individual 2D threads are described in Heinzel & Anzer (2001), so here we only give their brief outline (see Table 2). The parameter B_{x}(0) represents the magnetic field strength in the middle of the thread, M_{0} gives the maximum column density, and p_{0} is the boundary pressure. The temperature structure of the individual 2D thread is characterized by minimum central temperature T_{0}, maximum boundary temperature , and also by the shape of the PCTR with two different gradients of the temperature along and across the magnetic field lines. These temperature gradients can be described by exponents and , where prescribes very steep gradient across the field lines and represents gradual rise of the temperature along the field lines from the centre of the thread towards its boundaries (see Heinzel & Anzer 2001).
The layout of the threads of the multithread model is characterized by the total number of threads, by the distance between them, by their relative shifts in the direction along the magnetic field, and by the LOS velocities of the individual threads. In our present analysis we assume the total number of threads to be 10 and the distance between individual threads to be 1000 km for all models. For the maximum relative shift of individual threads with respect to the foremost thread, we assume that it is half of the length of the threads. Relative shifts of threads and their actual LOS velocities are generated randomly for each realization of given multithread structure. Thus the only input parameter describing the layout of the threads, which varies between individual multithread models, is the interval of the randomly generated LOS velocities. As was shown by Gunár et al. (2008), even relatively low LOS velocities (of the order of 10 km s^{1}) can have a considerable effect on the asymmetries of the Lyman line profiles. In the present investigation we consider three intervals of the LOS velocities, , , and km s^{1}.
3.2 Synthetic data sets
To obtain sufficiently extensive synthetic data sets of the Lyman line profiles, we computed one hundred realizations of each multithread model. Each realization has the same input parameters for individual 2D threads and the same total number of threads. We assume all threads of a given multithread model to be identical. The multithread realizations differ in the sets of the positional shifts of individual threads with respect to the foremost thread and in the actual LOS velocities of the individual threads. Then each synthetic data set consists of 83 equidistantly spaced positions along the foremost thread of the multithread model (one quasi slit spectrum) for each of one hundred realizations of the given multithread model, representing one hundred quasi observations at different times. Note that one pixel on the SUMER slit corresponds to approximately 700 km and the distance between two positions of our quasi slit is approximately 340 km.
All resulting synthetic Lyman line profiles are convolved with the instrumental profile of the SUMER A detector. This instrumental profile is a Gaussian with a standard deviation equivalent to one SUMER A detector spectral pixel, which corresponds to 0.043 Å. This agrees with the results of Chae et al. (1998), who measured the instrumental FWHM (full width at half maximum) of the SUMER A detector to be approximately 2.3 detector pixels. The FWHM for the Gaussian distribution is equal to approximately .
4 Statistical criteria
Figure 4: Integrated intensity histograms of the observed and synthetic (model V AR2, M ODEL1, and model V AR1) Lyman to Lyman lines. The xaxis gives the integrated intensities in erg cm^{2} s^{1} sr^{1}. The yaxis represents the number of counts normalized to 1000 points in the data set. The inscription in each plot gives the median and the MAD as a percentage of the median for the observed spectra and for models V AR2, M ODEL1, and V AR1, respectively. 

Open with DEXTER 
As shown by Heinzel et al. (2005) and Schmieder et al. (2007), the inclination of the LOS with respect to the orientation of the magnetic field plays an important role in the analysis of the Lyman line profiles. The line profiles obtained with the LOS perpendicular to the magnetic field exhibit significant reversals, even in the case of high members of the Lyman line series. On the other hand, the Lyman line profiles obtained with the LOS oriented along the magnetic field tend to be unreversed, especially the higher members of the Lyman series. This effect can be explained by a difference in the shape of the PCTR when observed across and along the magnetic field lines. The LOS perpendicular to the magnetic field intersects a very narrow PCTR with a steep gradient of the temperature. The LOS oriented along the magnetic field lines intersects much wider PCTR with a shallower decrease in the temperature from the prominence fine structure boundary towards its centre. Therefore, one might be able to infer the inclination of the LOS with respect to the orientation of the magnetic field from the shape of the observed Lyman line profiles.
The observed spectra we use in the present study mostly incorporates profiles with the significant reversals, which implies that the LOS during these observations was moreorless perpendicular to the predominant magnetic field inside the observed prominence. Therefore, we focus our statistical analyses on the multithread models with the LOS inclined no more than 60 from that perpendicular to the magnetic field.
The multithread models with the LOS close to the one parallel to the magnetic field show increasing predominance of purely emission Lyman line profiles, which is in a good agreement with the previous theoretical works. In the following analyses, we assume the LOS to be perpendicular to the magnetic field if not stated otherwise.
To assess the statistical properties of the observed and synthetic data, sets we use two robust statistical measures, the median and the median absolute deviation (MAD). The median is defined as a value in an ordered sample data set with an equal number of values above and below it. The MAD is defined as the median of the absolute deviations from the data set's median. We often give the MAD as a percentage of the median.
4.1 Integrated intensities
The most obvious criterion for statistical analyses of the spectral line profiles is the integrated intensity of the spectral lines. By comparing integrated intensities of the observed and synthetic spectra, we can easily distinguish the models with synthetic spectra considerably different from the observed one.
Figure 4 shows the comparison of the integrated intensity histograms of the observed (solid black lines) and synthetic Lyman to Lyman line profiles. The histograms for three different models are displayed. Green dashed lines represent model V AR2, blue dashdotted lines represent M ODEL1, and red dotted lines represent model V AR1. The xaxis gives the integrated intensities in erg cm^{2} s^{1} sr^{1}. The yaxis represents the number of counts normalized to 1000 points in the data set (Counts *1000/N, where N represents the total number of points in the data set). The inscription in each plot gives the median and the MAD as a percentage of the median for observed spectra and for models V AR2, M ODEL1, and V AR1.
Histograms of models V AR2 and M ODEL1 agree rather well with the histogram of the observed integrated intensities. In contrast, the histogram of the model V AR1, which differs from model V AR2 and M ODEL1 only in the shape of the PCTR described by parameters and , is considerably different from two other models and also from the observed spectra. We have also tested some additional models with different sets of the thread input parameters, and we have obtained even more pronounced deviations from the observed spectra than for model V AR1.
We have also analysed the effect of the choice of the LOS velocity intervals on the statistics of the integrated intensities. For this purpose we used the M ODEL1 with three different intervals of the LOS velocities, , , and km s^{1}. The effect on the integrated intensity statistics of the synthetic Lyman line profiles is negligible.
Furthermore, we studied the effect of the LOS inclination, which also proved to be relatively small for the inclinations of up to 60 from the LOS perpendicular to the magnetic field.
4.2 Lyman decrement
Figure 5: Lyman decrement histograms of the observed spectra and synthetic spectra of models V AR2, M ODEL1, and V AR1. Panels show Lyman to Lyman, Lyman to Lyman, and Lyman to Lyman ratios, respectively. The xaxis gives the ratio of integrated intensities of the given Lyman lines. The yaxis represents the number of counts normalized to 1000 points in the data set. The inscription in each plot gives the median and the MAD as a percentage of the median for the observed spectra and for models V AR2, M ODEL1, and V AR1. 

Open with DEXTER 
The Lyman to Lyman integrated intensity ratio is quite sensitive to the geometrical and thermodynamical properties of the prominence fine structures. Even more so is the whole series of the Lyman, Lyman, and Lyman to Lyman ratios, known as the Lyman decrement. Therefore the Lyman decrement represents a very useful criterion for statistical comparison of the observed and synthetic spectra.
Figure 5 shows the Lyman decrement histograms of the integrated intensity ratios of the Lyman, Lyman, and Lyman to the Lyman. Histograms of the observed and three sets of the synthetic spectra are shown. The xaxis gives the ratio of integrated intensities of the respective Lyman lines.
The histograms and their respective median values representing models V AR2 and M ODEL1 are in very good agreement with the observed spectra. Especially in case of Lyman to Lyman, and to some extent also Lyman to Lyman ratio, even the shape of the histograms represents the observed spectra rather well. In contrast, model V AR1 gives results that significantly differ from the observed data set throughout the whole Lyman decrement in both the shape of the histograms and their corresponding medians. Additional test models show even more pronounced deviations from the observed spectra than model V AR1.
The choice of the LOS velocity intervals has only a negligible effect on the statistics of the Lyman decrement. However, the inclination of the LOS with respect to that perpendicular to the magnetic field has a tangible effect on the Lyman decrement statistics, especially on the Lyman to Lyman ratio. This seems in contrast with the statistics of the integrated intensities where the effect of the LOS inclination is negligible for all Lyman lines (Sect. 4.1). However, one needs to consider that we compared the histograms of the integrated intensities in the first case, while in the second case we computed particular Lyman line ratios in each point of the spectra locally, and only then we did compare the histograms of the particular ratios. To study the effect of the LOS inclination on the Lyman decrement statistics we used the M ODEL1 with 0 (blue dashdotted lines in Fig. 5), 30, and 60 angle between the inclined LOS and that perpendicular to the magnetic field. Medians of the Lyman line integrated intensity ratios are given in Table 3. In this case the median values of the synthetic spectra with the best agreement with the observed spectra are those with the LOS perpendicular to the magnetic field. This might suggest the actual orientation of the LOS during the observations, but more profound investigation of the effects of the LOS inclination is needed.
Table 3: Median values of the Lyman line integrated intensity ratios for the M ODEL1 with different LOS inclinations.
4.3 Red to blue peak intensity ratio
Figure 6: Histograms of the ratios of the specific intensities of the red to dominant BLUE peak ( left panels) and blue to dominant RED peak ( right panels). The observed and the synthetic (model V AR2, M ODEL1, and model V AR1) Lyman and Lyman spectra are shown. The xaxis gives the ratio of the specific intensities in the given Lyman line peaks. The yaxis represents the number of counts normalized to 1000 points in the data set. The inscription in each plot gives the median, and the MAD as a percentage of the median for the observed spectra and for models V AR2, M ODEL1, and V AR1. 

Open with DEXTER 
Figure 7: Same as in Fig. 6, except that here we show the Lyman and Lyman lines. 

Open with DEXTER 
Figure 8: Same as in Fig. 6, except that here green dashed lines represent M ODEL1 with LOS velocities from the interval (marked M ODEL1 _V5), blue dashdotted lines represent M ODEL1 with LOS velocities from the interval (marked M ODEL1), and red dotted lines represent M ODEL1 with LOS velocities from the interval (marked M ODEL1 _V20). 

Open with DEXTER 
Figure 9: Same as in Fig. 6, except that here we show the Lyman and Lyman lines. Green dashed lines represent M ODEL1 with LOS velocities from the interval (marked M ODEL1 _V5), blue dashdotted lines represent M ODEL1 with LOS velocities from the interval (marked M ODEL1), and red dotted lines represent M ODEL1 with LOS velocities from the interval (marked M ODEL1 _V20). 

Open with DEXTER 
The significant asymmetries of the observed Lyman line profiles in prominences are a wellknown fact (e.g. Vial et al. 2007; Heinzel et al. 2001). Gunár et al. (2008) uses the multithread method for modelling the prominence fine structures to obtain the synthetic Lyman line profiles with asymmetries comparable to the observed ones. These authors show that even relatively low LOS velocities (of the order of 10 km s^{1}) can produce significantly asymmetrical Lyman line profiles. Such strong asymmetries cannot be attributed entirely to the effect of the Doppler shifts of the line profiles. Instead, as shown by Gunár et al. (2008), the asymmetrical profiles are produced by the combined effect of various Doppler shifts of the emission line profiles and the absorption coefficient profiles (optical thickness profiles) of the individual threads, together with relatively steep gradients of the absorption profiles and steep increase in intensity inside the line core. Gunár et al. (2008) also showed that such multithread fine structure models can produce Lyman line profiles with an asymmetry opposite to that of the higher Lyman lines at the same position. Such opposite asymmetries of different Lyman lines at the same position in the observed prominence were noticed by e.g. Vial et al. (2007).
Figures 6 and 7 show histograms of the ratios of the specific intensities in the peaks of the Lyman and Lyman, and Lyman and Lyman line profiles, respectively. The left panels of both figures display the specific intensity ratios of the red to dominant BLUE peaks, while the right panels show the ratios of the blue to dominant RED peaks. The xaxis gives the ratio of the specific intensities in the respective Lyman line profile peaks.
Figures 6 and 7 show that models V AR2, M ODEL1, and V AR1 represent the overall pattern of the observed spectra histograms reasonably well; however, they fail to reproduce the strong predominance of the red peaks over the blue peaks of the Lyman and Lyman line profiles. The histograms of the ratios of the specific intensities in the peaks are not affected by the shape of the PCTR of the individual threads of the multithread models.
On the other hand, the choice of the LOS velocity intervals has a significant effect on the statistics of the profile asymmetries. Figures 8 and 9 again show histograms of the ratios of the specific intensities in the peaks of the Lyman and Lyman, and Lyman and Lyman line profiles. The left hand panels display the specific intensity ratios of the red to dominant BLUE peaks and the right hand panels show the ratios of the blue to dominant RED peaks. Instead of V AR2 and V AR1 we show here two variations of the M ODEL1 with LOS velocities from the interval marked as M ODEL1 _V5 and from the interval marked as M ODEL1. M ODEL1 _V5 gives better agreement with the observed spectra histogram for Lyman and Lyman lines, but not as good agreement with higher Lyman lines. However, M ODEL1 _V20 differs from the observed spectra.
Table 4: Total number of profiles with the given asymmetry for the observed and synthetic spectra.
Table 4 gives the total number of profiles with the given asymmetry (dominant blue and red peak, respectively) for the observed and synthetic spectra. These numbers show the predominance of the red or blue asymmetry in the given data set. They can also indirectly indicate the occurrence of the opposite asymmetries of different Lyman lines at the same position in the prominence, when the predominance of one asymmetry for a given Lyman line changes to the predominance of the opposite asymmetry for another Lyman line. Such behaviour is in agreement with the findings of Vial et al. (2007) and Gunár et al. (2008) for observed and synthetic prominence Lyman spectra, respectively.
The effect of the inclination of the LOS with respect to the normal to the magnetic field is relatively small; however, the increasing angle between the LOS and the normal to the magnetic field produces a slightly lower number of asymmetrical profiles. Thus the increase in the angle of the LOS and the normal to the magnetic field has a complementary effect to that of a lowering of the LOS velocity interval, as one would expect.
4.4 Centre to meanpeak intensity ratio
Figure 10: Centre to mean peak intensity ratio histograms of the observed and synthetic (model V AR2, M ODEL1, and model V AR1) Lyman to Lyman line profiles. The xaxis gives the ratio of the centre to mean peak specific intensity. The yaxis represents the number of counts normalized to 1000 points in the data set. The inscription in each plot gives the median and the MAD as a percentage of the median for the observed spectra and for the models V AR2, M ODEL1, and V AR1. 

Open with DEXTER 
The last statistical criterion we employed in our statistical comparison of the observed and synthetic spectra is the depth of the reversal of the Lyman line profiles, characterized by the ratio of the specific intensity in the centre of the line to the averaged intensity of two peaks.
Figure 10 shows the histograms of the centre to the mean peak intensity ratios of the Lyman to Lyman lines. The xaxis gives the ratio of the centre to mean peak specific intensity.
For the Lyman and Lyman lines, the histograms of the observed and synthetic spectra are in relatively good agreement. However, the synthetic spectra are unable to reproduce the observed spectra histograms of the Lyman and Lyman lines. Therefore, the synthetic Lyman and Lyman line profile reversals are much deeper than the observed ones. Because the central intensities of the observed and synthetic line profiles agree rather well, the depth of the synthetic line profile reversals is due to very high peak intensities.
Unfortunately, we did not find any combination of the input parameters of the multithread models that would reproduce the observed spectra histograms. Nor has the LOS inclination any effect on the height of the peaks of the Lyman line profiles.
5 Discussion
We used four criteria wellsuited to statistical analyses of the Lyman spectra of the prominences. Using these statistical criteria, we analysed a unique, extensive data set of observed Lyman lines, including the Lyman line obtained outside the SUMER attenuator. We compared the observed spectra with synthetic Lyman spectra obtained from 2D multithread models of the prominence fine structure with randomly assigned LOS velocities (Gunár et al. 2008).
Although it was not the purpose of this paper to find a multithread model with the best match of the synthetic and observed Lyman spectra, we have demonstrated very good agreement between the observed spectra and the synthetic spectra resulting from M ODEL1. The set of input parameters of the individual 2D thread of the multithread M ODEL1 (see Table 2) was obtained by Gunár et al. (2007), who used a trailanderror method to find the model with the best agreement between the actual shape of the synthetic and observed Lyman line profiles. These authors used a part of the observed spectra also analysed in the present study. The approach of Gunár et al. (2007) of matching the shapes of individual synthetic Lyman line profiles with the observed profiles is significantly different from the analysis of the statistical criteria applied to all the data sets of synthetic and observed Lyman spectra.
The statistical criteria we have employed in this study are sensitive to the range of the input parameters of the 2D multithread models, so they represent a useful tool for assessing the properties of the prominence fine structures. The integrated intensities of the Lyman lines (Sect. 4.1) are suitable for analysis of the overall physical parameters of individual 2D finestructure threads and are especially sensitive to the shape of the PCTR (see Fig. 4) which is, in our case, described by parameters and .
The Lyman decrement criterion (Sect. 4.2) is very sensitive to all physical parameters of the individual threads and also to the inclination of the LOS with respect to the perpendicular to the magnetic field. Thus it is useful for more precise estimation of the physical parameters of the individual 2D threads within multithread models and for determining the inclination of the LOS. The median values of the Lyman decrement ratios (see Fig. 5) of the M ODEL1 are in very good agreement with observed spectra in all three studied ratios. This agreement has not been previously achieved (see e.g. Heinzel et al. 2001); however, no extensive statistical study has been done before. The Lyman to Lyman ratio has recently been considered by Vial et al. (2007), who analysed a part of the same prominence studied in the present paper. The values of the observed Lyman to Lyman ratio presented in Fig. 5 agree with values obtained by these authors and correspond to the earlier observations by the OSO8 (see Vial 1982). The synthetic Lyman and Lyman integrated intensities were studied by Heinzel et al. (1987) using 1D nonLTE models of individual prominence slabs. These authors obtained low integrated intensities of the synthetic Lyman line, which leads to large synthetic Lyman to Lyman ratios in comparison with the observations. Similar results were obtained also by Gouttebroze et al. (1993) (see Table 5). Inclusion of the PCTR to the 1D prominence modelling (Anzer & Heinzel 1999) led to higher Lyman integrated intensities and thus to lower Lyman to Lyman ratios. The Lyman to Lyman ratio was also studied by Vial et al. (1989) using 1D multislab models. Similar 1D multislab models were used by Fontenla et al. (1996), who also used the ambipolar diffusion. However, the resulting Lyman to Lyman ratio starkly contrasts with the observed values (see also conclusions of Vial et al. 2007).
Table 5: Lyman to Lyman integrated intensity ratios obtained by various prominence models.
The asymmetries of the Lyman line profiles are represented by the red to blue peak intensity ratio criterion (Sect. 4.3). This is mostly sensitive to the choice of the LOS velocity intervals (Figs. 8 and 9), so is useful for assessing of the LOS velocities of individual threads of the multithread models.
The depth of the reversed Lyman line profiles is described by the centretomeanpeak intensity ratio criterion (Sect. 4.4). The synthetic spectra obtained by multithread modelling are unable to reproduce the statistics of the observed Lyman spectra (Fig. 10). We did not succeed in finding a good fit to the observed spectra by any realistic combination of the input parameters of the multithread models. This comes from very high peak intensities compared to the observed profiles. The central intensities of the observed and synthetic profiles are in good agreement. Although the peak intensities of the synthetic line profiles are much higher than observed ones, the integrated intensities of the observed and synthetic lines are in good agreement, because of the sharpness of the synthetic line peaks and smaller linewidth of the synthetic Lyman lines, compared to the observed spectra.
6 Conclusions
Our 2D multithread models of the prominence fine structures are able to produce synthetic Lyman spectra in very good agreement with the observed spectra when evaluated by a number of statistical criteria. They are also able to reproduce the shape of the observed Lyman line profiles in particular parts of the spectra. In masse, however, our synthetic Lyman line profiles exhibit rather sharp peaks with high intensities, which are not commonly found in the observed spectra. This is, in fact, a rather common phenomenon in modelling not only of prominences but also of the quite chromosphere and other solar atmospheric features. This problem could be caused by effects of the redistribution in the radiative transfer computations, by some effects of the 2D radiative transfer that might be addressed by full 3D modelling or by the effect of the mutual radiative interaction of individual threads of the multithread model.
Our preliminary tests of the effect of mutual radiative interaction on the synthetic Lyman line profiles show a slight enhancement of the peak intensities in comparison with models without mutual radiative interaction. This even more emphasizes the discrepancy between the intensity and sharpness of line peaks of the observed and the synthetic spectra. However, these tests were only performed on ensembles of 3 and 5 finestructure threads without any LOS velocities. More realistic computations of the whole multithread models with LOS velocities and mutual radiative interaction are needed. However, such models are extremely timeconsuming, because the whole multithread model needs to be calculated within one large computational box, which is not the case for our method.
The solution of this phenomenon would represent an important step in modelling various solar atmospheric features. In achieving this, our set of statistical criteria could provide a useful test case because of its sensitivity to a number of different physical parameters.
The same method of employing a set of statistical criteria similar to our present study could also be used for analysing helium lines, as well as the spectral lines of other elements. In such a case one would need to use comprehensive atomic models of helium or any other species to obtain realistic synthetic spectra.
AcknowledgementsS.G. acknowledges support from grant 205/09/P554 of the Grant Agency of Czech Republic. S.G. and P.S. acknowledges support from grant 205/07/1100 of the Grant Agency of the Czech Republic. P.H. acknowledges the support from grant 1QS300120506 of the Grant Agency of the Academy of Sciences of the Czech Republic. This work was also supported by ESAPECS project No. 98030 and the institutional project AV0Z10030501.
References
 Anzer, U., & Heinzel, P. 1999, A&A, 349, 974 [NASA ADS] [Google Scholar]
 Chae, J., Schühle, U., & Lemaire, P. 1998, ApJ, 505, 957 [NASA ADS] [CrossRef] [Google Scholar]
 Curdt, W., Brekke, P., Feldman, U., et al. 2001, A&A, 375, 591 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fontenla, J. M., Rovira, M., Vial, J.C., & Gouttebroze, P. 1996, ApJ, 466, 496 [NASA ADS] [CrossRef] [Google Scholar]
 Gouttebroze, P., Heinzel, P., & Vial, J. C. 1993, A&AS, 99, 513 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
 Gunár, S., Heinzel, P., Schmieder, B., Schwartz, P., & Anzer, U. 2007, A&A, 472, 929 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gunár, S., Heinzel, P., Anzer, U., & Schmieder, B. 2008, A&A, 490, 307 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Harrison, R. A., Sawyer, E. C., Carter, M. K., et al. 1995, Sol. Phys., 162, 233 [NASA ADS] [CrossRef] [Google Scholar]
 Heasley, J. N., & Mihalas, D. 1976, ApJ, 205, 273 [NASA ADS] [CrossRef] [Google Scholar]
 Heinzel, P. 2007, in The Physics of Chromospheric Plasmas, ed. P. Heinzel, I. Dorotovic, & R. J. Rutten, ASP Conf. Ser., 368, 271 [Google Scholar]
 Heinzel, P., & Anzer, U. 2001, A&A, 375, 1082 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Heinzel, P., Anzer, U., & Gunár, S. 2005, A&A, 442, 331 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Heinzel, P., Gouttebroze, P., & Vial, J.C. 1987, A&A, 183, 351 [NASA ADS] [Google Scholar]
 Heinzel, P., Schmieder, B., Vial, J.C., & Kotrc, P. 2001, A&A, 370, 281 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kippenhahn, R., & Schlüter, A. 1957, ZAp, 43, 36 [NASA ADS] [Google Scholar]
 Mihalas, D., Auer, L. H., & Mihalas, B. R. 1978, ApJ, 220, 1001 [NASA ADS] [CrossRef] [Google Scholar]
 Patsourakos, S., & Vial, J.C. 2002, Sol. Phys., 208, 253 [NASA ADS] [CrossRef] [Google Scholar]
 Schmieder, B., Gunár, S., Heinzel, P., & Anzer, U. 2007, Sol. Phys., 241, 53 [NASA ADS] [CrossRef] [Google Scholar]
 Schühle, U. 2003, SUMER Data Cookbook, Published on internet at http://www.mps.mpg.de/projects/soho/sumer/text/cookbook.html [Google Scholar]
 Vial, J. C. 1982, ApJ, 253, 330 [NASA ADS] [CrossRef] [Google Scholar]
 Vial, J.C., Ebadi, H., & Ajabshirizadeh, A. 2007, Sol. Phys., 246, 327 [NASA ADS] [CrossRef] [Google Scholar]
 Vial, J. C., Rovira, M., Fontela, J., & Gouttebroze, P. 1989, Hvar Observatory Bulletin, 13, 347 [NASA ADS] [Google Scholar]
 Wilhelm, K., Curdt, W., Marsch, E., et al. 1995, Sol. Phys., 162, 189 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Table 1: Number of observations in each spectral window within each time block.
Table 2: List of input parameters of the 2D singlethread models.
Table 3: Median values of the Lyman line integrated intensity ratios for the M ODEL1 with different LOS inclinations.
Table 4: Total number of profiles with the given asymmetry for the observed and synthetic spectra.
Table 5: Lyman to Lyman integrated intensity ratios obtained by various prominence models.
All Figures
Figure 1: Raster obtained by SOHO/CDS (Coronal Diagnostic Spectrometer) (Harrison et al. 1995) in the He I 584.33 Å line on May 25, 2005 between 15:59:33 and 16:52:26 UT. The position of the SUMER slit is marked by the vertical bar crossing the southern leg of the prominence. Fullline part of the bar marks the working part of detector A. 

Open with DEXTER  
In the text 
Figure 2: Example of the observed spectra obtained on May 25, 2005 by SOHO/SUMER. Lyman line was observed at 16:22:45, Lyman at 16:20:27, and Lyman and Lyman at 16:18:18 UT. Dashed lines show the portion of the spectra we used in our analyses. 

Open with DEXTER  
In the text 
Figure 3: Scheme of the multithread model with N randomly shifted 2D threads and randomly assigned LOS velocities (). B is the horizontal component of the magnetic field. 

Open with DEXTER  
In the text 
Figure 4: Integrated intensity histograms of the observed and synthetic (model V AR2, M ODEL1, and model V AR1) Lyman to Lyman lines. The xaxis gives the integrated intensities in erg cm^{2} s^{1} sr^{1}. The yaxis represents the number of counts normalized to 1000 points in the data set. The inscription in each plot gives the median and the MAD as a percentage of the median for the observed spectra and for models V AR2, M ODEL1, and V AR1, respectively. 

Open with DEXTER  
In the text 
Figure 5: Lyman decrement histograms of the observed spectra and synthetic spectra of models V AR2, M ODEL1, and V AR1. Panels show Lyman to Lyman, Lyman to Lyman, and Lyman to Lyman ratios, respectively. The xaxis gives the ratio of integrated intensities of the given Lyman lines. The yaxis represents the number of counts normalized to 1000 points in the data set. The inscription in each plot gives the median and the MAD as a percentage of the median for the observed spectra and for models V AR2, M ODEL1, and V AR1. 

Open with DEXTER  
In the text 
Figure 6: Histograms of the ratios of the specific intensities of the red to dominant BLUE peak ( left panels) and blue to dominant RED peak ( right panels). The observed and the synthetic (model V AR2, M ODEL1, and model V AR1) Lyman and Lyman spectra are shown. The xaxis gives the ratio of the specific intensities in the given Lyman line peaks. The yaxis represents the number of counts normalized to 1000 points in the data set. The inscription in each plot gives the median, and the MAD as a percentage of the median for the observed spectra and for models V AR2, M ODEL1, and V AR1. 

Open with DEXTER  
In the text 
Figure 7: Same as in Fig. 6, except that here we show the Lyman and Lyman lines. 

Open with DEXTER  
In the text 
Figure 8: Same as in Fig. 6, except that here green dashed lines represent M ODEL1 with LOS velocities from the interval (marked M ODEL1 _V5), blue dashdotted lines represent M ODEL1 with LOS velocities from the interval (marked M ODEL1), and red dotted lines represent M ODEL1 with LOS velocities from the interval (marked M ODEL1 _V20). 

Open with DEXTER  
In the text 
Figure 9: Same as in Fig. 6, except that here we show the Lyman and Lyman lines. Green dashed lines represent M ODEL1 with LOS velocities from the interval (marked M ODEL1 _V5), blue dashdotted lines represent M ODEL1 with LOS velocities from the interval (marked M ODEL1), and red dotted lines represent M ODEL1 with LOS velocities from the interval (marked M ODEL1 _V20). 

Open with DEXTER  
In the text 
Figure 10: Centre to mean peak intensity ratio histograms of the observed and synthetic (model V AR2, M ODEL1, and model V AR1) Lyman to Lyman line profiles. The xaxis gives the ratio of the centre to mean peak specific intensity. The yaxis represents the number of counts normalized to 1000 points in the data set. The inscription in each plot gives the median and the MAD as a percentage of the median for the observed spectra and for the models V AR2, M ODEL1, and V AR1. 

Open with DEXTER  
In the text 
Copyright ESO 2010