There has been tremendous progress in the degree of realism of three-dimensional radiation magneto-hydrodynamic simulations of the solar atmosphere in the past decades. Four of the most frequently used numerical codes are Bifrost, CO5BOLD, MANCHA3D, and MURaM. Here we test and compare the wave propagation characteristics in model runs from these four codes by measuring the dispersion relation of acoustic-gravity waves at various heights. We find considerable differences between the various models. The height dependence of wave power, in particular of high-frequency waves, varies by up to two orders of magnitude between the models, and the phase difference spectra of several models show unexpected features, including $\pm180^\circ$ phase jumps.

The analyses of sunspot observations revealed a fundamental magnetic property of the umbral boundary, the invariance of the vertical component of the magnetic field. We aim to analyse the magnetic properties of the umbra-penumbra boundary in simulated sunspots and thus assess their similarity to observed sunspots. Also, we aim to investigate the role of plasma $\beta$ and the ratio of kinetic to magnetic energy in simulated sunspots on the convective motions. We use a set of non-grey simulation runs of sunspots with the MURaM code. These data are used to synthesise the Stokes profiles that are then degraded to the Hinode spectropolarimeter-like observations. Then, the data are treated like real Hinode observations of a sunspot and magnetic properties at the umbral boundaries are determined. Simulations with potential field extrapolation produce a realistic magnetic field configuration on their umbral boundaries. Two simulations with potential field upper boundary, but different subsurface magnetic field structures, differ significantly in the extent of their penumbrae. Increasing the penumbra width by forcing more horizontal magnetic fields at the upper boundary results in magnetic properties that are not consistent with observations. This implies that the size of the penumbra is given by the subsurface structure of the magnetic field. None of the sunspot simulations is consistent with observed properties of the magnetic field and direction of the Evershed flow at the same time. Strong outward directed Evershed flows are only found in setups with artificially enhanced horizontal component of the magnetic field at the top boundary that are not consistent with the observed magnetic field properties at the UP boundary. We want to stress out that the `photospheric' boundary of simulated sunspots is defined by a magnetic field strength of equipartition field value.

While the quiet Sun magnetic field shows only little variation with the solar cycle, long-term variations cannot be completely ruled out from first principles. We investigate the potential effect of quiet Sun magnetism on spectral solar irradiance through a series of small-scale dynamo simulations with zero vertical flux imbalance ($\langle B_z\rangle=0$) and varying levels of small-scale magnetic field strength, and one weak network case with an additional flux imbalance corresponding to a flux density of $\langle B_z\rangle=100$ G. From these setups we compute the dependence of the outgoing radiative energy flux on the mean vertical magnetic field strength in the photosphere at continuum optical depth $\tau=1$ ($\langle \vert B_z\vert\rangle_{\tau=1}$). We find that a quiet Sun setup with a mean vertical field strength of $\langle \vert B_z\vert\rangle_{\tau=1}=69$ G is about $0.6~\%$ brighter than a non-magnetic reference case. We find a linear dependence of the outgoing radiative energy flux on the mean field strength $\langle \vert B_z\vert\rangle_{\tau=1}$ with a relative slope of $1.4\cdot 10^{-4}$ G$^{-1}$. With this sensitivity, only a moderate change of the quiet Sun field strength by $10\%$ would lead to a total solar irradiance variation comparable to the observed solar cycle variation. While this does provide strong indirect constraints on possible quiet Sun variations during a regular solar cycle, it also emphasizes that potential variability over longer time scales could make a significant contribution to longer-term solar irradiance variations.

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

Inversion codes for the polarized radiative transfer equation can be used to infer the temperature $T$, line-of-sight velocity $v_{\rm los}$, and magnetic field $\rm{\bf B}$ as a function of the continuum optical-depth $\tau_{\rm c}$. However, they do not directly provide the gas pressure $P_{\rm g}$ or density $\rho$. In order to obtain these latter parameters, inversion codes rely instead on the assumption of hydrostatic equilibrium (HE) in addition to the equation of state (EOS). Unfortunately, the assumption of HE is rather unrealistic across magnetic field lines. This is because the role of the Lorentz force, among other factors, is neglected. This translates into an inaccurate conversion from optical depth $\tau_{\rm c}$ to geometrical height $z$. We aim at improving this conversion via the application of magneto-hydrostatic (MHS) equilibrium instead of HE. We develop a method to solve the momentum equation under MHS equilibrium (i.e., taking the Lorentz force into account) in three dimensions. The method is based on the solution of a Poisson-like equation. Considering the gas pressure $P_{\rm g}$ and density $\rho$ from three-dimensional magneto-hydrodynamic (MHD) simulations of sunspots as a benchmark, we compare the results from the application of HE and MHS equilibrium. We find that HE retrieves the gas pressure and density within an order of magnitude of the MHD values in only about 47 \% of the domain. This translates into an error of about $160-200$ km in the determination of the $z-\tau_{\rm c}$ conversion. On the other hand, the application of MHS equilibrium allows determination of $P_{\rm g}$ and $\rho$ within an order of magnitude in 84 \% of the domain. In this latter case, the $z-\tau_{\rm c}$ conversion is obtained with an accuracy of $30-70$ km.

This textbook gives a perspective of heliophysics in a way that emphasizes universal processes from a perspective that draws attention to what provides Earth (and similar (exo-)planets) with a relatively stable setting in which life as we know it can thrive. The book is intended for students in physical sciences in later years of their university training and for beginning graduate students in fields of solar, stellar, (exo-)planetary, and planetary-system sciences.

Recently, there have been some reports of unusually strong photospheric magnetic fields (which can reach values of over 7 kG) inferred from Hinode SOT/SP sunspot observations within penumbral regions. These superstrong penumbral fields are even larger than the strongest umbral fields on record and appear to be associated with supersonic downflows. The finding of such fields has been controversial since they seem to show up only when spatially coupled inversions are performed. Here, we investigate and discuss the reliability of those findings by studying in detail observed spectra associated with particularly strong magnetic fields at the inner edge of the penumbra of active region 10930. We apply classical diagnostic methods and various inversions with different model atmospheres to the observed Stokes profiles in two selected pixels with superstrong magnetic fields, and compare the results with a magnetohydrodynamic simulation of a sunspot whose penumbra contains localized regions with strong fields (nearly 5 kG at $\tau=1$) associated with supersonic downflows...

Recent numerical simulations and observations of sunspots show a significant amount of opposite polarity magnetic field within the sunspot penumbra. Most of the opposite polarity field is associated with convective downflows. We present an analysis of 3D MHD simulations through forward modeling of synthetic Stokes profiles of the Fe\sci 6301.5 \AA~ and Fe\sci 6302.5 \AA~ lines). The synthetic Stokes profiles are spatially and spectrally degraded considering typical instrument properties. Line bisector shifts of the Fe\sci 6301.5 \AA~ line are used to determine line-of-sight velocities. Far wing magnetograms are constructed from the Stokes V profiles of the Fe\sci 6302.5 \AA~ line. While we find an overall good agreement between observations and simulations, the fraction of opposite polarity magnetic field, the downflow filling factor and the opposite polarity-downflow association are strongly affected by spatial smearing and presence of strong gradients in the line-of-sight magnetic field and velocity. A significant fraction of opposite polarity magnetic field and downflows are hidden in the observations due to typical instrumental noise. Comparing simulations that differ by more than a factor of two in grid spacing we find that these quantities are robust within the simulations.

We observe the abrupt end of solar activity cycles at the Sun's equator by combining almost 140 years of observations from ground and space. These "terminator" events appear to be very closely related to the onset of magnetic activity belonging to the next sunspot cycle at mid-latitudes and the polar-reversal process at high-latitudes. Using multi-scale tracers of solar activity we examine the timing of these events in relation to the excitation of new activity and find that the time taken for the solar plasma to communicate this transition is of the order of one solar rotation, but could be shorter. Utilizing uniquely comprehensive solar observations from the Solar Terrestrial Relations Observatory (STEREO), and Solar Dynamics Observatory (SDO) we see that this transitional event is strongly longitudinal in nature. Combined, these characteristics imply that magnetic information is communicated through the solar interior rapidly. A range of possibilities exist to explain such behavior: the presence of magnetic reconnection in the deep interior, internal gravity waves on the solar tachocline, or that the magnetic fields present in the Sun's convection zone could be very large, with a poloidal field strengths reaching 50k - considerably larger than conventional explorations of solar and stellar dynamos estimate. Regardless of mechanism responsible, the rapid timescales demonstrated by the Sun's global magnetic field reconfiguration present strong constraints on first-principles numerical simulations of the solar interior and, by extension, other stars.

Because of the complex physics that governs the formation of chromospheric lines, interpretation of solar chromospheric observations is difficult. The origin and characteristics of many chromospheric features are, because of this, unresolved. We focus here on studying two prominent features: long fibrils and flare ribbons. To model them, we use a 3D MHD simulation of an active region which self-consistently reproduces both of them. We model the H$\alpha$, Mg II k, Ca II K, and Ca II 8542 {\AA} lines using the 3D non-LTE radiative transfer code Multi3D. This simulation reproduces long fibrils that span between the opposite-polarity sunspots and go up to 4 Mm in height. They can be traced in all lines due to density corrugation. Opposite to previous studies, H$\alpha$, Mg II h&k, and Ca II H&K, are formed at similar height in this model. Magnetic field lines are aligned with the H$\alpha$ fibrils, but the latter holds to a lesser extent for the Ca II 8542 {\AA} line. The simulation shows structures in the H$\alpha$ line core that look like flare ribbons. The emission in the ribbons is caused by a dense chromosphere and a transition region at high column mass. The ribbons are visible in all chromospheric lines, but least prominent in Ca II 8542 {\AA} line. In some pixels, broad asymmetric profiles with a single emission peak are produced, similar to the profiles observed in flare ribbons. They are caused by a deep onset of the chromospheric temperature rise and large velocity gradients. The simulation produces long fibrils similar to what is seen in observations. It also produces structures similar to flare ribbons despite the lack of non-thermal electrons in the simulation. The latter suggests that thermal conduction might be a significant agent in transporting flare energy to the chromosphere in addition to non-thermal electrons.

We show that at large magnetic Prandtl numbers, the Lorentz force does work on the flow at small scales and drives fluid motions, whose energy is dissipated viscously. This situation is opposite to that in a normal dynamo, where the flow does work against the Lorentz force. We compute the spectral conversion rates between kinetic and magnetic energies for several magnetic Prandtl numbers and show that normal (forward) dynamo action occurs on large scales over a progressively narrower range of wavenumbers as the magnetic Prandtl number is increased. At higher wavenumbers, reversed dynamo action occurs, i.e., magnetic energy is converted back into kinetic energy at small scales. We demonstrate this in both direct numerical simulations forced by volume stirring and in large eddy simulations of solar convectively driven small-scale dynamos. Low density plasmas such as stellar coronae tend to have large magnetic Prandtl numbers, i.e., the viscosity is large compared with the magnetic diffusivity. The regime in which viscous dissipation dominates over resistive dissipation for large magnetic Prandtl numbers was also previously found in large eddy simulations of the solar corona, i.e., our findings are a more fundamental property of MHD that is not just restricted to dynamos. Viscous energy dissipation is a consequence of positive Lorentz force work, which may partly correspond to particle acceleration in close-to-collisionless plasmas. This is, however, not modeled in the MHD approximation employed. By contrast, resistive energy dissipation on current sheets is expected to be unimportant in stellar coronae.

Contex: Solar magnetic activity exhibits variations with periods between 1.5--4 years, the so-called quasi-biennial oscillations (QBOs), in addition to the well-known 11-year Schwabe cycles. Solar dynamo is thought to be the responsible mechanism for generation of the QBOs. Aims: In this work, we analyse sunspot areas to investigate the spatial and temporal behaviour of the QBO signal and study the responsible physical mechanisms using simulations from fully nonlinear mean-field flux-transport dynamos. Methods: We investigated the behaviour of the QBOs in the sunspot area data in full disk, and northern and southern hemispheres, using wavelet and Fourier analyses. We also ran solar dynamos with two different approaches to generating a poloidal field from an existing toroidal field, Babcock-Leighton and turbulent $\alpha$ mechanisms. We then studied the simulated magnetic field strengths as well as meridional circulation and differential rotation rates using the same methods. Results: The results from the sunspot areas show that the QBOs are present in the full disk and hemispheric sunspot areas and they show slightly different spatial and temporal behaviours, indicating a slightly decoupled solar hemispheres. The QBO signal is generally intermittent and in-phase with the sunspot area data, surfacing when the solar activity is in maximum. The results from the BL-dynamos showed that they are neither capable of generating the slightly decoupled behaviour of solar hemispheres nor can they generate QBO-like signals. The turbulent $\alpha$-dynamos, on the other hand, generated decoupled hemispheres and some QBO-like shorter cycles. Conclusions: In conclusion, our simulations show that the turbulent $\alpha$-dynamos with the Lorentz force seems more efficient in generating the observed temporal and spatial behaviour of the QBO signal compared with those from the BL-dynamos.

We analyze recent high resolution photospheric small-scale dynamo simulations that were computed with the MURaM radiative MHD code. We focus the analysis on newly forming downflow lanes in exploding granules since they show how weakly magnetized regions in the photosphere (center of granules) evolve into strongly magnetized regions (downflow lanes). We find that newly formed downflow lanes exhibit initially mostly a laminar converging flow that amplifies the vertical magnetic field embedded in the granule from a few 10 G to field strengths exceeding 800 G. This results in extended magnetic sheets that have a length comparable to granular scales. Field amplification by turbulent shear happens first a few 100 km beneath the visible layers of the photosphere. Shallow recirculation transports the resulting turbulent field into the photosphere within minutes, after which the newly formed downflow lane shows a mix of strong magnetic sheets and turbulent field components. We stress in particular the role of shallow and deep recirculation for the organization and strength of magnetic field in the photosphere and discuss the photospheric and sub-photospheric energy conversion associated with the small-scale dynamo process. While the energy conversion through the Lorentz force depends only weakly on the saturation field strength (and therefore deep or shallow recirculation), it is strongly dependent on the magnetic Prandtl number. We discuss the potential of these findings for further constraining small-scale dynamo models through high resolution observations.

The motions of small-scale magnetic flux elements in the solar photosphere can provide some measure of the Lagrangian properties of the convective flow. Measurements of these motions have been critical in estimating the turbulent diffusion coefficient in flux-transport dynamo models and in determining the Alfven wave excitation spectrum for coronal heating models. We examine the motions of internetwork flux elements in a 24 hour long Hinode/NFI magnetogram sequence with 90 second cadence, and study both the scaling of their mean squared displacement and the shape of their displacement probability distribution as a function of time. We find that the mean squared displacement scales super-diffusively with a slope of about 1.48. Super-diffusive scaling has been observed in other studies for temporal increments as small as 5 seconds, increments over which ballistic scaling would be expected. Using high-cadence MURaM simulations, we show that the observed super-diffusive scaling at short temporal increments is a consequence of random changes in the barycenter positions due to flux evolution. We also find that for long temporal increments, beyond granular lifetimes, the observed displacement distribution deviates from that expected for a diffusive process, evolving from Rayleigh to Gaussian. This change in the distribution can be modeled analytically by accounting for supergranular advection along with motions due to granulation. These results complicate the interpretation of magnetic element motions as strictly advective or diffusive on short and long timescales and suggest that measurements of magnetic element motions must be used with caution in turbulent diffusion or wave excitation models. We propose that passive trace motions in measured photospheric flows may yield more robust transport statistics.

There have been a few reports in the literature of counter-Evershed flows observed in well developed sunspot penumbrae, i.e. flows directed towards the umbra along penumbral filaments. Here we investigate the driving forces of such counter-Evershed flows in a radiative magnetohydrodynamic simulation of a sunspot and compare them with the forces acting on the normal Evershed flow. The simulation covers a timespan of 100 solar hours and generates an Evershed outflow exceeding 8 km s$^{-1}$ in the penumbra along radially aligned filaments where the magnetic field is almost horizontal. Additionally, the simulation produces a fast counter-Evershed flow (i.e., an inflow near $\tau = 1$) in some regions within the penumbra, reaching peak flow speeds of $\sim$12 km s$^{-1}$. The counter-Evershed flows are transient and typically last a few hours before they turn into outflows again. By using the kinetic energy equation and evaluating its various terms in the simulation box, we found that the Evershed flow occurs due to overturning convection in a strongly inclined magnetic field while the counter-Evershed flows can be well described as siphon flows.

We aim to investigate the nature and occurrence characteristics of grand solar minimum and maximum periods, which are observed in the solar proxy records such as 10Be and 14C, using a fully non-linear Babcock-Leighton type flux-transport dynamo including momentum and entropy equations. The differential rotation and meridional circulation are generated from the effect of turbulent Reynolds stress and are subjected to back-reaction from the magnetic field. To generate grand minimum and maximum-like periods in our simulations, we used random fluctuations in the angular momentum transport process, namely the Lambda-mechanism, and in the Babcock-Leighton mechanism. To characterise the nature and occurrences of the identified grand minima and maxima in our simulations, we used the waiting time distribution analyses, which reflects whether the underlying distribution arises from a random or a memory-bearing process. The results show that, in majority of the cases, the distributions of grand minima and maxima reveal that the nature of these events originates from memoryless processes. We also found that in our simulations the meridional circulation speed tends to be smaller during grand maximum, while it is faster during grand minimum periods. The radial differential rotation tend to be larger during grand maxima, while it is smaller during grand minima. The latitudinal differential rotation on the other hand is found to be larger during grand minima.

We present a realistic numerical model of sunspot and active region formation based on the emergence of flux bundles generated in a solar convective dynamo. To this end we use the magnetic and velocity fields in a horizontal layer near the top boundary of the solar convective dynamo simulation to drive realistic radiative-magnetohydrodynamic simulations of the upper most layers of the convection zone. The main results are: (1) The emerging flux bundles rise with the mean speed of convective upflows, and fragment into small-scale magnetic elements that further rise to the photosphere, where bipolar sunspot pairs are formed through the coalescence of the small-scale magnetic elements. (2) Filamentary penumbral structures form when the sunspot is still growing through ongoing flux emergence. In contrast to the classical Evershed effect, the inflow seems to prevail over the outflow in a large part of the penumbra. (3) A well formed sunspot is a mostly monolithic magnetic structures that is anchored in a persistent deep-seated downdraft lane. The flow field outside the spot shows a giant vortex ring that comprises of an inflow below 15 Mm depth and an outflow above 15 Mm depth. (4) The sunspots successfully reproduce the fundamental properties of the observed solar active regions, including the more coherent leading spots with a stronger field strength, and the correct tilts of bipolar sunspot pairs. These asymmetries can be linked to the intrinsic asymmetries in the magnetic and flow fields adapted from the convective dynamo simulation.

We present the scientific case for a future space-based telescope aimed at very high spatial and temporal resolution imaging of the solar photosphere and chromosphere. Previous missions (e.g., HINODE, SUNRISE) have demonstrated the power of observing the solar photosphere and chromosphere at high spatial resolution without contamination from Earth's atmosphere. We argue here that increased spatial resolution (from currently 70 km to 25 km in the future) and high temporal cadence of the observations will vastly improve our understanding of the physical processes controlling solar magnetism and its characteristic scales. This is particularly important as the Sun's magnetic field drives solar activity and can significantly influence the Sun-Earth system. At the same time a better knowledge of solar magnetism can greatly improve our understanding of other astrophysical objects.

We present a new version of the MURaM radiative MHD code that allows for simulations spanning from the upper convection zone into the solar corona. We implemented the relevant coronal physics in terms of optically thin radiative loss, field aligned heat conduction and an equilibrium ionization equation of state. We artificially limit the coronal Alfv{\'e}n and heat conduction speeds to computationally manageable values using an approximation to semi-relativistic MHD with an artificially reduced speed of light (Boris correction). We present example solutions ranging from quiet to active Sun in order to verify the validity of our approach. We quantify the role of numerical diffusivity for the effective coronal heating. We find that the (numerical) magnetic Prandtl number determines the ratio of resistive to viscous heating and that owing to the very large magnetic Prandtl number of the solar corona, heating is expected to happen predominantly through viscous dissipation. We find that reasonable solutions can be obtained with values of the reduced speed of light just marginally larger than the maximum sound speed. Overall this leads to a fully explicit code that can compute the time evolution of the solar corona in response to photospheric driving using numerical time steps not much smaller than $0.1$ seconds. Numerical simulations of the coronal response to flux emergence covering a time span of a few days are well within reach using this approach.

We try to retrieve the power spectra with certainty to the highest spatial frequencies allowed by current instrumentation. For this, we use 2D inversion code that were able to recover information up to the instrumental diffraction limit. The retrieved power spectra have shallow slopes extending further down to much smaller scales than found before. They seem not to show any power law. The observed slopes at subgranular scales agree with those obtained from recent local dynamo simulations. Small differences are found for vertical component of kinetic energy that suggest that observations suffer from an instrumental effect that is not taken into account.

Magnetic field emerges at the surface of the Sun as sunspots and active regions. This process generates a poloidal magnetic field from a rising toroidal flux tube, it is a crucial but poorly understood aspect of the solar dynamo. The emergence of magnetic field is also important because it is a key driver of solar activity. We show that measurements of horizontal flows at the solar surface around emerging active regions, in combination with numerical simulations of solar magnetoconvection, can constrain the subsurface rise speed of emerging magnetic flux. The observed flows imply that the rise speed of the magnetic field is no larger than 150 m/s at a depth of 20 Mm, that is, well below the prediction of the (standard) thin flux tube model but in the range expected for convective velocities at this depth. We conclude that convective flows control the dynamics of rising flux tubes in the upper layers of the Sun and cannot be neglected in models of flux emergence.

Properties of magnetic field in the internetwork regions are still fairly unknown due to rather weak spectropolarimetric signals. We address the matter by using the 2D inversion code that is able to retrieve the information on smallest spatial scales, up to the diffraction limit, while being less susceptible to noise than most of the previous methods used. Performance of the code and the impact of the various effects on the retrieved field distribution is tested first on the realistic MHD simulations. The best inversion scenario is then applied to the real Hinode/SP data. Tests on simulations show: (1) the best choice of node position ensures a decent retrieval of all parameters, (2) code performs well for different configurations of magnetic field, (3) slightly different noise level or slightly different defocus included in the spatial PSF produces no significant effect on the results and (4) temporal integration shifts the field distribution to the stronger, more horizontally inclined field. Although the contribution of the weak field is slightly overestimated due to noise, the 2D inversions are able to recover well the overall distribution of the magnetic field strength. Application of the 2D inversion code on the Hinode/SP internetwork observations reveals a monotonic field strength distribution. The mean field strength at optical depth unity is $\sim 130$~G. At the higher layers, field strength drops as the field becomes more horizontal. Regarding the distribution of the field inclination, tests show that we cannot directly retrieve it with the observations/tools at hand, however the obtained distributions are consistent with those expected from simulations with a quasi-isotropic field inclination after accounting for observational effects.

We present a series of high-resolution sunspot simulations that cover a time span of up to 100 hours. The simulation domain extends about 18 Mm in depth beneath the photosphere and 98 Mm horizontally. We use open boundary conditions that do not maintain the initial field structure against decay driven by convective motions. We consider two setups: A sunspot simulation with penumbra, and a "naked-spot" simulation in which we removed the penumbra after 20 hours through a change in the magnetic top boundary condition. While the sunspot has an Evershed outflow of 3-4 km/s, the naked spot is surrounded by an inflow of 1-2 km/s in close proximity. However, both spots are surrounded by an outflow on larger scales with a few 100 m/s flow speed in the photosphere. While the sunspot has almost constant magnetic flux content for the simulated time span of 3-4 days, the naked spot decays steadily at a rate of $10^{21}$ Mx/day. A region with reduced downflow filling factor, which is more extended for the sunspot, surrounds both spots. The absence of downflows perturbs the upflow/downflow massflux balance and leads to a large-scale radially overturning flow system, the photospheric component of this flow is to the observable moat flow. The reduction of the downflow filling factor also inhibits submergence of magnetic field in the proximity of the spots, which stabilizes them against decay. While this effect is present for both spots, it is more pronounced for the sunspot and explains the almost stationary magnetic flux content.

We investigate small-scale dynamo action in the solar convection zone through a series of high resolution MHD simulations in a local Cartesian domain with 1$R_\odot$ (solar radius) of horizontal extent and a radial extent from 0.715 to 0.96$R_\odot$. The dependence of the solution on resolution and diffusivity is studied. For a grid spacing of less than 350 km, the root mean square magnetic field strength near the base of the convection zone reaches 95% of the equipartition field strength (i.e. magnetic and kinetic energy are comparable). For these solutions the Lorentz force feedback on the convection velocity is found to be significant. The velocity near the base of the convection zone is reduced to 50% of the hydrodynamic one. In spite of a significant decrease of the convection velocity, the reduction in the enthalpy flux is relatively small, since the magnetic field also suppresses the horizontal mixing of the entropy between up- and downflow regions. This effect increases the amplitude of the entropy perturbation and makes convective energy transport more efficient. We discuss potential implications of these results for solar global convection and dynamo simulations.

We present a high-resolution, highly stratified numerical simulation of rotating thermal convection in a spherical shell. Our aim is to study in detail the processes that can maintain a near surface shear layer (NSSL) as inferred from helioseismology. Using the reduced speed of sound technique we can extend our global convection simulation to $0.99\,R_{\odot}$ and include near the top of our domain small-scale convection with short time scales that is only weakly influenced by rotation. We find the formation of a NSSL preferentially in high latitudes in the depth range $r=0.95-0.975R_\odot$. The maintenance mechanisms are summarized as follows. Convection under weak influence of rotation leads to Reynolds stresses that transport angular momentum radially inward in all latitudes. This leads to the formation of a strong poleward directed meridional flow and a NSSL, which is balanced in the meridional plane by forces resulting from the $\langle v'_r v'_\theta\rangle$ correlation of turbulent velocities. The origin of the required correlations depends to some degree on latitude. In high latitudes a positive correlation $\langle v'_rv'_\theta\rangle$ is induced in the NSSL by the poleward meridional flow whose amplitude increases with the radius, while a negative correlation is generated by the Coriolis force in bulk of the convection zone. In low latitudes a positive correlation $\langle v'_rv'_\theta\rangle$ results from rotationally aligned convection cells ("banana cells"). The force caused by these Reynolds stresses is in balance with the Coriolis force in the NSSL.

We revisit an old question: what are the effects of observing stratified atmospheres on scales below a photon mean free path? The mean free path of photons emerging from the solar photosphere and chromosphere is near 100 km. Using current 1m-class telescopes, the mean free path is on the order of the angular resolution. But the Daniel K. Inoue Solar Telescope will have a diffraction limit of 0.020'' near the atmospheric cutoff at 310nm, corresponding to 14 km at the solar surface. Even a small amount of scattering in the source function leads to physical smearing due to this solar 'fog', with effects similar to a degradation of the telescope PSF. We discuss a unified picture that depends simply on the nature and amount of scattering in the source function. Scalings are derived from which the scattering in the solar atmosphere can be transcribed into an effective Strehl ratio, a quantity useful to observers. Observations in both permitted (e.g., Fe I 630.2 nm) and forbidden (Fe I 525.0 nm) lines will shed light on both instrumental performance as well as on small scale structures in the solar atmosphere.

Milne-Eddington (M-E) inversion codes for the radiative transfer equation are the most widely used tools to infer the magnetic field from observations of the polarization signals in photospheric and chromospheric spectral lines. Unfortunately, a comprehensive comparison between the different M-E codes available to the solar physics community is still missing, and so is a physical interpretation of their inferences. In this contribution we offer a comparison between three of those codes (VFISV, ASP/HAO, and HeLIx$^+$). These codes are used to invert synthetic Stokes profiles that were previously obtained from realistic non-grey three-dimensional magnetohydrodynamical (3D MHD) simulations. The results of the inversion are compared with each other and with those from the MHD simulations. In the first case, the M-E codes retrieve values for the magnetic field strength, inclination and line-of-sight velocity that agree with each other within $\sigma_B \leq 35$ (Gauss), $\sigma_\gamma \leq 1.2\deg$, and $\sigma_{\rm v} \leq 10$ ms$^{-1}$, respectively. Additionally, M-E inversion codes agree with the numerical simulations, when compared at a fixed optical depth, within $\sigma_B \leq 130$ (Gauss), $\sigma_\gamma \leq 5\deg$, and $\sigma_{\rm v} \leq 320$ ms$^{-1}$. Finally, we show that employing generalized response functions to determine the height at which M-E codes measure physical parameters is more meaningful than comparing at a fixed geometrical height or optical depth. In this case the differences between M-E inferences and the 3D MHD simulations decrease to $\sigma_B \leq 90$ (Gauss), $\sigma_\gamma \leq 3\deg$, and $\sigma_{\rm v} \leq 90$ ms$^{-1}$.

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

We model the solar horizontal velocity power spectrum at scales larger than granulation using a two-component approximation to the mass continuity equation. The model takes four times the density scale height as the integral (driving) scale of the vertical motions at each depth. Scales larger than this decay with height from the deeper layers. Those smaller are assumed to follow a Kolomogorov turbulent cascade, with the total power in the vertical convective motions matching that required to transport the solar luminosity in a mixing length formulation. These model components are validated using large scale radiative hydrodynamic simulations. We reach two primary conclusions: 1. The model predicts significantly more power at low wavenumbers than is observed in the solar photospheric horizontal velocity spectrum. 2. Ionization plays a minor role in shaping the observed solar velocity spectrum by reducing convective amplitudes in the regions of partial helium ionization. The excess low wavenumber power is also seen in the fully nonlinear three-dimensional radiative hydrodynamic simulations employing a realistic equation of state. This adds to other recent evidence suggesting that the amplitudes of large scale convective motions in the Sun are significantly lower than expected. Employing the same feature tracking algorithm used with observational data on the simulation output, we show that the observed low wavenumber power can be reproduced in hydrodynamic models if the amplitudes of large scale modes in the deep layers are artificially reduced. Since the large scale modes have reduced amplitudes, modes on the scale of supergranulation and smaller remain important to convective heat flux even in the deep layers, suggesting that small scale convective correlations are maintained through the bulk of the solar convection zone.

We present a series of radiative MHD simulations addressing the origin and distribution of mixed polarity magnetic field in the solar photosphere. To this end we consider numerical simulations that cover the uppermost 2-6 Mm of the solar convection zone and we explore scales ranging from 2 km to 25 Mm. We study how the strength and distribution of magnetic field in the photosphere and subsurface layers depend on resolution, domain size and boundary conditions. We find that 50% of the magnetic energy at the \tau=1 level comes from field with the less than 500 G strength and that 50% of the energy resides on scales smaller than about 100 km. While probability distribution functions are essentially independent of resolution, properly describing the spectral energy distribution requires grid spacings of 8 km or smaller. The formation of flux concentrations in the photosphere exceeding 1 kG requires a mean vertical field strength greater than 30-40 G at \tau=1. The filling factor of kG flux concentrations increases with overall domain size as magnetic field becomes organized by larger, longer lived flow structures. A solution with a mean vertical field strength of around 85 G at \tau=1 requires a subsurface RMS field strength increasing with depth at the same rate as the equipartition field strength. We consider this an upper limit for the quiet Sun field strength, which implies that most of the convection zone is magnetized close to equipartition. We discuss these findings in view of recent high-resolution spectropolarimetric observations of quiet Sun magnetism.

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

We carry out non-rotating high-resolution calculations of the solar global convection, which resolve convective scales of less than 10 Mm. To cope with the low Mach number conditions in the lower convection zone, we use the reduced speed of sound technique (RSST), which is simple to implement and requires only local communication in the parallel computation. In addition, the RSST allows us to expand the computational domain upward to about $0.99 R_{\odot}$ as it can also handle compressible flows. Using this approach, we study the solar convection zone on the global scale, including small-scale near-surface convection. In particular, we investigate the influence of the top boundary condition on the convective structure throughout the convection zone as well as on small-scale dynamo action. Our main conclusions are: 1. The small-scale downflows generated in the near-surface layer penetrate into deeper layers to some extent and excite small-scale turbulence in the region $>0.9R_\odot$, where $R_\odot$ is the solar radius. 2. In the deeper convection zone ($<0.9R_\odot$), the convection is not influenced by the location of the upper boundary. 3. Using an LES approach, we can achieve small-scale dynamo action and maintain a field of about $0.15-0.25 B_\mathrm{eq}$ throughout the convection zone, where $B_\mathrm{eq}$ is the equipartition magnetic field to the kinetic energy. 4. The overall dynamo efficiency varies significantly in the convection zone as a consequence of the downward directed Poynting flux and the depth variation of the intrinsic convective scales.

We present numerical simulations of active region scale flux emergence covering a time span of up to 6 days. Flux emergence is driven by a bottom boundary condition that advects a semi-torus of magnetic field with $1.7 \times 10^{22}$ Mx flux into the computational domain. The simulations show that, even in the absense of twist, the magnetic flux is able the rise through the upper $15.5$ Mm of the convection zone and emerge into the photosphere to form spots. We find that spot formation is sensitive to the persistence of upflows at the bottom boundary footpoints, i.e. a continuing upflow would prevent spot formation. In addition, the presence of a torus-aligned flow (such flow into the retrograde direction is expected from angular momentum conservation during the rise of flux ropes through the convection zone) leads to a significant asymmetry between the pair of spots, with the spot corresponding to the leading spot on the Sun being more axisymmetric and coherent, but also forming with a delay relative to the following spot. The spot formation phase transitions directly into a decay phase. Subsurface flows fragment the magnetic field and lead to intrusions of almost field free plasma underneath the photosphere. When such intrusions reach photospheric layers, the spot fragments. The time scale for spot decay is comparable to the longest convective time scales present in the simulation domain. We find that the dispersal of flux from a simulated spot in the first two days of the decay phase is consistent with self-similar decay by turbulent diffusion.

We investigate an intensification mechanism for the magnetic field near the base of the solar convection zone that does not rely on differential rotation. Such mechanism in addition to differential rotation has been suggested by studies of flux emergence, which typically require field strength in excess of those provided by differential rotation alone. We study here a process in which potential energy of the superadiabatically stratified convection zone is converted into magnetic energy. This mechanism, know as explosion of magnetic flux tubes, has been previously studied in the thin flux tube approximation as well as two-dimensional MHD simulations, we expand the investigation to three-dimensional MHD simulations. Our main result is that enough intensification can be achieved in a three-dimensional magnetic flux sheet as long as the spatial scale of the imposed perturbation normal to the magnetic field is sufficiently large. When this spatial scale is small, the flux sheet tends to rise toward the surface, resulting in a significant decrease of the magnetic field amplification.

We obtain lower limits on the amplitude of convective velocities in the deep solar convection zone based only on the observed properties of the differential rotation and meridional circulation together with simple and robust dynamical balances obtained from the fundamental MHD equations. The linchpin of the approach is the concept of gyroscopic pumping whereby the meridional circulation across isosurfaces of specific angular momentum is linked to the angular momentum transport by the convective Reynolds stress. We find that the amplitude of the convective velocity must be at least 30 m s$^{-1}$ in the upper CZ ($r \sim 0.95 R$) and at least 8 m s$^{-1}$ in the lower CZ ($r \sim 0.75 R$) in order to be consistent with the observed mean flows. Using the base of the near-surface shear layer as a probe of the rotational influence, we are further able to show that the characteristic length scale of deep convective motions must be no smaller than 5.5--30 Mm. These results are compatible with convection models but suggest that the efficiency of the turbulent transport assumed in advection-dominated flux-transport dynamo models is generally not consistent with the mean flows they employ.

The motion of dark striations across bright filaments in a sunspot penumbra has become an important new diagnostic of convective gas flows in penumbral filaments. The nature of these striations has, however, remained unclear. Here we present an analysis of small scale motions in penumbral filaments in both simulations and observations. The simulations, when viewed from above, show fine structure with dark lanes running outwards from the dark core of the penumbral filaments. The dark lanes either occur preferentially on one side or alternate between both sides of the filament. We identify this fine structure with transverse (kink) oscillations of the filament, corresponding to a sideways swaying of the filament. These oscillations have periods in the range of 5-7 min and propagate outward and downward along the filament. Similar features are found in observed G-band intensity time series of penumbral filaments in a sunspot located near disk center obtained by the Broadband Filter Imager (BFI) on board {\it Hinode}. We also find that some filaments show dark striations moving to both sides of the filaments. Based on the agreement between simulations and observations we conclude that the motions of these striations are caused by transverse oscillations of the underlying bright filaments.

MHD simulations of sunspots have successfully reproduced many aspects of sunspot fine structure as consequence of magneto convection in inclined magnetic field. We study how global sunspot properties and penumbral fine structure depend on the magnetic top boundary condition as well as on grid spacing. The overall radial extent of the penumbra is subject to the magnetic top boundary condition. All other aspects of sunspot structure and penumbral fine structure are resolved at an acceptable level starting from a grid resolution of 48 [24] km (horizontal [vertical]). We find that the amount of inverse polarity flux and the overall amount of overturning convective motions in the penumbra are robust with regard to both, resolution and boundary conditions. At photospheric levels Evershed flow channels are strongly magnetized. We discuss in detail the relation between velocity and magnetic field structure in the photosphere and point out observational consequences.

Context. The anelastic approximation is often adopted in numerical calculation with low Mach number, such as stellar internal convection. This approximation requires frequent global communication, because of an elliptic partial differential equation. Frequent global communication is negative factor for the parallel computing with a large number of CPUs. Aims. The main purpose of this paper is to test the validity of a method that artificially reduces the speed of sound for the compressible fluid equations in the context of stellar internal convection. The reduction of speed of sound allows for larger time steps in spite of low Mach number, while the numerical scheme remains fully explicit and the mathematical system is hyperbolic and thus does not require frequent global communication. Methods. Two and three dimensional compressible hydrodynamic equations are solved numerically. Some statistical quantities of solutions computed with different effective Mach numbers (due to reduction of speed of sound) are compared to test the validity of our approach. Results. Numerical simulations with artificially reduced speed of sound are a valid approach as long as the effective Mach number (based on the reduced speed of sound) remains less than 0.7.

We studied bright umbral dots (UDs) detected in a moderate size sunspot and compared their statistical properties to recent MHD models. The study is based on high resolution data recorded by the New Solar Telescope at the Big Bear Solar Observatory and 3D MHD simulations of sunspots. Observed UDs, living longer than 150 s, were detected and tracked in a 46 min long data set, using an automatic detection code. Total 1553 (620) UDs were detected in the photospheric (low chromospheric) data. Our main findings are: i) none of the analyzed UDs is precisely circular, ii) the diameter-intensity relationship only holds in bright umbral areas, and iii) UD velocities are inversely related to their lifetime. While nearly all photospheric UDs can be identified in the low chromospheric images, some small closely spaced UDs appear in the low chromosphere as a single cluster. Slow moving and long living UDs seem to exist in both the low chromosphere and photosphere, while fast moving and short living UDs are mainly detected in the photospheric images. Comparison to the 3D MHD simulations showed that both types of UDs display, on average, very similar statistical characteristics. However, i) the average number of observed UDs per unit area is smaller than that of the model UDs, and ii) on average, the diameter of model UDs is slightly larger than that of observed ones.

Numerical simulations indicate that the filamentation of sunspot penumbrae and the associated systematic outflow (the Evershed effect) are due to convectively driven fluid motions constrained by the inclined magnetic field. We investigate whether these motions, in particular the upflows in the bright filaments and the downflows at their edges can be reliably observed with existing instrumentation. We use a snapshot from a sunspot simulation to calculate 2D maps of synthetic line profiles for the spectral lines Fe\sci 7090.4 \AA ~ and C\sci 5380.34 \AA. The maps are spatially and spectrally degraded according to typical instrument properties. Line-of-sight velocities are determined from line bisector shifts. We find that the detectability of the convective flows is strongly affected by spatial smearing, particularly so for the downflows. Furthermore, the line-of-sight velocities are dominated by the Evershed flow unless the observation is made very near to disk center. These problems may have compromised recent attempts to detect overturning penumbral convection. Lines with a low formation height are best suited to detect the convective flows.

We present a series of numerical sunspot models addressing the subsurface field and flow structure in up to 16 Mm deep domains covering up to 2 days of temporal evolution. Changes in the photospheric appearance of the sunspots are driven by subsurface flows in several Mm depth. Most of magnetic field is pushed into a downflow vertex of the subsurface convection pattern, while some fraction of the flux separates from the main trunk of the spot. Flux separation in deeper layers is accompanied in the photosphere with light bridge formation in the early stages and formation of pores separating from the spot at later stages. Over a time scale of less than a day we see the development of a large scale flow pattern surrounding the sunspots, which is dominated by a radial outflow reaching about 50% of the convective rms velocity in amplitude. Several components of the large scale flow are found to be independent from the presence of a penumbra and the associated Evershed flow. While the simulated sunspots lead to blockage of heat flux in the near surface layers, we do not see compelling evidence for a brightness enhancement in their periphery. We further demonstrate that the influence of the bottom boundary condition on the stability and long-term evolution of the sunspot is significantly reduced in a 16 Mm deep domain compared to the shallower domains considered previously.

The stratification near the base of the Sun's convective envelope is governed by processes of convective overshooting and element diffusion, and the region is widely believed to play a key role in the solar dynamo. The stratification in that region gives rise to a characteristic signal in the frequencies of solar p modes, which has been used to determine the depth of the solar convection zone and to investigate the extent of convective overshoot. Previous helioseismic investigations have shown that the Sun's spherically symmetric stratification in this region is smoother than that in a standard solar model without overshooting, and have ruled out simple models incorporating overshooting, which extend the region of adiabatic stratification and have a more-or-less abrupt transition to subadiabatic stratification at the edge of the overshoot region. In this paper we consider physically motivated models which have a smooth transition in stratification bridging the region from the lower convection zone to the radiative interior beneath. We find that such a model is in better agreement with the helioseismic data than a standard solar model.

We analyze in detail the penumbral structure found in a recent radiative MHD simulation. Near tau=1, the simulation produces penumbral fine structure consistent with the observationally inferred interlocking comb structure. Fast outflows exceeding 8 km/s are present along almost horizontal stretches of the magnetic field; in the outer half of the penumbra, we see opposite polarity flux indicating flux returning beneath the surface. The bulk of the penumbral brightness is maintained by small-scale motions turning over on scales shorter than the length of a typical penumbral filament. The resulting vertical rms velocity at tau=1 is about half of that found in the quiet Sun. Radial outflows in the sunspot penumbra have two components. In the uppermost few 100 km, fast outflows are driven primarily through the horizontal component of the Lorentz force, which is confined to narrow boundary layers beneath tau=1, while the contribution from horizontal pressure gradients is reduced in comparison to granulation as a consequence of anisotropy. The resulting Evershed flow reaches its peak velocity near tau=1 and falls off rapidly with height. Outflows present in deeper layers result primarily from a preferred ring-like alignment of convection cells surrounding the sunspot. These flows reach amplitudes of about 50% of the convective rms velocity rather independent of depth. A preference for the outflow results from a combination of Lorentz force and pressure driving. While the Evershed flow dominates by velocity amplitude, most of the mass flux is present in deeper layers and likely related to a large-scale moat flow.

Sunspot fine structure has been modeled in the past by a combination of idealized magneto-convection simulations and simplified models that prescribe the magnetic field and flow structure to a large degree. Advancement in numerical methods and computing power has enabled recently 3D radiative MHD simulations of entire sunspots with sufficient resolution to address details of umbral dots and penumbral filaments. After a brief review of recent developments we focus on the magneto-convective processes responsible for the complicated magnetic structure of the penumbra and the mechanisms leading to the driving of strong horizontal outflows in the penumbra (Evershed effect). The bulk of energy and mass is transported on scales smaller than the radial extent of the penumbra. Strong horizontal outflows in the sunspot penumbra result from a redistribution of kinetic energy preferring flows along the filaments. This redistribution is facilitated primarily through the Lorentz force, while horizontal pressure gradients play only a minor role. The Evershed flow is strongly magnetized: While we see a strong reduction of the vertical field, the horizontal field component is enhanced within filaments.

A comprehensive understanding of the solar magnetic cycle requires detailed modeling of the solar interior including the maintenance and variation of large scale flows (differential rotation and meridional flow), the solar dynamo and the flux emergence process connecting the magnetic field in the solar convection zone with magnetic field in the photosphere and above. Due to the vast range of time and length scales encountered, a single model of the entire convection zone is still out of reach. However, a variety of aspects can be modeled through a combined approach of 3D MHD models and simplified descriptions. We will briefly review our current theoretical understanding of these processes based on numerical models of the solar interior.

The Mount Wilson Ca HK survey revealed magnetic activity variations in a large sample of solar-type stars with timescales ranging from 2.5 to 25 years. This broad range of cycle periods is thought to reflect differences in the rotational properties and the depths of the surface convection zones for stars with various masses and ages. In 2007 we initiated a long-term monitoring campaign of Ca II H and K emission for a sample of 57 southern solar-type stars to measure their magnetic activity cycles and their rotational properties when possible. We report the discovery of a 1.6-year magnetic activity cycle in the exoplanet host star iota Horologii, and we obtain an estimate of the rotation period that is consistent with Hyades membership. This is the shortest activity cycle so far measured for a solar-type star, and may be related to the short-timescale magnetic variations recently identified in the Sun and HD49933 from helio- and asteroseismic measurements. Future asteroseismic observations can be compared to those obtained near the magnetic minimum in 2006 to search for cycle-induced shifts in the oscillation frequencies. If such short activity cycles are common in F stars, then NASA's Kepler mission should observe their effects in many of its long-term asteroseismic targets.

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

We present a radiative magnetohydrodynamics simulation of the formation of an Active Region on the solar surface. The simulation models the rise of a buoyant magnetic flux bundle from a depth of 7.5 Mm in the convection zone up into the solar photosphere. The rise of the magnetic plasma in the convection zone is accompanied by predominantly horizontal expansion. Such an expansion leads to a scaling relation between the plasma density and the magnetic field strength such that $B\propto\varrho^{1/2}$. The emergence of magnetic flux into the photosphere appears as a complex magnetic pattern, which results from the interaction of the rising magnetic field with the turbulent convective flows. Small-scale magnetic elements at the surface first appear, followed by their gradual coalescence into larger magnetic concentrations, which eventually results in the formation of a pair of opposite polarity spots. Although the mean flow pattern in the vicinity of the developing spots is directed radially outward, correlations between the magnetic field and velocity field fluctuations allow the spots to accumulate flux. Such correlations result from the Lorentz-force driven, counter-streaming motion of opposite-polarity fragments. The formation of the simulated Active Region is accompanied by transient light bridges between umbrae and umbral dots. Together with recent sunspot modeling, this work highlights the common magnetoconvective origin of umbral dots, light bridges and penumbral filaments.

The Mount Wilson Ca HK survey revealed magnetic activity variations in a large sample of solar-type stars with timescales ranging from 2.5 to 25 years. This broad range of cycle periods is thought to reflect differences in the rotational properties and the depths of the surface convection zones for stars with various masses and ages. Asteroseismic data will soon provide direct measurements of these quantities for individual stars, but many of the most promising targets are in the southern sky (e.g., alpha Cen A & B, beta Hyi, mu Ara, tau Cet, nu Ind), while long-term magnetic activity cycle surveys are largely confined to the north. In 2007 we began using the SMARTS 1.5-m telescope to conduct a long-term monitoring campaign of Ca II H & K emission for a sample of 57 southern solar-type stars to measure their magnetic activity cycles and their rotational properties when possible. This sample includes the most likely southern asteroseismic targets to be observed by the Stellar Oscillations Network Group (SONG), currently scheduled to begin operations in 2012. We present selected results from the first two years of the survey, and from the longer time baseline sampled by a single-epoch survey conducted in 1992.

Sunspots are concentrations of magnetic field on the visible solar surface that strongly affect the convective energy transport in their interior and surroundings. The filamentary outer regions (penumbrae) of sunspots show systematic radial outward flows along channels of nearly horizontal magnetic field. These flows were discovered 100 years ago and are present in all fully developed sunspots. Using a comprehensive numerical simulation of a sunspot pair, we show that penumbral structures with such outflows form when the average magnetic field inclination to the vertical exceeds about 45 degrees. The systematic outflows are a component of the convective flows that provide the upward energy transport and result from anisotropy introduced by the presence of the inclined magnetic field.

Results of a 3D MHD simulation of a sunspot with a photospheric size of about 20 Mm are presented. The simulation has been carried out with the MURaM code, which includes a realistic equation of state with partial ionization and radiative transfer along many ray directions. The largely relaxed state of the sunspot shows a division in a central dark umbral region with bright dots and a penumbra showing bright filaments of about 2 to 3 Mm length with central dark lanes. By a process similar to the formation of umbral dots, the penumbral filaments result from magneto-convection in the form of upflow plumes, which become elongated by the presence of an inclined magnetic field: the upflow is deflected in the outward direction while the magnetic field is weakened and becomes almost horizontal in the upper part of the plume near the level of optical depth unity. A dark lane forms owing to the piling up of matter near the cusp-shaped top of the rising plume that leads to an upward bulging of the surfaces of constant optical depth. The simulated penumbral structure corresponds well to the observationally inferred interlocking-comb structure of the magnetic field with Evershed outflows along dark-laned filaments with nearly horizontal magnetic field and overturning perpendicular (`twisting') motion, which are embedded in a background of stronger and less inclined field. Photospheric spectral lines are formed at the very top and somewhat above the upflow plumes, so that they do not fully sense the strong flow as well as the large field inclination and significant field strength reduction in the upper part of the plume structures.

We present independent observations of the solar-cycle variation of flows near the solar surface and at a depth of about 60 Mm, in the latitude range $\pm45^\circ$. We show that the time-varying components of the meridional flow at these two depths have opposite sign, while the time-varying components of the zonal flow are in phase. This is in agreement with previous results. We then investigate whether the observations are consistent with a theoretical model of solar-cycle dependent meridional circulation based on a flux-transport dynamo combined with a geostrophic flow caused by increased radiative loss in the active region belt (the only existing quantitative model). We find that the model and the data are in qualitative agreement, although the amplitude of the solar-cycle variation of the meridional flow at 60 Mm is underestimated by the model.

Helioseismology has revealed many details of solar differential rotation and also its time variation, known as torsional oscillations. So far there is no generally accepted theoretical explanation for torsional oscillations, even though a close relation to the solar activity cycle is evident. On the theoretical side non-kinematic dynamo models (including the Lorentz force feedback on differential rotation) have been used to explain torsional oscillations. In this paper we use a slightly different approach by forcing torsional oscillations in a mean field differential rotation model. Our aim is not a fully self-consistent model but rather to point out a few general properties of torsional oscillations and their possible origin that are independent from a particular dynamo model. We find that the poleward propagating high latitude branch of the torsional oscillations can be explained as a response of the coupled differential rotation / meridional flow system to periodic forcing in mid-latitudes, of either mechanical (Lorentz force) or thermal nature. The speed of the poleward propagation sets constraints on the value of the turbulent viscosity in the solar convection zone to be less than 3x10^8 m^2/s. We also show that the equatorward propagating low latitude branch is very unlikely a consequence of mechanical forcing (Lorentz force) alone, but rather of thermal origin due to the Taylor-Proudman theorem.

In this paper we discuss the transport of toroidal magnetic field by a weak meridional flow at the base of the convection zone. We utilize the differential rotation and meridional flow model developed by Rempel and incorporate feedback of a purely toroidal magnetic field in two ways: directly through the Lorentz force (magnetic tension) and indirectly through quenching of the turbulent viscosity, which affects the parametrized turbulent angular momentum transport in the model. In the case of direct Lorentz force feedback we find that a meridional flow with an amplitude of around 2 m/s can transport a magnetic field with a strength of 20 to 30 kG. Quenching of turbulent viscosity leads to deflection of the meridional flow from the magnetized region and a significant reduction of the transport velocity if the magnetic field is above equipartition strength.

We present a mean field model based on the approach taken by Rempel (astro-ph/0604451) in order to investigate the influence of stochastic fluctuations in the Reynolds stresses on meridional flow and differential rotation. The stochastic fluctuations found in the meridional flow pattern directly resemble the stochastic fluctuations of the Reynolds stresses, while the stochastic fluctuations in the differential rotation are smaller by almost two orders of magnitude. It is further found that the correlation length and time scale of the stochastic fluctuations have only a weak influence on meridional flow, but a significant influence on the magnitude of variations in the differential rotation. We analyze the energy fluxes within the model to estimate time scales for the replenishment of differential rotation and meridional flow. We find that the time scale for the replenishment of differential rotation (~10 years) is nearly four orders of magnitude longer than the time scale for the replenishment of meridional flow, which explains the differences in the response to stochastic fluctuations of the Reynolds stress found for both flow fields.

In this paper we discuss a dynamic flux-transport dynamo model that includes the feedback of the induced magnetic field on differential rotation and meridional flow. We consider two different approaches for the feedback: meanfield Lorentz force and quenching of transport coefficients such as turbulent viscosity and heat conductivity. We find that even strong feedback on the meridional flow does not change the character of the flux-transport dynamo significantly; however it leads to a significant reduction of differential rotation. To a large degree independent from the dynamo parameters, the saturation takes place when the toroidal field at the base of the convection zone reaches between 1.2 an 1.5 T, the energy converted intomagnetic energy corresponds to about 0.1 to 0.2% of the solar luminosity. The torsional oscillations produced through Lorentz force feedback on differential rotation show a dominant poleward propagating branch with the correct phase relation to the magnetic cycle. We show that incorporating enhanced surface cooling of the active region belt (as proposed by Spruit) leads to an equatorward propagating branch in good agreement with observations.

We present a simple model for the solar differential rotation and meridional circulation based on a mean field parameterization of the Reynolds stresses that drive the differential rotation. We include the subadiabatic part of the tachocline and show that this, in conjunction with turbulent heat conductivity within the convection zone and overshoot region, provides the key physics to break the Taylor-Proudman constraint, which dictates differential rotation with contour lines parallel to the axis of rotation in case of an isentropic stratification. We show that differential rotation with contour lines inclined by 10 - 30 degrees with respect to the axis of rotation is a robust result of the model, which does not depend on the details of the Reynolds stress and the assumed viscosity, as long as the Reynolds stress transports angular momentum toward the equator. The meridional flow is more sensitive with respect to the details of the assumed Reynolds stress, but a flow cell, equatorward at the base of the convection zone and poleward in the upper half of the convection zone, is the preferred flow pattern.

The uncombed penumbral model explains the structure of the sunspot penumbra in terms of thick magnetic fibrils embedded in a magnetic surrounding atmosphere. This model has been successfully applied to explain the polarization signals emerging from the sunspot penumbra. Thick penumbral fibrils face some physical problems, however. In this contribution we will offer possible solutions to these shortcomings.

After a dynamically active emergence phase, magnetic flux at the solar surface soon ceases to show strong signs of the subsurface dynamics of its parent magnetic structure. This indicates that some kind of disconnection of the emerged flux from its roots in the deep convection zone should take place. We propose a mechanism for the dynamical disconnection of the surface flux based upon the buoyant upflow of plasma along the field lines. Such flows arise in the upper part of a rising flux loop during the final phases of its buoyant ascent towards the surface. The combination of the pressure buildup by the upflow and the cooling of the upper layers of an emerged flux tube by radiative losses at the surface lead to a progressive weakening of the magnetic field in several Mm depth. When the field strength has become sufficiently low, convective motions and the fluting instability disrupt the flux tube into thin, passively advected flux fragments, thus providing a dynamical disconnection of the emerged part from its roots. We substantiate this scenario by considering the quasi-static evolution of a sunspot model under the effects of radiative cooling, convective energy transport, and pressure buildup by a prescribed inflow at the bottom of the model. For inflow speeds in the range shown by simulations of thin flux tubes, we find that the disconnection takes place in a depth between 2 and 6 Mm for disconnection times up to 3 days.

We show that if the turbulent magnetic diffusivity used in solar dynamos is assumed to be 'quenched' by increasing toroidal fields, much larger amplitude and more concentrated toroidal fields can be induced by differential rotation from an assumed poloidal field than if there is no quenching. This amplification and concentration mechanism is weakened and bounded by j x B feedbacks on the differential rotation. Nevertheless, it is strong enough to contribute to the creation of ~100 kG toroidal fields near the base of the convection zone, perhaps in conjunction with the 'exploding flux tube' process. Such high fields are necessary for sunspots to occur in low solar latitudes.