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 Rossby waves (r modes) have recently been discovered in the near-surface horizontal flow field using the techniques of granulation-tracking and ring-diagram analysis applied to six years of SDO/HMI data. Aims. Here we apply time-distance helioseismology to the combined SOHO/MDI and SDO/HMI data sets, which cover 21 years of observations from May 1996 to April 2017. The goal of this study is to provide an independent confirmation over two solar cycles and in deeper layers of the Sun. Methods. We have measured south-north helioseismic travel times along the equator, which are sensitive to subsurface north-south flows. To reduce noise, the travel times were averaged over travel distances from 6$^\circ$ to 30$^\circ$; the mean distance corresponds to a p-mode lower turning point of 0.91 $R_\odot$. The 21-year time series of travel-time measurements was split into three seven-year subsets and transformed to obtain power spectra in a corotating frame. Results. The power spectra all show peaks near the frequencies of the classical sectoral Rossby waves for azimuthal wavenumbers in the range $3 \leq m \leq 15$. The mode frequencies and linewidths of the modes with $m \leq 9$ are consistent with a previous study whereas modes with $m \geq 10$ are shifted toward less negative frequencies by 10--20 nHz. While most of these modes have e-folding lifetimes on the order of a few months, the longest lived mode, $m=3$, has an e-folding lifetime of more than one year. For each mode, the rms velocity at the equator is in the range of 1--3 m s$^{-1}$ , with the largest values for $m\sim10$. No evidence for the $m=2$ sectoral mode is found in the power spectrum, implying that the rms velocity of this mode is below $\sim$0.5 m s$^{-1}$.

The south-north travel-time differences are measured by applying time-distance helioseismology to the MDI and HMI medium-degree Dopplergrams covering May 1996-April 2017. Our data analysis corrects for several sources of systematic effects: P-angle error, surface magnetic field effects, and center-to-limb variations. An interpretation of the travel-time measurements is obtained using a forward-modeling approach in the ray approximation. The travel-time differences are similar in the southern hemisphere for cycles 23 and 24. However, they differ in the northern hemisphere between cycles 23 and 24. Except for cycle 24's northern hemisphere, the measurements favor a single-cell meridional circulation model where the poleward flows persist down to $\sim$0.8 $R_\odot$, accompanied by local inflows toward the activity belts in the near-surface layers. Cycle 24's northern hemisphere is anomalous: travel-time differences are significantly smaller when travel distances are greater than 20$^\circ$. This asymmetry between northern and southern hemispheres during cycle 24 was not present in previous measurements (e.g., Rajaguru & Antia 2015), which assumed a different P-angle error correction where south-north travel-time differences are shifted to zero at the equator for all travel distances. In our measurements, the travel-time differences at the equator are zero for travel distances less than $\sim$30$^\circ$, but they do not vanish for larger travel distances. This equatorial offset for large travel distances need not be interpreted as a deep cross-equator flow; it could be due to the presence of asymmetrical local flows at the surface near the end points of the acoustic ray paths.

We present a new method for helioseismic diagnostics of the three-dimensional structure of sound speed, magnetic fields and flow velocities in the convection zone by inversion of acoustic travel-time data. The data are measurements of the time for acoustic waves to travel between points on the solar surface and surrounding annuli obtained from continuous observations at the South Pole in 1991 and from high-resolution observations from the Solar and Heliospheric Observatory (SOHO) in 1996. The travel time of the waves depends primarily on the sound speed perturbations and the velocity of flow along the ray paths. The effects of the sound speed perturbations and flows can be separated by measuring the travel time of waves propagating in opposite directions along the same ray paths. Magnetic fields result in anisotropy of the wave speed. A 3D inversion method based on Fermat's Principle and a regularized least-squares technique have been applied to infer the properties of convection in the quiet Sun and in active regions.

Previous helioseismology of sunspots has been sensitive to both the structural and magnetic aspects of sunspot structure. We aim to develop a technique that is insensitive to the magnetic component so the two aspects can be more readily separated. We study waves reflected almost vertically from the underside of a sunspot. Time-distance helioseismology was used to measure travel times for the waves. Ray theory and a detailed sunspot model were used to calculate travel times for comparison. It is shown that these large distance waves are insensitive to the magnetic field in the sunspot. The largest travel time differences for any solar phenomena are observed. With sufficient modeling effort, these should lead to better understanding of sunspot structure.

The Sun's complex dynamics is controlled by buoyancy and rotation in the convection zone and by magnetic forces in the atmosphere and corona. While small-scale solar convection is well understood, the dynamics of large-scale flows in the solar convection zone is not explained by theory or simulations. Waves of vorticity due to the Coriolis force, known as Rossby waves, are expected to remove energy out of convection at the largest scales. Here we unambiguously detect and characterize retrograde-propagating vorticity waves in the shallow subsurface layers of the Sun at angular wavenumbers below fifteen, with the dispersion relation of textbook sectoral Rossby waves. The waves have lifetimes of several months, well-defined mode frequencies below 200 nHz in a co-rotating frame, and eigenfunctions of vorticity that peak at the equator. Rossby waves have nearly as much vorticity as the convection at the same scales, thus they are an essential component of solar dynamics. We find a transition from turbulence-like to wave-like dynamics around the Rhines scale of angular wavenumber of twenty; this might provide an explanation for the puzzling deficit of kinetic energy at the largest spatial scales.

Context. The solar dynamo consists of a process that converts poloidal field to toroidal field followed by a process which creates new poloidal field from the toroidal field. Aims. Our aim is to observe the poloidal and toroidal fields relevant to the global solar dynamo and see if their evolution is captured by a Babcock-Leighton dynamo. Methods. We use synoptic maps of the surface radial field from the KPNSO/VT and SOLIS observatories to construct the poloidal field as a function of time and latitude, and Wilcox Solar Observatory and SOHO/MDI full disk images to infer the longitudinally averaged surface azimuthal field. We show that the latter is consistent with an estimate of that due to flux emergence and therefore closely related to the subsurface toroidal field. Results. We present maps of the poloidal and toroidal magnetic field of the global solar dynamo. The longitude-averaged azimuthal field observed at the surface results from flux emergence. At high latitudes this component follows the radial component of the polar fields with a short time lag (1-3 years). The lag increases at lower latitudes. The observed evolution of the poloidal and toroidal magnetic fields is described by the (updated) Babcock-Leighton dynamo model.

Context. Local helioseismology has detected spatially extended converging surface flows into solar active regions. These play an important role in flux-transport models of the solar dynamo. Aims. We aim to validate the existence of the inflows by deriving horizontal flow velocities around active regions with local correlation tracking of granulation. Methods. We generate a six-year long-time series of full-disk maps of the horizontal velocity at the solar surface by tracking granules in continuum intensity images provided by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). Results. On average, active regions are surrounded by inflows extending up to 10 deg from the center of the active region of magnitudes of 20-30 m/s, reaching locally up to 40 m/s, which is in agreement with results from local helioseismology. By computing an ensemble average consisting of 243 individual active regions, we show that the inflows are not azimuthally symmetric but converge predominantly towards the trailing polarity of the active region with respect to the longitudinally and temporally averaged flow field.

Time-distance helioseismology is one of the primary tools for studying the solar meridional circulation. However, travel-time measurements of the subsurface meridional flow suffer from a variety of systematic errors, such as a center-to-limb variation and an offset due to the P-angle uncertainty of solar images. Here we apply the time-distance technique to contemporaneous medium-degree Dopplergrams produced by SOHO/MDI and SDO/HMI to obtain the travel-time difference caused by meridional circulation throughout the solar convection zone. The P-angle offset in MDI images is measured by cross-correlating MDI and HMI images. The travel-time measurements in the south-north and east-west directions are averaged over the same observation period for the two data sets and then compared to examine the consistency of MDI and HMI travel times after correcting the systematic errors. The offsets in the south-north travel-time difference from MDI data induced by the P-angle error gradually diminish with increasing travel distance. However, these offsets become noisy for travel distances corresponding to waves that reach the base of the convection zone. This suggests that a careful treatment of the P-angle problem is required when studying a deep meridional flow. After correcting the P-angle and the removal of the center-to-limb effect, the travel-time measurements from MDI and HMI are consistent within the error bars for meridional circulation covering the entire convection zone. The fluctuations observed in both data sets are highly correlated and thus indicate their solar origin rather than an instrumental origin. Although our results demonstrate that the ad hoc correction is capable of reducing the wide discrepancy in the travel-time measurements from MDI and HMI, we cannot exclude the possibility that there exist other systematic effects acting on the two data sets in the same way.

NASA's Solar Dynamics Observatory (SDO) was launched 11 February 2010 with 3 instruments on board, including the Helioseismic and Magnetic Imager (HMI). Since beginning normal operations on 1 May 2010, HMI has observed the Sun's entire visible disk almost continuously. HMI collects sequences of polarized filtergrams taken at a fixed cadence with two 4096 x 4096 cameras from which are computed arcsecond-resolution maps of photospheric observables: the line-of-sight (LoS) velocity and magnetic field, continuum intensity, line width, line depth, and the Stokes polarization parameters, I Q U V, at 6 wavelengths. Two processing pipelines implemented at the SDO Joint Science Operations Center (JSOC) at Stanford University compute observables from calibrated Level-1 filtergrams. One generates LoS quantities every 45s, and the other, primarily for the vector magnetic field, computes averages on a 720s cadence. Corrections are made for static and temporally changing CCD characteristics, bad pixels, image alignment and distortion, polarization irregularities, filter-element uncertainty and non-uniformity, as well as Sun-spacecraft velocity. This report explains issues affecting the resulting physical quantities, describes the impact of regular updates to the instrument calibration, and shows how the computations are optimized for actual HMI observations. During the 5 years of the SDO prime mission, regular calibration sequences have been used to regularly improve and update the instrument calibration and to monitor instrument changes. The instrument more than satisfies the original specifications for data quality and continuity. The procedures described here still have significant room for improvement. The most significant remaining systematic errors are associated with the spacecraft orbital velocity.

Context. Local correlation tracking of granulation (LCT) is an important method for measuring horizontal flows in the photosphere. This method exhibits a systematic error that looks like a flow converging towards disk center, also known as the shrinking-Sun effect. Aims. We aim at studying the nature of the shrinking-Sun effect for continuum intensity data and at deriving a simple model that can explain its origin. Methods. We derived LCT flow maps by running the local correlation tracking code FLCT on tracked and remapped continuum intensity maps provided by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory. We also computed flow maps from synthetic continuum images generated from STAGGER code simulations of solar surface convection. We investigated the origin of the shrinking-Sun effect by generating an average granule from synthetic data from the simulations. Results. The LCT flow maps derived from HMI and from the simulations exhibit a shrinking-Sun effect of comparable magnitude. The origin of this effect is related to the apparent asymmetry of granulation originating from radiative transfer effects when observing with a viewing angle inclined from vertical. This causes, in combination with the expansion of the granules, an apparent motion towards disk center.

We review recent advances and results in enhancing and developing helioseismic analysis methods and in solar data assimilation. In the first part of this paper we will focus on selected developments in time-distance and global helioseismology. In the second part, we review the application of data assimilation methods on solar data. Relating solar surface observations as well as helioseismic proxies with solar dynamo models by means of the techniques from data assimilation is a promising new approach to explore and to predict the magnetic activity cycle of the Sun.

Context. Several upcoming and proposed space missions, such as Solar Orbiter, will be limited in telemetry and thus require data compression. Aims. We test the impact of data compression on local correlation tracking (LCT) of time-series of continuum intensity images. We evaluate the effect of several lossy compression methods (quantization, JPEG compression, and a reduced number of continuum images) on measurements of solar differential rotation with LCT. Methods. We apply the different compression methods to tracked and remapped continuum intensity maps obtained by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory. We derive 2D vector velocities using the local correlation tracking code FLCT and determine the additional bias and noise introduced by compression to differential rotation. Results. We find that probing differential rotation with LCT is very robust to lossy data compression when using quantization. Our results are severely affected by systematic errors of the LCT method and the HMI instrument. The sensitivity of LCT to systematic errors is a concern for Solar Orbiter.

Recently, Duvall and Hanasoge ({\it Solar Phys.} {\bf 287}, 71-83, 2013) found that large distance $[\Delta]$ separation travel-time differences from a center to an annulus $[\delta t_{\rm{oi}}]$ implied a model of the average supergranular cell that has a peak upflow of $240\rm{ms^{-1}}$ at a depth of $2.3\rm{Mm}$ and a corresponding peak outward horizontal flow of $700\rm{ms^{-1}}$ at a depth of $1.6\rm{Mm}$. In the present work, this effect is further studied by measuring and modeling center-to-quadrant travel-time differences $[\delta t_{\rm{qu}}]$, which roughly agree with this model. Simulations are analyzed that show that such a model flow would lead to the expected travel-time differences. As a check for possible systematic errors, the center-to-annulus travel-time differences $[\delta t_{\rm{oi}}]$ are found not to vary with heliocentric angle. A consistency check finds an increase of $\delta t_{\rm{oi}}$ with the temporal frequency $[\nu]$ by a factor of two, which is not predicted by the ray theory.

Convection in the solar interior is thought to comprise structures on a spectrum of scales. This conclusion emerges from phenomenological studies and numerical simulations, though neither covers the proper range of dynamical parameters of solar convection. Here, we analyze observations of the wavefield in the solar photosphere using techniques of time-distance helioseismology to image flows in the solar interior. We downsample and synthesize 900 billion wavefield observations to produce 3 billion cross-correlations, which we average and fit, measuring 5 million wave travel times. Using these travel times, we deduce the underlying flow systems and study their statistics to bound convective velocity magnitudes in the solar interior, as a function of depth and spherical-harmonic degree $\ell$. Within the wavenumber band $\ell<60$, Convective velocities are 20-100 times weaker than current theoretical estimates. This suggests the prevalence of a different paradigm of turbulence from that predicted by existing models, prompting the question: what mechanism transports the heat flux of a solar luminosity outwards? Advection is dominated by Coriolis forces for wavenumbers $\ell<60$, with Rossby numbers smaller than $\sim10^{-2}$ at $r/R_\odot=0.96$, suggesting that the Sun may be a much faster rotator than previously thought, and that large-scale convection may be quasi-geostrophic. The fact that iso-rotation contours in the Sun are not co-aligned with the axis of rotation suggests the presence of a latitudinal entropy gradient.

We report on a systematic center-to-limb variation in measured helioseismic travel times, which must be taken into account for an accurate determination of solar interior meridional flows. The systematic variation, found in time-distance helioseismology analysis using SDO/HMI and SDO/AIA observations, is different in both travel-time magnitude and variation trend for different observables. It is not clear what causes this systematic effect. Subtracting the longitude-dependent east-west travel times, obtained along the equatorial area, from the latitude-dependent north-south travel times, obtained along the central meridian area, gives remarkably similar results for different observables. We suggest this as an effective procedure for removing the systematic center-to-limb variation. The subsurface meridional flows obtained from inversion of the corrected travel times are approximately 10 m/s slower than those obtained without removing the systematic effect. The detected center-to-limb variation may have important implications in the derivation of meridional flows in the deep interior, and needs a better understanding.

The Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory (SDO/HMI) provides continuous full-disk observations of solar oscillations. We develop a data-analysis pipeline based on the time-distance helioseismology method to measure acoustic travel times using HMI Doppler-shift observations, and infer solar interior properties by inverting these measurements. The pipeline is used for routine production of near-real-time full-disk maps of subsurface wave-speed perturbations and horizontal flow velocities for depths ranging from 0 to 20 Mm, every eight hours. In addition, Carrington synoptic maps for the subsurface properties are made from these full-disk maps. The pipeline can also be used for selected target areas and time periods. We explain details of the pipeline organization and procedures, including processing of the HMI Doppler observations, measurements of the travel times, inversions, and constructions of the full-disk and synoptic maps. Some initial results from the pipeline, including full-disk flow maps, sunspot subsurface flow fields, and the interior rotation and meridional flow speeds, are presented.

Sunspot regions often form complexes of activity that may live for several solar rotations, and represent a major component of the Sun's magnetic activity. It had been suggested that the close appearance of active regions in space and time might be related to common subsurface roots, or "nests" of activity. EUV images show that the active regions are magnetically connected in the corona, but subsurface connections have not been established. We investigate the subsurface structure and dynamics of a large complex of activity, NOAA 10987-10989, observed during the SOHO/MDI Dynamics run in March-April 2008, which was a part of the Whole Heliospheric Interval (WHI) campaign. The active regions in this complex appeared in a narrow latitudinal range, probably representing a subsurface toroidal flux tube. We use the MDI full-disk Dopplergrams to measure perturbations of travel times of acoustic waves traveling to various depths by using time-distance helioseismology, and obtain sound-speed and flow maps by inversion of the travel times. The subsurface flow maps show an interesting dynamics of decaying active regions with persistent shearing flows, which may be important for driving the flaring and CME activity, observed during the WHI campaign. Our analysis, including the seismic sound-speed inversion results and the distribution of deep-focus travel-time anomalies, gave indications of diverging roots of the magnetic structures, as could be expected from $\Omega$-loop structures. However, no clear connection in the depth range of 0-48 Mm among the three active regions in this complex of activity was detected.

Local helioseismology provides unique information about the subsurface structure and dynamics of sunspots and active regions. However, because of complexity of sunspot regions local helioseismology diagnostics require careful analysis of systematic uncertainties and physical interpretation of the inversion results. We present new results of comparison of the ring-diagram analysis and time-distance helioseismology for active region NOAA 9787, for which a previous comparison showed significant differences in the subsurface sound-speed structure, and discuss systematic uncertainties of the measurements and inversions. Our results show that both the ring-diagram and time-distance techniques give qualitatively similar results, revealing a characteristic two-layer seismic sound-speed structure consistent with the results for other active regions. However, a quantitative comparison of the inversion results is not straightforward. It must take into account differences in the sensitivity, spatial resolution and the averaging kernels. In particular, because of the acoustic power suppression, the contribution of the sunspot seismic structure to the ring-diagram signal can be substantially reduced. We show that taking into account this effect reduces the difference in the depth of transition between the negative and positive sound-speed variations inferred by these methods. Further detailed analysis of the sensitivity, resolution and averaging properties of the local helioseismology methods is necessary for consolidation of the inversion results. It seems to be important that both methods indicate that the seismic structure of sunspots is rather deep and extends to at least 20 Mm below the surface, putting constraints on theoretical models of sunspots.

While sunspots are easily observed at the solar surface, determining their subsurface structure is not trivial. There are two main hypotheses for the subsurface structure of sunspots: the monolithic model and the cluster model. Local helioseismology is the only means by which we can investigate subphotospheric structure. However, as current linear inversion techniques do not yet allow helioseismology to probe the internal structure with sufficient confidence to distinguish between the monolith and cluster models, the development of physically realistic sunspot models are a priority for helioseismologists. This is because they are not only important indicators of the variety of physical effects that may influence helioseismic inferences in active regions, but they also enable detailed assessments of the validity of helioseismic interpretations through numerical forward modeling. In this paper, we provide a critical review of the existing sunspot models and an overview of numerical methods employed to model wave propagation through model sunspots. We then carry out an helioseismic analysis of the sunspot in Active Region 9787 and address the serious inconsistencies uncovered by \citeauthor{gizonetal2009}~(\citeyear{gizonetal2009,gizonetal2009a}). We find that this sunspot is most probably associated with a shallow, positive wave-speed perturbation (unlike the traditional two-layer model) and that travel-time measurements are consistent with a horizontal outflow in the surrounding moat.

Various methods of helioseismology are used to study the subsurface properties of the sunspot in NOAA Active Region 9787. This sunspot was chosen because it is axisymmetric, shows little evolution during 20-28 January 2002, and was observed continuously by the MDI/SOHO instrument. (...) Wave travel times and mode frequencies are affected by the sunspot. In most cases, wave packets that propagate through the sunspot have reduced travel times. At short travel distances, however, the sign of the travel-time shifts appears to depend sensitively on how the data are processed and, in particular, on filtering in frequency-wavenumber space. We carry out two linear inversions for wave speed: one using travel-times and phase-speed filters and the other one using mode frequencies from ring analysis. These two inversions give subsurface wave-speed profiles with opposite signs and different amplitudes. (...) From this study of AR9787, we conclude that we are currently unable to provide a unified description of the subsurface structure and dynamics of the sunspot.

We constrain the velocity spectral distribution of global-scale solar convective cells at depth using techniques of local helioseismology. We calibrate the sensitivity of helioseismic waves to large-scale convective cells in the interior by analyzing simulations of waves propagating through a velocity snapshot of global solar convection via methods of time-distance helioseismology. Applying identical analysis techniques to observations of the Sun, we are able to bound from above the magnitudes of solar convective cells as a function of spatial convective scale. We find that convection at a depth of $r/R_\odot = 0.95$ with spatial extent $\ell <20$, where $\ell$ is the spherical harmonic degree, comprise weak flow systems, on the order of 15 m/s or less. Convective features deeper than $r/R_\odot = 0.95$ are more difficult to image due to the rapidly decreasing sensitivity of helioseismic waves.

We derive expectations for signatures in the measured travel times of waves that interact with thermal anomalies and jets. A series of numerical experiments that involve the dynamic linear evolution of an acoustic wave field in a solar-like stratified spherical shell in the presence of fully 3D time-stationary perturbations are performed. The imprints of these interactions are observed as shifts in wave travel times, which are extracted from these data through methods of time-distance helioseismology \citep{duvall}. In situations where at least one of the spatial dimensions of the scatterer was smaller than a wavelength, oscillatory time shift signals were recovered from the analyses, pointing directly to a means of resolving sub-wavelength features. As evidence for this claim, we present analyses of simulations with spatially localized jets and sound-speed perturbations. We analyze 1 years' worth solar observations to estimate the noise level associated with the time differences. Based on theoretical estimates, Fresnel zone time shifts associated with the (possible) sharp rotation gradient at the base of the convection zone are of the order 0.01 - 0.1 s, well below the noise level that could be reached with the currently available amount of data ($\sim 0.15-0.2$ s with 10 yrs of data).

The propagation of solar waves through the sunspot of AR 9787 is observed using temporal cross-correlations of SOHO/MDI Dopplergrams. We then use three-dimensional MHD numerical simulations to compute the propagation of wave packets through self-similar magneto-hydrostatic sunspot models. The simulations are set up in such a way as to allow a comparison with observed cross-covariances (except in the immediate vicinity of the sunspot). We find that the simulation and the f-mode observations are in good agreement when the model sunspot has a peak field strength of 3 kG at the photosphere, less so for lower field strengths. Constraining the sunspot model with helioseismology is only possible because the direct effect of the magnetic field on the waves has been fully taken into account. Our work shows that the full-waveform modeling of sunspots is feasible.

Time-distance helioseismology has shown that f-mode travel times contain information about horizontal flows in the Sun. The purpose of this study is to provide a simple interpretation of these travel times. We study the interaction of surface-gravity waves with horizontal flows in an incompressible, plane-parallel solar atmosphere. We show that for uniform flows less than roughly 250 m s$^{-1}$, the travel-time shifts are linear in the flow amplitude. For stronger flows, perturbation theory up to third order is needed to model waveforms. The case of small-amplitude spatially-varying flows is treated using the first-order Born approximation. We derive two-dimensional Fr\'{e}chet kernels that give the sensitivity of travel-time shifts to local flows. We show that the effect of flows on travel times depends on wave damping and on the direction from which the observations are made. The main physical effect is the advection of the waves by the flow rather than the advection of wave sources or the effect of flows on wave damping. We compare the two-dimensional sensitivity kernels with simplified three-dimensional kernels that only account for wave advection and assume a vertical line of sight. We find that the three-dimensional f-mode kernels approximately separate in the horizontal and vertical coordinates, with the horizontal variations given by the simplified two-dimensional kernels. This consistency between quite different models gives us confidence in the usefulness of these kernels for interpreting quiet-Sun observations.

It is well known that the observed amplitude of solar oscillations is lower in sunspots than in quiet regions of the Sun. We show that this local reduction in oscillation amplitudes combined with the phase-speed filtering procedure in time-distance helioseismic analyses could be a source of systematic errors in the range of 5 - 40% in the measured travel-time anomalies of acoustic waves around sunspots. Removing these travel time artifacts is important for correctly inferring the subsurface structure of sunspots. We suggest an empirical correction procedure and illustrate its usage for a small sunspot. This work utilizes data from MDI/SOHO.

Supergranulation on the surface of the Sun is an organized cellular flow pattern with a characteristic scale of 30 Mm. It is superficially similar to the well understood granulation that operates at the 1.5 Mm natural scale of convection, which has led to the conventional view that supergranulation has its origin in the convective motion of cells of gas, though this does not explain the observation that the supergranulation pattern appears to move faster around the Sun than the bulk of its surface. A wave origin has been proposed for supergranulation that may explain the superrotation, but it has never had much support. Here we report that the supergranulation pattern has oscillatory components with periods of 5-10 days, for which the best explanation is a spectrum of traveling waves. We show that there is excess power in the prograde and equatorward directions, which explains the observation of superrotation.