the data is largest.

Figure \[figbicoh\] displays results as a function of $\bk = (k_x,k_y)$ and $|\bmm|$ for selected values of $|\bmm|$. The results are averaged in such a way that $\bmm$ (orange arrow) is rotated to be in the $x$ direction and thereby is fixed; see Appendix \[appBispecEst\]. Each data point is then plotted at the location corresponding to $\bk$ relative to the direction of $\bmm$. For each $\bk$, the third coupling Fourier component is then given by $\bl=\bmm-\bk$. The strongest coupling is present at a location where the triad $\bk+\bl = \bmm$ corresponds to a hexagonal pattern at the supergranular scale; see middle panel in Figure \[figbicoh\]. Figures \[figbicohall\] and \[figbicohKLall\] illustrate these results further.

The results for TD are very similar for $|\bmm|R_\sun \lesssim 300$ (Figure \[figbicohTD\] and second column in Table \[tabResults\]). For $|\bmm|R_\sun > 300$, the signal considerably weakens compared to LCT and is dominated by noise and leakage due to the spatial smoothing. This is the case for both TD datasets.

Estimates for the bispectrum are qualitatively very similar to the bicoherence; see Figure \[figbispec\]. As the sign of the bispectrum is predominantly positive, the dominant part of the quadratic interaction results in a 