Issue |
A&A
Volume 494, Number 3, February II 2009
|
|
---|---|---|
Page(s) | 1091 - 1106 | |
Section | The Sun | |
DOI | https://doi.org/10.1051/0004-6361:200810448 | |
Published online | 20 November 2008 |
Dynamics of small-scale magnetic fields on the Sun:
observations and numerical simulations![[*]](/icons/foot_motif.gif)
N. Bello González1 - L. Yelles Chaouche2 - O. Okunev1,3 - F. Kneer1
1 - Institut für Astrophysik, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
2 - Max-Planck-Institut für Sonnensystemforschung, Max-Planck-Str. 2, 37191 Katlenburg-Lindau, Germany
3 - Central Astronomical Observatory of the Russian Academy of Sciences, Pulkovskoye chaussee 65/1, 196140 St. Petersburg, Russia
Received 23 June 2008 / Accepted 18 September 2008
Abstract
Context. Small-scale magnetic fields play an important role in the structure and the dynamics of the solar atmosphere.
Aims. This study aims at revealing the evolution of magnetic fields, together with granular convection in the quiet Sun by means of observations with high spatial, spectral, and temporal resolution and of numerical MHD simulations.
Methods. Time sequences from quiet Sun disc centre were obtained with the upgraded ``Göttingen'' Fabry-Perot spectropolarimeter at the Vacuum Tower Telescope, Observatorio del Teide/Tenerife, in the Fe I 6173 Å line. The data were reconstructed with speckle methods. For comparison with the observations, numerical simulations of granular magnetoconvection were carried out with the MURaM code. The intensities and Stokes vectors emerging from the simulation box were degraded in wavelength, spatial co-ordinates, and noise to the quality of the observations.
Results. The noise in the observed magnetograms from the centre-of-gravity method is
G, yielding a polarimetric sensitivity of
Mx, at a cadence of 23 s with 0
33 spatial resolution in a field of view of
.
Many of the observed V profiles in network and internetwork (IN) areas exhibit strong asymmetries that indicate strong magnetoconvection. The temporal evolutions of IN structures and of a bright point (BP), as seen in broadband and line-minimum images, in Dopplergrams, and in magnetograms, are presented. The magnetic field structure in the numerical MHD simulations is even more complex than seen in the observations. Correspondingly, the emergent Stokes profiles are often very abnormal. The degradation yields a reduction of the intrinsic field strength to the ``observed'' one by a factor 4-5. The spectral resolution of the spectrometer is adequate, yet the limitation in spatial resolution and by noise filtering swamps the details seen in the non-degraded simulations. A BP was not found in the simulations, presumably because BP's are not so common and the simulated box is quiet small: only 6 Mm wide in both horizontal directions.
Conclusions. The combination of high-resolution observations with numerical simulation is a highly valuable means for studying small-scale magnetic fields on the Sun. Two-dimensional, low-noise data with spectral resolution as good as achieved here and with spatial resolution of 0
1 and better are needed to better understand this important part of the solar magnetism.
Key words: Sun: magnetic fields - techniques: high angular resolution - techniques: polarimetric - methods: numerical - methods: observational
1 Introduction
Small-scale fields on the Sun and their dynamics have become a lively field of research during the past decades. Since the discovery of ``magnetic knots'' by Beckers & Schröter (1968) and of kilo-Gauss fields in the network of the quiet Sun by Stenflo (1973), steady improvement in polarimetric sensitivity and spatial resolution has revealed that the internetwork (IN) magnetic fields play an important role in the solar magnetism. Lites et al. (1996) estimated the IN flux emergence to
Mx/day, which exceeds the estimated flux emerging in active regions and ephemeral magnetic regions of
1020 Mx/day and
1022 Mx/day, respectively, by orders of magnitude (cf. Durney et al. 1993; Petrovay & Szakály 1993).
The literature on results from ground-based observations of IN fields is rich (e.g., Domínguez Cerdeña et al. 2003; Socas-Navarro et al. 2004; Lin & Rimmele 1999; Martínez González et al. 2008; Khomenko et al. 2003; Lites et al. 1999; Keller & et al. 1994, and references therein). Recently, the high-resolution polarimetric data obtained with the Solar Optical Telescope (SOT; e.g. Suematsu et al. 2008) on board the HINODE satellite have uncovered a wealth of small-scale magnetic structures, both in active and in quiet regions, also of IN fields (e.g. Orozco Suárez et al. 2007; Lites et al. 2008; Rezaei et al. 2007; Centeno et al. 2007). While the above measurements of small magnetic flux are based on analyses using the Zeeman effect, the Hanle effect allowed important findings on weak and entangled, turbulent fields, which are hidden to Zeeman diagnostics (e.g. Trujillo Bueno et al. 2004; Bianda et al. 1999; Faurobert et al. 2001).
Table 1:
Observing parameters. The r0 is Fried parameter estimated by KAOS and
the scanning step.
On the side of numerical simulations of small-scale magnetic field dynamics, huge progress has been achieved by several groups during the last years. Magnetohydrodynamic (MHD) simulations of magnetoconvection have arrived at a high degree of realism on small scales in space and time (Stein & Nordlund 2006; Abbett 2007; Vögler et al. 2005, see also the review by Steiner 2007). The spatial scales resolved in the simulations are typically around 20 km, much smaller than the spatial resolution of 150-360 km (0
2-0
5) obtained in today's observations of magnetic fields. Thanks to their superb resolution, much can be learned from the MHD simulations by comparison with observations. E.g., Janßen et al. (2003) studied fractal dimensions of facular regions with observation and simulation. Khomenko et al. (2005a,b) used MHD simulations of granular magnetoconvection for an interpretation of observations of IN magnetic fields in visible and infrared spectral lines. De Pontieu et al. (2006) compared MHD modelling with high-resolution observations to study the rapid temporal variability of faculae. Schüssler & Vögler (2008) found horizontal magnetic fields in their simulations in agreement with observations, e.g. from HINODE (Lites et al. 2008). They concluded that a near-surface turbulent dynamo is a viable mechanism for generating dominating horizontal fields in the mid photosphere (see also Steiner et al. 2008).
We investigate small-scale magnetic field dynamics on the Sun with both high-resolution observations and MHD simulations. We observed quiet Sun and weak activity with the two-dimensional Göttingen Fabry-Perot interferometric spectropolarimeter (in short FPI) at the Vacuum Tower Telescope (VTT)/Observatorio del Teide/Tenerife (Koschinsky et al. 2001; Bendlin et al. 1992; Volkmer et al. 1995; Bendlin & Volkmer 1995). The upgrades of this instrument in 2005 (Puschmann et al. 2006) and in 2007 (Bello González & Kneer 2008) had been successful. The performance of the FPI in spatial resolution (0
25 in broadband images and 0
33 in magnetograms) and spectral resolution (
25 mÅ) is as good as that of the SOT spectrometer/polarimeter on board of HINODE. In the data taken for this work, the polarimetric sensitivity of the FPI is
(with
intensity), approximately a factor of two less than the usual HINODE SOT polarimetric data. The temporal resolution with cadences of 20-50 s, however, is much better, thanks to the superior throughput of FPIs (e.g. Kneer 2007). In addition, the computing capacities for MHD simulations have again increased considerably during the past few years. It has become feasible to follow the temporal evolution of magnetic fields in 3D simulations of magnetoconvection for many minutes (of real time on the Sun). This allows us to compare time sequences of observed, small-scale magnetic structures with theoretical expectation and to make predictions on observations with still higher spatial resolution.
This paper is organised as follows. In Sect. 2, we describe the observations and the data analysis. Section 3 presents the observational results. Section 4 deals with the MHD simulations and their degradation for comparison with the observations. Section 5 compares the observations and the simulations, and Sect. 6 concludes the paper.
2 Observations and data analysis
2.1 Observations
The observations for this study were taken on July 23 and on July 25, 2007, from disc centre of the Sun. We used the Göttingen FPI spectrometer/polarimeter (Bello González & Kneer 2008; Puschmann et al. 2006) at the VTT/Observatorio del Teide/Tenerife. The observations were supported by the Kiepenheuer Adaptive Optics System (KAOS, von der Lühe et al. 2003). On both days, the Fe I 6173 Å line, with Landé factor
,
was scanned with the spectrometer. The reduction of the area to be read out from the CCDs allowed reduction of the scanning time. The field of view (FOV) is still appreciable,
,
with a pixel size corresponding to 0
109. Table 1 gives the parameters of the observations.
The polarimeter is based on two ferroelectric liquid crystals (FLCs), a (nominal)
retarder plate followed by a (nominal)
retarder, and a modified Savart plate (for details see Bello González & Kneer 2008). This allows the measurement of the full Stokes vector with minimised seeing-induced crosstalk. By readjustment of the FLCs before the observing period, the efficiencies of the polarimeter could be improved compared to those during the first tests of the setup in 2007. They were then
,
,
and
.
The instrumental crosstalk by KAOS and the subsequent optics to the FPI spectrometer has been measured as described in Beck et al. (2005a), while the crosstalk from the telescope proper can be taken into account by means of the polarisation model by Beck et al. (2005b). The
crosstalk was removed by subtracting an appropriate fraction of
from
such that
at continuum wavelengths. Averaged over the FOV,
was then found to be less than
.
To further shorten the spectral scanning time, we used a polarimetric setup on July 23 that measures just Stokes I and V. For this, the
retarder FLC was removed, while in front of the polarimeter a
retarder was mounted with fixed orientation of its fast axis, under 45
with respect to the zero orientation of the
retarder. This, together with the Savart plate, provides alternating left and right circularly polarised components on either half of the detector. The efficiency was high,
.
For image reconstruction with speckle methods, broadband frames were recorded strictly simultaneously with the narrow-band spectropolarimetric frames. Furthermore, dark frames and flat field frames were taken with varying pointing of the telescope and with defocused telescope. Finally, images from a continuum light source (halogen lamp) were obtained to measure the wavelength dependence of the transmission of the spectropolarimeter.
In the present study, we concentrate on the Stokes I and V signals. The degree of linear polarisation from weak fields of the Sun is low for Fe I 6173 and the signals require a special analysis, as demonstrated by Lites et al. (2008). The analysis of the linear polarisation signals and the discussion of horizontal magnetic field components are certainly important for understanding the IN magnetic fields (e.g. Steiner et al. 2008; Lites et al. 2008; Abbett 2007; De Pontieu 2002; Martínez González et al. 2007b; Lites et al. 1996; Martínez González et al. 2007a; Centeno et al. 2007). Here, we restrict the analysis to circular polarisation and defer the investigation of the measured Stokes Q and U profiles to future work.
2.2 Data analysis
Dark corrections and flat fielding were applied. The broadband frames were reconstructed with the ``Göttingen'' speckle code (de Boer 1996), which uses the spectral ratio method (von der Lühe 1984) and speckle masking (Weigelt 1977). The speckle reconstruction takes into account that the image correction by KAOS varies with distance from the AO lock point. Thus, the variation in the Fried parameter across the field of view is fitted by a Gaussian centred on the lock point. Such distribution is then used to apply the proper speckle transfer function to each isoplanatic patch (cf. Bello González & Kneer 2008; Bello González et al. 2007).
Given the estimates of the broadband scenes and the instantaneous speckle images in broad band and narrow band, the speckle reconstruction of the narrow-band images is carried out as in Keller & von der Lühe (1992) and as repeatedly described in our earlier work (e.g. Krieg et al. 1999; Bello González et al. 2005). The polarimetric measurements of the modulated Stokes vector need special care, as outlined by Bello González & Kneer (2008). The entangled Stokes signals have to be demodulated, which yields
.
The Stokes components have then been corrected for instrumental crosstalk.
2.3 Filtering
The narrow-band data were noise-filtered in wavelength and in horizontal co-ordinates (x,y) with the filters shown in Fig. 1. The wavelength filter has a constant value of 1.0 up to a certain frequency, then continues with an exponential decay, and is finally smoothed with a Gaussian. The spatial filter has a transmission of 0.5 at a frequency corresponding to 0
33.
Finally, after the calculation of the physical parameters (cf. below), their maps were subjected to temporal, subsonic k-
filtering (k = horizontal wavenumber) with a limiting propagation speed of 7.5 km s-1.
![]() |
Figure 1:
Power spectrum of I profile from a pixel in the FOV (solid) and azimuthally averaged power spectrum of the continuum intensity fluctuations (dash-dotted), together with their low-pass filters for wavelength (dotted) and for intensity (dashed, radial cut through filter); arbitrary units for power spectra, filters are normalised to 1 at frequency f=0; abscissa in units of Nyquist frequency, i.e., 1/(
|
Open with DEXTER |
2.4 Physical parameters
We deal with time sequences of intensity images, Dopplergrams, and magnetograms. The intensities, as proxies for temperature fluctuations, are from (speckle reconstructed) broadband images, from images in the continuum
near the Fe I line, and from
at line minimum. The lattest were determined from parabolic fits from the I profiles at five wavelength points around the minimum. This calculation also yields the velocity
from the Doppler shift of the fitted parabola. Furthermore, the centre-of-gravity (COG) of the I profiles gives another measurement of the velocity. We calculated the wavelength positions of the COG of profiles with three limits of upper intensity
in the profiles for the integration over the profile,
such that
,
with a=0.9, 0.8, and 0.5. These will be denoted below by 90%, 80%, and 50% levels, respectively.
We also measured the magnetic fields with the COG method applied to the
and
profiles (Semel 1967; Rees & Semel 1979). According to these authors, this yields the line-of-sight (LOS) components of the magnetic fields averaged over the resolution element and over the formation height, and reduced by spurious light from scattering in Earth's atmosphere and by the telescope optics. We denote this as
where




The zero position in the wavelength axes of the V profiles presented below are shifted by +8.6 mÅ relative to the COG position of the average I profiles in the FOVs (from disc centre). This corresponds to the convective blueshift of the average I profile due to granular motion (Schröter 1957). As will be shown and discussed below, the V profiles become very asymmetric at high spatial resolution in both the observations and the simulations. Sometimes, only one V lobe is present, sometimes two positive/negative lobes, strongly shifted in wavelength (see also Sigwarth et al. 1999). We thus refrain from trying inversions of such profiles. We do not feel confident that a Milne-Eddington inversion assuming magnetic field and velocity independent of height, or a two-component inversion with one magnetic and one non-magnetic component, will lead to realistic descriptions. Recent work by Socas-Navarro et al. (2008) has cast doubts on the reliability of inversions from line pair observations to determine the field strengths. We compare the observed profiles and the physical parameters derived from them with those from simulated model atmospheres.
To obtain an idea of the formation height of the signals, we calculated response functions for velocity
and magnetic fields
in the VAL C model atmosphere (Vernazza et al. 1981). The following atomic parameters for the Fe I 6173 Å line were chosen: statistical weight of lower level g=3, oscillator strength
,
and abundance of iron
.
We employed the method by Eibe et al. (2001). There, small perturbations are introduced as step functions that are >0 and constant from bottom to a specific height and zero above that height. Differences in the emergent intensities, or COG values, from varying height of the step, yield the response functions. The perturbation amplitudes were chosen as 0.25 km s-1 in the case of velocities and 250 G for magnetic fields. Figure 2 shows the response functions
and
,
for various limits of the profile levels included in the determination of v and
with the COG method. Similar to Wunnenberg et al. (2002), we obtain a better height discrimination of the formation of the velocity signal by the linear combinations
Here, the subscripts ``90'' and ``min'' refer to the velocity determined with the limit at the 90% level and to the velocity from the line minimum shift, respectively, and the subscripts ``high'' and ``low'' refer to formation heights at high and low atmospheric layers, respectively. Such combinations with the same coefficients were also applied to the data. As seen from Fig. 2a, the low layer response function



![]() |
Figure 2:
Upper panel a): velocity response functions
|
Open with DEXTER |
We note that the temperature response function for Fe I 6173 Å in the VAL C model has its COG at a height of 270 km and an
of
300 km. The line minimum optical depth
is also located at approximately 270 km.
2.5 Noise
After filtering we obtained a standard deviation
in the continuum of the V profiles. We added Gaussian noise at each spatial and wavelength position of the proper amount and proportional to
.
With two different realisations of noise, this yielded test data
and
with the same noise amplitude in the continuum of V, at each point in the FOV. The profiles from inter-granular lanes are somewhat more sensitive to noise because they tend to be broader and shallower. However, after applying the above subsonic filter to the noise data, the granular pattern was no longer detected in the ``magnetograms''. The analysis of the noise-simulated data gave a standard deviation of the field strength of
G. At a 3
significance level and with a spatial resolution of 0
33, i.e. a radius of the resolution element
,
this translates into a magnetic sensitivity of
Mx. We note that this sensitivity at high spatial resolution is achieved within 23 s in an FOV of
.
3 Observational results
![]() |
Figure 3:
FOV of the 2007 July 23 time series. Upper left and middle: broadband and line minimum images, respectively, with intensities relative to the average continuum intensity; lower left: magnetogram of
|
Open with DEXTER |
We present in Fig. 3 maps from one of the time sequences (2007 July 23): broadband and line minimum images, magnetograms with original amplitude, as well as clipped to
G, and velocity maps in high and low layers (cf. Eqs. (2) and (3)). In the magnetograms, we show the field strengths multiplied by -1 for presentation purpose. The maps cover, in the upper left quadrant, part of an enhanced network region, as seen from the magnetograms, from the line minimum image, and from the abnormal granulation there. The lower half of the maps appears quiet, with some weak network in the right part and a very quiet internetwork in the lower left
area. The smallest structures in the magnetograms exhibit an
of 0
33, which is the spatial resolution of the spectropolarimetric data (see also Bello González & Kneer 2008).
At the positions with the co-ordinates (x,y) = (15
5, 9
3) and (20
8, 12
5), two especially bright points are present in the broadband image, which are also recognised in the line minimum image. We study the temporal evolution of the upper right one below. Apart from the bright points, the FOV is rich in small-scale magnetodynamic processes, and more examples of temporal developments of magnetic finestructures are presented in a separate contribution (Bello González et al. 2008).
In the Dopplergrams, the signal from the 5-min oscillations has been suppressed by subtracting velocity images that were boxcar-smoothed over
.
The velocity patterns from the upper and lower atmospheric layers clearly appear different. At high layers, the granular flow structures are roundish, while, in deep layers, they appear flat in the granules with narrow lanes of downflow in the inter-granular lanes. Such structure, with smooth upflow and narrow downflow channels had been found in numerical simulations by e.g. Nordlund et al. (1997). The measured amplitudes in deep layers are high, up to 2.8 km s-1. The rms velocities in the deep and high layers amount to 0.59 km s-1 and 0.47 km s-1, respectively, after subtraction of the 5-min signal.
3.1 Enhanced network and strong fields
The enhanced network region in the upper left corner of Fig. 3 possesses strong fields, we measure close to kilo-Gauss field strengths. The magnetogram, in which the fields are clipped to
G, shows a broadening of the areas with fields >50 G. Such expansion was attributed by Lites et al. (1999) to blooming, i.e. scattering of V signal from the strong field areas. We find for the present observations that blooming can only be part of the area broadening. Also, the clipped magnetogram shows finestructure at the borders and outskirts of the strong field areas. In addition, the V profiles change across the ``bloomed'' areas and can become strongly asymmetric. Four profiles are shown in Fig. 4, two of them exhibiting strong asymmetries. Such asymmetries may occur in the presence of a magnetic canopy (or magnetopause), with rapid changes in magnetic field and velocity along the LOS (cf. Steiner 2000; Sánchez Almeida et al. 1988; Grossmann-Doerth et al. 1988).
![]() |
Figure 4: Examples of V profiles from the enhanced network region in the upper left corner of Fig. 3, labelled there with the same letters as here. |
Open with DEXTER |
3.2 Temporal evolution of bright point
The temporal evolution of the bright point structure at (
and its surroundings in Fig. 3 is shown in Fig. 5. The velocities are represented as bright for downward motion and dark for upward motion. The position of the bright point is marked in the five rows by small circles corresponding to broadband intensities at 1.1
.
The dynamic process took place at the junction of three inter-granular lanes, cf. first row. The granules did not move much during the time sequence and, apparently, did not have any influence on the dynamics, e.g. via squeezing. The bright point in the broadband image attained an intensity of 1.44 of the average intensity in the FOV.
The magnetograms (second row of Fig. 5) exhibit a highly dynamic behaviour. The field structure gets pushed around, splits into parts, disperses, and, at the end of the sequence, has almost dissolved, apparently.
Right at the beginning of the sequence, a small-scale brightening in the line minimum image has started (third row of Fig. 5). Its maximum intensity occurred at approximately the same time as for the bright point in the broadband images.
The velocities from the I profiles (forth and fifth rows of Fig. 5) show downflows in the inter-granular lanes. They are, in high layers, somewhat displaced from the position of the bright point, and coincide better with it at deep layers. Also, the velocities from the I profiles are lower in the upper atmosphere than in deep layers. There they attain speeds downward of 2.4 km s-1.
![]() |
Figure 5:
Temporal evolution of the bright point structure and its surroundings at position (x,y) = (15
|
Open with DEXTER |
Figure 6 depicts the evolution of the V profiles from the position marked at image Nr. 8 by the circle in the magnetogram sequence. This was the position of the broadband bright point. The V profiles are averaged over pixel, i.e. over the spatial resolution element, and further smoothed in wavelength. The profiles are very broad.
At the beginning of the sequence, their zero crossings are redshifted by
80 mÅ which corresponds to a downflow velocity of 3.9 km s-1. At that time, the COG gives an LOS component of the field strength of
G. In the course of the evolution, the red V lobe has decreased and has essentially vanished at time position Nr. 8 (cf. first profile on the left of the bottom row of Fig. 6). At this time, the remaining lobe is strongly broadened. COG gives a small magnetic field strength of
G. Assuming that this is indeed the blue lobe, the position where it reaches zero is redshifted by approximately 120 mÅ corresponding to a downflow velocity of 5.8 km s-1. Proceeding in time, this lobe also decreased in amplitude until close to vanishing, cf. the third and second last profiles in the bottom row, and then a ``normal'', blue-asymmetric V profile, only a little shifted in wavelength apparently started to appear, meaning that the gas in the flux tube had almost come to rest. We again emphasise that the value of
is calculated by the COG method. If the magnetic structure is unresolved and varies in brightness with respect to the surrounding non-magnetic plasma then systematic effects in the determination of
are to be expected.
![]() |
Figure 6:
Temporal evolution of V profiles at the position within the bright point in Fig. 5 where the apparent field decrease occurred. Time runs from left to right, first the upper then the lower row,
|
Open with DEXTER |
Changes in the atmospheric states and to the resulting I, Q, and V profiles on time scales of 50 s were obtained in numerical simulations of magnetic flux sheets by Steiner et al. (1998). Nevertheless, the changes in the V profiles in the present study are different from those in the cited simulations. The behaviour of a bright point and of the measured I and V profiles is also different from what is described by Bellot Rubio et al. (2001) and interpreted as convective collapse and upward moving shock. In their observation, the I profiles are strongly blueshifted during the course of the event and a blueshifted V profile arises along the time. In our observations, the I profiles from the bright point are slightly redshifted, if at all, yielding downward velocities of few hundred m s-1.
Here, we tentatively suggest the scenario that we see a downflow within the magnetic feature caused by a convective collapse (Parker 1979). This has been simulated in the slender flux tube approximation by Hasan & Schüssler (1985). Their modelling shows that an adiabatic downflow heats high atmospheric layers. The flow velocities in the models are much lower than those measured from the shift in the zero crossing of the two-lobed V profiles in Fig. 6, though. Also, the bright point studied here is a transient phenomenon. The velocities measured from the I profiles are lower than those measured from the shift of the zero crossing of the V profiles, presumably because the I profiles contain much spurious light from outside the bright point area.
![]() |
Figure 7:
|
Open with DEXTER |
We speculate that the dynamic development comes from the change in the magnetic field structure in the inter-granular lanes and at their junction by subphotospheric magnetoconvection. The magnetograms in Fig. 5 suggest a reorganisation of the magnetic field with splitting, dispersion, and re-enforcement by convective collapse in an unstable situation. The collapse causes an evacuation, consequently a brightening of the broadband intensity originating in the deep interior of the flux structures due to decreased opacity. The brightening in the line minimum images is produced by the downward entrainment of chromospheric gas (Hasan & Schüssler 1985). The apparent decrease in the magnetic field signal measured with COG results from the heating and evacuation. The magnetic field structure contains little Fe I due to ionisation and strongly reduced gas density. Certainly, this scenario has to be tested by modelling.
The formation of a bright point has also been observed recently by Nagata et al. (2008) with HINODE's SOT. They too report a fast downflow prior to the rise in intensity. From Milne-Eddington inversions, they inferred an increase in the field strength from 500 G to 2 000 G during the evolution and concluded with the mechanism of convective collapse. The finally formed magnetic and intensity feature found by Nagata et al. appears to be much more persistent than the bright point discussed here.
3.3 IN magnetic fields and their temporal evolution
The IN magnetic fields have attained much attention during the past one to two decades, as outlined in Sect. 1 of the Introduction. We concentrate first on the data from July 23, 2007. In Fig. 7, we show a subfield from a very quiet IN region, the lower left corner in the FOV of Fig. 3 at time position 10 of the sequence, of size
.
From left to right we present the broadband image, the magnetogram, the line minimum image, and the velocity at deep layers. Comparing the broadband image with the magnetogram shows that the
fields
occur almost exclusively in the inter-granular lanes (see also Bello González & Kneer 2008, Fig. 8). In a few cases, a weak signal with
of 10-15 G is detected within granules (cf. Bello González et al. 2008). The line minimum image resembles the inverse granulation pattern, but only partly.
The broadband image shows dark dots within the granules, and an example is indicated by the arrow in Fig. 7. Many more are seen in the broadband image in Fig. 3. Their properties were discussed in detail on the basis of broadband time sequences by Hirzberger et al. (1993). They also are visible in numerical simulations (Vögler et al. 2005; Khomenko et al. 2005b, and Fig. 13 below from the MHD simulations). From an inspection of the temporal evolution of some dark dots we found that they harbour LOS fields of approximately 20 G during their formation (cf. Bello González et al. 2008). Sometimes they disappear while the field remains. Sometimes they develop into intergranular lanes. They exhibit moderate downflows of 0.3-1.0 km s-1.
![]() |
Figure 8:
Temporal evolution of V profiles in IN at the positions a) ( upper panel); b) ( middle panel); and c) ( lower panel) marked in Fig. 7. The numbers given in the panels are the temporal positions,
|
Open with DEXTER |
![]() |
Figure 9:
Overview of the data from July 25, 2007; images from 9
|
Open with DEXTER |
We draw attention to the V profiles in Fig. 8 with lobes that are persistent for several time steps. Tests with pure Gaussian noise of the appropriate amplitude have shown that the shapes of such noisy ``V'' profiles have amplitudes of 0.002-0.004 ,
but only for one time position; i.e., the shapes are not persistent. The changes seen in the V profiles may result both from local development of the magnetic structures and from horizontal motion of the features. We describe the temporal evolution:
- -
- the profiles from position ``a'' developed from no signal to profiles with one strongly blueshifted lobe and obtained at t=9 (second profile from left, lower row), a red lobe. At the end, the profile appeared unshifted with blue asymmetry;
- -
- the profiles in example ``b'' started with little signal and then became strongly asymmetric with only a red lobe (upper two rightmost and lower two leftmost profiles);
- -
- in example ``c'', the profile was very red-asymmetric at the beginning, and the blue lobe is not visible at times 1 and 2 (second and third profiles in the upper row). Later from t=3 on, an asymmetric three-lobed V profile developed (times 5-7) and became a two-lobed, red-asymmetric profile again (lower row, rightmost profiles).
![]() |
Figure 10:
Three examples of temporal evolution of IN structures from the 2007 July 25 data. For each example, three rows are shown: broadband images, magnetograms, and line minimum images. Only fields with
|
Open with DEXTER |
3.4 IN evolution, data from July 25, 2007
Figure 9 gives an overview of the second observed quiet area at disc centre at 9
3 after the start of the sequence. We show from left to right: broadband image, magnetogram clipped to
G, and line minimum image. The magnetograms are somewhat noisier,
G, compared to the magnetograms from the 2007 July 23 data. This was to be expected since the efficiency at measuring Stokes V profiles was lower on July 25, due to measuring the full Stokes vector. The broadband image shows two bright points near the centre of the FOV that also appear bright in the line-minimum image and are related to strong fields, up to 300 G from COG. The movies from the area show that the transient bright points move around and seem to become squeezed into elongated structures. Apart from these, the area appears to be free of network boundary structures. The white rectangles in the magnetogram mark three sub-areas. The time changes of the structures in these areas will be described in detail below.
When studying the temporal evolution, one obtains the impression that the small-scale magnetic field structures change shape more rapidly than the granular pattern and the line minimum intensity fluctuations. While the granules and the inter-granular lanes evolve and move around smoothly in movies with time steps of
s, this appears to be a long step for variations in magnetic field structure. To avoid large jumps and trace magnetic fields better along the time, the sequences from July 25, 2007, were smoothed with a subsonic filter with a limiting horizontal velocity of 4.5 km s-1.
We discuss the time variation of the structures in the areas 1-3 in Fig. 9 by means of the time sequences presented in Fig. 10. There, for each example, the broadband images, the magnetograms, and the line minimum intensities are shown.
- 1
- Upper three rows:
- -
- Magnetic field to the left of the short dash with ``a'' in the leftmost magnetogram. This field structure is located in the intergranular lane and appears like a chequerboard, or two rows from a chequerboard. Such an appearance is very common among IN fields; the magnetic field changes rapidly there, the pattern is clearly present at the first two snap shots and at later times, but re-organised and displaced; the line minimum intensity above the pattern is marginally enhanced.
- -
- Evolution in the area to the left and above the short dash labelled with ``b''. At the beginning of the sequence, a flux concentration within the junction of intergranular lanes s being moved upwards in the magnetograms. It is related with enhanced line-minimum intensity, and in its way, it interacts with other field structures of opposite polarity and seemingly ``eats'' them. At the same time, towards the upper right of this structure and near the upper border of this sub-FOV, magnetic field inhomogeneities appear together with brightening at line minimum. The two features then move together, producing a further brightening at line minimum.
- -
- Evolution to the left of the dash labelled with ``c''. Here, a positive polarity field is present close to negative polarity with substantial brightening in line minimum. After few time steps, both fields have dispersed.
- 2
- Middle three rows, negative magnetic polarity to the right of the dash labelled with `d' in the first magnetogram:
- -
- This example shows, with time steps
s, an intensification of the magnetic field along the time and then a field dilution. The amplitude of the measured field strength here is -150 G at the time of maximum field strength at time position 3. Inspection of the full granular movie shows proper motion of the granules converging towards the centre of the sub-FOV, i.e. towards the magnetic field structure, during the first phases and, during the final phases, inward motion from two opposite sides and outward flows into the two perpendicular directions. The motion advects fields of both magnetic polarities, and the line-minimum images exhibit intensity enhancement during the phases of maximum magnetic field signal.
- 3
- Lower three rows:
- -
- Magnetic fields with negative polarity to the left of the dash labelled with ``e''. The temporal evolution of this field structure is fast. One notices advection of magnetic flux with positive and negative polarity. The structure becomes distorted, coagulates with other field strands and splits again; meanwhile, the line minimum intensities are not remarkably bright.
- -
- Intergranular magnetic fields in the half ``circle'' above the label ``f'' and indicated with dashes. This example shows especially well the alignment of the fields in intergranular lanes and their development when following in time backward and forward. Again, the evolution of the line minimum intensities does not show especially bright features.
4 Numerical MHD simulations
4.1 MHD simulations
The simulations were carried out using the compressible 3D radiation MHD code MURaM (Vögler et al. 2005; Vögler 2003). This code solves the MHD equations taking energy exchange via radiative transfer into account as well as the effect of partial ionisation in the equation of state.
The horizontal size of the simulation domain is Mm,
spanned by
gridpoints. The height of the domain is
1.68 Mm, which is spanned by 120 gridpoints. The resulting
gridspacing in the horizontal and vertical directions is 20.83 and
14 km, respectively. Periodic boundary conditions are imposed at
the lateral boundaries of the simulation domain. The simulation box
extends between -920 km below to 750 km above the average
optical depth unity (
)
at 500 nm.
The magnetic field above the upper and lower boundaries is constrained to be vertical. The upper and lower boundaries are ``open'' to plasma motions. The choice of a tall simulation box (750 km above optical depth unity) minimises the influence of the upper boundary conditions on the regions where the considered iron line is formed.
Initially, the simulations begin with a plane-parallel solar-model atmosphere and no magnetic field. After the purely hydrodynamical convection has reached a statistically steady state, a mixed polarity vertical magnetic field was introduced. The initial configuration of the magnetic field was made such that the simulation domain was split into 4 parts, and the adjacent parts have opposite polarity fields with 200 G amplitude.
The flux is redistributed by convective dynamics, and the opposite
polarity fields cancel each other out resulting in an exponential
decrease in the overall magnetic flux. Since we are interested in
studying the weak internetwork magnetic field, we choose a time
sequence of snapshots starting with a mean unsigned field strength
of 25 G at
.
The selected time sequence of snapshots
covers a time span of 9.2 mn with a cadence of 23 s between
successive snapshots.
To study the observational properties of these simulated
small-scale flux concentrations, we calculate the Stokes profiles
for the iron
line Fe I 6173.3 Å with the same atomic parameters as given
above in Sect. 2.4. The computations are done in LTE using
the STOPRO code in the SPINOR package (Frutiger et al. 2000; Solanki 1987).
The emergent Stokes profiles and the neighbouring continuum
intensity are calculated at disc centre for each of the
pixels. The thus obtained Stokes profiles are
normalised to the local continuum.
![]() |
Figure 11:
Left: vertical magnetic field component at 150 km with bright indicating direction outward of the atmosphere and dark direction inward. Right: sections of vertical cuts through the atmosphere to the left of the lines indicated by ``1'' and ``2'' in the left panel; upper panels: vertical magnetic field components with bright indicating direction outward and dark inward; lower panels: vertical velocities with bright upward and dark inward motion. The sections extend from from -50 km to +750 km, with average
|
Open with DEXTER |
![]() |
Figure 12: Stokes V profiles from MHD simulations; thin: original profiles, thick: after degradation in wavelength (see Sect. 4.2). |
Open with DEXTER |
We obtained from the emergent Stokes I and V profiles a magnetogram with the COG method, i.e. from the wavelength separation of the COGs of
and
.
Dopplergrams at ``high'' and ``deep'' photospheric layers are calculated as well, according to the combinations in Eqs. (2) and (3).
Figure 11 shows the vertical magnetic field component at a height of 150 km in the left panel, where the maximum of the magnetic field response function is located. The right panels of this figure give vertical cuts extending from 50 km below average
to the top of the simulation box. The vertical magnetic field and the vertical velocities are presented.
Figure 12 depicts some Stokes V profiles. These exhibit exotic shapes, and such profiles are very common. One can generally state that the ``strangeness'' of the V profiles increases with decreasing amplitudes. Low-amplitude V profiles often possess four to five lobes. The structures with the largest magnetic field in the simulation box have two-lobed V profiles with amplitudes
,
which are often very asymmetric. We also note without showing that the zero crossings of the strong V signals exhibit rapid changes across the magnetic features fluctuations of velocities of
km s-1. On the one hand, V profiles with large amplitudes, thus rather simple shapes, can be subjected to inversion calculations from which trustworthy, approximate atmospheric parameters can be retrieved. On the other hand, V profiles exhibiting complex shapes like those observed (e.g. Fig. 4, two profiles on the right, and Fig. 8) and as those from MHD simulations would be difficult to interpret by inversions assuming a few atmospheric parameters, e.g. by a Milne-Eddington inversion. In view of the magnetic field structure in Fig. 11 and of the profiles in Fig. 12, the hypothesis of MIcro-Structured Atmosphere (MISMA, Domínguez Cerdeña et al. 2006, and references therein) appears a viable concept for retrieving the magnetic and velocity structure, yet with sufficiently resolved observations.
4.2 Degradation to resolution of observations
For comparison with the observations, the observable intensities and Stokes profiles that emerge from the simulation box had to be degraded to the resolution of the observations in both wavelength and spatial resolution. The degradation consisted in several steps. The first degradation occurs intrinsically, by the transfer of the Stokes vector through the photosphere. This is seen when comparing the spatial structure of the vertical magnetic field at a height of 150 km in Fig. 11 with the magnetogram in the upper right corner of Fig. 13. The small-scale alternation of magnetic polarity has mostly disappeared because many different atmospheric layers contribute with variable polarity fields and with their vertical velocities of opposite flow direction (cf. vertical cuts in Fig. 11). The Stokes V profiles from various positions in the (simulated) FOV in Fig. 12 are correspondingly very abnormal.
Next we scaled the wavelength resolution and step width to the spectrometer's resolution and step width. Obviously, the COG magnetic field measurement is very robust against spectral degradation since the COG method itself is an averaging method. Broad V profiles with large amplitude from strong field structures are little affected by the limitation in spectral resolution. Some degraded V profiles are shown with the original profiles in Fig. 12. The wavelength degradation preserves the shape of the profiles well, but decreases them in amplitude.
![]() |
Figure 13:
Effects of degradation. Upper row from original simulation data, lower row after full degradation to quality of observations, from left to right: intensity in continuum near Fe I line, intensity at local line minimum,
|
Open with DEXTER |
The main effect comes from the limited spatial resolution. To account for the latter we convolved the simulated data in spatial direction with a normalised point spread function that contains the telescopic Airy function and a 5% contribution from a Lorentzian with
on the Sun. This Lorentzian with its extended wings is chosen to account for scattered light by telescope and Earth's atmosphere from regions with few arcsec distance. Additionally a ``wide-angle'' scattering has been applied with an amplitude of
,
where
is the intensity at this specific wavelength
averaged over the FOV. Finally, the images were filtered with the spatial low-pass filter to scales
0
33 (dashed filter in Fig. 1). The image resolution of 20.83 km was then binned to the pixel size of the detectors in the observations corresponding to 79 km. The spatial degradation brings the simulated rms continuum fluctuations from 0.17
to the observed fluctuations of 0.08-0.09
.
At the same time, the average unsigned LOS component decreases from 16-17 G at 150 km to 5.1-5.2 G from COG measurement.
The last degradation deals with the noise. We added Gaussian noise proportional to the intensities in the (I+V) and (I-V) profiles and filtered again in space and wavelength with the filters depicted in Fig. 1. The noise amplitude was chosen to obtain noise in the magnetic field measurements of 2 G.
Some of the effects are shown in Fig. 13 for the degradation of the intensity images and magnetograms. This figure stems from the first temporal position of the simulations with calculated profiles. The bright small-scale structure at the centre of the FOV of the continuum images has a maximum intensity of
,
so not higher than the brightest granular intensities. Its intensity in the image degraded to the resolution of the observed, speckle reconstructed broadband images (0
25) is only
above the dark background, so it cannot be considered as a typical bright point.
Figure 14, from the 11th temporal position in the simulations, shows examples of V profiles from simulations (thin curves) and from their full degradation (thick). The amplitudes of the degraded profiles are multiplied with a factor of 2, for the presentation. We comment without showing that the large-amplitude, original profiles from the strong fields are generally reproduced well, apart from some broadening. The main reason for this broadening is the rapid fluctuation of field strength and velocity across the magnetic structure. Noise filtering has a minor effect. Similarly for low-amplitude V profiles, a slow variation of the atmospheric structure including magnetic fields yields slowly varying V profiles across the resolution element of the degraded data.
The V profile ``a'' in Fig. 14 is redshifted and two-lobed. The degraded profile almost has a factor of 10 smaller amplitude than the original one and is less redshifted. The degraded profiles ``b'' and ``c'' allow the strong asymmetry to be seen. Profile ``d'' has only a tiny red lobe that is buried in the noise after degradation. The original V profile ``e'' is reproduced well in the degraded profile, apart from the amplitude and the broadening. These last two profiles ``d'' and ``e'' stem from regions with variation that it is not too fast across the resolution element. In contrast to them, the complex profiles ``f'' and ``g'' are from areas with rapid variation. Thus the degraded profiles, albeit not consisting of pure noise, exhibit no resemblance with the original ones but are apparently a mixture of very different profiles.
![]() |
Figure 14: Comparison of V profiles from original MHD simulations (thin) and from the degraded data set thick, latter multiplied with a factor of 2. The profiles stem from the eleventh temporal position of the simulated time sequence. |
Open with DEXTER |
Test calculations have shown that one needs low-noise data with a spatial resolution of 0
1 or better to preserve the Stokes profiles, at least the shapes of the V profiles, while it will be extremely difficult to preserve their amplitudes.
Table 2: Values of various parameters from observations and from simulations at different states of degradation.
In Table 2 we summarise the effects of the various degrees of degradation in different parameters, i.e.
the line-of-sight (LOS) components of unsigned magnetic field
averaged over field of view, minimum and maximum of
,
area coverage of fields
(Cov), and rms values of LOS velocities at deep and high layers. Index lim in the third column indicates limitation to
G in the original COG ``measurements'' and elimination of those pixels with larger field strengths. The first two rows in the table refer to the simulation box, and the other rows to Stokes profiles.
One notices a decrease in
by a factor 4-5, from the intrinsic field to the detected field. The effect of limited spatial resolution is the strongest. We have given two values for the average unsigned magnetic field strengths on the following grounds. The area contains some parts with kilo-G fields. Although their number of pixels in the FOV is small, a few tens, as seen from the histograms in Fig. 17, they contribute to the average unsigned field. To eliminate their effect, we chose to include in the averaging only areas with
smaller than 700 G in the COG measurements of the non-degraded data. Also, since the observations contain noise of
G, the fields of those pixels with
were set to zero in the averaging, in both the simulated and the observed data sets.
The entries for the observations from this work stem from an IN area with little to no strong magnetic flux. The average LOS magnetic field strengths of 3 G is from the observing day 2007 June 23. The initial conditions for numerical MHD simulations were actually chosen to obtain
G after the full degradation process.
The velocities are included in Table 2 to show the overshooting granular flows. The 5-min oscillations do not contribute since they are removed from the observed data by subtraction of velocities smoothed in spatial co-ordinates. The velocity amplitudes determined from the formulae in Eqs. (2) and (3) can be high, in deep layers up to
km s-1, and in high layers in the range of
km s-1. The degradation in wavelength decreases the rms velocity at deep layers and increases it at high layers due to increased mixing of contributions from different heights of formation.
![]() |
Figure 15:
Temporal evolution in numerically simulated data set. Three sequences from the fully degraded data are shown, the first two from the first half of the series, the third from its end. The (sub)columns show from left to right: continuum intensity, line minimum intensity, and
|
Open with DEXTER |
4.3 Temporal evolution in simulations
The degraded data sets, including the V profiles at each wavelength position, were subjected to a subsonic filter with a cutoff of 7.5 km s-1, as the observations. Figure 15 gives three examples of the temporal evolution of the (degraded) numerically simulated continuum intensities, line minimum intensities, and COG magnetograms in Fe I 6173 Å. Time runs from top to bottom, and the initial times are given in the upper left corners of the continuum images. The second set of columns are continuations of the first set.
The coarse distribution of the strong magnetic fields stays the same during the full sequence, and yet the finestructures, especially those of the weaker fields, change from one temporal position to the next. This is seen in the accompanying movie (see footnote). It illustrates the need for high temporal resolution. A cadence of 46 s gives some jitter of the small-scale fields from one time step to the next, which was seen also in the observations.
Generally, the magnetic fields detectable in the magnetograms from the degraded MHD simulations are located in the intergranular structures seen in the continuum images. The original data also show weak fields in the granules (cf. Fig. 13). Some apparently ``diffuse'' fields, such as in the upper left corner at the initial time of the sequence, are related to small convection elements of sub-arcsec size.
The strong positive polarity field close to the centre of the FOV has its counterpart as brightening in the line minimum images. We mentioned above that it does not have the properties of a bright point discussed in the observational part, together with Fig. 5. It is more like a line gap feature with reduced line depression. In fact, within this structure, the Stokes I profiles exhibit clear Zeeman splitting.
What is missing in the magnetograms are the many structures in the intergranular lanes with opposite magnetic polarity next to each other. They were seen in the observations (cf. Figs. 7 and 10) with field strengths from COG of 10-50 G, well above the noise. More recent numerical MHD simulations with the MURaM code, using a very high spatial resolution and correspondingly high magnetic Reynold's number, exhibits dynamo action and has many more of these features with close opposite polarity (Danilovic et al. 2008; Schüssler & Vögler 2008).
Figure 16 depicts two examples of V profile evolution. Here also the degraded profiles were multiplied by a factor of 2. The examples are selected outside the strong flux regions in the FOV. In the upper sequence of this figure, the V profiles develop, starting with two lobes and rapidly varying amplitude and asymmetry. During the first nine time positions, the degraded profiles follow the original ones closely, apart from strong amplitude reduction of the former. The positions of the zero crossings, which indicate downflows with up to 1.5 km s-1, are reproduced well. Nevertheless, towards the end of the sequence, starting from time position 10 (second row, third panel), the original profiles weaken and become very abnormal. The degraded profiles cannot follow this evolution. They represent some averages over the resolution element plus noise.
![]() |
Figure 16:
Two examples of temporal evolution of V profiles from MHD simulations, with 23 s cadence and at fixed positions in the FOV. Time runs from left to right, first upper rows then lower rows. Thin profiles: from original data; thick profiles: from fully degraded data, multiplied with a factor of 2 and averaged over
|
Open with DEXTER |
The second example in the lower part of Fig. 16 shows rapid development of the original one-lobed V profile and then a reversal of magnetic polarity within 23 s, between time positions 4 and 5. At these positions, the negative V lobe changes from being on the blue side to the red side of the zero wavelength position. The degraded profiles can reproduce the polarity change, although with noise and much broader than the original profiles. In view of the abnormal, nearly single-lobed V profiles, it is better to state that one sees a change in dominating polarity along the LOS. During the second half of the sequence (lower row), the V profiles weaken and increase again, reappearing first with a one-lobed profile, then with two lobes. As in the first example above, the degraded V profiles do not reproduce the very abnormal, low-amplitude original profiles, but find the correct polarity, at least for original profiles with sufficient amplitude.
We note from inspection of the movies that the temporal changes in the profiles result from both local evolution and lateral convective flux transport, as expected.
5 Discussion
In this section we first compare the observational results with those from the numerical MHD simulations. Then we compare our results on IN fields with other work.
5.1 Comparison of observations and simulations
A comparison of the observations with the numerical simulations has two aspects. First, one can discuss general similarities and dissimilarities. Second, one has to degrade the simulation outcome to the quality of the observations. This allows to see the processes that are behind the observations, but hidden by the lack of resolution. Thus the simulations are extremely valuable for the interpretation.
Without attempting completeness, we list several points in this comparison:
- 1.
- As known from earlier studies of MHD processes in the granulation, the magnetic fields are mainly located in the intergranular lanes as observed. They are distributed by the convection processes. The reaction of the magnetic fields is fast. Therefore, it turns out that time steps of
40 s and longer are too long for studying processes on scales of 0
3 and smaller.
- 2.
- The simulation results show that the MHD processes are much more complex than observations suggest, even at the best spatial resolution. The often very abnormal V profiles result from rapid changes in magnetic fields, in their polarity, and in velocities of several km s-1 along the LOS.
- 3.
- The non-degraded, as well as the degraded, simulated ``measurements'' of the magnetic fields do not exhibit chains of positive and negative polarity fields next to each other in intergranular lanes. Possibly a much higher magnetic Reynolds number is needed for a more adequate concept for IN fields showing more alternating polarities (Danilovic et al. 2008; Schüssler & Vögler 2008).
- 4.
- Related to this last point is that, after we had degraded the simulation output to the best of our knowledge, the interganular lanes show less field coverage in the simulations, 17% according to Table 2, than in the observations (30% coverage).
- 5.
- Some of the areas in the simulation surface appear more like network with kilo-G fields even after degradation as described above, although the transition from network to internetwork fields is not sharp. In the observational results, we have excluded such areas from the study of IN fields. Possibly, this is again an outcome of the simulation initialisation, and a numerical simulation with a surface dynamo action starting from a low-amplitude seed field could yield less network features.
- 6.
- In the observations, bright points were found with high broadband intensities well above the brightest granular level. They represent transient processes with fast downflows on small scales. We did not find any such features in the simulation - possibly because the bright points are relatively rare and the simulation box was small (6 Mm by 6 Mm in the horizontal directions) and of limited duration (9.2 mn). A longer time series or a bigger box are obvious targets for future simulations.




![]() |
Figure 17:
Histograms of
|
Open with DEXTER |
Certainly, fields exist with strengths above 250 G, but very few compared to those that entered the histograms of Fig. 17. The interesting part here is the one with low apparent field strengths, below 250 G. Both,
and
exhibit very similar distributions at fields larger than 100 G. These fields and COG ``measurements'' are from areas with one polarity strongly dominating any mixed polarity situation along the LOS. Below
20 G, both distributions are also similar. These data are from granules, with many pixels exhibiting very weak fields. Between 20 G and 100 G, the original ``measurements'' show less fields than at 150 km, and they are intrinsically degraded by radiative transfer where mixed polarities are playing a role.
An important result in Fig. 17 is the difference between the histograms of the observed COG fields and the fully degraded COG fields. In the range of 10...100 G, the latter histogram is lower than the one from the observations by a factor 2 to 3. This explains the lower area coverage of the degraded fields compared to the observed fields in Table 2 and also the lower appearance of fields in the
10...100 G range in Figs. 13 and 15 from the simulations compared with Figs. 7 and 10 from the observations.
5.2 Comparison with other work on IN fields
In Table 3, we collect values from the literature and from this work for the unsigned LOS component of the internetwork magnetic field averaged over the FOV.
Table 3:
Line-of-sight components of unsigned IN magnetic field
,
averaged over field of view, at disc centre.
Khomenko et al. (2003) have used full Stokes measurements in the Fe I lines at 1.56
to determine the field strengths B, inclination angles
and filling factors f (in their notation). With their
in approximately 50% of the FOV and B=250-300 G, one obtains
G. The authors themselves give
G averaged over approximately 50% of the IN area with a significant signal. Domínguez Cerdeña et al. (2003) performed speckle reconstruction of polarimetric images in the Fe 1 6302 Å line pair and conclude with many strong fields from the line ratio. With ASP measurements in the same line pair, Lites & Socas-Navarro (2004) performed a normalised integration over the Stokes V profile,
and calibrated this quantity. The value given by Orozco Suárez et al. (2007) is from inversion of SOT-HINODE measurements, again in the 6302 Å line pair. IN measurements from other SOT-
data and the calibration with fields in strong magnetic flux regions lead Lites et al. (2008) to an unsigned average
G. Martínez González et al. (2008) simultaneously measured the full Stokes parameters in the IR and the 6302 Å line pairs and also carried out inversions simultaneously in all four lines. In the work of Bello González & Kneer (2008), the COG method was also applied to speckle reconstructed data to measure the LOS components of the magnetic field. From the fields with
,
the data from July 23, 2007, give an average IN field
of 3 G, those from July 25, 2007, give 4.0 G. By means of the Hanle effect, Trujillo Bueno et al. (2004) found unresolved, mixed polarity (``hidden'') fields, which are possibly entangled on micro-scales.
The average unsigned LOS field component found in the present investigation from observation compares satisfactorily with values found by Khomenko et al. (2003), Orozco Suárez et al. (2007), the value in Lites et al. (2008), and the one in Martínez González et al. (2008). As noted by Lites & Socas-Navarro (2004), the field strength obtained in our observations and also in the Orozco Suárez et al. (2007) measurements is not larger than from other data despite the high spatial spatial resolution and the high sensitivity.
The entries in the lines ``simulated'' in Table 3 refer to the MHD simulations. The values of
from Table 2 are repeated here. After degradation in spectral and spatial resolution and after adding an appropriate amount of noise, the simulations were designed to give approximately 3 G for the average LOS field component measured with COG. At an atmospheric height of 150 km, from where we receive the magnetic signal with the COG method (cf. Fig. 2b), we obtain
G. After degradation in wavelength and spatial resolution, this value has decreased by a factor of 3. Then, by adding noise and analysing only the V profiles with amplitude 3
above noise, the average unsigned field strength has again decreased by a factor of 1.7.
6 Conclusions
The observations presented in the first part of this study revealed highly dynamic processes in the solar photosphere. The fast observational cadences allowed study of important new phenomena and dynamics. Some examples are given in this study, and more can be found in Bello González et al. (2008). The efforts to combine 2D spectropolarimetry of high sensitivity with efficient data gathering and image reconstruction were successful. At present, the results from the observations and data analysis are among the best to our knowledge. The data will be subjected to further analyses for more in-depth studies of the small-scale magnetic field dynamics of the quiet solar photosphere. We demonstrated that especially the investigation of temporal evolution in photospheric magnetoconvection has become possible.
Likewise, the numerical MHD simulations used here contain many more details than could be presented in this study. They contain a wealth of information on MHD processes in a stratified hot plasma to be uncovered in future investigations.
The combination of high spatial resolution observations and MHD simulations turned out to be very fruitful. We have learned much from the comparison of observational results and simulations. Much we have seen in the brightness, magnetic field, and velocity evolutions is also reflected in the simulations. And not all observational results are reproduced by the simulations, such as bright points and changes in polarity along inter-granular lanes. Possible reasons for these discrepancies include the relatively low magnetic Reynolds number in the simulations and the smallness of the simulation box that leads to limited statistics.
Most importantly, the simulations have taught that the actual properties of the magnetic photosphere is much more complex at small scales than the observations, although of high quality, suggest. The many abnormal Stokes V profiles turn out to be the normal case in simulations. The observation of such profiles is best performed by comparing with the simulations and the physical conditions producing such profiles.
In view of the results from the numerical MHD simulations and in comparison with the observations, there is still a long way to go for observational improvement. The spectral resolution of the ``Göttingen'' FPI spectropolarimeter appears to be adequate. As the degradation of the simulated V profiles have shown, the spectrometer by itself can even reproduce the shape of very abnormal V profiles with four to five lobes, although with reduced amplitude. Sacrificing the spectral purity for the sake of decreasing the noise would hide much of the dynamics reflected in the spectral profiles.
The combination of high-resolution observations with simulations of much higher spatial resolution has demonstrated the need for high spatial resolution of 0
1 and better, for high spectral resolution, for high temporal resolution, and for very low noise observation. Present detector technology will improve again and allow fast cadences of images, which are to be analysed with improved computer facilities. Perhaps the 1.5 m GREGOR telescope, which will be operative in the near future, will allow 0
1 resolution in magnetograms. This will certainly be a new breakthrough in the research into the small-scale magnetism of the Sun. Still, it is important to have many photons to reduce noise. This is a typical case for solar telescopes of the 4 m class like EST and ATST.
Acknowledgements
We are greatly indebted to Robert Cameron for his invaluable and unselfish support with the numerical MHD simulations and for his encouragement with our English. NBG acknowledges financial support by the Deutsche Forchungsgemeinschaft through grant KN 152/31-1. L.Y.C. is grateful to the Max-Planck-Institut für Sonnensystemforschung, Katlenburg-Lindau, for a stipend of the International Max Planck Research School on ``Physical Processes in the Solar System and Beyond''. OO thanks for support by Deutsche Forschungsgemeinschaft through grant KN 152/29-3 and by the Institut für Astrophysik of the Georg-August-Universität Göttingen. The Vacuum Tower Telescope is operated by the Kiepenheuer-Institut für Sonnenphysik, Freiburg, at the Spanish Observatorio del Teide of the Instituto de Astrofísica de Canarias.
References
- Abbett, W. P. 2007, ApJ, 665, 1469 [NASA ADS] [CrossRef]
- Beckers, J. M., & Schröter, E. H. 1968, Sol. Phys., 4, 142 [NASA ADS] [CrossRef] (In the text)
- Beck, C., Schmidt, W., Kentischer, T., & Elmore, D. 2005a, A&A, 437, 1159 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Beck, C., Schlichenmaier, R., Collados, M., Bellot Rubio, L., & Kentischer, T. 2005b, A&A, 443, 1047 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Bello González, N., & Kneer, F. 2008, A&A, 480, 265 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Bello González, N., Okunev, O. V., Domínguez Cerdeña, I., Kneer, F., & Puschmann, K. G. 2005, A&A, 434, 317 [NASA ADS] [CrossRef] [EDP Sciences]
- Bello González, N., Kneer, F., & Puschmann, K. G. 2007, in Modern Solar Facilities - Advanced Solar Science, ed. F. Kneer, K. G. Puschmann, & A. D. Wittmann, Universitätsverlag Göttingen, 217
- Bello González, N., Okunev, O., & Kneer, F. 2008, A&A, 490, L23 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Bellot Rubio, L. R., Rodríguez Hidalgo, I., Collados, M., Khomenko, E., & Ruiz Cobo, B. 2001, ApJ, 560, 1010 [NASA ADS] [CrossRef] (In the text)
- Bendlin, C., & Volkmer, R. 1995, A&AS, 112, 371 [NASA ADS]
- Bendlin, C., Volkmer, R., & Kneer, F. 1992, A&A, 257, 817 [NASA ADS]
- Bianda, M., Stenflo, J. O., & Solanki, S. K. 1999, A&A, 350, 1060
- Carlsson, M., Stein, R. F., Nordlund, Å., & Scharmer, G. B. 2004, ApJ, 610, L137 [NASA ADS] [CrossRef]
- Centeno, R., Socas-Navarro, H., Lites, B., et al. 2007, ApJ, 666, L137 [NASA ADS] [CrossRef]
- Danilovic, S., and the MHD team at Max-Planck-Institut für Sonnensystemforschung, Katlenburg-Lindau, Germany 2008, private communication
- de Boer, C. R. 1996, A&AS, 120, 195 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- De Pontieu, B. 2002, ApJ, 569, 474 [NASA ADS] [CrossRef]
- De Pontieu, B., Carlsson, M., Stein, R., et al. 2006, ApJ, 646, 1405 [NASA ADS] [CrossRef] (In the text)
- Domínguez Cerdeña, I., Sánchez Almeida, J., & Kneer, F. 2003, A&A, 407, 741 [NASA ADS] [CrossRef] [EDP Sciences]
- Domínguez Cerdeña, I., Sánchez Almeida, J., & Kneer, F. 2006, ApJ, 646, 1421 [NASA ADS] [CrossRef] (In the text)
- Durney, B. R., De Young, D. S., & Roxburgh, I. W. 1993, Sol. Phys., 145, 207 [NASA ADS] [CrossRef]
- Eibe, M. T., Mein, P., Roudier, Th., & Faurobert, M. 2001, A&A, 371, 1128 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Faurobert, M., Arnaud, J., Vigneau, J., & Frisch, H. 2001, A&A, 378, 627 [NASA ADS] [CrossRef] [EDP Sciences]
- Frutiger, C., Solanki, S. K., Fligge, M., & Bruls, J. H. M. J. 2000, A&A, 358, 1109 [NASA ADS]
- Grossmann-Doerth, U., Schüssler, M., & Solanki, S. K. 1988, A&A, 206, L37 [NASA ADS]
- Hasan, S. S., & Schüssler, M. 1985, A&A, 151, 69 (In the text)
- Hirzberger, J., Bonet, J. A., Vázquez, M., & Hanslmeier, A. 1999, ApJ, 527, 405 [CrossRef] (In the text)
- Janßen, K., Vögler, A., & Kneer, F. 2003, A&A, 409, 1127 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Keller, C. U., & von der Lühe, O. 1992, A&A, 261, 321 [NASA ADS] (In the text)
- Keller, C. U., Deubner, F.-L., Egger, U., Fleck, B., & Povel, H. P. 1994, A&A, 286, 626 [NASA ADS]
- Khomenko, E. V., Collados, M., Solanki, S. K., Lagg, A., & Trujillo Bueno, J. 2003, A&A, 408, 1115 [NASA ADS] [CrossRef] [EDP Sciences]
- Khomenko, E. V., Martínez González, M. J., Collados, M., et al. 2005a, A&A, 436, L27 [NASA ADS] [CrossRef] [EDP Sciences]
- Khomenko, E. V., Shelyag, S., Solanki, S. K., & Vögler, A. 2005b, A&A, 442, 1059 [NASA ADS] [CrossRef] [EDP Sciences]
- Kneer, F. 2007, in Modern Solar Facilities - Advanced Solar Science, ed. F. Kneer, K. G. Puschmann, & A. D. Wittmann, Universitätsverlag Göttingen, 67 (In the text)
- Koschinsky, M., Kneer, F., & Hirzberger, J. 2001, A&A, 365, 588 [NASA ADS] [CrossRef] [EDP Sciences]
- Krieg, J., Wunnenberg, M., Kneer, F., Koschinsky, M., & Ritter, C. 1999, A&A, 343, 983
- Lin, H., & Rimmele, T. 1999, ApJ, 514, 448 [CrossRef]
- Lites, B. W., & Socas-Navarro, H. 2004, ApJ, 613, 600 [NASA ADS] [CrossRef] (In the text)
- Lites, B. W., Leka, K. D., Skumanich, A., Martínez Pillet, V., & Shimizu, T. 1996, ApJ, 460, 1019 [NASA ADS] [CrossRef] (In the text)
- Lites, B. W., Rutten, R. J., & Berger, T. E. 1999, ApJ, 517, 1013 [CrossRef]
- Lites, B. W., Kubo, M., Socas-Navarro, H., et al. 2008, ApJ, 672, 1237 [NASA ADS] [CrossRef]
- Martínez González, M. J., Collados, M., Ruiz Cobo, B., & Solanki, S. K. 2007a, A&A, 469, L39 [NASA ADS] [CrossRef] [EDP Sciences]
- Martínez González, M. J., Asensio Ramos, A., López Ariste, A., & Manso Sainz, R. 2007b, A&A, 479, 229 [NASA ADS]
- Martínez González, M. J., Collados, M., Ruiz Cobo, B., & Beck, C. 2008, A&A, 477, 953 [NASA ADS] [CrossRef] [EDP Sciences]
- Martínez Pillet, V., Lites, B. W., & Skumanich, A. 1997, ApJ, 474, 810 [CrossRef]
- Nagata, S., Tsuneta, S., Suematsu, Y., et al. 2008, ApJ, 677, L145 [NASA ADS] [CrossRef] (In the text)
- Nordlund, Å., Spruit, H. C., Ludwig, H.-G., & Trampedach, R. 1997, A&A, 328, 229 (In the text)
- Orozco Suárez, D., Bellot Rubio, L. R., del Toro Iniesta, J. C., et al. 2007, ApJ, 670, L61 [NASA ADS] [CrossRef]
- Parker, E. N. 1979, Cosmical Magnetic Fields (Oxford: Clarendon) (In the text)
- Petrovay, K., & Szakály, G. 1993, A&A, 274, 543 [NASA ADS]
- Puschmann, K. G., Kneer, F., Seelemann, T., & Wittmann, A. D. 2006, A&A, 445, 337 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Rees, D. E., & Semel, M. D. 1979, A&A, 74, 1 [NASA ADS]
- Rezaei, R., Steiner, O., Wedemeyer-Böhm, S., Schlichenmaier, R., et al. 2007, A&A, 476, L33 [NASA ADS] [CrossRef] [EDP Sciences]
- Sánchez Almeida, J., Collados, M., & del Toro Iniesta, J. C. 1988, A&A, 201, L37 [NASA ADS]
- Sánchez Almeida, J., Landi Degl'Innocenti, E., Martínez Pillet, V., & Lites, B. W. 1996, ApJ, 466, 537 [NASA ADS] [CrossRef]
- Schröter, E. H. 1957, ZAp, 41, 141 [NASA ADS] (In the text)
- Schüssler, M., & Vögler, A. 2008, A&A, 481, L5 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Semel, M. D. 1967, Ann. Astrophys., 30, 513 [NASA ADS]
- Sigwarth, M., Balasubramaniam, K. S., Knölker, M., & Schmidt, W. 1999, A&A, 349, 941 (In the text)
- Socas-Navarro, H., Martínez Pillet, V., & Lites, B. W. 2004, ApJ, 611, 1139 [NASA ADS] [CrossRef]
- Socas-Navarro, H., Borrero, J. M., Asensio Ramos, A., et al. 2008, ApJ, 674, 596 [NASA ADS] [CrossRef] (In the text)
- Solanki, S. K. 1987, Ph.D. Thesis, ETH Zürich
- Stein, R. F., & Nordlund, Å 2006, ApJ, 642, 1246 [NASA ADS] [CrossRef]
- Steiner, O. 2000, Sol. Phys., 196, 245 [NASA ADS] [CrossRef]
- Steiner, O. 2007, in Modern Solar Facilities - Advanced Solar Science, ed. F. Kneer, K. G. Puschmann, & A. D. Wittmann, Universitätsverlag Göttingen, 321 (In the text)
- Steiner, O., Grossmann-Doerth, U., Knölker, M., & Schüssler, M. 1998, ApJ, 495, 468 [NASA ADS] [CrossRef] (In the text)
- Steiner, O., Rezaei, R., Schaffenberger, W., & Wedemeyer-Böhm, S. 2008, ApJ, 680, L85 [NASA ADS] [CrossRef] (In the text)
- Stenflo, J. O. 1973, Sol. Phys., 32, 41 [NASA ADS] [CrossRef] (In the text)
- Suematsu, Y., Tsuneta, S., Ichimoto, K., et al. 2008, Sol. Phys., 249, 197 [NASA ADS] [CrossRef] (In the text)
- Trujillo Bueno, J., Shchukina, N., & Asensio Ramos, A. 2004, Nature, 430, 326 [NASA ADS] [CrossRef]
- Vernazza, J. E., Avrett, E. H., Loeser, R. 1981, ApJS, 45, 635 [NASA ADS] [CrossRef] (In the text)
- Vögler, A. 2003, Ph.D. Thesis, Göttingen university
- Vögler, A., Shelyag, S., Schüssler, M., et al. 2005, A&A, 429, 335 [NASA ADS] [CrossRef] [EDP Sciences]
- Volkmer, R., Kneer, F., & Bendlin, C. 1995, A&A, 304, L1 [NASA ADS]
- von der Lühe, O. 1984, J. Opt. Soc. Am. A1, 510 [NASA ADS] [CrossRef] (In the text)
- von der Lühe, O., Soltau, D., Berkefeld, T., & Schelenz, T. 2003, SPIE, 4853, 187 [NASA ADS] (In the text)
- Weigelt, G. P. 1977, Optics Comm., 21, 55 [CrossRef] (In the text)
- Wunnenberg, M., Kneer, F., & Hirzberger, J. 2002, A&A, 395, L51 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
Online Material
Footnotes
- ... simulations
- A movie is only available in electronic form at http://www.aanda.org
- ...Å
- The full time sequence, containing the original and degraded simulated continuum intensities, line minimum intensities, magnetograms, and Dopplergrams, is available as a gif movie as online material.
All Tables
Table 1:
Observing parameters. The r0 is Fried parameter estimated by KAOS and
the scanning step.
Table 2: Values of various parameters from observations and from simulations at different states of degradation.
Table 3:
Line-of-sight components of unsigned IN magnetic field
,
averaged over field of view, at disc centre.
All Figures
![]() |
Figure 1:
Power spectrum of I profile from a pixel in the FOV (solid) and azimuthally averaged power spectrum of the continuum intensity fluctuations (dash-dotted), together with their low-pass filters for wavelength (dotted) and for intensity (dashed, radial cut through filter); arbitrary units for power spectra, filters are normalised to 1 at frequency f=0; abscissa in units of Nyquist frequency, i.e., 1/(
|
Open with DEXTER | |
In the text |
![]() |
Figure 2:
Upper panel a): velocity response functions
|
Open with DEXTER | |
In the text |
![]() |
Figure 3:
FOV of the 2007 July 23 time series. Upper left and middle: broadband and line minimum images, respectively, with intensities relative to the average continuum intensity; lower left: magnetogram of
|
Open with DEXTER | |
In the text |
![]() |
Figure 4: Examples of V profiles from the enhanced network region in the upper left corner of Fig. 3, labelled there with the same letters as here. |
Open with DEXTER | |
In the text |
![]() |
Figure 5:
Temporal evolution of the bright point structure and its surroundings at position (x,y) = (15
|
Open with DEXTER | |
In the text |
![]() |
Figure 6:
Temporal evolution of V profiles at the position within the bright point in Fig. 5 where the apparent field decrease occurred. Time runs from left to right, first the upper then the lower row,
|
Open with DEXTER | |
In the text |
![]() |
Figure 7:
|
Open with DEXTER | |
In the text |
![]() |
Figure 8:
Temporal evolution of V profiles in IN at the positions a) ( upper panel); b) ( middle panel); and c) ( lower panel) marked in Fig. 7. The numbers given in the panels are the temporal positions,
|
Open with DEXTER | |
In the text |
![]() |
Figure 9:
Overview of the data from July 25, 2007; images from 9
|
Open with DEXTER | |
In the text |
![]() |
Figure 10:
Three examples of temporal evolution of IN structures from the 2007 July 25 data. For each example, three rows are shown: broadband images, magnetograms, and line minimum images. Only fields with
|
Open with DEXTER | |
In the text |
![]() |
Figure 11:
Left: vertical magnetic field component at 150 km with bright indicating direction outward of the atmosphere and dark direction inward. Right: sections of vertical cuts through the atmosphere to the left of the lines indicated by ``1'' and ``2'' in the left panel; upper panels: vertical magnetic field components with bright indicating direction outward and dark inward; lower panels: vertical velocities with bright upward and dark inward motion. The sections extend from from -50 km to +750 km, with average
|
Open with DEXTER | |
In the text |
![]() |
Figure 12: Stokes V profiles from MHD simulations; thin: original profiles, thick: after degradation in wavelength (see Sect. 4.2). |
Open with DEXTER | |
In the text |
![]() |
Figure 13:
Effects of degradation. Upper row from original simulation data, lower row after full degradation to quality of observations, from left to right: intensity in continuum near Fe I line, intensity at local line minimum,
|
Open with DEXTER | |
In the text |
![]() |
Figure 14: Comparison of V profiles from original MHD simulations (thin) and from the degraded data set thick, latter multiplied with a factor of 2. The profiles stem from the eleventh temporal position of the simulated time sequence. |
Open with DEXTER | |
In the text |
![]() |
Figure 15:
Temporal evolution in numerically simulated data set. Three sequences from the fully degraded data are shown, the first two from the first half of the series, the third from its end. The (sub)columns show from left to right: continuum intensity, line minimum intensity, and
|
Open with DEXTER | |
In the text |
![]() |
Figure 16:
Two examples of temporal evolution of V profiles from MHD simulations, with 23 s cadence and at fixed positions in the FOV. Time runs from left to right, first upper rows then lower rows. Thin profiles: from original data; thick profiles: from fully degraded data, multiplied with a factor of 2 and averaged over
|
Open with DEXTER | |
In the text |
![]() |
Figure 17:
Histograms of
|
Open with DEXTER | |
In the text |
Copyright ESO 2009
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.