Travel-time Sensitivity Kernels Research Articles

Testing the volume integrals of travel-time sensitivity kernels for flows

Context. Helioseismic inversions largely rely on sensitivity kernels, in which 3D spatial functions describe how the changes in the solar interior translate into the change in helioseismic observables. These sensitivity kernels in most cases come from the forward modelling that is used in the most advanced solar models. Aims. We aim to test the sensitivity kernels by comparing their volume integrals with measured values from helioseismic travel times. Methods. By manipulating the tracking rate, we mimicked the additional zonal velocity in the Dopplergram datacubes. These datacubes were then processed by a standard travel-time measurements pipeline. We investigated the dependence of the east-west travel time averaged over a box around the disc centre on the implanted tracking velocity. The slope of this dependence is directly proportional to the total volume integral of the sensitivity kernel that corresponds to the travel-time geometry that is used. Results. The agreement between measurements and models for travel times that are computed with a ridge filtering is very good to acceptable. The dependence we sought to determine indeed resembles a linear function, and its slope agrees with the expected volume integral from the forward-modelled sensitivity kernel. The agreement is poorer for the phase-speed filtered datacubes. The disagreement is particularly strong for the slowest phase speeds (filters td1–td4). For the higher phase speeds, our result indicates that the measured kernel integrals are systematically larger than expected from the forward modelling. We admit our testing procedure may not be appropriate for high phase speeds and higher radial modes.

Time-lapse imaging of seismic scattering property and velocity in the northeastern Japan

Investigating spatio-temporal changes in crustal structure is a key to understanding the dynamic processes in the earth's crust associated with large earthquakes and volcanic activities. The dense seismic network in the northeast part of Japan (Tohoku region) allows us to apply seismic interferometry to continuous records of ambient noise data for studying crustal changes related to the 2008 Mw 6.9 Iwate-Miyagi Nairiku earthquake and the 2011 Mw 9.0 Tohoku-oki earthquake. We track changes in both seismic scattering properties and velocities over five years and spatially localize them based on the coda wave decorrelation and travel time sensitivity kernels calculated from the solution of the radiative transfer equation. The 2008 Iwate-Miyagi Nairiku earthquake caused significant changes in both seismic scattering properties and seismic velocities in the vicinity of the epicenter subjected to dynamic or/and static stress perturbations. Changes in scattering properties associated with the 2011 Tohoku-Oki earthquake around some active volcanoes in the Tohoku region are likely related to geofluid migrations triggered by strong vibrations during the large earthquake. The decrease in seismic velocity exists mainly along the eastern coastal region and extends to the central Quaternary volcanic chains under the combined dynamic and static stress changes by the earthquake. The waveform correlation is possibly sensitive to the fluid content and might be a useful indicator of fluid migration in the crust. This is the first comprehensive demonstration and comparison of the spatio-temporal changes of seismic scattering properties and velocities at long-term and regional scales under the influence of two large earthquakes.

Decomposition and Properties of the Traveltime Sensitivity Kernels in Transversely Isotropic Elastic Media

This study discusses the properties and decompositions of the traveltime sensitivity kernels in elastic isotopic and transversely isotropic media. Analytical solutions in homogeneous isotropic elastic media are derived and help better understand their properties. One key component is the kernel phase function which determines the frequency components of the sensitivity kernels, the geometrical setting of the Fresnel zones, and the resolving ability for velocity inversion. For the purpose of exploration seismology, the sensitivity kernels for finite-frequency traveltime inversion in transversely isotropic elastic media are expressed in terms of Thomsen’s anisotropic parameters. The sensitivity kernels in anisotropic media can be decomposed into the summation of the isotropic terms constructed by the elastic properties in the symmetry direction and the correction terms due to the presence of anisotropy. Sensitivity kernels computed from different types of acquisitions contain information for different parameters. We compare the kernels computed using finite-difference with those by dynamic ray tracing. We conclude that the full wave numerical method is able to provide high-resolution kernels, which form the basis for 3D finite-frequency traveltime inversion for velocity and anisotropic parameters.

High-Accuracy Location of Microseismic Events in a Strong Inhomogeneous Mining Environment by Optimized Global Full Waveform Inversion

High-accuracy determination of a microseismic (MS) location is the core task in MS monitoring. In this study, a 3D multi-scale grid Green’s function database, depending on recording wavefield frequency band for the target mining area, is pre-generated based on the reciprocity theorem and 3D spectral element method (SEM). Then, a multi-scale global grid search strategy is performed based on this pre-stored Green’s function database, which can be effectively and hierarchically processed by searching for the spatial location. Numerical wavefield modeling by SEM effectively overcomes difficulties in traditional and simplified ray tracing modeling, such as difficult wavefield amplitude and multi-path modeling in 3D focusing and defusing velocity regions. In addition, as a key step for broadband waveform simulation, the source-time function estimated from a new data-driven singular value decomposition averaged fractional derivative based wavelet function (DD-SVD-FD wavelet) was proposed to generate high-precision synthetic waveforms for better fitting observed broadband waveform than those by simple and traditional source-time function. Combining these sophisticated processing procedures, a new robust grid search and waveform inversion-based location (GSWI location) approach is integrated. In the synthetic test, we discuss and demonstrate the importance of 3D velocity model accuracy to waveform inversion-based location results for a practical MS monitoring configuration. Furthermore, the average location error of the 3D GSWI location for eight real blasting events is only 15.0 m, which is smaller than error from 3D ray tracing-based location (26.2 m) under the same velocity model. These synthetic and field application investigations prove the crucial role of 3D velocity model, finite-frequency travel-time sensitivity kernel characteristics and accurate numerical 3D broadband wavefield modeling for successful MS location in a strong heterogeneous velocity model that are induced by the presence of ore body, host rocks, complex tunnels, and large excavations.

A General Formulation for Computing Spherical Helioseismic Sensitivity Kernels while Incorporating Systematical Effects

Abstract As helioseismology matures and turns into a precision science, modeling finite-frequency, geometric, and systematical effects is becoming increasingly important. Here we introduce a general formulation for treating perturbations of arbitrary tensor rank in spherical geometry using fundamental ideas of quantum mechanics and their extensions in geophysics. We include line-of-sight projections and center-to-limb differences in line formation heights in our analysis. We demonstrate the technique by computing a travel-time sensitivity kernel for sound-speed perturbations. The analysis produces the spherical harmonic coefficients of the sensitivity kernels, which leads to better-posed and computationally efficient inverse problems.

Rytov-approximation-based wave-equation traveltime tomography

Most finite-frequency traveltime tomography methods are based on the Born approximation, which requires that the scale of the velocity heterogeneity and the magnitude of the velocity perturbation should be small enough to satisfy the Born approximation. On the contrary, the Rytov approximation works well for large-scale velocity heterogeneity. Typically, the Rytov-approximation-based finite-frequency traveltime sensitivity kernel (Rytov-FFTSK) can be obtained by integrating the phase-delay sensitivity kernels with a normalized weighting function, in which the calculation of sensitivity kernels requires the numerical solution of Green’s function. However, solving the Green’s function explicitly is quite computationally demanding, especially for 3D problems. To avoid explicit calculation of the Green’s function, we show that the Rytov-FFTSK can be obtained by crosscorrelating a forward-propagated incident wavefield and reverse-propagated adjoint wavefield in the time domain. In addition, we find that the action of the Rytov-FFTSK on a model-space vector, e.g., the product of the sensitivity kernel and a vector, can be computed by calculating the inner product of two time-domain fields. Consequently, the Hessian-vector product can be computed in a matrix-free fashion (i.e., first calculate the product of the sensitivity kernel and a model-space vector and then calculate the product of the transposed sensitivity kernel and a data-space vector), without forming the Hessian matrix or the sensitivity kernels explicitly. We solve the traveltime inverse problem with the Gauss-Newton method, in which the Gauss-Newton equation is approximately solved by the conjugate gradient using our matrix-free Hessian-vector product method. An example with a perfect acquisition geometry found that our Rytov-approximation-based traveltime inversion method can produce a high-quality inversion result with a very fast convergence rate. An overthrust synthetic data test demonstrates that large- to intermediate-scale model perturbations can be recovered by diving waves if long-offset acquisition is available.

Finite-frequency traveltime tomography using the Generalized Rytov approximation

SUMMARY Wave-equation-based traveltime tomography has been extensively applied in both global tomography and seismic exploration. Typically, the traveltime Fréchet derivative is obtained using the first-order Born approximation, which is only satisfied for weak velocity perturbations and small phase shifts (i.e. the weak-scattering assumption). Although the small phase-shift restriction can be handled with the Rytov approximation, the weak velocity-perturbation assumption is still a major limitation. The recently developed generalized Rytov approximation (GRA) method can achieve an improved phase accuracy of the forward-scattered wavefield, in the presence of large-scale and strong velocity perturbations. In this paper, we combine GRA with the classical finite-frequency theory and propose a GRA-based traveltime sensitivity kernel (GRA-TSK), which overcomes the weak-scattering limitation of the conventional finite-frequency methods. Numerical examples demonstrate that the accumulated time delay of forward-scattered waves caused by large-scale smooth perturbations can be correctly handled by the GRA-TSK, regardless of the magnitude of the velocity perturbations. Then, we apply the new sensitivity kernel to solve the traveltime inverse problem, and we propose a matrix-free Gauss–Newton method that has a faster convergence rate compared with the gradient-based method. Numerical tests show that, compared with the conventional adjoint traveltime tomography, the proposed GRA-based traveltime tomography can obtain a more accurate model with a faster convergence rate, making it more suited for recovering the large-intermediate scale of the velocity model, even for strong-perturbation and complex subsurface structures.

Combined helioseismic inversions for 3D vector flows and sound-speed perturbations

Context. Time–distance helioseismology is the method of the study of the propagation of waves through the solar interior via the travel times of those waves. The travel times of wave packets contain information about the conditions in the interior integrated along the propagation path of the wave. The travel times are sensitive to perturbations of a variety of quantities. The usual task is to invert for the vector of plasma flows or the sound–speed perturbations separately. The separate inversions may be polluted by systematic bias, for instance, originating in the leakage of vector flows into the sound–speed perturbations and vice versa (called a cross-talk). Information about the cross-talk is necessary for a proper interpretation of results. Aims. We introduce an improved methodology of the time-distance helioseismology which allows us to invert for a full 3D vector of plasma flows and the sound–speed perturbations at once. Using this methodology one can also derive the mean value of the vertical component of plasma flows and the cross-talk between the plasma flows and the sound–speed perturbations. Methods. We used the Subtractive Optimally Localised Averaging method with a minimisation of the cross-talk as a tool for inverse modelling. In the forward model, we use Born approximation travel-time sensitivity kernels with the Model S as a background. The methodology was validated using forward-modelled travel times with both mean and difference point-to-annulus averaging geometries applied to a snapshot of fully self-consistent simulation of the convection. Results. We tested the methodology on synthetic data. We demonstrate that we are able to recover flows and sound–speed perturbations in the near-surface layers. We have taken the advantage of the sensitivity of our methodology to entire vertical velocity, and not only to its variations as in other available methodologies. The cross-talk from both the vertical flow component and the sound–speed perturbation has only a negligible effect for inversions for the horizontal flow components. Furthermore, this cross-talk can be minimised if needed. The inversions for the vertical component of the vector flows or for the sound–speed perturbations are affected by the cross-talk from the horizontal components, which needs to be minimised in order to provide valid results. It seems that there is a nearly constant cross-talk between the vertical component of the vector flows and the sound–speed perturbations.

Sensitivity kernels for time-distance helioseismology

Context. The interpretation of helioseismic measurements, such as wave travel-time, is based on the computation of kernels that give the sensitivity of the measurements to localized changes in the solar interior. These kernels are computed using the ray or the Born approximation. The Born approximation is preferable as it takes finite-wavelength effects into account, although it can be computationally expensive. Aims. We propose a fast algorithm to compute travel-time sensitivity kernels under the assumption that the background solar medium is spherically symmetric. Methods. Kernels are typically expressed as products of Green’s functions that depend upon depth, latitude, and longitude. Here, we compute the spherical harmonic decomposition of the kernels and show that the integrals in latitude and longitude can be performed analytically. In particular, the integrals of the product of three associated Legendre polynomials can be computed. Results. The computations are fast and accurate and only require the knowledge of the Green’s function where the source is at the pole. The computation time is reduced by two orders of magnitude compared to other recent computational frameworks. Conclusions. This new method allows flexible and computationally efficient calculations of a large number of kernels, required in addressing key helioseismic problems. For example, the computation of all the kernels required for meridional flow inversion takes less than two hours on 100 cores.

Calculating a sensitivity kernel for the final cutoff of received energy in a broadband channel impulse response

The final cutoff of a deep-water ocean acoustic tomography reception is often interpreted through a modal analysis. An alternative approach is presented here. First, a stable travel-time is generated by using an estimator-correlator to account for the statistical variability of the final low-angle arrival with the assumption of an appropriate coherent bandwidth and depth coherence. A time-series of these travel-times can then be tracked using the Viterbi algorithm. Then, the resolvable, trackable peaks are interpreted using the travel-time sensitivity kernel (TSK). The TSK has been previously used to account for the full-wave effects of ray-like arrivals. But the TSK machinery can be adapted to any processing strategy with the application of the chain-rule. The question of how to interpret an arrival produced by an Estimator-Correlator vs. an arrival produced by standard processing is examined. The feasibility of this approach is first analytically determined with the Pekeris waveguide and then an example...

Computational helioseismology in the frequency domain: acoustic waves in axisymmetric solar models with flows

Local helioseismology has so far relied on semi-analytical methods to compute the spatial sensitivity of wave travel times to perturbations in the solar interior. These methods are cumbersome and lack flexibility. Here we propose a convenient framework for numerically solving the forward problem of time-distance helioseismology in the frequency domain. The fundamental quantity to be computed is the cross-covariance of the seismic wavefield. We choose sources of wave excitation that enable us to relate the cross-covariance of the oscillations to the Green's function in a straightforward manner. We illustrate the method by considering the 3D acoustic wave equation in an axisymmetric reference solar model, ignoring the effects of gravity on the waves. The symmetry of the background model around the rotation axis implies that the Green's function can be written as a sum of longitudinal Fourier modes, leading to a set of independent 2D problems. We use a high-order finite-element method to solve the 2D wave equation in frequency space. The computation is `embarrassingly parallel', with each frequency and each azimuthal order solved independently on a computer cluster. We compute travel-time sensitivity kernels in spherical geometry for flows, sound speed, and density perturbations under the first Born approximation. Convergence tests show that travel times can be computed with a numerical precision better than one millisecond, as required by the most precise travel-time measurements. The method presented here is computationally efficient and will be used to interpret travel-time measurements in order to infer, e.g., the large-scale meridional flow in the solar convection zone. It allows the implementation of (full-waveform) iterative inversions, whereby the axisymmetric background model is updated at each iteration.

Resolution, identification, and stability of broadband acoustic arrivals in Fram Strait.

An ocean acoustic tomography system consisting of three moorings with low frequency, broadband transceivers and a moored receiver located approximately in the center of the triangle formed by the transceivers was installed in the central, deep-water part of Fram Strait during 2010-2012. Comparisons of the acoustic receptions with predictions based on hydrographic sections show that the oceanographic conditions in Fram Strait result in complex arrival patterns in which it is difficult to resolve and identify individual arrivals. In addition, the early arrivals are unstable, with the arrival structures changing significantly over time. The stability parameter α suggests that the instability is likely not due to small-scale variability, but rather points toward strong mesoscale variability in the presence of a relatively weak sound channel as being largely responsible. The estimator-correlator [Dzieciuch, J. Acoust. Soc. Am. 136, 2512-2522 (2014)] is shown to provide an objective formalism for generating travel-time series given the complex propagation conditions. Because travel times obtained from the estimator-correlator are not associated with resolved, identified ray arrivals, inverse methods are needed that do not use sampling kernels constructed from geometric ray paths. One possible approach would be to use travel-time sensitivity kernels constructed for the estimator-correlator outputs.

Analysis of the traveltime sensitivity kernels for an acoustic transversely isotropic medium with a vertical axis of symmetry

ABSTRACTIn anisotropic media, several parameters govern the propagation of the compressional waves. To correctly invert surface recorded seismic data in anisotropic media, a multi‐parameter inversion is required. However, a tradeoff between parameters exists because several models can explain the same dataset. To understand these tradeoffs, diffraction/reflection and transmission‐type sensitivity‐kernels analyses are carried out. Such analyses can help us to choose the appropriate parameterization for inversion. In tomography, the sensitivity kernels represent the effect of a parameter along the wave path between a source and a receiver. At a given illumination angle, similarities between sensitivity kernels highlight the tradeoff between the parameters. To discuss the parameterization choice in the context of finite‐frequency tomography, we compute the sensitivity kernels of the instantaneous traveltimes derived from the seismic data traces. We consider the transmission case with no encounter of an interface between a source and a receiver; with surface seismic data, this corresponds to a diving wave path. We also consider the diffraction/reflection case when the wave path is formed by two parts: one from the source to a sub‐surface point and the other from the sub‐surface point to the receiver. We illustrate the different parameter sensitivities for an acoustic transversely isotropic medium with a vertical axis of symmetry. The sensitivity kernels depend on the parameterization choice. By comparing different parameterizations, we explain why the parameterization with the normal moveout velocity, the anellipitic parameter η, and the δ parameter is attractive when we invert diving and reflected events recorded in an active surface seismic experiment.

Interpretation of Helioseismic Travel Times

Time-distance helioseismology uses cross-covariances of wave motions on the solar surface to determine the travel times of wave packets moving from one surface location to another. We review the methodology to interpret travel-time measurements in terms of small, localized perturbations to a horizontally homogeneous reference solar model. Using the first Born approximation, we derive and compute 3D travel-time sensitivity (Fr\'echet) kernels for perturbations in sound-speed, density, pressure, and vector flows. While kernels for sound speed and flows had been computed previously, here we extend the calculation to kernels for density and pressure, hence providing a complete description of the effects of solar dynamics and structure on travel times. We treat three thermodynamic quantities as independent and do not assume hydrostatic equilibrium. We present a convenient approach to computing damped Green's functions using a normal-mode summation. The Green's function must be computed on a wavenumber grid that has sufficient resolution to resolve the longest lived modes. The typical kernel calculations used in this paper are computer intensive and require on the order of 600 CPU hours per kernel. Kernels are validated by computing the travel-time perturbation that results from horizontally-invariant perturbations using two independent approaches. At fixed sound-speed, the density and pressure kernels are approximately related through a negative multiplicative factor, therefore implying that perturbations in density and pressure are difficult to disentangle. Mean travel-times are not only sensitive to sound-speed, density and pressure perturbations, but also to flows, especially vertical flows. Accurate sensitivity kernels are needed to interpret complex flow patterns such as convection.

Wave-equation-based travel-time seismic tomography – Part 1: Method

Abstract. In this paper, we propose a wave-equation-based travel-time seismic tomography method with a detailed description of its step-by-step process. First, a linear relationship between the travel-time residual Δt = Tobs–Tsyn and the relative velocity perturbation δ c(x)/c(x) connected by a finite-frequency travel-time sensitivity kernel K(x) is theoretically derived using the adjoint method. To accurately calculate the travel-time residual Δt, two automatic arrival-time picking techniques including the envelop energy ratio method and the combined ray and cross-correlation method are then developed to compute the arrival times Tsyn for synthetic seismograms. The arrival times Tobs of observed seismograms are usually determined by manual hand picking in real applications. Travel-time sensitivity kernel K(x) is constructed by convolving a~forward wavefield u(t,x) with an adjoint wavefield q(t,x). The calculations of synthetic seismograms and sensitivity kernels rely on forward modeling. To make it computationally feasible for tomographic problems involving a large number of seismic records, the forward problem is solved in the two-dimensional (2-D) vertical plane passing through the source and the receiver by a high-order central difference method. The final model is parameterized on 3-D regular grid (inversion) nodes with variable spacings, while model values on each 2-D forward modeling node are linearly interpolated by the values at its eight surrounding 3-D inversion grid nodes. Finally, the tomographic inverse problem is formulated as a regularized optimization problem, which can be iteratively solved by either the LSQR solver or a~nonlinear conjugate-gradient method. To provide some insights into future 3-D tomographic inversions, Fréchet kernels for different seismic phases are also demonstrated in this study.

Structure and stability of wave-theoretic kernels in the ocean

Wave-theoretic modeling can be applied to obtain travel-time sensitivity kernels (TSKs) representing the amount ray travel times are affected by sound-speed variations anywhere in the medium. This work explores the spatial frequency content of the TSK compared to expected ocean variability. It also examines the stability of the TSK in environments that produce strong sensitivity of ray paths to initial conditions. The conclusion is that the linear TSK model is an effective predictor of travel-time changes and that the rays perform nearly as well as the full-wave kernel. The TSK is examined in physical space and in wavenumber space, and it is found that this is the key to understanding how the travel time reacts to ocean perturbations. There are minimum vertical and horizontal length scales of ocean perturbations that are required for the travel time to be affected. The result is that the correspondence between true travel times and those calculated from the kernel is high for large-scale perturbations and somewhat less for the small scales. This demonstrates the validity of ray-based inversion of travel time observations for the cases under study.

Long-range asymptotic behavior of vertical travel-time sensitivity kernels

Vertical travel-time sensitivity kernels (VTSKs) describe the effect of horizontally uniform sound-speed changes on travel times in range-independent ocean environments. Wave-theoretic VTSKs can be obtained either analytically, through perturbation of the normal-mode representation, or numerically, as horizontal marginals of the corresponding two-dimensional and three-dimensional travel-time sensitivity kernels. In previous works, it has been observed that wave-theoretic finite-frequency VTSKs approach the corresponding ray-theoretic sensitivity kernels as the propagation range increases. The present work is an attempt to explain this behavior. A stationary-phase approach is used to obtain a long-range asymptotic expression for the wave-theoretic VTSKs. The resulting asymptotic VTSKs are very close to the corresponding ray-theoretic ones. The smoothness condition, required for the stationary-phase approximation to hold, is used to obtain an estimate for the range beyond which the asymptotic behavior sets in.

Wave equation tomographic velocity inversion method based on the Born/Rytov approximation

This paper discusses Born/Rytov approximation tomographic velocity inversion methods constrained by the Fresnel zone. Calculations of the sensitivity kernel function and traveltime residuals are critical in tomographic velocity inversion. Based on the Born/Rytov approximation of the frequency-domain wave equation, we derive the traveltime sensitivity kernels of the wave equation on the band-limited wave field and simultaneously obtain the traveltime residuals based on the Rytov approximation. In contrast to single-ray tomography, the modified velocity inversion method improves the inversion stability. Tests of the near-surface velocity model and field data prove that the proposed method has higher accuracy and Computational efficiency than ray theory tomography and full waveform inversion methods.

Asymptotic behavior of vertical travel-time sensitivity kernels with increasing propagation range

Vertical travel-time sensitivity kernels (VTSKs) describe the effect of horizontally uniform sound-speed changes on travel-times in range-independent environments. Wave-theoretic VTSKs can be obtained either analytically, through perturbation of the normal-mode representation, or numerically, as horizontal marginals of the corresponding 2D and 3D travel-time sensitivity kernels. In previous works it has been observed that, as the propagation range increases, wave-theoretic VTSKs approach the corresponding ray-theoretic sensitivity kernels even for low frequencies. In the present work an asymptotic expression of the wave-theoretic finite-frequency VTSKs is obtained, using a stationary-phase approach. Numerical results show that wave-theoretic VTSKs converge with increasing range toward the asymptotic form, which in turn lies very close to the ray-theoretic VTSK. [Work supported by ONR.]

Validated helioseismic inversions for 3D vector flows

According to time-distance helioseismology, information about internal fluid motions is encoded in the travel times of solar waves. The inverse problem consists of inferring 3-D vector flows from a set of travel-time measurements. Here we investigate the potential of time-distance helioseismology to infer 3-D convective velocities in the near-surface layers of the Sun. We developed a new Subtractive Optimally Localised Averaging (SOLA) code suitable for pipeline pseudo-automatic processing. Compared to its predecessor, the code was improved by accounting for additional constraints in order to get the right answer within a given noise level. The main aim of this study is to validate results obtained by our inversion code. We simulate travel-time maps using a snapshot from a numerical simulation of solar convective flows, realistic Born travel-time sensitivity kernels, and a realistic model of travel-time noise. These synthetic travel times are inverted for flows and the results compared with the known input flow field. Additional constraints are implemented in the inversion: cross-talk minimization between flow components and spatial localization of inversion coefficients. Using modes f, p1 through p4, we show that horizontal convective flow velocities can be inferred without bias, at a signal-to-noise ratio greater than one in the top 3.5 Mm, provided that observations span at least four days. The vertical component of velocity (v_z), if it were to be weak, is more difficult to infer and is seriously affected by cross-talk from horizontal velocity components. We emphasise that this cross-talk must be explicitly minimised in order to retrieve v_z in the top 1 Mm. We also show that statistical averaging over many different areas of the Sun allows for reliably measuring of average properties of all three flow components in the top 5.5 Mm of the convection zone.