Context. The spatial power spectrum of supergranulation does not fully characterize the underlying physics of turbulent convection. For example, it does not describe the non-Gaussianity in the horizontal flow divergence. Aims. Our aim is to statistically characterize the spatial pattern of solar supergranulation beyond the power spectrum. The next-order statistic is the bispectrum. It measures correlations of three Fourier components and is related to the nonlinearities in the underlying physics. Methods. We estimated the bispectrum of supergranular horizontal surface divergence maps that were obtained using local correlation tracking (LCT) and time-distance helioseismology (TD) from one year of data from the Helioseismic and Magnetic Imager on-board the Solar Dynamics Observatory starting in May 2010. Results. We find significantly nonzero and consistent estimates for the bispectrum. The strongest nonlinearity is present when the three coupling wave vectors are at the supergranular scale. These are the same wave vectors that are present in regular hexagons, which were used in analytical studies of solar convection. At these Fourier components, the bispectrum is positive, consistent with the positive skewness in the data and with supergranules preferentially consisting of outflows surrounded by a network of inflows. We use the bispectrum to generate synthetic divergence maps that are very similar to the data by a model that consists of a Gaussian term and a weaker quadratic nonlinear component. Thereby, we estimate the fraction of the variance in the divergence maps from the nonlinear component to be of the order of 4-6%. Conclusions. We propose that bispectral analysis is useful for understanding solar turbulent convection, for example for comparing observations and numerical models of supergranular flows. This analysis may also be useful to generate synthetic flow fields.

Context. Solar gravity modes (g modes) are buoyancy waves trapped in the solar radiative zone that have been very difficult to detect at the surface. Solar g modes would complement solar pressure modes (p modes) in probing the central regions of the Sun, for example the core rotation rate. Aims. A detection of g modes using changes in the large frequency separation of p modes has recently been reported. However, it is unclear how p and g modes interact. The aim of this study is to evaluate to what extent g modes can perturb the frequencies of p modes. Methods. We computed the first-order perturbation to global p-mode frequencies due to a flow field and perturbations to solar structure caused by a g mode. We focused on long-period g modes and assumed that the g-mode perturbations are constant in time. The surface amplitude of g modes is assumed to be $1$ mm s$^{-1}$, which is close to the observational limit. Results. Gravity modes do perturb p-mode frequencies to first order if the harmonic degree of the g mode is even and if its azimuthal order is zero. The effect is extremely small. For dipole and quadrupole p modes, all frequency shifts are smaller than $0.1$ nHz, or $2\times10^{-8}$ in relative numbers. This is because the relative perturbation to solar structure quantities caused by a g mode of realistic amplitude is of the order of $10^{-6}$ to $10^{-5}$. We find that structural changes dominate over advection. Surprisingly, the interaction of g and p modes takes place to a large part near the surface, where p modes spend most of their propagation times and g modes generate the largest relative changes to solar structure. Conclusions. It appears to be impossible to detect g modes solely through their signature in p-mode frequency shifts. Whether g modes leave a detectable signature in p-mode travel times under a given observational setup remains an open question.

The solar meridional flow is a crucial ingredient in modern dynamo theory. Seismic estimates of this flow have, however, been contradictory in deeper layers below about $0.9\,R_\odot$. Results from time-distance helioseismology have so far been obtained using the ray approximation. Here, we perform inversions using the Born approximation. The initial result is similar to the result previously obtained by Jackiewicz et al. (2015) using ray kernels while using the same set of GONG data and the SOLA inversion technique. However, we show that the assumption of uncorrelated measurements used in earlier studies may lead to inversion errors being underestimated by a factor of about two to four. In a second step, refined inversions are performed using the full covariance matrix and a regularization for cross-talk. As the results are found to depend on the threshold used in the singular value decomposition, they were obtained for a medium threshold ($10^{-7}-10^{-5}$, about 50% of the values used) and a threshold lower by a factor of 10 (about 70% of the values used). The result obtained with the medium threshold is again similar to the original, with less latitudinal variation. However, using the lower threshold, the inverted flow in the southern hemisphere shows two or three cells stacked radially depending on the associated radial flows. Both the single-cell and the multi-cell profiles are consistent with the measured travel times. All our results confirm a shallow return flow at about $0.9\,R_\odot$.

Accurate measurements of deep solar meridional flow are of vital interest for understanding the solar dynamo. In this paper, we validate a recently developed method for obtaining sensitivity functions (kernels) for travel-time measurements to solar interior flows using the Born approximation in spherical geometry, which is expected to be more accurate than the classical ray approximation. Furthermore, we develop a numerical approach to efficiently compute a large number of kernels based on the separability of the eigenfunctions into their horizontal and radial dependence. The validation is performed using a hydrodynamic simulation of linear wave propagation in the Sun, which includes a standard single-cell meridional flow profile. We show that, using the Born approximation, it is possible to accurately model observational quantities relevant for time-distance helioseismology such as the mean power spectrum, disc-averaged cross-covariance functions, and travel times in the presence of a flow field. In order to closely match the model to observations, we show that it is beneficial to use mode frequencies and damping rates which were extracted from the measured power spectrum. Furthermore, the contribution of the radial flow to the total travel time is found to reach 20% of the contribution of the horizontal flow at travel distances over $40^\circ$. Using the Born kernels and a 2D SOLA inversion of travel times, we can recover most features of the input meridional flow profile. The Born approximation is thus a promising method for inferring large-scale solar interior flows.

We extend an existing Born approximation method for calculating the linear sensitivity of helioseismic travel times to flows from Cartesian to spherical geometry. This development is necessary for using the Born approximation for inferring large-scale flows in the deep solar interior. In a first sanity check, we compare two $f-$mode kernels from our spherical method and from an existing Cartesian method. The horizontal and total integrals agree to within 0.3 %. As a second consistency test, we consider a uniformly rotating Sun and a travel distance of 42 degrees. The analytical travel-time difference agrees with the forward-modelled travel-time difference to within 2 %. In addition, we evaluate the impact of different choices of filter functions on the kernels for a meridional travel distance of 42 degrees. For all filters, the sensitivity is found to be distributed over a large fraction of the convection zone. We show that the kernels depend on the filter function employed in the data analysis process. If modes of higher harmonic degree ($90\lesssim l \lesssim 170$) are permitted, a noisy pattern of a spatial scale corresponding to $l\approx 260$ appears near the surface. When mainly low-degree modes are used ($l\lesssim70$), the sensitivity is concentrated in the deepest regions and it visually resembles a ray-path-like structure. Among the different low-degree filters used, we find the kernel for phase-speed filtered measurements to be best localized in depth.