next previous
Up: Active region effects

2. Methods

Our approach was originally motivated by the aim to simulate the possible detection of g - modes by a Na I resonance spectrometer, like the one proposed for the GOLF experiment on board of SOHO, or those in use by the IRIS network. Therefore, we planned to follow this scheme:

The interaction with the GOLF community has stimulated us to afford first of all the problem of solar noise, and, for this work, we have modified the scheme outlined above in the first two items that became

We assume that active regions consist of two components, a warmer structure or plage and a cooler structure or spot. Let us discuss how the daily distribution of spots and plages over the solar disk has been determined.

2.1. Active region data

For the period from June 10, 1990 to November 29, 1991, the daily distribution of spots and plages over the disk has been determined from the Ca  II K line full disk archive of the Big Bear Solar Observatory.

For example, in Fig. 1 (click here), we show the K line image of the Sun for February 1, 1991 from BBSO, and the distribution we have inferred for plages and spots, after removing the center-to-limb variations. To get the active region areas, we used as contrast thresholds an upper bound of 0.8 for spots, and a lower bound of 1.3 for plages. However, because the fraction of solar disk covered by plages depends strongly on the latter threshold, in Sect. 4 we report also on results obtained with a plage lower bound of 1.15.

Note that eventually we used only the images from August 1 to November 29, 1991 for comparison with IRIS, because in this period the IRIS offset velocities have the minimum scatter, and also the BBSO archive has a gap from August to December 1990 as a result of the well-known earthquake.

2.2. Atmospheric models and sodium line formation

As model atmosphere for the quiet Sun we have selected model C of Fontenla et al. (1993). The spot model was the version of the Maltby et al. (1986) umbral model modified by Severino et al. (1994). Finally, the plage model is one of the active region models developed by Andretta & Giampapa (1995). Following the approach of Cram & Mullan (1979), these models have been built by shifting inward the temperature structure, by a fixed amount tex2html_wrap_inline1261, where m is the mass column density (in tex2html_wrap_inline1265) of the quiet atmosphere. In particular we used the model corresponding to tex2html_wrap_inline1267 (hereafter called PLAGE07). Moreover, we shifted inward of the same amount also the microturbulence structure. As shown by Severino et al. (1993), we selected this particular plage model, because, in this way, it is possible to reproduce the "antiplage'' behaviour observed in the sodium tex2html_wrap_inline1213 line at about 100 mÅ from the line center, (Ulrich et al. 1993). Otherwise, when the sodium line is computed with the Fontenla et al. model atmosphere for plage, the plage appears bright at the working points of a sodium cell (e.g. Fig. 3 of Paper I). However we remark that the agreement we got between computed and observed plage contrasts is qualitative, and in Sect. 4 we will allow for a number of different contrasts (see Figs. 6 (click here) and 7 (click here)).

Finally, we used the version 2.0 of the code MULTI of Mats Carlsson (1986) for the NLTE sodium line synthesis; more details on this are given in Sect. 2 of Paper I.

2.3. Flux determination

Since the velocity measures refer to low degree modes, we have to integrate the emergent intensity over the entire disk to find the flux in the sodium line, accounting for the three-component time dependent distribution of the intensity over the disk. We performed the flux integration analytically, which allowed us to have more insight into results. However, we checked that for a number of different distributions of spots over the disk the analytical method produces velocity fluctuations which agree within 5% with the results of the direct numerical integration.

  figure267
Figure 1: Distribution of active regions over the solar disk on February 1, 1991. From left to right, and top to bottom: CaII K line image from the Big Bear Solar Observatory; the same image with the center-to-limb variations removed; inferred plage distributions obtained with a lower bound for the plage contrast of 1.30 and 1.15 respectively; inferred spot distribution used for the simulation. The total projected area covered by spots is 0.27 percent of disk, and that covered by plages is 1.5% and 5.7% respectively

The analytical formulation of the irradiance variations produced by active regions was extensively developed in Sect. 3 of Paper I. The final result of that section was that, with the assumptions that the disk is covered by a discrete distribution of active regions with a contrast depending only upon the heliocentric angle, and that the emergent intensity can be adequately represented by the linear terms in the solar rotational velocity, we can write the emergent flux at wavelength tex2html_wrap_inline1273 as
equation272

where tex2html_wrap_inline1275 is the flux emergent from the quiet disk, broadened by solar rotation and shifted by tex2html_wrap_inline1277 because of the relative velocity V0, which includes the Sun-Earth relative motions and the gravitational redshift; the other two terms represent the modulations due to the active regions, with tex2html_wrap_inline1281 independent of the solar rotation, and tex2html_wrap_inline1283 at the first order in the rotation velocity. Both these terms can be expressed by summations over the active regions effectively present on the disk as (see also Livingston et al. 1991)
equation285

equation288

displaymath1271
where Ai is the area (in units of solar hemisphere), tex2html_wrap_inline1287 is the emergent intensity of the quiet Sun, Ci denotes the contrast of the tex2html_wrap_inline1291 active region located at distance tex2html_wrap_inline1293 from the disk center and with azimuth tex2html_wrap_inline1295, and tex2html_wrap_inline1297 is the line-of-sight component of the rotational velocity.

Once the active region distribution is known from the K line daily images, and the relative velocity is determined according to the ephemerides, we can compute the fluxes in the blue and red flanks of the sodium tex2html_wrap_inline1213 line from Eqs. (2), (3) and (4), and then we can determine the photometric ratio r according to Eq. (1).


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.