Search

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., 16 pages, 11 figures, accepted for publication in A&A, title and affiliation updated

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., 10 pages, 10 Figures, Accepted by A&A

Most future surveys designed to discover transiting exoplanets, including TESS and PLATO, will target bright (V3.0 subgiants. The relatively low amount of contamination (defined as the fraction of false positives, 14 page, 6 figures, 2 tables. Accepted in MNRAS

Retrograde-propagating waves of vertical vorticity with longitudinal wavenumbers between 3 and 15 have been observed on the Sun with a dispersion relation close to that of classical sectoral Rossby waves. The observed vorticity eigenfunctions are symmetric in latitude, peak at the equator, switch sign near $20^\circ$-$30^\circ$, and decrease at higher latitudes. We search for an explanation that takes into account solar latitudinal differential rotation. In the equatorial $\beta$ plane, we study the propagation of linear Rossby waves (phase speed $c, Comment: Submitted to Astron. Astrophys. on 29 May 2020

Context. Large-scale equatorial Rossby modes have been observed on the Sun over the last two solar cycles. Aims. We investigate the impact of the time-varying zonal flows on the frequencies of Rossby modes. Methods. A first-order perturbation theory approach is used to obtain an expression for the expected shift in the mode frequencies due to perturbations in the internal rotation rate. Results. Using the time-varying rotation from helioseismic inversions we predict the changes in Rossby mode frequencies with azimuthal orders from m = 1 to m = 15 over the last two solar cycles. The peak-to-peak frequency change is less than 1 nHz for the m = 1 mode, grows with m, and reaches 25 nHz for m = 15. Conclusions. Given the observational uncertainties on mode frequencies due to the finite mode lifetimes, we find that the predicted frequency shifts are near the limit of detectability.

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., 16 pages, 12 figures, accepted for publication by A&A

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., 6 pages, 4 figures, 1 table, Accepted in A&A

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., Submitted to and accepted by Astronomy and Astrophysics; includes minor changes from language editing report

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., 45 pages, 11 figures, 4 tables, published in Science

Context. 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 kernels are computed using the ray or the Born approximation. The Born approximation is preferable as it takes finite-wavelength effects into account, although it can be computationally expensive. Aims. We propose a fast algorithm to compute travel-time sensitivity kernels under the assumption that the background solar medium is spherically symmetric. Methods. 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. Results. 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. Conclusions. This new method allows 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.

International audience; This work offers some contributions to the numerical study of acoustic waves propagating in the Sun and its atmosphere. The main goal is to provide boundary conditions for outgoing waves in the solar atmosphere where it is assumed that the sound speed is constant and the density decays exponentially with radius. Outgoing waves are governed by a Dirichlet-to-Neumann map which is obtained from the factorization of the Helmholtz equation expressed in spherical coordinates. For the purpose of extending the outgoing wave equation to axisymmetric or 3D cases, different approximations are implemented by using the frequency and/or the angle of incidence as parameters of interest. This results in boundary conditions called Atmospheric Radiation Boundary Conditions (ARBC) which are tested in ideal and realistic configurations. These ARBCs deliver accurate results and reduce the computational burden by a factor of two in helioseismology applications.; Ce travail apporte quelques contributions à l'étude numérique des ondes acoustiques se propageant dans le Soleil et son atmosphère. Il se base sur la caractérisa-tion des ondes sortantes dans l'atmosphère représentée par une vitesse constante et une densité décroissant exponentiellement. Les ondes sortantes sont régies par un opérateur Dirichlet-to-Neumann qui est obtenu par la factorisation de l'équation de Helmholtz formulée dans les coordonnées sphériques. Afin d'étendre l'équation des ondes sortantes à des géométries axisymétriques ou 3D, différentes approximations sont menées en util-isant la fréquence et/ou l'angle d'incidence comme paramètres d'intérêt. Ceci mène à des conditions de frontière que nous appelons Conditions de Radiation Atmosphériques (ARBC) et qui sont testées en configuration idéalisées et réalistes. Ces conditions ARBC offrent des résultats précis et réduisent le coût de calcul d'un facteur deux pour le cas du Soleil.

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., Comment: 8 pages, 7 Figures, 1 Table, Accepted to ApJ

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., 15 pages, 11 Figures

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., 18 pages, 10 figures

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., Comment: This is the submitted version of the paper published in Nature Astronomy. 23 pages, 8 figures, 1 table

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., Submitted to Astronomy and Astrophysics

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}$., Comment: accepted for publication in A&A

Context.Previous helioseismology of sunspots has been sensitive to both the structural and magnetic aspects of sunspot structure.Aims.We aim to develop a technique that is insensitive to the magnetic component so the two aspects can be more readily separated.Methods.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.Results.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.Conclusions.With sufficient modeling effort, these should lead to better understanding of sunspot structure.

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., Comment: 8 pages, 4 figures, accepted in A&A Letters

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., accepted for publication in A&A

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., Accepted for publication in Astronomy & Astrophysics. 8 pages, 4 figures, 1 table

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., 24 pages, 8 Figures; MNRAS

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., Accepted for publication in Astronomy and Astrophysics. 16 pages, 10 figures and 8 tables

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., Comment: Accepted on 3 Nov 2016 for publication in Astronomy and Astrophysics

S. Jafarzadeh et. al., ©2017 The American Astronomical Society. All rights reserved.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., The German contribution to Sunrise and its reflight was funded by the Max Planck Foundation, the Strategic Innovations Fund of the President of the Max Planck Society (MPG), DLR, and private donations by supporting members of the Max Planck Society, which is gratefully acknowledged. The Spanish contribution was funded by the Ministerio de Economía y Competitividad under Projects ESP2013-47349-C6 and ESP2014-56169-C6, partially using European FEDER funds. The HAO contribution was partly funded through NASA grant number NNX13AE95G. This work was partly supported by the BK21 plus program through the National Research Foundation (NRF) funded by the Ministry of Education of Korea. S.J. receives support from the Research Council of Norway. T.W. acknowledges support by DFG-grant WI 3211/4-1. M.S. acknowledges support from the European Research Council under the European Union's Seventh Framework Programme / ERC Grant agreement nr. 291058.

©2017 The American Astronomical Society. All rights reserved.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 polarities 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., This work has been partially funded by the Spanish Ministerio de Economía y Competitividad, through Projects No. ESP2013-47349-C6 and ESP2014-56169-C6, including a percentage from European FEDER funds. The German contribution to sunrise and its reflight was funded by the Max Planck Foundation, the Strategic Innovations Fund of the President of the Max Planck Society (MPG), DLR, and private donations by supporting members of the Max Planck Society, which is gratefully acknowledged. The National Solar Observatory (NSO) is operated by the Association of Universities for Research in Astronomy (AURA) Inc. under a cooperative agreement with the National Science Foundation. The HAO contribution was partly funded through NASA grant number NNX13AE95G. This work was partly supported by the BK21 plus program through the National Research Foundation (NRF) funded by the Ministry of Education of Korea.

T. Wiegelmann et. al., ©2017 The American Astronomical Society. All rights reserved.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 2013 June as boundary conditions 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 was applied earlier 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 us to model intrinsic nonlinear structures like strongly localized electric currents., The German contribution to Sunrise and its reflight was funded by the Max Planck Foundation, the Strategic Innovations Fund of the President of the Max Planck Society (MPG), DLR, and private donations by supporting members of the Max Planck Society, which is gratefully acknowledged. The Spanish contribution was funded by the Ministerio de Economía y Competitividad under Projects ESP2013-47349-C6 and ESP2014-56169-C6, partially using European FEDER funds. The HAO contribution was partly funded through NASA grant number NNX13AE95G. This work was partly supported by the BK21 plus program through the National Research Foundation (NRF) funded by the Ministry of Education of Korea. The used HMI-data are courtesy of NASA/SDO and the HMI science team. TW acknowledges DLR-grant 50 OC 1301 and DFG-grant WI 3211/4-1. T.N. acknowledges support by the UK's Science and Technology Facilities Council via Consolidated Grants ST/K000950/1 and ST/N000609/1. D.N. was supported from GA ČR under grant numbers 16-05011S and 16-13277S. The Astronomical Institute Ondřejov is supported by the project RVO:67985815.

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., Comment: 10 pages, 9 figures. Submitted to Astronomy & Astrophysics

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. ..., Comment: Accepted for publication in Astronomy & Astrophysics. 13 pages, 6 figures, 5 tables. Abstract slightly abridged to meet arXiv's 1920 character limit

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., Comment: accepted for publication in A&A

S. Danilovic et. al., ©2017 The American Astronomical Society. All rights reserved.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 analyze Sunrise/IMaX data, along 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. However, the information on the magnetic field vector and velocity at this spatial resolution is extremely valuable because it shows what numerical models miss and how they can be improved., The German contribution to Sunrise and its reflight was funded by the Max Planck Foundation, the Strategic Innovations Fund of the President of the Max Planck Society (MPG), DLR, and private donations by supporting members of the Max Planck Society, which are all gratefully acknowledged. This work has benefited from the discussions at the meeting "Solar UV bursts—a new insight to magnetic reconnection" at the International Space Science Institute (ISSI) in Bern. The Spanish contribution was funded by the Ministerio de Economia y Competitividad under Projects ESP2013-47349-C6 and ESP2014-56169-C6, partially using European FEDER funds. The HAO contribution was partly funded through NASA grant number NNX13AE95G. The National Solar Observatory (NSO) is operated by the Association of Universities for Research in Astronomy (AURA) Inc. under a cooperative agreement with the National Science Foundation. This work was also partly supported by the BK21 plus program through the National Research Foundation (NRF), funded by the Ministry of Education of Korea.

We present results of modelling of α Cen A. In order to estimate the physical parameters of this star, we modelled spectroscopic, interferometric, astrometric and asteroseismic data. To that effect we chose to use a Bayesian approach to parameter estimation, which allowed us, in particular, to define our prior knowledge on the parameters. An important question we wanted to address was to assess whether or not α Cen A has a convective core. We found that the data we used give indecisive results on this issue. If the star has a convective core, and provided that overshooting is taken into account, there is a possibility for the star to be in the peculiar state in which the ppII chain is the main driver of nuclear energy generation. We also found a non-negligible probability for α Cen A to be a very early subgiant.

S. Jafarzadeh et. al., ©2017 The American Astronomical Society. All rights reserved. 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 1.1 Å wide passband covering the line core of Ca ii H 3969 Å 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−1 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−1. 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 Alfvénic 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., The German contribution to Sunrise and its reflight was funded by the Max Planck Foundation, the Strategic Innovations Fund of the President of the Max Planck Society (MPG), DLR, and private donations by supporting members of the Max Planck Society, which is gratefully acknowledged. The Spanish contribution was funded by the Ministerio de Economía y Competitividad under Projects ESP2013-47349-C6 and ESP2014-56169-C6, partially using European FEDER funds. The HAO contribution was partly funded through NASA grant number NNX13AE95G. This work was partly supported by the BK21 plus program through the National Research Foundation (NRF) funded by the Ministry of Education of Korea.

International audience; 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 cutoff frequency (∼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.

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., 11 pages, 5 figures, The Astrophysical Journal Supplement Series (accepted)

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., Comment: accepted for publication in Solar Physics

The level of solar activity varies from cycle to cycle. This variability is probably caused by a combination of nonlinear and random effects. Based on surface flux transport simulations, we show that the observed inflows into active regions and toward the activity belts provide an important nonlinearity in the framework of Babcock-Leighton model for the solar dynamo. Inclusion of these inflows also leads to a reproduction of the observed proportionality between the open heliospheric flux during activity minima and the maximum sunspot number of the following cycle. A substantial component of the random variability of the cycle strength is associated with the cross-equatorial flux plumes that occur when large, highly tilted sunspot groups emerge close to the equator. We show that the flux transported by these events is important for the amplitude of the polar fields and open flux during activity minima. The combined action of inflows and cross-equatorial flux plumes provides an explanation for the weakness of the polar fields at the end of solar cycle 23 (and hence for the relative weakness of solar cycle 24).

Context. The emergence of solar active regions is an important but poorly understood aspect of the solar dynamo.Aims. Knowledge of the flows associated with the rise of active-region-forming magnetic concentrations through the near-surface layers will help determine the mechanisms of active region formation.Methods. We used helioseismic holography and granulation tracking to measure the horizontal flows at the surface that precede the emergence of active regions. We then averaged these flows over about sixty emerging active regions to reduce the noise, selecting active regions that emerge into relatively quiet Sun. To help interpret the results, we constructed a simple model flow field by generating synthetic “emergence locations” that are probabilistically related to the locations of supergranulation-scale convergence regions in the quiet Sun.Results. The flow maps obtained from helioseismology and granulation tracking are very similar (correlation coefficients for single maps around 0.96). We find that active region emergence is, on average, preceded by converging horizontal flows of amplitude about 40 m s−1. The convergence region extends over about 40 Mm in the east-west direction and about 20 Mm in the north-south direction and is centered in the retrograde direction relative to the emergence location. This flow pattern is largely reproduced by a model in which active region emergence occurs preferentially in the prograde direction relative to supergranulation inflows.Conclusions. Averaging over many active regions reveals a statistically significant pattern of near-surface flows prior to emergence. The qualitative success of our simple model suggests that rising flux concentrations and supergranule-scale flows interact during the emergence process.

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., accepted A&A

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., Comment: Accepted in A&A

Context. 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 since starspot latitudes are generally unknown. The related changes in starspot rotation rates caused by latitudinal differential rotation can, however, be measured.Aims. Using the set of 3093Keplerstars with measured activity cycles, we aim to study the temporal change in starspot rotation rates over magnetic activity cycles, and how this relates to the activity level, the mean rotation rate of the star, and its effective temperature.Methods. We measured the photometric variability as a proxy for the magnetic activity and the spot rotation rate in each quarter over the duration of theKeplermission. We phase-folded these measurements with the cycle period. To reduce random errors, we performed averages over stars with comparable mean rotation rates and effective temperature at fixed activity-cycle phases.Results. 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 (rotation 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 observed around rotation periods of 10 − 15 days, where stars with an effective temperature above 4200 K instead show positive correlation.Conclusions. Our measurements can be interpreted in terms of a stellar “butterfly diagram”, but these appear different from that of the Sun since the starspot rotation rates are either in phase or anti-phase with the activity level. Alternatively, the activity cycle periods observed byKeplerare short (around 2.5 years) and may therefore be secondary cycles, perhaps analogous to the solar quasi-biennial oscillations.

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., accepted for publication in A&A

©2017 The American Astronomical Society. All rights reserved. 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 best match the observed profiles. 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 an 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 toward inversions that give results consistent with the MHD equations., The German contribution to Sunrise and its reflight was funded by the Max Planck Foundation, the Strategic Innovations Fund of the President of the Max Planck Society (MPG), DLR, and private donations by supporting members of the Max Planck Society, which is gratefully acknowledged. The Spanish contribution was funded by the Ministerio de Economía y Competitividad under Projects ESP2013-47349-C6 and ESP2014-56169-C6, partially using European FEDER funds. The National Center for Atmospheric Research is sponsored by the National Science Foundation. The HAO contribution was partly funded through NASA grant number NNX13AE95G. The National Solar Observatory (NSO) is operated by the Association of Universities for Research in Astronomy (AURA) Inc. under a cooperative agreement with the National Science Foundation. This work was partly supported by the BK21 plus program through the National Research Foundation (NRF) funded by the Ministry of Education of Korea.

S. Jafarzadeh et. al., ©2017 The American Astronomical Society. All rights reserved.Convective flows are known as the prime means of transporting magnetic fields on the solar surface. Thus, small magnetic structures are good tracers of 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\mbox{--}2.0$), and decelerating movers in the vicinity of flux emergence and/or within active regions ($\gamma =1.4\mbox{--}1.5$). The three types of regions host MBFs with mean diffusion coefficients of 130 km2 s−1, 80–90 km2 s−1, and 25–70 km2 s−1, 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%., The German contribution to Sunrise and its reflight was funded by the Max Planck Foundation, the Strategic Innovations Fund of the President of the Max Planck Society (MPG), DLR, and private donations by supporting members of the Max Planck Society, which is gratefully acknowledged. The Spanish contribution was funded by the Ministerio de Economía y Competitividad under Projects ESP2013-47349-C6 and ESP2014-56169-C6, partially using European FEDER funds. The HAO contribution was partly funded through NASA grant number NNX13AE95G. This work was partly supported by the BK21 plus program through the National Research Foundation (NRF) funded by the Ministry of Education of Korea. S.J. receives support from the Research Council of Norway.

L. P. Chitta et. al., ©2017 The American Astronomical Society. All rights reserved. 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., L.P.C. acknowledges funding by the Max-Planck-Princeton Center for Plasma Physics and funding from the European Union’s Horizon 2020 research and innovation program under Marie Skłodowska-Curie grant agreement No. 707837. The German contribution to Sunrise and its reflight was funded by the Max Planck Foundation, the Strategic Innovations Fund of the President of the Max Planck Society (MPG), DLR, and private donations by supporting members of the Max Planck Society, which is gratefully acknowledged. The Spanish contribution was funded by the Ministerio de Economía y Competitividad under Projects ESP2013-47349-C6 and ESP2014-56169-C6, partially using European FEDER funds. The HAO contribution was partly funded through NASA grant number NNX13AE95G. SDO data are the courtesy of NASA/SDO and the AIA, and HMI science teams. This work was partly supported by the BK21 plus program through the National Research Foundation (NRF) funded by the Ministry of Education of Korea

A. J. Kaithakkal et. al., ©2017 The American Astronomical Society. All rights reserved.Spectropolarimetric observations from Sunrise/IMaX, obtained in 2013 June, 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 blueshifted, whereas 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 ~1.2 × 1017 Mx., The German contribution to Sunrise and its reflight was funded by the Max Planck Foundation, the Strategic Innovations Fund of the President of the Max Planck Society (MPG), DLR, and private donations by supporting members of the Max Planck Society, which is gratefully acknowledged. The Spanish contribution was funded by the Ministerio de Economía y Competitividad under Projects ESP2013-47349-C6 and ESP2014-56169-C6, partially using European FEDER funds. The HAO contribution was partly funded through NASA grant number NNX13AE95G. This work was partly supported by the BK21 plus program through the National Research Foundation (NRF) funded by the Ministry of Education of Korea.