Meridional flow in the solar interior plays an important role in redistributing angular momentum and transporting magnetic flux inside the Sun. Although it has long been recognized that the meridional flow is predominantly poleward at the Sun's surface and in its shallow interior, the location of the equatorward return flow and the meridional flow profile in the deeper interior remain unclear. Using the first two years of continuous helioseismology observations from the Solar Dynamics Observatory / Helioseismic Magnetic Imager, we analyze travel times of acoustic waves that propagate through different depths of the solar interior carrying information about the solar interior dynamics. After removing a systematic center-to-limb effect in the helioseismic measurements and performing inversions for flow speed, we find that the poleward meridional flow of a speed of 15 m/s extends in depth from the photosphere to about 0.91 R_sun. An equatorward flow of a speed of 10 m/s is found between 0.82 to 0.91 R_sun in the middle of the convection zone. Our analysis also shows evidence of that the meridional flow turns poleward again below 0.82 R_sun, indicating an existence of a second meridional circulation cell below the shallower one. This double-cell meridional circulation profile with an equatorward flow shallower than previously thought suggests a rethinking of how magnetic field is generated and redistributed inside the Sun.

As large--distance rays (say, 10\,-\,$24 ^\circ$) approach the solar surface approximately vertically, travel times measured from surface pairs for these large separations are mostly sensitive to vertical flows, at least for shallow flows within a few Mm of the solar surface. All previous analyses of supergranulation have used smaller separations and have been hampered by the difficulty of separating the horizontal and vertical flow components. We find that the large separation travel times associated with supergranulation cannot be studied using the standard phase-speed filters of time-distance helioseismology. These filters, whose use is based upon a refractive model of the perturbations, reduce the resultant travel time signal by at least an order of magnitude at some distances. More effective filters are derived. Modeling suggests that the center--annulus travel time difference $[\delta t_{\rm{oi}}]$ in the separation range $\Delta=10$\,-\,$24 ^\circ$ is insensitive to the horizontally diverging flow from the centers of the supergranules and should lead to a constant signal from the vertical flow. Our measurement of this quantity, $5.1 \pm 0.1\secs$, is constant over the distance range. This magnitude of signal cannot be caused by the level of upflow at cell centers seen at the photosphere of $10\ms$ extended in depth. It requires the vertical flow to increase with depth. A simple Gaussian model of the increase with depth implies a peak upward flow of $240\ms$ at a depth of $2.3\Mm$ and a peak horizontal flow of $700\ms$ at a depth of $1.6\Mm$.

Measuring the depth variation of the meridional flows is important for understanding the solar cycle, at least according to a number of dynamo models. While attempting to extend the early observations of \citet{giles_thesis} of time-distance measurements of flow, we have stumbled upon some systematic errors that can affect these measurements: 1) the additional distance traveled by radiation coming from points away from disk center causes an apparent `shrinking' Sun, that is an apparent flow towards the disk center and 2) in measurements away from the central longitude, the rotation signal can leak into meridional flow signals. Attempts to understand and overcome these systematic problems will be presented. Forward models based on ray theory have been applied in order to test the sensitivity of travel times to various models.

In the recent papers, we introduced a method utilised to measure the flow field. The method is based on the tracking of supergranular structures. We did not precisely know, whether its results represent the flow field in the photosphere or in some sub-photospheric layers. In this paper, in combination with helioseismic data, we are able to estimate the depths in the solar convection envelope, where the detected large-scale flow field is well represented by the surface measurements. We got a clear answer to question what kind of structures we track in full-disc Dopplergrams. It seems that in the quiet Sun regions the supergranular structures are tracked, while in the regions with the magnetic field the structures of the magnetic field are dominant. This observation seems obvious, because the nature of Doppler structures is different in the magnetic regions and in the quiet Sun. We show that the large-scale flow detected by our method represents the motion of plasma in layers down to ~10 Mm. The supergranules may therefore be treated as the objects carried by the underlying large-scale velocity field.