Issue |
A&A
Volume 513, April 2010
|
|
---|---|---|
Article Number | L6 | |
Number of page(s) | 4 | |
Section | Letters | |
DOI | https://doi.org/10.1051/0004-6361/200913584 | |
Published online | 20 April 2010 |
LETTER TO THE EDITOR
Evidence of small-scale magnetic concentrations dragged by vortex motion of solar photospheric plasma
L. Balmaceda1,2 - S. Vargas Domínguez3 - J. Palacios1 - I. Cabello1 - V. Domingo1
1 - Image Processing Laboratory, University of Valencia, PO Box 22085, 46980 Paterna, Valencia, Spain
2 - Instituto de Ciencias Astronómicas, de la Tierra y el Espacio, ICATE-CONICET, San Juan, Argentina
3 - Mullard Space Science Laboratory, University College London, Holmbury St Mary, Dorking, Surrey, RH5 6NT, UK
Received 2 November 2009 / Accepted 18 March 2010
Abstract
Vortex-type motions have been measured by tracking bright points in
high-resolution observations of the solar photosphere. These
small-scale motions are thought to be determinant in the evolution of
magnetic footpoints and their interaction with plasma and therefore
likely to play a role in heating the upper solar atmosphere by twisting
magnetic flux tubes. We report the observation of magnetic
concentrations being dragged towards the center of a convective vortex
motion in the solar photosphere from high-resolution ground-based and
space-borne data. We describe this event by analyzing a series of
images at different solar atmospheric layers. By computing
horizontal proper motions, we detect a vortex whose center appears to
be the draining point for the magnetic concentrations detected in
magnetograms and well-correlated with the locations of bright points
seen in G-band and CN images.
Key words: convection - Sun: granulation - Sun: photosphere - magnetic fields
1 Introduction
Convection is an important mechanism of energy exchange responsible for
the formation of the extensively studied solar granulation. The
interaction between the convective plasma flows and solar magnetic
fields causes the appearance of solar structures on many spatial
scales, sunspots and pores being the most conspicuous ones. Magnetism
is often understood to be the main process controlling the whole Sun's
configuration and behavior, even in the so-called quiet regions where
small magnetic elements such as bright points (hereafter BPs) are
detected mainly along dark lanes in-between granules. Small-scale (0.5 Mm) convectively driven vortex-type motions were discovered by Bonet et al. (2008). These motions were suggested by Sturrock & Uchida (1981)
to account for about half of the photospheric kinetic energy, but this
could not be confirmed due to the limitations of early observations. On
larger spatial scales, i.e., granular sizes, vortical motions were
observed by Brandt et al. (1988), Tarbell et al. (1991), and Title et al. (1992). Signatures of photospheric vortex flows have also been observed on even larger scales such as supergranular junctions (Attie et al. 2009). Vortical motions in the solar photosphere were studied by simulations (Zirker 1993,
and references herein) and are understood to play an important role in
the evolution of magnetic footpoints and heating in the upper solar
atmosphere (Brandt et al. 1988; Wang et al. 1995).
Here we describe a vortex motion in the solar photosphere that persists
for about 20 min. It seems to affect the small-scale magnetic
concentrations, which appear to be dragged towards the vortex center.
![]() |
Figure 1:
Quasi-simultaneous observations showing the evolution of a quiet sun region on 29 September 2007. From top to bottom:
Mg I magnetogram, G-band, CN, Ca II H. Every panel
displays the average image over 4-min intervals and the time stamps
correspond to the initial time for each interval, respectively.
BPs are detected in all filtergrams and are co-spatial with the
magnetic concentrations observed in magnetograms. Spatial scales in
both axis are in arcsec and correspond to solar coordinates (FOV |
Open with DEXTER |
2 Observations and data processing
An international campaign at the Canary Islands observatories was
performed during September-October 2007 and corresponds to an
important collaboration among several European and Japanese
institutions. Coordinated observations involving all the ground-based
facilities at Canary Islands and the space solar telescope Hinode (Kosugi et al. 2007) were performed as part of the Hinode Operation Program 14.
In the present work, we use data recorded during a particular
observing run on September 29, 2007, whose target was a quiet
sun region north of the active region NOAA 10971, close to the
solar disk center ( = 0.99).
2.1 Ground-based observations
We used the Swedish 1-m Solar Telescope (SST, Scharmer et al. 2003) to acquire images in G-band (
nm) with an effective field-of-view (hereafter FOV) of 68
54
68
54 and a sampling of 0.034''/px. Image post-processing included the Multi-Object Multi-Frame Blind-Deconvolution (MOMFBD, Van Noort et al. 2005)
restoration technique to correct for atmospheric and instrumental
aberrations. The final products were two time series of images (s1: 08:47-09:07 UT and s2:
09:14-09:46 UT) with a cadence of 15 s. Additional data
post-processing steps were: compensation for diurnal field rotation,
rigid alignment of the images, correction for distortion, and subsonic
filtering to eliminate residual jittering (Title et al. 1986).
2.2 Satellite data
The Solar Optical Telescope (SOT, Tsuneta et al. 2008) onboard Hinode acquired filtegrams in CN (
FWHM 0.8 nm) and the core of Ca II H line (
0.3 nm) using the Broadband Filter Imager (BFI) with a cadence of
35 s. The observed FOV covered 19
18
74
09. Although Hinode's FOV is smaller than that of SST, they still overlap. Magnetograms in the Mg I line (
nm) were also obtained at a cadence of
20 s with the Narrowband Filter Imager (NFI). Data from the spectropolarimeter (SP, Ichimoto et al. 2008) were also acquired. This dataset consisted of Stokes parameters I, Q, U, and V measured along a slit of 256 pixels in raster scan mode of 18 scans from 08:20-09:44, with a sampling of 0.15
/px and FOV of
.
The operation mode was set to dynamic mode,
i.e., the exposure time per slit position was 1.6 s to study
the evolution of highly dynamic events. The noise level was 1.6
in Stokes I and 1.8
Stokes Q and U. The observed lines were the two FeI lines,
nm and
nm.
The SOT images were corrected for dark current, flat-field, and
cosmic rays with standard procedures. A subsonic filter was
applied to remove high frequency oscillations in filtergrams.
SP data were inverted using the full atmosphere inversion code
LILIA (Socas Navarro 2001), based on the SIR code (Ruiz Cobo & Del Toro Iniesta 1992).
![]() |
Figure 2:
From top to bottom: magnetic field strength and LOS velocity single maps obtained by inverting SP data. The FOV covered is shown in the left lower panel. Positive/negative
|
Open with DEXTER |
3 Data analysis and results
3.1 Multiwavelength description of the event
Quasi-simultaneous data taken by Hinode and SST were used. The maximum time-lag was 2 s and
12 s
between CN and Ca images and between filtergrams and
MgI magnetograms, respectively. Images were rotated, co-aligned,
and trimmed to maintain the same FOV at all wavelengths. Figure 1
shows from bottom to top, time-averaged Ca II H, CN, and
G-band filtergrams and Mg I magnetograms from 08:40 to
09:20 UT. Each frame corresponds to an average over 4 min,
obtained from
7 images in the case of filtergrams and
12 MgI magnetograms.
The averaging of the consecutive images allows us to easily identify
the structures present in successive images. We superimpose on the
corresponding filtergrams magnetogram contours, which are useful for
visualizing the co-spatiality of the features observed at different
wavelengths. The grayscale bars indicate the intensities normalized to
the mean intensity of the entire data set for the respective
wavelengths. The LOS magnetograms shown here consist of masks
obtained after applying a threshold of 3
to remove the noise following Krivova & Solanki (2004).
The good agreement between the magnetic features detected in these
images and the bright elements present in the filtergrams gives us
confidence to justify its use. On first inspection, the magnetic
feature appears to undergo an apparent rotation over the 40 min
elapsed. We can identify a magnetic structure with two main lobes
(denoted by L1 and L2 in the first frame, top row of
Fig. 1) one of which rotates around the other as indicated by the plus black and white plus-signs marking the 4-min averaged location of the centroid for L1 and L2, respectively
.
This movement is also detected at the various wavelengths,
i.e., different heights in the solar atmosphere. Both G-band and
CN data provide information about the photospheric level, while
Ca II H images and Mg I magnetograms retrieve
information from the low chromosphere. During the interval
between 09:00 and 09:08, the L2 lobe seems to
dissolve. The magnetic flux density in this region appears to decrease,
but increases again when the whole structure has rotated and L2
approaches the lower part of the FOV. An increase in intensity is
observed in Ca images and is accompanied by the appearance of BPs
in G-band, while in CN the signal is more diffuse.
As a consequence of the rotation, small-scale processes of
fragmentation and coalescence of BPs clearly discernible in G-band
filtergrams take place along the intergranular lane. The BPs are
also detectable in CN images, although their appearance is not so
well defined as in the SST data, especially towards the end of the
sequence. The good agreement between the small elements detected
in Ca II H and the other wavelengths is also remarkable.
A forthcoming work (Vargas Domínguez et al. 2010; hereafter
Paper I) will focus on the evolution of BPs in the sequence.
![]() |
Figure 3:
Map of horizontal velocities ( FWHM 1
|
Open with DEXTER |
3.2 Characterization of the magnetic configuration of the region
Figure 2 displays the evolution of the magnetic structure in the sequence of single maps of magnetic field strength |B| retrieved by LILIA (first row) and LOS-velocity,
,
obtained by estimating the Doppler shift in the line-center-wavelength using the center-of-gravity method to the Stokes I profiles (second row). For the inversion, only pixels with Stokes Q, U, or V amplitudes
larger than 3 times their noise levels are taken into
consideration to exclude profiles that cannot be inverted reliably.
As observed in MgI magnetograms, the two magnetic lobes
can be clearly distinguished in the |B| maps displayed in
the first row. SOT-SP data start earlier than the filtergrams and
magnetograms described in the previous section. The first 5 frames
display the evolution of the magnetic structure prior to 08:40,
where we can see that it is already rotating. The maximum value of the
magnetic field strength was close to 600 G at 08:42, though
not included in the displayed images. The rms of the errors in |B|
is 194 G. We defined the errors to be the differences between the
inferred from the model and the real parameters at optical depth
log
= -1.
To check the reliability of the inversion, we also compared these
values with those obtained by assuming a weak-field approximation (Landi degl'Innocenti 1992)
and achieved similar results. The longitudinal magnetic flux decreases
by about 30% during the whole time interval. The area of the
structure seems to shrink with time. In particular,
for 09:00-09:20 the magnetic patches evolve from being bound
together, when the two structures lose their identity and become
elongated. The Y-shape structure that is evident in the
MgI magnetogram (at about 09:12 UT) is not visible here.
One possible explanation is that Fe and MgI lines provide
information about different heights in the solar atmosphere. The
missing component of ``Y'' exhibits in MgI magnetograms a
very weak and variable signal that is co-spatial with bright regions in
CaII images. However, the component is absent in photospheric
filtergrams and SP maps. SP maps are available
until 09:44 (though not shown in the figure) and show that the
structure continues rotating after 09:20 when the SOT-NFI and
BFI time series end, L1 and L2 approaching each other
towards the end of the observing interval. The inclination
(not shown here) remains stable during the whole observing period.
Magnetic flux tubes are vertically oriented with inclination values
close to 180 degrees. However, changes in the inclination can be
seen at 08:44, when some BPs appear in-between the magnetic lobes.
In spite of the high noise level in Stokes Q and U,
the good agreement between the magnetic structure observed in the
LOS magnetogram and that of SP maps is indicative of nearly
vertical magnetic fields. In the
maps (second row of Fig. 2), the granules are clearly distinguished as upflows with velocities from 1 to 2 km s-1. The intergranular lanes are the areas showing downflows with velocities from -1.3 to -2 km s-1. A strong downflow region appears at 08:44 UT (encircled area in Fig. 2) that denotes the location where the photospheric vortex forms (Sect. 3.3).
It remains there for several minutes (until 08:56) and
at 09:20 another downflow region is observed at the location
of L2 at the end of the sequence, though much less intense. The
magnetic patches are located in the lanes, with a couple of exceptions
where the magnetic areas do excursions into granules (at 09:00
and 09:07 UT). During the whole sequence, numerous
small-scale processes take place leading to the formation of BPs
(Paper I).
3.3 Photospheric plasma flows
The G-band series from SST were used to analyze the horizontal
proper motions of structures in the FOV. Proper motions were
measured by using the local correlation tracking technique (LCT, November & Simon 1988) implemented by Molowny-Horas & Yi (1994). Maps of horizontal velocities are calculated for the time series by using a Gaussian tracking window of FWHM
,
i.e., roughly half of the typical granular size. Figure 3 shows the flow map computed from the first of the SST time series (s1). This map of horizontal velocities is calculated by averaging over the total duration of the series (
min). The background image represents in false-color the normalized (factor of 0.086 s-1) divergence field. Horizontal velocity magnitudes are averaging-dependent values (Vargas Domínguez 2008),
though the flow patterns yielded by the analysis of proper motions are
commonly assumed to represent plasma motions thus enabling the
detection of general trends, which is our main aim here
(e.g., sinks or granulation downdrafts where the cold plasma
returns to the solar interior). The map of horizontal displacements is
dominated by flows coming from granular explosive events and commonly
associated with mesogranulation (Bonet et al. 2005; Roudier & Muller 2004).
We identified only a pair of examples where strong sinks are present at
intergranular lanes (coordinates [6.5, 3.5]
and [8.5, 7.5] in Fig. 3)
but, in contrast to other areas displaying downflows, these are
characterized by being the draining point where all horizontal velocity
vectors in the neighborhood converge. This behavior was observed in
numerical simulations by Nordlund (1985), who showed that the plasma on granular scales in strong downflow regions tend to rotate around the center, like a bathtub
effect. Both vortex-type events last for several minutes
(>20 min). Nevertheless, the one framed in the two boxes
in Fig. 3 is not only less stable and has a defined onset in our data (at
08:48 UT,
Paper I) but also appears to be the final destination of the
magnetic concentrations being dragged from an upper location
(at coordinates [5, 7] in the same figure),
as described in Sect. 3.1. The latter vortex remains rather stable but disappears abruptly in the map of horizontal velocities computed for s2, which is dominated by an organized flow coming from the lower right part of the FOV (Fig. 4, right). Unfortunately, as already mentioned, there is a 7-min gap between s1 and s2 meaning that it is not possible to define exactly when and how the vortex disappears.
![]() |
Figure 4: Temporal evolution of the magnetic centroid for both lobes (L1 and L2) in Fig. 1. The trajectories followed by the centroids are plotted independently and superimposed on the flow maps computed for SST/G-band images in time intervals: s1 ( left) and s2 ( right). The background false-color image is the divergence field computed from each respective flow map. Alphabet letters correspond to UT times as follows: A-08:40, B-08:48, C-08:56, D-09:00, E-09:04, F-09:08, G-09:12, H-09:21, I-09:24, J-09:28, K-09:32, L-09:36, M-09:40, N-09:44. The FOV is the same boxed in solid line in Fig. 3. |
Open with DEXTER |
4 Discussion
In Fig. 4, we plot the trajectories followed by the centroid of the distribution of magnetic flux density for lobes L1 and L2. The centroids
denoted with plus-signs and asterisks for the respective lobes are now
calculated for the individual maps obtained from SP data with a
cadence of 36 s using the magnetic flux density M for weighting. They are defined by
=
and
=
,
where (xi, yi) are the coordinates of the ith pixel
for each lobe. In spite of its lower resolution, we use the
SP data set because its time series extends until 09:44 UT,
allowing us to track the dynamics of L1 and L2 over a longer
period. The trajectories are superimposed on both the respective flow
maps computed for the SST time series (s1 and s2)
and the divergence field for each series. The comparison with the
centroid positions calculated using the MgI magnetograms (top row
in Fig. 1) allowed us to estimate the errors to be
,
for L1 and
,
for L2. From the centroid positions, it was possible to
estimate the horizontal velocities for each lobe. While L1 motion
is confined to a radius of approximately 1
,
L2 travels a much larger distance. In the time interval prior
to 09:14, its motion seems to be strongly influenced by the
vortex formed in the coordinates [3.5, 1.5] in the left panel
of Fig. 4, increasing its velocity from
1.5 to 3 km s-1
as it approaches this location. After the vortex has disappeared,
L2 continues its movement with decreasing velocity probably
influenced by the continuous evolution of the surrounding granules
towards a new sink denoted by the region of negative divergence whose
center is located at coordinates [1.0, 2.5] in Fig. 4 (right). Because of the lack of space here, additional analysis of centroids velocities will be deferred to Paper I.
A possible explanation of the observed magnetic concentrations being dragged towards the vortex center could be that this is the result of the stochastic evolution of granules that allow the BPs to approach the vortex influence, increase their velocities, and eventually fall into the vortex (surface scenario). On the other hand, long-lasting (>1 h) vertical magnetic flux tubes anchored beneath the surface might be directly influenced by the vortex action due to some interaction mechanism taking place underneath. Although presenting additional evidence to support this scenario is beyond the scope of this paper, this possibility should not be ruled out.
AcknowledgementsThe Swedish 1-m Solar Telescope is operated by the Institute of Solar Physics of the Royal Swedish Academy of Sciences at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias. Hinode is a Japanese mission developed and launched by ISAS/JAXA, with NAOJ as domestic partner and NASA and STFC (UK) as international partners. It is operated by these agencies in co-operation with ESA and NSC (Norway). The authors would like to thank the anonymous referee for helpful comments on the manuscript.
References
- Attie, R., Innes, D. E., & Potts, H. E. 2009, A&A, 493, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bonet, J. A., Márquez, I., Muller, R., Sobotka, M., & Roudier, Th. 2005, A&A, 430, 1089 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bonet, J. A., Márquez, I., Sánchez Almeida, J., Cabello, I., & Domingo, V. 2008, ApJ, 687, L131 [NASA ADS] [CrossRef] [Google Scholar]
- Brandt, P. N., Scharmer, G. B., Ferguson, S., Shine, R. A., & Tarbell, T. D. 1988, Nature, 335, 238 [NASA ADS] [CrossRef] [Google Scholar]
- Ichimoto, K., Lites, B., Elmore, D., et al. 2008, Sol. Phys., 249, 233 [NASA ADS] [CrossRef] [Google Scholar]
- Kosugi, T., Matsukaki, K., Sakao, T., et al. 2007, Sol. Phys., 243, 3 [NASA ADS] [CrossRef] [Google Scholar]
- Krivova, N. A., & Solanki, S. K. 2004, A&A, 417, 1125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Landi degl'Innocenti, E. 1992, in Solar Observations: Techniques and Interpretation, ed. F. Sánchez, M. Collados, & M. Vázquez (Cambridge: Cambridge University Press), 73 [Google Scholar]
- Molowny-Horas, R., & Yi, Z. 1994, ITA (Oslo) Internal Rep., 31 [Google Scholar]
- Nordlund, A. 1985, Sol. Phys., 100, 209 [NASA ADS] [CrossRef] [Google Scholar]
- November, L. J., & Simon, G. W. 1988, ApJ, 333, 427 [Google Scholar]
- Roudier, Th., & Muller, R. 2004, A&A, 419, 757 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ruiz Cobo, B., & Del Toro Iniesta, J. C. 1992, ApJ, 398, 375 [NASA ADS] [CrossRef] [Google Scholar]
- Scharmer, G. B., Bjelksjö, K., Korhonen, T. K., Lindberg, B., & Petterson, B. 2003, Proc. SPIE., 4853, 341 [Google Scholar]
- Socas Navarro, H. 2001, in ASP Conf. Proc. 236, ed. M. Sigwarth (San Francisco, CA: ASP), 487 [Google Scholar]
- Sturrock, P. A., & Uchida, Y. 1981, ApJ, 246, 331 [NASA ADS] [CrossRef] [Google Scholar]
- Tarbell, T. D., Slater, G. L., Frank, Z. A., et al. 1991, BAAS, 23, 1048 [NASA ADS] [Google Scholar]
- Title, A. M., Tarbell, T. D., Simon, G., & the SOUP team 1986, Adv. Space Res., 6, 253 [NASA ADS] [CrossRef] [Google Scholar]
- Title, A. M., Topka, K. P., Tarbell, T. D., et al. 1992, ApJ, 393, 782 [NASA ADS] [CrossRef] [Google Scholar]
- Tsuneta, S., Ichimoto, K., Katsukawa, Y., et al. 2008, Sol. Phys., 249, 167 [NASA ADS] [CrossRef] [Google Scholar]
- Van Noort, M., Rouppe van der Voort, L., & Löfdahl, M. G. 2005, Sol. Phys., 228, 191 [NASA ADS] [CrossRef] [Google Scholar]
- Vargas Domínguez, S. 2008, Ph.D. Thesis, University of La Laguna - Instituto de Astrofísica de Canarias, Tenerife, Spain [Google Scholar]
- Wang, Y., Noyes, R. W., Tarbell, T. D., & Title, A. M. 1995, ApJ, 447, 419 [NASA ADS] [CrossRef] [Google Scholar]
- Zirker, J. B. 1993, Sol. Phys., 147, 47 [NASA ADS] [CrossRef] [Google Scholar]
Footnotes
- ... respectively
- To identify the lobes, we use the line orthogonal to the linear fit of the pixels belonging to the structure and that passes through the center-of-mass of the whole data cube.
All Figures
![]() |
Figure 1:
Quasi-simultaneous observations showing the evolution of a quiet sun region on 29 September 2007. From top to bottom:
Mg I magnetogram, G-band, CN, Ca II H. Every panel
displays the average image over 4-min intervals and the time stamps
correspond to the initial time for each interval, respectively.
BPs are detected in all filtergrams and are co-spatial with the
magnetic concentrations observed in magnetograms. Spatial scales in
both axis are in arcsec and correspond to solar coordinates (FOV |
Open with DEXTER | |
In the text |
![]() |
Figure 2:
From top to bottom: magnetic field strength and LOS velocity single maps obtained by inverting SP data. The FOV covered is shown in the left lower panel. Positive/negative
|
Open with DEXTER | |
In the text |
![]() |
Figure 3:
Map of horizontal velocities ( FWHM 1
|
Open with DEXTER | |
In the text |
![]() |
Figure 4: Temporal evolution of the magnetic centroid for both lobes (L1 and L2) in Fig. 1. The trajectories followed by the centroids are plotted independently and superimposed on the flow maps computed for SST/G-band images in time intervals: s1 ( left) and s2 ( right). The background false-color image is the divergence field computed from each respective flow map. Alphabet letters correspond to UT times as follows: A-08:40, B-08:48, C-08:56, D-09:00, E-09:04, F-09:08, G-09:12, H-09:21, I-09:24, J-09:28, K-09:32, L-09:36, M-09:40, N-09:44. The FOV is the same boxed in solid line in Fig. 3. |
Open with DEXTER | |
In the text |
Copyright ESO 2010
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.