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.

Prediction of solar activity cycles is challenging because physical processes inside the Sun involve a broad range of multiscale dynamics that no model can reproduce and because the available observations are highly limited and cover mostly surface layers. Helioseismology makes it possible to probe solar dynamics in the convective zone, but variations in differential rotation and meridional circulation are currently available for only two solar activity cycles. It has been demonstrated that sunspot observations, which cover over 400 years, can be used to calibrate the Parker-Kleeorin-Ruzmaikin dynamo model, and that the Ensemble Kalman Filter (EnKF) method can be used to link the modeled magnetic fields to sunspot observations and make reliable predictions of a following activity cycle. However, for more accurate predictions, it is necessary to use actual observations of the solar magnetic fields, which are available only for the last four solar cycles. In this paper I briefly discuss the influence of the limited number of available observations on the accuracy of EnKF estimates of solar cycle parameters, the criteria to evaluate the predictions, and application of synoptic magnetograms to the prediction of solar activity.

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.

Building a reliable forecast of solar activity is a long-standing problem that requires to accurately describe past and current global dynamics. However, synoptic observations of magnetic fields and subsurface flows became available relatively recently. In this paper, we present an investigation of effects of short observational data series on accuracy of solar cycle prediction. This analysis is performed using the annual sunspot number time-series applied to the Parker-Kleeorin-Ruzmaikin dynamo model and employing the Ensemble Kalman Filter (EnKF) data assimilation method. The testing of the cycle prediction accuracy is performed for the last six cycles (from Solar Cycle 19 to 24) by sequentially shortening the observational data series that are used for prediction of a target cycle, and evaluating the prediction accuracy according to specified criteria. According to the analysis, reliable activity predictions can be made using relatively short time-series of the sunspot number. It demonstrated that even two cycles of available observations allow us to obtain reasonable forecasts.

Synoptic magnetograms provide us with knowledge about the evolution of magnetic fields on the solar surface and present important information for forecasting future solar activity. In this work, poloidal and toroidal magnetic field components derived from synoptic magnetograms are assimilated, using the Ensemble Kalman Filter method, into a mean-field dynamo model based on Parker's migratory dynamo theory complemented by magnetic helicity conservation. It was found that the predicted toroidal field is in good agreement with observations for almost the entire following solar cycle. However, poloidal field predictions agree with observations only for the first 2 - 3 years of the predicted cycle. The results indicate that the upcoming Solar Maximum of Cycle 25 (SC25) is expected to be weaker than the current Cycle 24. The model results show that a deep extended solar activity minimum is expected during 2019 - 2021, and that the next solar maximum will occur in 2024 - 2025. The sunspot number at the maximum will be about 50 with an error estimate of 15 - 30 %. The maximum will likely have a double peak or show extended periods (for 2 - 2.5 years) of high activity. According to the hemispheric prediction results, SC25 will start in 2020 in the Southern hemisphere, and will have a maximum in 2024 with a sunspot number of about 28. In the Northern hemisphere the cycle will be delayed for about 1 year (with an error of $\pm 0.5$ year), and reach a maximum in 2025 with a sunspot number of about 23.

We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize their atmospheric structure, composition, and circulation, from gas giants to rocky planets. However, exoplanet atmospheric models capable of interpreting the upcoming observations are often limited by insufficiencies in the laboratory and theoretical data that serve as critical inputs to atmospheric physical and chemical tools. Here we provide an up-to-date and condensed description of areas where laboratory and/or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry, building off a larger 2016 white paper, and endorsed by the NAS Exoplanet Science Strategy report. Now is the ideal time for progress in these areas, but this progress requires better access to, understanding of, and training in the production of spectroscopic data as well as a better insight into chemical reaction kinetics both thermal and radiation-induced at a broad range of temperatures. Given that most published efforts have emphasized relatively Earth-like conditions, we can expect significant and enlightening discoveries as emphasis moves to the exotic atmospheres of exoplanets.

It is a truism within the exoplanet field that "to know the planet, you must know the star." This pertains to the physical properties of the star (i.e. mass, radius, luminosity, age, multiplicity), the activity and magnetic fields, as well as the stellar elemental abundances which can be used as a proxy for planetary composition. In this white paper, we discuss important stellar characteristics that require attention in upcoming ground- and space-based missions, such that their processes can be understood and either detangled from that of the planet, correlated with the presence of a planet, or utilized in lieu of direct planetary observations.

Surveys reveal that terrestrial- to Neptune-sized planets (1 $< R <$ 4 R$_{\rm{Earth}}$) are the most common type of planets in our galaxy. Detecting and characterizing such small planets around nearby stars holds the key to understanding the diversity of exoplanets and will ultimately address the ubiquitousness of life in the universe. The following fundamental questions will drive research in the next decade and beyond: (1) how common are terrestrial to Neptune-sized planets within a few AU of their host star, as a function of stellar mass? (2) How does planet composition depend on planet mass, orbital radius, and host star properties? (3) What are the energy budgets, atmospheric dynamics, and climates of the nearest worlds? Addressing these questions requires: a) diffraction-limited spatial resolution; b) stability and achievable contrast delivered by adaptive optics; and c) the light-gathering power of extremely large telescopes (ELTs), as well as multi-wavelength observations and all-sky coverage enabled by a comprehensive US ELT Program. Here we provide an overview of the challenge, and promise of success, in detecting and comprehensively characterizing small worlds around the very nearest stars to the Sun with ELTs. This white paper extends and complements the material presented in the findings and recommendations published in the National Academy reports on Exoplanet Science Strategy and Astrobiology Strategy for the Search for Life in the Universe.

The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditions favorable for the origin, development and sustainment of life as we know it. This requires the understanding of global (astrospheric) and local (atmospheric, surface and internal) environments of exoplanets in the framework of the physical processes of the interaction between evolving planet-hosting stars along with exoplanetary evolution over geological timescales, and the resulting impact on climate and habitability of exoplanets. Feedbacks between astrophysical, physico-chemical atmospheric and geological processes can only be understood through interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary, Earth sciences, astrobiology, and the origin of life communities. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets and potential exomoons around them may significantly modify the extent and the location of the habitable zone and provide new directions for searching for signatures of life. Thus, characterization of stellar ionizing outputs becomes an important task for further understanding the extent of habitability in the universe. The goal of this white paper is to identify and describe promising key research goals to aid the theoretical characterization and observational detection of ionizing radiation from quiescent and flaring upper atmospheres of planet hosts as well as properties of stellar coronal mass ejections and stellar energetic particle events.

Transmission spectra probe the atmospheres of transiting exoplanets, but these observations are also subject to signals introduced by magnetic active regions on host stars. Here we outline scientific opportunities in the next decade for providing useful constraints on stellar photospheres and inform interpretations of transmission spectra of the smallest ($R<4\,R_{\odot}$) exoplanets. We identify and discuss four primary opportunities: (1) refining stellar magnetic active region properties through exoplanet crossing events; (2) spectral decomposition of active exoplanet host stars; (3) joint retrievals of stellar photospheric and planetary atmospheric properties with studies of transmission spectra; and (4) continued visual transmission spectroscopy studies to complement longer-wavelength studies from $\textit{JWST}$. We make five recommendations to the Astro2020 Decadal Survey Committee: (1) identify the transit light source (TLS) effect as a challenge to precise exoplanet transmission spectroscopy and an opportunity ripe for scientific advancement in the coming decade; (2) include characterization of host star photospheric heterogeneity as part of a comprehensive research strategy for studying transiting exoplanets; (3) support the construction of ground-based extremely large telescopes (ELTs); (4) support multi-disciplinary research teams that bring together the heliophysics, stellar physics, and exoplanet communities to further exploit transiting exoplanets as spatial probes of stellar photospheres; and (5) support visual transmission spectroscopy efforts as complements to longer-wavelength observational campaigns with $\textit{JWST}$.

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.

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

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.

Spectra of the cellular photospheric flows are determined from full-disk Doppler velocity observations acquired by the Helioseismic and Magnetic Imager (HMI) instrument on the Solar Dynamics Observatory (SDO) spacecraft. Three different analysis methods are used to separately determine spectral coefficients representing the poloidal flows, the toroidal flows, and the radial flows. The amplitudes of these spectral coefficients are constrained by simulated data analyzed with the same procedures as the HMI data. We find that the total velocity spectrum rises smoothly to a peak at a wavenumber of about 120 (wavelength of about 35 Mm), which is typical of supergranules. The spectrum levels off out to wavenumbers of about 400, and then rises again to a peak at a wavenumber of about 3500 (wavelength of about 1200 km), which is typical of granules. The velocity spectrum is dominated by the poloidal flow component (horizontal flows with divergence but no curl) at wavenumbers above 30. The toroidal flow component (horizontal flows with curl but no divergence) dominates at wavenumbers less than 30. The radial flow velocity is only about 3\% of the total flow velocity at the lowest wavenumbers, but increases in strength to become about 50\% at wavenumbers near 4000. The spectrum compares well with the spectrum of giant cell flows at the lowest wavenumbers and with the spectrum of granulation from a 3D radiative-hydrodynamic simulation at the highest wavenumbers.

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.

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.

The solar atmosphere is extremely dynamic, and many important phenomena develop on small scales that are unresolved in observations with the Helioseismic and Magnetic Imager (HMI) instrument on the Solar Dynamics Observatory (SDO). For correct calibration and interpretation of the observations, it is very important to investigate the effects of small-scale structures and dynamics on the HMI observables, such as Doppler shift, continuum intensity, spectral line depth, and width. We use 3D radiative hydrodynamics simulations of the upper turbulent convective layer and the atmosphere of the Sun, and a spectro-polarimetric radiative transfer code to study observational characteristics of the Fe I 6173A line observed by HMI in quiet-Sun regions. We use the modeling results to investigate the sensitivity of the line Doppler shift to plasma velocity, and also sensitivities of the line parameters to plasma temperature and density, and determine effective line formation heights for observations of solar regions located at different distances from the disc center. These estimates are important for the interpretation of helioseismology measurements. In addition, we consider various center-to-limb effects, such as convective blue-shift, variations of helioseismic travel-times, and the 'concave' Sun effect, and show that the simulations can qualitatively reproduce the observed phenomena, indicating that these effects are related to a complex interaction of the solar dynamics and radiative transfer.

Highly turbulent nature of convection on the Sun causes strong multi-scale interaction of subsurface layers with the photosphere and chromosphere. According to realistic 3D radiative MHD numerical simulations ubiquitous small-scale vortex tubes are generated by turbulent flows below the visible surface and concentrated in the intergranular lanes. The vortex tubes can capture and amplify magnetic field, penetrate into chromospheric layers and initiate quasi-periodic flow eruptions that generates Alfv\'enic waves, transport mass and energy into the solar atmosphere. The simulations revealed high-speed flow patterns, and complicated thermodynamic and magnetic structures in the erupting vortex tubes. The spontaneous eruptions are initiated and driven by strong pressure gradients in the near-surface layers, and accelerated by the Lorentz force in the low chromosphere. In this paper, the simulation data are used to further investigate the dynamics of the eruptions, their spectro-polarimetric characteristics for the Fe I 6301.5 and 6302.5 A spectral lines, and demonstrate expected signatures of the eruptions in the Hinode SP data. We found that the complex dynamical structure of vortex tubes (downflows in the vortex core and upflows on periphery) can be captured by the Stokes I profiles. During an eruption, the ratio of down and upflows can suddenly change, and this effect can be observed in the Stokes V profile. Also, during the eruption the linear polarization signal increases, and this also can be detected with Hinode SP.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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

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.

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.

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.

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.