Issue 
A&A
Volume 570, October 2014



Article Number  A90  
Number of page(s)  7  
Section  The Sun  
DOI  https://doi.org/10.1051/00046361/201424201  
Published online  23 October 2014 
Timedistance helioseismology: A new averaging scheme for measuring flow vorticity
^{1}
GeorgAugustUniversität, Institut für Astrophysik,
FriedrichHundPlatz 1,
37077
Göttingen, Germany
email:
jlangfe@astro.physik.unigoettingen.de
^{2}
MaxPlanckInstitut für Sonnensystemforschung,
JustusvonLiebigWeg 3,
37077
Göttingen,
Germany
Received:
13
May
2014
Accepted:
12
August
2014
Context. Timedistance helioseismology provides information about vector flows in the nearsurface layers of the Sun by measuring wave travel times between points on the solar surface. Specific spatial averages of travel times have been proposed for distinguishing between flows in the eastwest and northsouth directions and measuring the horizontal divergence of the flows. No specific measurement technique has, however, been developed to measure flow vorticity.
Aims. Here we propose a new measurement technique tailored to measuring the vertical component of vorticity. Fluid vorticity is a fundamental property of solar convection zone dynamics and of rotating turbulent convection in particular.
Methods. The method consists of measuring the travel time of waves along a closed contour on the solar surface in order to approximate the circulation of the flow along this contour. Vertical vorticity is related to the difference between clockwise and anticlockwise travel times.
Results. We applied the method to characterize the vortical motions of solar convection using helioseismic data from the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory (SDO/HMI) and from the Michelson Doppler Imager onboard the Solar and Heliospheric Observatory (SOHO/MDI). Away from the equator, a clear correlation between vertical vorticity and horizontal divergence is detected. Horizontal outflows are associated with negative vorticity in the northern hemisphere and positive vorticity in the southern hemisphere. The signal is much stronger for HMI than for MDI observations. We characterize the spatial power spectrum of the signal by comparison with a noise model. Vertical vorticity at horizontal wavenumbers below 250 /R_{⊙} can be probed with this helioseismic technique.
Key words: Sun: helioseismology / Sun: oscillations / convection
© ESO, 2014
1. Introduction
The Sun exhibits complex flow patterns in the convection zone, such as turbulent convection, differential rotation, and meridional circulation. These flows are important ingredients for understanding global solar dynamics and the dynamo responsible for the solar 22year magnetic cycle (cf. Toomre 2002). Fluid vorticity is a fundamental characteristic of fluid dynamics. The interplay between turbulent convection and rotation can generate cyclonic motions with a net kinetic helicity that depends on solar latitude (Duvall & Gizon 2000). These motions may convert the toroidal magnetic field into a poloidal field (Parker 1979). Vortices are not confined to convective motions. Hindman et al. (2009) detected that inflows into active regions (Gizon et al. 2001) have a cyclonic component that is presumably caused by solar rotation.
Duvall et al. (1993) showed that nearsurface solar flows can be measured using timedistance helioseismology. The idea is to measure the time it takes for solar waves to travel between two surface locations from the temporal crosscovariance of the observable measured at these locations. Typically, the observable is a series of lineofsight Doppler velocity images, φ(r,t), which has been filtered in the Fourier domain to select particular wave packets. We consider a pair of points r_{1} and r_{2} (“pointtopoint geometry”). The crosscovariance C at time lag t is (1)where h_{t} is the temporal cadence, T = (2N + 1)h_{t} the observation time, and t_{i} = ih_{t} with i = −N, − N +1,...,N are the times at which the observable is sampled. From the crosscovariance, the travel time can be measured by fitting a wavelet (e.g., Duvall et al. 1997) or a sliding reference crosscovariance (Gizon & Birch 2004).
Waves are advected by the flow field v(r), and travel times are sensitive to flows in the vicinity of the ray connecting the points r_{1} and r_{2}. If the flow has a component in the direction r_{2} − r_{1}, then the travel time from r_{1} to r_{2} (denoted by τ^{+}(r_{1},r_{2})) is reduced, while the travel time from r_{2} to r_{1} (denoted by τ^{−}(r_{1},r_{2})) is increased.
To obtain a measurement that is particularly sensitive to the horizontal flow divergence div_{h}, travel times are measured between a central point r and a surrounding annulus with radius Δ (Duvall et al. 1993). This “pointtoannulus geometry” is displayed in Fig. 1a. The flow divergence is related to the difference between inward and outward travel times.
Duvall et al. (1997) proposed to break the annulus into four quadrants pointing in the east, west, north, and south directions, respectively. Here we remind the reader that the solar convention is that west is in the prograde direction of solar rotation. The travel time measured between r and the west (or the east) quadrant (“pointtoquadrant geometry”) is sensitive to the component of the flow velocity in the west direction, v_{x}. In practice, the difference of the quadrants is used. In the same fashion, the north component of the flow velocity, v_{y}, can be obtained using the north and south quadrants.
There is no specific measurement geometry, however, that is directly sensitive to the flow vorticity. So far, the vertical component of flow vorticity, ω_{z} = ∂_{x}v_{y} − ∂_{y}v_{x}, has been estimated by taking spatial derivatives of the westeast and northsouth travel times (see, e.g., Gizon et al. 2000). Alternatively, one could take the spatial derivatives of inverted flow velocities. We would like though a traveltime measurement that is close to the vorticity before performing any inversion. Furthermore, taking derivatives of noisy quantities (as in both cases above) is a dangerous operation. Thus it is desirable to have a traveltime measurement geometry that is explicitly tailored to measure vorticity and that avoids numerical derivatives.
Fig. 1 Traveltime measurement geometries. a) Measurement geometry sensitive to the horizontal component of the flow divergence. Travel times are measured between a central point r and the average over a surrounding annulus with radius Δ as introduced by Duvall et al. (1993). b) Proposed new measurement geometry sensitive to the vertical component of flow vorticity. Travel times are measured sequentially along neighboring pairs of points r_{i} and r_{i + 1} located on a closed contour. In this example, n = 6 points are used, forming a regular hexagon. c) As b), but the hexagon is rotated by an angle β around r. 
2. Measuring vortical flows along a closed contour
In this paper, we implement a measurement technique that is used in ocean acoustic tomography (Munk et al. 1995), where wave travel times are measured along a closed contour . This measurement returns the flow circulation along the contour. The flow circulation is related to the vertical component of vorticity ω_{z} (averaged over the area A enclosed by the contour) by Stokes’ theorem: (2)where v is the flow velocity vector on the surface A. The vector dA is normal to the solar surface (upward) and the contour integral runs anticlockwise with dl tangential to the contour.
2.1. Geometry for anticlockwise travel times
We approximate the contour integral as follows. We select n points r_{1}, r_{2}, ..., r_{n} along a circular contour and measure the travel times τ^{+} pairwise in the anticlockwise direction. Neighboring points are each separated by an equal distance Δ. The points form the vertices of a regular polygon (Fig. 1b). Averaging over the τ^{+} measurements yields what we call the “anticlockwise travel time” , (3)with the notation r_{n + 1} = r_{1}. With this definition, is reduced when there is a flow velocity v_{φ} tangential to the circle of radius R = Δ / [ 2sin(π/n) ] in anticlockwise direction. In order to provide a simplified description of the relationship between and ω_{z}, we may write the perturbation to the anticlockwise travel time caused by v_{φ} as (4)where τ_{0} is the unperturbed travel time and V_{ref} is the reference wave speed. This description only provides a ruleofthumb connection between and ω_{z}. The proper relationship is described by 3D sensitivity kernels (Birch & Gizon 2007). Additionally, these quantities are functions of all three spatial dimensions, which has not been accounted for.
We note that the distance Δ between the points must be greater than the wavelength (e.g., about 5 Mm for the f mode at 3 mHz), in order to distinguish waves propagating from r_{i} to r_{i + 1} from waves propagating in the opposite direction. Also there is some freedom in selecting the contour C. For example, active regions are often shaped irregularly. To measure the vorticity around active regions, it may be useful to adapt the contour to the shape of the active region.
2.2. Reducing the noise level
The definition of the anticlockwise travel times (Eq. (3)) assumes that pairwise travel times can be measured, irrespective of the noise level. This is not a problem in the quiet Sun using the traveltime definition of Gizon & Birch (2004), which is very robust with respect to noise. However, a wavelet fit to the crosscovariance (as in Duvall et al. 1997) is only possible when the noise is sufficiently low. This is not the case for a single pair of points (see Fig. 2a for an example using fmodefiltered SDO/HMI data with Δ = 10 Mm).
One option is to average C before performing the wavelet fit. At fixed Δ and n, the measurement polygon can be rotated by angles β around r (Fig. 1c). Since plane waves are only weakly correlated for different propagation directions, taking the average over various angles β will lower the noise level. In Fig. 2c, an average over eight angles β is shown for n = 6. Furthermore, C can be averaged over several annulus radii R (several n at fixed distance Δ). In Fig. 2d, the crosscovariance is averaged over three different annuli (n = 4, 6, and 8) and angles β. An additional 4 × 4 averaging over the centers of annuli (Fig. 2e) gives a crosscovariance function that has a sufficiently low level of noise to be analyzed by a wavelet fit. Such averaging procedures are often used for measuring outward−inward travel times. Any spatial averaging must be properly taken into account when traveltime inversions are performed later.
Fig. 2 Fmode crosscovariance C as a function of time lag for HMI data (see Sect. 3.1). All curves except e) and f) are for a single center point r near disk center (8 h averaging, 1 May 2010 dataset). The distance Δ between consecutive points is 10 Mm. a) C for a single pair. b) Averaging over annulus with n = 6. c) Further averaging over eight angles β. d) Further averaging over two additional radii: n = 4 with four angles β and n = 8 with eight angles β. e) Further averaging over annuli centers (4 × 4 pixels). f) For reference, n = 6 annulus averaged over an entire map (about 180 × 180 Mm^{2}). 
2.3. Decoupling from isotropic wavespeed perturbations
Since the purpose of is to measure the vorticity ω_{z}, it should ideally not be sensitive to any other solar perturbation. However, a wavespeed perturbation along the contour due to fluctuations in, for instance, temperature or density will also perturb . In order to remove traveltime perturbations that are not caused by the flow, we also measure the clockwise travel time (5)For example, the travel times and are affected in the same way by a temperature perturbation. We thus introduce the difference between the anticlockwise and the clockwise travel times (denoted with the superscript “ac”), (6)Note that , where is given by Eq. (4). The travel time τ^{ac} should be largely independent of perturbations other than a vortical flow. This approach is similar to the one proposed by Duvall et al. (1997) to measure the flow divergence from pointtoannulus travel times.
3. Proof of concept using SDO/HMI and SOHO/MDI observations
In order to test if the proposed averaging scheme is able to measure flow vorticity, we have carried out two simple experiments using f modes. The first experiment (Sect. 3.3) consists of making maps of using SDO/HMI observations (Schou et al. 2012), computing the spatial power spectrum of these maps, and comparing with the predicted power spectrum of pure realization noise. The second experiment (Sect. 3.4) is to look for a correlation between vertical vorticity and horizontal divergence. The sign and the amplitude of this correlation are expected to scale like the local Coriolis number.
3.1. Observations
We used 112 × 24 h series of SDO/HMI lineofsight velocity images. The Dopplergrams were taken from 1 May to 28 August 2010 when the Sun was relatively quiet. Regions of the size 180 × 180 Mm^{2} at solar latitudes from −60° to +60° in steps of 20° were tracked for one day as they crossed the central meridian. Images were remapped using Postel’s projection and tracked at the local surface rotation rate from Snodgrass (1984). The resulting data cubes were cut into three 8 h datasets. A ridge filter was applied to select f modes.
We also used 56 × 24 h series of SOHO/MDI (Scherrer et al. 1995) fulldisk lineofsight velocity images from 8 May through 11 July 2010, thus overlapping with the HMI observations. The MDI data were processed in the same way as for HMI, however the spatial sampling of MDI is lower by a factor of four (2.0 arcsec px^{1} instead of 0.5 arcsec px^{1}) and the temporal cadence is 60 s instead of 45 s. The spatial resolution is given by the instrumental point spread function (PSF), which can be approximated by a Gaussian with a full width at half maximum (FWHM) of about 3.3 arcsec (2.4 Mm) for MDI (Korzennik et al. 2004, 2012) and about 1.0 arcsec (0.7 Mm) for HMI (Yeo et al. 2014).
3.2. Traveltime maps
From the fmodefiltered Dopplergrams, we computed the crosscovariance C in Fourier space, (7)where ω is the angular frequency and h_{ω} = 2π/T is the frequency resolution. The symbol φ_{T} denotes the Fourier transform of φ(r,t) multiplied by the temporal window function. This way of computing C is equivalent to Eq. (1), but is much faster. To measure the travel times and , we used the linearized definition of travel times as defined by Eq. (3) in Gizon & Birch (2004) with W given by Eq. (4) in that paper. We obtained the reference crosscovariance C^{ref} by spatially averaging C over the whole map. We used Δ = 10 Mm and n = 6, which corresponds to an annulus radius of 10 Mm, and used four different values for β (0°, 15°, 30°, and 45°). We computed τ^{ac} as defined in Eq. (6). Additionally, we computed outward−inward mean traveltime maps (τ^{oi}) between an annulus radius of 10 Mm and the central point. Again, we used the traveltime definition from Gizon & Birch (2004).
Figure 3a shows a τ^{oi} traveltime map for an example 8 h dataset at the solar equator. The bluish features of size 20 to 30 Mm are areas of positive divergence. They represent supergranular outflow regions. Conversely, the reddish areas show the supergranular network of converging flows. For the same dataset, a τ^{ac} map that was averaged over the four angles β is depicted in Fig. 3b. There is no evidence of excess power at the scales of supergranulation.
Fig. 3 a) Example HMI τ^{oi} traveltime map for one fmodefiltered 8 h dataset centered at the solar equator and annulus radius Δ = 10 Mm. The coordinate x is in the west direction and y in the north direction. The travel times have been measured using the linearized definition in Gizon & Birch (2004). The color scale has been truncated to increase the contrast. The minimum and maximum values are −152.4 s and 117.3 s. b) Example τ^{ac} traveltime map for the same dataset as in a) with Δ = 10 Mm and n = 6, averaged over the four angles β. The minimum and maximum values are − 149.7 s and 144.8 s. c) Power spectrum of τ^{oi} for both HMI data and noise model (Gizon & Birch 2004), averaged over azimuth and 336 × 8 h datasets (size about 180 × 180 Mm^{2}) centered at the solar equator, plotted versus the product of horizontal wavenumber k and solar radius R_{⊙}. The thickness of the lines denotes the 1σ error. d) Power spectrum of τ^{ac} for both HMI data and noise model, averaged over azimuth and 336 × 8 h datasets. The τ^{ac} maps were averaged over four angles β before computing the power. 
3.3. Test 1: evidence of a vorticity signal in τ^{ac} as a function of wavenumber
In order to evaluate the signaltonoise ratio (S/N) in the τ^{ac} maps, we compare the spatial power of the traveltime maps with a noise model. For the noise model, we use the recipe of Gizon & Birch (2004) to construct artificial datasets as follows. In 3D Fourier space the observable is modeled by a Gaussian complex random variable with zero mean and variance given by the expected power spectrum (estimated from the observations). In the noise model, wavenumbers and frequencies are uncorrelated to exclude wave scattering by flows and heterogeneities (signal). The expectation value of the power spectrum is chosen to match that of HMI observations. A detailed study and a validation of the noise model is provided by Fournier et al. (2014).
For each HMI dataset, one realization of the noise model was generated, based on the corresponding power spectrum. From these noise datasets, we computed τ^{oi} and τ^{ac} traveltime maps in the same manner as for the HMI observations (cf. Sect. 3.2).
The τ^{oi} power spectrum averaged over azimuth and all datasets is shown in Fig. 3c. There is signal above noise level for kR_{⊙}< 500, with a maximum S/N of 50 at kR_{⊙} = 120. This is the wellknown supergranulation peak. For high wavenumbers (kR_{⊙}> 500), the traveltime maps are dominated by noise. The contribution from convection features that are much shorter lived than the observation time of 8 h is very small. At kR_{⊙} ~ 2000 there is a cutoff in power corresponding to a wavelength of 2.5 Mm (half the fmode wavelength at 3 mHz). For MDI (not shown here), the S/N at the supergranulation peak is about 16 and the S/N vanishes at kR_{⊙} = 400.
For τ^{ac}, we averaged the maps over four angles β, computed the power spectrum for each resulting map and averaged the power spectra over azimuth and all datasets. The result for HMI is shown in Fig. 3d. The power for the HMI travel times significantly exceeds that of the noise model for kR_{⊙} < 250, with a S/N increasing toward larger scales (S/N about 1.5 at kR_{⊙} = 120) and reaching a S/N of about 2.6 for kR_{⊙} = 25. This is qualitatively different than for the τ^{oi} power. Again, there is a power cutoff at kR_{⊙} ~ 2000. In the MDI case, the S/N is similar but lower than for the HMI data (S/N = 0.7 at kR_{⊙} = 120, S/N = 1.3 at kR_{⊙} = 25).
For the results presented in Fig. 3d as much as four months of HMI data were averaged. What is the minimum number of days of observations needed to achieve a clear detection of the vorticity signal in the τ^{ac} maps (180 Mm on the side)? By clear detection we mean that, at fixed wavenumber, the power in the observed travel times and the power in the noise model are separated by at least 2 × 3σ. This requirement is somewhat arbitrary but safe. Figure 4a shows the spatial power of τ^{ac} at kR_{⊙} = 109 as a function of observation duration, for τ^{ac} maps averaged over four angles β. Overplotted are the 3σ error estimates (filled areas). Since the data cubes are not correlated from one day to the next (different longitudes), the variance of the noise decreases like 1 over the number of days. The criterion for clear signal detection is fulfilled after two days of observations.
Fig. 4 a) Azimuthally averaged power of HMI and noisemodel τ^{ac} maps (average maps over four angles β) at kR_{⊙} = 109, as a function of days of data averaged over (the same data as in Fig. 3 were used). The power was averaged over three 8 h datasets per day. The filled areas denote the 3σ error estimates. b) Variance of τ^{ac} measurements computed from both HMI and noisemodel fmodefiltered data at the solar equator. The observation time of the datasets is 8 h and their spatial size is about 180 × 180 Mm^{2}. The variance is shown as a function of the number of angles β over which the data is averaged. For reference, the black line shows Var(τ^{ac}) ∝ 1 /n_{β}. 
This detection is also a function of the number of angles β over which the contours are rotated and averaged. In Fig. 4b, the variance of the τ^{ac} travel times is plotted versus this number, n_{β}. For Δ = 10 Mm, we find that the variance decreases almost like 1 /n_{β} for small n_{β} and reaches a plateau for n_{β}> 8. This is because the τ^{ac} measurements are highly correlated for small rotation angles. In this case, n_{β} = 4 is a good compromise between efficient data use and computation time.
3.4. Test 2: effect of rotation on vorticity in supergranules
Here we compute ⟨ τ^{oi}τ^{ac} ⟩, where the angle brackets denote an average over the solar surface and over all datasets. Since τ^{oi} ∝ − div_{h} and τ^{ac} ∝ − ω_{z}, the product ⟨ τ^{oi}τ^{ac} ⟩ serves as a proxy for ⟨ div_{h}ω_{z} ⟩, which is a component of the kinetic helicity and is sensitive to the effect of the Coriolis force on convection (e.g., Zeldovich et al. 1990; Rüdiger et al. 1999).
The latitudinal dependencies of and for HMI data are plotted in Fig. 5a. Note that we filtered the traveltime maps spatially by removing the power for kR_{⊙}> 300 since there is no significant signal in τ^{ac} at high wavenumbers. is negative in the northern hemisphere and positive in the southern hemisphere. For the pattern is reversed. At the equator, both quantities have a low positive value of 2 s^{2}, which is probably due to wavespeed perturbations associated with the magnetic network.
Fig. 5 a) Horizontal averages and for fmodefiltered HMI observations as functions of solar latitude, averaged over 336 × 8 h datasets of about 180 × 180 Mm^{2} and four angles β. b) Horizontal average ⟨ τ^{oi}τ^{ac} ⟩ for different data: HMI at full resolution (0.5 arcsec px^{1}), MDI fulldisk data (2.0 arcsec px^{1}), and HMI data spatially averaged over 4 × 4 px (after remapping) and convolved with a Gaussian with 2.4 Mm FWHM to match the MDI sampling and PSF. 
In Fig. 5b, ⟨ τ^{oi}τ^{ac} ⟩ is plotted versus solar latitude for both HMI and MDI data. This allows to compare the sensitivity of both instruments. Additionally, ⟨ τ^{oi}τ^{ac} ⟩ was computed from HMI Dopplergrams that were binned over 4 × 4 pixels and convolved with a Gaussian with 2.4 Mm FWHM (black curve) such that the spatial resolution matches approximately the resolution of MDI fulldisk data. For both instruments, ⟨ τ^{oi}τ^{ac} ⟩ is negative in the northern hemisphere and positive in the southern hemisphere. As expected, ⟨ τ^{oi}τ^{ac} ⟩ is consistent with zero at the equator (no Coriolis force acts on horizontal flows). Our measurements are qualitatively consistent with the ⟨ div_{h}ω_{z} ⟩ estimates by Duvall & Gizon (2000) and Gizon & Duvall (2003).
Away from the equator, we see that the amplitude of ⟨ τ^{oi}τ^{ac} ⟩ is higher for HMI than for MDI by a factor of five. This means HMI has a much higher sensitivity to ⟨ τ^{oi}τ^{ac} ⟩ than MDI. Only when the HMI Dopplergrams are degraded to match the MDI sampling and PSF, a good agreement of the curves is achieved. Hence the difference in sensitivity to vertical vorticity between HMI and MDI is directly related to the spatial resolution of the instruments.
As in Sect. 3.3, we investigate how much data is needed for a clear nonzero detection (3σ level) of ⟨ τ^{oi}τ^{ac} ⟩ away from the equator. Figure 6 shows ⟨ τ^{oi}τ^{ac} ⟩ versus the duration of the HMI observations at solar latitudes ± 40°. Again the 3σ error estimates are plotted as shaded areas. The curves for observations at 40° and −40° latitude do not overlap after averaging over one or more days of data. This means a difference in ⟨ τ^{oi}τ^{ac} ⟩ between the hemispheres can be detected at the 3σ level after one day of averaging. To distinguish between 40° observations and the noisemodel data takes roughly three days of data. Note that these results are only valid for the given map size of about 180 × 180 Mm^{2}. This corresponds to averaging over roughly 40 supergranules per day.
Fig. 6 Product ⟨ τ^{oi}τ^{ac} ⟩ as a function of days of data averaged over, for fmodefiltered HMI observations at − 40° and 40° latitude as well as τ^{ac} computed for the noise model. The filled areas denote the 3σ estimates. All data have been averaged over space (about 180 × 180 Mm^{2}), four angles β and three 8 h datasets per day. 
4. Conclusion
We have presented a new averaging scheme for timedistance helioseismology, which has direct sensitivity to the vertical component of the flow vorticity. The anticlockwise minus clockwise HMI traveltime maps for f modes show power above the noise level. Unlike the divergence signal, the vorticity signal does not peak at supergranular scales but increases continuously toward larger spatial scales. Furthermore, the latitudinal dependence of the correlation between the vorticity and the divergence signals is consistent with the effect of the Coriolis force on turbulent convection. We find that HMI has a much higher sensitivity to this correlation than MDI.
Acknowledgments
We acknowledge research funding by Deutsche Forschungsgemeinschaft (DFG) under grant SFB 963/1 “Astrophysical flow instabilities and turbulence” (Project A1). The HMI data used are courtesy of NASA/SDO and the HMI science team. The data were processed at the German Data Center for SDO (GDCSDO), funded by the German Aerospace Center (DLR). We are grateful to C. Lindsey, H. Schunker and R. Burston for providing help with tracking and mapping. We acknowledge the workflow management system Pegasus (funded by The National Science Foundation under OCI SI2SSI program grant #1148515 and the OCI SDCI program grant #0722019).
References
 Birch, A. C., & Gizon, L. 2007, Astron. Nachr., 328, 228 [NASA ADS] [CrossRef] [Google Scholar]
 Duvall, Jr., T. L., & Gizon, L. 2000, Sol. Phys., 192, 177 [NASA ADS] [CrossRef] [Google Scholar]
 Duvall, Jr., T. L., Jefferies, S. M., Harvey, J. W., & Pomerantz, M. A. 1993, Nature, 362, 430 [NASA ADS] [CrossRef] [Google Scholar]
 Duvall, Jr., T. L., Kosovichev, A. G., Scherrer, P. H., et al. 1997, Sol. Phys., 170, 63 [NASA ADS] [CrossRef] [Google Scholar]
 Fournier, D., Gizon, L., Hohage, T., & Birch, A. C. 2014, A&A, 567, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gizon, L., & Birch, A. C. 2004, ApJ, 614, 472 [NASA ADS] [CrossRef] [Google Scholar]
 Gizon, L., & Duvall, Jr., T. L. 2003, in GONG+ 2002. Local and Global Helioseismology: the Present and Future, ed. H. SawayaLacoste, ESA SP, 517, 43 [Google Scholar]
 Gizon, L., Duvall, Jr., T. L., & Larsen, R. M. 2000, J. Astrophys. Astron., 21, 339 [Google Scholar]
 Gizon, L., Duvall, Jr., T. L., & Larsen, R. M. 2001, in Recent Insights into the Physics of the Sun and Heliosphere: Highlights from SOHO and Other Space Missions, eds. P. Brekke, B. Fleck, & J. B. Gurman, IAU Symp., 203, 189 [Google Scholar]
 Hindman, B. W., Haber, D. A., & Toomre, J. 2009, ApJ, 698, 1749 [NASA ADS] [CrossRef] [Google Scholar]
 Korzennik, S. G., RabelloSoares, M. C., & Schou, J. 2004, ApJ, 602, 481 [NASA ADS] [CrossRef] [Google Scholar]
 Korzennik, S. G., RabelloSoares, M. C., & Schou, J. 2012, ApJ, 760, 156 [NASA ADS] [CrossRef] [Google Scholar]
 Munk, W. H., Worcester, P., & Wunsch, C. 1995, Ocean acoustic tomography (Cambridge University Press) [Google Scholar]
 Parker, E. N. 1979, Cosmic Magnetic Fields (Oxford) [Google Scholar]
 Rüdiger, G., Brandenburg, A., & Pipin, V. V. 1999, Astron. Nachr., 320, 135 [Google Scholar]
 Scherrer, P. H., Bogart, R. S., Bush, R. I., et al. 1995, Sol. Phys., 162, 129 [Google Scholar]
 Schou, J., Scherrer, P. H., Bush, R. I., et al. 2012, Sol. Phys., 275, 229 [Google Scholar]
 Snodgrass, H. B. 1984, Sol. Phys., 94, 13 [NASA ADS] [CrossRef] [Google Scholar]
 Toomre, J. 2002, Science, 296, 64 [CrossRef] [Google Scholar]
 Yeo, K. L., Feller, A., Solanki, S. K., et al. 2014, A&A, 561, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zeldovich, Y. B., Ruzmaikin, A. A., & Sokoloff, D. D. 1990, Magnetic Fields in Astrophysics (Gordon and Breach) [Google Scholar]
All Figures
Fig. 1 Traveltime measurement geometries. a) Measurement geometry sensitive to the horizontal component of the flow divergence. Travel times are measured between a central point r and the average over a surrounding annulus with radius Δ as introduced by Duvall et al. (1993). b) Proposed new measurement geometry sensitive to the vertical component of flow vorticity. Travel times are measured sequentially along neighboring pairs of points r_{i} and r_{i + 1} located on a closed contour. In this example, n = 6 points are used, forming a regular hexagon. c) As b), but the hexagon is rotated by an angle β around r. 

In the text 
Fig. 2 Fmode crosscovariance C as a function of time lag for HMI data (see Sect. 3.1). All curves except e) and f) are for a single center point r near disk center (8 h averaging, 1 May 2010 dataset). The distance Δ between consecutive points is 10 Mm. a) C for a single pair. b) Averaging over annulus with n = 6. c) Further averaging over eight angles β. d) Further averaging over two additional radii: n = 4 with four angles β and n = 8 with eight angles β. e) Further averaging over annuli centers (4 × 4 pixels). f) For reference, n = 6 annulus averaged over an entire map (about 180 × 180 Mm^{2}). 

In the text 
Fig. 3 a) Example HMI τ^{oi} traveltime map for one fmodefiltered 8 h dataset centered at the solar equator and annulus radius Δ = 10 Mm. The coordinate x is in the west direction and y in the north direction. The travel times have been measured using the linearized definition in Gizon & Birch (2004). The color scale has been truncated to increase the contrast. The minimum and maximum values are −152.4 s and 117.3 s. b) Example τ^{ac} traveltime map for the same dataset as in a) with Δ = 10 Mm and n = 6, averaged over the four angles β. The minimum and maximum values are − 149.7 s and 144.8 s. c) Power spectrum of τ^{oi} for both HMI data and noise model (Gizon & Birch 2004), averaged over azimuth and 336 × 8 h datasets (size about 180 × 180 Mm^{2}) centered at the solar equator, plotted versus the product of horizontal wavenumber k and solar radius R_{⊙}. The thickness of the lines denotes the 1σ error. d) Power spectrum of τ^{ac} for both HMI data and noise model, averaged over azimuth and 336 × 8 h datasets. The τ^{ac} maps were averaged over four angles β before computing the power. 

In the text 
Fig. 4 a) Azimuthally averaged power of HMI and noisemodel τ^{ac} maps (average maps over four angles β) at kR_{⊙} = 109, as a function of days of data averaged over (the same data as in Fig. 3 were used). The power was averaged over three 8 h datasets per day. The filled areas denote the 3σ error estimates. b) Variance of τ^{ac} measurements computed from both HMI and noisemodel fmodefiltered data at the solar equator. The observation time of the datasets is 8 h and their spatial size is about 180 × 180 Mm^{2}. The variance is shown as a function of the number of angles β over which the data is averaged. For reference, the black line shows Var(τ^{ac}) ∝ 1 /n_{β}. 

In the text 
Fig. 5 a) Horizontal averages and for fmodefiltered HMI observations as functions of solar latitude, averaged over 336 × 8 h datasets of about 180 × 180 Mm^{2} and four angles β. b) Horizontal average ⟨ τ^{oi}τ^{ac} ⟩ for different data: HMI at full resolution (0.5 arcsec px^{1}), MDI fulldisk data (2.0 arcsec px^{1}), and HMI data spatially averaged over 4 × 4 px (after remapping) and convolved with a Gaussian with 2.4 Mm FWHM to match the MDI sampling and PSF. 

In the text 
Fig. 6 Product ⟨ τ^{oi}τ^{ac} ⟩ as a function of days of data averaged over, for fmodefiltered HMI observations at − 40° and 40° latitude as well as τ^{ac} computed for the noise model. The filled areas denote the 3σ estimates. All data have been averaged over space (about 180 × 180 Mm^{2}), four angles β and three 8 h datasets per day. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext 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 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.