# Supergranules as Probes of Solar Convection Zone Dynamics

###### Abstract

Supergranules are convection cells seen at the Sun’s surface as a space filling pattern of horizontal flows. While typical supergranules have diameters of about 35 Mm, they exhibit a broad spectrum of sizes from Mm to Mm. Here we show that supergranules of different sizes can be used to probe the rotation rate in the Sun’s outer convection zone. We find that the equatorial rotation rate as a function of depth as measured by global helioseismology matches the equatorial rotation as a function of wavelength for the supergranules. This suggests that supergranules are advected by flows at depths equal to their wavelengths and thus can be used to probe flows at those depths. The supergranule rotation profiles show that the surface shear layer, through which the rotation rate increases inward, extends to depths of Mm and to latitudes of at least . Typical supergranules are well observed at high latitudes and have a range of sizes that extend to greater depths than those typically available for measuring subsurface flows with local helioseismology. These characteristics indicate that probing the solar convection zone dynamics with supergranules can complement the results of helioseismology.

## 1 Introduction

Supergranules were discovered in the 1950s by Hart (1954) but were best characterized in the 1960s by Leighton et al. (1962) who gave them their name. Leighton et al. (1962) showed that this cellular pattern of horizontal flows covers the solar surface and that the boundaries of the cells coincide with the chromospheric/magnetic network. Typical supergranules have diameters of Mm and maximum flow speeds of m s. While the kinetic energy spectrum has a distinct peak at wavelengths of Mm, the spectrum includes cells at least three times larger and extends to much smaller cells where the supergranule spectrum blends into the granulation spectrum (Hathaway et al., 2000).

Larger cells live longer than smaller cells. Typical supergranules with diameters of Mm live for hr (Simon & Leighton, 1964; Wang & Zirin, 1989). Typical granules with diameters of Mm only live for min (Title et al., 1989). Cells of intermediate size ( Mm) have intermediate lifetimes of hr (November et al., 1981).

The rotation rate of the supergranule pattern was first measured by Duvall (1980) who cross-correlated the Doppler velocity pattern from equatorial spectral scans obtained over several days. He found that the pattern rotates about 3% faster than the photospheric plasma and faster rates are found for the 24-hr time lags from day-to-day than for the 8-hr time lags from the beginning to end of an observing day. He concluded that larger cells dominate the longer time lags and that the observations are consistent with supergranules embedded in a surface shear layer in which the rotation rate increases with depth.

The presence of this shear layer was first suggested by Foukal & Jokipii (1975) as a consequence of the conservation of angular momentum by convective elements moving inward and outward in the near surface layers. Global helioseismology inversions (Thompson et al., 1996; Schou et al., 1998) indicate that the shear layer extends to depths of 35-50 Mm near the equator but may disappear or reverse at latitudes above . While local helioseismology (Basu et al., 1999; Corbard & Thompson, 2002) does not probe as deeply, it produces similar results which suggest that the rotation rate may not continue to increase inward at higher latitudes.

Beck & Schou (2000) measured the rotation of the supergranule pattern using a Fourier technique with space-based Doppler data from the ESA/NASA Solar and Heliospheric Observatory (SOHO) Michelson Doppler Imager (MDI) (Scherrer et al., 1995). They mapped the data onto heliographic coordinates, took the Fourier transform in longitude of data from equatorial latitudes, and then the Fourier transform in time of those spectral coefficients over 6 10-day intervals in 1996 which had continuous coverage at a 15-min cadence. They found that the larger cells do indeed rotate more rapidly than the smaller cells, but with rotation rates that exceeded the peak internal rotation rate at the base (50 Mm depth) of the surface shear layer as determined from global helioseismology (Schou et al., 1998). This discovery led them to conclude that supergranules must have wave-like properties in order to rotate faster than the flows they are embedded in.

However, Hathaway et al. (2006) showed that line-of-sight projection effects on a rigidly rotating fixed velocity pattern could reproduce the excess rotation velocities found by Beck & Schou (2000). In projecting the vector velocities onto the line-of-sight the function (where is the heliographic longitude relative to the central meridian) multiplies the longitudinal flow velocities near the equator. Since the flows are largely horizontal the multiplier effectively pushes the peaks in the Doppler pattern away from the central meridian and makes the pattern appear to rotate more rapidly.

In fact, Schou (2003) largely removed the line-of-sight projection effects from the equatorial Doppler data by dividing the data by a function that approximated the function and found that rotation velocities were much more in line with those from global helioseismology (but noted that there were still motions relative to this that suggested wave-like properties for supergranules). Recently, Hathaway et al. (2010) found that the rotation profiles as functions of latitude determined by the cross-correlation technique used by Duvall (1980) for time lags from 2-hr to 16-hr could be reproduced by cellular patterns that are advected by a differential rotation with a peak velocity consistent with that found in the Sun’s surface shear layer by global helioseismology.

Here we measure the rotation of the pattern of supergranules by analyzing the same data used by Beck & Schou (2000) and by Schou (2003). We execute a series of 2D Fourier transform analyses. We remove the line-of-sight projection effects near the equator as was done by Schou (2003) and repeat the Fourier analysis done by Beck & Schou (2000) to show that the equatorial rotation rate of supergranules as a function of wavelength matches the equatorial rotation rate as a function of depth determined from global helioseismology (Schou et al., 1998). This “de-projection” can only be done at the equator and is only approximate since the flows are not purely horizontal. We determine the rotation rates at other latitudes by repeating the Fourier transform analysis on the raw Doppler data (without the removal of projections effects) and using a data simulation to support our conclusions.

## 2 Data Preparation

The data consist of pixel images of the line-of-sight velocity determined from the Doppler shift of a spectral line due to the trace element nickel in the solar atmosphere. The images are acquired at a 1 min cadence and cover the full visible disk of the Sun. We average the data over 31 min with a Gaussian weighting function which filters out variations on time scales less than about 16 min, and sample that data at 15 min intervals. We then map these temporally filtered images onto a grid in heliographic latitude from pole to pole and in longitude from the central meridian (Figure 1). This mapping accounts for the position angle of the Sun’s rotation axis relative to the imaging CCD and the tilt angle of the Sun’s rotation axis toward or away from the spacecraft. Both of these angles include modification in line with the most recent determinations of the orientation of the Sun’s rotation axis (Beck & Giles, 2005; Hathaway & Rightmire, 2010). We analyze data obtained during a 60 day period of continuous coverage in 1996 from May 24 to July 22.

Line-of-sight projection effects influence the results so we also generate and analyze simulated data to assist in our determination of the actual rotation as a function of latitude and depth in the Sun’s outer convection zone. We construct the simulated data from an evolving spectrum of vector spherical harmonics in such a manner as to reproduce the spatial, spectral, and temporal behavior of the observed cellular flows. The amplitudes of the spectral coefficients are constrained by matching the observed velocity spectrum (Hathaway et al., 2000) with the radial flow component constrained by the disk center to limb variation in the RMS Doppler signal(Hathaway et al., 2002).

The cells are given finite lifetimes and made to rotate by adding changes to the phases of the spectral coefficients (Hathaway et al., 2010). The rotation rates are constrained by matching the observed rotation rates as functions of latitude and wavelength.

The rotation rates of the cells in the simulated data are given by a fairly simple function with the latitudinal, , variation separated from the wavelength, , variation such that

(1) |

with

(2) |

and

(3) |

where the wavelength, , is given in Mm.

## 3 Equatorial Rotation Rate

We determine the equatorial rotation rate of the supergranules by 2D Fourier transforms with and without removing line-of-sight projection effects. The line-of-sight projection effects can be minimized using the method described by Schou (2003). The mapped Doppler velocities near the equator are divided by a function, which approximates the geometric factor, , that multiplies the longitudinal velocity in producing the Doppler signal. This “de-projected” signal and the raw Doppler signal are both then apodized near the limb and then Fourier transformed over longitude for the 50 latitude positions that straddle the equator. These spectral coefficients are then Fourier transformed in time over six 10-day intervals.

The rotation rate as a function of wavelength is determined by first finding the temporal frequency of the centroid of the spectral power for that wavelength. This temporal frequency is then divided by the longitudinal wavenumber to give the synodic rotation rate which is then converted to a sidereal rotation rate by adding a correction based on the rate of change of the ecliptic longitude during the observations ().

The results of these analyses are shown in Fig. 2 along with the rotation rate with depth from a global helioseismology analysis by Schou et al. (1998). The sidereal rotation rate of the de-projected supergranules as a function of longitudinal wavelength matches the rotation rate as a function of depth through the outer half of the solar convection zone and the simulation matches the MDI observations. The raw Doppler data give faster rotation rates than the de-projected data at all wavelengths (in both the MDI and the simulated data) with larger increases for larger cells.

The rotation rates seen with both de-projected datasets fall slightly below both the helioseismology results and the input profile for the simulated data. Experiments with the simulated data suggest that this can be attributed to the presence of a small radial flow component. This component is projected into the line-of-sight along the equator by multiplying by a projection factor . When the Doppler signal is de-projected by dividing by this small signal can become large and influence the results in this manner.

We also see that the rotation rates for the raw simulated data fall systematically below the MDI results for wavelengths Mm. This too may be a result of the radial flows.

These results do however lead to a key conclusion - that supergranules with sizes from 10 Mm to 100 Mm are advected by flows within the convection zone at depths equal to their widths.

## 4 Rotation Profiles

We determine the rotation rate of the supergranules as functions of latitude and wavelength (depth) by repeating the 2D Fourier transforms on the raw Doppler data for a series of latitude strips. Each strip is 11 pixels or high in latitude and offset from the previous strip by 5 pixels. We repeat the procedure for 172 positions between latitude and for each of six 10-day intervals for both MDI and simulated data. Latitudinal rotation profiles are obtained for a series of cell wavelengths by averaging the profiles for all waveumbers that produce wavelengths within 5 Mm of the target wavelength. These profiles are then averaged between hemispheres and smoothed with a 9-point binomial smoothing kernal which then limits the data to latitude. The results for averages from the MDI datasets and the from simulated data are shown in Fig. 3.

The rotation velocity appears to increase at all latitudes with increasing wavelength. However, the small increase at 70 Mm wavelength (in particular near the equator) can be attributed to projection effects. The actual rotation rate of the convection cells in the simulation is given by equations 1-3 and the solid lines in Fig 2. which give a rotation rate that increases with wavelength to a maximum at Mm at all latitudes.

The smallest cells, Mm, are difficult to resolve at high latitudes but clearly show slower rotation than larger cells over their observed latitude range. Cells with wavelengths of Mm can be resolved at all latitudes and have significantly lower rotation rates than the larger cells.

While this method with the raw Doppler data is subject to systematic offsets due to line-of-sight projection effects, the same offsets are present in the simulated data and are much smaller for the smaller cells. We conclude that the surface shear layer, in which the rotation rate increases inward, extends to latitudes of at least .

## 5 Conclusions

We conclude that supergranules are anchored or steered in the subsurface flows at depths equal to their wavelengths. This is a simple explanation for the match with the rotation rate with depth from global helioiseismology (Figure 2 lower panel) and is based on well known physics. This is consistent with numerical simulations of convection in the outermost 16 Mm of the Sun by Stein et al. (2011) who show that flow structures at different depths have diameters about equal to the depth itself.

This conclusion is, however, somewhat surprising given the much smaller estimates for the depth of typical supergranules previously determined from the visibility of their internal flows using local helioseismology. Duvall (1998) estimated a depth of 8 Mm for typical supergranules while Zhao & Kosovichev (2003) estimated a depth of 15 Mm. This suggests that local helioseismology is less sensitive to these deeper (and slower) flows and that this new method of probing the convection zone with supergranules can probe flows at greater depths.

We also conclude that the surface shear layer extends to a depth of Mm at all latitudes. The increase in rotation rate with depth has long been suggested by observations and is attributed to the conservation of angular momentum for fluid elements moving inward and outward in the near surface layers (Foukal & Jokipii, 1975; Hathaway, 1982). Measurements of this rotation rate increase from helioseismology (Schou et al., 1998; Corbard & Thompson, 2002) indicate that it follows this critical gradient to depths of 10-15 Mm and reaches a maximum rotation rate at depths of 35-50 Mm. However, many helioseismology results also suggest that the shear layer disappears at latitudes above about . The results reported here, using supergranules, indicate that the shear layer extends to the highest latitudes probed with this MDI data - .

The rotation increase with depth given by Equation 3 follows the critical gradient (with ) given by angular momentum mixing near the surface but then drops below that value at greater depths as the cell turn-over times become longer and the convective flows adjust to the solar rotation.

Global helioseismology can measure the internal rotation rate to great depths but it gives less reliable results at high latitudes. Local helioseismology can measure the non-axisymmetric flows (as well as the axisymmetric meridional flow) but only in the near surface layers. Using supergranules of different sizes to probe the flows in the Sun’s convection zone extends these measurements to greater depths and higher latitudes. This new method of probing solar convection zone dynamics should provide information complementary to that obtained with helioseismology.

## References

- Basu et al. (1999) Basu, S., Antia, H. M., & Tripathy, S. C. 1999, ApJ 512, 458
- Beck & Giles (2005) Beck, J. G., & Giles, P. 2005, ApJ 621, L153
- Beck & Schou (2000) Beck, J. G., & Schou, J. 2000, Sol. Phys. 193, 333
- Corbard & Thompson (2002) Corbard, T., & Thompson, M. J. 2002, Sol. Phys. 205, 211
- Duvall (1980) Duvall, T. L., Jr. 1980, Sol. Phys. 66, 213
- Duvall (1998) Duvall, T. L., Jr. 1998, in Structure and Dynamics of the Interior of the Sun and Sun-like Stars ed. S. Korzennik (Noordwijk, Netherlands: ESA SP 418), 581
- Foukal & Jokipii (1975) Foukal, P., & Jokipii, R. 1975, ApJ 199, L71
- Hart (1954) Hart, A. B. 1954, MNRAS 114, 17
- Hathaway (1982) Hathaway, D. H. 1982, Sol. Phys. 77, 341
- Hathaway et al. (2000) Hathaway, D. H., Beck, J. G., Bogart, R. S., Bachmann, K. T., Khatri, G., Petitto, J. M., Han, S., & Raymond, J. 2000, Sol. Phys. 193, 299
- Hathaway et al. (2002) Hathaway, D. H., Beck, J. G., Han, S., & Raymond, J. 2002, Sol. Phys. 205, 25
- Hathaway & Rightmire (2010) Hathaway, D. H. & Rightmire, L. 2010, Science 327, 1350
- Hathaway et al. (2006) Hathaway, D. H., Williams, & Cuntz, M. 2006, ApJ 644, 598
- Hathaway et al. (2010) Hathaway, D. H., Williams, P. E., Dela Rosa, K., & Cuntz, M. 2010, ApJ 725, 1082
- Leighton et al. (1962) Leighton, R. B., Noyes, R. W., & Simon, G. W. 1962, ApJ 135, 474
- November et al. (1981) November, L. J., Toomre, J., Gebbie, K. B. & Simon, G. W. 1981, ApJ 245, L123
- Scherrer et al. (1995) Scherrer, P. H., et al. 1995, Sol. Phys. 162, 129
- Schou et al. (1998) Schou, J. et al. 1998, ApJ 505, 390
- Schou (2003) Schou, J. 2003, ApJ 596, L259
- Simon & Leighton (1964) Simon, G. W. & Leighton, R. B. 1964, ApJ 140, 1120
- Stein et al. (2011) Stein, R. F., LagerFjärd, A., Nordlund, Å., & Georgobiani, D. 2011, Sol. Phys. 268, 271
- Thompson et al. (1996) Thompson, M. J., et al. 1996, Science 272, 1300
- Title et al. (1989) Title, A. M., Tarbell, T. D., Topka, K. P., Ferguson, S. H., Shine, R. A. & The SOUP Team 1989, ApJ 336, 475
- Wang & Zirin (1989) Wang, H. & Zirin, H. 1989, Sol. Phys. 120, 1
- Zhao & Kosovichev (2003) Zhao, J. & Kosovichev, A. G. 2003, in GONG+ 2002. Local and Global Helioseismology: the Present and Future ed. H. Sawaya-Lacoste (Noordwijk, Netherlands: ESA SP 517), 417