Using a nonlinear mean-field solar dynamo model, we study relationships between the amplitude of the `extended' mode of migrating zonal flows (`torsional oscillations') and magnetic cycles, and investigate whether properties the torsional oscillations in subsurface layers and in the deep convection zone can provide information about the future solar cycles. We consider two types of dynamo models: models with regular variations of the alpha-effect, and models with stochastic fluctuations, simulating `long'- and 'short-memory' types of magnetic activity variations. It is found that torsional oscillation parameters, such the zonal acceleration, show a considerable correlation with the magnitude of the subsequent cycles with a time lag of 11-20 yr. The sign of the correlation and the time-lag parameters can depend on the depth and latitude of the torsional oscillations as well as on the properties of long-term (`centennial') variations of the dynamo cycles. The strongest correlations are found for the zonal acceleration at high latitudes at the base of the convection zone. The model results demonstrate that helioseismic observations of the torsional oscillations can be useful for advanced prediction of the solar cycles, one-two sunspot cycles ahead.

The authors of this white paper met on 16-17 January 2020 at the New Jersey Institute of Technology, Newark, NJ, for a 2-day workshop that brought together a group of heliophysicists, data providers, expert modelers, and computer/data scientists. Their objective was to discuss critical developments and prospects of the application of machine and/or deep learning techniques for data analysis, modeling and forecasting in Heliophysics, and to shape a strategy for further developments in the field. The workshop combined a set of plenary sessions featuring invited introductory talks interleaved with a set of open discussion sessions. The outcome of the discussion is encapsulated in this white paper that also features a top-level list of recommendations agreed by participants.

Sunquakes are one of the more distinct secondary phenomena related to solar flares, where energy deposition in the lower layers of the Sun's atmosphere excites acoustic waves easily visible in Helioseismic and Magnetic Imager (HMI) dopplergrams. We explore two possible sources of sunquakes in the context of the electron beam hypothesis: an instantaneous transfer of momentum and a gradual applied force due to flare eruption. We model the sunquake excitation and compare with five observed sunquake events using a cross-correlation analysis. We find that at least half the events studied are consistent with the electron beam hypothesis and estimate the energy required to excite the sunquakes to be within the range determined by previous studies.

Identification of solar coronal holes (CHs) provides information both for operational space weather forecasting and long-term investigation of solar activity. Source data for the first problem are typically most recent solar disk observations, while for the second problem it is convenient to consider solar synoptic maps. Motivated by the idea that the concept of CHs should be similar for both cases we investigate universal models that can learn a CHs segmentation in disk images and reproduce the same segmentation in synoptic maps. We demonstrate that Convolutional Neural Networks (CNN) trained on daily disk images provide an accurate CHs segmentation in synoptic maps and their pole-centric projections. Using this approach we construct a catalog of synoptic maps for the period of 2010-20 based on SDO/AIA observations in the 193 Angstrom wavelength. The obtained CHs synoptic maps are compared with magnetic synoptic maps in the time-latitude and time-longitude diagrams. The initial results demonstrate that while in some cases the CHs are associated with magnetic flux transport events there are other mechanisms contributing to the CHs formation and evolution. To stimulate further investigations the catalog of synoptic maps is published in open access.

The paper presents results of a search for helioseismic events (sunquakes) produced by M-X class solar flares during Solar Cycle 24. The search is performed by analyzing photospheric Dopplergrams from Helioseismic Magnetic Imager (HMI). Among the total number of 500 M-X class flares, 94 helioseismic events were detected. Our analysis has shown that many strong sunquakes were produced by solar flares of low M class (M1-M5), while in some powerful X-class flares helioseismic waves were not observed or were weak. Our study also revealed that only several active regions were characterized by the most efficient generation of helioseismic waves during flares. We found that the sunquake power correlates with the maximum value of the soft X-ray flux time derivative better than with the X-ray class, indicating that the sunquake mechanism is associated with high-energy particles. We also show that the seismically active flares are more impulsive than the flares without helioseismic perturbations. We present a new catalog of helioseismic solar flares, which opens opportunities for performing statistical studies to better understand the physics of sunquakes as well as the flare energy release and transport.

We perform an analysis of the line-of-sight (LOS) observables of the SDO/HMI for models of the solar atmosphere heated by precipitating high-energy electrons during solar flares. The radiative hydrodynamic (RADYN) flare models are obtained from the F-CHROMA database. The Stokes profiles for the Fe 6173 A line observed by SDO/HMI are calculated using the radiative transfer code RH1.5D assuming statistical equilibrium for atomic level populations and imposing uniform background vertical magnetic field of various strength. The SDO/HMI observing sequence and LOS data processing pipeline algorithm are applied to derive the observables (continuum intensity, line depth, Doppler velocity, LOS magnetic field). Our results reveal that the strongest deviations of the observables from the actual spectroscopic line parameters are found for the model with a total energy deposited of $E_{total}=1.0\times{}10^{12}$ erg cm$^{-2}$, injected with power-law spectral index of $\delta=3$ above a low-energy cutoff of $E_{c}=25$ keV. The magnitudes of the velocity and magnetic field deviations depend on the imposed magnetic field, and can reach 0.35 km/s for LOS velocities, 90 G for LOS magnetic field, and 3\% for continuum enhancement for the 1000 G imposed LOS magnetic field setup. For $E_{total}\geq{}3.0\times{}10^{11}$ erg cm$^{-2}$ models, the velocity and magnetic field deviations are most strongly correlated with the energy flux carried by $\sim$50 keV electrons, and the continuum enhancement is correlated with the synthesized $\sim$55-60 keV hard X-ray photon flux. The relatively low magnitudes of perturbations of the observables and absence of magnetic field sign reversals suggest that the considered radiative hydrodynamic beam heating models augmented with the uniform vertical magnetic field setups cannot explain strong transient changes found in the SDO/HMI observations.

The Kepler Asteroseismic Legacy Project provided frequencies, separation ratios, error estimates, and covariance matrices for 66 Kepler main sequence targets. Most of the previous analysis of these data was focused on fitting standard stellar models. We present results of direct asteroseismic inversions using the method of optimally localized averages (OLA), which effectively eliminates the surface effects and attempts to resolve the stellar core structure. The inversions are presented for various structure properties, including the density stratification and sound speed. The results show that the mixed modes observed in post-main sequence F-type stars allow us to resolve the stellar core structure and reveal significant deviations from the evolutionary models obtained by the grid-fitting procedure to match the observed oscillation frequencies.

Various models of solar subsurface stratification are tested in the global EULAG-MHD solver to simulate diverse regimes of near-surface convective transport. Sub- and superadiabacity are altered at the surface of the model ($ r > 0.95~R_{\odot}$) to either suppress or enhance convective flow speeds in an effort to investigate the impact of the near-surface layer on global dynamics. A major consequence of increasing surface convection rates appears to be a significant alteration of the distribution of angular momentum, especially below the tachocline where the rotational frequency predominantly increases at higher latitudes. These hydrodynamic changes correspond to large shifts in the development of the current helicity in this stable layer ($r<0.72R_{\odot}$), significantly altering its impact on the generation of poloidal and toroidal fields at the tachocline and below, acting as a major contributor towards transitions in the dynamo cycle. The enhanced near-surface flow speed manifests in a global shift of the toroidal field ($B_{\phi}$) in the butterfly diagram - from a North-South symmetric pattern to a staggered anti-symmetric emergence.

We present a nonlinear mean-field model of the solar interior dynamics and dynamo, which reproduces the observed cyclic variations of the global magnetic field of the Sun, as well as the differential rotation and meridional circulation. Using this model, we explain, for the first time, the extended 22-year pattern of the solar torsional oscillations, observed as propagation of zonal variations of the angular velocity from high latitudes to the equator during the time equal to the full dynamo cycle. In the literature, this effect is usually attributed to the so-called "extended solar cycle". In agreement with the commonly accepted idea our model shows that the torsional oscillations can be driven by a combinations of magnetic field effects acting on turbulent angular momentum transport, and the large-scale Lorentz force. We find that the 22-year pattern of the torsional oscillations can result from a combined effect of an overlap of subsequent magnetic cycles and magnetic quenching of the convective heat transport. The latter effect results in cyclic variations of the meridional circulation in the sunspot formation zone, in agreement with helioseismology results. The variations of the meridional circulation together with other drivers of the torsional oscillations maintain their migration to the equator during the 22-year magnetic cycle, resulting in the observed extended pattern of the torsional oscillations.

Various models of solar sub-surface stratification are tested in the global EULAG-MHD solver to simulate diverse regimes of near-surface convective transport. Sub- and superadiabacity are altered at the surface of the model ($ r > 0.95~R_{\odot}$) to either suppress or enhance convective flow speeds in an effort to investigate the impact of the near-surface layer on global dynamics. A major consequence of increasing surface convection rates appears to be a significant alteration of the distribution of angular momentum, especially below the tachocline where the rotational frequency predominantly increases at higher latitudes. These hydrodynamic changes correspond to large shifts in the development of the current helicity in this stable layer ($r<0.72R_{\odot}$), significantly altering its impact on the generation of poloidal and toroidal fields at the tachocline and below, acting as a major contributor towards transitions in the dynamo cycle. The enhanced near-surface flow speed manifests in a global shift of the toroidal field ($B_{\phi}$) in the butterfly diagram -- from a North-South symmetric pattern to a staggered anti-symmetric emergence.

Two fundamental properties of stellar magnetic fields have been determined by observations for solar-like stars with different Rossby numbers (Ro), namely, the magnetic field strength and the magnetic cycle period. The field strength exhibits two regimes: 1) for fast rotation it is independent of Ro, 2) for slow rotation it decays with Ro following a power law. For the magnetic cycle period two regimes of activity, the active and inactive branches, also have been identified. For both of them, the longer the rotation period, the longer the activity cycle. Using global dynamo simulations of solar like stars with Rossby numbers between ~0.4 and ~2, this paper explores the relevance of rotational shear layers in determining these observational properties. Our results, consistent with non-linear alpha^2-Omega dynamos, show that the total magnetic field strength is independent of the rotation period. Yet at surface levels, the origin of the magnetic field is determined by Ro. While for Ro<1 it is generated in the convection zone, for Ro>1 strong toroidal fields are generated at the tachocline and rapidly emerge towards the surface. In agreement with the observations, the magnetic cycle period increases with the rotational period. However, a bifurcation is observed for Ro~1, separating a regime where oscillatory dynamos operate mainly in the convection zone, from the regime where the tachocline has a predominant role. In the latter the cycles are believed to result from the periodic energy exchange between the dynamo and the magneto-shear instabilities developing in the tachocline and the radiative interior.

We propose a novel approach to reconstruct the surface magnetic helicity density on the Sun or sun-like stars. The magnetic vector potential is determined via decomposition of vector magnetic field measurements into toroidal and poloidal components. The method is verified using data from a non-axisymmetric dynamo model. We apply the method to vector field synoptic maps from Helioseismic and Magnetic Imager (HMI) onboard of Solar Dynamics Observatory (SDO) to study evolution of the magnetic helicity density during solar cycle 24. It is found that the mean helicity density of the non-axisymmetric magnetic field of the Sun evolves in a way which is similar to that reported for the current helicity density of the solar active regions. It has predominantly the negative sign in the northern hemisphere, and it is positive in the southern hemisphere. Also, the hemispheric helicity rule for the non-axisymmetric magnetic field showed the sign inversion at the end of cycle 24. Evolution of magnetic helicity density of large-scale axisymmetric magnetic field is different from that expected in dynamo theory. On one hand, the mean large- and small-scale components of magnetic helicity density display the hemispheric helicity rule of opposite sign at the beginning of cycle 24. However, later in the cycle, the two helicities exhibit the same sign in contrast with the theoretical expectations.

We present a statistical analysis of properties of Soft X-Ray (SXR) emission, plasma temperature (T), and emission measure (EM), derived from GOES observations of flares in 2002-2017. The temperature and emission measures are obtained using the TEBBS algorithm (Ryan et al. 2012), which delivers reliable results together with uncertainties even for weak B-class flare events. More than 96% of flares demonstrate a sequential appearance of T, SXR, and EM maxima, in agreement with the expected behavior of the chromospheric evaporation process. The relative number of such flares increases with increasing the SXR flux maximum. The SXR maximum is closer in time to the T maximum for B-class flares than for >=C-class flares, while it is very close to the EM maximum for M- and X-class flares. We define flares as "T-controlled" if the time interval between the SXR and T maxima is at least two times shorter than the interval between the EM and SXR maxima, and as "EM-controlled" if the time interval between the EM and SXR maxima is at least two times shorter than the interval between the SXR and T maxima. For any considered flare class range, the T-controlled events compared to EM-controlled events have: a) higher EM but lower T; b) longer durations and shorter relative growth times; c) longer FWHM and characteristic decay times. Interpretation of these statistical results based on analysis of a single loop dynamics suggests that for flares of the same class range, the T-controlled events can be developed in longer loops than the EM-controlled events.

Torsional oscillations represent bands of fast and slow zonal flows around the Sun, which extend deep into the convection zone and migrate during solar cycles towards the equator following the sunspot butterfly diagram. Analysis of helioseismology data obtained in 1996-2018 for almost two solar cycles reveals zones of deceleration of the torsional oscillations inside the Sun due to dynamo-generated magnetic field. The zonal deceleration originates near the bottom of the convection zone at high latitudes, and migrates to the surface revealing patterns of magnetic dynamo waves predicted by the Parker's dynamo theory. The analysis reveals that the primary seat of the solar dynamo is located in a high-latitude zone of the tachocline. It suggests a dynamo scenario that can explain 'extended solar cycles' previously observed in the evolving shape of the solar corona. The results show a substantial decrease of the zonal acceleration in the current solar cycle and indicate further decline of activity in the next solar cycle. Although the relationship between the magnitude of zonal deceleration and the amount of emerged toroidal field that leads to formation of sunspots is not yet established, the results open a new perspective for solar cycle modeling and prediction using helioseismology data.

We present a statistical study of chromospheric evaporation in solar flares using simultaneous observations by the RHESSI X-ray telescope and the IRIS UV spectrograph. The results are compared with radiation hydrodynamic flare models from the F-CHROMA RADYN database. For each event, we study spatially-resolved Doppler shifts of spectral lines formed in the transition region (C\,II\,1334.5\,{\AA}) and hot coronal plasma (Fe\,XXI\,1354.1\,{\AA}) to investigate the dynamics of the solar atmosphere during the flare impulsive phase. We estimate the energy fluxes deposited by high-energy electrons using X-ray imaging spectroscopy and assuming the standard thick-target model. Using the RADYN flare models, the RH 1.5D radiative transfer code and the Chianti atomic line database, we calculate C\,II and Fe\,XXI line profiles and compare with the observations. While the RADYN models predict a correlation between the Doppler shifts and deposited energy flux for both lines, this was only observed in the C II data. Several quantitative discrepancies are found between the observations and models: the Fe\,XXI Doppler shifts are substantially stronger in the models than in the data; the C\,II mean blueshifts are absent in the observations but predicted by the models. The transition energy between `gentle' and `explosive' evaporation regimes estimated from the observations ($2-8\times{}10^{9}$~erg~cm$^{-2}$~s$^{-1}$) and derived from the models ($2.2-10.1\times{}10^{9}$~erg~cm$^{-2}$~s$^{-1}$) are comparable with each other. The results illustrate relationships among the processes of chromospheric evaporation, response of the colder layers, and the flare energy flux deposited by high-energy electrons, although demonstrating discrepancy between analyzed observations and RADYN models.

The authors of this report met on 28-30 March 2018 at the New Jersey Institute of Technology, Newark, New Jersey, for a 3-day workshop that brought together a group of data providers, expert modelers, and computer and data scientists, in the solar discipline. Their objective was to identify challenges in the path towards building an effective framework to achieve transformative advances in the understanding and forecasting of the Sun-Earth system from the upper convection zone of the Sun to the Earth's magnetosphere. The workshop aimed to develop a research roadmap that targets the scientific challenge of coupling observations and modeling with emerging data-science research to extract knowledge from the large volumes of data (observed and simulated) while stimulating computer science with new research applications. The desire among the attendees was to promote future trans-disciplinary collaborations and identify areas of convergence across disciplines. The workshop combined a set of plenary sessions featuring invited introductory talks and workshop progress reports, interleaved with a set of breakout sessions focused on specific topics of interest. Each breakout group generated short documents, listing the challenges identified during their discussions in addition to possible ways of attacking them collectively. These documents were combined into this report-wherein a list of prioritized activities have been collated, shared and endorsed.

It is generally accepted that solar acoustic (p) modes are excited by near-surface turbulent motions, in particular, by downdrafts and interacting vortices in intergranular lanes. Recent analysis of Solar Dynamics Observatory data by (Zhao et al., 2015) revealed fast-moving waves around sunspots, which are consistent with magnetoacoustic waves excited approximately 5 Mm beneath the sunspot. We analyzed 3D radiative MHD simulations of solar magnetoconvection with a self-organized pore-like magnetic structure, and identified more than 600 individual acoustic events both inside and outside of this structure. By performing a case-by-case study, we found that surrounding the magnetic structure, acoustic sources are associated with downdrafts. Their depth correlates with downdraft speed and magnetic fields. The sources often can be transported into deeper layers by downdrafts. The wave front shape, in the case of a strong or inclined downdraft, can be stretched along the downdraft. Inside the magnetic structure, excitation of acoustic waves is driven by converging flows. Frequently, strong converging plasma streams hit the structure boundaries, causing compressions in its interior that excite acoustic waves. Analysis of the depth distribution of acoustic events shows the strongest concentration at 0.2 - 1 Mm beneath the surface for the outside sources and 2.5 - 3 Mm for the excitation event inside the structure.

We explore effects of random non-axisymmetric perturbations of kinetic helicity (the $\alpha$ effect) and diffusive decay of bipolar magnetic regions on generation and evolution of large-scale non-axisymmetric magnetic fields on the Sun. Using a reduced 2D nonlinear mean-field dynamo model and assuming that bipolar regions emerge due to magnetic buoyancy in situ of the large-scale dynamo action, we show that fluctuations of the $\alpha$ effect can maintain the non-axisymmetric magnetic fields through a solar-type $\alpha^{2}\Omega$ dynamo process. It is found that diffusive decay of bipolar active regions is likely to be the primary source of the non-axisymmetric magnetic fields observed on the Sun. Our results show that the non-axisymmetric dynamo model with stochastic perturbations of the $\alpha$ effect can explain periods of extremely high activity (`super-cycle' events) as well as periods of deep decline of magnetic activity. We compare the models with synoptic observations of solar magnetic fields for the last four activity cycles, and discuss implications of our results for interpretation of observations of stellar magnetic activity.

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

We analyze the relationship between the flare X-ray peak flux, and characteristics of the Polarity Inversion Line (PIL) and Active Regions (AR), derived from line-of-sight (LOS) magnetograms. The PIL detection algorithm based on a magnetogram segmentation procedure is applied for each AR with 1 hour cadence. The PIL and AR characteristics are associated with the AR flare history and divided into flaring and non-flaring cases. Effectiveness of the derived characteristics for flare forecasting is determined by the number of non-flaring cases separated from flaring cases by a certain threshold, and by their Fisher ranking score. The Support Vector Machine (SVM) classifier trained only on the PIL characteristics is used for the flare prediction. We have obtained the following results: (1) the PIL characteristics are more effective than global characteristics of ARs, (2) the highest True Skill Statistics (TSS) values of 0.76$\pm$0.03 for $\geq$M1.0 flares and 0.84$\pm$0.07 for $\geq$X1.0 flares are obtained using the "Sigmoid" SVM kernel, (3) the TSS scores obtained using only the LOS magnetograms are slightly lower than the scores obtained using vector magnetograms, but significantly better than current expert-based predictions, (4) for prediction of $\geq$M1.0 class flares 74.4 % of all cases, and 91.2 % for $\geq$X1.0 class, can be pre-classified as negative with no significant effect on the results, (5) the inclusion of global AR characteristics does not improve the forecast. The study confirms the unique role of the PIL region characteristics in the flare initiation process, and demonstrate possibilities of flare forecasting using only the line-of-sight magnetograms.

The questions whether the Sun shrinks with the solar activity and what causes this have been a subject of debate. Helioseismology provides means to measure with high precision the radial displacement of subsurface layers, co-called "seismic radius", through analysis of oscillation frequencies of surface gravity (f) modes. Here, we present results of a new analysis of twenty one years of helioseismology data from two space missions, Solar and Heliospheric Observatory (SoHO) and Solar Dynamics Observatory (SDO), which allow us to resolve previous uncertainties and compare variations of the seismic radius in two solar cycles. After removing the f-mode frequency changes associated with the surface activity we find that the mean seismic radius is reduced by 1-2 km during the solar maxima, and that most significant variations of the solar radius occur beneath the visible surface of the Sun at the depth of about 5 Mm, where the radius is reduced by 5-8 km. These variations can be interpreted as changes in the solar subsurface structure caused by predominately vertical ~10 kG magnetic field.

The measurement of the Sun's diameter has been first tackled by the Greek astronomers from a geometric point of view. Their estimation of ~1800", although incorrect, was not truly called into question for several centuries. The first pioneer works for measuring the Sun's diameter with an astrometric precision were made around the year 1660 by Gabriel Mouton, then by Picard and La Hire. A canonical value of the solar radius of 959".63 was adopted by Auwers in 1891. Despite considerable efforts during the second half of the XXth century, involving dedicated space instruments, no consensus was reached on this issue. However, with the advent of high sensitivity instruments on board satellites, such as the Michelson Doppler Imager (MDI) on Solar and Heliospheric Observatory (SoHO) and the Helioseismic and Magnetic Imager (HMI) aboard NASA's Solar Dynamics Observatory (SDO), it was possible to extract with an unprecedented accuracy the surface gravity oscillation f modes, over nearly two solar cycles, from 1996 to 2017. Their analysis in the range of angular degree l=140-300 shows that the so-called "seismic radius" exhibits a temporal variability in anti-phase with the solar activity. Even if the link between the two radii (photospheric and seismic) can be made only through modeling, such measurements provide an interesting alternative which led to a revision of the standard solar radius by the International Astronomical Union in 2015. This new look on such modern measurements of the Sun's global changes from 1996 to 2017 gives a new way for peering into the solar interior, mainly to better understand the subsurface fields which play an important role in the implementation of the solar cycles.

The X9.3 flare of September 6, 2017, was the most powerful flare of Solar Cycle 24. It generated strong white-light emission and multiple helioseismic waves (sunquakes). By using data from Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) as well as hard X-ray data from KONUS instrument onboard WIND spacecraft, and Anti-Coincidence System (ACS) onboard the INTERGRAL space observatory, we investigate spatio-temporal dynamics of photospheric emission sources, identify sources of helioseismic waves and compare the flare photospheric dynamics with the hard X-ray (HXR) temporal profiles. The results show that the photospheric flare impacts started to develop in compact regions in close vicinity of the magnetic polarity inversion line (PIL) in the pre-impulsive phase before detection of the HXR emission. The initial photospheric disturbances were localized in the region of strong horizontal magnetic field of the PIL, and, thus, are likely associated with a compact sheared magnetic structure elongated along the PIL. The acoustic egression power maps revealed two primary sources of generation of sunquakes, which were associated with places of the strongest photospheric impacts in the pre-impulsive phase and the early impulsive phase. This can explain the two types of helioseismic waves observed in this flare. Analysis of the high-cadence HMI filtergrams suggests that the flare energy release developed in the form of sequential involvement of compact low-lying magnetic loops that were sheared along the PIL.

The precise shape of the Sun is sensitive to the influence of gravity, differential rotation, local turbulence and magnetic fields. It has been previously shown that the solar shape exhibits asphericity that evolves with the 11-year cycle. Thanks to the capability of the SoHO/MDI and SDO/HMI instruments to observe with an unprecedented accuracy the surface gravity oscillation (f) modes, it is possible to extract information concerning the coefficients of rotational frequency splitting, a1, a3 and a5, that measure the differential rotation, together with the a2, a4 and a6 asphericity coefficients. Analysis of these helioseismology data for almost two solar cycles, from 1996 to 2017, reveals a close correlation of the a1 and a5 coefficients with the solar activity, whilst a3 exhibits a long-term trend and a weak correlation in the current cycle indicating a substantial change of the global rotation, potentially associated with a long-term evolution of the solar cycles. Looking in more details, the asphericity coefficients, a2, a4 and a6 are more strongly associated with the solar cycle when applying a time lag of respectively 0.1, 1.6 and -1.6 years. The magnitude of a6-coefficient varies in phase with the sunspot number (SN), but its amplitude is ahead of the SN variation. The last measurements made in mid 2017 indicate that the magnitude of a6-coefficient has probably reached its minimum; therefore, the next solar minimum can be expected by the end of 2018 or in the beginning of 2019. The so-called seismic radius in the range of f-mode angular degree: l=137-299 exhibits a temporal variability in anti-phase with the solar activity; its relative value decreased by 2.3E-05 in Solar Cycle 23 and 1.7E-05 in Cycle 24. Such results will be useful for better understanding the physical mechanisms which act inside the Sun, and so, better constrain dynamo models for forecasting the solar cycles.

Recent advances in helioseismology, numerical simulations and mean-field theory of solar differential rotation have shown that the meridional circulation pattern may consist of two or more cells in each hemisphere of the convection zone. According to the mean-field theory the double-cell circulation pattern can result from the sign inversion of a nondiffusive part of the radial angular momentum transport (the so-called $\Lambda$-effect) in the lower part of the solar convection zone. Here, we show that this phenomenon {can result} from the radial inhomogeneity of the Coriolis number, which depends on the convective turnover time. We demonstrate that if this effect is taken into account then the solar-like differential rotation and the double-cell meridional circulation are both reproduced by the mean-field model. The model is consistent with the distribution of turbulent velocity correlations determined from observations by tracing motions of sunspots and large-scale magnetic fields, indicating that these tracers are rooted just below the shear layer.

We study flare processes in the solar atmosphere using observational data for a M1-class flare of June 12, 2014, obtained by New Solar Telescope (NST/BBSO) and Helioseismic Magnetic Imager (HMI/SDO). The main goal is to understand triggers and manifestations of the flare energy release in the photosphere and chromosphere using high-resolution optical observations and magnetic field measurements. We analyze optical images, HMI Dopplergrams and vector magnetograms, and use Non-Linear Force-Free Field (NLFFF) extrapolations for reconstruction of the magnetic topology and electric currents. The NLFFF modelling reveals interaction of two magnetic flux ropes with oppositely directed magnetic field in the PIL. These flux ropes are observed as a compact sheared arcade along the PIL in the high-resolution broad-band continuum images from NST. In the vicinity of PIL, the NST H alpha observations reveal formation of a thin three-ribbon structure corresponding to a small-scale photospheric magnetic arcade. The observational results evidence in favor of location of the primary energy release site in the chromospheric plasma with strong electric currents concentrated near the polarity inversion line. In this case, magnetic reconnection is triggered by the interacting magnetic flux ropes forming a current sheet elongated along the PIL.

Recent observations of the magnetic field in pre-main sequence stars suggest that the magnetic field topology changes as a function of age. The presence of a tachocline could be an important factor in the development of magnetic field with higher multipolar modes. In this work we performed MHD simulations using the EULAG-MHD code to study the magnetic field generation and evolution in models that mimic stars at two evolutionary stages. The stratification for both stellar phases was computed by fitting stellar structure profiles obtained with the ATON stellar evolution code. The first stage is at $1.1 \mathrm{Myr}$, when the star is completely convective. The second stage is at $14 \mathrm{Myrs}$, when the star is partly convective, with a radiative core developed up to $30\%$ of the stellar radius. In this proceedings we present a preliminary analysis of the resulting mean-flows and magnetic field. The mean-flow analysis shown that the star rotate almost rigidly on the fully convective phase, whereas at the partially convective phase there is differential rotation with conical contours of iso-rotation. As for the mean magnetic field both simulations show similarities with respect to the field evolution. However, the topology of the magnetic field is different.

The present paper discusses results of a statistical study of the characteristics of coronal hole (CH) rotation in order to find connections to the internal rotation of the Sun. The goal is to measure CH rotation rates and study their distribution over latitude and their area sizes. In addition, the CH rotation rates are compared with the solar photospheric and inner layer rotational profiles. We study coronal holes observed within $\pm 60$ latitude and longitude degrees from the solar disc centre during the time span from the 1 January 2013 to 20 April 2015, which includes the extended peak of solar cycle 24.We used data created by the Spatial Possibilistic Clustering Algorithm (SPoCA), which provides the exact location and characterisation of solar coronal holes using SDO=AIA 193 {\AA} channel images. The CH rotation rates are measured with four-hour cadence data to track variable positions of the CH geometric centre. North-south asymmetry was found in the distribution of coronal holes: about 60 percent were observed in the northern hemisphere and 40 percent were observed in the southern hemisphere. The smallest and largest CHs were present only at high latitudes. The average sidereal rotation rate for 540 examined CHs is $13:86 (\pm 0:05)$ degrees/d. Conclusions. The latitudinal characteristics of CH rotation do not match any known photospheric rotation profile. The CH angular velocities exceed the photospheric angular velocities at latitudes higher than 35-40 degrees. According to our results, the CH rotation profile perfectly coincides with tachocline and the lower layers of convection zone at around 0.71 $R_{\odot}$; this indicates that CHs may be linked to the solar global magnetic field, which originates in the tachocline region.

Solar flares are complicated physical phenomena that are observable in a broad range of the electromagnetic spectrum, from radiowaves to $\gamma$-rays. For a more comprehensive understanding of flares, it is necessary to perform a combined multi-wavelength analysis using observations from many satellites and ground-based observatories. For efficient data search, integration of different flare lists and representation of observational data, we have developed an Interactive Multi-Instrument Database of Solar Flares (https://solarflare.njit.edu/). The web accessible database is fully functional and allows the user to search for uniquely-identified flare events based on their physical descriptors and availability of observations by a particular set of instruments. Currently, the data from three primary flare lists (GOES, RHESSI and HEK) and a variety of other event catalogs (Hinode, Fermi GBM, Konus-Wind, OVSA flare catalogs, CACTus CME catalog, Filament eruption catalog) and observing logs (IRIS and Nobeyama coverage) are integrated, and an additional set of physical descriptors (temperature and emission measure) is provided along with an observing summary, data links, and multi-wavelength light curves for each flare event since January, 2002. We envision that this new tool will allow researchers to significantly speed up the search of events of interest for statistical and case studies.

The X-class solar flare of October 23, 2012, generated continuum photospheric emission and a strong helioseismic wave ("sunquake") that points to an intensive energy release in the dense part of the solar atmosphere. We study properties of the energy release with high temporal and spatial resolutions, using photospheric data from the Helioseismic Magnetic Imager (HMI) onboard Solar Dynamics Observatory (SDO), and hard X-ray observations made by the Ramaty High-Energy Solar Spectroscopic Imager (RHESSI). For this analysis we use level-1 HMI data (filtergrams), obtained by scanning the Fe I line (6731~\AA) with the time cadence of $\sim 3.6$ s and spatial resolution of $\sim 0.5^{\prime\prime}$ per pixel. It is found that the photospheric disturbances caused by the flare spatially coincide with the region of hard X-ray emission, but are delayed by $\lesssim 4$ seconds. This delay is consistent with predictions of the flare hydrodynamics RADYN models. However, the models fail to explain the magnitude of variations observed by the HMI. The data indicate that the photospheric impact and helioseismic wave might be caused by the electron energy flux substantially higher than that in the current flare radiative hydrodynamic models.

The size of the diameter of the Sun has been debated for a very long time. First tackled by the Greek astronomers from a geometric point of view, an estimate, although incorrect, has been determined, not truly called into question for several centuries. The French school of astronomy, under the impetus of Mouton and Picard in the XVIIth century can be considered as a pioneer in this issue. It was followed by the German school at the end of the XIXth century whose works led to a canonical value established at 959".63 (second of arc). A number of ground-based observations has been made in the second half of the XIXth century leading to controversial results mainly due to the difficulty to disentangle between the solar and atmospheric effects. Dedicated space measurements yield to a very faint dependence of the solar diameter with time. New studies over the entire radiation spectrum lead to a clear relationship between the solar diameter and the wavelength, reflecting the height at which the lines are formed. Thus the absolute value of the solar diameter, which is a reference for many astrophysical applications, must be stated according to the wavelength. Furthermore, notable features of the Near Sub-Surface Layer (NSSL), called the leptocline, can be established in relation to the solar limb variations, mainly through the shape asphericities coefficients. The exact relationship has not been established yet, but recent studies encourage further in-depth investigations of the solar subsurface dynamics, both observationally and by numerical MHD simulations.

The phenomenon of solar "torsional oscillations" (TO) represents migratory zonal flows associated with the solar cycle. These flows are observed on the solar surface and, according to helioseismology, extend through the convection zone. We study the origin of the TO using results from a global MHD simulation of the solar interior that reproduces several of the observed characteristics of the mean-flows and magnetic fields. Our results indicate that the magnetic tension (MT) in the tachocline region is a key factor for the periodic changes in the angular momentum transport that causes the TO. The torque induced by the MT at the base of the convection zone is positive at the poles and negative at the equator. A rising MT torque at higher latitudes causes the poles to speed-up, whereas a declining negative MT torque at the lower latitudes causes the equator to slow-down. These changes in the zonal flows propagate through the convection zone up to the surface. Additionally, our results suggest that it is the magnetic field at the tachocline that modulates the amplitude of the surface meridional flow rather than the opposite as assumed by flux-transport dynamo models of the solar cycle.

Local helioseismology has opened new frontiers in our quest for understanding of the internal dynamics and dynamo on the Sun. Local helioseismology reconstructs subsurface structures and flows by extracting coherent signals of acoustic waves traveling through the interior and carrying information about subsurface perturbations and flows, from stochastic oscillations observed on the surface. The initial analysis of the subsurface flow maps reconstructed from the 5 years of SDO/HMI data by time-distance helioseismology reveals the great potential for studying and understanding of the dynamics of the quiet Sun and active regions, and the evolution with the solar cycle. In particular, our results show that the emergence and evolution of active regions are accompanied by multi-scale flow patterns, and that the meridional flows display the North-South asymmetry closely correlating with the magnetic activity. The latitudinal variations of the meridional circulation speed, which are probably related to the large-scale converging flows, are mostly confined in shallow subsurface layers. Therefore, these variations do not necessarily affect the magnetic flux transport. The North-South asymmetry is also pronounced in the variations of the differential rotation ("torsional oscillations"). The calculations of a proxy of the subsurface kinetic helicity density show that the helicity does not vary during the solar cycle, and that supergranulation is a likely source of the near-surface helicity.

Local helioseismology provides a unique opportunity to investigate the subsurface structure and dynamics of active regions and their effect on the large-scale flows and global circulation of the Sun. We use measurements of plasma flows in the upper convection zone, provided by the Time-Distance Helioseismology Pipeline developed for analysis of solar oscillation data obtained by Helioseismic and Magnetic Imager (HMI) on Solar Dynamics Observatory (SDO), to investigate the subsurface dynamics of emerging active region NOAA 11726. The active region emergence was detected in deep layers of the convection zone about 12 hours before the first bipolar magnetic structure appeared on the surface, and 2 days before the emergence of most of the magnetic flux. The speed of emergence determined by tracking the flow divergence with depth is about 1.4 km/s, very close to the emergence speed in the deep layers. As the emerging magnetic flux becomes concentrated in sunspots local converging flows are observed beneath the forming sunspots. These flows are most prominent in the depth range 1-3 Mm, and remain converging after the formation process is completed. On the larger scale converging flows around active region appear as a diversion of the zonal shearing flows towards the active region, accompanied by formation of a large-scale vortex structure. This process occurs when a substantial amount of the magnetic flux emerged on the surface, and the converging flow pattern remains stable during the following evolution of the active region. The Carrington synoptic flow maps show that the large-scale subsurface inflows are typical for active regions. In the deeper layers (10-13 Mm) the flows become diverging, and surprisingly strong beneath some active regions. In addition, the synoptic maps reveal a complex evolving pattern of large-scale flows on the scale much larger than supergranulation.

We present multidimensional modeling of convection and oscillations in main-sequence stars somewhat more massive than the Sun, using three separate approaches: 1) Using the 3-D planar StellarBox radiation hydrodynamics code to model the envelope convection zone and part of the radiative zone. Our goals are to examine the interaction of stellar pulsations with turbulent convection in the envelope, excitation of acoustic modes, and the role of convective overshooting; 2) Applying the spherical 3-D MHD ASH (Anelastic Spherical Harmonics) code to simulate the core convection and radiative zone. Our goal is to determine whether core convection can excite low-frequency gravity modes, and thereby explain the presence of low frequencies for some hybrid gamma Doradus/delta Scuti variables for which the envelope convection zone is too shallow for the convective blocking mechanism to drive gravity modes; 3) Applying the ROTORC 2-D stellar evolution and dynamics code to calculate evolution with a variety of initial rotation rates and extents of core convective overshooting. The nonradial adiabatic pulsation frequencies of these nonspherical models are calculated using the 2-D pulsation code NRO. We present new insights into pulsations for stars of one to two solar masses gained by multidimensional modeling.

Chromospheric evaporation is observed as Doppler blueshift during solar flares. It plays one of key roles in dynamics and energetics of solar flares, however, its mechanism is still unknown. In this paper we present a detailed analysis of spatially-resolved multi-wavelength observations of chromospheric evaporation during an M 1.0 class solar flare (SOL2014-06-12T21:12) using data from the NASA's IRIS (Interface Region Imaging Spectrograph) and HMI/SDO (Helioseismic and Magnetic Imager onboard Solar Dynamics Observatory) telescopes, and VIS/NST (Visible Imaging Spectrometer at New Solar Telescope) high-resolution observations, covering the temperature range from 10^4 K to 10^7 K. The results show that the averaged over the region Fe XXI blueshift of the hot evaporating plasma is delayed relative to the C II redshift of the relatively cold chromospheric plasma by about 1 min. The spatial distribution of the delays is not uniform across the region and can be as long as 2 min in several zones. Using vector magnetograms from HMI we reconstruct the magnetic field topology and the quasi-separatrix layer (QSL) and find that the blueshift delay regions as well as the H-alpha flare ribbons are connected to the region of magnetic polarity inversion line (PIL) and an expanding flux rope via a system of low-lying loop arcades with height < ~4.5 Mm. This allows us to propose an interpretation of the chromospheric evaporation based on the geometry of local magnetic fields, and the primary energy source associated with the PIL.

We study flare processes in the lower solar atmosphere using observational data for a M1-class flare of June 12, 2014, obtained by New Solar Telescope (NST/BBSO) and Helioseismic Magnetic Imager (HMI/SDO). The main goal is to understand triggers and manifestations of the flare energy release in the lower layers of the solar atmosphere (the photosphere and chromosphere) using high-resolution optical observations and magnetic field measurements. We analyze optical images, HMI Dopplergrams and vector magnetograms, and use Non-Linear Force-Free Field (NLFFF) extrapolations for reconstruction of the magnetic topology. The NLFFF modelling reveals interaction of oppositely directed magnetic flux-tubes in the PIL. These two interacting magnetic flux tubes are observed as a compact sheared arcade along the PIL in the high-resolution broad-band continuum images from NST. In the vicinity of the PIL, the NST H alpha observations reveal formation of a thin three-ribbon structure corresponding to the small-scale photospheric magnetic arcade. Presented observational results evidence in favor of location of the primary energy release site in the dense chromosphere where plasma is partially ionized in the region of strong electric currents concentrated near the polarity inversion line. Magnetic reconnection may be triggered by two interacting magnetic flux tubes with forming current sheet elongated along the PIL.

We study turbulent generation of large-scale magnetic fields using nonlinear dynamo models for solar-type stars in the range of rotational periods from 14 to 30 days. Our models take into account non-linear effects of dynamical quenching of magnetic helicity, and escape of magnetic field from the dynamo region due to magnetic buoyancy. The results show that the observed correlation between the period of rotation and the duration of activity cycles can be explained in the framework of a distributed dynamo model with a dynamical magnetic feedback acting on the turbulent generation either from magnetic buoyancy or magnetic helicity. We discuss implications of our findings for the understanding of dynamo processes operating in solar-like stars.

Using mean-field hydrodynamic models of the solar angular momentum balance we show that the non-monotonic latitudinal dependence of the radial angular momentum fluxes caused by Lambda-effect can affect the number of the meridional circulation cells stacking in radial direction in the solar convection zone. In particular, our results show the possibility of a complicated triple-cell meridional circulation structure. This pattern consists of two large counterclockwise circulation cells (the N-hemisphere) and a smaller clockwise cell located at low latitudes at the bottom of the convection zone.

Continued progress in observational stellar astrophysics requires a deep understanding of the underlying convection dynamics. We present results of realistic 3D radiative hydrodynamic simulations of the outer layers of a moderate mass star (1.47 Msun), including the full convection zone, the overshoot region, and the top layers of the radiative zone. The simulation results show that the surface granulation has a broad range of scales, from 2 to 12 Mm, and that large granules are organized in well-defined clusters, consisting of several granules. Comparison of the mean structure profiles from 3D simulations with the corresponding 1D standard stellar model shows an increase of the stellar radius by ~800 km, as well as significant changes in the thermodynamic structure and turbulent properties of the ionization zones. Convective downdrafts in the intergranular lanes between granulation clusters reach speeds of more than 20 km/s, penetrate through the whole convection zone, hit the radiative zone, and form a 8 Mm thick overshoot layer. Contrary to semi-empirical overshooting models, our results show that the 3D dynamic overshoot region consists of two layers: a nearly adiabatic extension of the convection zone and a deeper layer of enhanced subadiabatic stratification. This layer is formed because of heating caused by the braking of the overshooting convective plumes. This effect has to be taken into account in stellar modeling and the interpretation of asteroseismology data. In particular, we demonstrate that the 3D model can qualitatively explain deviations from the standard solar model found by helioseismology.

Rotational shear layers at the boundary between radiative and convective zones, tachoclines, play a key role in the process of magnetic field generation in solar-like stars. We present two sets of global simulations of rotating turbulent convection and dynamo. The first set considers a stellar convective envelope only; the second one, aiming at the formation of a tachocline, considers also the upper part of the radiative zone. Our results indicate that the resulting mean-flows and dynamo properties like the growth rate, saturation energy and mode depend on the Rossby (Ro) number. For the first set of models either oscillatory (with ~2 yr period) or steady dynamo solutions are obtained. The models in the second set naturally develop a tachocline which, in turn, leads to the generation of strong mean magnetic field. Since the field is also deposited into the stable deeper layer, its evolutionary time-scale is much longer than in the models without a tachocline. Surprisingly, the magnetic field in the upper turbulent convection zone evolves in the same time scale as the deep field. These models result in either an oscillatory dynamo with ~30 yr period or in a steady dynamo depending on Ro. In terms of the mean-field dynamo coefficients computed using FOSA, the field evolution in the oscillatory models without a tachocline seems to be consistent with dynamo waves propagating according to the Parker-Yoshimura sign rule. In the models with tachoclines the dynamics is more complex involving other transport mechanisms as well as tachocline instabilities.

Numerical MHD simulations play increasingly important role for understanding mechanisms of stellar magnetism. We present simulations of convection and dynamos in density-stratified rotating spherical fluid shells. We employ a new 3D simulation code for the solution of a physically consistent anelastic model of the process with a minimum number of parameters. The reported dynamo simulations extend into a "buoyancy-dominated" regime where the buoyancy forcing is dominant while the Coriolis force is no longer balanced by pressure gradients and strong anti-solar differential rotation develops as a result. We find that the self-generated magnetic fields, despite being relatively weak, are able to reverse the direction of differential rotation from anti-solar to solar-like. We also find that convection flows in this regime are significantly stronger in the polar regions than in the equatorial region, leading to non-oscillatory dipole-dominated dynamo solutions, and to concentration of magnetic field in the polar regions. We observe that convection has different morphology in the inner and at the outer part of the convection zone simultaneously such that organized geostrophic convection columns are hidden below a near-surface layer of well-mixed highly-chaotic convection. While we focus the attention on the buoyancy-dominated regime, we also demonstrate that conical differential rotation profiles and persistent regular dynamo oscillations can be obtained in the parameter space of the rotation-dominated regime even within this minimal model.

We explore a response of a non-linear non-axisymmetric mean-field solar dynamo model to shallow non-axisymmetric perturbations. After a relaxation period the amplitude of the non-axisymmetric field depends on the initial condition, helicity conservation, and the depth of perturbation. It is found that a perturbation which is anchored at r=0.9R has a profound effect on the dynamo process, producing a transient magnetic cycle of the axisymmetric magnetic field, if it is initiated at the growing phase of the cycle. The non-symmetric with respect to the equator perturbation results in a hemispheric asymmetry of the magnetic activity. The evolution of the axisymmetric and non-axisymmetric field depends on the turbulent magnetic Reynolds number R_m. In the range of R_m=10^{4-6} the evolution returns to the normal course in the next cycle, in which the non-axisymmetric field is generated due to a non-linear alpha-effect and magnetic buoyancy. In the stationary state the large-scale magnetic field demonstrates a phenomenon of "active longitudes" with cyclic 180 degree "flip-flop" changes of the large-scale magnetic field orientation. The flip-flop effect is known from observations of solar and stellar magnetic cycles. However this effect disappears in the model which includes the meridional circulation pattern determined by helioseismology. The rotation rate of the non-axisymmetric field components varies during the relaxation period, and carries important information about the dynamo process.

Realistic numerical simulations, i.e., those that make minimal use of ad hoc modeling, are essential for understanding the complex turbulent dynamics of the interiors and atmospheres of the Sun and other stars and the basic mechanisms of their magnetic activity and variability. The goal of this paper is to present a detailed description and test results of a compressible radiative MHD code, `StellarBox', specifically developed for simulating the convection zones, surface, and atmospheres of the Sun and moderate-mass stars. The code solves the three-dimensional, fully coupled compressible MHD equations using a fourth-order Pad\'e spatial differentiation scheme and a fourth-order Runge-Kutta scheme for time integration. The radiative transfer equation is solved using the Feautrier method for bi-directional ray tracing and an opacity-binning technique. A specific feature of the code is the implementation of subgrid-scale MHD turbulence models. The data structures are automatically configured, depending on the computational grid and the number of available processors, to achieve good load balancing. We present test results and illustrate the code's capabilities for simulating the granular convection on the Sun and a set of main-sequence stars. The results reveal substantial changes in the near-surface turbulent convection in these stars, which in turn affect properties of the surface magnetic fields. For example, in the solar case initially uniform vertical magnetic fields tend to self-organize into compact (pore-like) magnetic structures, while in more massive stars such structures are not formed and the magnetic field is distributed more-or-less uniformly in the intergranular lanes.

Helioseismic and magnetohydrodynamic waves are abundant in and above sunspots. Through cross-correlating oscillation signals in the photosphere observed by the SDO/HMI, we reconstruct how waves propagate away from virtual wave sources located inside a sunspot. In addition to the usual helioseismic wave, a fast-moving wave is detected traveling along the sunspot's radial direction from the umbra to about 15 Mm beyond the sunspot boundary. The wave has a frequency range of 2.5 - 4.0 mHz with a phase velocity of 45.3 km/s, substantially faster than the typical speeds of Alfven and magnetoacoustic waves in the photosphere. The observed phenomenon is consistent with a scenario of that a magnetoacoustic wave is excited at approximately 5 Mm beneath the sunspot, and its wavefront travels to and sweeps across the photosphere with a speed higher than the local magnetoacoustic speed. The fast-moving wave, if truly excited beneath the sunspot's surface, will help open a new window to study the internal structure and dynamics of sunspots.

Magnetic fields are usually observed in the quiet Sun as small-scale elements that cover the entire solar surface (the `salt and pepper' patterns in line-of-sight magnetograms). By using 3D radiative MHD numerical simulations we find that these fields result from a local dynamo action in the top layers of the convection zone, where extremely weak 'seed' magnetic fields (e.g., from a $10^{-6}$ G) can locally grow above the mean equipartition field, to a stronger than 2000~G field localized in magnetic structures. Our results reveal that the magnetic flux is predominantly generated in regions of small-scale helical downflows. We find that the local dynamo action takes place mostly in a shallow, about 500~km deep, subsurface layer, from which the generated field is transported into the deeper layers by convective downdrafts. We demonstrate that the observed dominance of vertical magnetic fields at the photosphere and horizontal fields above the photosphere can be explained by small-scale magnetic loops produced by the dynamo. Such small-scale loops play an important role in the structure and dynamics of the solar atmosphere and that their detection in observations is critical for understanding the local dynamo action on the Sun.

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.

Dynamics of hot chromospheric plasma of solar flares is a key to understanding of mechanisms of flare energy release and particle acceleration. A moderate M1.0 class flare of 12 June, 2014 (SOL2014-06-12T21:12) was simultaneously observed by NASA's Interface Region Imaging Spectrograph (IRIS), other spacecraft, and also by New Solar Telescope (NST) at the BBSO. This paper presents the first part of our investigation focused on analysis of the IRIS data. Our analysis of the IRIS data in different spectral lines reveals strong redshifted jet-like flow with the speed of ~100 km/s of the chromospheric material before the flare. Strong nonthermal emission of the C II k 1334.5 A line, formed in the chromosphere-corona transition region, is observed at the beginning of the impulsive phase in several small (with a size of ~1 arcsec) points. It is also found that the C II k line is redshifted across the flaring region before, during and after the impulsive phase. A peak of integrated emission of the hot (1.1 MK) plasma in the Fe XXI 1354.1 A line is detected approximately 5 minutes after the integrated emission peak of the lower temperature C II k. A strong blueshift of the Fe XXI line across the flaring region corresponds to evaporation flows of the hot chromospheric plasma with a speed of 50 km/s. Additional analysis of the Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI) data supports the idea that the upper chromospheric dynamics observed by IRIS has features of "gentle" evaporation driven by heating of the solar chromosphere by accelerated electrons and by a heat flux from the flare energy release site.

Two different mechanisms may act to induce quasi-periodic pulsations (QPP) in whole-disk observations of stellar flares. One mechanism may be magneto-hydromagnetic (MHD) forces and other processes acting on flare loops as seen in the Sun. The other mechanism may be forced local acoustic oscillations due to the high-energy particle impulse generated by the flare (known as `sunquakes' in the Sun). We analyze short-cadence Kepler data of 257 flares in 75 stars to search for QPP in the flare decay branch or post-flare oscillations which may be attributed to either of these two mechanisms. About 18 percent of stellar flares show a distinct bump in the flare decay branch of unknown origin. The bump does not seem to be a highly-damped global oscillation because the periods of the bumps derived from wavelet analysis do not correlate with any stellar parameter. We detected damped oscillations covering several cycles (QPP), in seven flares on five stars. The periods of these oscillations also do not correlate with any stellar parameter, suggesting that these may be a due to flare loop oscillations. We searched for forced global oscillations which might result after a strong flare. To this end, we investigated the behaviour of the amplitudes of solar-like oscillations in eight stars before and after a flare. However, no clear amplitude change could be detected. We also analyzed the amplitudes of the self-excited pulsations in two delta Scuti stars and one gamma Doradus star before and after a flare. Again, no clear amplitude changes were found. Our conclusions are that a new process needs to be found to explain the high incidence of bumps in stellar flare light curves, that flare loop oscillations may have been detected in a few stars and that no conclusive evidence exists as yet for flare induced global acoustic oscillations (starquakes).

The solar flare on July 30, 2011 was of a modest X-ray class (M9.3), but it made a strong photospheric impact and produced a "sunquake," observed with the Helioseismic and Magnetic Imager (HMI) on NASA's Solar Dynamics Observatory (SDO). In addition to the helioseismic waves (also observed with the SDO/AIA instrument), the flare caused a large expanding area of white-light emission and was accompanied by substantial restructuring of magnetic fields, leading to the rapid formation of a sunspot structure in the flare region. The flare produced no significant hard X-ray emission and no coronal mass ejection. This indicates that the flare energy release was mostly confined to the lower atmosphere. The absence of significant coronal mass ejection rules out magnetic rope eruption as a mechanism of helioseismic waves. We discuss the connectivity of the flare energy release with the electric currents dynamics and show the potential importance of high-speed plasma flows in the lower solar atmosphere during the flare energy release.

We present results of non-LTE modeling of H-alpha 6563 A and Mg II k&h 2796 A and 2803 A lines. This modeling is important for interpretation of coordinated observations from the recently launched NASA's IRIS mission and from the New Solar Telescope at Big Bear Solar Observatory. Among available codes for the non-LTE modeling, the RH code is chosen as the most appropriate for modeling of the line profiles. The most suitable Hydrogen and Magnesium atomic models are selected by performing several tests of the code. The influence of the ionization degree on the line profiles is also studied. Radiative-MHD simulations of the solar atmosphere, obtained with the Bifrost code, are used as input data for calculation of synthetic spectra of the H-alpha and Mg II h&k lines for particular locations evolving with time. The spectral line variations reveal the presence of MHD waves in the simulation results. We construct oscillation power spectra of the line intensity for different wavelength, and compare these with the corresponding height-dependent power spectra of atmospheric parameters from the simulations. We find correlations between the power spectra of intensities of the line profiles at certain wavelengths and the power spectra of the atmospheric parameters at the tau-unity heights for these wavelengths. These results provide a new diagnostic method of chromospheric oscillations; however, larger amounts of data are needed to confirm these correlations.

We present analysis of C7.0 solar flare of Febrary 17, 2013, revealing a strong helioseismic response (sunquake) caused by a very compact impact in the photosphere. This is the weakest known C-class flare generating a sunquake event. To investigate possible mechanisms of this event, and to understand the role of accelerated charged particles and photospheric electric currents, we use data from three space observatories: Ramaty High Energy Solar Spectroscopic Imager (RHESSI), Solar Dynamics Observatory (SDO) and Geostationary Operational Environmental Satellite (GOES). We find that the photospheric flare impact does not spatially correspond to the strongest HXR emission source, but both of these events are parts of the same energy release. Our analysis reveals a close association of the flare energy release with a rapid increase of the electric currents, and suggests that the sunquake initiation is unlikely to be explained by the impact of high-energy electrons but may be associated with a rapid current dissipation or a localized impulsive Lorentz force.

Solar subsurface zonal and meridional-flow profiles during the rising phase of solar cycle 24 are studied using time-distance helioseismology technique. The faster zonal bands in the torsional-oscillation pattern show strong hemispheric asymmetries and temporal variations in both width and speed. The faster band in the northern hemisphere is located closer to the equator than the band in the southern hemisphere, and migrates past the equator when the magnetic activity in the southern hemisphere is reaching maximum. The meridional-flow speed decreases substantially with the increase of magnetic activity, and the flow profile shows two zonal structures in each hemisphere. The residual meridional flow, after subtracting a mean meridional-flow profile, converges toward the activity belts and shows faster and slower bands like the torsional-oscillation pattern. More interestingly, the meridional-flow speed above latitude 30 degree shows an anti-correlation with the poleward-transporting magnetic flux, slower when the following-polarity flux is transported and faster when the leading-polarity flux is transported. It is expected that this phenomenon slows the process of magnetic cancellation and polarity reversal in the high-latitude areas.

Emission of solar flares across the electromagnetic spectrum is often observed in the form of two expanding ribbons. The standard flare model explains the flare ribbons as footpoints of magnetic arcades, emitting due to interaction of energetic particles with the chromospheric plasma. However, the physics of this interaction and properties of the accelerated particles are still unknown. We present results of multiwavelength observations of C2.1 flare of August 15, 2013, observed with New Solar Telescope (NST) of Big Bear Solar Observatory, Solar Dynamics Observatory (SDO), GOES and FERMI spacecraft. The observations reveal previously unresolved sub-arcsecond structure of the flare ribbons in regions of strong magnetic field consisting from numerous small-scale bright knots. We observe red-blue asymmetry of H alpha flare ribbons with a width as small as 100 km. We discuss the relationship between the ribbons and vertical electric currents estimated from vector magnetograms, and show that Joule heating can be responsible for energization of H alpha knots in the ribbons.

State-of-the-art solar instrumentation is revealing magnetic activity of the Sun with unprecedented resolution. Observations with the 1.6m New Solar Telescope of the Big Bear Solar Observatory are making next steps in our understanding of the solar surface structure. Granular-scale magnetic flux emergence and the response of the solar atmosphere are among the key research. As part of a joint observing program with NASA's IRIS mission, the NST observed active region NOAA 11810 in photospheric and chromospheric wavelengths. Complimentary data are provided by SDO and Hinode space-based telescopes. The region displayed a group of solar pores, in the vicinity of which we detect a small-scale buoyant horizontal magnetic flux tube causing abnormal granulation and interacting with the pre-existing ambient field in upper atmospheric layers. Following the expansion of distorted granules at the emergence site, we observed a sudden appearance of an extended surge in the HeI data. IRIS catched ejection of a hot plasma jet associated with the HeI-surge. The SDO/HMI data reveal emerging magnetic looplike structures. Hinode/Ca II H and IRIS filtergrams detail the connectivities of the newly emerged magnetic field in the lower solar chromosphere. From these data we fid that the orientation of the emerging flux tube was almost perpendicular to the overlying ambient field. Nevertheless the interaction of emerging magnetic field lines with the pre-existing overlying field generates high-temperature emission regions and boosts the surge/jet production. The localized heating is detected before and after the first signs of the surge/jet ejection. We compare the results with previous observations and theoretical models, and propose a scenario for the activation of plasma jet/surges and confined heating. Such process may play a significant role in the mass and energy flow from the interior to the corona.

Helioseismology provides important constraints for the solar dynamo problem. However, the basic properties and even the depth of the dynamo process, which operates also in other stars, are unknown. Most of the dynamo models suggest that the toroidal magnetic field that emerges on the surface and forms sunspots is generated near the bottom of the convection zone, in the tachocline. However, there is a number of theoretical and observational problems with justifying the deep-seated dynamo models. This leads to the idea that the subsurface angular velocity shear may play an important role in the solar dynamo. Using helioseismology measurements of the internal rotation and meridional circulation, we investigate a mean-field MHD model of dynamo distributed in the bulk of the convection zone but shaped in a near-surface layer. We show that if the boundary conditions at the top of the dynamo region allow the large-scale toroidal magnetic fields to penetrate into the surface, then the dynamo wave propagates along the isosurface of angular velocity in the subsurface shear layer, forming the butterfly diagram in agreement with the Parker-Yoshimura rule and solar-cycle observations. Unlike the flux-transport dynamo models, this model does not depend on the transport of magnetic field by meridional circulation at the bottom of the convection zone, and works well when the meridional circulation forms two cells in radius, as recently indicated by deep-focus time-distance helioseismology analysis of the SDO/HMI and SOHO/MDI data. We compare the new dynamo model with various characteristics if the solar magnetic cycles, including the cycle asymmetry (Waldmeier's relations) and magnetic `butterfly' diagrams.

Sunquakes observed in the form of expanding wave ripples on the surface of the Sun during solar flares represent packets of acoustic waves excited by flare impacts and traveling through the solar interior. The excitation impacts strongly correlate with the impulsive flare phase, and are caused by the energy and momentum transported from the energy release sites. The flare energy is released in the form of energetic particles, waves, mass motions, and radiation. However, the exact mechanism of the localized hydrodynamic impacts which generate sunquakes is unknown. Solving the problem of the sunquake mechanism will substantially improve our understanding of the flare physics. In addition, sunquakes offer a unique opportunity for studying the interaction of acoustic waves with magnetic fields and flows in flaring active regions, and for developing new approaches to helioseismic acoustic tomography.

In addition to well-known mechanisms of excitation of solar and stellar oscillations by turbulent convection and instabilities, the oscillations can be excited by an impulsive localized force caused by the energy release in solar and stellar flares. Such oscillations have been observed on the Sun (`sunquakes'), and created a lot of interesting discussions about physical mechanisms of the impulsive excitation and their relationship to the flare physics. The observation and theory have shown that most of a sunquake's energy is released in high-degree, high-frequency p modes. In addition, there have been reports on helioseismic observations of low-degree modes excited by strong solar flares. Much more powerful flares observed on other stars can cause `starquakes' of substantially higher amplitude. Observations of such oscillations can provide new asteroseismic information and also constraints on mechanisms of stellar flares. The basic properties of sunquakes and initial attempts to detect flare-excited oscillations in Kepler short-cadence data are discussed.

In the quiet Sun, magnetic fields are usually observed as small-scale magnetic elements, `salt and pepper', covering the entire solar surface. By using 3D radiative MHD numerical simulations we demonstrate that these fields are a result of local dynamo action in the top layers of the convection zone, where extremely weak `seed' magnetic fields can locally grow above the mean equipartition field (e.g., from a $10^{-6}$ G `seed' field to more than 1000 G magnetic structures). We find that the local dynamo action takes place only in a shallow, about 500 km deep, subsurface layer, from which the generated field is transported into deeper layers by convection downdrafts. We demonstrate that the observed dominance of vertical magnetic fields at the photosphere and the horizontal fields above the photosphere can be explained by multi-scale magnetic loops produced by the dynamo.

We study how anisotropy of turbulent convection affects diffusion of large-scale magnetic fields and the dynamo process on the Sun. The effect of anisotropy is calculated in a mean-field magneto-hydrodynamics framework using the minimal $\tau$-approximation. We examine two types of mean-field dynamo models: the well-known benchmark flux-transport model, and a distributed-dynamo model with the subsurface rotational shear layer. For both models we investigate effects of the double-cell meridional circulation, recently suggested by helioseismology. We introduce a parameter of anisotropy as a ratio of the radial and horizontal intensity of turbulent mixing, to characterize the anisotropy effects. It is found that the anisotropy of turbulent convection affects the distribution of magnetic fields inside the convection zone. The concentration of the magnetic flux near the bottom and top boundaries of the convection zone is greater when the anisotropy is stronger. It is shown that the critical dynamo number and the dynamo period approach to constant values for the large anisotropy parameter. The anisotropy reduces the overlap of the toroidal fields of subsequent cycles in the time-latitude ``butterfly' diagram. If we assume that sunspots are formed in the vicinity of the subsurface shear layer, then the distributed dynamo model with anisotropic diffusivity satisfies the observational constraints from heloseismology results and is consistent with the value of effective turbulent diffusion, estimated from the dynamics of surface magnetic fields.

To explore the physics of large-scale flows in solar-like stars, we perform 3D anelastic simulations of rotating convection for global models with stratification resembling the solar interior. The numerical method is based on an implicit large-eddy simulation approach designed to capture effects from non-resolved small scales. We obtain two regimes of differential rotation, with equatorial zonal flows accelerated either in the direction of rotation (solar-like) or in the opposite direction (anti-solar). While the models with the solar-like differential rotation tend to produce multiple cells of meridional circulation, the models with anti-solar differential rotation result in only one or two meridional cells. Our simulations indicate that the rotation and large-scale flow patterns critically depend on the ratio between buoyancy and Coriolis forces. By including a subadiabatic layer at the bottom of the domain, corresponding to the stratification of a radiative zone, we reproduce a layer of strong radial shear similar to the solar tachocline. Similarly, enhanced superadiabaticity at the top results in a near-surface shear layer located mainly at lower latitudes. The models reveal a latitudinal entropy gradient localized at the base of the convection zone and in the stable region, which however does not propagate across the convection zone. In consequence, baroclinicity effects remain small and the rotation iso-contours align in cylinders along the rotation axis. Our results confirm the alignment of large convective cells along the rotation axis in the deep convection zone, and suggest that such "banana-cell" pattern can be hidden beneath the supergranulation layer.

Meridional flow in the solar interior plays an important role in redistributing angular momentum and transporting magnetic flux inside the Sun. Although it has long been recognized that the meridional flow is predominantly poleward at the Sun's surface and in its shallow interior, the location of the equatorward return flow and the meridional flow profile in the deeper interior remain unclear. Using the first two years of continuous helioseismology observations from the Solar Dynamics Observatory / Helioseismic Magnetic Imager, we analyze travel times of acoustic waves that propagate through different depths of the solar interior carrying information about the solar interior dynamics. After removing a systematic center-to-limb effect in the helioseismic measurements and performing inversions for flow speed, we find that the poleward meridional flow of a speed of 15 m/s extends in depth from the photosphere to about 0.91 R_sun. An equatorward flow of a speed of 10 m/s is found between 0.82 to 0.91 R_sun in the middle of the convection zone. Our analysis also shows evidence of that the meridional flow turns poleward again below 0.82 R_sun, indicating an existence of a second meridional circulation cell below the shallower one. This double-cell meridional circulation profile with an equatorward flow shallower than previously thought suggests a rethinking of how magnetic field is generated and redistributed inside the Sun.

Recent findings of helioseismology as well as advances in direct numerical simulations of global dynamics of the Sun have indicated that in each solar hemisphere the meridional circulation forms the two cells along the in the convection zone. We investigate properties of a mean-field solar dynamo with such double-cell meridional circulation. The dynamo model also includes the realistic profile of solar differential rotation (including the tachocline and subsurface shear layer), and takes into account effects of turbulent pumping, anisotropic turbulent diffusivity, and conservation of magnetic helicity. Contrary to previous flux-transport dynamo models, it is found that the dynamo model can robustly reproduce the basic properties of the solar magnetic cycles for a wide range of model parameters and the circulation speed. The best agreement with observations is achieved when the surface speed of meridional circulation is about 12 m/s. For this circulation speed the simulated sunspot activity shows good synchronization with the polar magnetic fields. Such synchronization was indeed observed during the past sunspot cycles 21 and 22. We compare theoretical and observed phase diagrams of the sunspot number and the polar field strength and discuss the peculiar properties of Cycle 23.

Convective turbulent motions in the solar interior, as well as the mean flows resulting from them, determine the evolution of the solar magnetic field. With the aim to get a better understanding of these flows we study anelastic rotating convection in a spherical shell whose stratification resembles that of the solar interior. This study is done through numerical simulations performed with the EULAG code. Due to the numerical formulation, these simulations are known as implicit large eddy simulations (ILES), since they intrinsically capture the contribution of, non-resolved, small scales at the same time maximizing the effective Reynolds number. We reproduce some previous results and find a transition between buoyancy and rotation dominated regimes which results in anti-solar or solar like rotation patterns. Even thought the rotation profiles are dominated by Taylor-Proudman columnar rotation, we are able to reproduce the tachocline and a low latitude near-surface shear layer. We find that simulations results depend on the grid resolution as a consequence of a different sub-grid scale contribution.

The solar surface is covered by high-speed jets transporting mass and energy into the solar corona and feeding the solar wind. The most prominent of these jets have been known as spicules. However, the mechanism initiating these eruptions events is still unknown. Using realistic numerical simulations we find that small-scale eruptions are produced by ubiquitous magnetized vortex tubes generated by the Sun's turbulent convection in subsurface layers. The swirling vortex tubes (resembling tornadoes) penetrate into the solar atmosphere, capture and stretch background magnetic field, and push surrounding material up, generating quasiperiodic shocks. Our simulations reveal a complicated high-speed flow patterns, and thermodynamic and magnetic structure in the erupting vortex tubes. We found that the eruptions are initiated in the subsurface layers and are driven by the high-pressure gradients in the subphotosphere and photosphere, and by the Lorentz force in the higher atmosphere layers.

This article presents a comparative analysis of solar activity data, Mt Wilson diameter data, Super-Kamiokande solar neutrino data, and nuclear decay data acquired at the Lomonosov Moscow State University (LMSU). We propose that salient periodicities in all of these datasets may be attributed to r-mode oscillations. Periodicities in the solar activity data and in Super-Kamiokande solar neutrino data may be attributed to r-mode oscillations in the known tachocline, with normalized radius in the range 0.66 to 0.74, where the sidereal rotation rate is in the range 13.7 to 14.6 year-1. We propose that periodicities in the Mt Wilson and LMSU data may be attributed to similar r-mode oscillations where the sidereal rotation rate is approximately which we attribute to a hypothetical inner tachocline separating a slowly rotating core from the radiative zone. We also discuss the possible role of the RSFP (Resonant Spin Flavor Precession) process, which leads to estimates of the neutrino magnetic moment and of the magnetic field strength in or near the solar core.

Solar five-minute oscillations have been detected in the power spectra of two six-day time intervals from soft X-ray measurements of the Sun observed as a star using the Extreme Ultraviolet Spectrophotometer (ESP) onboard the Solar Dynamics Observatory (SDO) Extreme Ultraviolet Variability Experiment (EVE). The frequencies of the largest amplitude peaks were found matching within 3.7 microHz the known low-degree (l = 0--3) modes of global acoustic oscillations, and can be explained by a leakage of the global modes into the corona. Due to strong variability of the solar atmosphere between the photosphere and the corona the frequencies and amplitudes of the coronal oscillations are likely to vary with time. We investigate the variations in the power spectra for individual days and their association with changes of solar activity, e.g. with the mean level of the EUV irradiance, and its short-term variations due to evolving active regions. Our analysis of samples of one-day oscillation power spectra for a 49-day period of low and intermediate solar activity showed little correlation with the mean EUV irradiance and the short-term variability of the irradiance. We suggest that some other changes in the solar atmosphere, e.g. magnetic fields and/or inter-network configuration may affect the mode leakage to the corona.

We performed 3D numerical simulations of the solar surface wave field for the quiet Sun and for three models with different localized sound-speed variations in the interior with: (i) deep, (ii) shallow, and (iii) two-layer structures. We used simulated data generated by two different codes which use the same standard solar model as a background model, but utilize two different integration techniques and use different models of stochastic wave excitation. Acoustic travel times were measured from all data sets using the time-distance helioseismology technique and compared with the ray theory predictions, frequently used for helioseismic travel-time inversions. It is found that the measured travel-time shifts agree well with the ray theory in both cases with and without phase-speed filtering for the shallow and deep perturbations. This testing verifies the whole measuring-filtering-inversion procedure for sound-speed anomalies inside the Sun. It is shown, that the phase-speed filtering, frequently used to improve the signal-to-noise ratio does not introduce significant systematic errors. Results of the sound-speed inversion procedure show good agreement with the background sound-speed profiles in all cases. Due to its smoothing nature, the inversion procedure overestimates sound speed variations in areas with sharp gradients of the sound-speed profile.

The meridional flow in the Sun is an axisymmetric flow that is generally poleward directed at the surface, and is presumed to be of fundamental importance in the generation and transport of magnetic fields. Its true shape and strength, however, is debated. We present a numerical simulation of helioseismic wave propagation in the whole solar interior in the presence of a prescribed, stationary, single-cell, deep meridional circulation serving as a test-bed for helioseismic measurement techniques. A deep-focusing time-distance helioseismology technique is applied to the artificial data showing that it can in fact be used to measure the effects of the meridional flow very deep in the solar convection zone. It is shown that the ray-approximation which is commonly used for interpretation of helioseismology measurements remains a reasonable approximation even for the very long distances between 12 and 42 degrees corresponding to depths between 52 and 195 Mm considered here. From the measurement noise we extrapolate that on the order of a full solar cycle may be needed to probe the flow all the way to the base of the convection zone.

Results of a statistical analysis of solar granulation are presented. A data set of 36 images of a quiet Sun area on the solar disk center was used. The data were obtained with the 1.6 m clear aperture New Solar Telescope (NST) at Big Bear Solar Observatory (BBSO) and with a broad-band filter centered at the TiO (705.7 nm) spectral line. The very high spatial resolution of the data (diffraction limit of 77 km and pixel scale of 0.$"$0375) augmented by the very high image contrast (15.5$\pm$0.6%) allowed us to detect for the first time a distinct subpopulation of mini-granular structures. These structures are dominant on spatial scales below 600 km. Their size is distributed as a power law with an index of -1.8 (which is close to the Kolmogorov's -5/3 law) and no predominant scale. The regular granules display a Gaussian (normal) size distribution with a mean diameter of 1050 km. Mini-granular structures contribute significantly to the total granular area. They are predominantly confined to the wide dark lanes between regular granules and often form chains and clusters, but different from magnetic bright points. A multi-fractality test reveals that the structures smaller than 600 km represent a multi-fractal, whereas on larger scales the granulation pattern shows no multi-fractality and can be considered as a Gaussian random field. The origin, properties and role of the newly discovered population of mini-granular structures in the solar magneto-convection are yet to be explored.

Turbulent properties of the quiet Sun represent the basic state of surface conditions, and a background for various processes of solar activity. Therefore understanding of properties and dynamics of this `basic' state is important for investigation of more complex phenomena, formation and development of observed phenomena in the photosphere and atmosphere. For characterization of the turbulent properties we compare kinetic energy spectra on granular and sub-granular scales obtained from infrared TiO observations with the New Solar Telescope (Big Bear Solar Observatory) and from 3D radiative MHD numerical simulations ('SolarBox' code). We find that the numerical simulations require a high spatial resolution with 10 - 25 km grid-step in order to reproduce the inertial (Kolmogorov) turbulence range. The observational data require an averaging procedure to remove noise and potential instrumental artifacts. The resulting kinetic energy spectra show a good agreement between the simulations and observations, opening new perspectives for detailed joint analysis of more complex turbulent phenomena on the Sun, and possibly on other stars. In addition, using the simulations and observations we investigate effects of background magnetic field, which is concentrated in self-organized complicated structures in intergranular lanes, and find an increase of the small-scale turbulence energy and its decrease at larger scales due to magnetic field effects.

Investigating the properties of magnetic flux emergence is one of the most important problems of solar physics. In this study we present a newly developed deep-focus time-distance measurement scheme which is able to detect strong emerging flux events in the deep solar interior, before the flux becomes visible on the surface. We discuss in detail the differences between our method and previous methods, and demonstrate step-by-step how the signal-to-noise (S/N) ratio is increased. The method is based on detection of perturbations in acoustic phase travel times determined from cross-covariances of solar oscillations observed on the surface. We detect strong acoustic travel-time reductions of an order of 12 - 16 seconds at a depth of 42 - 75 Mm. These acoustic anomalies are detected 1 - 2 days before high peaks in the photospheric magnetic flux rate implying that the average emerging speed is 0.3 - 0.6 km/s. The results of this work contribute to our understanding of solar magnetism and benefit space weather forecasting.

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

We use 3D radiative MHD simulations to investigate the formation and dynamics of small-scale (less than 0.5 Mm in diameter) vortex tubes spontaneously generated by turbulent convection in quiet-Sun regions with initially weak mean magnetic fields. The results show that the vortex tubes penetrate into the chromosphere and substantially affect the structure and dynamics of the solar atmosphere. The vortex tubes are mostly concentrated in intergranular lanes and are characterized by strong (near sonic) downflows and swirling motions that capture and twist magnetic field lines, forming magnetic flux tubes that expand with height and which attain magnetic field strengths ranging from 200 G in the chromosphere to more than 1 kG in the photosphere. We investigate in detail the physical properties of these vortex tubes, including thermodynamic properties, flow dynamics, and kinetic and current helicities, and conclude that magnetized vortex tubes provide an important path for energy and momentum transfer from the convection zone into the chromosphere.

Investigation of the turbulent properties of solar convection is extremely important for understanding the multi-scale dynamics observed on the solar surface. In particular, recent high-resolution observations have revealed ubiquitous vortical structures, and numerical simulations have demonstrated links between vortex tube dynamics and magnetic field organization and have shown the importance of vortex tube interactions in the mechanisms of acoustic wave excitation on the Sun. In this paper we investigate the mechanisms of the formation of vortex tubes in highly-turbulent convective flows near the solar surface by using realistic radiative hydrodynamic LES simulations. Analysis of data from the simulations indicates two basic processes of vortex tube formation: 1) development of small-scale convective instability inside convective granules, and 2) a Kelvin-Helmholtz type instability of shearing flows in intergranular lanes. Our analysis shows that vortex stretching during these processes is a primary source of generation of small-scale vorticity on the Sun.

The magnetic flux that is generated by dynamo inside the Sun emerges in the form of bipolar magnetic regions. We have analyzed the whole set of solar magnetograms obtained with the SOHO/MDI instrument in 1995-2011, and automatically identified 160,079 bipolar magnetic regions that span a range of scale sizes across nearly four orders of magnitude. Their properties have been statistically analyzed, in particular with respect to the polarity orientations of the bipolar regions, including their tilt angle distributions. The latitude variation of the average tilt angles (with respect to the E-W direction), known as Joy's law, is found to closely follow the relation 32.1*sin(latitude)[deg]. There is no indication of a dependence on region size that one may expect if the tilts were produced by the Coriolis force during the buoyant rise of flux loops from the tachocline region. A few percent of all regions have orientations that violate Hale's polarity law. We show examples, from different phases of the solar cycle, where well defined medium-size bipolar regions with opposite polarity orientations occur side by side in the same latitude zone. Such oppositely oriented large bipolar regions cannot be part of the same toroidal flux system, but different flux systems must coexist in the same latitude zones. These examples are incompatible with the paradigm of coherent, subsurface toroidal flux ropes as the source of sunspots, and instead show that fluctuations must play a major role at all scales for the turbulent dynamo. We see no observational support for a separation of scales or a division between a global and a local dynamo, since also the smallest scales in the data set retain a non-random component that significantly contributes to the accumulated emergence of a N-S dipole moment that leads to the replacement of the old global poloidal field with a new one that has the opposite orientation.

Helioseismology can be used to place new constraints on the properties of dark matter, allowing solar observations to complement more conventional dark matter searches currently in operation. During the course of its lifetime, the Sun accretes a sizeable amount of dark matter. This accreted matter affects the heat transport of the stellar core in ways that helioseismology can detect, or at least constrain. We modify the CESAM stellar evolution code to take account of dark matter and determine the effect of WIMP models on the stellar structure and normal-mode oscillation frequencies.

We study the dynamical and statistical properties of turbulent cross-helicity (correlation of the aligned fluctuating velocity and magnetic field components). We derive an equation governing generation and evolution of the turbulent cross-helicity and discuss its meaning for the dynamo. Using symmetry properties of the problem we suggest a general expression for the turbulent cross-helicity pseudo-scalar and compute the turbulent coefficients in this expression. Effects of the density stratification, large-scale magnetic fields, differential rotation and turbulent convection are taken into account. We investigate the relative contribution of these effects to the cross-helicity evolution for two kinds of dynamo models of the solar cycle including a distributed mean-field model and a flux-transport dynamo model. We show that the contribution from the density stratification follows the evolution of the radial magnetic field, while large-scale electric currents produce a more complicated pattern of the cross-helicity of the comparable magnitude. We suggest that the results of observational analysis of the cross-helicity will depend on the averaging scales. Our results show that the pattern of the cross-helicity evolution strongly depends on details of the dynamo mechanism. Thus, we anticipate that direct observations of the cross-helicity on the Sun may serve for the diagnostic purpose of the solar dynamo process.

The paper presents a study of a solar dynamo model operating in the bulk of the convection zone with the toroidal magnetic field flux concentrated in the subsurface rotational shear layer. We explore how this type of dynamo may depend on spatial variations of turbulent parameters and on the differential rotation near the surface. The mean-field dynamo model takes into account the evolution of magnetic helicity and describes its nonlinear feedback on the generation of large-scale magnetic field by the $\alpha$-effect. We compare the magnetic cycle characteristics predicted by the model, including the cycle asymmetry (associated with the growth and decay times) and the duration - amplitude relation (Waldmeier's effects), with the observed sunspot cycle properties. We show that the model qualitatively reproduces the basic properties of the solar cycles.

Seismic traveltime measurements are a crucial tool in the investigation of the solar interior, particularly in the examination of fine-scale structure. Traditional analysis of traveltimes relies on a geometrical ray picture of acoustic wave propagation, which assumes high frequencies. However, it is well-known that traveltimes obtained from finite-frequency waves are sensitive to variations of medium parameters in a wide Fresnel zone around the ray path. To address this problem, Fr\'{e}chet traveltime sensitivity kernels have previously been developed. These kernels use a more realistic approximation of the wave propagation to obtain a linear relationship between traveltimes and variations in medium parameters. Fr\'{e}chet kernels take into account the actual frequency content of the measured waves and, thus, reproduce the Fresnel zone. Kernel theory has been well-developed in previous work on plane-parallel models of the Sun for use in local helioseismology. Our primary purpose is to apply kernel theory to much larger scales and in a spherical geometry. We also present kernel theory in a different way, using basic functional analytic methods, in the hope that this approach provides an even clearer understanding of the theory, as well as a set of tools for calculating kernels. Our results are very general and can be used to develop kernels for sensitivity to sound speed, density, magnetic fields, fluid flows, and any other medium parameter which can affect wave propagation.

The X2.2-class solar flare of February 15, 2011, produced a powerful sunquake event, representing a helioseismic response to the flare impact in the solar photosphere, which was observed with the HMI instrument on the Solar Dynamics Observatory (SDO). The impulsively excited acoustic waves formed a compact wavepacket traveling through the solar interior and appearing on the surface as expanding wave ripples. The initial flare impacts were observed in the form of compact and rapid variations of the Doppler velocity, line-of-sight magnetic field and continuum intensity. These variations formed a typical two-ribbon flare structure, and are believed to be associated with thermal and hydrodynamic effects of high-energy particles heating the lower atmosphere. The analysis of the SDO/HMI and X-ray data from the Ramaty High Energy Solar Spectroscopic Imager (RHESSI) shows that the helioseismic waves were initiated by the photospheric impact in the early impulsive phase, observed prior to the hard X-ray (50-100 keV) impulse, and were probably associated with atmospheric heating by relatively low-energy electrons (~6-50 keV) and heat flux transport. The impact caused a short motion in the sunspot penumbra prior to the appearance of the helioseismic wave. It is found that the helioseismic wave front traveling through a sunspot had a lower amplitude and was significantly delayed relative to the front traveling outside the spot. These observations open new perspectives for studying the flare photospheric impacts and for using the flare-excited waves for sunspot seismology.

The paper presents a study of kinematic axisymmetric mean-field dynamo models for a case of the meridional circulation with a deep-seated stagnation point and a strong return flow at the bottom of the convection zone. This kind of circulation follows from mean-field models of the angular momentum balance in the solar convection zone. We show that it is possible for this types of meridional circulation to construct kinematic dynamo models that resemble in some aspects the sunspot magnetic activity cycle. The dynamo model includes turbulent sources of the large-scale poloidal magnetic field production, due to kinetic helicity and a combined effect due to Coriolis force and the large-scale current. In these models the toroidal magnetic field, which is responsible for the sunspot production, is concentrated at the bottom of the convection zone, and is transported to low-latitude regions by the meridional flow. The meridional component of the poloidal field is also concentrated at the bottom of the convection zone while the radial component is concentrated in near polar regions. There are some issues which, perhaps, are resulted from the given meridional circulation pattern and the distribution of the magnetic diffusivity inside convection zone. In particular, in the near-equatorial regions the phase relations between the toroidal and poloidal components disagree with observations. Also, we show that the period of the magnetic cycle may not always monotonically decrease with the increase of the meridional flow speed. Thus, for the further progress it is important to determine the structure of the meridional circulation, which is one of the critical properties, from helioseismology observations.

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

Helioseismology studies the structure and dynamics of the Sun's interior by observing oscillations on the surface. These studies provide information about the physical processes that control the evolution and magnetic activity of the Sun. In recent years, helioseismology has made substantial progress towards the understanding of the physics of solar oscillations and the physical processes inside the Sun, thanks to observational, theoretical and modeling efforts. In addition to the global seismology of the Sun based on measurements of global oscillation modes, a new field of local helioseismology, which studies oscillation travel times and local frequency shifts, has been developed. It is capable of providing 3D images of the subsurface structures and flows. The basic principles, recent advances and perspectives of global and local helioseismology are reviewed in this article.

The X2.2-class solar flare of February 15, 2011, produced a powerful `sunquake' event, representing a seismic response to the flare impact. The impulsively excited seismic waves formed a compact wavepacket traveling through the solar interior and appeared on the surface as expanding wave ripples. The Helioseismic and Magnetic Imager (HMI), instrument on SDO, observes variations of intensity, magnetic field and plasma velocity (Dopplergrams) on the surface of Sun almost uninterruptedly with high resolution (0.5 arcsec/pixel) and high cadence (45 sec). The flare impact on the solar surface was observed in the form of compact and rapid variations of the HMI observables (Doppler velocity, line-of-sight magnetic field and continuum intensity). These variations, caused by the impact of high-energy particles in the photosphere, formed a typical two-ribbon flare structure. The sunquake can be easily seen in the raw Dopplergram differences without any special data processing. The source of this quake was located near the outer boundary of a very complicated complicated sunspot region, NOAA 1158, in a sunspot penumbra and at the penumbra boundary. This caused an interesting plasma dynamics in the impact region. I present some preliminary results of analysis of the near-real-time data from HMI, and discuss properties of the sunquake and the flare impact sources.

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

Rotation of the leading sunspot of active region NOAA 10421 was investigated using magnetograms and Dopplergrams from the MDI instrument of the Solar and Heliospheric Observatory (SOHO), and white-light images from the Transition Region and Coronal Explorer (TRACE). The vertical, radial, and azimuthal axisymmetrical components of both magnetic and velocity field vectors were reconstructed for the sunspot umbra and penumbra. All three components of both vectors in the umbra and penumbra show torsional oscillations with the same rotational period of about 3.8 days. The TRACE white-light data also show that the sunspot umbra and penumbra are torsionally rotating with the same period. Possible mechanisms of sunspot torsional motions are discussed.

Mechanisms of the formation and stability of sunspots are among the longest-standing and intriguing puzzles of solar physics and astrophysics. Sunspots are controlled by subsurface dynamics hidden from direct observations. Recently, substantial progress in our understanding of the physics of the turbulent magnetized plasma in strong-field regions has been made by using numerical simulations and local helioseismology. Both the simulations and helioseismic measurements are extremely challenging, but it becomes clear that the key to understanding the enigma of sunspots is a synergy between models and observations. Recent observations and radiative MHD numerical models have provided a convincing explanation to the Evershed flows in sunspot penumbrae. Also, they lead to the understanding of sunspots as self-organized magnetic structures in the turbulent plasma of the upper convection zone, which are maintained by a large-scale dynamics. Local helioseismic diagnostics of sunspots still have many uncertainties, some of which are discussed in this review. However, there have been significant achievements in resolving these uncertainties, verifying the basic results by new high-resolution observations, testing the helioseismic techniques by numerical simulations, and comparing results obtained by different methods. For instance, a recent analysis of helioseismology data from the Hinode space mission has successfully resolved several uncertainties and concerns (such as the inclined-field and phase-speed filtering effects) that might affect the inferences of the subsurface wave-speed structure of sunspots and the flow pattern. It becomes clear that for the understanding of the phenomenon of sunspots it is important to further improve the helioseismology methods and investigate the whole life cycle of active regions, from magnetic-flux emergence to dissipation.

In this study we investigate the effects of turbulent convection on formation of large-scale inhomogeneous magnetic structures by means of Large-Eddy Simulation (LES) for convection in solar-type stars. The main idea of this study is the implementation of a new subgrid-scale model for the effective Lorentz force in a three-dimensional nonlinear radiative magnetohydrodynamics (MHD) code developed for simulating the upper solar convection zone and lower atmosphere. To this end we derived the energy budget equations, which include the effects of the subgrid-scale turbulence on the Lorentz-force, and implemented the new subgrid-scale turbulence model (TELF-Model) in a three-dimensional nonlinear MHD LES code. Using imposed initial vertical and horizontal uniform magnetic fields in LES with the TELF-Model, we have shown that the magnetic flux tubes formation is started when the initial mean magnetic field is larger than a threshold value (about 100 G). This is in agreement with the theoretical studies by Rogachevskii and Kleeorin (2007). We have determined the vertical profiles of the velocity and magnetic fluctuations, total MHD energy and anisotropy of turbulent magneto-convection, kinetic and current and cross helicities.

We propose a solar dynamo model distributed in the bulk of the convection zone with the toroidal magnetic field the flux concentrated in the near-surface layer. We show that if the boundary conditions at the top of the dynamo region allow the large-scale toroidal magnetic fields to penetrate closer to the surface, then the pattern of the modeled butterfly diagram for the toroidal magnetic fields in the upper part of the convection zone is formed by the surface rotational shear layer. The model is in agreement with observed properties of the magnetic solar cycle.

We report on a new phenomenon of `alignment' of supergranulation cells in the polar regions of the Sun. Recent high-resolution datasets obtained by the Solar Optical Telescope onboard the Hinode satellite enabled us to investigate supergranular structures in high-latitude regions of the Sun. We have carried out a local helioseismology time-distance analysis of the data, and detected acoustic travel-time variations due to the supergranular flows. The supergranulation cells in both the north and south polar regions show systematic alignment patterns in the north-south direction. The south-pole datasets obtained in a month-long Hinode campaign indicate that the supergranulation alignment property may be quite common in the polar regions. We also discuss the latitudinal dependence of the supergranulation cell sizes; the data show that the east-west cell size decreases towards higher latitudes.

The sunspot penumbra is a transition zone between the strong vertical magnetic field area (sunspot umbra) and the quiet Sun. The penumbra has a fine filamentary structure that is characterized by magnetic field lines inclined toward the surface. Numerical simulations of solar convection in inclined magnetic field regions have provided an explanation of the filamentary structure and the Evershed outflow in the penumbra. In this paper, we use radiative MHD simulations to investigate the influence of the magnetic field inclination on the power spectrum of vertical velocity oscillations. The results reveal a strong shift of the resonance mode peaks to higher frequencies in the case of a highly inclined magnetic field. The frequency shift for the inclined field is significantly greater than that in vertical field regions of similar strength. This is consistent with the behavior of fast MHD waves.

We investigate how helioseismic waves that originate from effective point sources interact with a sunspot. These waves are reconstructed from observed stochastic wavefields on the Sun by cross-correlating photospheric Doppler-velocity signals. We select the wave sources at different locations relative to the sunspot, and investigate the p- and f-mode waves separately. The results reveal a complicated picture of waveform perturbations caused by the wave interaction with the sunspot. In particular, it is found that for waves originating from outside of the sunspot, p-mode waves travel across the sunspot with a small amplitude reduction and slightly higher speed, and wave amplitude and phase get mostly restored to the quiet-Sun values after passing the sunspot. The f-mode wave experiences some amplitude reduction passing through the sunspot, and the reduced amplitude is not recovered after that. The wave-propagation speed does not change before encountering the sunspot and inside the sunspot, but the wavefront becomes faster than the reference wave after passing through the sunspot. For waves originating from inside the sunspot umbra, both f- and p-mode waves show significant amplitude reductions and faster speed during all courses of propagation. A comparison of positive and negative time lags of cross-correlation functions shows an apparent asymmetry in the waveform changes for both the f- and p-mode waves. We suggest that the waveform variations of the helioseismic waves interacting with a sunspot found in this article can be used for developing a method of waveform heliotomography, similar to the waveform tomography of the Earth.

Five-minutes oscillations is one of the basic properties of solar convection. Observations show mixture of a large number of acoustic wave fronts propagating from their sources. We investigate the process of acoustic waves excitation from the point of view of individual events, by using realistic 3D radiative hydrodynamic simulation of the quiet Sun. The results show that the excitation events are related to dynamics vortex tubes (or swirls) in the intergranular lanes. These whirlpool-like flows are characterized by very strong horizontal velocities (7 - 11 km/s) and downflows (~ 7 km/s), and are accompanied by strong decreases of the temperature, density and pressure at the surface and in a ~ 0.5-1 Mm deep layer below the surface. High-speed whirlpool flows can attract and capture other vortices. According to our simulation results, the processes of the vortex interaction, such as vortex annihilation, can cause the excitation of acoustic waves.

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

Accurate subgrid-scale turbulence models are needed to perform realistic numerical magnetohydrodynamic (MHD) simulations of the subsurface flows of the Sun. To perform large-eddy simulations (LES) of turbulent MHD flows, three unknown terms have to be modeled. As a first step, this work proposes to use a priori tests to measure the accuracy of various models proposed to predict the SGS term appearing in the transport equation of the filtered magnetic field. It is proposed to evaluate the SGS model accuracy in term of "structural" and "functional" performance, i.e. the model capacity to locally approximate the unknown term and to reproduce its energetic action, respectively. From our tests, it appears that a mixed model based on the scale-similarity model has better performance.

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 results on modeling solar pores and sunspots using 2D axisymmetric magneto-hydrodynamic (MHD) simulations. These models are helpful for understanding the mechanisms of magnetic field concentration in sunspots, and the large-scale flow patterns associated with them. The simulations provide consistent, self-maintained, although not fully realistic, models of concentrated magnetic field near the solar surface. In this paper, we explore under which conditions the associated flows are converging or diverging near the surface. We find that in most cases in which a stable, pore-like concentration of magnetic field forms, a configuration with converging over diverging flow is established.

One of the puzzling features of solar magnetism is formation of long-living compact magnetic structures; such as sunspots and pores, in the highly turbulent upper layer of the solar convective zone. We use realistic radiative 3D MHD simulations to investigate the interaction between magnetic field and turbulent convection. In the simulations, a weak vertical uniform magnetic field is imposed in a region of fully developed granular convection; and the total magnetic flux through the top and bottom boundaries is kept constant. The simulation results reveal a process of spontaneous formation of stable magnetic structures, which may be a key to understanding of the magnetic self-organization on the Sun and formation of pores and sunspots. This process consists of two basic steps: 1) formation of small-scale filamentary magnetic structures associated with concentrations of vorticity and whirlpool-type motions, and 2) merging of these structures due to the vortex attraction, caused by converging downdrafts around magnetic concentration below the surface. In the resulting large-scale structure maintained by the converging plasma motions, the magnetic field strength reaches ~1.5 kG at the surface and ~6 kG in the interior; and the surface structure resembles solar pores. The magnetic structure remains stable for the whole simulation run of several hours with no sign of decay.

We show that under certain conditions, subsurface structures in the solar interior can alter the average acoustic power observed at the photosphere above them. By using numerical simulations of wave propagation, we show that this effect is large enough for it to be potentially used for detecting emerging active regions before they appear on the surface. In our simulations, simplified subsurface structures are modeled as regions with enhanced or reduced acoustic wave speed. We investigate the dependence of the acoustic power above a subsurface region on the sign, depth, and strength of the wave speed perturbation. Observations from the Solar and Heliospheric Observatory/Michelson Doppler Imager (SOHO/MDI) prior and during the emergence of NOAA active region 10488 are used to test the use of acoustic power as a potential precursor of magnetic flux emergence.

Recent spectro-polarimetric observations of a sunspot showed the formation of bipolar magnetic patches in the mid penumbra and their propagation toward the outer penumbral boundary. The observations were interpreted as being caused by sea-serpent magnetic fields near the solar surface (Sainz Dalda & Bellot Rubio 2008). In this Letter, we develop a 3D radiative MHD numerical model to explain the sea-serpent structure and the wave-like behavior of the penumbral magnetic field lines. The simulations reproduce the observed behavior, suggesting that the sea-serpent phenomenon is a consequence of magnetoconvection in a strongly inclined magnetic field. It involves several physical processes: filamentary structurization, high-speed overturning convective motions in strong, almost horizontal magnetic fields with partially frozen field lines, and traveling convective waves. The results demonstrate a correlation of the bipolar magnetic patches with high-speed Evershed downflows in the penumbra. This is the first time that a 3D numerical model of the penumbra results in downward directed magnetic fields, an essential ingredient of sunspot penumbrae that has eluded explanation until now.

We present comparison of numerical simulations of propagation of MHD waves,excited by subphotospheric perturbations, in two different ("deep" and "shallow") magnetostatic models of the sunspots. The "deep" sunspot model distorts both the shape of the wavefront and its amplitude stronger than the "shallow" model. For both sunspot models, the surface gravity waves (f-mode) are affected by the sunspots stronger than the acoustic p-modes. The wave amplitude inside the sunspot depends on the photospheric strength of the magnetic field and the distance of the source from the sunspot axis. For the source located at 9 Mm from the center of the sunspot, the wave amplitude increases when the wavefront passes through the central part of the sunspot. For the source distance of 12 Mm, the wave amplitude inside the sunspot is always smaller than outside. For the same source distance from the sunspot center but for the models with different strength of the magnetic field, the wave amplitude inside the sunspot increases with the strength of the magnetic field. The simulations show that unlike the case of the uniform inclined background magnetic field, the p- and f-mode waves are not spatially separated inside the sunspot where the magnetic field is strongly non-uniform. These properties have to be taken into account for interpretation of observations of MHD waves traveling through sunspot regions.

In recent years solar oscillations have been studied in great detail, both observationally and theoretically; so, perhaps, the Sun currently is the best understood pulsating star. The observational studies include long, almost uninterrupted series of oscillation data from the SOHO spacecraft and ground-based networks, GONG and BiSON, and more recently, extremely high-resolution observations from the Hinode mission. These observational data cover the whole oscillation spectrum, and have been extensively used for helioseismology studies, providing frequencies and travel times for diagnostics of the internal stratification, differential rotation, zonal and meridional flows, subsurface convection and sunspots. Together with realistic numerical simulations they lead to better understanding of the excitation mechanism and interactions of the oscillations with turbulence and magnetic fields. However, many problems remain unsolved. In particular, the precision of the helioseismology measurements is still insufficient for detecting the dynamo zone and deep routes of sunspots. Our knowledge of the oscillation physics in strong magnetic field regions is inadequate for interpretation of MHD waves in sunspots and for sunspot seismology. A new significant progress in studying the solar oscillations is expected from the Solar Dynamics Observatory scheduled for launch in 2010.

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.

Correct interpretation of acoustic travel times measured by time-distance helioseismology is essential to get an accurate understanding of the solar properties that are inferred from them. It has long been observed that sunspots suppress p-mode amplitude, but its implications on travel times has not been fully investigated so far. It has been found in test measurements using a 'masking' procedure, in which the solar Doppler signal in a localized quiet region of the Sun is artificially suppressed by a spatial function, and using numerical simulations that the amplitude modulations in combination with the phase-speed filtering may cause systematic shifts of acoustic travel times. To understand the properties of this procedure, we derive an analytical expression for the cross-covariance of a signal that has been modulated locally by a spatial function that has azimuthal symmetry, and then filtered by a phase speed filter typically used in time-distance helioseismology. Comparing this expression to the Gabor wavelet fitting formula without this effect, we find that there is a shift in the travel times, that is introduced by the amplitude modulation. The analytical model presented in this paper can be useful also for interpretation of travel time measurements for non-uniform distribution of oscillation amplitude due to observational effects.

We analyze a solar active region observed by the Hinode CaII H line using the time-distance helioseismology technique, and infer wave-speed perturbation structures and flow fields beneath the active region with a high spatial resolution. The general subsurface wave-speed structure is similar to the previous results obtained from SOHO/MDI observations. The general subsurface flow structure is also similar, and the downward flows beneath the sunspot and the mass circulations around the sunspot are clearly resolved. Below the sunspot, some organized divergent flow cells are observed, and these structures may indicate the existence of mesoscale convective motions. Near the light bridge inside the sunspot, hotter plasma is found beneath, and flows divergent from this area are observed. The Hinode data also allow us to investigate potential uncertainties caused by the use of phase-speed filter for short travel distances. Comparing the measurements with and without the phase-speed filtering, we find out that inside the sunspot, mean acoustic travel times are in basic agreement, but the values are underestimated by a factor of 20-40% inside the sunspot umbra for measurements with the filtering. The initial acoustic tomography results from Hinode show a great potential of using high-resolution observations for probing the internal structure and dynamics of sunspots.

Solar gravity modes (or g modes) -- oscillations of the solar interior for which buoyancy acts as the restoring force -- have the potential to provide unprecedented inference on the structure and dynamics of the solar core, inference that is not possible with the well observed acoustic modes (or p modes). The high amplitude of the g-mode eigenfunctions in the core and the evanesence of the modes in the convection zone make the modes particularly sensitive to the physical and dynamical conditions in the core. Owing to the existence of the convection zone, the g modes have very low amplitudes at photospheric levels, which makes the modes extremely hard to detect. In this paper, we review the current state of play regarding attempts to detect g modes. We review the theory of g modes, including theoretical estimation of the g-mode frequencies, amplitudes and damping rates. Then we go on to discuss the techniques that have been used to try to detect g modes. We review results in the literature, and finish by looking to the future, and the potential advances that can be made -- from both data and data-analysis perspectives -- to give unambiguous detections of individual g modes. The review ends by concluding that, at the time of writing, there is indeed a consensus amongst the authors that there is currently no undisputed detection of solar g modes.

The solar tachocline at the bottom of the convection zone is an important region for the dynamics of the Sun and the solar dynamo. In this region, the sound speed inferred by global helioseismology exhibits a bump of approximately 0.4% relative to the standard solar model. Global helioseismology does not provide any information on possible latitudinal variations or asymmetries between the Northern and Southern hemisphere. Here, we develop a time-distance helioseismology technique, including surface- and deep-focusing measurement schemes and a combination of both, for two-dimensional tomographic imaging of the solar tachocline that infers radial and latitudinal variations in the sound speed. We test the technique using artificial solar oscillation data obtained from numerical simulations. The technique successfully recovers major features of the simplified tachocline models. The technique is then applied to SOHO/MDI medium-l data and provides for the first time a full two-dimensional sound-speed perturbation image of the solar tachocline. The one-dimensional radial profile obtained by latitudinal averaging of the image is in good agreement with the previous global helioseismology result. It is found that the amplitude of the sound-speed perturbation at the tachocline varies with latitude, but it is not clear whether this is in part or fully an effect of instrumental distortion. Our initial results demonstrate that time-distance helioseismology can be used to probe the deep interior structure of the Sun, including the solar tachocline.

Discovered in 1909 the Evershed effect represents strong mass outflows in sunspot penumbra, where the magnetic field of sunspots is filamentary and almost horizontal. These flows play important role in sunspots and have been studied in detail using large ground-based and space telescopes, but the basic understanding of its mechanism is still missing. We present results of realistic numerical simulations of the Sun's subsurface dynamics, and argue that the key mechanism of this effect is in non-linear magnetoconvection that has properties of traveling waves in the presence of strong, highly inclined magnetic field. The simulations reproduce many observed features of the Evershed effect, including the high-speed "Evershed clouds", the filamentary structure of the flows, and the non-stationary quasi-periodic behavior. The results provide a synergy of previous theoretical models and lead to an interesting prediction of a large-scale organization of the outflows.

We report on a signature of chromospheric downflows in two emerging-flux regions detected by time-distance helioseismology analysis. We use both chromospheric intensity oscillation data in the Ca II H line and photospheric Dopplergrams in the Fe I 557.6nm line obtained by Hinode/SOT for our analyses. By cross-correlating the Ca II oscillation signals, we have detected a travel-time anomaly in the plage regions; outward travel times are shorter than inward travel times by 0.5-1 minute. However, such an anomaly is absent in the Fe I data. These results can be interpreted as evidence of downflows in the lower chromosphere. The downflow speed is estimated to be below 10 km/s. This result demonstrates a new possibility of studying chromospheric flows by time-distance analysis.

Variations of the sunspot number are important indicators of the solar activity cycles. The sunspot formation is a result of a dynamo process inside the Sun, which is far from being understood. We use simple dynamical models of the dynamo process to simulate the magnetic field evolution and investigate general properties of the sunspot number variations during the solar cycles. We have found that the classical Parker's model with a standard kinetic helicity quenching cannot represent the typical profiles of the solar-cycle variations of the sunspot number, and also does not give chaotic solutions. For modeling of the solar cycle properties we use a nonlinear dynamo model of Kleeorin and Ruzmaikin (1982), which takes into account dynamics of the turbulent magnetic helicity. We have obtained a series of periodic and chaotic solutions for different layers of the convective zone. The solutions qualitatively reproduce some basic observational features of the solar cycle properties, in particular, the relationship between the growth time and the cycle amplitude. Also, on the longer time scale the dynamo model with the magnetic helicity has intermittent solutions, which may be important for modeling long-term variations of the solar cycles.

Despite the known general properties of the solar cycles, a reliable forecast of the 11-year sunspot number variations is still a problem. The difficulties are caused by the apparent chaotic behavior of the sunspot numbers from cycle to cycle and by the influence of various turbulent dynamo processes, which are far from understanding. For predicting the solar cycle properties we make an initial attempt to use the Ensemble Kalman Filter (EnKF), a data assimilation method, which takes into account uncertainties of a dynamo model and measurements, and allows to estimate future observational data. We present the results of forecasting of the solar cycles obtained by the EnKF method in application to a low-mode nonlinear dynamical system modeling the solar $\alpha\Omega$-dynamo process with variable magnetic helicity. Calculations of the predictions for the previous sunspot cycles show a reasonable agreement with the actual data. This forecast model predicts that the next sunspot cycle will be significantly weaker (by $\sim 30%$) than the previous cycle, continuing the trend of low solar activity.

We have developed 3D, compressible, non-linear radiative MHD simulations to study the influence of the magnetic field of various strength and geometry on the turbulent convective cells and on the excitation mechanisms of the acoustic oscillations. The results reveal substantial changes of the granulation structure with increased magnetic field, and a frequency-dependent reduction in the oscillation power. These simulation results reproduce the enhanced high-frequency acoustic emission observed at the boundaries of active region ("acoustic halo" phenomenon). In the presence of inclined magnetic field the solar convection develops filamentary structure with flows concentrated along magnetic filaments, and also exhibits behavior of running magnetoconvective waves, resembling recent observations of the sunspot penumbra dynamics from Hinode/SOT.

The complex dynamics that lead to the emergence of active regions on the sun are poorly understood. One possibility is that magnetic structures (flux tubes, etc.) rise from below the surface by self induction and convection that lead to the formation of active regions and sunspots on the solar surface. For space weather forecasting, one would like to detect the subsurface structures before they reach the surface. The goal of this study is to investigate whether sound speed perturbations associated with subsurface structures could affect the acoustic power observed at the solar surface above them. Possible mechanisms for this effect are wave reflection, scattering or diffraction. By using numerical simulations of wave propagation in the solar interior, we investigate whether observations of the acoustic power can be used to detect emerging active regions before they appear on the surface. In the simulations, subsurface structures are modeled as regions with enhanced or reduced acoustic wavespeed. We show how the acoustic power above a subsurface region depends on the sign, depth and strength of the wavespeed perturbation. For comparison, we analyze observations from SOHO/MDI of the emergence of solar active region NOAA 10488.

Magnetic fields emerging from the Sun's interior carry information about physical processes of magnetic field generation and transport in the convection zone. Soon after appearance on the solar surface the magnetic flux gets concentrated in sunspot regions and causes numerous active phenomena on the Sun. This paper discusses some properties of the emerging magnetic flux observed on the solar surface and in the interior. A statistical analysis of variations of the tilt angle of bipolar magnetic regions during the emergence shows that the systematic tilt with respect to the equator (the Joy's law) is most likely established below the surface. However, no evidence of the dependence of the tilt angle on the amount of emerging magnetic flux, predicted by the rising magnetic flux rope theories, is found. Analysis of surface plasma flows in a large emerging active region reveals strong localized upflows and downflows at the initial phase of emergence but finds no evidence for large-scale flows indicating future appearance a large-scale magnetic structure. Local helioseismology provides important tools for mapping perturbations of the wave speed and mass flows below the surface. Initial results from SOHO/MDI and GONG reveal strong diverging flows during the flux emergence, and also localized converging flows around stable sunspots. The wave speed images obtained during the process of formation of a large active region, NOAA 10488, indicate that the magnetic flux gets concentrated in strong field structures just below the surface. Further studies of magnetic flux emergence require systematic helioseismic observations from the ground and space, and realistic MHD simulations of the subsurface dynamics.

In the recent papers, we introduced a method utilised to measure the flow field. The method is based on the tracking of supergranular structures. We did not precisely know, whether its results represent the flow field in the photosphere or in some sub-photospheric layers. In this paper, in combination with helioseismic data, we are able to estimate the depths in the solar convection envelope, where the detected large-scale flow field is well represented by the surface measurements. We got a clear answer to question what kind of structures we track in full-disc Dopplergrams. It seems that in the quiet Sun regions the supergranular structures are tracked, while in the regions with the magnetic field the structures of the magnetic field are dominant. This observation seems obvious, because the nature of Doppler structures is different in the magnetic regions and in the quiet Sun. We show that the large-scale flow detected by our method represents the motion of plasma in layers down to ~10 Mm. The supergranules may therefore be treated as the objects carried by the underlying large-scale velocity field.

We use 2D numerical simulations and eikonal approximation, to study properties of MHD waves traveling below the solar surface through the magnetic structure of sunspots. We consider a series of magnetostatic models of sunspots of different magnetic field strengths, from 10 Mm below the photosphere to the low chromosphere. The purpose of these studies is to quantify the effect of the magnetic field on local helioseismology measurements by modeling waves excited by sub-photospheric sources. Time-distance propagation diagrams and wave travel times are calculated for models of various field strength and compared to the non-magnetic case. The results clearly indicate that the observed time-distance helioseismology signals in sunspot regions correspond to fast MHD waves. The slow MHD waves form a distinctly different pattern in the time-distance diagram, which has not been detected in observations. The numerical results are in good agreement with the solution in the short-wavelength (eikonal) approximation, providing its validation. The frequency dependence of the travel times is in a good qualitative agreement with observations.

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

Investigation of propagation, conversion, and scattering of MHD waves in the Sun is very important for understanding the mechanisms of observed oscillations and waves in sunspots and active regions. We have developed 3D linear MHD numerical model to investigate influence of the magnetic field on excitation and properties of the MHD waves. The results show that the magnetic field can substantially change the properties of the surface gravity waves (f-mode), but their influence on the acoustic-type waves (p-modes) is rather moderate. Comparison our simulations with the time-distance helioseismology results from SOHO/MDI shows that the travel time variations caused by the inclined magnetic field do not exceed 25% of the observed amplitude even for strong fields of 1400-1900 G. This can be an indication that other effects (e.g. background flows and non-uniform distribution of magnetic field) can contribute to the observed travel time variations. The travel time variations caused by the wave interaction with magnetic field are in phase with the observations for strong fields of 1400-1900 G if Doppler velocities are taken at the height of 300 km above the photosphere where plasma parameter beta<<1. The simulations show that the travel times only weakly depend on the height of velocity observation. For the photospheric level the travel times are systematically smaller on approximately 0.12 min then for the hight of 300 km above the photosphere for all studied ranges of the magnetic field strength and inclination angles. The numerical MHD wave modeling and new data from the HMI instrument of the Solar Dynamics Observatory will substantially advance our knowledge of the wave interaction with strong magnetic fields on the Sun and improve the local helioseismology diagnostics.

The goal of this research is to investigate how magnetic field affects the dynamics of granular convection and excitation of solar oscillations by means of realistic numerical simulations. We have used a 3D, compressible, non-linear radiative magnetohydrodynamics code developed at the NASA Ames Research Center. This code takes into account several physical phenomena: compressible fluid flow in a highly stratified medium, sub-grid scale turbulence models, radiative energy transfer between the fluid elements, and a real-gas equation of state. We have studied the influence of the magnetic field of various strength on the convective cells and on the excitation mechanisms of the acoustic oscillations by calculating spectral properties of the convective motions and oscillations. The results reveal substantial changes of the granulation structure with increased magnetic field, and a frequency-dependent reduction in the oscillation power in a good agreement with solar observations. These simulations suggest that the enhanced high-frequency acoustic emission at the boundaries of active region ("acoustic halo" phenomenon) is caused by the changes of the spatial-temporal spectrum of the turbulent convection in magnetic field, resulting in turbulent motions of smaller scales and higher frequencies than in quiet Sun regions.

The aim of this paper is to extend our previous study of the solar-cycle variations of the meridional flows and to investigate their latitudinal and longitudinal structure in the subphotospheric layer, especially their variations in magnetic regions. Helioseismology observations indicate that mass flows around active regions are dominated by inflows into those regions. On average, those local flows are more important around leading magnetic polarities of active regions than around the following polarities, and depend on the evolutionary stage of particular active regions. We present a statistical study based on MDI/SOHO observations of 1996-2002 and show that this effect explains a significant part of the cyclic change of meridional flows in near-equatorial regions, but not at higher latitudes. A different mechanism driving solar-cycle variations of the meridional flow probably operates.

Far-side images of solar active regions have become one of the routine products of helioseismic observations, and are of importance for space weather forecasting by allowing the detection of sunspot regions before they become visible on the Earth side of the Sun. An accurate assessment of the quality of the far-side maps is difficult, because there are no direct observations of the solar far side to verify the detections. In this paper we assess far-side imaging based on the time-distance helioseismology method, by using numerical simulations of solar oscillations in a spherical solar model. Localized variations in the speed of sound in the surface and subsurface layers are used to model the perturbations associated with sunspots and active regions. We examine how the accuracy of the resulting far-side maps of acoustic travel times depends on the size and location of active regions. We investigate potential artifacts in the far-side imaging procedure, such as those caused by the presence of active regions on the solar near side, and suggest how these artifacts can be identified in the real Sun far-side images obtained from SOHO/MDI and GONG data.

By analyzing numerically simulated solar oscillation data, we study the influence of non-uniform distribution of acoustic wave amplitude, acoustic source strength, and perturbations of the sound speed on the shifts of acoustic travel times measured by the time-distance helioseismology method. It is found that for short distances, the contribution to the mean travel time shift caused by non-uniform distribution of acoustic sources in sunspots may be comparable to (but smaller than) the contribution from the sound speed perturbation in sunspots, and that it has the opposite sign to the sound-speed effect. This effect may cause some underestimation of the negative sound-speed perturbations in sunspots just below the surface, that was found in previous time-distance helioseismology inferences. This effect cannot be corrected by artificially increasing the amplitude of oscillations in sunspots. For large time-distance annuli, the non-uniform distribution of wavefields does not have significant effects on the mean travel times, and thus the sound-speed inversion results. The measured travel time differences, which are used to determine the mass flows beneath sunspots, can also be systematically shifted by this effect, but only in an insignificant magnitude.

Aims. With the Solar Optical Telescope on Hinode, we investigate the basic properties of high-degree solar oscillations observed at two levels in the solar atmosphere, in the G-band (formed in the photosphere) and in the Ca II H line (chromospheric emission). Methods. We analyzed the data by calculating the individual power spectra as well as the cross-spectral properties, i.e., coherence and phase shift. The observational properties are compared with a simple theoretical model, which includes the effects of correlated noise. Results. The results reveal significant frequency shifts between the Ca II H and G-band spectra, in particular above the acoustic cut-off frequency for pseudo-modes. The cross-spectrum phase shows peaks associated with the acoustic oscillation (p-mode) lines, and begins to increase with frequency around the acoustic cut-off. However, we find no phase shift for the (surface gravity wave) f-mode. The observed properties for the p-modes are qualitatively reproduced in a simple model with a correlated background if the correlated noise level in the Ca II H data is higher than in the G-band data. These results suggest that multi-wavelength observations of solar oscillations, in combination with the traditional intensity-velocity observations, may help to determine the level of the correlated background noise and to determine the type of wave excitation sources on the Sun.

Two methods are used to observationally determine the solar radius: One is the observation of the intensity profile at the limb, the other one uses f-mode frequencies to derive a 'seismic' solar radius which is then corrected to optical depth unity. The two methods are inconsistent and lead to a difference in the solar radius of $\sim$0.3 Mm. Because of the geometrical extention of the solar photosphere and the increased path lengths of tangential rays the Sun appears to be larger to an observer who measures the extent of the solar disk. Based on radiative transfer calculations we show that this discrepancy can be explained by the difference between the height at disk center where $\tau_{\mathrm{5000}}=1$ ($\tau_{\mathrm{Ross}}=2/3$) and the inflection point of the intensity profile on the limb. We calculate the intensity profile of the limb for the MDI continuum and the continuum at 5000 {\AA} for two atmosphere structures and compare the position of the inflection points with the radius at $\tau_{\mathrm{5000}}=1$ ($\tau_{\mathrm{Ross}}=2/3$). The calculated difference between the 'seismic' radius and the inflection point is $0.347\pm 0.06$ Mm with respect to $\tau_{\mathrm{5000}}=1$ and $0.333\pm 0.08$ Mm with respect to $\tau_{\mathrm{Ross}}=2/3$. We conclude that the standard solar radius in evolutionary models has to be lowered by $0.333\pm 0.08$ Mm and is 695.66 Mm. Furthermore, this correction reconciles inflection point measurements and the 'seismic' radii within the uncertainty.

The goal of this research is to investigate how well various turbulence models can describe physical properties of the upper convective boundary layer of the Sun. An accurate modeling of the turbulence motions is necessary for understanding the excitation mechanisms of solar oscillation modes. We have carried out realistic numerical simulations using several different physical Large Eddy Simulation (LES) models (Hyperviscosity approach, Smagorinsky, and dynamic models) to investigate how the differences in turbulence modeling affect the damping and excitation of the oscillations and their spectral properties and compare with observations. We have first calculated the oscillation power spectra of radial and non-radial modes supported by the computational box with the different turbulence models. Then we have calculated the work integral input to the modes to estimate the influence of the turbulence model on the depth and strength of the oscillation sources. We have compared these results with previous studies and with the observed properties of solar oscillations. We find that the dynamic turbulence model provides the best agreement with the helioseismic observations.

We use the magnetic butterfly diagram to determine the speed of the magnetic flux transport on the solar surface towards the poles. The manifestation of the flux transport is clearly visible as elongated structures extended from the sunspot belt to the polar regions. The slopes of these structures are measured and interpreted as meridional magnetic flux transport speed. Comparison with the time-distance helioseismology measurements of the mean speed of the meridional flows at the depth of 3.5--12 Mm shows a generally good agreement, but the speeds of the flux transport and the meridional flow are significantly different in areas occupied by the magnetic field. The local circulation flows around active regions, especially the strong equatorward flows on the equatorial side of active regions affect the mean velocity profile derived by helioseismology, but do not influence the magnetic flux transport. The results show that the mean longitudinally averaged meridional flow measurements by helioseismology may not be used directly in solar dynamo models for describing the magnetic flux transport, and that it is necessary to take into account the longitudinal structure of these flows.

Observations of a large solar flare of December 13, 2006, using Solar Optical Telescope (SOT) on Hinode spacecraft revealed high-frequency oscillations excited by the flare in the sunspot chromosphere. These oscillations are observed in the region of strong magnetic field of the sunspot umbra, and may provide a new diagnostic tool for probing the structure of sunspots and understanding physical processes in solar flares.

Analysis of the hydrodynamic and helioseismic effects in the photosphere during the solar flare of July 23, 2002, observed by Michelson Doppler Imager (MDI) on SOHO, and high-energy images from RHESSI shows that these effects are closely associated with sources of the hard X-ray emission, and that there are no such effects in the centroid region of the flare gamma-ray emission. These results demonstrate that contrary to expectations the hydrodynamic and helioseismic responses (''sunquakes") are more likely to be caused by accelerated electrons than by high-energy protons. A series of multiple impulses of high-energy electrons forms a hydrodynamic source moving in the photosphere with a supersonic speed. The moving source plays a critical role in the formation of the anisotropic wave front of sunquakes.

Results from initial helioseismic observations by Solar Optical Telescope onboard Hinode are reported. It has been demonstrated that intensity oscillation data from Broadband Filter Imager can be used for various helioseismic analyses. The k-omega power spectra, as well as corresponding time-distance cross-correlation function that promises high-resolution time-distance analysis below 6-Mm travelling distance, were obtained for G-band and CaII-H data. Subsurface supergranular patterns have been observed from our first time-distance analysis. The results show that the solar oscillation spectrum is extended to much higher frequencies and wavenumbers, and the time-distance diagram is extended to much shorter travel distances and times than they were observed before, thus revealing great potential for high-resolution helioseismic observations from Hinode.

Exploiting high-resolution observations made by the Solar Optical Telescope onboard Hinode, we investigate the spatial distribution of power spectral density of oscillatory signal in and around NOAA active region 10935. The G-band data show that in the umbra the oscillatory power is suppressed in all frequency ranges. On the other hand, in Ca II H intensity maps oscillations in the umbra, so-called umbral flashes, are clearly seen with the power peaking around 5.5 mHz. The Ca II H power distribution shows the enhanced elements with the spatial scale of the umbral flashes over most of the umbra but there is a region with suppressed power at the center of the umbra. The origin and property of this node-like feature remain unexplained.

Five-minute oscillations on the Sun (acoustic and surface gravity waves) are excited by subsurface turbulent convection. However, in sunspots the excitation is suppressed because strong magnetic field inhibits convection. We use 3D simulations to investigate how the suppression of excitation sources affects the distribution of the oscillation power in sunspot regions. The amplitude of random acoustic sources was reduced in circular-shaped regions to simulate the suppression in sunspots. The simulation results show that the amplitude of the oscillations can be approximately 2-4 times lower in the sunspot regions in comparison to the quiet Sun, just because of the suppressed sources. Using SOHO/MDI data we measured the amplitude ratio for the same frequency bands outside and inside sunspots, and found that this ratio is approximately 3-4. Hence, the absence of excitation sources inside sunspots makes a significant contribution (about 50% or higher) to the observed amplitude ratio and must be taken into account in sunspot seismology.

In time-distance helioseismology, the time signals (Doppler shifts) at two points on the solar surface, separated by a fixed angular distance are cross-correlated, and this leads to a wave packet signal. Accurately measuring the travel times of these wave packets is crucial for inferring the sub-surface properties in the Sun. The observed signal is quite noisy, and to improve the signal-to-noise ratio and make the cross-correlation more robust, the temporal oscillation signal is phase-speed filtered at the two points in order to select waves that travel a fixed horizontal distance. Hence a new formula to estimate the travel times is derived, in the presence of a phase speed filter, and it includes both the radial and horizontal component of the oscillation displacement signal. It generalizes the previously used Gabor wavelet that was derived without a phase speed filter and included only the radial component of the displacement. This is important since it will be consistent with the observed cross-correlation that is computed using a phase speed filter, and also it accounts for both the components of the displacement. The new formula depends on the location of the two points on the solar surface that are being cross correlated and accounts for the travel time shifts at different locations on the solar surface.

Recent models of variations of the Sun's structure with the 11-year activity cycle by Sofia et al (2005) predict strong non-homologous changes of the radius of subsurface layers, due to subsurface magnetic fields and field-modulated turbulence. According to their best model the changes of the surface radius may be 1000 times larger than those at the depth of 5 Mm. We use f-mode oscillation frequency data from the MDI instrument of Solar and Heliospheric Observatory (SOHO) and measurements of the solar surface radius variations from SOHO and ground-based observatories during solar cycle 23 (1996-2005) to put constraints on the radius changes. The results show that the above model overestimates the change of the radius at the surface relative to the change at 5Mm.

We apply time-distance helioseismology, local correlation tracking and Fourier spatial-temporal filtering methods to realistic supergranule scale simulations of solar convection and compare the results with high-resolution observations from the SOHO Michelson Doppler Imager (MDI). Our objective is to investigate the surface and sub-surface convective structures and test helioseismic measurements. The size and grid of the computational domain are sufficient to resolve various convective scales from granulation to supergranulation. The spatial velocity spectrum is approximately a power law for scales larger than granules, with a continuous decrease in velocity amplitude with increasing size. Aside from granulation no special scales exist, although a small enhancement in power at supergranulation scales can be seen. We calculate the time-distance diagram for f- and p-modes and show that it is consistent with the SOHO/MDI observations. From the simulation data we calculate travel time maps for surface gravity waves (f-mode). We also apply correlation tracking to the simulated vertical velocity in the photosphere to calculate the corresponding horizontal flows. We compare both of these to the actual large-scale (filtered) simulation velocities. All three methods reveal similar large scale convective patterns and provide an initial test of time-distance methods.

We present a direct comparison between two different techniques time-distance helioseismology and a local correlation tracking method for measuring mass flows in the solar photosphere and in a near-surface layer: We applied both methods to the same dataset (MDI high-cadence Dopplergrams covering almost the entire Carrington rotation 1974) and compared the results. We found that after necessary corrections, the vector flow fields obtained by these techniques are very similar. The median difference between directions of corresponding vectors is 24 degrees, and the correlation coefficients of the results for mean zonal and meridional flows are 0.98 and 0.88 respectively. The largest discrepancies are found in areas of small velocities where the inaccuracies of the computed vectors play a significant role. The good agreement of these two methods increases confidence in the reliability of large-scale synoptic maps obtained by them.

Recent progress in realistic simulations of solar convection have given us an unprecedented opportunity to evaluate the robustness of solar interior structures and dynamics obtained by methods of local helioseismology. We present results of testing the time-distance method using realistic simulations. By computing acoustic wave propagation time and distance relations for different depths of the simulated data, we confirm that acoustic waves propagate into the interior and then turn back to the photosphere. This demonstrates that in the numerical simulations properties of acoustic waves (p-modes) are similar to the solar conditions, and that these properties can be analyzed by the time-distance technique. For the surface gravity waves (f-mode), we calculate perturbations of their travel times, caused by localized downdrafts, and demonstrate that the spatial pattern of these perturbations (representing so-called sensitivity kernels) is similar to the patterns obtained from the real Sun, displaying characteristic hyperbolic structures. We then test the time-distance measurements and inversions by calculating acoustic travel times from a sequence of vertical velocities at the photosphere of the simulated data, and inferring a mean 3D flow fields by performing inversion based on the ray approximation. The inverted horizontal flow fields agree very well with the simulated data in subsurface areas up to 3 Mm deep, but differ in deeper areas. Due to the cross-talk effects between the horizontal divergence and downward flows, the inverted vertical velocities are significantly different from the mean convection velocities of the simulation dataset.

Results of numerical 3D simulations of propagation of acoustic waves inside the Sun are presented. A linear 3D code which utilizes realistic OPAL equation of state was developed by authors. Modified convectively stable standard solar model with smoothly joined chromosphere was used as a background model. High order dispersion relation preserving numerical scheme was used to calculate spatial derivatives. The top non-reflecting boundary condition established in the chromosphere absorbs waves with frequencies greater than the acoustic cut-off frequency which pass to the chromosphere, simulating a realistic situation. The acoustic power spectra obtained from the wave field generated by sources randomly distributed below the photosphere are in good agreement with observations. The influence of the height of the top boundary on results of simulation was studied. It was shown that the energy leakage through the acoustic potential barrier damps all modes uniformly and does not change the shape of the acoustic spectrum. So the height of the top boundary can be used for controlling a damping rate without distortion of the acoustic spectrum. The developed simulations provide an important tool for testing local helioseismology.

Recent observations of helioseismic holography revealed that magnetic fields inclined to the line-of-sight direction could cause systematic variations in measured acoustic phase shifts, and that the surface magnetic field may shift the phases and impair the coherence of acoustic waves (known as "showerglass effect"). We examine how these affect time-distance helioseismology measurements. It is confirmed that the inclined magnetic field could cause variations in acoustic travel times inside sunspot penumbra, however, inversions of the measured times show that this effect only slightly shifts the location of negative sound-speed variations near the solar surface, but does not change the inverted deeper interior structures. Regarding to the showerglass effect, we find that outgoing and incoming travel time perturbations through sunspots are significantly smaller than those reported from helioseismic holography. In addition, our second-skip cross-correlation experiments demonstrate that inside sunspots, the half of the double-skip travel times are very similar to the mean single-skip travel times. We conclude that these surface magnetism effects do not cause considerable systematic errors in time-distance helioseismology of active regions.

The solar seismic waves excited by solar flares (``sunquakes'') are observed as circular expanding waves on the Sun's surface. The first sunquake was observed for a flare of July 9, 1996, from the Solar and Heliospheric Observatory (SOHO) space mission. However, when the new solar cycle started in 1997, the observations of solar flares from SOHO did not show the seismic waves, similar to the 1996 event, even for large X-class flares during the solar maximum in 2000-2002. The first evidence of the seismic flare signal in this solar cycle was obtained for the 2003 ``Halloween'' events, through acoustic ``egression power'' by Donea and Lindsey. After these several other strong sunquakes have been observed. Here, I present a detailed analysis of the basic properties of the helioseismic waves generated by three solar flares in 2003-2005. For two of these flares, X17 flare of October 28, 2003, and X1.2 flare of January 15, 2005, the helioseismology observations are compared with simultaneous observations of flare X-ray fluxes measured from the RHESSI satellite. These observations show a close association between the flare seismic waves and the hard X-ray source, indicating that high-energy electrons accelerated during the flare impulsive phase produced strong compression waves in the photosphere, causing the sunquake. The results also reveal new physical properties such as strong anisotropy of the seismic waves, the amplitude of which varies significantly with the direction of propagation. The waves travel through surrounding sunspot regions to large distances, up to 120 Mm, without significant decay. These observations open new perspectives for helioseismic diagnostics of flaring active regions on the Sun and for understanding the mechanisms of the energy release and transport in solar flares.

We report on the changes of the Sun's subsurface stratification inferred from helioseismology data. Using SOHO/MDI (SOlar and Heliospheric Observatory/Michelson Doppler Imager) data for the last 9 years and, more precisely, the temporal variation of f-mode frequencies, we have computed the variation of the radius of subsurface layers of the Sun by applying helioseismic inversions. We have found a variability of the ``helioseismic'' radius in antiphase with the solar activity, with the strongest variations of the stratification being just below the surface around 0.995$R_{\odot}$. Besides, the radius of the deeper layers of the Sun, between 0.975$R_{\odot}$ and 0.99$R_{\odot}$ changes in phase with the 11-year cycle.

The SoHO spacecraft made astrophysicists achieve a major breakthrough in the knowledge of the Sun. In helioseismology, both GOLF and MDI experiments aboard SoHO greatly improve the accuracy of seismic data. More specifically, the detection of an enhanced number of low degree low order modes improves the accuracy on the sound speed and density profiles in the solar core. After a description of these profiles, we build solar seismic models. Different models are considered and enable us to derive precise emitted neutrino fluxes. These ones are validated by the seismic data and are in agreement with the recent detected neutrinos, assuming 3 neutrino flavors. The seismic models are also used to put limits on large scale magnetic fields in the solar interior. This analysis puts some upper bounds of about 3x10^7 G in the radiative zone. Such a field could slightly improve the emitted neutrino flux, which remains in agreement with the Sudbury Neutrino Observatory result of 2001. From the models we deduce gravity mode predictions, and the electron and neutron radial densities that are useful to calculate the neutrino oscillations. We also begin to discuss how the external magnetic field may influence such quantities.

We suggest a solution to an important problem of observational helioseismology of the separation of lines of solar acoustic (p) modes of low angular degree in oscillation power spectra by constructing optimal masks for Doppler images of the Sun. Accurate measurements of oscillation frequencies of low-degree modes are essential for the determination of the structure and rotation of the solar core. However, these measurements for a particular mode are often affected by leakage of other p modes arising when the Doppler images are projected on to spherical-harmonics masks. The leakage results in overlaping peaks corresponding to different oscillation modes in the power spectra. In this paper we present a method for calculating optimal masks for a given (target) mode by minimizing the signals of other modes appearing in its vicinity. We apply this method to time series of 2 years obtained from Michelson Doppler Imager (MDI) instrument on board SOHO space mission and demonstrate its ability to reduce efficiently the mode leakage.

We use the three-dimensional hydrodynamic code of Stein and Nordlund to realistically simulate the upper layers of the solar convection zone in order to study physical characteristics of solar oscillations. Our first result is that the properties of oscillation modes in the simulation closely match the observed properties. Recent observations from SOHO/MDI and GONG have confirmed the asymmetry of solar oscillation line profiles, initially discovered by Duvall et al. In this paper we compare the line profiles in the power spectra of the Doppler velocity and continuum intensity oscillations from the SOHO/MDI observations with the simulation. We also compare the phase differences between the velocity and intensity data. We have found that the simulated line profiles are asymmetric and have the same asymmetry reversal between velocity and intensity as observed. The phase difference between the velocity and intensity signals is negative at low frequencies and jumps in the vicinity of modes as is also observed. Thus, our numerical model reproduces the basic observed properties of solar oscillations, and allows us to study the physical properties which are not observed.

The 11-year cycle of solar activity follows Hale's law by reversing the magnetic polarity of leading and following sunspots in bipolar regions during the minima of activity. In the 1996-97 solar minimum, most solar activity emerged in narrow longitudinal zones - `active longitudes' but over a range in latitude. Investigating the distribution of solar magnetic flux, we have found that the Hale sunspot polarity reversal first occurred in these active zones. We have estimated the rotation rates of the magnetic flux in the active zones before and after the polarity reversal. Comparing these rotation rates with the internal rotation inferred by helioseismology, we suggest that both `old' and `new' magnetic fluxes were probably generated in a low-latitude zone near the base of the solar convection zone. The reversal of active region polarity observed in certain longitudes at the beginning of a new solar cycle suggests that the phenomenon of active longitudes may give fundamental information about the mechanism of the solar cycle. The non-random distribution of old-cycle and new-cycle fluxes presents a challenge for dynamo theories, most of which assume a uniform longitudinal distribution of solar magnetic fields.

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.