As the amount of asteroseismic data available continues to grow, the inability to accurately model observed oscillation frequencies is becoming a critical problem for interpreting these frequencies. A major component of this problem is the modeling of the near-surface layers. Our aim is to develop a method to estimate the effect of the near-surface layers on oscillation frequencies. In the proposed method we numerically estimate eigenfunctions in 3D hydrodynamic simulations. We match those to the eigenfunctions calculated from the classic equations applied to the horizontal averages of the structure variables. We use this procedure to calculate the frequency perturbation resulting from the dynamical part of the interaction of the oscillations with near-surface convection. As the last step we scale the numbers to the Sun. To provide a qualitative test of our method we performed a series of simulations, calculated the perturbations using our procedure, and compared them to previously reported residuals relative to solar models. We find that we can largely reproduce the observed frequency residuals without resorting to poorly justified theoretical models. We find that, while the calculations of Houdek et al. (2017, MNRAS, 464, L124) produce similar frequency perturbations, the density-pressure phase differences computed here do not match those of that work.

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.

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.

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.

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.

The extension of on-board data processing capabilities is an attractive option to reduce telemetry for scientific instruments on deep space missions. The challenges that this presents, however, require a comprehensive software system, which operates on the limited resources a data processing unit in space allows. We implemented such a system for the Polarimetric and Helioseismic Imager (PHI) on-board the Solar Orbiter (SO) spacecraft. It ensures autonomous operation to handle long command-response times, easy changing of the processes after new lessons have been learned and meticulous book-keeping of all operations to ensure scientific accuracy. This contribution presents the requirements and main aspects of the software implementation, followed by an example of a task implemented in the software frame, and results from running it on SO/PHI. The presented example shows that the different parts of the software framework work well together, and that the system processes data as we expect. The flexibility of the framework makes it possible to use it as a baseline for future applications with similar needs and limitations as SO/PHI.

The pattern of migrating zonal flow bands associated with the solar cycle, known as the torsional oscillation, has been monitored with continuous global helioseismic observations by the Global Oscillations Network Group, together with those made by the Michelson Doppler Imager onboard the Solar and Heliosepheric Observatory and its successor the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory, since 1995, giving us nearly two full solar cycles of observations. We report that the flows now show traces of the mid-latitude acceleration that is expected to become the main equatorward-moving branch of the zonal flow pattern for Cycle 25. Based on the current position of this branch, we speculate that the onset of widespread activity for Cycle 25 is unlikely to be earlier than the middle of 2019.

Context. Recently our ability to study stars using asteroseismic techniques has increased dramatically, largely through the use of space based photometric observations. Work has also been done using ground based spectroscopic observations and more is expected in the near future from the SONG network. Unfortunately, the intensity observations have an inferior signal to noise ratio and details of the observations do not agree with theory, while the data analysis used in the spectroscopic method has often been based on overly simple models of the spectra. Aims. The aim is to improve the reliability of measurements of the parameters of stellar oscillations using spectroscopic observations and to enable the optimal use of the observations. Methods. While previous investigations have used 1D models, I use realistic magnetohydrodynamic simulations, combined with radiative transfer calculations, to estimate the effects of the oscillations on the spectra. I then calculate the visibility of the oscillation modes for a variety of stellar parameters and using various fitting methods. In addition to the methods used in previous investigations, I use a Singular Value Decomposition technique, which allows one to determine the information content available from the spectral perturbations and how that can be expressed most compactly. Finally I describe how the time series obtained may be analyzed. Results. It is shown that it is important to model the visibilities carefully and that the results deviate substantially from models previously used, especially in the presence of rotation. Detailed spectral modeling may be exploited to measure the properties of a larger number of modes than possible using the commonly used cross-correlation method. With moderate rotation, there is as much information in the line shape changes as in the Doppler shift and an outline of how to extract this is given.

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 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 Helioseismic and Magnetic Imager (HMI) instrument is a major component of NASA's Solar Dynamics Observatory (SDO) spacecraft. Since beginning normal science operations on 1 May 2010, HMI has operated with remarkable continuity, e.g. during the more than five years of the SDO prime mission that ended 30 September 2015, HMI collected 98.4% of all possible 45-second velocity maps; minimizing gaps in these full-disk Dopplergrams is crucial for helioseismology. HMI velocity, intensity, and magnetic-field measurements are used in numerous investigations, so understanding the quality of the data is important. We describe the calibration measurements used to track HMI performance and detail trends in important instrument parameters during the mission. Regular calibration sequences provide information used to improve and update the HMI data calibration. The set-point temperature of the instrument front window and optical bench is adjusted regularly to maintain instrument focus, and changes in the temperature-control scheme have been made to improve stability in the observable quantities. The exposure time has been changed to compensate for a 15% decrease in instrument throughput. Measurements of the performance of the shutter and tuning mechanisms show that they are aging as expected and continue to perform according to specification. Parameters of the tunable-optical-filter elements are regularly adjusted to account for drifts in the central wavelength. Frequent measurements of changing CCD-camera characteristics, such as gain and flat field, are used to calibrate the observations. Infrequent expected events, such as eclipses, transits, and spacecraft off-points, interrupt regular instrument operations and provide the opportunity to perform additional calibration. Onboard instrument anomalies are rare and seem to occur quite uniformly in time. The instrument continues to perform very well.

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.

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.

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.

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

In this brief communication we provide the rationale for, and the outcome of the International Astronomical Union (IAU) resolution vote at the XXIX-th General Assembly in Honolulu, Hawaii, in 2015, on recommended nominal conversion constants for selected solar and planetary properties. The problem addressed by the resolution is a lack of established conversion constants between solar and planetary values and SI units: a missing standard has caused a proliferation of solar values (e.g., solar radius, solar irradiance, solar luminosity, solar effective temperature and solar mass parameter) in the literature, with cited solar values typically based on best estimates at the time of paper writing. As precision of observations increases, a set of consistent values becomes increasingly important. To address this, an IAU Working Group on Nominal Units for Stellar and Planetary Astronomy formed in 2011, uniting experts from the solar, stellar, planetary, exoplanetary and fundamental astronomy, as well as from general standards fields to converge on optimal values for nominal conversion constants. The effort resulted in the IAU 2015 Resolution B3, passed at the IAU General Assembly by a large majority. The resolution recommends the use of nominal solar and planetary values, which are by definition exact and are expressed in SI units. These nominal values should be understood as conversion factors only, not as the true solar/planetary properties or current best estimates. Authors and journal editors are urged to join in using the standard values set forth by this resolution in future work and publications to help minimize further confusion.

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.

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.

We quantify the effect of observational spectroscopic and asteroseismic uncertainties on regularised least squares (RLS) inversions for the radial differential rotation of Sun-like and subgiant stars. We first solved the forward problem to model rotational splittings plus the observed uncertainties for models of a Sun-like star, HD 52265, and a subgiant star, KIC 7341231. We randomly perturbed the parameters of the stellar models within the uncertainties of the spectroscopic and asteroseismic constraints and used these perturbed stellar models to compute rotational splittings. We experimented with three rotation profiles: solid body rotation, a step function, and a smooth rotation profile decreasing with radius. We then solved the inverse problem to infer the radial differential rotation profile using a RLS inversion and kernels from the best-fit stellar model.We found that the inversions for Sun-like stars with solar-like radial differential rotation profiles are insensitive to the uncertainties in the stellar models. We found that when the rotation rate below the convection zone is increased to six times that of the surface rotation rate the inferred rotation profile excluded solid body rotation. With the current observational uncertainties, we found that inversions of subgiant stars are sensitive to the uncertainties in the stellar model. Our findings suggest that inversions for the radial differential rotation of subgiant stars would benefit from more tightly constrained stellar models. In Sun-like stars, the insensitivity of the inversions to stellar model uncertainties suggests that it may be possible to perform ensemble inversions for the average radial differential rotation of many stars with a range of stellar types to better constrain the inversions.

We present a comprehensive study of one method for measuring various parameters of global modes of oscillation of the Sun. Using velocity data taken by the Michelson Doppler Imager (MDI), we analyze spherical harmonic degrees l <= 300. Both current and historical methodologies are explained, and the various differences between the two are investigated to determine their effects on global-mode parameters and systematic errors in the analysis. These differences include a number of geometric corrections made during spherical harmonic decomposition; updated routines for generating window functions, detrending timeseries, and filling gaps; and consideration of physical effects such as mode profile asymmetry, horizontal displacement at the solar surface, and distortion of eigenfunctions by differential rotation. We apply these changes one by one to three years of data, and then reanalyze the entire MDI mission applying all of them, using both the original 72-day long timeseries and 360-day long timeseries. We find significant changes in mode parameters, both as a result of the various changes to the processing, as well as between the 72-day and 360-day analyses. We find reduced residuals of inversions for internal rotation, but seeming artifacts remain, such as the peak in the rotation rate near the surface at high latitudes. An annual periodicity in the f-mode frequencies is also investigated.

The XXIXth IAU General Assembly in Honolulu adopted IAU 2015 Resolution B2 on recommended zero points for the absolute and apparent bolometric magnitude scales. The resolution was proposed by the IAU Inter-Division A-G Working Group on Nominal Units for Stellar and Planetary Astronomy after consulting with a broad spectrum of researchers from the astronomical community. Resolution B2 resolves the long-standing absence of an internationally-adopted zero point for the absolute and apparent bolometric magnitude scales. Resolution B2 defines the zero point of the absolute bolometric magnitude scale such that a radiation source with $M_{\rm Bol}$ = 0 has luminosity L$_{\circ}$ = 3.0128e28 W. The zero point of the apparent bolometric magnitude scale ($m_{\rm Bol}$ = 0) corresponds to irradiance $f_{\circ}$ = 2.518021002e-8 W/m$^2$. The zero points were chosen so that the nominal solar luminosity (3.828e26 W) adopted by IAU 2015 Resolution B3 corresponds approximately to $M_{\rm Bol}$(Sun) = 4.74, the value most commonly adopted in recent literature. The nominal total solar irradiance (1361 W/m$^2$) adopted in IAU 2015 Resolution B3 corresponds approximately to apparent bolometric magnitude $m_{\rm bol}$(Sun) = -26.832. Implicit in the IAU 2015 Resolution B2 definition of the apparent bolometric magnitude scale is an exact definition for the parsec (648000/$\pi$ au) based on the IAU 2012 Resolution B2 definition of the astronomical unit.

Astronomers commonly quote the properties of celestial objects in units of parameters for the Sun, Jupiter, or the Earth. The resolution presented here was proposed by the IAU Inter-Division Working Group on Nominal Units for Stellar and Planetary Astronomy and passed by the XXIXth IAU General Assembly in Honolulu. IAU 2015 Resolution B3 adopts a set of nominal solar, terrestrial, and jovian conversion constants for stellar and (exo)planetary astronomy which are defined to be exact SI values. While the nominal constants are based on current best estimates (CBEs; which have uncertainties, are not secularly constant, and are updated regularly using new observations), they should be interpreted as standard values and not as CBEs. IAU 2015 Resolution B3 adopts five solar conversion constants (nominal solar radius, nominal total solar irradiance, nominal solar luminosity, nominal solar effective temperature, and nominal solar mass parameter) and six planetary conversion constants (nominal terrestrial equatorial radius, nominal terrestrial polar radius, nominal jovian equatorial radius, nominal jovian polar radius, nominal terrestrial mass parameter, and nominal jovian mass parameter).

Context. The interaction of solar oscillations with near surface convection is poorly understood. These interactions are likely the cause of several problems in helio- and astero-seismology, including the so-called surface effect and apparently unphysical travel time shifts as a function of center to limb distance. There is thus a clear need for further theoretical understanding and observational tests. Aims. The aim is to determine how the observed modes are affected by the convection. Methods. I use HMI velocity and intensity images to construct k-$\omega$ diagrams showing how the oscillation amplitude and phase depend on the local granulation intensity. Results. There is a clear and significant dependence of the observed properties of the oscillations on the local convection state.

A new fitting methodology is presented which is equally well suited for the estimation of low-, medium-, and high-degree mode parameters from $m$-averaged solar oscillation power spectra of widely differing spectral resolution. This method, which we call the "Windowed, MuLTiple-Peak, averaged spectrum", or WMLTP Method, constructs a theoretical profile by convolving the weighted sum of the profiles of the modes appearing in the fitting box with the power spectrum of the window function of the observing run using weights from a leakage matrix that takes into account both observational and physical effects, such as the distortion of modes by solar latitudinal differential rotation. We demonstrate that the WMLTP Method makes substantial improvements in the inferences of the properties of the solar oscillations in comparison with a previous method that employed a single profile to represent each spectral peak. We also present an inversion for the internal solar structure which is based upon 6,366 modes that we have computed using the WMLTP method on the 66-day long 2010 SOHO/MDI Dynamics Run. To improve both the numerical stability and reliability of the inversion we developed a new procedure for the identification and correction of outliers in a frequency data set. We present evidence for a pronounced departure of the sound speed in the outer half of the solar convection zone and in the subsurface shear layer from the radial sound speed profile contained in Model~S of Christensen-Dalsgaard and his collaborators that existed in the rising phase of Solar Cycle~24 during mid-2010.

We discuss recent observational, theoretical and modeling progress made in understanding the Sun's internal dynamics, including its rotation, meridional flow, convection and overshoot. Over the past few decades, substantial theoretical and observational effort has gone into appreciating these aspects of solar dynamics. A review of these observations, related helioseismic methodology and inference and computational results in relation to these problems is undertaken here.

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.

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 Helioseismic and Magnetic Imager (HMI) began near-continuous full-disk solar measurements on 1 May 2010 from the Solar Dynamics Observatory (SDO). An automated processing pipeline keeps pace with observations to produce observable quantities, including the photospheric vector magnetic field, from sequences of filtergrams. The primary 720s observables were released in mid 2010, including Stokes polarization parameters measured at six wavelengths as well as intensity, Doppler velocity, and the line-of-sight magnetic field. More advanced products, including the full vector magnetic field, are now available. Automatically identified HMI Active Region Patches (HARPs) track the location and shape of magnetic regions throughout their lifetime. The vector field is computed using the Very Fast Inversion of the Stokes Vector (VFISV) code optimized for the HMI pipeline; the remaining 180 degree azimuth ambiguity is resolved with the Minimum Energy (ME0) code. The Milne-Eddington inversion is performed on all full-disk HMI observations. The disambiguation, until recently run only on HARP regions, is now implemented for the full disk. Vector and scalar quantities in the patches are used to derive active region indices potentially useful for forecasting; the data maps and indices are collected in the SHARP data series, hmi.sharp_720s. Patches are provided in both CCD and heliographic coordinates. HMI provides continuous coverage of the vector field, but has modest spatial, spectral, and temporal resolution. Coupled with limitations of the analysis and interpretation techniques, effects of the orbital velocity, and instrument performance, the resulting measurements have a certain dynamic range and sensitivity and are subject to systematic errors and uncertainties that are characterized in this report.

The Very Fast Inversion of the Stokes Vector (VFISV) is a Milne-Eddington spectral line inversion code used to determine the magnetic and thermodynamic parameters of the solar photosphere from observations of the Stokes vector in the 6173 A Fe I line by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). We report on the modifications made to the original VFISV inversion code in order to optimize its operation within the HMI data pipeline and provide the smoothest solution in active regions. The changes either sped up the computation or reduced the frequency with which the algorithm failed to converge to a satisfactory solution. Additionally, coding bugs which were detected and fixed in the original VFISV release, are reported here.

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.

White-light observations by the Solar Dynamics Observatory's Helioseismic and Magnetic Imager of a loop-prominence system occurring in the aftermath of an X-class flare on 2013 May 13 near the eastern solar limb show a linearly polarized component, reaching up to $\sim$20% at an altitude of $\sim$33 Mm, about the maximal amount expected if the emission were due solely to Thomson scattering of photospheric light by the coronal material. The mass associated with the polarized component was 8.2$\times$10$^{14}$ g. At 15 Mm altitude, the brightest part of the loop was 3(+/-0.5)% linearly polarized, only about 20% of that expected from pure Thomson scattering, indicating the presence of an additional unpolarized component at wavelengths near Fe I (617.33 nm), probably thermal emission. We estimated the free electron density of the white-light loop system to possibly be as high as 1.8$\times$10$^{12}$ cm$^{-3}$.

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 report observations of white-light ejecta in the low corona, for two X-class flares on the 2013 May 13, using data from the Helioseismic and Magnetic Imager (HMI) of the Solar Dynamics Observatory. At least two distinct kinds of sources appeared (chromospheric and coronal), in the early and later phases of flare development, in addition to the white-light footpoint sources commonly observed in the lower atmosphere. The gradual emissions have a clear identification with the classical loop-prominence system, but are brighter than expected and possibly seen here in the continuum rather than line emission. We find the HMI flux exceeds the radio/X-ray interpolation of the bremsstrahlung produced in the flare soft X-ray sources by at least one order of magnitude. This implies the participation of cooler sources that can produce free-bound continua and possibly line emission detectable by HMI. One of the early sources dynamically resembles "coronal rain", appearing at a maximum apparent height and moving toward the photosphere at an apparent constant projected speed of 134 $\pm$ 8 $\mathrm{km s^{-1}}$. Not much literature exists on the detection of optical continuum sources above the limb of the Sun by non-coronagraphic instruments, and these observations have potential implications for our basic understanding of flare development, since visible observations can in principle provide high spatial and temporal resolution.

We propose to observe with Kepler an age sequence of nearby uncrowded open clusters. The current Kepler field contains very few (only 4) somewhat distant and/or old clusters. Nearby open clusters are already well characterized from ground- and many space-based observations. Our proposal focuses on mid and upper main-sequence variables (main sequence and pre-main sequence gamma Dor or delta Scuti stars, SPB/beta Cep stars, Cepheids, or yellow supergiants), having periods of hours to days accessible by longer cadence observations. Asteroseismology for these objects to date is limited by the number of modes observable from the ground, difficulty in obtaining spectroscopic or photometric mode ID for the fainter targets that have been observed by Kepler, uncertainties in interior and initial abundances, especially for stars with abundance anomalies, uncertainties in distance/luminosity, and lack of knowledge of prior evolution history. The additional constraints of common age, distance and initial abundances in clusters will place these variables in their evolutionary context, and help unlock some of the science and reduce uncertainties for understanding star formation and stellar evolution.

With the increasing quality of asteroseismic observations it is important to minimize the random and systematic errors in mode parameter estimates. To this end it is important to understand how the oscillations relate to the directly observed quantities, such as intensities and spectra, and to derived quantities, such as Doppler velocity. Here I list some of the effects we need to take into account and show an example of the impact of some of them.

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.

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.

As a response to the white paper call, we propose to turn Kepler to the South Ecliptic Pole (SEP) and observe thousands of large amplitude variables for years with high cadence in the frame of the Kepler-SEP Mission. The degraded pointing stability will still allow observing these stars with reasonable (probably better than mmag) accuracy. Long-term continuous monitoring already proved to be extremely helpful to investigate several areas of stellar astrophysics. Space-based missions opened a new window to the dynamics of pulsation in several class of pulsating variable stars and facilitated detailed studies of eclipsing binaries. The main aim of this mission is to better understand the fascinating dynamics behind various stellar pulsational phenomena (resonances, mode coupling, chaos, mode selection) and interior physics (turbulent convection, opacities). This will also improve the applicability of these astrophysical tools for distance measurements, population and stellar evolution studies. We investigated the pragmatic details of such a mission and found a number of advantages: minimal reprogramming of the flight software, a favorable field of view, access to both galactic and LMC objects. However, the main advantage of the SEP field comes from the large sample of well classified targets, mainly through OGLE. Synergies and significant overlap (spatial, temporal and in brightness) with both ground- (OGLE, LSST) and space-based missions (GAIA, TESS) will greatly enhance the scientific value of the Kepler-SEP mission. GAIA will allow full characterization of the distance indicators. TESS will continuously monitor this field for at least one year, and together with the proposed mission provide long time series that cannot be obtained by other means. If Kepler-SEP program is successful, there is a possibility to place one of the so-called LSST "deep-drilling" fields in this region.

Many helioseismic measurements suffer from substantial systematic errors. A particularly frustrating one is that time-distance measurements suffer from a large center to limb effect which looks very similar to the finite light travel time, except that the magnitude depends on the observable used and can have the opposite sign Zhao et al. (2012). This has frustrated attempts to determine the deep meridional flow in the solar convection zone, with Zhao et al. (2012) applying an ad hoc correction with little physical basis to correct the data. In this letter we propose that part of this effect can be explained by the highly asymmetrical nature of the solar granulation which results in what appears to the oscillation modes as a net radial flow, thereby imparting a phase shift on the modes as a function of observing height and thus heliocentric angle.

We present the first accurate characterization of high-degree modes, derived using the best MDI full-disk full-resolution data set available. A ninety day long time series of full-disk two arc-second per pixel resolution dopplergrams was acquired in 2001. These dopplergrams were spatially decomposed using our best estimate of the image scale and the known components of MDI's image distortion. A multi-taper power spectrum estimator was used to generate power spectra up to l = 1000, with a large number of tapers to reduce the realization noise, while the blending at high degrees negates the need for high spectral resolution. These power spectra were fitted for all degrees and all azimuthal orders, between l = 100 and l = 1000. This fitting generated in excess of 6x10^6 individual estimates of ridge frequencies, line-widths, amplitudes and asymmetries (singlets: l,n,m), corresponding to some 6,000 multiplets (l,n). Fitting at high degrees generates ridge characteristics, characteristics that do not correspond to the underlying mode characteristics. We used a sophisticated forward modeling to recover the best possible estimate of the underlying mode characteristics (mode frequencies, as well as line-widths, amplitudes and asymmetries). We describe in detail this modeling and its validation. The modeling has been extensively reviewed and refined, by including an iterative process to improve its input parameters to better match the observations. Also, the contribution of the leakage matrix on the accuracy of the procedure has been carefully assessed. We present the derived set of corrected mode characteristics, discuss their uncertainties and the precision of the ridge to mode correction schemes. In our conclusions, we address how to further improve these estimates, and the implications for other data sets, like GONG+ and HMI.

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.

We describe observations of a white-light flare (SOL2011-02-24T07:35:00, M3.5) close to the limb of the Sun, from which we obtain estimates of the heights of the optical continuum sources and those of the associated hard X-ray sources.For this purpose we use hard X-ray images from the Reuven Ramaty High Energy Spectroscopic Imager (RHESSI), and optical images at 6173 \AA from the Solar Dynamics Observatory (SDO). We find that the centroids of the impulsive-phase emissions in white light and hard X-rays (30-80 keV) match closely in central distance (angular displacement from Sun center), within uncertainties of order 0.2". This directly implies a common source height for these radiations, strengthening the connection between visible flare continuum formation and the accelerated electrons. We also estimate the absolute heights of these emissions, as vertical distances from Sun center. Such a direct estimation has not been done previously, to our knowledge. Using a simultaneous 195 \AA image from the Solar-Terrestrial RElations Observatory (STEREO-B) spacecraft to identify the heliographic coordinates of the flare footpoints, we determine mean heights above the photosphere (as normally defined; \tau = 1 at 5000 \AA) of 305 \pm 170 km and 195 \pm 70 km, respectively, for the centroids of the hard X-ray (HXR) and white light (WL) footpoint sources of the flare. These heights are unexpectedly low in the atmosphere, and are consistent with the expected locations of \tau = 1 for the 6173 \AA and the ~40 keV photons observed, respectively.

Coronal mass ejections (CMEs) and corotating interaction regions (CIRs) as well as their source regions are important because of their space weather consequences. The current understanding of CMEs primarily comes from the Solar and Heliospheric Observatory (SOHO) and the Solar Terrestrial Relations Observatory (STEREO) missions, but these missions lacked some key measurements: STEREO did not have a magnetograph; SOHO did not have in-situ magnetometer. SOHO and other imagers such as the Solar Mass Ejection Imager (SMEI) located on the Sun-Earth line are also not well-suited to measure Earth-directed CMEs. The Earth-Affecting Solar Causes Observatory (EASCO) is a proposed mission to be located at the Sun-Earth L5 that overcomes these deficiencies. The mission concept was recently studied at the Mission Design Laboratory (MDL), NASA Goddard Space Flight Center, to see how the mission can be implemented. The study found that the scientific payload (seven remote-sensing and three in-situ instruments) can be readily accommodated and can be launched using an intermediate size vehicle; a hybrid propulsion system consisting of a Xenon ion thruster and hydrazine has been found to be adequate to place the payload at L5. Following a 2-year transfer time, a 4-year operation is considered around the next solar maximum in 2025.

The Helioseismic and Magnetic Imager (HMI) has just started producing data that will help determine what the sources and mechanisms of variability in the Sun's interior are. The instrument measures the Doppler shift and the polarization of the Fe I 6173 A line, on the entire solar disk at a relatively-high cadence, in order to study the oscillations and the evolution of the full vector magnetic field of the solar Photosphere. After the data are properly calibrated, they are given to a Milne-Eddington inversion code (VFISV, Borrero et al. 2010) whose purpose is to infer certain aspects of the physical conditions in the Sun's Photosphere, such as the full 3-D topology of the magnetic field and the line-of-sight velocity at the solar surface. We will briefly describe the characteristics of the inversion code, its advantages and limitations --both in the context of the model atmosphere and the actual nature of the data--, and other aspects of its performance on such a remarkable data load. Also, a cross-comparison with near-simultaneous maps from the Spectro-Polarimeter (SP) onboard Hinode will be made.

We report observations of a white-light solar flare (SOL2010-06-12T00:57, M2.0) observed by the Helioseismic Magnetic Imager (HMI) on the Solar Dynamics Observatory (SDO) and the Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI). The HMI data give us the first space-based high-resolution imaging spectroscopy of a white-light flare, including continuum, Doppler, and magnetic signatures for the photospheric FeI line at 6173.34{\AA} and its neighboring continuum. In the impulsive phase of the flare, a bright white-light kernel appears in each of the two magnetic footpoints. When the flare occurred, the spectral coverage of the HMI filtergrams (six equidistant samples spanning \pm172m{\AA} around nominal line center) encompassed the line core and the blue continuum sufficiently far from the core to eliminate significant Doppler crosstalk in the latter, which is otherwise a possibility for the extreme conditions in a white-light flare. RHESSI obtained complete hard X-ray and \Upsilon-ray spectra (this was the first \Upsilon-ray flare of Cycle 24). The FeI line appears to be shifted to the blue during the flare but does not go into emission; the contrast is nearly constant across the line profile. We did not detect a seismic wave from this event. The HMI data suggest stepwise changes of the line-of-sight magnetic field in the white-light footpoints.

The primary science goal of the Kepler Mission is to provide a census of exoplanets in the solar neighborhood, including the identification and characterization of habitable Earth-like planets. The asteroseismic capabilities of the mission are being used to determine precise radii and ages for the target stars from their solar-like oscillations. Chaplin et al. (2010) published observations of three bright G-type stars, which were monitored during the first 33.5 days of science operations. One of these stars, the subgiant KIC 11026764, exhibits a characteristic pattern of oscillation frequencies suggesting that it has evolved significantly. We have derived asteroseismic estimates of the properties of KIC 11026764 from Kepler photometry combined with ground-based spectroscopic data. We present the results of detailed modeling for this star, employing a variety of independent codes and analyses that attempt to match the asteroseismic and spectroscopic constraints simultaneously. We determine both the radius and the age of KIC 11026764 with a precision near 1%, and an accuracy near 2% for the radius and 15% for the age. Continued observations of this star promise to reveal additional oscillation frequencies that will further improve the determination of its fundamental properties.

We present preliminary asteroseismic results from Kepler on three G-type stars. The observations, made at one-minute cadence during the first 33.5d of science operations, reveal high signal-to-noise solar-like oscillation spectra in all three stars: About 20 modes of oscillation may be clearly distinguished in each star. We discuss the appearance of the oscillation spectra, use the frequencies and frequency separations to provide first results on the radii, masses and ages of the stars, and comment in the light of these results on prospects for inference on other solar-type stars that Kepler will observe.

We examine the evolution of the zonal flow pattern in the upper solar convection zone during the current extended solar minimum, and compare it with that during the previous minimum. The results suggest that a configuration matching that at the previous minimum was reached during 2008, but that the flow band corresponding to the new cycle has been moving more slowly towards the equator than was observed in the previous cycle, resulting in a gradual increase in the apparent length of the cycle during the 2007 -- 2008 period. The current position of the lower-latitude fast-rotating belt corresponds to that seen around the onset of activity in the previous cycle.

In this paper we describe in detail the implementation and main properties of a new inversion code for the polarized radiative transfer equation (VFISV: Very Fast inversion of the Stokes vector). VFISV will routinely analyze pipeline data from the Helioseismic and Magnetic Imager (HMI) on-board of the Solar Dynamics Observatory (SDO). It will provide full-disk maps (4096$\times$4096 pixels) of the magnetic field vector on the Solar Photosphere every 10 minutes. For this reason VFISV is optimized to achieve an inversion speed that will allow it to invert 16 million pixels every 10 minutes with a modest number (approx. 50) of CPUs. Here we focus on describing a number of important details, simplifications and tweaks that have allowed us to significantly speed up the inversion process. We also give details on tests performed with data from the spectropolarimeter on-board of the Hinode spacecraft.

Here we present a detailed analysis of solar acoustic mode frequencies and their rotational splittings for modes with degree up to 900. They were obtained by applying spherical harmonic decomposition to full-disk solar images observed by the Michelson Doppler Imager onboard the Solar and Heliospheric Observatory spacecraft. Global helioseismology analysis of high-degree modes is complicated by the fact that the individual modes cannot be isolated, which has limited so far the use of high-degree data for structure inversion of the near-surface layers (r > 0.97 R). In this work, we took great care to recover the actual mode characteristics using a physically motivated model which included a complete leakage matrix. We included in our analysis the following instrumental characteristics: the correct instantaneous image scale, the radial and non-radial image distortions, the effective position angle of the solar rotation axis and a correction to the Carrington elements. We also present variations of the mode frequencies caused by the solar activity cycle. We have analyzed seven observational periods from 1999 to 2005 and correlated their frequency shift with four different solar indices. The frequency shift scaled by the relative mode inertia is a function of frequency alone and follows a simple power law, where the exponent obtained for the p modes is twice the value obtained for the f modes. The different solar indices present the same result.

Using full-disk observations obtained with the Michelson Doppler Imager (MDI) on board the Solar and Heliospheric Observatory (SOHO) spacecraft, we present variations of the solar acoustic mode frequencies caused by the solar activity cycle. High-degree (100 < l < 900) solar acoustic modes were analyzed using global helioseismology analysis techniques over most of solar cycle 23. We followed the methodology described in details in Korzennik, Rabello-Soares and Schou (2004) to infer unbiased estimates of high-degree mode parameters (see also Rabello-Soares, Korzennik and Schou, 2006). We have removed most of the known instrumental and observational effects that affect specifically high-degree modes. We show that the high-degree changes are in good agreement with the medium-degree results, except for years when the instrument was highly defocused. We analyzed and discuss the effect of defocusing on high degree estimation. Our results for high-degree modes confirm that the frequency shift scaled by the relative mode inertia is a function of frequency and it is independent of degree.

We investigate the accuracy to which we can retrieve the solar photospheric magnetic field vector using the Helioseismic and Magnetic Imager (HMI) that will fly onboard of the Solar Dynamics Observatory (SDO) by inverting simulated HMI profiles. The simulated profiles realistically take into account the effects of the photon noise, limited spectral resolution, instrumental polarization modulation, solar p modes and temporal averaging. The accuracy of the determination of the magnetic field vector is studied considering the different operational modes of the instrument.

We present a study of two spectral lines, Fe I 6173 Angstroms and Ni I 6768 Angstroms, that were candidates to be used in the Helioseismic and Magnetic Imager (HMI) for observing Doppler velocity and the vector magnetic field. The line profiles were studied using the Mt. Wilson Observatory, the Advanced Stokes Polarimeter and the Kitt Peak McMath telescope and one meter Fourier transform spectrometer atlas. Both Fe I and Ni I profiles have clean continua and no blends that threaten instrument performance. The Fe I line is 2% deeper, 15% narrower and has a 6% smaller equivalent width than the Ni I line. The potential of each spectral line to recover pre-assigned solar conditions is tested using a least-squares minimization technique to fit Milne-Eddington models to tens of thousands of line profiles that have been sampled at five spectral positions across the line. Overall, the Fe I line has a better performance than the Ni I line for vector magnetic field retrieval. We selected the Fe I spectral line for use in HMI due to its better performance for magnetic diagnostics while not sacrificing velocity information.

We present results of fitting the 50-day time series of photometry of alpha Cen A taken by the WIRE satellite in 1999. Both power spectrum and autocovariance function (ACF) fitting techniques were used in an attempt to determine mode frequencies, rotational splittings, lifetimes and amplitudes of low-l p-modes. In all, using both techniques, we managed to fit 18 modes (seven l = 0, eight l = 1 and three l = 2) with frequencies determined to within 1 - 2 micro-Hz. These estimates are shown to be 0.6 +/- 0.3 micro-Hz lower, on average, than the frequencies determined from two other more recent studies (Bouchy & Carrier, 2002; Bedding et al. 2004) which used data gathered about 19 months after the WIRE observations. This could be indicative of an activity cycle, although due to the large uncertainty, more data would be needed to confirm this. Over a range of 1700 to 2650 micro-Hz we were also able to use the ACF fitting to determine an average lifetime of 3.9 +/- 1.4 days, and an average rotational splitting of 0.54 +/- 0.22 micro-Hz, which is the first ever reliable estimate of this parameter. In contrast to the ACF, the power spectrum fitting was shown to return significantly biased results for these parameters.

The Helioseismic and Magnetic Imager (HMI), on board the Solar Dynamics Observatory (SDO), will begin data acquisition in 2008. It will provide the first full disk, high temporal cadence observations of the full Stokes vector with a 0.5 arc sec pixel size. This will allow for a continuous monitoring of the Solar magnetic field vector. HMI data will advance our understanding of the small and large-scale magnetic field evolution, its relation to the solar and global dynamic processes, coronal field extrapolations, flux emergence, magnetic helicity and the nature of the polar magnetic fields. We summarize HMI's expected operation modes, focusing on the polarization cross-talk induced by the solar oscillations and how this affects the magnetic field vector determinations.

We argue that the non-detection of oscillations in Procyon by the MOST satellite reported by Matthews et al. (2004) is fully consistent with published ground-based velocity observations of this star. We also examine the claims that the MOST observations represent the best photometric precision so far reported in the literature by about an order of magnitude and are the most sensitive data set for asteroseismology available for any star other than the Sun. These statements are not correct, with the most notable exceptions being observations of oscillations in alpha Cen A that are far superior. We further disagree that the hump of excess power seen repeatedly from velocity observations of Procyon can be explained as an artefact caused by gaps in the data. The MOST observations failed to reveal oscillations clearly because their noise level is too high, possibly from scattered Earthlight in the instrument. We did find an excess of strong peaks in the MOST amplitude spectrum that is inconsistent with a simple noise source such as granulation, and may perhaps indicate oscillations at roughly the expected level.

The characteristic of the solar acoustic spectrum is such that mode lifetimes get shorter and spatial leaks get closer in frequency as the degree of a mode increases for a given order. A direct consequence of this property is that individual p-modes are only resolved at low and intermediate degrees, and that at high degrees, individual modes blend into ridges. Once modes have blended into ridges, the power distribution of the ridge defines the ridge central frequency and it will mask the true underlying mode frequency. An accurate model of the amplitude of the peaks that contribute to the ridge power distribution is needed to recover the underlying mode frequency from fitting the ridge. We present the results of fitting high degree power ridges (up to l = 900) computed from several two to three-month-long time-series of full-disk observations taken with the Michelson Doppler Imager (MDI) on-board the Solar and Heliospheric Observatory between 1996 and 1999. We also present a detailed discussion of the modeling of the ridge power distribution, and the contribution of the various observational and instrumental effects on the spatial leakage, in the context of the MDI instrument. We have constructed a physically motivated model (rather than some ad hoc correction scheme) resulting in a methodology that can produce an unbiased determination of high-degree modes, once the instrumental characteristics are well understood. Finally, we present changes in high degree mode parameters with epoch and thus solar activity level and discuss their significance.

Recently Gizon, Duvall and Schou (2002) suggested that supergranulation has a wave-like component. In this paper I show that the same phenomenon can be observed using surface Doppler shift data, thereby confirming their observations. I am also able to measure the dispersion relation to lower wavenumbers and to extend the results for rotation and meridional flows beyond +/-70 degrees latitude.

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.

We study the variation of the frequency splitting coefficients describing the solar asphericity in both GONG and MDI data, and use these data to investigate temporal sound-speed variations as a function of both depth and latitude during the period from 1995-2000 and a little beyond. The temporal variations in even splitting coefficients are found to be correlated to the corresponding component of magnetic flux at the solar surface. We confirm that the sound-speed variations associated with the surface magnetic field are superficial. Temporally averaged results show a significant excess in sound speed around 0.92 solar radii and latitude of 60 degrees.

We have analyzed the full set of SOHO/MDI f- and p-mode oscillation frequencies from 1996 to date in a search for evidence of solar radius evolution during the rising phase of the current activity cycle. Like Antia et al. (2000), we find that a significant fraction of the f-mode frequency changes scale with frequency; and that if these are interpreted in terms of a radius change, it implies a shrinking sun. Our inferred rate of shrinkage is about 1.5 km/y, which is somewhat smaller than found by Antia et al. We argue that this rate does not refer to the surface, but rather to a layer extending roughly from 4 to 8 Mm beneath the visible surface. The rate of shrinking may be accounted for by an increasing radial component of the rms random magnetic field at a rate that depends on its radial distribution. If it were uniform, the required field would be ~7 kG. However, if it were inwardly increasing, then a 1 kG field at 8 Mm would suffice. To assess contribution to the solar radius change arising above 4Mm, we analyzed the p-mode data. The evolution of the p-mode frequencies may be explained by a magnetic^M field growing with activity. The implications of the near-surface magnetic field changes depend on the anisotropy of the random magnetic field. If the field change is predominantly radial, then we infer an additional shrinking at a rate between 1.1-1.3 km/y at the photosphere. If on the other hand the increase is isotropic, we find a competing expansion at a rate of 2.3 km/y. In any case, variations in the sun's radius in the activity cycle are at the level of 10^{-5} or less, hence have a negligible contribution to the irradiance variations.

Helioseismic observations have detected small temporal variations of the rotation rate below the solar surface corresponding to the so-called `torsional oscillations' known from Doppler measurements of the surface. These appear as bands of slower and faster than average rotation moving equatorward. Here we establish, using complementary helioseismic observations over four years from the GONG network and from the MDI instrument on board SOHO, that the banded flows are not merely a near-surface phenomenon: rather they extend downward at least 60 Mm (some 8% of the total solar radius) and thus are evident over a significant fraction of the nearly 200 Mm depth of the solar convection zone.

Advances in hydrodynamical simulations have provided new insights into the effects of convection on the frequencies of solar oscillations. As more accurate observations become available, this may lead to an improved understanding of the dynamics of convection and the interaction between convection and pulsation (\cite{Rosenthal+99}). Recent high-resolution observations from the SOI/MDI instrument on the SOHO spacecraft have provided the so-far most-detailed observations of high-degree modes of solar oscillations, which are particularly sensitive to the near-surface properties of the Sun. Here we present preliminary results of a comparison between these observations and frequencies computed for models based on realistic simulations of near-surface convection. Such comparisons may be expected to help in identifying the causes of the remaining differences between the observed frequencies and those of solar models.

Observations of the Sun with the LOWL instrument provide a homogeneous set of solar p-mode frequencies from low to intermediate degree which allows one to determine the structure of much of the solar interior avoiding systematic errors that are introduced when different data sets are combined, i.e., principally the effects of solar cycle changes on the frequencies. Unfortunately, the LOWL data set contains very few of the lowest-degree modes, which are essential for determining reliably the structure of the solar core - in addition, these lowest-degree data have fairly large associated uncertainties. However, observations made by the Birmingham Solar-Oscillations Network (BiSON) in integrated sunlight provide high-accuracy measurements of a large number of low-degree modes. In this paper we demonstrate that the low-degree mode set of the LOWL data can be successfully combined with the more accurate BiSON data, provided the observations are contemporaneous for those frequencies where the solar-cycle-induced effects are important. We show that this leads to a factor-of-two decrease in the error on the inferred sound speed in the solar core. We find that the solar sound speed is higher than in solar models for r < 0.2Rsun. The density of the solar core is, however, lower than that in solar models.