Geological evidence suggests liquid water near the Earth's surface as early as 4.4 gigayears ago when the faint young Sun only radiated about 70 % of its modern power output. At this point, the Earth should have been a global snowball. An extreme atmospheric greenhouse effect, an initially more massive Sun, release of heat acquired during the accretion process of protoplanetary material, and radioactivity of the early Earth material have been proposed as alternative reservoirs or traps for heat. For now, the faint-young-sun paradox persists as one of the most important unsolved problems in our understanding of the origin of life on Earth. Here we use astrophysical models to explore the possibility that the new-born Moon, which formed about 69 million years (Myr) after the ignition of the Sun, generated extreme tidal friction - and therefore heat - in the Hadean and possibly the Archean Earth. We show that the Earth-Moon system has lost about 3e31 J, (99 % of its initial mechanical energy budget) as tidal heat. Tidal heating of roughly 10 W/m^2 through the surface on a time scale of 100 Myr could have accounted for a temperature increase of up to 5 degrees Celsius on the early Earth. This heating effect alone does not solve the faint-young-sun paradox but it could have played a key role in combination with other effects. Future studies of the interplay of tidal heating, the evolution of the solar power output, and the atmospheric (greenhouse) effects on the early Earth could help in solving the faint-young-sun paradox.

The tilt of solar active regions described by Joy's law is essential for converting a toroidal field to a poloidal field in Babcock-Leighton dynamo models. In thin flux tube models the Coriolis force causes Joy's law, acting on east-west flows as they rise towards the surface. Our goal is to measure the evolution of the average tilt angle of hundreds of active regions as they emerge, so that we can constrain the origins of Joy's law. We measured the tilt angle of the primary bipoles in 153 emerging active regions in the Solar Dynamics Observatory Helioseismic Emerging Active Region survey. We used line-of-sight magnetic field measurements averaged over 6 hours to define the polarities and measure the tilt angle up to four days after emergence. We find that at the time of emergence the polarities are on average aligned east-west, and that neither the separation nor the tilt depends on latitude. We do find, however, that ARs at higher latitudes have a faster north-south separation speed than those closer to the equator at the emergence time. After emergence, the tilt angle increases and Joy's law is evident about two days later. The scatter in the tilt angle is independent of flux until about one day after emergence, when higher-flux regions have a smaller scatter in tilt angle than lower-flux regions. Our finding that active regions emerge with an east-west alignment is surprising since thin flux tube models predict that tilt angles of rising flux tubes are generated below the surface. Previously reported tilt angle relaxation of deeply anchored flux tubes can be largely explained by the change in east-west separation. We conclude that Joy's law is caused by an inherent north-south separation speed present when the flux first reaches the surface, and that the scatter in the tilt angle is consistent with buffeting of the polarities by supergranulation.

Many moons have been detected around planets in our Solar System, but none has been detected unambiguously around any of the confirmed extrasolar planets. We test the feasibility of a supervised convolutional neural network to classify photometric transit light curves of planet-host stars and identify exomoon transits, while avoiding false positives caused by stellar variability or instrumental noise. Convolutional neural networks are known to have contributed to improving the accuracy of classification tasks. The network optimization is typically performed without studying the effect of noise on the training process. Here we design and optimize a 1D convolutional neural network to classify photometric transit light curves. We regularize the network by the total variation loss in order to remove unwanted variations in the data features. Using numerical experiments, we demonstrate the benefits of our network, which produces results comparable to or better than the standard network solutions. Most importantly, our network clearly outperforms a classical method used in exoplanet science to identify moon-like signals. Thus the proposed network is a promising approach for analyzing real transit light curves in the future.

While the solar system contains about 20 times more moons than planets, no moon has been confirmed around any of the thousands of extrasolar planets known so far. Tools for an uncomplicated identification of the most promising exomoon candidates could be beneficial to streamline follow-up studies.} Here we study three exomoon indicators that emerge if well-established planet-only models are fitted to a planet-moon transit light curve: transit timing variations (TTVs), transit duration variations (TDVs), and apparent planetary transit radius variations (TRVs). We re-evaluate under realistic conditions the previously proposed exomoon signatures in the TTV and TDV series. We simulate light curves of a transiting exoplanet with a single moon. These model light curves are then fitted with a planet-only transit model, pretending there were no moon, and we explore the resulting TTV, TDV, and TRV series for evidence of the moon. The previously described ellipse in the TTV-TDV diagram of an exoplanet with a moon emerges only for high-density moons. Low-density moons distort the sinusoidal shapes of the TTV and the TDV series due to their photometric contribution to the combined planet-moon transit. Sufficiently large moons can produce periodic apparent TRVs of their host planets that could be observable. We find that Kepler and PLATO have similar performances in detecting the exomoon-induced TRV effect around simulated bright ($m_V=8$) stars. These stars, however, are rare in the Kepler sample but will be abundant in the PLATO sample. Moreover, PLATO's higher cadence yields a stronger TTV signal. The periodogram of the sequence of transit radius measurements can indicate the presence of a moon. The TTV and TDV series of exoplanets with moons can be more complex than previously assumed. We propose that TRVs could be a more promising means to identify exomoons in large exoplanet surveys.

Large-scale Rossby waves have recently been discovered from measurements of horizontal surface and near-surface solar flows (L\"optien at al. 2018). We are interested in understanding why only the sectoral modes are seen in the observations and also in modelling the radial structure of the observed modes. To do so, we characterise here the radial eigenfunctions of r modes for slowly-rotating polytropes in uniform rotation. We find that for free-surface boundary conditions on a spheroid of non-vanishing surface density, r modes can only exist for $\ell=m$ spherical harmonics in the inviscid case, and we compute their depth dependence and frequencies to leading order. For quasi-adiabatic stratification the sectoral modes with no radial nodes are the only modes which are almost toroidal and the depth dependence of the corresponding horizontal motion scales as $r^m$. For all r modes except the zero radial order sectoral ones, non-adiabatic stratification plays a crucial role in the radial force balance. The lack of quasi-toroidal solutions when stratification is close to neutral, except for the sectoral modes without nodes in radius, follows from the statement that the system needs to be in both horizontal and radial force balance. In the absence of super- or subadiabatic stratification and viscosity, both the horizontal and radial force balances independently determine the pressure perturbation. The only quasi-toroidal cases in which the two determinations of the pressure perturbation are consistent are the special cases where $\ell=m$, and the horizontal displacement scales with $r^m$.

The variability of the spectral solar irradiance (SSI) over the course of the 11-year solar cycle is one of the manifestations of solar magnetic activity. There is a strong evidence that the SSI variability has an effect on the Earth's atmosphere. The faster rotation of the Sun in the past lead to a more vigorous action of solar dynamo and thus potentially to larger amplitude of the SSI variability on the timescale of the solar activity cycle. This could led to a stronger response of the Earth's atmosphere as well as other solar system planets' atmospheres to the solar activity cycle. We calculate the amplitude of the SSI and TSI variability over the course of the solar activity cycle as a function of solar age. We employ the relationship between the stellar magnetic activity and the age based on observations of solar twins. Using this relation we reconstruct solar magnetic activity and the corresponding solar disk area coverages by magnetic features (i.e. spots and faculae) over the last four billion years. These disk coverages are then used to calculate the amplitude of the solar-cycle SSI variability as a function of wavelength and solar age. Our calculations show that the young Sun was significantly more variable than the present Sun. The amplitude of the solar-cycle Total Solar Irradiance (TSI) variability of the 600 Myr old Sun was about 10 times larger than that of the present Sun. Furthermore, the variability of the young Sun was spot-dominated (the Sun being brighter at the activity minimum than in the maximum), i.e. the Sun was overall brighter at activity minima than at maxima. The amplitude of the TSI variability decreased with solar age until it reached a minimum value at 2.8 Gyr. After this point, the TSI variability is faculae-dominated (the Sun is brighter at the activity maximum) and its amplitude increases with age.

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 sectoral Rossby waves have only recently been unambiguously identified in Helioseimsic and Magnetic Imager (HMI) and Michelson Doppler Imager (MDI) maps of flows near the solar surface. So far this has not been done with the Global Oscillation Network Group (GONG) ground-based observations, which have different noise properties. Aims: We utilize 17 years of GONG++ data, to identify and characterize solar Rossby waves using ring-diagram helioseismology. We compare directly with HMI ring-diagram analysis. Methods: Maps of the radial vorticity are obtained for flows within the top 2 Mm of the surface for 17 years of GONG++. The data is corrected for systematic effects including the annual periodicity related to the $B_0$ angle. We then compute the Fourier components of the radial vorticity of the flows in the co-rotating frame. We perform the same analysis on the HMI data that overlap in time. Results: We find that the solar Rossby waves have measurable amplitudes in the GONG++ sectoral power spectra for azimuthal orders between $m=3$ and $m=15$. The measured mode characteristics (frequencies, lifetimes and amplitudes) from GONG++ are consistent with the HMI measurements in the overlap period from 2010 to 2018 for $m\le9$. For higher-$m$ modes the amplitudes and frequencies agree within two sigmas. The signal-to-noise ratio of modes in GONG++ power spectra is comparable to HMI for $8\le m\le11$, but is lower by a factor of two for other modes. Conclusions: The GONG++ data provide a long and uniform data set to study solar global-scale Rossby waves from 2001.

While long period variables (LPVs) have been extensively investigated, especially with MACHO and OGLE data for the Magellanic Clouds, there still exist open questions in their pulsations regarding the excitation mechanisms, radial order and angular degree assignment. Here, we perform asteroseismic analyses on LPVs observed by the 4-year Kepler mission. Using a cross-correlation method, we detect unambiguous pulsation ridges associated with radial fundamental modes ($n=1$) and overtones ($n\geqslant2$), where the radial order assignment is made by using theoretical frequencies and observed frequencies. Our results confirm that the amplitude variability seen in semiregulars is consistent with oscillations being solar-like. We identify that the dipole modes, $l=1$, are dominant in the radial orders of $3\leq n \leq6$, and that quadrupole modes, $l=2$, are dominant in the first overtone $n=2$. A test of seismic scaling relations using Gaia DR2 parallaxes reveals the possibility that the relations break down when $\nu_{\rm max}$ $\lesssim$ 3 $\mu$Hz (R $\gtrsim$ 40 R$_{\odot}$, or log $\rm L/L_{\odot}$ $\gtrsim$ 2.6). Our homogeneous measurements of pulsation amplitude and period for 3213 LPVs will be very valuable for probing effects of pulsation on mass loss, in particular in those stars with periods around 60 days, which has been argued as a threshold of substantial pulsation-triggered mass loss.

Global-scale Rossby waves have recently been unambiguously identified on the Sun. Here we study the latitude and depth dependence of the Rossby wave eigenfunctions. By applying helioseismic ring-diagram analysis and granulation tracking to SDO/HMI observations, we compute maps of the radial vorticity of flows in the upper solar convection zone (down to depths of more than $16$ Mm). We use a Fourier transform in longitude to separate the different azimuthal orders m in the range $3 \le m \le 15$. At each $m$ we obtain the phase and amplitude of the Rossby waves as a function of depth using the helioseismic data. At each $m$ we also measure the latitude dependence of the eigenfunctions by calculating the covariance between the equator and other latitudes. We then study the horizontal and radial dependences of the radial vorticity eigenfunctions. The horizontal eigenfunctions are complex. As observed previously, the real part peaks at the equator and switches sign near $\pm 30^\circ$, thus the eigenfunctions show significant non-sectoral contributions. The imaginary part is smaller than the real part. The phase of the radial eigenfunctions varies by only roughly $\pm 5^\circ$ over the top $15$ Mm. The amplitude of the radial eigenfunctions decreases by about $10\%$ from the surface down to $8$ Mm (the region where ring-diagram analysis is most reliable, as seen by comparing with the rotation rate measured by global-mode seismology). The radial dependence of the radial vorticity eigenfunctions deduced from ring-diagram analysis is consistent with a power-law down to $8$ Mm and is unreliable at larger depths. However, the observations provide only weak constraints on the power-law exponents. For the real part, the latitude dependence of the eigenfunctions is consistent with previous work (using granulation tracking). The imaginary part is smaller than the real part but significantly nonzero.

Axisymmetric magnetic activity on the Sun and Sun-like stars increases the frequencies of the modes of acoustic oscillation. However, it is unclear how a corotating patch of activity affects the oscillations, since such a perturbation is unsteady in the frame of the observer. In this paper we qualitatively describe the asteroseismic signature of a large active region in the power spectrum of the dipole and quadrupole p modes. In the corotating frame of the active region, the perturbations due to (differential) rotation and the active region completely lift the $(2\ell + 1)$-fold azimuthal degeneracy of the frequency spectrum of modes with harmonic degree $\ell$. In the frame of the observer, the unsteady nature of the perturbation leads to the appearance of $(2\ell+1)^2$ peaks in the power spectrum of a multiplet. These peaks blend into each other to form asymmetric line profiles. In the limit of a small active region, we approximate the power spectrum of a multiplet in terms of $2\times(2\ell+1)$ peaks, whose amplitudes and frequencies depend on the latitude of the active region and the inclination angle of the star's rotation axis. In order to check the results and to explore the nonlinear regime, we also perform numerical simulations using the 3D time-domain pseudo-spectral linear pulsation code GLASS.

Context: There is a wide discrepancy in current estimates of the strength of convection flows in the solar interior obtained using different helioseismic methods applied to observations from SDO/HMI. The cause for these disparities is not known. Aims: As one step in the effort to resolve this discrepancy, we aim to characterize the multi-ridge fitting code for ring-diagram helioseismic analysis that is used to obtain flow estimates from local power spectra of solar oscillations. Methods: We updated the multi-ridge fitting code developed by Greer et al.(2014) to solve several problems we identified through our inspection of the code. In particular, we changed the merit function to account for the smoothing of the power spectra, model for the power spectrum, and noise estimates. We used Monte Carlo simulations to generate synthetic data and to characterize the noise and bias of the updated code by fitting these synthetic data. Results: The bias in the output fit parameters, apart from the parameter describing the amplitude of the p-mode resonances in the power spectrum, is below what can be measured from the Monte-Carlo simulations. The amplitude parameters are underestimated; this is a consequence of choosing to fit the logarithm of the averaged power. We defer fixing this problem as it is well understood and not significant for measuring flows in the solar interior. The scatter in the fit parameters from the Monte-Carlo simulations is well-modeled by the formal error estimates from the code. Conclusions: We document and demonstrate a reliable multi-ridge fitting method for ring-diagram analysis. The differences between the updated fitting results and the original results are less than one order of magnitude and therefore we suspect that the changes will not eliminate the aforementioned orders-of-magnitude discrepancy in the amplitude of convective flows in the solar interior.

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.

Our goal is to constrain models of active region formation by tracking the average motion of active region polarity pairs as they emerge onto the surface. We measured the motion of the two main opposite polarities in 153 emerging active regions (EARs) using line-of-sight magnetic field observations from the Solar Dynamics Observatory Helioseismic Emerging Active Region (SDO/HEAR) survey (Schunker et al. 2016). We first measured the position of each of the polarities eight hours after emergence and tracked their location forwards and backwards in time. We find that, on average, the polarities emerge with an east-west orientation and the separation speed between the polarities increases. At about 0.1 days after emergence, the average separation speed reaches a peak value of 229 +/- 11 m/s, and then starts to decrease, and about 2.5 days after emergence the polarities stop separating. We also find that the separation and the separation speed in the east-west direction are systematically larger for active regions with higher flux. Our results reveal two phases of the emergence process defined by the rate of change of the separation speed as the polarities move apart. Phase 1 begins when the opposite polarity pairs first appear at the surface, with an east-west alignment and an increasing separation speed. We define Phase 2 to begin when the separation speed starts to decrease, and ends when the polarities have stopped separating. This is consistent with the picture of Chen, Rempel, & Fan (2017): the peak of a flux tube breaks through the surface during Phase 1. During Phase 2 the magnetic field lines are straightened by magnetic tension, so that the polarities continue to move apart, until they eventually lie directly above their anchored subsurface footpoints.

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}$.

This paper describes the Polarimetric and Helioseismic Imager on the Solar Orbiter mission (SO/PHI), the first magnetograph and helioseismology instrument to observe the Sun from outside the Sun-Earth line. It is the key instrument meant to address the top-level science question: How does the solar dynamo work and drive connections between the Sun and the heliosphere? SO/PHI will also play an important role in answering the other top-level science questions of Solar Orbiter, as well as hosting the potential of a rich return in further science. SO/PHI measures the Zeeman effect and the Doppler shift in the FeI 617.3nm spectral line. To this end, the instrument carries out narrow-band imaging spectro-polarimetry using a tunable LiNbO_3 Fabry-Perot etalon, while the polarisation modulation is done with liquid crystal variable retarders (LCVRs). The line and the nearby continuum are sampled at six wavelength points and the data are recorded by a 2kx2k CMOS detector. To save valuable telemetry, the raw data are reduced on board, including being inverted under the assumption of a Milne-Eddington atmosphere, although simpler reduction methods are also available on board. SO/PHI is composed of two telescopes; one, the Full Disc Telescope (FDT), covers the full solar disc at all phases of the orbit, while the other, the High Resolution Telescope (HRT), can resolve structures as small as 200km on the Sun at closest perihelion. The high heat load generated through proximity to the Sun is greatly reduced by the multilayer-coated entrance windows to the two telescopes that allow less than 4% of the total sunlight to enter the instrument, most of it in a narrow wavelength band around the chosen spectral line.

Asteroseismology has undergone a profound transformation as a scientific field following the CoRoT and Kepler space missions. The latter is now yielding the first measurements of latitudinal differential rotation obtained directly from oscillation frequencies. Differential rotation is a fundamental mechanism of the stellar dynamo effect. Our goal is to measure the amount of differential rotation in the solar analogues 16 Cyg A and B, which are the components of a binary system. These stars are the brightest observed by Kepler and have therefore been extensively observed, with exquisite precision on their oscillation frequencies. We modelled the acoustic power spectrum of 16 Cyg A and B using a model that takes into account the contribution of differential rotation to the rotational frequency splitting. The estimation was carried out in a Bayesian setting. We then inverted these results to obtain the rotation profile of both stars under the assumption of a solar-like functional form. We observe that the magnitude of latitudinal differential rotation has a strong chance of being solar-like for both stars, their rotation rates being higher at the equator than at the pole. The measured latitudinal differential rotation, defined as the difference of rotation rate between the equator and the pole, is $320\pm269$ nHz and $440^{+363}_{-383}$ nHz for 16 Cyg A and B, respectively, confirming that the rotation rates of these stars are almost solar-like. Their equatorial rotation rates are $535\pm75$ nHz and $565_{-129}^{+150}$ nHz. Our results are in good agreement with measurements obtained from spectropolarimetry, spectroscopy, and photometry. We present the first conclusive measurement of latitudinal differential rotation for solar analogues. Their rotational profiles are very close to those of the Sun. These results depend weakly on the uncertainties of the stellar parameters.

Radial velocity (RV) measurements are used to search for planets orbiting late-type main-sequence stars and confirm the transiting planets. The most advanced spectrometers are approaching a precision of $\sim 10$ cm/s that implies the need to identify and correct for all possible sources of RV oscillations intrinsic to the star down to this level and possibly beyond. The recent discovery of global-scale equatorial Rossby waves in the Sun, also called r modes, prompted us to investigate their possible signature in stellar RV measurements. R modes are toroidal modes of oscillation whose restoring force is the Coriolis force and propagate in the retrograde direction in a frame that corotates with the star. The solar r modes with azimuthal orders $3 \leq m \lesssim 15$ were identified unambiguously because of their dispersion relation and their long e-folding lifetimes of hundreds of days. Here we simulate the RV oscillations produced by sectoral r modes with $2 \leq m \leq 5$ assuming a stellar rotation period of 25.54 days and a maximum amplitude of the surface velocity of each mode of 2 m/s. This amplitude is representative of the solar measurements, except for the $m=2$ mode which has not yet been observed. Sectoral r modes with azimuthal orders $m=2$ and $3$ would produce RV oscillations with amplitudes of 76.4 and 19.6 cm/s and periods of 19.16 and 10.22 days, respectively, for a star with an inclination of the rotation axis $i=60^{\circ}$. Therefore, they may produce rather sharp peaks in the Fourier spectrum of the radial velocity time series that could lead to spurious planetary detections. Sectoral r~modes may represent a source of confusion in the case of slowly rotating inactive stars that are preferential targets for RV planet search. The main limitation of the present investigation is the lack of observational constraint on the amplitude of the $m=2$ mode on the Sun.

The inversion of ring fit parameters to obtain subsurface flow maps in ring-diagram analysis for 8 years of SDO observations is computationally expensive, requiring ~3200 CPU hours. In this paper we apply machine learning techniques to the inversion in order to speed up calculations. Specifically, we train a predictor for subsurface flows using the mode fit parameters and the previous inversion results, to replace future inversion requirements. We utilize Artificial Neural Networks as a supervised learning method for predicting the flows in 15 degree ring tiles. To demonstrate that the machine learning results still contain the subtle signatures key to local helioseismic studies, we use the machine learning results to study the recently discovered solar equatorial Rossby waves. The Artificial Neural Network is computationally efficient, able to make future flow predictions of an entire Carrington rotation in a matter of seconds, which is much faster than the current ~31 CPU hours. Initial training of the networks requires ~3 CPU hours. The trained Artificial Neural Network can achieve a root mean-square error equal to approximately half that reported for the velocity inversions, demonstrating the accuracy of the machine learning (and perhaps the overestimation of the original errors from the ring-diagram pipeline). We find the signature of equatorial Rossby waves in the machine learning flows covering six years of data, demonstrating that small-amplitude signals are maintained. The recovery of Rossby waves in the machine learning flow maps can be achieved with only one Carrington rotation (27.275 days) of training data. We have shown that machine learning can be applied to, and perform more efficiently than the current ring-diagram inversion. The computation burden of the machine learning includes 3 CPU hours for initial training, then around 0.0001 CPU hours for future predictions.

During the solar magnetic activity cycle the emergence latitudes of sunspots change, leading to the well-known butterfly diagram. This phenomenon is poorly understood for other stars as starspot latitudes are generally unknown. The related changes in starspot rotation rates caused by latitudinal differential rotation can however be measured. Using the set of 3093 Kepler stars with activity cycles identified by Reinhold et al. (2017), we aim to study the temporal change in starspot rotation rates over magnetic activity cycles, and how this relates to the activity level, mean rotation rate, and effective temperature of the star. We measure the photometric variability as a proxy for the magnetic activity and the spot rotation rate in each quarter over the duration of the Kepler mission. We phase-fold these measurements with the cycle period. We perform averages over stars with comparable mean rotation rates and effective temperature at fixed activity-cycle phases. We detect a clear correlation between the variation of activity level and the variation of the starspot rotation rate. The sign and amplitude of this correlation depends on the mean stellar rotation and, to a lesser extent, on the effective temperature. For slowly rotating stars (with periods between 15-28 days) the starspot rotation rates are clearly anti-correlated with the level of activity during the activity cycles. A transition is seen at periods of 10-15 days, where stars with effective temperature above 4200K instead show positive correlation. Our results can be interpreted in terms of a stellar "butterfly diagram", but these appear different from the Sun's as the starspot rotation rates are either in phase or anti-phase with the activity level. Alternatively, the activity cycles seen by Kepler are short (around 2.5 years) and may therefore be secondary cycles, perhaps analogous to the solar quasi-biennial oscillations.

Stellar magnetic fields are poorly understood but are known to be important for stellar evolution and exoplanet habitability. They drive stellar activity, which is the main observational constraint on theoretical models for magnetic field generation and evolution. Starspots are the main manifestation of the magnetic fields at the stellar surface. In this study we measure the variation of their latitude with time, called a butterfly diagram in the solar case, for the solar analogue HD 173701 (KIC 8006161). To that effect, we use Kepler data, to combine starspot rotation rates at different epochs and the asteroseismically determined latitudinal variation of the stellar rotation rates. We observe a clear variation of the latitude of the starspots. It is the first time such a diagram is constructed using asteroseismic data.

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.

Helioseismic holography is an imaging technique used to study heterogeneities and flows in the solar interior from observations of solar oscillations at the surface. Holograms contain noise due to the stochastic nature of solar oscillations. We provide a theoretical framework for modeling signal and noise in Porter-Bojarski helioseismic holography. The wave equation may be recast into a Helmholtz-like equation, so as to connect with the acoustics literature and define the holography Green's function in a meaningful way. Sources of wave excitation are assumed to be stationary, horizontally homogeneous, and spatially uncorrelated. Using the first Born approximation we calculate holograms in the presence of perturbations in sound-speed, density, flows, and source covariance, as well as the noise level as a function of position. This work is a direct extension of the methods used in time-distance helioseismology to model signal and noise. To illustrate the theory, we compute the hologram intensity numerically for a buried sound-speed perturbation at different depths in the solar interior. The reference Green's function is obtained for a spherically-symmetric solar model using a finite-element solver in the frequency domain. Below the pupil area on the surface, we find that the spatial resolution of the hologram intensity is very close to half the local wavelength. For a sound-speed perturbation of size comparable to the local spatial resolution, the signal-to-noise ratio is approximately constant with depth. Averaging the hologram intensity over a number $N$ of frequencies above 3 mHz increases the signal-to-noise ratio by a factor nearly equal to the square root of $N$. This may not be the case at lower frequencies, where large variations in the holographic signal are due to the individual contributions of the long-lived modes of oscillation.

The differentially rotating outer layers of stars are thought to play a role in driving their magnetic activity, but the underlying mechanisms that generate and sustain differential rotation are poorly understood. We report the measurement of latitudinal differential rotation in the convection zones of 40 Sun-like stars using asteroseismology. For the most significant detections, the stars' equators rotate approximately twice as fast as their mid-latitudes. The latitudinal shear inferred from asteroseismology is much larger than predictions from numerical simulations.

Transit photometry of the exoplanet candidate Kepler-1625b has recently been interpreted to show hints of a moon. We aim to clarify whether the exomoon-like signal is really caused by a large object in orbit around Kepler-1625b. We explore several detrending procedures, i.e. polynomials and the Cosine Filtering with Autocorrelation Minimization (CoFiAM). We then supply a light curve simulator with the co-planar orbital dynamics of the system and fit the resulting planet-moon transit light curves to the Kepler data. We employ the Bayesian Information Criterion (BIC) to assess whether a single planet or a planet-moon system is a more likely interpretation of the light curve variations. We carry out a blind hare-and-hounds exercise using many noise realizations by injecting simulated transits into different out-of-transit parts of the original Kepler-1625 data: 100 sequences with 3 synthetic transits of a Kepler-1625b-like planet and 100 sequences with 3 synthetic transits of this planet with a Neptune-sized moon. The statistical significance and characteristics of the exomoon-like signal strongly depend on the detrending method, and the data chosen for detrending, and on the treatment of gaps in the light curve. Our injection-retrieval experiment shows evidence for moons in about 10% of those light curves that do not contain an injected moon. Strikingly, many of these false-positive moons resemble the exomoon candidate. We recover up to about half of the injected moons, depending on the detrending method, with radii and orbital distances broadly corresponding to the injected values. A $\Delta$BIC of -4.9 for the CoFiAM-based detrending indicates an exomoon around Kepler-1625b. This solution, however, is only one out of many and we find very different solutions depending on the details of the detrending method. It is worrying that the detrending is key to the interpretation of the data.

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.

Solar supergranulation presents us with many mysteries. For example, previous studies in spectral space found that supergranulation has wave-like properties. Here we study, in real space, the wave-like evolution of the average supergranule over a range of spatial scales (from 10 to 80 Mm). We complement this by characterizing the evolution of the associated network magnetic field. We use one year of data from the Helioseismic and Magnetic Imager (HMI) to measure horizontal near-surface flows near the solar equator by applying time-distance helioseismology on Dopplergrams and granulation tracking on intensity images. The average supergranule outflow (or inflow) is constructed by averaging over 10000 individual outflows (or inflows). The contemporaneous evolution of the magnetic field is studied with HMI line-of-sight observations. We confirm and extend previous measurements of the supergranular wave dispersion relation to angular wavenumbers in the range 50<kR<270. We find a plateau for kR>120. In real space, larger supergranules undergo oscillations with longer periods and lifetimes than smaller cells. We find excellent agreement between TD and LCT and obtain wave properties that are independent of the tracking rate. The observed network magnetic field follows the oscillations of the supergranular flows with a six-hour time lag. This behavior can be explained by computing the motions of corks carried by the supergranular flows. Signatures of supergranular waves in surface horizontal flows near the solar equator can be observed in real space. These oscillatory flows control the evolution of the network magnetic field, in particular they explain the recently discovered east-west anisotropy of the magnetic field around the average supergranule. Background flow measurements that we obtain from Doppler frequency shifts do not favor shallow models of supergranulation.

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.

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 are computed using the ray or the Born approximation. The Born approximation is preferable as it takes finite-wavelength effects into account, but can be computationally expensive. We propose a fast algorithm to compute travel-time sensitivity kernels under the assumption that the background solar medium is spherically symmetric. 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 thanks to the algorithm of Dong and Lemus (2002). 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. This new method allows for 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.

Oscillation properties are usually measured by fitting symmetric Lorentzian profiles to the power spectra of Sun-like stars. However the line profiles of solar oscillations have been observed to be asymmetrical for the Sun. The physical origin of this line asymmetry is not fully understood, although it should depend on the depth dependence of the source of wave excitation (convective turbulence) and details of the observable (velocity or intensity). For oscillations of the Sun, it has been shown that neglecting the asymmetry leads to systematic errors in the frequency determination. This could subsequently affects the results of seismic inferences of the solar internal structure. Using light curves from the {\it Kepler} spacecraft we have measured mode asymmetries in 43 stars. We confirm that neglecting the asymmetry leads to systematic errors that can exceed the $1\sigma$ confidence intervals for seismic observations longer than one year. Therefore, the application of an asymmetric Lorentzian profile is to be favoured to improve the accuracy of the internal stellar structure and stellar fundamental parameters. We also show that the asymmetry changes sign between cool Sun-like stars and hotter stars. This provides the best constraints to date on the location of the excitation sources across the Hertzsprung-Russel diagram.

The internal gravity modes of the Sun are notoriously difficult to detect, and the claimed detection of gravity modes presented in Fossat et al. 2017 is thus very exciting. Given the importance of these modes for understanding solar structure and dynamics, the results must be robust. While Fossat et al. 2017 described their method and parameter choices in detail, the sensitivity of their results to several parameters were not presented. Therefore, we test the sensitivity to a selection of them. The most concerning result is that the detection vanishes when we adjust the start time of the 16.5 year velocity time series by a few hours. We conclude that this reported detection of gravity modes is extremely fragile and should be treated with utmost caution.

The purpose of deep-focusing time--distance helioseismology is to construct seismic measurements that have a high sensitivity to the physical conditions at a desired target point in the solar interior. With this technique, pairs of points on the solar surface are chosen such that acoustic ray paths intersect at this target (focus) point. Considering acoustic waves in a homogeneous medium, we compare travel-time and amplitude measurements extracted from the deep-focusing cross-covariance functions. Using a single-scattering approximation, we find that the spatial sensitivity of deep-focusing travel times to sound-speed perturbations is zero at the target location and maximum in a surrounding shell. This is unlike the deep-focusing amplitude measurements, which have maximum sensitivity at the target point. We compare the signal-to-noise ratio for travel-time and amplitude measurements for different types of sound-speed perturbations, under the assumption that noise is solely due to the random excitation of the waves. We find that, for highly localized perturbations in sound speed, the signal-to-noise ratio is higher for amplitude measurements than for travel-time measurements. We conclude that amplitude measurements are a useful complement to travel-time measurements in time--distance helioseismology.

The observed magnetic field on the solar surface is characterized by a very complex spatial and temporal behavior. Although feature-tracking algorithms have allowed us to deepen our understanding of this behavior, subjectivity plays an important role in the identification and tracking of such features. In this paper, we continue studies Gorobets, A. Y., Borrero, J. M., & Berdyugina, S. 2016, ApJL, 825, L18 of the temporal stochasticity of the magnetic field on the solar surface without relying either on the concept of magnetic features or on subjective assumptions about their identification and interaction. We propose a data analysis method to quantify fluctuations of the line-of-sight magnetic field by means of reducing the temporal field's evolution to the regular Markov process. We build a representative model of fluctuations converging to the unique stationary (equilibrium) distribution in the long time limit with maximum entropy. We obtained different rates of convergence to the equilibrium at fixed noise cutoff for two sets of data. This indicates a strong influence of the data spatial resolution and mixing-polarity fluctuations on the relaxation process. The analysis is applied to observations of magnetic fields of the relatively quiet areas around an active region carried out during the second flight of the Sunrise/IMaX and quiet Sun areas at the disk center from the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory satellite.

The temporal covariance between seismic waves measured at two locations on the solar surface is the fundamental observable in time-distance helioseismology. Above the acoustic cut-off frequency ($\sim$5.3~mHz), waves are not trapped in the solar interior and the covariance function can be used to probe the upper atmosphere. We wish to implement appropriate radiative boundary conditions for computing the propagation of high-frequency waves in the solar atmosphere. We consider the radiative boundary conditions recently developed by Barucq et al. (2017) for atmospheres in which sound-speed is constant and density decreases exponentially with radius. We compute the cross-covariance function using a finite element method in spherical geometry and in the frequency domain. The ratio between first- and second-skip amplitudes in the time-distance diagram is used as a diagnostic to compare boundary conditions and to compare with observations. We find that a boundary condition applied 500 km above the photosphere and derived under the approximation of small angles of incidence accurately reproduces the `infinite atmosphere' solution for high-frequency waves. When the radiative boundary condition is applied 2 Mm above the photosphere, we find that the choice of atmospheric model affects the time-distance diagram. In particular, the time-distance diagram exhibits double-ridge structure when using a VAL atmospheric model.

We develop a helioseismic inversion algorithm that can be used to recover sub-surface vertical profiles of 2-dimensional supergranular flows from surface measurements of synthetic wave travel times. We carry out seismic wave-propagation simulations through a 2-dimensional section of a flow profile that resembles an averaged supergranule, and a starting model that has flows only at the surface. We assume that the wave measurements are entirely without realization noise for the purpose of our test. We expand the vertical profile of the supergranule stream function on a basis of B-splines. We iteratively update the B-spline coefficients of the supergranule model to reduce the travel-times differences observed between the two simulations. We carry out the exercise for four different vertical profiles peaking at different depths below the solar surface. We are able to accurately recover depth profiles of four supergranule models at depths up to $8-10\,\text{Mm}$ below the solar surface using $f-p_4$ modes, under the assumption that there is no realization noise. We are able to obtain the peak depth and the depth of the return flow for each model. A basis-resolved inversion performs significantly better than one where the flow field is inverted for at each point in the radial grid. This is an encouraging result and might act as a guide in developing more realistic inversion strategies that can be applied to supergranular flows in the Sun.

In this paper, we detail the scientific objectives and outline a strawman payload of the SOLAR sail Investigation of the Sun (SOLARIS). The science objectives are to study the 3D structure of the solar magnetic and velocity field, the variation of total solar irradiance with latitude, and the structure of the corona. We show how we can meet these science objective using solar-sail technologies currently under development. We provide a tentative mission profile considering several trade-off approaches. We also provide a tentative mass budget breakdown and a perspective for a programmatic implementation.

We discuss current advances in forward and inverse modeling for local helioseismology. We report theoretical uniqueness results, in particular the Novikov-Agaltsov reconstruction algorithm, which is relevant to solving the non-linear inverse problem of time-distance helioseismology (finite amplitude pertubations to the medium). Numerical experiments were conducted to determine the number of frequencies required to reconstruct density and sound speed in the solar interior.

Deciphering the assembly history of the Milky Way is a formidable task, which becomes possible only if one can produce high-resolution chrono-chemo-kinematical maps of the Galaxy. Data from large-scale astrometric and spectroscopic surveys will soon provide us with a well-defined view of the current chemo-kinematical structure of the Milky Way, but will only enable a blurred view on the temporal sequence that led to the present-day Galaxy. As demonstrated by the (ongoing) exploitation of data from the pioneering photometric missions CoRoT, Kepler, and K2, asteroseismology provides the way forward: solar-like oscillating giants are excellent evolutionary clocks thanks to the availability of seismic constraints on their mass and to the tight age-initial-mass relation they adhere to. In this paper we identify five key outstanding questions relating to the formation and evolution of the Milky Way that will need precise and accurate ages for large samples of stars to be addressed, and we identify the requirements in terms of number of targets and the precision on the stellar properties that are needed to tackle such questions. By quantifying the asteroseismic yields expected from PLATO for red-giant stars, we demonstrate that these requirements are within the capabilities of the current instrument design, provided that observations are sufficiently long to identify the evolutionary state and allow robust and precise determination of acoustic-mode frequencies. This will allow us to harvest data of sufficient quality to reach a 10% precision in age. This is a fundamental pre-requisite to then reach the more ambitious goal of a similar level of accuracy, which will only be possible if we have to hand a careful appraisal of systematic uncertainties on age deriving from our limited understanding of stellar physics, a goal which conveniently falls within the main aims of PLATO's core science.

Rotational shear in Sun-like stars is thought to be an important ingredient in models of stellar dynamos. Thanks to helioseismology, rotation in the Sun is characterized well, but the interior rotation profiles of other Sun-like stars are not so well constrained. Until recently, measurements of rotation in Sun-like stars have focused on the mean rotation, but little progress has been made on measuring or even placing limits on differential rotation. Using asteroseismic measurements of rotation we aim to constrain the radial shear in five Sun-like stars observed by the NASA Kepler mission: KIC004914923, KIC005184732, KIC006116048, KIC006933899, and KIC010963065. We used stellar structure models for these five stars from previous works. These models provide the mass density, mode eigenfunctions, and the convection zone depth, which we used to compute the sensitivity kernels for the rotational frequency splitting of the modes. We used these kernels as weights in a parametric model of the stellar rotation profile of each star, where we allowed different rotation rates for the radiative interior and the convective envelope. This parametric model was incorporated into a fit to the oscillation power spectrum of each of the five Kepler stars. This fit included a prior on the rotation of the envelope, estimated from the rotation of surface magnetic activity measured from the photometric variability. The asteroseismic measurements without the application of priors are unable to place meaningful limits on the radial shear. Using a prior on the envelope rotation enables us to constrain the interior rotation rate and thus the radial shear. In the five cases that we studied, the interior rotation rate does not differ from the envelope by more than approximately +/-30%. Uncertainties in the rotational splittings are too large to unambiguously determine the sign of the radial shear.

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.

In this article, we derive and compute the sensitivity of measurements of coupling between normal modes of oscillation in the Sun to underlying flows. The theory is based on first-Born perturbation theory, and the analysis is carried out using the formalism described by \citet{lavely92}. Albeit tedious, we detail the derivation and compute the sensitivity of specific pairs of coupled normal modes to anomalies in the interior. Indeed, these kernels are critical for the accurate inference of convective flow amplitudes and large-scale circulations in the solar interior. We resolve some inconsistencies in the derivation of \citet{lavely92} and reformulate the fluid-continuity condition. We also derive and compute sound-speed kernels, paving the way for inverting for thermal anomalies alongside flows.

In recent years it has been claimed that the length of stellar activity cycles is determined by the stellar rotation rate. It is observed that the cycle period increases with rotation period along the so-called active and inactive sequences. In this picture the Sun occupies a solitary position in between the two sequences. Our goal is to measure cyclic variations of the stellar light curve amplitude and the rotation period using four years of Kepler data. Periodic changes of the light curve amplitude or the stellar rotation period are associated with an underlying activity cycle. Using the McQuillan et al. 2014 sample we compute the rotation period and the variability amplitude for each Kepler quarter and search for periodic variations of both time series. To test for periodicity in each stellar time series we consider Lomb-Scargle periodograms and use a selection based on a False Alarm Probability (FAP). We detect amplitude periodicities in 3203 stars between 0.5-6 years covering rotation periods between 1-40 days. Given our sample size of 23,601 stars and our selection criteria that the FAP is less than 5%, this number is almost three times higher than that expected from pure noise. We do not detect periodicities in the rotation period beyond those expected from noise. Our measurements reveal that the cycle period shows a weak dependence on rotation rate, slightly increasing for longer rotation period. We further show that the shape of the variability deviates from a pure sine curve, consistent with observations of the solar cycle. Our measurements do not support the existence of distinct sequences in the P_rot-P_cyc plane, although there is some evidence for the inactive sequence for rotation periods between 5-25 days. Unfortunately, the total observing time is too short to draw sound conclusions on activity cycles with similar length as the solar cycle.

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.

Convective flows are known as the prime means of transporting magnetic fields on the solar surface. Thus, small magnetic structures are good tracers of the turbulent flows. We study the migration and dispersal of magnetic bright features (MBFs) in intergranular areas observed at high spatial resolution with Sunrise/IMaX. We describe the flux dispersal of individual MBFs as a diffusion process whose parameters are computed for various areas in the quiet Sun and the vicinity of active regions from seeing-free data. We find that magnetic concentrations are best described as random walkers close to network areas (diffusion index, gamma=1.0), travelers with constant speeds over a supergranule (gamma=1.9-2.0), and decelerating movers in the vicinity of flux emergence and/or within active regions (gamma=1.4-1.5). The three types of regions host MBFs with mean diffusion coefficients of 130 km^2/s, 80-90 km^2/s, and 25-70 km^2/s, respectively. The MBFs in these three types of regions are found to display a distinct kinematic behavior at a confidence level in excess of 95%.

We present observations of transverse oscillations in slender Ca II H fibrils (SCFs) in the lower solar chromosphere. We use a 1 hr long time series of high- (spatial and temporal-) resolution seeing-free observations in a 0.11 nm wide passband covering the line core of Ca II H 396.9 nm from the second flight of the Sunrise balloon-borne solar observatory. The entire field of view, spanning the polarity inversion line of an active region close to the solar disk center, is covered with bright, thin, and very dynamic fine structures. Our analysis reveals the prevalence of transverse waves in SCFs with median amplitudes and periods on the order of 2.4+-0.8 km/s and 83+-29 s, respectively (with standard deviations given as uncertainties). We find that the transverse waves often propagate along (parts of) the SCFs with median phase speeds of 9+-14 km/s. While the propagation is only in one direction along the axis in some of the SCFs, propagating waves in both directions, as well as standing waves are also observed. The transverse oscillations are likely Alfvenic and are thought to be representative of magnetohydrodynamic kink waves. The wave propagation suggests that the rapid high-frequency transverse waves, often produced in the lower photosphere, can penetrate into the chromosphere with an estimated energy flux of ~ 15 kW/m^2. Characteristics of these waves differ from those reported for other fibrillar structures, which, however, were observed mainly in the upper solar chromosphere.

A dense forest of slender bright fibrils near a small solar active region is seen in high-quality narrowband Ca II H images from the SuFI instrument onboard the Sunrise balloon-borne solar observatory. The orientation of these slender Ca II H fibrils (SCF) overlaps with the magnetic field configuration in the low solar chromosphere derived by magnetostatic extrapolation of the photospheric field observed with Sunrise/IMaX and SDO/HMI. In addition, many observed SCFs are qualitatively aligned with small-scale loops computed from a novel inversion approach based on best-fit numerical MHD simulation. Such loops are organized in canopy-like arches over quiet areas that differ in height depending on the field strength near their roots.

How and where are coronal loops rooted in the solar lower atmosphere? The details of the magnetic environment and its evolution at the footpoints of coronal loops are crucial to understanding the processes of mass and energy supply to the solar corona. To address the above question, we use high-resolution line-of-sight magnetic field data from the Imaging Magnetograph eXperiment instrument on the SUNRISE balloon-borne observatory and coronal observations from the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory of an emerging active region. We find that the coronal loops are often rooted at the locations with minor small-scale but persistent opposite-polarity magnetic elements very close to the larger dominant polarity. These opposite-polarity small-scale elements continually interact with the dominant polarity underlying the coronal loop through flux cancellation. At these locations we detect small inverse Y-shaped jets in chromospheric Ca II H images obtained from the SUNRISE Filter Imager during the flux cancellation. Our results indicate that magnetic flux cancellation and reconnection at the base of coronal loops due to mixed polarity fields might be a crucial feature for the supply of mass and energy into the corona.

Accurate modelling of solar-like oscillators requires that modelled mode frequencies are corrected for the systematic shift caused by improper modelling of the near-surface layers, known as the surface effect. ... We investigate how much additional uncertainty is introduced to stellar model parameters by our uncertainty about the functional form of the surface effect. At the same time, we test whether any of the parametrizations is significantly better or worse at modelling observed subgiants and low-luminosity red giants. We model six stars observed by Kepler that show clear mixed modes. We fix the input physics of the stellar models and vary the choice of surface correction ... Models using a solar-calibrated power law correction consistently fit the observations more poorly than the other four corrections. Models with the remaining four corrections generally fit ... about equally well, with the combined surface correction by Ball & Gizon perhaps being marginally superior. The fits broadly agree on the model parameters within about the $2\sigma$ uncertainties, with discrepancies between the modified Lorentzian and free power law corrections occasionally exceeding the $3\sigma$ level. Relative to the best-fitting values, the total uncertainties on the masses, radii and ages of the stars are all less than 2, 1 and 6 per cent, respectively. A solar-calibrated power law ... appears unsuitable for use with more evolved solar-like oscillators. Among the remaining surface corrections, the uncertainty in the model parameters introduced by the surface effects is about twice as large as the uncertainty in the individual fits for these six stars. Though the fits are thus somewhat less certain because of our uncertainty of how to manage the surface effect, these results also demonstrate that it is feasible to model the individual mode frequencies of subgiants and low-luminosity red giants. ...

Magneto-static models may overcome some of the issues facing force-free magnetic field extrapolations. So far they have seen limited use and have faced problems when applied to quiet-Sun data. Here we present a first application to an active region. We use solar vector magnetic field measurements gathered by the IMaX polarimeter during the flight of the \sunrise{} balloon-borne solar observatory in June 2013 as boundary condition for a magneto-static model of the higher solar atmosphere above an active region. The IMaX data are embedded in active region vector magnetograms observed with SDO/HMI. This work continues our magneto-static extrapolation approach, which has been applied earlier ({\it Paper I}) to a quiet Sun region observed with \sunrise{} I. In an active region the signal-to-noise-ratio in the measured Stokes parameters is considerably higher than in the quiet Sun and consequently the IMaX measurements of the horizontal photospheric magnetic field allow us to specify the free parameters of the model in a special class of linear magneto-static equilibria. The high spatial resolution of IMaX (110-130 km, pixel size 40 km) enables us to model the non-force-free layer between the photosphere and the mid chromosphere vertically by about 50 grid points. In our approach we can incorporate some aspects of the mixed beta layer of photosphere and chromosphere, e.g., taking a finite Lorentz force into account, which was not possible with lower resolution photospheric measurements in the past. The linear model does not, however, permit to model intrinsic nonlinear structures like strongly localized electric currents.

We report the detection of oscillations in slender Ca II H fibrils (SCFs) from high-resolution observations acquired with the Sunrise balloon-borne solar observatory. The SCFs show obvious oscillations in their intensity, but also their width. The oscillatory behaviors are investigated at several positions along the axes of the SCFs. A large majority of fibrils show signs of oscillations in intensity. Their periods and phase speeds are analyzed using a wavelet analysis. The width and intensity perturbations have overlapping distributions of the wave period. The obtained distributions have median values of the period of $32\pm17$s and $36\pm25$s, respectively. We find that the fluctuations of both parameters propagate in the SCFs with speeds of ${11}^{+49}_{-11}$ km/s and ${15}^{+34}_{-15}$ km/s, respectively. Furthermore, the width and intensity oscillations have a strong tendency to be either in anti-phase, or, to a smaller extent, in phase. This suggests that the oscillations of both parameters are caused by the same wave mode and that the waves are likely propagating. Taking all the evidence together, the most likely wave mode to explain all measurements and criteria is the fast sausage mode.

The SUNRISE balloon-borne solar observatory, consisting of a 1~m aperture telescope that provided a stabilized image to a UV filter imager and an imaging vector polarimeter, carried out its second science flight in June 2013. It provided observations of parts of active regions at high spatial resolution, including the first high-resolution images in the Mg~{\sc ii}~k line. The obtained data are of very high quality, with the best UV images reaching the diffraction limit of the telescope at 3000~\AA\ after Multi-Frame Blind Deconvolution reconstruction accounting for phase-diversity information. Here a brief update is given of the instruments and the data reduction techniques, which includes an inversion of the polarimetric data. Mainly those aspects that evolved compared with the first flight are described. A tabular overview of the observations is given. In addition, an example time series of a part of the emerging active region NOAA AR~11768 observed relatively close to disk centre is described and discussed in some detail. The observations cover the pores in the trailing polarity of the active region, as well as the polarity inversion line where flux emergence was ongoing and a small flare-like brightening occurred in the course of the time series. The pores are found to contain magnetic field strengths ranging up to 2500~G and, while large pores are clearly darker and cooler than the quiet Sun in all layers of the photosphere, the temperature and brightness of small pores approach or even exceed those of the quiet Sun in the upper photosphere.

Context. In time-distance helioseismology, wave travel times are measured from the two-point cross-covariance function of solar oscillations and are used to image the solar convection zone in three dimensions. There is, however, also information in the amplitude of the cross-covariance function, for example about seismic wave attenuation. Aims. Here we develop a convenient procedure to measure the amplitude of the cross-covariance function of solar oscillations. Methods. In this procedure, the amplitude of the cross-covariance function is linearly related to the cross-covariance function and can be measured even for high levels of noise. Results. As an example application, we measure the amplitude perturbations of the seismic waves that propagate through the sunspot in active region NOAA 9787. We can recover the amplitude variations due to the scattering and attenuation of the waves by the sunspot and associated finite-wavelength effects. Conclusions. The proposed definition of cross-covariance amplitude is robust to noise, can be used to relate measured amplitudes to 3D perturbations in the solar interior under the Born approximation, and will provide independent information from the travel times.

We use seeing-free high spatial resolution Ca II H data obtained by the SUNRISE observatory to determine properties of slender fibrils in the lower solar chromosphere. In this work we use intensity images taken with the SUFI instrument in the Ca II H line during the second scientific flight of the SUNRISE observatory to identify and track elongated bright structures. After the identification, we analyze theses structures in order to extract their morphological properties. We identify 598 slender Ca II H fibrils (SCFs) with an average width of around 180 km, a length between 500 km and 4000 km, an average lifetime of ~400 s, and an average curvature of 0.002 arcsec^-1. The maximum lifetime of the SCFs within our time series of 57 minutes is ~2000 s. We discuss similarities and differences of the SCFs with other small-scale, chromospheric structures such as spicules of type I and II, or Ca II K fibrils.

We study the dynamics and topology of an emerging magnetic flux concentration using high spatial resolution spectropolarimetric data acquired with the Imaging Magnetograph eXperiment on board the Sunrise balloon-borne solar observatory. We obtain the full vector magnetic field and the line-of-sight (LOS) velocity through inversions of the Fe I line at 525.02 nm with the SPINOR code. The derived vector magnetic field is used to trace magnetic field lines. Two magnetic flux concentrations with different polarity and LOS velocities are found to be connected by a group of arch-shaped magnetic field lines. The positive polarity footpoint is weaker (1100 G) and displays an upflow, while the negative polarity footpoint is stronger (2200 G) and shows a downflow. This configuration is naturally interpreted as a siphon flow along an arched magnetic flux tube.

Stars are not perfectly spherically symmetric. They are deformed by rotation and magnetic fields. Until now, the study of stellar shapes has only been possible with optical interferometry for a few of the fastest-rotating nearby stars. We report an asteroseismic measurement, with much better precision than interferometry, of the asphericity of an A-type star with a rotation period of 100 days. Using the fact that different modes of oscillation probe different stellar latitudes, we infer a tiny but significant flattening of the star's shape of $\Delta R/R = (1.8 \pm 0.6) \times 10^{-6}$. For a stellar radius $R$ that is $2.24$ times the solar radius, the difference in radius between the equator and the poles is $\Delta R = 3 \pm 1$ km. Because the observed $\Delta R/R$ is only one-third of the expected rotational oblateness, we conjecture the presence of a weak magnetic field on a star that does not have an extended convective envelope. This calls to question the origin of the magnetic field.

We present a new method of Stokes inversion of spectropolarimetric data and evaluate it by taking the example of a SUNRISE/IMaX observation. An archive of synthetic Stokes profiles is obtained by the spectral synthesis of state-of-the-art magnetohydrodynamics (MHD) simulations and a realistic degradation to the level of the observed data. The definition of a merit function allows the archive to be searched for the synthetic Stokes profiles that match the observed profiles best. In contrast to traditional Stokes inversion codes, which solve the Unno-Rachkovsky equations for the polarized radiative transfer numerically and fit the Stokes profiles iteratively, the new technique provides the full set of atmospheric parameters. This gives us the ability to start an MHD simulation that takes the inversion result as initial condition. After a relaxation process of half an hour solar time we obtain physically consistent MHD data sets with a target similar to the observation. The new MHD simulation is used to repeat the method in a second iteration, which further improves the match between observation and simulation, resulting in a factor of 2.2 lower mean $\chi^2$ value. One advantage of the new technique is that it provides the physical parameters on a geometrical height scale. It constitutes a first step towards inversions giving results consistent with the MHD equations.

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.

In June 2013, the two scientific instruments onboard the second Sunrise mission witnessed, in detail, a small-scale magnetic flux emergence event as part of the birth of an active region. The Imaging Magnetograph Experiment (IMaX) recorded two small (~5 arcsec) emerging flux patches in the polarized filtergrams of a photospheric Fe I spectral line. Meanwhile, the Sunrise Filter Imager (SuFI) captured the highly dynamic chromospheric response to the magnetic fields pushing their way through the lower solar atmosphere. The serendipitous capture of this event offers a closer look at the inner workings of active region emergence sites. In particular, it reveals in meticulous detail how the rising magnetic fields interact with the granulation as they push through the Sun's surface, dragging photospheric plasma in their upward travel. The plasma that is burdening the rising field slides along the field lines, creating fast downflowing channels at the footpoints. The weight of this material anchors this field to the surface at semi-regular spatial intervals, shaping it in an undulatory fashion. Finally, magnetic reconnection enables the field to release itself from its photospheric anchors, allowing it to continue its voyage up to higher layers. This process releases energy that lights up the arch-filament systems and heats the surrounding chromosphere.

Binary star systems are important for understanding stellar structure and evolution, and are especially useful when oscillations can be detected and analysed with asteroseismology. However, only four systems are known in which solar-like oscillations are detected in both components. Here, we analyse the fifth such system, HD 176465, which was observed by Kepler. We carefully analysed the system's power spectrum to measure individual mode frequencies, adapting our methods where necessary to accommodate the fact that both stars oscillate in a similar frequency range. We also modelled the two stars independently by fitting stellar models to the frequencies and complementary parameters. We are able to cleanly separate the oscillation modes in both systems. The stellar models produce compatible ages and initial compositions for the stars, as is expected from their common and contemporaneous origin. Combining the individual ages, the system is about 3.0$\pm$0.5 Gyr old. The two components of HD 176465 are young physically-similar oscillating solar analogues, the first such system to be found, and provide important constraints for stellar evolution and asteroseismology.

While the velocity fluctuations of supergranulation dominate the spectrum of solar convection at the solar surface, very little is known about the fluctuations in other physical quantities like temperature or density at supergranulation scale. Using SDO/HMI observations, we characterize the intensity contrast of solar supergranulation at the solar surface. We identify the positions of ${\sim}10^4$ outflow and inflow regions at supergranulation scales, from which we construct average flow maps and co-aligned intensity and magnetic field maps. In the average outflow center, the maximum intensity contrast is $(7.8\pm0.6)\times10^{-4}$ (there is no corresponding feature in the line-of-sight magnetic field). This corresponds to a temperature perturbation of about $1.1\pm0.1$ K, in agreement with previous studies. We discover an east-west anisotropy, with a slightly deeper intensity minimum east of the outflow center. The evolution is asymmetric in time: the intensity excess is larger 8 hours before the reference time (the time of maximum outflow), while it has almost disappeared 8 hours after the reference time. In the average inflow region, the intensity contrast mostly follows the magnetic field distribution, except for an east-west anisotropic component that dominates 8 hours before the reference time. We suggest that the east-west anisotropy in the intensity is related to the wave-like properties of supergranulation.

Spectropolarimetric observations from Sunrise II/IMaX obtained in June 2013 are used for a statistical analysis to determine the physical properties of moving magnetic features (MMFs) observed near a pore. MMFs of the same and opposite polarity with respect to the pore are found to stream from its border at an average speed of 1.3 km s$^{-1}$ and 1.2 km s$^{-1}$ respectively, with mainly same-polarity MMFs found further away from the pore. MMFs of both polarities are found to harbor rather weak, inclined magnetic fields. Opposite-polarity MMFs are blue-shifted, while same-polarity MMFs do not show any preference for up- or downflows. Most of the MMFs are found to be of sub-arcsecond size and carry a mean flux of $\sim$ 1.2$\times 10^{17}$ Mx.

Ellerman Bombs are signatures of magnetic reconnection, which is an important physical process in the solar atmosphere. How and where they occur is a subject of debate. In this paper we analyse Sunrise/IMaX data together with 3D MHD simulations that aim to reproduce the exact scenario proposed for the formation of these features. Although the observed event seems to be more dynamic and violent than the simulated one, simulations clearly confirm the basic scenario for the production of EBs. The simulations also reveal the full complexity of the underlying process. The simulated observations show that the Fe I 525.02 nm line gives no information on the height where reconnection takes place. It can only give clues about the heating in the aftermath of the reconnection. The information on the magnetic field vector and velocity at this spatial resolution is, however, extremely valuable because it shows what numerical models miss and how they can be improved.

Most future surveys designed to discover transiting exoplanets, including TESS and PLATO, will target bright (V<13) and nearby solar-type stars having a spectral type later than F5. In order to enhance the probability of identifying transits, these surveys must cover a very large area on the sky, because of the intrinsically low areal density of bright targets. Unfortunately, no existing catalog of stellar parameters is both deep and wide enough to provide a homogeneous input list. As the first Gaia data release exploitable for this purpose is expected to be released not earlier than late 2017, we have devised an improved reduced-proper-motion method to discriminate late field dwarfs and giants by combining UCAC4 proper motions with APASS DR6 photometry, and relying on RAVE DR4 as an external calibrator. The output, named UCAC4-RPM, is a publicly-available, complete all-sky catalog of solar-type dwarfs down to V<13.5, plus an extension to logg>3.0 subgiants. The relatively low amount of contamination (defined as the fraction of false positives; <30%) also makes UCAC4-RPM a useful tool for the past and ongoing ground-based transit surveys, which need to discard candidate signals originating from early-type or giant stars. As an application, we show how UCAC4-RPM may support the preparation of the TESS (that will map almost the entire sky) input catalog and the input catalog of PLATO, planned to survey more than half of the whole sky with exquisite photometric precision.

Observations from the Solar Dynamics Observatory (SDO) have the potential for allowing the helioseismic study of the formation of hundreds of active regions, which would enable us to perform statistical analyses. Our goal is to collate a uniform data set of emerging active regions observed by the SDO/HMI instrument suitable for helioseismic analysis up to seven days before emergence. We restricted the sample to active regions that were visible in the continuum and emerged into quiet Sun largely avoiding pre-existing magnetic regions. As a reference data set we paired a control region (CR), with the same latitude and distance from central meridian, with each emerging active region (EAR). We call this data set, which is currently comprised of 105 emerging active regions observed between May 2010 and November 2012, the SDO Helioseismic Emerging Active Region (SDO/HEAR) survey. To demonstrate the utility of a data set of a large number of emerging active regions, we measure the relative east-west velocity of the leading and trailing polarities from the line-of-sight magnetogram maps during the first day after emergence. The latitudinally averaged line-of-sight magnetic field of all the EARs shows that, on average, the leading (trailing) polarity moves in a prograde (retrograde) direction with a speed of 121 +/- 22 m/s (-70 +/- 13 m/s) relative to the Carrington rotation rate in the first day. However, relative to the differential rotation of the surface plasma, the east-west velocity is symmetric, with a mean of 95 +/- 13 m/s. The SDO/HEAR data set will not only be useful for helioseismic studies, but will also be useful to study other features such as the surface magnetic field evolution of a large sample of EARs.

Helioseismology has revealed that the angular velocity of the Sun increases with depth in the outermost 35 Mm of the Sun. Recently, we have shown that the logarithmic radial gradient ($\rm d\ln\Omega/\rm d\ln r $) in the upper 10~Mm is close to $-1$ from the equator to $60^\circ$ latitude.We aim to measure the temporal variation of the rotational shear over solar cycle 23 and the rising phase of cycle 24 (1996-2015). We used f mode frequency splitting data spanning 1996 to 2011 from the Michelson Doppler Imager (MDI) and 2010 to 2015 from the Helioseismic Magnetic Imager (HMI). In a first for such studies, the f mode frequency splitting data were obtained from 360-day time series. We used the same method as in our previous work for measuring $\rm d\ln\Omega/d\ln r $ from the equator to $80^\circ$ latitude in the outer 13~Mm of the Sun. Then, we calculated the variation of the gradient at annual cadence relative to the average over 1996 to 2015. We found the rotational shear at low latitudes ($0^\circ$ to $30^\circ$) to vary in-phase with the solar activity, varying by $\sim \pm 10$\% over the period 1996 to 2015. At high latitudes ($60^\circ$ to $80^\circ$), we found rotational shear to vary in anti-phase with the solar activity. By comparing the radial gradient obtained from the splittings of the 360-day and the corresponding 72-day time series of HMI and MDI data, we suggest that the splittings obtained from the 72-day HMI time series suffer from systematic errors. We provide a quantitative measurement of the temporal variation of the outer part of the near surface shear layer which may provide useful constraints on dynamo models and differential rotation theory.

Magnetic field emerges at the surface of the Sun as sunspots and active regions. This process generates a poloidal magnetic field from a rising toroidal flux tube, it is a crucial but poorly understood aspect of the solar dynamo. The emergence of magnetic field is also important because it is a key driver of solar activity. We show that measurements of horizontal flows at the solar surface around emerging active regions, in combination with numerical simulations of solar magnetoconvection, can constrain the subsurface rise speed of emerging magnetic flux. The observed flows imply that the rise speed of the magnetic field is no larger than 150 m/s at a depth of 20 Mm, that is, well below the prediction of the (standard) thin flux tube model but in the range expected for convective velocities at this depth. We conclude that convective flows control the dynamics of rising flux tubes in the upper layers of the Sun and cannot be neglected in models of flux emergence.

Context: The cross-covariance of solar oscillations observed at pairs of points on the solar surface is a fundamental ingredient in time-distance helioseismology. Wave travel times are extracted from the cross-covariance function and are used to infer the physical conditions in the solar interior. Aims: Understanding the statistics of the two-point cross-covariance function is a necessary step towards optimizing the measurement of travel times. Methods: By modeling stochastic solar oscillations, we evaluate the variance of the cross-covariance function as function of time-lag and distance between the two points. Results: We show that the variance of the cross-covariance is independent of both time-lag and distance in the far field, i.e., when they are large compared to the coherence scales of the solar oscillations. Conclusions: The constant noise level for the cross-covariance means that the signal-to-noise ratio for the cross-covariance is proportional to the amplitude of the expectation value of the cross-covariance. This observation is important for planning data analysis efforts.

... [C]urrent stellar models predict oscillation frequencies that are systematically affected by simplified modelling of the near-surface layers. We use three-dimensional radiation hydrodynamics simulations to better model the near-surface equilibrium structure of dwarfs with spectral types F3, G2, K0 and K5, and examine the differences between oscillation mode frequencies. ... We precisely match stellar models to the simulations' gravities and effective temperatures at the surface, and to the temporally- and horizontally-averaged densities and pressures at their deepest points. We then replace the near-surface structure with that of the averaged simulation and compute the change in the oscillation mode frequencies. We also fit the differences using several parametric models currently available in the literature. The surface effect in the stars of solar-type and later is qualitatively similar and changes steadily with decreasing effective temperature. In particular, the point of greatest frequency difference decreases slightly as a fraction of the acoustic cut-off frequency and the overall scale of the surface effect decreases. The surface effect in the hot, F3-type star follows the same trend in scale (i.e. it is larger in magnitude) but shows a different overall variation with mode frequency. We find that the two-term fit by Ball & Gizon (2014) is best able to reproduce the surface terms across all four spectral types, although the scaled solar term and a modified Lorentzian function also match the three cooler simulations reasonably well. ... Our simplified results suggest that the surface effect is generally larger in hotter stars (and correspondingly smaller in cooler stars) and of similar shape in stars of solar type and cooler. However, we cannot presently predict whether this will remain so when other components of the surface effect are included.

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 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.

High-quality astrometric, spectroscopic, interferometric and, importantly, asteroseismic observations are available for $\alpha$ Cen A, which is the closest binary star system to earth. Taking all these constraints into account, we study the internal structure of the star by means of theoretical modelling. Using the Aarhus STellar Evolution Code (ASTEC) and the tools of Computational Bayesian Statistics, in particular a Markov chain Monte Carlo algorithm, we perform statistical inferences for the physical characteristics of the star. We find that $\alpha$ Cen A has a probability of approximately 40% of having a convective core. This probability drops to few percents if one considers reduced rates for the $^{14}$N(p,$\gamma$)$^{15}$O reaction. These convective cores have fractional radii less than 8% when overshoot is neglected. Including overshooting also leads to the possibility of a convective core mostly sustained by the ppII chain energy output. We finally show that roughly 30% of the stellar models describing $\alpha$ Cen A are in the subgiant regime.

Radial differential rotation is an important parameter for stellar dynamo theory and for understanding angular momentum transport. We investigate the potential of using a large number of similar stars simultaneously to constrain their average radial differential rotation gradient: we call this 'ensemble fitting'. We use a range of stellar models along the main sequence, each with a synthetic rotation profile. The rotation profiles are step functions with a step of -350 nHz, which is located at the base of the convection zone. These models are used to compute the rotational splittings of the p modes and to model their uncertainties. We then fit an ensemble of stars to infer the average step size. All the uncertainties on the inferred step size for individual stars are of the order 1 micro Hz. Using 15 stellar models in an ensemble fit, we show that the uncertainty on the average step size is reduced to less than the input step size, which allows us to constrain the sign of the radial differential rotation. We show that a solar-like step size (approximately 30 nHz) can be constrained by an ensemble fit of thousands of main-sequence stars. Observing the number of stars required to successfully exploit the ensemble fitting method will be possible with future asteroseismology missions, such as PLATO. We demonstrate the potential of ensemble fitting by showing that any systematic differences in the average step size between F, G, and K-type stars larger than 100 nHz can be detected.

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.

In previous work we identified six Sun-like stars observed by Kepler with exceptionally clear asteroseismic signatures of rotation. Here, we show that five of these stars exhibit surface variability suitable for measuring rotation. In order to further constrain differential rotation, we compare the rotation periods obtained from light-curve variability with those from asteroseismology. The two rotation measurement methods are found to agree within uncertainties, suggesting that radial differential rotation is weak, as is the case for the Sun. Furthermore, we find significant discrepancies between ages from asteroseismology and from three different gyrochronology relations, implying that stellar age estimation is problematic even for Sun-like stars.

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.

The high-precision photometry from the CoRoT and Kepler satellites has led to measurements of surface rotation periods for tens of thousands of stars. Our main goal is to derive ages of thousands of field stars using consistent rotation period measurements in different gyrochronology relations. Multiple rotation periods are interpreted as surface differential rotation (DR). We re-analyze the sample of 24,124 Kepler stars from Reinhold et al. (2013) using different approaches based on the Lomb-Scargle periodogram. Each quarter (Q1-Q14) is treated individually using a prewhitening approach. Additionally, the full time series, and different segments thereof are analyzed. For more than 18,500 stars our results are consistent with the rotation periods from McQuillan et al. (2014). Thereof, more than 12,300 stars show multiple significant peaks, which we interpret as DR. Gyrochronology ages between 100 Myr and 10 Gyr were derived for more than 17,000 stars using different gyrochronology relations. We find a bimodal age distribution for Teff between 3200-4700 K. The derived ages reveal an empirical activity-age relation using photometric variability as stellar activity proxy. Additionally, we found 1079 stars with extremely stable (mostly short) periods. Half of these periods may be associated with rotation stabilized by non-eclipsing companions, the other half might be due to pulsations. The derived gyrochronology ages are well constrained since more than 93.0 % of the stars seem to be younger than the Sun where calibration is most reliable. Explaining the bimodality in the age distribution is challenging, and limits accurate stellar age predictions. The existence of cool stars with almost constant rotation period over more than three years of observation might be explained by synchronization with stellar companions, or a dynamo mechanism keeping the spot configurations extremely stable.

We present the first detections by the NASA K2 Mission of oscillations in solar-type stars, using short-cadence data collected during K2 Campaign\,1 (C1). We understand the asteroseismic detection thresholds for C1-like levels of photometric performance, and we can detect oscillations in subgiants having dominant oscillation frequencies around $1000\,\rm \mu Hz$. Changes to the operation of the fine-guidance sensors are expected to give significant improvements in the high-frequency performance from C3 onwards. A reduction in the excess high-frequency noise by a factor of two-and-a-half in amplitude would bring main-sequence stars with dominant oscillation frequencies as high as ${\simeq 2500}\,\rm \mu Hz$ into play as potential asteroseismic targets for K2.

Flow vorticity is a fundamental property of turbulent convection in rotating systems. Solar supergranules exhibit a preferred sense of rotation, which depends on the hemisphere. This is due to the Coriolis force acting on the diverging horizontal flows. We aim to spatially resolve the vertical flow vorticity of the average supergranule at different latitudes, both for outflow and inflow regions. To measure the vertical vorticity, we use two independent techniques: time-distance helioseismology (TD) and local correlation tracking of granules in intensity images (LCT) using data from the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). Both maps are corrected for center-to-limb systematic errors. We find that 8-h TD and LCT maps of vertical vorticity are highly correlated at large spatial scales. Associated with the average supergranule outflow, we find tangential (vortical) flows that reach about 10 m/s in the clockwise direction at 40{\deg} latitude. In average inflow regions, the tangential flow reaches the same magnitude, but in the anti-clockwise direction. These tangential velocities are much smaller than the radial (diverging) flow component (300 m/s for the average outflow and 200 m/s for the average inflow). The results for TD and LCT as measured from HMI are in excellent agreement for latitudes between $-$60{\deg} and 60{\deg}. From HMI LCT, we measure the vorticity peak of the average supergranule to have a full width at half maximum of about 13 Mm for outflows and 8 Mm for inflows. This is larger than the spatial resolution of the LCT measurements (about 3 Mm). On the other hand, the vorticity peak in outflows is about half the value measured at inflows (e.g. 4/(10^6 s) clockwise compared to 8/(10^6 s) anti-clockwise at 40{\deg} latitude). Results from MDI/SOHO obtained in 2010 are biased compared to the HMI/SDO results for the same period.

Supergranules in the quiet Sun are outlined by a web-like structure of enhanced magnetic field strength, the so-called magnetic network. We aim to map the magnetic network field around the average supergranule near disk center. We use observations of the line-of-sight component of the magnetic field from the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). The average supergranule is constructed by coaligning and averaging over 3000 individual supergranules. We determine the positions of the supergranules with an image segmentation algorithm that we apply on maps of the horizontal flow divergence measured using time-distance helioseismology. In the center of the average supergranule the magnetic (intranetwork) field is weaker by about 2.2 Gauss than the background value (3.5 Gauss), whereas it is enhanced in the surrounding ring of horizontal inflows (by about 0.6 Gauss on average). We find that this network field is significantly stronger west (prograde) of the average supergranule than in the east (by about 0.3 Gauss). With time-distance helioseismology, we find a similar anisotropy. The observed anisotropy of the magnetic field adds to the mysterious dynamical properties of solar supergranulation.

Acoustic modes of oscillation are affected by stellar activity, however it is unclear how starspots contribute to these changes. Here we investigate the nonmagnetic effects of starspots on global modes with angular degree $\ell \leq 2$ in highly active stars, and characterize the spot seismic signature on synthetic light curves. We perform 3D time-domain simulations of linear acoustic waves to study their interaction with a model starspot. We model the spot as a 3D change in the sound speed stratification with respect to a convectively stable stellar background, built from solar Model S. We perform a parametric study by considering different depths and perturbation amplitudes. Exact numerical simulations allow the investigation of the wavefield-spot interaction beyond first order perturbation theory. The interaction of the axisymmetric modes with the starspot is strongly nonlinear. As mode frequency increases, the frequency shifts for radial modes exceed the value predicted by linear theory, while the shifts for the $\ell=2, m=0$ modes are smaller than predicted by linear theory, with avoided-crossing-like patterns forming between the $m=0$ and $m=1$ mode frequencies. The nonlinear behavior increases with increasing spot amplitude and/or decreasing depth. Linear theory still reproduces the correct shifts for nonaxisymmetric modes. In the nonlinear regime the mode eigenfunctions are not pure spherical harmonics, but rather a mixture of different spherical harmonics. This mode mixing, together with the frequency changes, may lead to misidentification of the modes in the observed acoustic power spectra.

We measure rotation periods for 12151 stars in the Kepler field, based on the photometric variability caused by stellar activity. Our analysis returns stable rotation periods over at least six out of eight quarters of Kepler data. This large sample of stars enables us to study the rotation periods as a function of spectral type. We find good agreement with previous studies and vsini measurements for F, G and K stars. Combining rotation periods, B-V color, and gyrochronology relations, we find that the cool stars in our sample are predominantly younger than ~1Gyr.

Our Sun, primarily composed of ionized hydrogen and helium, has a surface temperature of 5777~K and a radius $R_\odot \approx 696,000$ km. In the outer $R_\odot/3$, energy transport is accomplished primarily by convection. Using typical convective velocities $u\sim100\,\rm{m\,s^{-1}}$ and kinematic viscosities of order $10^{-4}$ m$^{2}$s$^{-1}$, we obtain a Reynolds number $Re \sim 10^{14}$. Convection is thus turbulent, causing a vast range of scales to be excited. The Prandtl number, $Pr$, of the convecting fluid is very low, of order $10^{-7}$\,--\,$10^{-4}$, so that the Rayleigh number ($\sim Re^2 Pr$) is on the order of $10^{21}\,-\,10^{24}$. Solar convection thus lies in extraordinary regime of dynamical parameters, highly untypical of fluid flows on Earth. Convective processes in the Sun drive global fluid circulations and magnetic fields, which in turn affect its visible outer layers ("solar activity") and, more broadly, the heliosphere ("space weather"). The precise determination of the depth of solar convection zone, departures from adiabaticity of the temperature gradient, and the internal rotation rate as a function of latitude and depth are among the seminal contributions of helioseismology towards understanding convection in the Sun. Contemporary helioseismology, which is focused on inferring the properties of three-dimensional convective features, suggests that transport velocities are substantially smaller than theoretical predictions. Furthermore, helioseismology provides important constraints on the anisotropic Reynolds stresses that control the global dynamics of the solar convection zone. This review discusses the state of our understanding of convection in the Sun, with a focus on helioseismic diagnostics. We present our considerations with the interests of fluid dynamicists in mind.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

Time-distance helioseismology provides information about vector flows in the near-surface layers of the Sun by measuring wave travel times between points on the solar surface. Specific spatial averages of travel times have been proposed for distinguishing between flows in the east-west and north-south directions and measuring the horizontal divergence of the flows. No specific measurement technique has, however, been developed to measure flow vorticity. Here we propose a new measurement technique tailored to measuring the vertical component of vorticity. Fluid vorticity is a fundamental property of solar convection zone dynamics and of rotating turbulent convection in particular. The method consists of measuring the travel time of waves along a closed contour on the solar surface in order to approximate the circulation of the flow along this contour. Vertical vorticity is related to the difference between clockwise and counter-clockwise travel times. We applied the method to characterize the vortical motions of solar convection using helioseismic data from the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory (SDO/HMI) and from the Michelson Doppler Imager onboard the Solar and Heliospheric Observatory (SOHO/MDI). Away from the equator, a clear correlation between vertical vorticity and horizontal divergence is detected. Horizontal outflows are associated with negative vorticity in the northern hemisphere and positive vorticity in the southern hemisphere. The signal is much stronger for HMI than for MDI observations. We characterize the spatial power spectrum of the signal by comparison with a noise model. Vertical vorticity at horizontal wavenumbers below 250/R_Sun can be probed with this helioseismic technique.

Context. Several upcoming helioseismology space missions are very limited in telemetry and will have to perform extensive data compression. This requires the development of new methods of data compression. Aims. We give an overview of the influence of lossy data compression on local helioseismology. We investigate the effects of several lossy compression methods (quantization, JPEG compression, and smoothing and subsampling) on power spectra and time-distance measurements of supergranulation flows at disk center. Methods. We applied different compression methods to tracked and remapped Dopplergrams obtained by the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory. We determined the signal-to-noise ratio of the travel times computed from the compressed data as a function of the compression efficiency. Results. The basic helioseismic measurements that we consider are very robust to lossy data compression. Even if only the sign of the velocity is used, time-distance helioseismology is still possible. We achieve the best results by applying JPEG compression on spatially subsampled data. However, our conclusions are only valid for supergranulation flows at disk center and may not be valid for all helioseismology applications.

Helioseismology has provided unprecedented information about the internal rotation of the Sun. One of the important achievements was the discovery of two radial shear layers: one near the bottom of the convection zone (the tachocline) and one near the surface. These shear layers may be important ingredients for explaining the magnetic cycle of the Sun. We measure the logarithmic radial gradient of the rotation rate ($\rm d\ln\Omega/\rm d\ln r$) near the surface of the Sun using 15 years of f mode rotational frequency splittings from the Michelson Doppler Imager (MDI) and four years of data from the Helioseismic and Magnetic Imager (HMI). We model the angular velocity of the Sun in the upper $\sim 10$ Mm as changing linearly with depth and use a multiplicative optimally localized averaging inversion to infer the gradient of the rotation rate as a function of latitude. Both the MDI and HMI data show that $\rm d\ln\Omega/\rm d\ln r$ is close to $-1$ from the equator to 60$^{\circ}$ latitude and stays negative up to 75$^{\circ}$ latitude. However, the value of the gradient is different for MDI and HMI for latitudes above $60^{\circ}$. Additionally, there is a significant difference between the value of $d\ln\Omega/d\ln r$ using an older and recently reprocessed MDI data for latitudes above $30^\circ$. We could reliably infer the value of $\rm d\ln\Omega/\rm d\ln r$ up to 60$^{\circ}$, but not above this latitude, which will hopefully constrain theories of the near-surface shear layer and dynamo. Furthermore, the recently reprocessed MDI splitting data are more reliable than the older versions which contained clear systematic errors in the high degree f modes.

Asteroseismology offers the prospect of constraining differential rotation in Sun-like stars. Here we have identified six high signal-to-noise main-sequence Sun-like stars in the Kepler field, which all have visible signs of rotational splitting of their p-mode frequencies. For each star, we extract the rotational frequency splitting and inclination angle from separate mode sets (adjacent modes with l=2, 0, and 1) spanning the p-mode envelope. We use a Markov chain Monte Carlo method to obtain the best fit and errors associated with each parameter. We are able to make independent measurements of rotational splittings of ~8 radial orders for each star. For all six stars, the measured splittings are consistent with uniform rotation, allowing us to exclude large radial differential rotation. This work opens the possibility of constraining internal rotation of Sun-like stars.

After emerging to the solar surface, the Sun's magnetic field displays a complex and intricate evolution. The evolution of the surface field is important for several reasons. One is that the surface field, and its dynamics, sets the boundary condition for the coronal and heliospheric magnetic fields. Another is that the surface evolution gives us insight into the dynamo process. In particular, it plays an essential role in the Babcock-Leighton model of the solar dynamo. Describing this evolution is the aim of the surface flux transport model. The model starts from the emergence of magnetic bipoles. Thereafter, the model is based on the induction equation and the fact that after emergence the magnetic field is observed to evolve as if it were purely radial. The induction equation then describes how the surface flows -- differential rotation, meridional circulation, granular, supergranular flows, and active region inflows -- determine the evolution of the field (now taken to be purely radial). In this paper, we review the modeling of the various processes that determine the evolution of the surface field. We restrict our attention to their role in the surface flux transport model. We also discuss the success of the model and some of the results that have been obtained using this model.

Space-based observations of solar-like oscillations present an opportunity to constrain stellar models using individual mode frequencies. However, current stellar models are inaccurate near the surface, which introduces a systematic difference that must be corrected. We introduce and evaluate two parametrizations of the surface corrections based on formulae given by Gough (1990). The first we call a cubic term proportional to $\nu^3/\mathcal{I}$ and the second has an additional inverse term proportional to $\nu^{-1}/\mathcal{I}$, where $\nu$ and $\mathcal{I}$ are the frequency and inertia of an oscillation mode. We first show that these formulae accurately correct model frequencies of two different solar models ... We then incorporate the parametrizations into a modelling pipeline that simultaneously fits the surface effects and the underlying stellar model parameters. We apply this pipeline to synthetic observations of a Sun-like stellar model, solar observations degraded to typical asteroseismic uncertainties, and observations of the well-studied CoRoT target HD52265. For comparison, we also run the pipeline with the scaled power-law correction proposed by Kjeldsen et al. (2008). The fits to synthetic and degraded solar data show that the method is unbiased and produces best-fit parameters that are consistent with the input models and known parameters of the Sun. Our results for HD52265 are consistent with previous modelling efforts and the magnitude of the surface correction is similar to that of the Sun. The fit using a scaled power-law correction is significantly worse but yields consistent parameters ... We find that the cubic term alone is suitable for asteroseismic applications and it is easy to implement in an existing pipeline. ... This parametrization is thus a useful new way to correct model frequencies so that observations of individual mode frequencies can be exploited.

Context: Studying the motions on the solar surface is fundamental for understanding how turbulent convection transports energy and how magnetic fields are distributed across the solar surface. Aims: From horizontal velocity measurements all over the visible disc of the Sun and using data from the Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI), we investigate the structure and evolution of solar supergranulation. Methods: Horizontal velocity fields were measured by following the proper motions of solar granules using a newly developed version of the coherent structure tracking (CST) code. With this tool, maps of horizontal divergence were computed. We then segmented and identified supergranular cells and followed their histories by using spatio-temporal labelling. With this dataset we derived the fundamental properties of supergranulation, including their motion. Results: We find values of the fundamental parameters of supergranulation similar to previous studies: a mean lifetime of 1.5 days and a mean diameter of 25~Mm. The tracking of individual supergranular cells reveals the solar differential rotation and a poleward circulation trend of the meridional flow. The shape of the derived differential rotation and meridional flow does not depend on the cell size. If there is a background magnetic field, the diverging flows in supergranules are weaker. Conclusions: This study confirms that supergranules are suitable tracers that may be used to investigate the large-scale flows of the solar convection as long as they are detectable enough on the surface.

In time-distance helioseismology, information about the solar interior is encoded in measurements of travel times between pairs of points on the solar surface. Travel times are deduced from the cross-covariance of the random wave field. Here we consider travel times and also products of travel times as observables. They contain information about e.g. the statistical properties of convection in the Sun. The basic assumption of the model is that noise is the result of the stochastic excitation of solar waves, a random process which is stationary and Gaussian. We generalize the existing noise model (Gizon and Birch 2004) by dropping the assumption of horizontal spatial homogeneity. Using a recurrence relation, we calculate the noise covariance matrices for the moments of order 4, 6, and 8 of the observed wave field, for the moments of order 2, 3 and 4 of the cross-covariance, and for the moments of order 2, 3 and 4 of the travel times. All noise covariance matrices depend only on the expectation value of the cross-covariance of the observed wave field. For products of travel times, the noise covariance matrix consists of three terms proportional to $1/T$, $1/T^2$, and $1/T^3$, where $T$ is the duration of the observations. For typical observation times of a few hours, the term proportional to $1/T^2$ dominates and $Cov[\tau_1 \tau_2, \tau_3 \tau_4] \approx Cov[\tau_1, \tau_3] Cov[\tau_2, \tau_4] + Cov[\tau_1, \tau_4] Cov[\tau_2, \tau_3]$, where the $\tau_i$ are arbitrary travel times. This result is confirmed for $p_1$ travel times by Monte Carlo simulations and comparisons with SDO/HMI observations. General and accurate formulae have been derived to model the noise covariance matrix of helioseismic travel times and products of travel times. These results could easily be generalized to other methods of local helioseismology, such as helioseismic holography and ring diagram analysis.

The Solar Orbiter mission, to be launched in July 2017, will carry a suite of remote sensing and in-situ instruments, including the Polarimetric and Helioseismic Imager (PHI). PHI will deliver high-cadence images of the Sun in intensity and Doppler velocity suitable for carrying out novel helioseismic studies. The orbit of the Solar Orbiter spacecraft will reach a solar latitude of up to 21 deg (up to 34 deg by the end of the extended mission) and thus will enable the first local helioseismology studies of the polar regions. Here we consider an array of science objectives to be addressed by helioseismology within the baseline telemetry allocation (51 Gbit per orbit, current baseline) and within the science observing windows (baseline 3 x 10 days per orbit). A particularly important objective is the measurement of large-scale flows at high latitudes (rotation and meridional flow), which are largely unknown but play an important role in flux transport dynamos. The full range of Earth-Sun-spacecraft angles provided by the orbit will enable helioseismology from two vantage points by combining PHI with another instrument: stereoscopic helioseismology will allow the study of the deep solar interior and a better understanding of the physics of solar oscillations in both quiet Sun and sunspots. We have used a model of the PHI instrument to study its performance for helioseismology applications. As input we used a 6 hr time-series of realistic solar magneto-convection simulation (Stagger code) and the SPINOR radiative transfer code to synthesize the observables. The simulated power spectra of solar oscillations show that the instrument is suitable for helioseismology. In particular, the specified point spread function, image jitter, and photon noise are no obstacle to a successful mission.

The Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI) filtergrams, taken at six wavelengths around the Fe I 6173.3 {\AA} line, contain information about the line-of-sight velocity over a range of heights in the solar atmosphere. Multi-height velocity inferences from these observations can be exploited to study wave motions and energy transport in the atmosphere. Using realistic convection simulation datasets provided by the STAGGER and MURaM codes, we generate synthetic filtergrams and explore several methods for estimating Dopplergrams. We investigate at which height each synthetic Dopplergram correlates most strongly with the vertical velocity in the model atmospheres. On the basis of the investigation, we propose two Dopplergrams other than the standard HMI-algorithm Dopplergram produced from HMI filtergrams: a line-center Dopplergram and an average-wing Dopplergram. These two Dopplergrams correlate most strongly with vertical velocities at the heights of 30 - 40 km above (line-center) and 30 - 40 km below (average-wing) the effective height of the HMI-algorithm Dopplergram. Therefore, we can obtain velocity information from two layers separated by about a half of a scale height in the atmosphere, at best. The phase shifts between these multi-height Dopplergrams from observational data as well as those from the simulated data are also consistent with the height-difference estimates in the frequency range above the photospheric acoustic cutoff frequency.

PLATO 2.0 has recently been selected for ESA's M3 launch opportunity (2022/24). Providing accurate key planet parameters (radius, mass, density and age) in statistical numbers, it addresses fundamental questions such as: How do planetary systems form and evolve? Are there other systems with planets like ours, including potentially habitable planets? The PLATO 2.0 instrument consists of 34 small aperture telescopes (32 with 25 sec readout cadence and 2 with 2.5 sec candence) providing a wide field-of-view (2232 deg2) and a large photometric magnitude range (4-16 mag). It focusses on bright (4-11 mag) stars in wide fields to detect and characterize planets down to Earth-size by photometric transits, whose masses can then be determined by ground-based radial-velocity follow-up measurements. Asteroseismology will be performed for these bright stars to obtain highly accurate stellar parameters, including masses and ages. The combination of bright targets and asteroseismology results in high accuracy for the bulk planet parameters: 2%, 4-10% and 10% for planet radii, masses and ages, respectively. The planned baseline observing strategy includes two long pointings (2-3 years) to detect and bulk characterize planets reaching into the habitable zone (HZ) of solar-like stars and an additional step-and-stare phase to cover in total about 50% of the sky. PLATO 2.0 will observe up to 1,000,000 stars and detect and characterize hundreds of small planets, and thousands of planets in the Neptune to gas giant regime out to the HZ. It will therefore provide the first large-scale catalogue of bulk characterized planets with accurate radii, masses, mean densities and ages. This catalogue will include terrestrial planets at intermediate orbital distances, where surface temperatures are moderate. Coverage of this parameter range with statistical numbers of bulk characterized planets is unique to PLATO 2.0.

Context. The object HD 43587Aa is a G0V star observed during the 145-day LRa03 run of the COnvection, ROtation and planetary Transits space mission (CoRoT), for which complementary High Accuracy Radial velocity Planet Searcher (HARPS) spectra with S/N>300 were also obtained. Its visual magnitude is 5.71, and its effective temperature is close to 5950 K. It has a known companion in a highly eccentric orbit and is also coupled with two more distant companions. Aims. We undertake a preliminary investigation of the internal structure of HD 43587Aa. Methods. We carried out a seismic analysis of the star, using maximum likelihood estimators and Markov Chain Monte Carlo methods. Results. We established the first table of the eigenmode frequencies, widths, and heights for HD 43587Aa. The star appears to have a mass and a radius slightly larger than the Sun, and is slightly older (5.6 Gyr). Two scenarios are suggested for the geometry of the star: either its inclination angle is very low, or the rotation velocity of the star is very low. Conclusions. A more detailed study of the rotation and of the magnetic and chromospheric activity for this star is needed, and will be the subject of a further study. New high resolution spectrometric observations should be performed for at least several months in duration.

Context : We still do not know which mechanisms are responsible for the transport of angular momentum inside stars. The recent detection of mixed modes that contain the signature of rotation in the spectra of Kepler subgiants and red giants gives us the opportunity to make progress on this issue. Aims: Our aim is to probe the radial dependance of the rotation profiles for a sample of Kepler targets. For this purpose, subgiants and early red giants are particularly interesting targets because their rotational splittings are more sensitive to the rotation outside the deeper core than is the case for their more evolved counterparts. Methods: We first extract the rotational splittings and frequencies of the modes for six young Kepler red giants. We then perform a seismic modeling of these stars using the evolutionary codes CESAM2k and ASTEC. By using the observed splittings and the rotational kernels of the optimal models, we perform inversions of the internal rotation profiles of the six stars. Results: We obtain estimates of the mean rotation rate in the core and in the convective envelope of these stars. We show that the rotation contrast between the core and the envelope increases during the subgiant branch. Our results also suggest that the core of subgiants spins up with time, contrary to the RGB stars whose core has been shown to spin down. For two of the stars, we show that a discontinuous rotation profile with a deep discontinuity reproduces the observed splittings significantly better than a smooth rotation profile. Interestingly, the depths that are found most probable for the discontinuities roughly coincide with the location of the H-burning shell, which separates the layers that contract from those that expand. These results will bring observational constraints to the scenarios of angular momentum transport in stars.

We study Doppler velocity measurements at multiple heights in the solar atmosphere using a set of six filtergrams obtained by the Helioseismic magnetic Imager on board the Solar Dynamics Observatory. There are clear and significant phase differences between core and wing Dopplergrams in the frequency range above the photospheric acoustic cutoff frequency, which indicates that these are really "multi-height" datasets.

Linear time-domain simulations of acoustic oscillations are unstable in the stellar convection zone. To overcome this problem it is customary to compute the oscillations of a stabilized background stellar model. The stabilization, however, affects the result. Here we propose to use a perturbative approach (running the simulation twice) to approximately recover the acoustic wave field, while preserving seismic reciprocity. To test the method we considered a 1D standard solar model. We found that the mode frequencies of the (unstable) standard solar model are well approximated by the perturbative approach within $1~\mu$Hz for low-degree modes with frequencies near $3~\mu$Hz. We also show that the perturbative approach is appropriate for correcting rotational-frequency kernels. Finally, we comment that the method can be generalized to wave propagation in 3D magnetized stellar interiors because the magnetic fields have stabilizing effects on convection.

We use asteroseismic data obtained by the NASA Kepler Mission to estimate the fundamental properties of more than 500 main-sequence and sub-giant stars. Data obtained during the first 10 months of Kepler science operations were used for this work, when these solar-type targets were observed for one month each in a survey mode. Stellar properties have been estimated using two global asteroseismic parameters and complementary photometric and spectroscopic data. Homogeneous sets of effective temperatures were available for the entire ensemble from complementary photometry; spectroscopic estimates of T_eff and [Fe/H] were available from a homogeneous analysis of ground-based data on a subset of 87 stars. [Abbreviated version... see paper for full abstract.]

We present the first high-resolution solar images in the Mg II k 2796 \AA{} line. The images, taken through a 4.8 \AA{} broad interference filter, were obtained during the second science flight of SUNRISE in June 2013 by the SuFI instrument. The Mg II k images display structures that look qualitatively very similar to images taken in the core of Ca II H. The Mg II images exhibit reversed granulation (or shock waves) in the internetwork regions of the quiet Sun, at intensity contrasts that are similar to those found in Ca II H. Very prominent in Mg II are bright points, both in the quiet Sun and in plage regions, particularly near disk center. These are much brighter than at other wavelengths sampled at similar resolution. Furthermore, Mg II k images also show fibril structures associated with plage regions. Again, the fibrils are similar to those seen in Ca II H images, but tend to be more pronounced, particularly in weak plage.

High-quality data from space-based observatories present an opportunity to fit stellar models to observations of individually-identified oscillation frequencies, not just the large and small frequency separations. But such fits require the evaluation of a large number of accurate stellar models, which remains expensive. Here, we show that global-mode oscillation frequencies interpolated in a grid of stellar models are precise and accurate, at least in the neighbourhood of a solar model.

Aims: We aim to measure the starspot rotation periods of active stars in the Kepler field as a function of spectral type and to extend reliable rotation measurements from F-, G-, and K-type to M-type stars. Methods: Using the Lomb-Scargle periodogram we searched more than 150 000 stellar light curves for periodic brightness variations. We analyzed periods between 1 and 30 days in eight consecutive Kepler quarters, where 30 days is an estimated maximum for the validity of the PDC_MAP data correction pipeline. We selected stable rotation periods, i.e., periods that do not vary from the median by more than one day in at least six of the eight quarters. We averaged the periods for each stellar spectral class according to B - V color and compared the results to archival vsini data, using stellar radii estimates from the Kepler Input Catalog. Results: We report on the stable starspot rotation periods of 12 151 Kepler stars. We find good agreement between starspot velocities and vsini data for all F-, G- and early K-type stars. The 795 M-type stars in our sample have a median rotation period of 15.4 days. We find an excess of M-type stars with periods less than 7.5 days that are potentially fast-rotating and fully convective. Measuring photometric variability in multiple Kepler quarters appears to be a straightforward and reliable way to determine the rotation periods of a large sample of active stars, including late-type stars.

We comment on the potential for continuing asteroseismology of solar-type and red-giant stars in a 2-wheel Kepler Mission. Our main conclusion is that by targeting stars in the ecliptic it should be possible to perform high-quality asteroseismology, as long as favorable scenarios for 2-wheel pointing performance are met. Targeting the ecliptic would potentially facilitate unique science that was not possible in the nominal Mission, notably from the study of clusters that are significantly brighter than those in the Kepler field. Our conclusions are based on predictions of 2-wheel observations made by a space photometry simulator, with information provided by the Kepler Project used as input to describe the degraded pointing scenarios. We find that elevated levels of frequency-dependent noise, consistent with the above scenarios, would have a significant negative impact on our ability to continue asteroseismic studies of solar-like oscillators in the Kepler field. However, the situation may be much more optimistic for observations in the ecliptic, provided that pointing resets of the spacecraft during regular desaturations of the two functioning reaction wheels are accurate at the < 1 arcsec level. This would make it possible to apply a post-hoc analysis that would recover most of the lost photometric precision. Without this post-hoc correction---and the accurate re-pointing it requires---the performance would probably be as poor as in the Kepler-field case. Critical to our conclusions for both fields is the assumed level of pointing noise (in the short-term jitter and the longer-term drift). We suggest that further tests will be needed to clarify our results once more detail and data on the expected pointing performance becomes available, and we offer our assistance in this work.

Rotation is thought to drive cyclic magnetic activity in the Sun and Sun-like stars. Stellar dynamos, however, are poorly understood owing to the scarcity of observations of rotation and magnetic fields in stars. Here, inferences are drawn on the internal rotation of a distant Sun-like star by studying its global modes of oscillation. We report asteroseismic constraints imposed on the rotation rate and the inclination of the spin axis of the Sun-like star HD 52265, a principal target observed by the CoRoT satellite that is known to host a planetary companion. These seismic inferences are remarkably consistent with an independent spectroscopic observation (rotational line broadening) and with the observed rotation period of star spots. Furthermore, asteroseismology constrains the mass of exoplanet HD 52265b. Under the standard assumption that the stellar spin axis and the axis of the planetary orbit coincide, the minimum spectroscopic mass of the planet can be converted into a true mass of 1.85 (+0.52,-0.42) M_Jupiter, which implies that it is a planet, not a brown dwarf.

Measurements of seismic wave travel times at the photosphere of the Sun have enabled inferences of its interior structure and dynamics. In interpreting these measurements, the simplifying assumption that waves propagate through a temporally stationary medium is almost universally invoked. However, the Sun is in a constant state of evolution, on a broad range of spatio-temporal scales. At the zero wavelength limit, i.e., when the wavelength is much shorter than the scale over which the medium varies, the WKBJ (ray) approximation may be applied. Here, we address the other asymptotic end of the spectrum, the infinite wavelength limit, using the technique of homogenization. We apply homogenization to scenarios where waves are propagating through rapidly varying media (spatially and temporally), and derive effective models for the media. One consequence is that a scalar sound speed becomes a tensorial wavespeed in the effective model and anisotropies can be induced depending on the nature of the perturbation. The second term in this asymptotic two-scale expansion, the so-called corrector, contains contributions due to higher-order scattering, leading to the decoherence of the wavefield. This decoherence may be causally linked to the observed wave attenuation in the Sun. Although the examples we consider here consist of periodic arrays of perturbations to the background, homogenization may be extended to ergodic and stationary random media. This method may have broad implications for the manner in which we interpret seismic measurements in the Sun and for modeling the effects of granulation on the scattering of waves and distortion of normal-mode eigenfunctions.

We compare measurements of horizontal flows on the surface of the Sun using helioseismic time--distance inversions and coherent structure tracking of solar granules. Tracking provides 2D horizontal flows on the solar surface, whereas the time--distance inversions estimate the full 3-D velocity flows in the shallow near-surface layers. Both techniques use HMI observations as an input. We find good correlations between the various measurements resulting from the two techniques. Further, we find a good agreement between these measurements and the time-averaged Doppler line-of-sight velocity, and also perform sanity checks on the vertical flow that resulted from the 3-D time--distance inversion.

In order to assess the ability of helioseismology to probe the subsurface structure and magnetic field of sunspots, we need to determine how helioseismic travel times depend on perturbations to sunspot models. Here we numerically simulate the propagation of f, p1, and p2 wave packets through magnetic sunspot models. Among the models we considered, a ~50 km change in the height of the Wilson depression and a change in the subsurface magnetic field geometry can both be detected above the observational noise level. We also find that the travel-time shifts due to changes in a sunspot model must be modeled by computing the effects of changing the reference sunspot model, and not by computing the effects of changing the subsurface structure in the quiet-Sun model. For p1 modes the latter is wrong by a factor of four. In conclusion, numerical modeling of MHD wave propagation is an essential tool for the interpretation of the effects of sunspots on seismic waveforms.

The measurement of the Sun's surface motions with a high spatial and temporal resolution is still a challenge. We wish to validate horizontal velocity measurements all over the visible disk of the Sun from Solar Dynamics Observatory/ Helioseismic and Magnetic Imager (SDO/HMI) data. Horizontal velocity fields are measured by following the proper motions of solar granules using a newly developed version of the Coherent Structure Tracking (CST) code. The comparison of the surface flows measured at high spatial resolution (Hinode, 0.1 arcsec) and low resolution (SDO/HMI, 0.5 arcsec) allows us to determine corrections to be applied to the horizontal velocity measured from HMI white light data. We derive horizontal velocity maps with spatial and temporal resolutions of respectively 2.5 Mm and 30 min. From the two components of the horizontal velocity Vx and Vy measured in the sky plane and the simultaneous line of sight component from SDO/HMI dopplergrams v_D, we derive the spherical velocity components (Vr, Vtheta, Vphi). The azimuthal component Vphi gives the solar differential rotation with a high precision (+-0.037km/s) from a temporal sequence of only three hours. By following the proper motions of the solar granules, we can revisit the dynamics of the solar surface at high spatial and temporal resolutions from hours to months and years with the SDO data.

The Square Kilometre Array (SKA) is the most ambitious radio telescope ever planned. With a collecting area of about a square kilometre, the SKA will be far superior in sensitivity and observing speed to all current radio facilities. The scientific capability promised by the SKA and its technological challenges provide an ideal base for interdisciplinary research, technology transfer, and collaboration between universities, research centres and industry. The SKA in the radio regime and the European Extreme Large Telescope (E-ELT) in the optical band are on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI) and have been recognised as the essential facilities for European research in astronomy. This "White Paper" outlines the German science and R&D interests in the SKA project and will provide the basis for future funding applications to secure German involvement in the Square Kilometre Array.

Convection is the mechanism by which energy is transported through the outermost 30% of the Sun. Solar turbulent convection is notoriously difficult to model across the entire convection zone where the density spans many orders of magnitude. In this issue of PNAS, Hanasoge et al. (2012) employ recent helioseismic observations to derive stringent empirical constraints on the amplitude of large-scale convective velocities in the solar interior. They report an upper limit that is far smaller than predicted by a popular hydrodynamic numerical simulation.

Sunspots are prominent manifestations of solar magnetoconvection and imaging their subsurface structure is an outstanding problem of wide physical importance. Travel times of seismic waves that propagate through these structures are typically used as inputs to inversions. Despite the presence of strongly anisotropic magnetic waveguides, these measurements have always been interpreted in terms of changes to isotropic wavespeeds and flow-advection related Doppler shifts. Here, we employ PDE-constrained optimization to determine the appropriate parameterization of the structural properties of the magnetic interior. Seven different wavespeeds fully characterize helioseismic wave propagation: the isotropic sound speed, a Doppler-shifting flow-advection velocity and an anisotropic magnetic velocity. The structure of magnetic media is sensed by magnetoacoustic slow and fast modes and Alfv\'{e}n waves, each of which propagates at a different wavespeed. We show that even in the case of weak magnetic fields, significant errors may be incurred if these anisotropies are not accounted for in inversions. Translation invariance is demonstrably lost. These developments render plausible the accurate seismic imaging of magnetoconvection in the Sun.

Rotation is expected to have an important influence on the structure and the evolution of stars. However, the mechanisms of angular momentum transport in stars remain theoretically uncertain and very complex to take into account in stellar models. To achieve a better understanding of these processes, we desperately need observational constraints on the internal rotation of stars, which until very recently were restricted to the Sun. In this paper, we report the detection of mixed modes - i.e. modes that behave both as g modes in the core and as p modes in the envelope - in the spectrum of the early red giant KIC7341231, which was observed during one year with the Kepler spacecraft. By performing an analysis of the oscillation spectrum of the star, we show that its non-radial modes are clearly split by stellar rotation and we are able to determine precisely the rotational splittings of 18 modes. We then find a stellar model that reproduces very well the observed atmospheric and seismic properties of the star. We use this model to perform inversions of the internal rotation profile of the star, which enables us to show that the core of the star is rotating at least five times faster than the envelope. This will shed new light on the processes of transport of angular momentum in stars. In particular, this result can be used to place constraints on the angular momentum coupling between the core and the envelope of early red giants, which could help us discriminate between the theories that have been proposed over the last decades.

This paper presents a detailed and precise study of the characteristics of the Exoplanet Host Star and CoRoT main target HD 52265, as derived from asteroseismic studies. The results are compared with previous estimates, with a comprehensive summary and discussion. The basic method is similar to that previously used by the Toulouse group for solar-type stars. Models are computed with various initial chemical compositions and the computed p-mode frequencies are compared with the observed ones. All models include atomic diffusion and the importance of radiative accelerations is discussed. Several tests are used, including the usual frequency combinations and the fits of the \'echelle diagrams. The possible surface effects are introduced and discussed. Automatic codes are also used to find the best model for this star (SEEK, AMP) and their results are compared with that obtained with the detailed method. We find precise results for the mass, radius and age of this star, as well as its effective temperature and luminosity. We also give an estimate of the initial helium abundance. These results are important for the characterization of the star-planet system.

For the first time, the motion of granules (solar plasma on the surface on scales larger than 2.5 Mm) has been followed over the entire visible surface of the Sun, using SDO/HMI white-light data. Horizontal velocity fields are derived from image correlation tracking using a new version of the coherent structure tracking algorithm.The spatial and temporal resolutions of the horizontal velocity map are 2.5 Mm and 30 min respectively . From this reconstruction, using the multi-resolution analysis, one can obtain to the velocity field at different scales with its derivatives such as the horizontal divergence or the vertical component of the vorticity. The intrinsic error on the velocity is ~0.25 km/s for a time sequence of 30 minutes and a mesh size of 2.5 Mm.This is acceptable compared to the granule velocities, which range between 0.3 km/s and 1.8 km/s. A high correlation between velocities computed from Hinode and SDO/HMI has been found (85%). From the data we derive the power spectrum of the supergranulation horizontal velocity field, the solar differential rotation, and the meridional velocity.

The evolved solar-type stars 16 Cyg A & B have long been studied as solar analogs, yielding a glimpse into the future of our own Sun. The orbital period of the binary system is too long to provide meaningful dynamical constraints on the stellar properties, but asteroseismology can help because the stars are among the brightest in the Kepler field. We present an analysis of three months of nearly uninterrupted photometry of 16 Cyg A & B from the Kepler space telescope. We extract a total of 46 and 41 oscillation frequencies for the two components respectively, including a clear detection of octupole (l=3) modes in both stars. We derive the properties of each star independently using the Asteroseismic Modeling Portal, fitting the individual oscillation frequencies and other observational constraints simultaneously. We evaluate the systematic uncertainties from an ensemble of results generated by a variety of stellar evolution codes and fitting methods. The optimal models derived by fitting each component individually yield a common age (t=6.8+/-0.4 Gyr) and initial composition (Z_i=0.024+/-0.002, Y_i=0.25+/-0.01) within the uncertainties, as expected for the components of a binary system, bolstering our confidence in the reliability of asteroseismic techniques. The longer data sets that will ultimately become available will allow future studies of differential rotation, convection zone depths, and long-term changes due to stellar activity cycles.

Inversions for local helioseismology are an important and necessary step for obtaining three-dimensional maps of various physical quantities in the solar interior. Frequently, the full inverse problems that one would like to solve prove intractable because of computational constraints. Due to the enormous seismic data sets that already exist and those forthcoming, this is a problem that needs to be addressed. To this end, we present a very efficient linear inversion algorithm for local helioseismology. It is based on a subtractive optimally localized averaging (SOLA) scheme in the Fourier domain, utilizing the horizontal-translation invariance of the sensitivity kernels. In Fourier space the problem decouples into many small problems, one for each horizontal wave vector. This multi-channel SOLA method is demonstrated for an example problem in time-distance helioseismology that is small enough to be solved both in real and Fourier space. We find that both approaches are successful in solving the inverse problem. However, the multi-channel SOLA algorithm is much faster and can easily be parallelized.

In this paper, we construct background solar models that are stable against convection, by modifying the vertical pressure gradient of Model S (Christensen-Dalsgaard et al., 1996, Science, 272, 1286) relinquishing hydrostatic equilibrium. However, the stabilisation affects the eigenmodes that we wish to remain as close to Model S as possible. In a bid to recover the Model S eigenmodes, we choose to make additional corrections to the sound speed of Model S before stabilisation. No stabilised model can be perfectly solar-like, so we present three stabilised models with slightly different eigenmodes. The models are appropriate to study the f and p1 to p4 modes with spherical harmonic degrees in the range from 400 to 900. Background model CSM has a modified pressure gradient for stabilisation and has eigenfrequencies within 2% of Model S. Model CSM_A has an additional 10% increase in sound speed in the top 1 Mm resulting in eigenfrequencies within 2% of Model S and eigenfunctions that are, in comparison with CSM, closest to those of Model S. Model CSM_B has a 3% decrease in sound speed in the top 5 Mm resulting in eigenfrequencies within 1% of Model S and eigenfunctions that are only marginally adversely affected. These models are useful to study the interaction of solar waves with embedded three-dimensional heterogeneities, such as convective flows and model sunspots. We have also calculated the response of the stabilised models to excitation by random near-surface sources, using simulations of the propagation of linear waves. We find that the simulated power spectra of wave motion are in good agreement with an observed SOHO/MDI power spectrum. Overall, our convectively stabilised background models provide a good basis for quantitative numerical local helioseismology. The models are available for download from http://www.mps.mpg.de/projects/seismo/NA4/.

For a given {\it misfit function}, a specified optimality measure of a model, its gradient describes the manner in which one may alter properties of the system to march towards a stationary point. The adjoint method, arising from partial-differential-equation-constrained optimization, describes a means of extracting derivatives of a misfit function with respect to model parameters through finite computation. It relies on the accurate calculation of wavefields that are driven by two types of sources, namely the average wave-excitation spectrum, resulting in the {\it forward wavefield}, and differences between predictions and observations, resulting in an {\it adjoint wavefield}. All sensitivity kernels relevant to a given measurement emerge directly from the evaluation of an interaction integral involving these wavefields. The technique facilitates computation of sensitivity kernels (Fr\'{e}chet derivatives) relative to three-dimensional heterogeneous background models, thereby paving the way for non-linear iterative inversions. An algorithm to perform such inversions using as many observations as desired is discussed.

The star HD 52265 is a G0V metal-rich exoplanet-host star observed in the seismology field of the CoRoT space telescope from November 2008 to March 2009. The satellite collected 117 days of high-precision photometric data on this star, showing that it presents solar-like oscillations. HD 52265 was also observed in spectroscopy with the Narval spectrograph at the same epoch. We characterise HD 52265 using both spectroscopic and seismic data. The fundamental stellar parameters of HD 52265 were derived with the semi-automatic software VWA, and the projected rotational velocity was estimated by fitting synthetic profiles to isolated lines in the observed spectrum. The parameters of the observed p modes were determined with a maximum-likelihood estimation. We performed a global fit of the oscillation spectrum, over about ten radial orders, for degrees l=0 to 2. We also derived the properties of the granulation, and analysed a signature of the rotation induced by the photospheric magnetic activity. Precise determinations of fundamental parameters have been obtained: Teff = 6100 +- 60 K, log g = 4.35 +- 0.09, [M/H] = 0.19 +- 0.05, as well as vsini = 3.6 +0.3 -1.0 km/s. We have measured a mean rotation period P_rot = 12.3 +- 0.15 days, and find a signature of differential rotation. The frequencies of 31 modes are reported in the range 1500-2550 micro-Hz. The large separation exhibits a clear modulation around the mean value <Dnu> = 98.3 +- 0.1 micro-Hz. Mode widths vary with frequency along an S-shape with a clear local maximum around 1800 micro-Hz. We deduce lifetimes ranging between 0.5 and 3 days for these modes. Finally, we find a maximal bolometric amplitude of about 3.96 +- 0.24 ppm for radial modes.

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.

Three-dimensional numerical simulations have been used to study the scattering of a surface-gravity wave packet by vertical magnetic flux tubes, with radii from 200 km to 3 Mm, embedded in stratified polytropic atmosphere. The scattered wave was found to consist primarily of m=0 (axisymmetric) and m=1 modes. It was found that the ratio of the amplitude of these two modes is strongly dependant on the radius of the flux tube: The kink mode is the dominant mode excited in tubes with a small radius while the sausage mode is dominant for large tubes. Simulations of this type provide a simple, efficient and robust way to start understanding the seismic signature of flux tubes, which have recently began to be observed.

We observe and characterize the scattering of acoustic wave packets by a sunspot, in a regime where the wavelength is comparable to the size of the sunspot. Spatial maps of wave traveltimes and amplitudes are measured from the cross-covariance function of the random wave field. The averaging procedure is such that incoming wave packets are plane wave packets. Observations show that the magnitude of the traveltime perturbation caused by the sunspot diminishes as waves propagate away from the sunspot -- a finite-wavelength phenomenon known as wavefront healing. Observations also show a reduction of the amplitude of the waves after their passage through the sunspot. A significant fraction of this amplitude reduction is due to the defocusing of wave energy by the fast wave-speed perturbation introduced by the sunspot. This ``geometrical attenuation'' will contribute to the wave amplitude reduction in addition to the physical absorption of waves. In addition, we observe an enhancement of wave amplitude away from the central path: diffracted rays intersect with unperturbed rays (caustics) and wavefronts fold and triplicate. Thus we find that ray tracing is useful to interpret these phenomena, although it cannot explain wavefront healing.

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.

Perfectly matched layers are a very efficient and accurate way to absorb waves in media. We present a stable convolutional unsplit perfectly matched formulation designed for the linearized stratified Euler equations. However, the technique as applied to the Magneto-hydrodynamic (MHD) equations requires the use of a sponge, which, despite placing the perfectly matched status in question, is still highly efficient at absorbing outgoing waves. We study solutions of the equations in the backdrop of models of linearized wave propagation in the Sun. We test the numerical stability of the schemes by integrating the equations over a large number of wave periods.

A data and computation center for helioseismology has been set up at the Max Planck Institute for Solar System Research in Germany to prepare for the SDO mission. Here we present the system infrastructure and the scientific aims of this project, which is funded through grants from the German Aerospace Center and the European Research Council.

One goal of helioseismology is to determine the subsurface structure of sunspots. In order to do so, it is important to understand first the near-surface effects of sunspots on solar waves, which are dominant. Here we construct simplified, cylindrically-symmetric sunspot models, which are designed to capture the magnetic and thermodynamics effects coming from about 500 km below the quiet-Sun $\tau_{5000}=1$ level to the lower chromosphere. We use a combination of existing semi-empirical models of sunspot thermodynamic structure (density, temperature, pressure): the umbral model of Maltby et al. (1986) and the penumbral model of Ding and Fang (1989). The OPAL equation of state tables are used to derive the sound speed profile. We smoothly merge the near-surface properties to the quiet-Sun values about 1mm below the surface. The umbral and penumbral radii are free parameters. The magnetic field is added to the thermodynamic structure, without requiring magnetostatic equilibrium. The vertical component of the magnetic field is assumed to have a Gaussian horizontal profile, with a maximum surface field strength fixed by surface observations. The full magnetic field vector is solenoidal and determined by the on-axis vertical field, which, at the surface, is chosen such that the field inclination is 45$^\circ$ at the umbral-penumbral boundary. We construct a particular sunspot model based on SOHO/MDI observations of the sunspot in active region NOAA 9787. The helioseismic signature of the model sunspot is studied using numerical simulations of the propagation of f, p$_1$, and p$_2$ wave packets. These simulations are compared against cross-covariances of the observed wave field. We find that the sunspot model gives a helioseismic signature that is similar to the observations.

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.

Time-distance helioseismology and related techniques show great promise for probing the structure and dynamics of the subphotospheric layers of the Sun. Indeed time-distance helioseismology has already been applied to make inferences about structures and flows under sunspots and active regions, to map long-lived convective flow patterns, and so on. Yet certainly there are still many inadequacies in the current approaches and, as the data get better and the questions we seek to address get more subtle, methods that were previously regarded as adequate are no longer acceptable. Here we give a short and partial description of outstanding problems in local helioseismology, using time-distance helioseismology as a guiding example.

In this paper we describe the semi-spectral linear MHD (SLiM) code which we have written to follow the interaction of linear waves through an inhomogeneous three-dimensional solar atmosphere. The background model allows almost arbitrary perturbations of density, temperature, sound speed as well as magnetic and velocity fields. We give details of several of the tests we have used to check the code. The code will be useful in understanding the helioseismic signatures of various solar features, including sunspots.

Time-distance helioseismology is a technique for measuring the time for waves to travel from one point on the solar surface to another. These wave travel times are affected by advection by subsurface flows. Inferences of plasma flows based on observed travel times depend critically on the ability to accurately model the effects of subsurface flows on time-distance measurements. We present a Born approximation based computation of the sensitivity of time distance travel times to weak, steady, inhomogeneous subsurface flows. Three sensitivity functions are obtained, one for each component of the 3D vector flow. We show that the depth sensitivity of travel times to horizontally uniform flows is given approximately by the kinetic energy density of the oscillation modes which contribute to the travel times. For flows with strong depth dependence, the Born approximation can give substantially different results than the ray approximation.

Solar Orbiter is intended to become ESA's next solar mission in heritage of the successful SOHO project. The scientific objectives of the mission, its design, and its scientific payload are reviewed. Specific emphasis is given to the perspectives of Solar Orbiter with respect to helioseismology.

In local helioseismology numerical simulations of wave propagation are useful to model the interaction of solar waves with perturbations to a background solar model. However, the solution to the equations of motions include convective modes that can swamp the waves we are interested in. For this reason, we choose to first stabilise the background solar model against convection by altering the vertical pressure gradient. Here we compare the eigenmodes of our convectively stabilised model with a standard solar model (Model S) and find a good agreement.

The Sun supports a rich spectrum of internal waves that are continuously excited by turbulent convection. The GONG network and the MDI/SOHO space instrument provide an exceptional data base of spatially-resolved observations of solar oscillations, covering an entire sunspot cycle (11 years). Local helioseismology is a set of tools for probing the solar interior in three dimensions using measurements of wave travel times and local mode frequencies. Local helioseismology has discovered (i) near-surface vector flows associated with convection (ii) 250 m/s subsurface horizontal outflows around sunspots (iii) ~50 m/s extended horizontal flows around active regions (converging near the surface and diverging below), (iv) the effect of the Coriolis force on convective flows and active region flows (v) the subsurface signature of the 15 m/s poleward meridional flow, (vi) a +/-5 m/s time-varying depth-dependent component of the meridional circulation around the mean latitude of activity, and (vii) magnetic activity on the far side of the Sun.

ASTROD I is a planned interplanetary space mission with multiple goals. The primary aims are: to test General Relativity with an improvement in sensitivity of over 3 orders of magnitude, improving our understanding of gravity and aiding the development of a new quantum gravity theory; to measure key solar system parameters with increased accuracy, advancing solar physics and our knowledge of the solar system and to measure the time rate of change of the gravitational constant with an order of magnitude improvement and the anomalous Pioneer acceleration, thereby probing dark matter and dark energy gravitationally. It is an international project, with major contributions from Europe and China and is envisaged as the first in a series of ASTROD missions. ASTROD I will consist of one spacecraft carrying a telescope, four lasers, two event timers and a clock. Two-way, two-wavelength laser pulse ranging will be used between the spacecraft in a solar orbit and deep space laser stations on Earth, to achieve the ASTROD I goals. A second mission, ASTROD II is envisaged as a three-spacecraft mission which would test General Relativity to one part per billion, enable detection of solar g-modes, measure the solar Lense-Thirring effect to 10 parts per million, and probe gravitational waves at frequencies below the LISA bandwidth. In the third phase (ASTROD III or Super-ASTROD), larger orbits could be implemented to map the outer solar system and to probe primordial gravitational-waves at frequencies below the ASTROD II bandwidth.

The main goals of the HELAS local helioseismology network activity are to consolidate this field of research in Europe, to organise scientific workshops, and to facilitate the distribution of observations and data analysis software. Most of this is currently accomplished via a dedicated website at http://www.mps.mpg.de/projects/seismo/NA4/ . In this paper we list the outreach material, observational data, analysis tools and modelling tools currently available from the website and describe the focus of the scientific workshops and their proceedings.

The POLAR Investigation of the Sun (POLARIS) mission uses a combination of a gravity assist and solar sail propulsion to place a spacecraft in a 0.48 AU circular orbit around the Sun with an inclination of 75 degrees with respect to solar equator. This challenging orbit is made possible by the challenging development of solar sail propulsion. This first extended view of the high-latitude regions of the Sun will enable crucial observations not possible from the ecliptic viewpoint or from Solar Orbiter. While Solar Orbiter would give the first glimpse of the high latitude magnetic field and flows to probe the solar dynamo, it does not have sufficient viewing of the polar regions to achieve POLARIS' primary objective : determining the relation between the magnetism and dynamics of the Sun's polar regions and the solar cycle.

With the aim of studying magnetic effects in time-distance helioseismology, we use the first-order Born approximation to compute the scattering of acoustic plane waves by a magnetic cylinder embedded in a uniform medium. We show, by comparison with the exact solution, that the travel-time shifts computed in the Born approximation are everywhere valid to first order in the ratio of the magnetic to the gas pressures. We also show that, for arbitrary magnetic field strength, the Born approximation is not valid in the limit where the radius of the magnetic cylinder tends to zero.

Solar oscillations consist of a rich spectrum of internal acoustic waves and surface gravity waves, stochastically excited by turbulent convection. They have been monitored almost continuously over the last ten years with high-precision Doppler images of the solar surface. The purpose of helioseismology is to retrieve information about the structure and the dynamics of the solar interior from the frequencies, phases, and amplitudes of solar waves. Methods of analysis are being developed to make three-dimensional images of subsurface motions and temperature inhomogeneities in order to study convective structures and regions of magnetic activity, like sunspots.

Quantitative helio- and asteroseismology require very precise measurements of the frequencies, amplitudes, and lifetimes of the global modes of stellar oscillation. It is common knowledge that the precision of these measurements depends on the total length (T), quality, and completeness of the observations. Except in a few simple cases, the effect of gaps in the data on measurement precision is poorly understood, in particular in Fourier space where the convolution of the observable with the observation window introduces correlations between different frequencies. Here we describe and implement a rather general method to retrieve maximum likelihood estimates of the oscillation parameters, taking into account the proper statistics of the observations. Our fitting method applies in complex Fourier space and exploits the phase information. We consider both solar-like stochastic oscillations and long-lived harmonic oscillations, plus random noise. Using numerical simulations, we demonstrate the existence of cases for which our improved fitting method is less biased and has a greater precision than when the frequency correlations are ignored. This is especially true of low signal-to-noise solar-like oscillations. For example, we discuss a case where the precision on the mode frequency estimate is increased by a factor of five, for a duty cycle of 15%. In the case of long-lived sinusoidal oscillations, a proper treatment of the frequency correlations does not provide any significant improvement; nevertheless we confirm that the mode frequency can be measured from gapped data at a much better precision than the 1/T Rayleigh resolution.

We present independent observations of the solar-cycle variation of flows near the solar surface and at a depth of about 60 Mm, in the latitude range $\pm45^\circ$. We show that the time-varying components of the meridional flow at these two depths have opposite sign, while the time-varying components of the zonal flow are in phase. This is in agreement with previous results. We then investigate whether the observations are consistent with a theoretical model of solar-cycle dependent meridional circulation based on a flux-transport dynamo combined with a geostrophic flow caused by increased radiative loss in the active region belt (the only existing quantitative model). We find that the model and the data are in qualitative agreement, although the amplitude of the solar-cycle variation of the meridional flow at 60 Mm is underestimated by the model.

We calculate the scattering effects associated with the interaction of a surface gravity or {\it f} mode with a thin magnetic flux tube embedded in a realistically stratified medium. We find that the dominant scattered wave is an $f$ mode with amplitude and phase of 1.17% and around $49^\circ$ relative to the incident wave, compared to the values of 0.13% and $40^\circ$ estimated from observations. The extent of scattering into high-order acoustic $p$ modes is too weak to be accurately characterized. We recover the result that the degree of scattering is enhanced as (a) the frequency of the incident wave increases and (b) the flux tube becomes magnetically dominated.

To recover the flow information encoded in travel-time data of time-distance helioseismology, accurate forward modeling and a robust inversion of the travel times are required. We accomplish this using three-dimensional finite-frequency travel-time sensitivity kernels for flows along with a 2+1 dimensional (2+1D) optimally localized averaging (OLA) inversion scheme. Travel times are measured by ridge filtering MDI full-disk Doppler data and the corresponding Born sensitivity kernels are computed for these particular travel times. We also utilize the full noise covariance properties of the travel times which allow us to accurately estimate the errors for all inversions. The whole procedure is thus fully consistent. Due to ridge filtering, the kernel functions separate in the horizontal and vertical directions, motivating our choice of a 2+1D inversion implementation. The inversion procedure also minimizes cross-talk effects among the three flow components, and the averaging kernels resulting from the inversion show very small amounts of cross-talk. We obtain three-dimensional maps of vector solar flows in the quiet Sun at spatial resolutions of 7-10 Mm using generally 24 h of data. For all of the flow maps we provide averaging kernels and the noise estimates. We present examples to test the inferred flows, such as a comparison with Doppler data, in which we find a correlation of 0.9. We also present results for quiet-Sun supergranular flows at different depths in the upper convection zone.

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.

We present an up-to-date estimate for the prospect of using the Astrodynamical Space Test of Relativity using Optical Devices (ASTROD) for an unambiguous detection of solar g modes (f < 400 micro Hertz) through their gravitational signature. There are currently two major efforts to detect low-frequency gravitational effects, ASTROD and the Laser Interferometer Space Antenna (LISA). Using the most recent g mode surface amplitude estimates, both observational and theoretical, it is unclear whether LISA will be capable of successfully detecting these modes. The ASTROD project may be better suited for detection as its sensitivity curve is shifted towards lower frequencies with the best sensitivity occurring in the range 100-300 micro Hertz.

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.

We perform a two-dimensional inversion of f-mode travel times to determine near-surface solar flows. The inversion is based on optimally localized averaging of travel times. We use finite-wavelength travel-time sensitivity functions and a realistic model of the data errors. We find that it is possible to obtain a spatial resolution of 2 Mm. The error in the resulting flow estimate ultimately depends on the observation time and the number of travel distances used in the inversion.

We compute f-mode sensitivity kernels for flows. Using a two-dimensional model, the scattered wavefield is calculated in the first Born approximation. We test the correctness of the kernels by comparing an exact solution (constant flow), a solution linearized in the flow, and the total integral of the kernel. In practice, the linear approximation is acceptable for flows as large as about 400 m/s.

We compute f-mode travel-time sensitivity kernels for flows. Using a two-dimensional model, we show that it is important to account for several systematic effects, such as the foreshortening and the projection of the velocity vector onto the line of sight. Correcting for these effects is necessary before any data inversion is attempted away from the center of the solar disk.

The solar magnetism is no more considered as a purely superficial phenomenon. The SoHO community has shown that the length of the solar cycle depends on the transition region between radiation and convection. Nevertheless, the internal solar (stellar) magnetism stays poorly known. Starting in 2008, the American instrument HMI/SDO and the European microsatellite PICARD will enrich our view of the Sun-Earth relationship. Thus obtaining a complete MHD solar picture is a clear objective for the next decades and it requires complementary observations of the dynamics of the radiative zone. For that ambitious goal, space prototypes are being developed to improve gravity mode detection. The Sun is unique to progress on the topology of deep internal magnetic fields and to understand the complex mechanisms which provoke photospheric and coronal magnetic changes and possible longer cycles important for human life. We propose the following roadmap in Europe to contribute to this "impressive" revolution in Astronomy and in our Sun-Earth relationship: SoHO (1995-2007), PICARD (2008-2010), DynaMICS (2009-2017) in parallel to SDO (2008-2017) then a world-class mission located at the L1 orbit or above the solar poles.

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.

In Part I we have developed a theory for fitting p-mode Fourier spectra assuming that these spectra have a multi-normal distribution. We showed, using Monte-Carlo simulations, how one can obtain p-mode parameters using 'Maximum Likelihood Estimators'. In this article, hereafter Part II, we show how to use the theory developed in Part I for fitting real data. We introduce 4 new diagnostics in helioseismology: the $(m,\nu)$ echelle diagramme, the cross echelle diagramme, the inter echelle diagramme, and the ratio cross spectrum. These diagnostics are extremely powerful to visualize and understand the covariance matrices of the Fourier spectra, and also to find bugs in the data analysis code. These diagrammes can also be used to derive quantitative information on the mode leakage and noise covariance matrices. Numerous examples using the LOI/SOHO and GONG data are given.

In this article we present our state of the art of fitting helioseismic p-mode spectra. We give a step by step recipe for fitting the spectra: statistics of the spectra both for spatially unresolved and resolved data, the use of Maximum Likelihood estimates, the statistics of the p-mode parameters, the use of Monte-Carlo simulation and the significance of fitted parameters. The recipe is applied to synthetic low-resolution data, similar to those of the LOI, using Monte-Carlo simulations. For such spatially resolved data, the statistics of the Fourier spectrum is assumed to be a multi-normal distribution; the statistics of the power spectrum is \emph{not} a $\chi^{2}$ with 2 degrees of freedom. Results for $l=1$ shows that all parameters describing the p modes can be obtained without bias and with minimum variance provided that the leakage matrix is known. Systematic errors due to an imperfect knowledge of the leakage matrix are derived for all the p-mode parameters.