One-sided arc averaging geometries in time-distance local helioseismology

The study of solar oscillations (helioseismology) has been a very successful method of researching the Sun. Helioseismology teaches us about the structure and mean properties of the Sun. Together with mid-resolution data, the local properties were uncovered in quiet-Sun regions. However, magnetic fields affect the oscillations and prevent us from studying the properties of magnetically active regions with helioseismology. We aim to create a new methodology to suppress the negative effects of magnetic fields on solar oscillations and measure plasma properties close to active regions. The methodology consists of new averaging geometries, a non-linear approach to travel-time measurements, and a consistent inversion method that combines plasma flows and sound-speed perturbations. We constructed the one-sided arc averaging geometries and applied them to the non-linear approach of travel-time measurements. Using the one-sided arc travel times, we reconstructed the annulus travel times in a quiet-Sun region. We tested the methodology against the validated helioseismic inversion pipeline. We applied the new methodology for an inversion for surface horizontal flows in a region with a circular H-type sunspot. The inverted surface horizontal flows are comparable with the output of the coherent structure tracking, which is not strongly affected by the presence of the magnetic field. We show that the new methodology suppresses the negative effects of magnetic fields up to outer penumbra. We measure divergent flows with properties comparable to the moat flow. The new methodology can teach us about the depth structure of active regions and physical conditions that contribute to the evolution of the active regions.


Introduction
Helioseismology studies the solar interior via properties of the solar oscillations caused by the sound and the surface gravity waves that propagate through the Sun. The waves are affected by plasma conditions along their trajectory. Typical methods of local helioseismology utilised for inferring the plasma properties in mid resolution are the ring-diagram method (Hill 1988) and the time-distance method (Duvall et al. 1993). The ring-diagram method is based on the direct study of an average power spectrum and its resonant frequencies. The resonant frequencies are also sensitive to the plasma flows and the sound speed via the Doppler effect. Time-distance helioseismology estimates the properties of the plasma flows and the sound-speed perturbations via the travel times of the waves between two points on the solar surface. The travel times are naturally affected by the plasma flows and sound speed. These methods have been utilised in many studies (e.g. Basu et al. 1999;Kosovichev et al. 2000;Brüggen & Spruit 2000;Haber et al. Send offprint requests to: David Korda, e-mail: korda@sirrah.troja.mff.cuni.cz A&A proofs: manuscript no. main 2014; Braun 2019). To avoid issues with the magnetic fields, researchers have tried to mask out the active regions (e.g. Liang & Chou 2015). Zhao & Kosovichev (2003) inverted for plasma flows and sound-speed perturbations beneath a sunspot. They have also compared results with and without a maskedout active region but did not find significant differences between the corresponding results. Korzennik (2006) inverted for a set of models of sound-speed perturbations. For each model, he masked out a circular region around a sunspot with different radius and found out that 'plume' sound-speed perturbations below an active region are probably caused by the surface contamination.
A different strategy was chosen by Duvall et al. (2018). They observed points far from a sunspot and selected waves with two skips, where the first skip is inside the sunspot. Hughes et al. (2005) compared the inverted sound-speed perturbations using single-skip and double-skip geometries. They found a good agreement at depths larger than about 5 Mm below a sunspot. The waves with two skips cannot be separated from those with one skip if the observed points were closer to each other. For this reason, they distinguished a lack of waves propagating through the shallow layers in the double-skip geometry. For this reason, this type of measurement is not suitable in shallow layers.
Taking the arguments given above into account, the subsurface structure and the dynamics of the active regions are still unknown. This is in contrast with the fact that the large active regions are potentially dangerous for a technical civilisation via solar flares.

Motivation
Around an evolved sunspot, divergent horizontal flows called moat flows are observed. The amplitude of a moat flow is about 500 m s −1 (see e.g. Švanda et al. (2014)). Numerical simulations (Rempel et al. 2009) predicted a uniform outflow from the centre of the sunspot. The helioseismic inversions (e.g. Zhao et al. 2010) resulted in a multi-layer structure, in which the direction of the flows and the sign of the sound-speed perturbations changed repeatedly. Moreover, helioseismic results depend on the specific approach (Gizon et al. 2009). Therefore, the depth profile of the moat flow is unknown unless it has a large amplitude.
In the presence of the magnetic field, the conversion between sound waves and magnetoacoustic waves appears. The magnetoacoustic waves are not trapped inside the Sun and escape into the atmosphere (Schunker et al. 2013). In this case, we can observe the wave at only one of the points between which we measure the travel times. Therefore, the amplitudes of the detected waves are lower in active regions. This effect is usually misinterpreted as plasma flows or sound-speed perturbations. The amplitude of the detected waves (the amplitude of the cross-correlation function C defined in Sect. 3) explicitly appears in the linear approach of the travel-time measurements (GB04; Gizon & Birch 2004), while it does not enter non-linear approach (GB02; Gizon & Birch 2002). Therefore, the GB02 travel times are more suitable for studies of the active regions. The comparison between GB02 and GB04 travel times in an active region was done by DeGrave et al. (2014) who showed the advantages of the GB02 approach in active regions.
Time-distance helioseismology uses the averaging geometries. These increase the signal-to-noise ratio. The most often used averaging geometries assume quasi-uniform properties of the observed data since the cross-correlation is computed between the observed data at a point and the data averaged over Examples of the annulus (the outer annulus) and the arc (the inner arc) averaging geometries. The black dot is the point where the travel time is measured. The colours correspond to the weights of the averaging geometries. The black circle delimits the artificial active region. The annulus geometry averages the quiet and the active regions because of its two-sided character; the observed point in the one-sided arc geometry is not affected by the active region. For display purposes, the radii of the averaging geometries differ, the full width at half maxima (FWHMs) of the Gaussian annuli were enlarged, and the weights of the arc geometry were divided by 12.
the annulus or quadrant around the point (e.g. Duvall et al. 1997;Kosovichev et al. 2000;Švanda 2013). Then, the final travel times are computed from a combination of the individual traveltime measurements. For example, in quadrant geometry the travel times sensitive to flows in the longitudinal direction are computed as the difference between the travel times measured in and against the direction of the solar rotation. This effectively makes the final averaging geometries (and the travel times) sensitive to the conditions all around the observed point. Therefore, these averaging geometries are very useful in quiet-Sun regions, but not in the vicinity of active regions, because they average the active and the quiet regions together. An effect of an active region to the cross-correlation function was demonstrated, for example, by Cameron et al. (2008). They computed the crosscorrelation function and showed that an active region absorbs non-magnetic waves (see their Fig. 4). Similar results were obtained by Gizon et al. (2009), who measured the absorption of surface gravity waves and sound waves with different radial orders. The absorption coefficients were between 49% and 57%. Moreover, the two-sided character of the averaging geometries causes omnidirectional disturbances even in quiet-Sun regions close to an active region, as was pointed out by Hughes et al. (2005).
As mentioned in the introduction, the possible way to correctly measure plasma parameters below a sunspot is to utilise the double-skip geometry. As a trade-off, the user cannot reliably invert for the parameters in the first few megametres (Hughes et al. 2005). The way to avoid these issues is to abandon the approaches of the double-skip geometry and the two-sided averaging geometries.
In this paper, we introduce a new type of averaging geometry: one-sided arc geometries. As seen in Fig. 1, the one-sided geometries do not average the quiet and the active regions together. Naturally, the one-sided averaging geometries do not need to use the double-skip geometry. Therefore, the one-sided averaging geometries can be effective in the vicinity of the active regions, where the moat flow appears. In Fig. 1, we plotted examples of the annulus (the outer annulus) and the arc (the inner arc) averaging geometries. In this case, both geometries are sensitive to flows in the x direction. The black dot at the coordinates [0, 0] is the point at which the travel time is measured, the colours indicate the weights of the averaging geometries, and the black circle corresponds to the edge of an artificial active region. In the case of the annulus geometries, the central point is affected by the active region, even though the point itself is not inside the active region (part of the annulus is inside it). This is not the case of the one-sided arc geometries. Therefore, we can obtain unspoiled information about a much closer vicinity of the active region if we combine the one-sided arc geometries and the non-linear GB02 travel-time approach.

Travel times
The travel times τ of the sound and the surface gravity waves, used in the time-distance helioseismology, are measured between the two points r 1 and r 2 on the solar surface. The travel time τ (r 1 , r 2 ) is defined as the time lag that maximises the temporal cross-correlation C (r 1 , r 2 , t) of the observed signal ψ (r, t) at the points r 1 and r 2 (Gizon & Birch 2002): where T is the duration of the observation. In the time-distance helioseismology, authors often utilised the perturbed travel times δτ (r 1 , r 2 ) , defined as the difference between the measured travel time τ msm and the modelled travel time τ model : δτ (r 1 , r 2 ) ≡ τ msm (r 1 , r 2 ) − τ model (r 1 , r 2 ) .
3.1. Point-to-annulus averaging geometries The perturbed travel times are usually measured under the given averaging geometries, which reduce the realisation noise component. In applications widely used in literature, four pointto-annulus (PtA) averaging geometries are distinguished. The temporal cross-correlation is always computed between the observed signal ψ at the given point and the observed signal averaged over the normalised annulus with the radius ∆. The final perturbed travel times are constructed from the combinations of the travel times measured from the central point to the rim of the annulus and the travel times measured in the opposite direction. If the annulus is symmetrical with respect to its centre, we refer to the outflow-inflow (o-i; the difference between δτ (r = r 1 , ∆ = r 1 − r 2 ) and δτ (∆, r)) and the mean (mn; the mean value of the two travel times) averaging geometries. The annulus values can be weighted with the cosine (east-west geometry; e-w) or the sine (north-south geometry; n-s) of the polar angle. These geometries are sensitive to the direction of the wave trajectory. Only the differences between the two travel times are computed for the e-w and the n-s geometries.

Point-to-arc averaging geometries
The PtA averaging geometries are very useful in quiet-Sun regions because the statistical properties of waves with the same phase speed and generated at the same point are comparable regardless of the direction of propagation. If this is not true, the PtA averaging geometries average regions that are not similar. Hence, the PtA geometries do not fit the real physical conditions in such regions well. For this reason, we introduce new averaging geometries: one-sided point-to-arc (PtArc) averaging geometries. They represent a natural transition between noisy one-sided point-to-point travel-time geometries and omnidirectional PtA travel-time geometries with a reasonable signal-to-noise ratio. These geometries do not average the observed signal along all the angles but only over small arcs. Therefore, the travel time measured at the edge of a region with different properties is correctly split into travel times inside and outside this region (see Fig. 1). With the PtArc averaging geometries we can, in principle, invert for plasma properties around active regions as a function of depth, which is highly important. Theoretically, the ideal solution is to use very narrow arcs, or, equivalently, point-to-point travel times defined by Eq. (3); however, the user must deal with increasing contribution from the realisation noise. Additionally, many independent observations are needed to obtain the same quantity of information, and the inversion process is slowed down.

Construction
One particular arc geometry is defined by three numbers; the angular width of the arc α, the orientation of the arc axis ϕ, and the distance ∆ (simplified draw is plotted in Fig. 2). The geometry is constructed in the following way: 1. select the three numbers, 2. construct the unit Gaussian annulus with the given FWHM at the distance ∆ from the central point, 3. construct a function that is equal to one inside the angle interval ϕ − α/2; ϕ + α/2 and zero elsewhere, 4. weight the function with cos (θπ/α), where θ is the inclination of a point of the function with respect to the reference direction defined by ϕ, 5. multiply the Gaussian annulus and the weighted function, and then 6. the arc geometry is attained via normalisation of the previous multiplication to unity.
For the given α and ∆, we constructed a set of arcs directing to a set of polar angles ϕ, usually with a critical sampling given the α.
The construction of the arc geometry defined by α = π/4, ϕ = 0, and ∆ ≈ 22 Mm is shown in Fig. 3. We note that the annulus geometries can be reconstructed from the arc geometries only if α = π; for instance, the annulus geometry plotted in Fig. 1 corresponds to the combination of two arc geometries defined by α 1 = π, ϕ 1 = 0, ∆ 1 ≈ 29 Mm and α 2 = π, ϕ 2 = π, ∆ 2 ≈ 29 Mm. We tested the methodology with the following set of numbers: α = π/4, ϕ from 0 to 2π − α with the step α (eight geometries overall), the FWHM of the Gaussian annulus of approximately 1.46 Mm, and ∆ from about 7.3 Mm to about 29.2 Mm with the step of about 1.46 Mm. The selected angular width and orientations covered all directions. The step α corresponded to the smallest value which enabled to enforce the given symmetries of a reference model without an interpolation.

Forward and inverse models
Using an approximate solution of the wave equation in a perturbed environment we are able to compute the sensitivity kernels K. The sensitivity kernels quantify the changes in travel times under the given changes in the reference model. In general, the sensitivity kernels differ by reference model and with a level of approximation. We utilised the model S (Christensen-Dalsgaard et al. 1996) as the reference model and the Born approximation (Gizon & Birch 2002). The flow and sound-speed sensitivity kernels were computed according to Birch & Gizon (2007) and Birch et al. (2004).
The perturbations of travel times δτ can be computed as where r and r ′ are the horizontal positions, z the vertical position, the upper index a denotes the independent measurements, P is the number of perturbers, and n a the realisation noise.
The inverse methods deconvolve Eq. (4) and teach us about the perturbers δq β . Due to the noise term, the deconvolution leads to the approximate solution δq inv α . The multi-channel optimally localised averaging method (MC-SOLA; Jackiewicz et al. 2012) applies an approximate inverse method and searches for the approximate solution in the form where r i and r 0 are the horizontal positions, z 0 the target depth, N the number of horizontal positions, and w α a the unknown weight functions (inverse filters) that must be determined. In our inverse models, the weight functions minimise the χ 2 functional in the form: where µ, ν and ǫ are the trade-off parameters, K α β the averaging kernel, T α β the target function, Λ ab the noise covariance matrices, λ β the Lagrange multipliers, and δ α β the Kronecker delta function. We invite the reader to consult  for more details about the implementation.

Consistency check
The validity regions for both PtA and PtArc travel-time measurements overlap in quiet-Sun regions. We utilised the aforementioned inverse method and perform consistency checks. We tested the consistency of the new one-sided arc geometries in two ways. First, we approximated the PtA travel times from the specific combination of the PtArc travel times. This reconstruction was a consistency check of the new PtArc travel times. Second, we compared the inversions for horizontal flows in a quiet-Sun region based on PtA and PtArc travel times. In quiet-Sun regions, the results should be very similar, hence, this was a check of the whole PtArc pipeline. We note that if the travel times were consistent but the pipeline was not, the issue must have been in the sensitivity kernels.

Travel times
Due to our selected α π, the PtA travel times cannot be fully reconstructed from the PtArc travel times. For the reconstruction, we utilised the Dopplergrams measured with SDO/HMI Schou et al. 2012) on 8 January 2011. The Dopplergrams were averaged over 24 hours, and the central part of the observed map was at the Carrington coordinates of 156.6 • longitude and 0.0 • latitude (the disc centre was at the coordinates of 156.6 • longitude and −3.8 • latitude). From the Dopplergrams, we measured the travel times in the PtA and the PtArc averaging geometries. In order to compare the travel times, we combined the PtArc travel times in the following way: 1. measure the PtArc travel times from the central point to the rim of the arc and in the opposite direction, 2. for the PtArc travel times with the given ϕ, multiply the travel times with the following: one in the case of o-i and mn geometries; cos (ϕ) in the case of the e-w geometry; sin (ϕ) in the case of the n-s geometry, and sum up the results, 3. normalise the result so that the corresponding averaging geometries have the same integrals for positive and negative parts, 4. subtract the 'opposite-direction' travel times (o-i, e-w, ns) and compute the mean of the 'opposite-direction' travel times (mn).
In Figs. 4 and 5, we show the exact GB04 annulus travel times (top rows) and the reconstructed annulus travel times from the GB02 one-sided arc travel times (bottom rows). We note that we computed the f -mode travel times with ∆ ≈ 22 Mm. The corresponding sub-plots have the same colour bars. In the right panels, there are the corresponding averaging geometries. As one can see, the reconstruction could not be perfect because the averaging geometries were slightly different. Because of the cosine and the sine weighting, only six arc geometries were used in the cases of the e-w and the n-s travel times, while all eight in the cases of the o-i and the mn travel times. The statistical properties of the comparison are summarised in Table 1, namely the correlation coefficients between the corresponding travel times and the root mean square (rms) of the differences between the corresponding travel times.
Histograms of the arc-to-annulus and the corresponding annulus travel times are plotted in Figs. 6 and 7. The bin widths are 5 seconds in all histograms. In the histograms in Fig. 6, the black bars indicate the arc-to-annulus travel times, and the yellow bars represent the annulus travel times. The larger widths of the histograms of the arc-to-annulus travel times are probably caused by lower averaging. The black lines in Fig. 7 have unit slopes. The slopes of the linear regressions are 0.89 (o-i), 0.82 (mn), 0.89 (e-w), and 0.90 (n-s). Except the histograms of the mean travel times, the other histograms are centred close to δτ = 0 s. The histogram of the mean travel times is centred at about δτ = 12 s. This shift is most probably due to the difference of the plasma properties in the near-surface layers of the Sun (roughly the last 1 Mm, where the f -mode is significant) as compared to the model S, which is a reference for our traveltime measurement. This systematic offset does not affect the difference travel times, because it is automatically subtracted. It is important to note that the offset in the mean travel times is the same for both the PtA travel times and the equivalent travel times derived from our PtArc measurements.
In the cases of the o-i and the mn geometries, the reconstructed travel times are very similar to the exact annulus travel times. Not only are the structures of the travel times nearly the same, but the amplitudes are also comparable. The correlation coefficient of the two averaging geometries (right panels in Fig. 4) is 0.90; therefore, the correlation coefficients between the corresponding travel times should not be larger than this value because of the different realisation noise contributions.  The e-w and the n-s averaging geometries and travel times are plotted in Fig. 5. The differences between the corresponding travel times are slightly larger because of the larger differences between the corresponding averaging geometries. However, the

Horizontal flows
In the past, our 'PtA' horizontal flows in quiet-Sun regions were compared with various types of independent methods (see e.g. Švanda et al. 2013;). Therefore, we can validate the new arc geometries via a comparison of the inverted horizontal flows in quiet-Sun regions with our PtA inversions. We performed two sets of surface ( f -mode) inversions for horizontal flows. The first set was based on the PtA averaging geometries, while the other one was based on the new PtArc averaging geometries. Similarly to the approach taken in the previous section, we utilised the Dopplergrams measured at the quiet disc centre on 8 January 2011. The inverted horizontal flows are plotted in Fig. 8. In the top row, we show the inversions for longitudinal flows, and in the bottom row we show the inversions for the latitudinal flows. In the left column, one can see the 'GB04, annulus' horizontal flows, while in the right column the new onesided 'GB02, arc' horizontal flows are documented. Even though a different data averaging and slightly different averaging kernels and realisation noises were used, the results in both rows are nearly the same. The larger amplitude of the PtArc horizontal flows corresponds to the larger amount of longer PtArc travel times (see Fig. 6) and is likely caused by lower averaging of the input data and differences between the averaging kernels. In Table 2, one can find correlation coefficients, rms and mean of the 'arc' minus 'annulus', and the slope s of a linear fit with an arc model as an independent variable (so PtArc × s ≈ PtA).

Test on a circular active region
In the previous section, we show that the new arc travel-time geometries provide us with very similar results to the annulus  geometries in the quiet-Sun regions. The advantage of the arc geometries emerges close to active regions.
For simplicity, we assume a circular sunspot with the radius ρ and the centre at the position r s , the averaging geometry with the distance ∆, and the position of the measurements r m . Then, in the case of the annulus averaging geometries, the affected measurements are inside the disc defined by r s − r m ≤ ρ + ∆. In the case of the one-sided arc geometries, the affected points are inside the ellipse-like figure, whose 'major axis' is in the ϕ direction. The lengths of the 'semi-major axes' differ in the direction of ϕ and ϕ + π. In the ϕ + π direction, its length is ρ + ∆, similarly to the annulus geometries, but in the direction ϕ, the length is ρ. The distance ∆ is often greater than 20 Mm; therefore, the combination of the arc geometries with different ϕ parameters should provide us with the unspoiled information about the plasma properties much closer to the sunspot. The damaged parts of the travel-time measurements are visualised in Fig. 9. In the left panel, we show the affected region for the outflow-inflow annulus averaging geometry, while in the right panel we show the affected region in the case of one of the arc geometries. The yellow central contour corresponds to the edge of an active region. By selecting a set of ϕs, we can 'scan' the close vicinity of the spot. Outflow-inflow annulus geometry. Right: One-sided arc geometry with α = π/4 and ϕ = π. The black contours correspond to umbra (full) and penumbra (dashed). The red wedges show the orientation and the width of the arc geometry.

Travel times
For the demonstration, we selected a roundish H-type active region NOAA 11084 observed with SDO/HMI on 2 July 2010 at Carrington coordinates of 144.5 • longitude and −19.1 • latitude (the Dopplergram centre was at the coordinates 143.0 • longitude and −20.0 • latitude; the disc centre was at 137.7 • longitude and 3.1 • latitude). We measured the surface-mode travel times with ∆ ≈ 29 Mm around this active region in both the annulus and the one-sided arc averaging geometries. Both travel times were measured using the GB02 approach. The measured travel times are plotted in Fig. 10. In the left panel, we show the annulus travel time, and in the right panel the arc travel time (α = π/4, ϕ = π) is shown. The full central contour roughly surrounds the umbra and the outer dashed contour surrounds the penumbra. The red wedges show the orientation and the width of the arc geometry. In accordance with the estimation, the arc travel time is spoiled, simply said, in the region outside the red wedge. Inside the wedge in the penumbra, the amplitude of the arc travel time does not seem to be reduced by the presence of the magnetic field. Therefore, the one-sided arc travel times in this specific geometry with combination of the GB02 travel-time approach contain undisturbed information about plasma flows in the eastern direction.

Horizontal flows
The surface inversions for the horizontal flows can be directly compared with other methods independent on helioseismology. One of the suitable methods is local correlation tracking, which does not suffer from the magnetic fields. The advanced coherent structure tracking (CST; Roudier et al. 1999;Rieutord et al. 2007) is particularly very convenient for the direct comparison of the horizontal flows. These methods are based on tracking granules. Therefore, these methods fail at the edges and within the active regions where the granules are not observed.
To compare the horizontal flows from the CST technique and helioseismic inversions, we need to decrease the spatial resolution of the CST results. The spatial resolution of our helioseismic inversions is quantified by the averaging kernels K α β . Therefore, we convolved the CST horizontal flows with corresponding components of the averaging kernels. The CST provides us with surface flows only; therefore, we only utilised the surface cut through the averaging kernel. The convolution with the averaging kernels has one side effect. The CST method utilised the segmentation of granules and their tracking in the co-aligned time-lapse datacubes. In the umbra and penumbra, there are no granules. The structures segmented by the segmentation (such as penumbral grains) have different velocity properties, not to mention that their velocities are only apparent (see e.g. Schlichenmaier et al. 1998). Hence, the magnitudes of the CST-derived horizontal flows are much smaller (by almost an order of magnitude) within the umbra and penumbra regions and this decrease of the flow magnitude is abrupt on the quiet Sunpenumbra boundary. When smoothing by the averaging kernel, in the moat region the kernel also accumulates the values from the penumbra, which leads to a decrease in the CST-derived flow magnitude.
The results are plotted in Fig. 11. In each panel, one can see the inversion model under the given approach. The arrows visualise the inverted horizontal flows. The green contours surround the umbra and the penumbra. The background images are the corresponding intensitygrams. The agreement between the middle (PtArc) and the right (CST * K) panels in regions where both methods return valid results provides strong evidence that the new one-sided arc averaging geometries can teach us about the properties of active regions. Using PtArc geometries, we detect a divergent flow whose location and amplitude is comparable with the moat flow.
In the left panel, we show the horizontal flows derived from the usual PtA inversion applied to the GB02 travel times. A similar approach was also used, for example, by Švanda et al. (2014). In the quiet-Sun regions close to the active region, one can see the moat flow with the amplitude of about 500 m s −1 . The measured moat flow is about one-third more extended in comparison to both the PtAr inversion model and the CST model(see Fig. 12). Hence, the divergent moat flow is measured in places where both the PtArc and the CST models changed their signs and no longer detect the moat flow. The larger extension of the moat flow in the PtA flow inversion is likely due to the two-sided travel-time averaging geometries.
In this regard, we note that Hughes et al. (2005), for example, also found that the travel times of the waves are affected by the nearby magnetised region in a region much larger than the size of the active region. It is in agreement with our finding that the PtA inversions are affected in a region larger than  Radial profiles of azimuthally averaged speeds. The local minima of the speed are at distances 37.3, 28.1, and 28.5 Mm for PtA, PtArc, and CST respectively. The green vertical lines indicate the edges of umbra and penumbra. The distances correspond to minima of local parabolic fits. We note that to estimate the extent of the moat flow we used raw CST data; otherwise, the edge is systematically shifted by 1 Mm.
the magnetised region itself. The PtArc model, unlike the PtA model, at least provides us with correct information about the extent of the moat flow and perhaps also about the flows up to the edge of penumbra. As a further test, we performed a similar comparison for even simpler flow inversions than PtA presented here in Appendix B. Those inversions are very straightforward in the interpretation. However, they once again involve the difference (two-sided) travel-time geometries and thus lead to the derived moat-like region, which is more extended than the moat flow measured in PtArc or CST models. We note that in the more distant quiet-Sun regions all the discussed models are comparable.
In both travel-time inversions, the magnitude of the moat flow seems about 50% larger than that derived by the CST granulation tracking. We note that the surface feature-tracking methods and time-distance inversions do not represent exactly the Table 3. Correlation coefficients between the CST * K and the 'arc' horizontal flows and rms and mean of the arc minus CST * K. We note that the parameters were computed from the active region and its quiet surrounding plotted in Fig. 11.
same depths, as the time-distance inversions aggregate some signal from sub-surface layers. Some numerical simulations (e.g. Rempel 2011) indicate that there might be an increase of magnitude of the flows around the simulated sunspot just under the surface. The decrease of the CST-derived flow magnitude as a side effect of the resolution matching described above also plays a role here.
The statistical comparison between the PtArc and the CST horizontal flows is in Table 3. We note that larger rms and mean values were mostly due to the active region and slightly different spatial resolution.
In Tables 4 and 5, we summarise the comparison between the helioseismic inversions for surface horizontal flows and the surface CST measurements. The correlation coefficients between the models are listed above the diagonals. The parameter s corresponds to the slopes of the linear fits of the scatter plots between the given models. The model listed in the row was considered to be the independent variable. Due to slightly different smearing of the models, the fits contained a non-negligible constant parameter. We computed two sets of the s parameters. The first set was computed from the area visualised in Fig. 11 after masking out the active region. The slopes are about unity in this case due to most of the area is in quiet Sun. In Tables 4 and 5, these values are in parentheses. The second set was computed from an annulus around the active region. The inner and outer radii of the annulus were about 10 Mm and 37 Mm, respectively. The annulus masked out the imprecise measurements inside the umbra and the inner penumbra and the quiet-Sun region around it. These slopes are outside the parentheses in the two tables. Table 4. Summary table of the comparisons between PtA and PtArc helioseismic inversions and the surface CST * K measurements. Above the diagonal, the correlation coefficients are given. Below the diagonal, the scaling coefficient s corresponds to the slope of a fitted line. The value of s represents the coefficient with which we have to multiply the row in order to obtain the column. The values outside and inside the parentheses were computed from a closer or wider surrounding of the active region, respectively. P P P P P P P s

Conclusions
We introduce the one-sided arc averaging geometries in the timedistance local helioseismology. The one-sided geometries allow for inhomogeneities in the observed signal, which can be useful especially in the vicinity of active regions. The one-sided arc geometries are a superset of the annulus geometries. We selected the specific arc geometries with the width of the arc α = π/4 and the arc orientations ϕ from 0 to 2π − α with the step α. With these arc geometries, we successfully reconstruct the corresponding annulus travel times. The remaining differences are caused by the different realisation noise components. From the histograms and scatter plots of the corresponding travel times, we note the travel-time histograms are narrower (see Fig. 6) in the case of the annulus geometries. This may be related to the higher level of averaging in the case of the annulus geometries.
We observe the trend in the arc travel times. Similar to the mean travel times, the trend is caused by the foreshortening. The radial centre-to-limb trend in the case of the mean travel times is changed to the band-like trend with the same orientation ϕ as the given arc geometry. We reconstruct the radial trend with the combination of the arc geometries, and we propose and test a trend-subtraction method (see Appendix A).
We compared the horizontal flows inverted under the annulus and the new one-sided arc averaging geometries in the quiet-Sun region, where the inversions of the annulus travel times were tested against various independent measurements. Both inverted results are similar (see Fig. 8); therefore, the new one-sided arc geometries give us the correct results in the quiet-Sun regions.
We demonstrate the advantages of the new one-sided arc geometries in active regions. We tested the inverted horizontal flows around the H-type active region NOAA 11084 using the CST technique. Furthermore, we show that by using this methodology we can successfully measure flows up to the outer penumbra, including the moat flow. Again, both sets of horizontal flows (CST-based and PtArc) are comparable (see Fig. 11). Therefore, we suggest applying the one-sided arc averaging geometries and GB02 travel-time approximation in the vicinity of active regions in order to learn about the depth structure of flows around them.

Appendix B: Point-to-quadrant-like inversions
To further demonstrate the improvement of our PtArc inversions over the PtA inversions in the magnetised regions, we constructed very simple and very basic point-to-quadrant-like inversions, essentially following the strategy used in the advent of time-distance A&A proofs: manuscript no. main helioseismology. We limited the set of travel-time geometries to the simplest possible case, essentially following the findings of Burston et al. (2015). According to that work, for instance, the east-west component of the horizontal velocity (v x in our notation) has the strongest impact on e-w travel times (following the notation described in Section 3.1). Similarly, the north-south component of velocity (v y ) has the strongest impact on n-s travel times. These inversions are very simple to interpret. Therefore, we constructed inversions involving only those travel-time geometries. Such inversions are, in principle, equivalent to point-to-quadrant flow inversions used, for example, by Hindman et al. (2004) (their Fig. 6 and related text) and introduced much earlier by . Our approach differs from that of , Hindman et al. (2004), and others, by the fact that our averaging geometry is a smooth function of the azimuthal angle (sine or cosine), whereas they used a stepwise function of the azimuthal angle to separate the east, west, north, and south quadrants. Despite the differences, the principles are similar, and hence we can denote these simplistic inversions as point-to-quadrant (PtQ-1).
For simplicity and to obtain results directly comparable to the inversions shown in Fig. 11, we once again focused only on the f -mode travel times. In fact, our very simplistic PtQ-1 inversions considered a combination of only 16 independent traveltime measurements with varying sizes of the surrounding annulus. This is significantly less than in the case of the PtA inversions presented in Fig.11 (48 travel-time measurements combined) and PtArc (128 independent measurements combined).
The results are shown in Fig. B.1. Due to the lower number of travel-time measurements considered, the signal-to-noise ratio S/N ≡ rms 2 (inverted) − rms 2 (noise)/rms (noise) is slightly worse in PtQ-1 inversions (the S/Ns are 6.4 for PtQ-1, 6.6 for PtA, and 20.0 for PtArc inversions, respectively). The averaging kernel (Fig. B.2) shows that there is a considerable cross-talk from the perpendicular horizontal velocity component (about 25%), of which the realisation depends on a particular configuration of the flow. For pure divergent flows, these cross-talk contributions act as convergent flows. In order to minimise the cross-talk between horizontal flows, we performed another set of PtQ inversion (hereafter referred to as PtQ-2). In these inversions, we combined both horizontal-flow components and e-w and n-w averaging geometries in one inversion setup; therefore, the cross-talk can be suppressed by the technique introduced by Švanda et al. (2011) that is an inseparable part of our inversion pipelines. The averaging kernel for v x component is shown in Fig. B.3. We note that S/N = 7.8 in this case.
The correlation coefficients between PtQ-1 and CST * K in the same format as in Tables 4 and 5 are 0.81 (0.69) for v x and 0.74 (0.68) for v y , and for the PtQ-2 numbers these are 0.87 (0.80) for v x and 0.80 (0.76) for v y . These are lower than correlation coefficients between PtArc and CST * K presented in Tables 4 and 5. The slopes s corresponding to the values in parentheses in Tables 4 and 5 between PtQ-1, PtQ-2, and CST * K are 0.88 for PtQ-1 (both components) and 1.17 (PtQ-2, v x ) and 1.18 (PtQ-2, v y ). All these values are farther from unity than the corresponding values between PtArc and CST * K. The spatial extent of the moat The red and the green contours correspond to the half-maximum of the target function and the averaging kernel, respectively. The full and the dotted blue contours correspond to +5% and −5% of the maximum of the averaging kernel. The red and the green full and dotted contours correspond to the half-maximum of the target function and the half-maximum and the minus half-maximum the of averaging kernel, respectively. The full and the dotted blue contours correspond to +5% and −5% of the maximum of the averaging kernel. We note that the crosstalk contribution from the perpendicular horizontal component (second column) is greatly minimised as compared to Fig. B.2.
still seems larger in the case of the PtQ models as seen in Fig. B.4. The PtQ-1 and the PtQ-2 moat flows are about 25% and 10% more extended in comparison with PtArc and CST, respectively. In the moat region, the PtQ-derived horizontal flow magnitude is about a factor of 2 smaller than the magnitude of the flows inverted using PtA and PtArc approaches. This magnitude decrease is probably due to the affection of the travel times by the presence of the magnetic field nearby. In the travel-time maps, the east-west and north-south travel times are greatly affected in the region of the sunspot, which inevitably leads to the suppression of the derived flow magnitude. On the other hand, the out-in travel time has a very strong magnitude in the sunspot region (see Fig. 10; leaving it proper interpretation aside), which explains why the PtA inversion (which combined not only the directional e-w and n-s but also the omnidirectional o-i geometry) yields a much larger flow amplitude. Therefore, we must conclude that even the very simple PtQ inversions are affected by the presence of the magnetic field in the sunspot in a much larger region, and as a consequence the topology of the flows is different from the PtArc one-sided inversions. Fig. B.4. Radial profiles of azimuthally averaged speeds. The local minima of the speed are at distances 37.3, 34.6, 30.9, 28.1, and 28.5 Mm for PtA, PtQ-1, PtQ-2, PtArc, and CST respectively. The green vertical lines indicate the edges of umbra and penumbra.