Time-distance helioseismology is a set of powerful tools to study features below the Sun's surface. Inverse methods are needed to interpret time-distance measurements, with many examples in the literature. However, techniques that utilize a more statistical approach to inferences, and broadly used in the astronomical community, are less-commonly found in helioseismology. This article aims to introduce a potentially powerful inversion scheme based on Bayesian probability theory and Monte Carlo sampling that is suitable for local helioseismology. We describe the probabilistic method and how it is conceptually different from standard inversions used in local helioseismology. Several example calculations are carried out to compare and contrast the setup of the problems and the results that are obtained. The examples focus on two important phenomena studied with helioseismology: meridional circulation and supergranulation. Numerical models are used to compute synthetic observations, providing the added benefit of knowing the solution against which the results can be tested. For demonstration purposes, the problems are formulated in two and three dimensions, using both ray- and Born-theoretical approaches. The results seem to indicate that the probabilistic inversions not only find a better solution with much more realistic estimation of the uncertainties, but they also provide a broader view of the range of solutions possible for any given model, making the interpretation of the inversion more quantitative in nature. Unlike the progress being made in fundamental measurement schemes in local helioseismology that image the far side of the Sun, or have detected signatures of global Rossby waves, among many others, inversions of those measurements have had significantly less success. Such statistical methods may help overcome some of these barriers to move the field forward.

The objective of this work is to determine what fraction of red-giant (RG) stars shows photometric rotational modulation, and understand its origin. One of the underlying questions is the role of close binarity in this population, standing upon the fact that RGs in short-period binary systems (<150 days or so) have been observed to display strong rotational modulation. We select a sample of about 4500 relatively bright RGs observed by Kepler, and show that 370 of them (8%) display rotational modulation. Almost all have oscillation amplitudes below the median of the sample, while 30 of them are not oscillating at all. Of the 85 of these RGs with rotational modulation chosen for follow-up radial-velocity observation and analysis, 34 show clear evidence of spectroscopic binarity. Surprisingly, 26 of the 30 non-oscillators are in this group of binaries. To the contrary, about 85% of the active RGs with detectable oscillations are not part of close binaries. With the help of stellar masses and evolutionary states computed from the oscillation properties, it appears that low-mass red-giant branch stars tend to be magnetically inactive, while intermediate-mass ones tend to be highly active. The opposite trends are true for helium-core burning (red clump) stars, whereby the lower-mass clump stars are comparatively more active and the higher-mass ones less so. In other words, we find that low-mass red-giant branch stars gain angular momentum as they evolve to clump stars, while higher-mass ones lose angular momentum. The trend observed with low-mass stars leads to possible scenarios of planet engulfment or other merging events during the shell-burning phase. Regarding intermediate-mass stars, the rotation periods are long with respect to theoretical expectations reported in the literature, which reinforces the existence of an unidentified sink of angular momentum after the main sequence.

Asteroseismology offers the possibility of probing stellar interiors and testing evolutionary and seismic models. Precise photometry and spectroscopy obtained during multi-site campaigns on young open clusters allows discovering rich samples of pulsating stars and using them in a simultaneous seismic modelling called ensemble asteroseismology. The aim of this study is to obtain the age of the open cluster NGC 6910 by means of ensemble asteroseismology of the early-type pulsating members, to derive their stellar parameters, and to classify the excited modes. We used time-series analysis, performed photometric and spectroscopic mode identification, and calculated grids of evolutionary and seismic models to apply the procedure of ensemble asteroseismology for nine pulsating members of NGC 6910. With two iterations of the procedure of ensemble asteroseismology applied to nine pulsating stars we derived an age of 10.6$^{+0.9}_{-0.8}$ Myr for NGC 6910. Of the nine pulsating stars examined in the paper, eight are $\beta$ Cep stars, including three that are hybrid $\beta$ Cep and slowly pulsating B-type (SPB) pulsators, and one is an SPB star. Interestingly, the least massive $\beta$ Cep star, NGC 6910-38, has a mass of about 5.6 M$_\odot$. The present theory does not predict unstable $p$ modes in B-type stars with such a low mass. The $g$ modes with relatively high frequencies ($>3.5$ d$^{-1}$), observed in three members of the cluster, are also stable according to seismic modelling. Both findings pose a challenge for theoretical calculations and prompt a revision of the opacities. The procedure of ensemble asteroseismology was found to be successful for NGC 6910 and $\chi$ Per on the basis of pulsating B-type stars and can therefore be applied to other young open clusters that are rich in such stars.

Spectroscopic Eclipsing Binaries (SEBs) are fundamental benchmarks in stellar astrophysics and today are observed in breathtaking detail by missions like TESS, Kepler, and APOGEE. We develop a methodology for simultaneous analysis of high precision Kepler light curves and high resolution near-IR spectra from APOGEE and present orbital solutions and evolutionary histories for a subset of SEBs within this overlap. Radial velocities extracted from APOGEE spectra using the Broadening Function technique are combined with Kepler light curves and to yield binary orbital solutions. The Broadening Function approach yields more precise radial velocities than the standard Cross-Correlation Function, which in turn yields more precise orbital parameters and enables the identification of tertiary stars. The orbital periods of these seven SEBs range from 4 to 40 days. Four of the systems (KIC 5285607, KIC 6864859, KIC 6778289, and KIC 4285087) are well-detached binaries. The remaining three systems have apparent tertiary companions, but each exhibits two eclipses along with at least one spectroscopically varying component (KIC 6449358, KIC 6131659, and KIC 6781535). Gaia distances are available for four targets which we use to estimate temperatures of both members of these SEBs. We explore evolutionary histories in H-R diagram space and estimate ages for this subset of our sample. Finally, we consider the implications for the formation pathways of close binary systems via interactions with tertiary companions. Our methodology combined with the era of big data and observation overlap opens up the possibility of discovering and analyzing large numbers of diverse SEBs, including those with high flux ratios and those in triple systems.

Asteroseismology is the only observational tool in astronomy that can probe the interiors of stars, and is a benchmark method for deriving fundamental properties of stars and exoplanets. Over the coming decade, space-based and ground-based observations will provide a several order of magnitude increase of solar-like oscillators, as well as a dramatic increase in the number and quality of classical pulsator observations, providing unprecedented possibilities to study stellar physics and galactic stellar populations. In this white paper, we describe key science questions and necessary facilities to continue the asteroseismology revolution into the 2020's.

In this second paper of the series on internal gravity waves (IGWs), we present a study of the generation and propagation of IGWs in a model solar atmosphere with diverse magnetic conditions. A magnetic field free, and three magnetic models that start with an initial, vertical, homogeneous field of 10 G, 50 G, and 100 G magnetic flux density, are simulated using the CO$^{\rm 5}$BOLD code. We find that the IGWs are generated in similar manner in all four models in spite of the differences in the magnetic environment. The mechanical energy carried by IGWs is significantly larger than that of the acoustic waves in the lower part of the atmosphere, making them an important component of the total wave energy budget. The mechanical energy flux (10$^{6}$ - 10$^{3}$ W m$^{-2}$) is few orders of magnitude larger than the Poynting flux (10$^{3}$ - 10$^{1}$ W m$^{-2}$). The Poynting fluxes show a downward component in the frequency range corresponding to the IGWs, which confirm that these waves do not propagate upwards in the atmosphere when the fields are predominantly vertical and strong. We conclude that, in the upper photosphere, the propagation properties of IGWs depend on the average magnetic field strength and therefore these waves can be potential candidate for magnetic field diagnostics of these layers. However, their subsequent coupling to Alfv\'{e}nic waves are unlikely in a magnetic environment permeated with predominantly vertical fields and therefore they may not directly or indirectly contribute to the heating of layers above plasma-$\beta$ less than 1.

NGC 6910 is the northern hemisphere open cluster known to be rich in $\beta$ Cephei-type stars. Using four-season photometry obtained in Bia\l k\'ow (Poland) and Xinglong (China) observatories, we performed variability survey of NGC 6910. As the result, we found over 100 variable stars in the field of the cluster, including many stars showing variability due to pulsations and binarity. Thanks to the spectroscopic observations, we also detected changes in the profiles of spectral lines of $\beta$ Cep stars, caused by pulsations.

Recent observations of Jupiter have suggested the existence of global oscillatory modes at millihertz frequencies, yet the source mechanism responsible for driving these modes is still unknown. However, the energies necessary to produce observable surface oscillations have been predicted. Here we investigate if moist convection in Jupiter's upper atmosphere can be responsible for driving the global oscillations and what moist convective energy requirements are necessary to achieve these theoretical mode energies and surface amplitudes. We begin by creating a one-dimensional moist convective cloud model and find that the available kinetic energy of the rising cloud column falls below theoretical estimates of oscillations energies. That is, mode excitation cannot occur with a single storm eruption. We then explore stochastic excitation scenarios of the oscillations by moist convective storms. We find that mode energies and amplitudes can reach theoretical estimates if the storm energy available to the modes is more than just kinetic. In order for the modes to be excited, we find that they require $5 \times 10^{27}$ to 10$^{28}$ erg per day. However, even for a large storm eruption each day, the available kinetic energy from the storms falls two orders of magnitude short of the required driving energy. Although our models may oversimplify the true complexity of the coupling between Jovian storms and global oscillations, our findings reveal that enough thermal energy is associated with moist convection to drive the modes, should it be available to them.

Asteroseismology is a powerful tool to measure the fundamental properties of stars and probe their interiors. This is particularly efficient for red giants because their modes are well detectable and give information on their deep layers. However, the seismic relations used to infer the mass and radius of a star have been calibrated on the Sun. Therefore, it is crucial to assess their accuracy for red giants which are not perfectly homologous to it. We study eclipsing binaries with a giant component to test their validity. We identified 16 systems for which we intend to compare the dynamical masses and radii obtained by combined photometry and spectroscopy to the values obtained from asteroseismology. In the present work, we illustrate our approach on a system from our sample.

The solar meridional flow is a crucial ingredient in modern dynamo theory. Seismic estimates of this flow have, however, been contradictory in deeper layers below about $0.9\,R_\odot$. Results from time-distance helioseismology have so far been obtained using the ray approximation. Here, we perform inversions using the Born approximation. The initial result is similar to the result previously obtained by Jackiewicz et al. (2015) using ray kernels while using the same set of GONG data and the SOLA inversion technique. However, we show that the assumption of uncorrelated measurements used in earlier studies may lead to inversion errors being underestimated by a factor of about two to four. In a second step, refined inversions are performed using the full covariance matrix and a regularization for cross-talk. As the results are found to depend on the threshold used in the singular value decomposition, they were obtained for a medium threshold ($10^{-7}-10^{-5}$, about 50% of the values used) and a threshold lower by a factor of 10 (about 70% of the values used). The result obtained with the medium threshold is again similar to the original, with less latitudinal variation. However, using the lower threshold, the inverted flow in the southern hemisphere shows two or three cells stacked radially depending on the associated radial flows. Both the single-cell and the multi-cell profiles are consistent with the measured travel times. All our results confirm a shallow return flow at about $0.9\,R_\odot$.

Sequential chromospheric brightenings (SCBs) are often observed in the immediate vicinity of erupting flares and are associated with coronal mass ejections. Since their initial discovery in 2005, there have been several subsequent investigations of SCBs. These studies have used differing detection and analysis techniques, making it difficult to compare results between studies. This work employs the automated detection algorithm of Kirk et al. (Solar Phys. 283, 97, 2013) to extract the physical characteristics of SCBs in 11 flares of varying size and intensity. We demonstrate that the magnetic substructure within the SCB appears to have a significantly smaller area than the corresponding H-alpha emission. We conclude that SCBs originate in the lower corona around 0.1 R_sun above the photosphere, propagate away from the flare center at speeds of 35 - 85 km/s, and have peak photosphere magnetic intensities of 148 +/- 2.9 G. In light of these measurements, we infer SCBs to be distinctive chromospheric signatures of erupting coronal mass ejections.

The chromosphere is a complex region that acts as an intermediary between the magnetic flux emergence in the photosphere and the magnetic features seen in the corona. Large eruptions in the chromosphere of flares and filaments are often accompanied by ejections of coronal mass off the sun. Several studies have observed fast-moving progressive trains of compact bright points (called Sequential Chromospheric Brightenings or SCBs) streaming away from chromospheric flares that also produce a coronal mass ejection (CME). In this work, we review studies of SCBs and search for commonalties between them. We place these findings into a larger context with contemporary chromospheric and coronal observations. SCBs are fleeting indicators of the solar atmospheric environment as it existed before their associated eruption. Since they appear at the very outset of a flare eruption, SCBs are good early indication of a CME measured in the chromosphere.

Accurate measurements of deep solar meridional flow are of vital interest for understanding the solar dynamo. In this paper, we validate a recently developed method for obtaining sensitivity functions (kernels) for travel-time measurements to solar interior flows using the Born approximation in spherical geometry, which is expected to be more accurate than the classical ray approximation. Furthermore, we develop a numerical approach to efficiently compute a large number of kernels based on the separability of the eigenfunctions into their horizontal and radial dependence. The validation is performed using a hydrodynamic simulation of linear wave propagation in the Sun, which includes a standard single-cell meridional flow profile. We show that, using the Born approximation, it is possible to accurately model observational quantities relevant for time-distance helioseismology such as the mean power spectrum, disc-averaged cross-covariance functions, and travel times in the presence of a flow field. In order to closely match the model to observations, we show that it is beneficial to use mode frequencies and damping rates which were extracted from the measured power spectrum. Furthermore, the contribution of the radial flow to the total travel time is found to reach 20% of the contribution of the horizontal flow at travel distances over $40^\circ$. Using the Born kernels and a 2D SOLA inversion of travel times, we can recover most features of the input meridional flow profile. The Born approximation is thus a promising method for inferring large-scale solar interior flows.

The $\kappa$-mechanism has been successful in explaining the origin of observed oscillations of many types of "classical" pulsating variable stars. Here we examine quantitatively if that same process is prominent enough to excite the potential global oscillations within Jupiter, whose energy flux is powered by gravitational collapse rather than nuclear fusion. Additionally, we examine whether external radiative forcing, i.e. starlight, could be a driver for global oscillations in hot Jupiters orbiting various main-sequence stars at defined orbital semimajor axes. Using planetary models generated by the Modules for Experiments in Stellar Astrophysics (MESA) and nonadiabatic oscillation calculations, we confirm that Jovian oscillations cannot be driven via the $\kappa$-mechanism. However, we do show that in hot Jupiters oscillations can likely be excited via the suppression of radiative cooling due to external radiation given a large enough stellar flux and the absence of a significant oscillatory damping zone within the planet. This trend seems to not be dependent on the planetary mass. In future observations we can thus expect that such planets may be pulsating, thereby giving greater insight into the internal structure of these bodies.

Observations of the solar atmosphere show that internal gravity waves are generated by overshooting convection, but are suppressed at locations of magnetic flux, which is thought to be the result of mode conversion into magneto-acoustic waves. Here, we present a study of the acoustic-gravity wave spectrum emerging from a realistic, self-consistent simulation of solar (magneto-)convection. A magnetic field free, hydrodynamic simulation and a magneto-hydrodynamic (MHD) simulation with an initial, vertical, homogeneous field of 50 G flux density were carried out and compared with each other to highlight the effect of magnetic fields on the internal gravity wave propagation in the Sun's atmosphere. We find that the internal gravity waves are absent or partially reflected back into the lower layers in the presence of magnetic fields and argue that the suppression is due to the coupling of internal gravity waves to slow magneto acoustic waves still within the high-beta region of the upper photosphere. The conversion to Alfven waves is highly unlikely in our model because there is no strongly inclined magnetic field present. We argue that the suppression of internal waves observed within magnetic flux-concentrations may also be due to non-linear breaking of internal waves due to vortex flows that are ubiquitously present in the upper photosphere and the chromosphere.

Starting in December 2014, Kepler K2 observed Neptune continuously for 49 days at a 1-minute cadence. The goals consisted of studying its atmospheric dynamics (Simon et al. 2016), detecting its global acoustic oscillations (Rowe et al., submitted), and those of the Sun, which we report on here. We present the first indirect detection of solar oscillations in intensity measurements. Beyond the remarkable technical performance, it indicates how Kepler would see a star like the Sun. The result from the global asteroseismic approach, which consists of measuring the oscillation frequency at maximum amplitude "nu_max" and the mean frequency separation between mode overtones "Delta nu", is surprising as the nu_max measured from Neptune photometry is larger than the accepted value. Compared to the usual reference nu_max_sun = 3100 muHz, the asteroseismic scaling relations therefore make the solar mass and radius appear larger by 13.8 +/- 5.8 % and 4.3 +/- 1.9 % respectively. The higher nu_max is caused by a combination of the value of nu_max_sun, being larger at the time of observations than the usual reference from SOHO/VIRGO/SPM data (3160 +/- 10 muHz), and the noise level of the K2 time series, being ten times larger than VIRGO's. The peak-bagging method provides more consistent results: despite a low signal-to-noise ratio (SNR), we model ten overtones for degrees l=0,1,2. We compare the K2 data with simultaneous SOHO/VIRGO/SPM photometry and BiSON velocity measurements. The individual frequencies, widths, and amplitudes mostly match those from VIRGO and BiSON within 1 sigma, except for the few peaks with lowest SNR.

We present phase curves of Jupiter from 0-140 degrees as measured in multiple optical bandpasses by Cassini/ISS during the Millennium flyby of Jupiter in late 2000 to early 2001. Phase curves are of interest for studying the energy balance of Jupiter and understanding the scattering behavior of Jupiter as an exoplanet analog. We find that Jupiter is significantly darker at partial phases than an idealized Lambertian planet by roughly 25% and is not well fit by Jupiter-like exoplanet atmospheric models across all wavelengths. We provide analytic fits to Jupiter's phase function in several Cassini/ISS imaging filter bandpasses. In addition, these observations show that Jupiter's color is more variable with phase angle than predicted by models. Therefore, the color of even a near Jupiter-twin planet observed at a partial phase cannot be assumed to be comparable to that of Jupiter at full phase. We discuss how WFIRST and other future direct-imaging missions can enhance the study of cool giants.

Given the potential of ensemble asteroseismology for understanding fundamental properties of large numbers of stars, it is critical to determine the accuracy of the scaling relations on which these measurements are based. From several powerful validation techniques, all indications so far show that stellar radius estimates from the asteroseismic scaling relations are accurate to within a few percent. Eclipsing binary systems hosting at least one star with detectable solar-like oscillations constitute the ideal test objects for validating asteroseismic radius and mass inferences. By combining radial-velocity measurements and photometric time series of eclipses, it is possible to determine the masses and radii of each component of a double-lined spectroscopic binary. We report the results of a four-year radial-velocity survey performed with the \'echelle spectrometer of the Astrophysical Research Consortium's 3.5-m telescope and the APOGEE spectrometer at Apache Point Observatory. We compare the masses and radii of 10 red giants obtained by combining radial velocities and eclipse photometry with the estimates from the asteroseismic scaling relations. We find that the asteroseismic scaling relations overestimate red-giant radii by about 5% on average and masses by about 15% for stars at various stages of red-giant evolution. Systematic overestimation of mass leads to underestimation of stellar age, which can have important implications for ensemble asteroseismology used for Galactic studies. As part of a second objective, where asteroseismology is used for understanding binary systems, we confirm that oscillations of red giants in close binaries can be suppressed enough to be undetectable, an hypothesis that was proposed in a previous work.

The NASA Kepler spacecraft data revealed a large number of multimode nonradially pulsating gamma Dor and delta Sct variable star candidates. The high precision long time-series photometry makes it possible to study amplitude variations of the frequencies. We summarize recent literature on amplitude and frequency variations in pulsating variables. We are searching for amplitude variability in several dozen faint gamma Doradus or delta Scuti variable-star candidates observed as part of the Kepler Guest Observer program. We apply several methods, including a Matlab-script wavelet analysis developed by J. Jackiewicz, and the wavelet technique of the VSTAR software (http://www.aavso.org/vstar-overview). Here we show results for two stars, KIC 2167444 and KIC 2301163. We discuss the magnitude and timescale of the amplitude variations, and the presence or absence of correlations between amplitude variations for different frequencies of a given star. Amplitude variations may be detectable using Kepler data even for stars with Kepler magnitude >14 with low-amplitude frequencies (around 100 ppm) using only one or a few quarters of long-cadence data. We discuss proposed causes of amplitude spectrum variability that will require further investigation.

We combine Kepler photometry with ground-based spectra to present a comprehensive dynamical model of the double red giant eclipsing binary KIC 9246715. While the two stars are very similar in mass (M1 = 2.171 [+0.006 / -0.008], M2 = 2.149 [+0.006 / -0.008] Msun) and radius (R1 = 8.37 [+0.03 / -0.07], R2 = 8.30 [+0.04 / -0.03] Rsun), an asteroseismic analysis finds one main set of solar-like oscillations with unusually low-amplitude, wide modes. A second set of oscillations from the other star may exist, but this marginal detection is extremely faint. Because the two stars are nearly twins, KIC 9246715 is a difficult target for a precise test of the asteroseismic scaling relations, which yield M = 2.17 +/- 0.14 Msun and R = 8.26 +/- 0.18 Rsun. Both stars are consistent with the inferred asteroseismic properties, but we suspect the main oscillator is Star 2 because it is less active than Star 1. We find evidence for stellar activity and modest tidal forces acting over the 171-day eccentric orbit, which are likely responsible for the essential lack of solar-like oscillations in one star and weak oscillations in the other. Mixed modes indicate the main oscillating star is on the secondary red clump (a core-He-burning star), and stellar evolution modeling supports this with a coeval history for a pair of red clump stars. This system is a useful case study and paves the way for a detailed analysis of more red giants in eclipsing binaries, an important benchmark for asteroseismology.

Meridional flow is thought to play a very important role in the dynamics of the solar convection zone; however, because of its relatively small amplitude, precisely measuring it poses a significant challenge. Here we present a complete time-distance helioseismic analysis of about two years of ground-based GONG Doppler data to retrieve the meridional circulation profile for modest latitudes, in an attempt to corroborate results from other studies. We use an empirical correction to the travel times due to an unknown center-to-limb systematic effect. The helioseismic inversion procedure is first tested and reasonably validated on artificial data from a large-scale numerical simulation, followed by a test to broadly recover the solar differential rotation found from global seismology. From GONG data, we measure poleward photospheric flows at all latitudes with properties that are comparable with earlier studies, and a shallow equatorward flow about $65$\,Mm beneath the surface, in agreement with recent findings from HMI data. No strong evidence of multiple circulation cells in depth nor latitude is found, yet the whole phase space has not yet been explored. Tests of mass flux conservation are then carried out on the inferred GONG and HMI flows and compared to a fiducial numerical baseline from models, and we find that the continuity equation is poorly satisfied. While the two disparate data sets do give similar results for about the outer $15\%$ of the interior radius, the total inverted circulation pattern appears to be unphysical in terms of mass conservation when interpreted over modest time scales. We can likely attribute this to both the influence of realization noise and subtle effects in the data and measurement procedure.

The subsurface structure of an "average" supergranule is derived from existing HMI pipeline time-distance data products and compared to the best helioseismic flow model detailed in Duvall and Hanasoge (2013). We find that significant differences exist between them. Unlike the shallow structure predicted by the model, the average HMI supergranule is very extended in depth, exhibiting horizontal outflow down to $7$--$10$~Mm, followed by a weak inflow reaching a depth of $\sim20$~Mm below the photosphere. The maximal velocities in the horizontal direction for the average supergranule are much smaller than the model, and its near-surface flow field RMS value is about an order of magnitude smaller than the often-quoted values of $\sim250-350$~$\rm{m\,s^{-1}}$ for supergranulation. Much of the overall HMI supergranule structure and its weak flow amplitudes can be explained by examining the HMI pipeline averaging kernels for the near-surface inversions, which are found to be very broad in depth, and nearly identical to one another in terms of sensitivity along the $z$-direction. We also show that forward-modeled travel times in the Born approximation using the model (derived from a ray theory approach) are inconsistent with measured travel times for an average supergranule at any distance. Our findings suggest systematic inaccuracies in the typical techniques used to study supergranulation, confirming some of the results in Duvall and Hanasoge (2013).

The light curves of 2768 stars with effective temperatures and surface gravities placing them near the gamma Doradus/delta Scuti instability region were observed as part of the Kepler Guest Observer program from Cycles 1 through 5. The light curves were analyzed in a uniform manner to search for gamma Doradus, delta Scuti, and hybrid star pulsations. The gamma Doradus, delta Scuti, and hybrid star pulsations extend asteroseismology to stars slightly more massive (1.4 to 2.5 solar masses) than our Sun. We find 207 gamma Doradus, 84 delta Scuti, and 32 hybrid candidate stars. Many of these stars are cooler than the red edge of the gamma Doradus instability strip as determined from ground-based observations made before Kepler. A few of our gamma Doradus candidate stars lie on the hot side of the ground-based gamma Doradus instability strip. The hybrid candidate stars cover the entire region between 6200 K and the blue edge of the ground-based delta Scuti instability strip. None of our candidate stars are hotter than the hot edge of the ground-based delta Scuti instability strip. Our discoveries, coupled with the work of others, shows that Kepler has discovered over 2000 gamma Doradus, delta Scuti, and hybrid star candidates in the 116 square degree Kepler field of view. We found relatively few variable stars fainter than magnitude 15, which may be because they are far enough away to lie between spiral arms in our Galaxy, where there would be fewer stars.

In our 2013 Astronomical Review article, we discussed the statistics of variability for 633 faint spectral type A-F stars observed by the Kepler spacecraft during Quarters 6-13. We found six stars that showed no variability with amplitude 20 ppm or greater in the range 0.2 to 24.4 cycles/day, but whose positions in the log g--Teff diagram place them in the delta Sct or gamma Dor pulsation instability regions established from pre-Kepler ground-based observations. Here we present results for 2137 additional stars observed during Quarters 14-17, and find 34 stars that lie within the instability regions. In Paper I, we included a +229 K offset to the Kepler Input Catalog Teff to take into account an average systematic difference between the KIC values and the Teff derived from SDSS color photometry for main-sequence F stars (Pinsonneault et al. 2012). Here we compare the KIC Teff value and the Teff derived from spectroscopy taken by the LAMOST instrument (Molenda-Zakowicz et al. 2013, 2014) for 54 stars common to both samples. We find no trend to support applying the offset, but instead find that a small average temperature decrease relative to the KIC Teff may be more appropriate for the stars in our spectral-type range. If the offset is omitted, only 17 of our 34 `constant' stars fall within the instability regions. For the two `constant' stars also observed by LAMOST, the LAMOST Teff values are cooler than the KIC Teff by several hundred K, and would move these stars out of the instability regions. It is possible that a more accurate determination of their Teff and log g would move some of the other `constant' stars out of the instability regions. However, if average (random) errors in Teff are taken into account, 15 to 52 stars may still persist within the instability regions. Explanations for these `constant' stars, both theoretical and observational, remain to be investigated.

Chromospheric flare ribbons observed in H-alpha appear well-organized when first examined: ribbons impulsively brighten, morphologically evolve, and exponentially decay back to pre-flare levels. Upon closer inspection of H-alpha flares, there is often a significant number of compact areas brightening in concert with the flare eruption but are spatially separated from the evolving flare ribbon. One class of these brightenings is known as sequential chromospheric brightenings (SCBs). SCBs are often observed in the intimidate vicinity of erupting flares and are associated with coronal mass ejections. In the past decade there have been several previous investigations of SCBs. These studies have exclusively relied upon H-alphaimages to discover and analyze these ephemeral brightenings. This work employs the automated detection algorithm of Kirk et al. (2011) to extract the physical qualities of SCBs in observations of ground-based H-alpha images and complementary AIA images in HeII, Civ, and 1700 \AA. The meta-data produced in this tracking process are then culled using complementary Doppler velocities to isolate three distinguishable types of SCBs. From a statistical analysis, we find that the SCBs at the chromospheric H-alpha layer appear earlier, and last longer than their corresponding signatures measured in AIA. From this multi-layer analysis, we infer that SCBs are spatially constrained to the mid-chromosphere. We also derive an energy budget to explain SCBs in which SCBs have a postulated energy of not more than 0.01% of the total flare energy.

Seismology applied to giant planets could drastically change our understanding of their deep interiors, as it has happened with the Earth, the Sun, and many main-sequence and evolved stars. The study of giant planets' composition is important for understanding both the mechanisms enabling their formation and the origins of planetary systems, in particular our own. Unfortunately, its determination is complicated by the fact that their interior is thought not to be homogeneous, so that spectroscopic determinations of atmospheric abundances are probably not representative of the planet as a whole. Instead, the determination of their composition and structure must rely on indirect measurements and interior models. Giant planets are mostly fluid and convective, which makes their seismology much closer to that of solar-like stars than that of terrestrial planets. Hence, helioseismology techniques naturally transfer to giant planets. In addition, two alternative methods can be used: photometry of the solar light reflected by planetary atmospheres, and ring seismology in the specific case of Saturn. The current decade has been promising thanks to the detection of Jupiter's acoustic oscillations with the ground-based imaging-spectrometer SYMPA and indirect detection of Saturn's f-modes in its rings by the NASA Cassini orbiter. This has motivated new projects of ground-based and space-borne instruments that are under development. The purpose of this paper is to support the chapter Seismology of Giant Planets which will appear in May 2015 in the book Extraterrestrial Seismology, published in Cambridge University Press. Here, we review the key questions about jovian interiors as well as the most recent observational results, and refer the reader to the book chapter above-mentioned for further details.

Linear time-distance helioseismic inversions are carried out using several filtering schemes to determine vector flow velocities within two $\sim100^2\,{\rm Mm^2}\times 20\,{\rm Mm}$ realistic magnetohydrodynamic sunspot simulations of 25~hr. One simulation domain contains a model of a full sunspot (i.e. one with both an umbra and penumbra), while the other contains a pore (i.e. a spot without a penumbra). The goal is to test current helioseismic methods using these state-of-the-art simulations of magnetic structures. We find that horizontal flow correlations between inversion and simulation flow maps are reasonably high ($\sim0.5$--0.8) in the upper 3~Mm at distances exceeding 25--30~Mm from spot center, but are substantially lower at smaller distances and larger depths. Inversions of forward-modeled travel times consistently outperform those of our measured travel times in terms of horizontal flow correlations, suggesting that our inability to recover flow structure near these active regions is largely due to the fact that we are unable to accurately measure travel times near strong magnetic features. In many cases the velocity amplitudes from the inversions underestimate those of the simulations by up to 50\%, possibly indicating nonlinearity of the forward problem. In every case, we find that our inversions are unable to recover the vertical flow structure of the simulations at any depth.

Linear time-distance helioseismic inversions are carried out for vector flow velocities using travel times measured from two $\sim 100^2\,{\rm Mm^2}\times 20\,{\rm Mm}$ realistic magnetohydrodynamic quiet-Sun simulations of about 20 hr. The goal is to test current seismic methods on these state-of-the-art simulations. Using recent three-dimensional inversion schemes, we find that inverted horizontal flow maps correlate well with the simulations in the upper $\sim 3$ Mm of the domains for several filtering schemes, including phase-speed, ridge, and combined phase-speed and ridge measurements. In several cases, however, the velocity amplitudes from the inversions severely underestimate those of the simulations, possibly indicating nonlinearity of the forward problem. We also find that, while near-surface inversions of the vertical velocites are best using phase-speed filters, in almost all other example cases these flows are irretrievable due to noise, suggesting a need for statistical averaging to obtain better inferences.

As part of the NASA Kepler Guest Observer program, we requested and obtained long-cadence data on about 2700 faint (magnitude 14-16) Kepler stars with effective temperatures and surface gravities that lie near or within the pulsation instability region for main-sequence gamma Doradus and delta Scuti pulsating variables. These variables are of spectral type A-F with masses of 1.4 to 2.5 solar masses. The delta Scuti stars pulsate in radial and non-radial acoustic modes, with periods of a few hours (frequencies around 10 cycles/day), while gamma Doradus variables pulsate in nonradial gravity modes with periods 0.3 to 3 days (frequencies around 1 cycle/day). Here we consider the light curves and Fourier transforms of 633 stars in an unbiased sample observed by Kepler in Quarters 6-13 (June 2010-June 2012). We show the location of these stars in the log surface gravity--effective temperature diagram compared to the instability region limits established from ground-based observations, and taking into account uncertainties and biases in the Kepler Input Catalog T_eff values. While hundreds of variables have been discovered in the Kepler data, about 60% of the stars in our sample do not show any frequencies between 0.2 and 24.4 cycles per day with amplitude above 20 parts per million. We find that six of these apparently constant stars lie within the pulsation instability region. We discuss some possible reasons that these stars do not show photometric variability in the Kepler data. We also comment on the non-constant stars, and on 26 variable-star candidates, many of which also do not lie within the expected instability regions.

Large-scale plasma flows in the Sun's convection zone likely play a major role in solar dynamics on decadal timescales. In particular, quantifying meridional motions is a critical ingredient for understanding the solar cycle and the transport of magnetic flux. Because the signal of such features can be quite small in deep solar layers and be buried in systematics or noise, the true meridional velocity profile has remained elusive. We perform time-distance helioseismology measurements on several years worth of GONG Doppler data. A spherical harmonic decomposition technique is applied to a subset of acoustic modes to measure travel-time differences to try to obtain signatures of meridional flows throughout the solar convection zone. Center-to-limb systematics are taken into account in an intuitive, yet ad hoc manner. Travel-time differences near the surface that are consistent with a poleward flow in each hemisphere and are similar to previous work are measured. Additionally, measurements in deep layers near the base of the convection zone suggest a possible equatorward flow, as well as partial evidence of a sign change in the travel-time differences at mid-convection zone depths. This analysis on an independent data set using different measurement techniques strengthens recent conclusions that the convection zone may have multiple "cells" of meridional flow. The results may challenge the common understanding of one large conveyor belt operating in the solar convection zone. Further work with helioseismic inversions and a careful study of systematic effects are needed before firm conclusions of these large-scale flow structures can be made.

Among 19 red-giant stars belonging to eclipsing binary systems that have been identified in Kepler data, 15 display solar-like oscillations. We study whether the absence of mode detection in the remaining 4 is an observational bias or possibly evidence of mode damping that originates from tidal interactions. A careful analysis of the corresponding Kepler light curves shows that modes with amplitudes that are usually observed in red giants would have been detected if they were present. We observe that mode depletion is strongly associated with short-period systems, in which stellar radii account for 16-24 % of the semi-major axis, and where red-giant surface activity is detected. We suggest that when the rotational and orbital periods synchronize in close binaries, the red-giant component is spun up, so that a dynamo mechanism starts and generates a magnetic field, leading to observable stellar activity. Pressure modes would then be damped as acoustic waves dissipate in these fields.

Context: OB stars are important in the chemistry and evolution of the Universe, but the sample of targets well understood from an asteroseismological point of view is still too limited to provide feedback on the current evolutionary models. Our study extends this sample with two spectroscopic binary systems. AIMS. Our goal is to provide orbital solutions, fundamental parameters and abundances from disentangled high-resolution high signal-to-noise spectra, as well as to analyse and interpret the variations in the Kepler light curve of these carefully selected targets. This way we continue our efforts to map the instability strips of beta Cep and SPB stars using the combination of high-resolution ground-based spectroscopy and uninterrupted space-based photometry. Methods: We fit Keplerian orbits to radial velocities measured from selected absorption lines of high-resolution spectroscopy using synthetic composite spectra to obtain orbital solutions. We use revised masks to obtain optimal light curves from the original pixel-data from the Kepler satellite, which provided better long term stability compared to the pipeline processed light curves. We use various time-series analysis tools to explore and describe the nature of variations present in the light curve. Results: We find two eccentric double-lined spectroscopic binary systems containing a total of three main sequence B-type stars (and one F-type component) of which at least one in each system exhibits light variations. The light curve analysis (combined with spectroscopy) of the system of two B stars points towards the presence of tidally excited g modes in the primary component. We interpret the variations seen in the second system as classical g mode pulsations driven by the kappa mechanism in the B type primary, and explain the unexpected power in the p mode region as a result of nonlinear resonant mode excitation.

Red-giant stars are an incredible source of information for testing models of stellar evolution, as asteroseismology has opened up a window into their interiors. Such insights are a direct result of the unprecedented data from space missions Kepler and CoRoT as well as recent theoretical advances. Eclipsing binaries are also fundamental astrophysical objects, and when coupled with asteroseismology, binaries would provide two independent methods to obtain masses and radii and exciting opportunities to develop highly constrained stellar models. The possibility of discovering pulsating red giants in eclipsing binary systems is therefore an important goal that could potentially offer very robust characterization of these systems. Hitherto, only 1 case has been discovered with Kepler. We cross-correlated the detected red-giant and eclipsing-binary catalogs from Kepler data to find candidate systems. Light-curve modeling and mean asteroseismic properties are combined to yield measurements of periods, masses, radii, temperatures, eclipse timing variations, core rotation rates, and red-giant evolutionary state. After elimination of false positives, out of the 70 systems common to the red-giant and eclipsing-binary catalogs we find 13 new strong candidates that are pulsating red giants in eclipsing binaries, and 10 more that are pulsating red giants in hierarchical triple systems. The systems span a range of orbital eccentricities, periods, and spectral types for the red-giant companion. One case even suggests a system composed of two red-giants. The discovery of pulsating red giants in multiple systems provides exciting test beds for precise astrophysical modeling, and follow-up spectroscopic observations of the candidate systems are encouraged.

Chromospheric wave activity around flares and filaments has been a research focus for years, and could provide indirect measurements of local conditions that are not otherwise accessible. One interesting observed phenomenon is oscillations in filaments, activated by distant flares and the large-scale waves they produce. Characteristics of these oscillations, such as periods, amplitudes, and lifetimes, can provide unique information about the filament. We measure oscillation properties in flares and filaments from H-alpha chromospheric data using a new method that provides important spatial and frequency content of the dynamics. We apply the method to two flare events where filaments are observed to oscillate and determine their properties. We find strong oscillatory signal in flaring active regions in the chromosphere over a range of frequencies. Two filaments are found to oscillate without any detectable chromospheric wave acting as an activation mechanism. We find that filaments oscillate with periods of tens of minutes, but variations are significant at small spatial scales along the filamentary region. The results suggest there is a frequency dependence of the oscillation amplitude, as well as a spatial dependence along single filaments that is more difficult to quantify. It also appears that the strength of the oscillations do not necessarily depend on the strength of the trigger, although there are other possible effects that make this conclusion preliminary. Applications of this technique to other events and different data sets will provide important new insights into the local energy densities and magnetic fields associated with dynamic chromospheric structures.

Jupiter is expected to pulsate in a spectrum of acoustic modes and recent re-analysis of a spectroscopic time series has identified a regular pattern in the spacing of the frequencies \citep{gaulme2011}. This exciting result can provide constraints on gross Jovian properties and warrants a more in-depth theoretical study of the seismic structure of Jupiter. With current instrumentation, such as the SYMPA instrument \citep{schmider2007} used for the \citet{gaulme2011} analysis, we assume that, at minimum, a set of global frequencies extending up to angular degree $\ell=25$ could be observed. In order to identify which modes would best constrain models of Jupiter's interior and thus help motivate the next generation of observations, we explore the sensitivity of derived parameters to this mode set. Three different models of the Jovian interior are computed and the theoretical pulsation spectrum from these models for $\ell\leq 25$ is obtained. We compute sensitivity kernels and perform linear inversions to infer details of the expected discontinuities in the profiles in the Jovian interior. We find that the amplitude of the sound-speed jump of a few percent in the inner/outer envelope boundary seen in two of the applied models should be reasonably inferred with these particular modes. Near the core boundary where models predict large density discontinuities, the location of such features can be accurately measured, while their amplitudes have more uncertainty. These results suggest that this mode set would be sufficient to infer the radial location and strength of expected discontinuities in Jupiter's interior, and place strong constraints on the core size and mass. We encourage new observations to detect these Jovian oscillations.

We make a comparison between small scale chromospheric brightenings and energy release processes through examining the temporal evolution of sequential chromospheric brightenings (SCBs), derive propagation velocities, and propose a connection of the small-scale features to solar flares. Our automated routine detects and distinguishes three separate types of brightening regularly observed in the chromosphere: plage, flare ribbon, and point brightenings. By studying their distinct dynamics, we separate out the flare-associated bright points commonly known as SCBs and identify a propagating Moreton wave. Superimposing our detections on complementary off-band images, we extract a Doppler velocity measurement beneath the point brightening locations. Using these dynamic measurements, we put forward a connection between point brightenings, the erupting flare, and overarching magnetic loops. A destabilization of the pre-flare loop topology by the erupting flare directly leads to the SCBs observed.

We report on the physical properties of solar sequential chromospheric brightenings (SCBs) observed in conjunction with moderate-sized chromospheric flares with associated CMEs. To characterize these ephemeral events, we developed automated procedures to identify and track subsections (kernels) of solar flares and associated SCBs using high resolution H-alpha images. Following the algorithmic identification and a statistical analysis, we compare and find the following: SCBs are distinctly different from flare kernels in their temporal characteristics of intensity, Doppler structure, duration, and location properties. We demonstrate that flare ribbons are themselves made up of subsections exhibiting differing characteristics. Flare kernels are measured to have a mean propagation speed of 0.2 km/s and a maximum speed of 2.3 km/s over a mean distance of 5 x 10^3 km. Within the studied population of SCBs, different classes of characteristics are observed with coincident negative, positive, or both negative and positive Doppler shifts of a few km/s. The appearance of SCBs precede peak flare intensity by ~12 minutes and decay ~1 hour later. They are also found to propagate laterally away from flare center in clusters at 41 km/s or 89 km/s. Given SCBs distinctive nature compared to flares, we suggest a different physical mechanism relating to their origin than the associated flare. We present a heuristic model of the origin of SCBs.

The light curves of 252 B-star candidates in the Kepler data base are analyzed in a similar fashion to that done by Balona et al. (2011) to further characterize B star variability, increase the sample of variable B stars for future study, and to identify stars whose power spectra include particularly interesting features such as frequency groupings. Stars are classified as either constant light emitters, $\beta$ Cep stars, slowly pulsating B stars, hybrid pulsators, binaries or stars whose light curves are dominated by rotation (Bin/Rot), hot subdwarfs, or white dwarfs. One-hundred stars in our sample were found to be either light contants or to be variable at a level of less than 0.02 mmag. We increase the number of candidate B-star variables found in the Kepler data base by Balona et al. (2011) in the following fashion: $\beta$ Cep stars from 0 to 10, slowly pulsating B stars from 8 to 54, hybrid pulsators from 7 to 21, and Bin/Rot stars from 23 to 82. For comparison purposes, approximately 51 SPBs and 6 hybrids had been known prior to 2007. The number of $\beta$ Cep stars known prior to 2004 was 93. A secondary result of this study is the identification of an additional 11 pulsating white dwarf candidates, four of which possess frequency groupings.

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

The Kepler spacecraft is providing time series of photometric data with micromagnitude precision for hundreds of A-F type stars. We present a first general characterization of the pulsational behaviour of A-F type stars as observed in the Kepler light curves of a sample of 750 candidate A-F type stars. We propose three main groups to describe the observed variety in pulsating A-F type stars: gamma Dor, delta Sct, and hybrid stars. We assign 63% of our sample to one of the three groups, and identify the remaining part as rotationally modulated/active stars, binaries, stars of different spectral type, or stars that show no clear periodic variability. 23% of the stars (171 stars) are hybrid stars, which is a much larger fraction than what has been observed before. We characterize for the first time a large number of A-F type stars (475 stars) in terms of number of detected frequencies, frequency range, and typical pulsation amplitudes. The majority of hybrid stars show frequencies with all kinds of periodicities within the gamma Dor and delta Sct range, also between 5 and 10 c/d, which is a challenge for the current models. We find indications for the existence of delta Sct and gamma Dor stars beyond the edges of the current observational instability strips. The hybrid stars occupy the entire region within the delta Sct and gamma Dor instability strips, and beyond. Non-variable stars seem to exist within the instability strips. The location of gamma Dor and delta Sct classes in the (Teff,logg)-diagram has been extended. We investigate two newly constructed variables 'efficiency' and 'energy' as a means to explore the relation between gamma Dor and delta Sct stars. Our results suggest a revision of the current observational instability strips, and imply an investigation of other pulsation mechanisms to supplement the kappa mechanism and convective blocking effect to drive hybrid pulsations.

We present a new automated algorithm to identify, track, and characterize small-scale brightening associated with solar eruptive phenomena observed in H{\alpha}. The temporal spatially-localized changes in chromospheric intensities can be separated into two categories: flare ribbons and sequential chromospheric brightenings (SCBs). Within each category of brightening we determine the smallest resolvable locus of pixels, a kernel, and track the temporal evolution of the position and intensity of each kernel. This tracking is accomplished by isolating the eruptive features, identifying kernels, and linking detections between frames into trajectories of kernels. We fully characterize the evolving intensity and morphology of the flare ribbons by observing the tracked flare kernels in aggregate. With the location of SCB and flare kernels identified, they can easily be overlaid on top of complementary data sets to extract Doppler velocities and magnetic field intensities underlying the kernels. This algorithm is adaptable to any dataset to identify and track solar features.

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

We present the ground-based activities within the different working groups of the Kepler Asteroseismic Science Consortium (KASC). The activities aim at the systematic characterization of the 5000+ KASC targets, and at the collection of ground-based follow-up time-series data of selected promising Kepler pulsators. So far, 36 different instruments at 31 telescopes on 23 different observatories in 12 countries are in use, and a total of more than 530 observing nights has been awarded. (Based on observations made with the Isaac Newton Telescope, William Herschel Telescope, Nordic Optical Telescope, Telescopio Nazionale Galileo, Mercator Telescope (La Palma, Spain), and IAC-80 (Tenerife, Spain). Also based on observations taken at the observatories of Sierra Nevada, San Pedro Martir, Vienna, Xinglong, Apache Point, Lulin, Tautenburg, Loiano, Serra la Nave, Asiago, McDonald, Skinakas, Pic du Midi, Mauna Kea, Steward Observatory, Mt Wilson, Bialkow Observatory of the Wroclaw University, Piszkesteto Mountain Station, Observatoire de Haute Provence, and Centro Astronomico Hispano Aleman at Calar Alto. Based on data from the AAVSO International Database.)

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

Observations of the pulsations of stars can be used to infer their interior structure and test theoretical models. The main sequence $\gamma$ Doradus (Dor) and $\delta$ Scuti (Sct) stars with masses 1.2-2.5 $M_{\sun}$ are particularly useful for these studies. The $\gamma$ Dor stars pulsate in high-order $g$ modes with periods of order 1 day, driven by convective blocking at the base of their envelope convection zone. The $\delta$ Sct stars pulsate in low-order $g$ and $p$ modes with periods of order 2 hours, driven by the $\kappa$ mechanism operating in the Heii ionization zone. Theory predicts an overlap region in the Hertzsprung-Russell diagram between instability regions, where 'hybrid' stars pulsating in both types of modes should exist. The two types of modes with properties governed by different portions of the stellar interior provide complementary model constraints. Among the known $\gamma$ Dor and $\delta$ Sct stars, only four have been confirmed as hybrids. Now, analysis of combined Quarter 0 and Quarter 1 Kepler data for hundreds of variable stars shows that the frequency spectra are so rich that there are practically no pure $\delta$ Sct or $\gamma$ Dor pulsators, i.e. essentially all of the stars show frequencies in both the $\delta$ Sct and $\gamma$ Dor frequency range. A new observational classification scheme is proposed that takes into account the amplitude as well as the frequency, and is applied to categorize 234 stars as $\delta$ Sct, $\gamma$ Dor, $\delta$ Sct/$\gamma$ Dor or $\gamma$ Dor/$\delta$ Sct hybrids.

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

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

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

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

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