Nonlinear dynamo models using quasi-biennial oscillations constrained by sunspot area data

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

Sunspot modeling: From simplified models to radiative MHD simulations

We review our current understanding of sunspots from the scales of their fine structure to their large scale (global) structure including the processes of their formation and decay. Recently, sunspot models have undergone a dramatic change. In the past, several aspects of sunspot structure have been addressed by static MHD models with parametrized energy transport. Models of sunspot fine structure have been relying heavily on strong assumptions about flow and field geometry (e.g., flux-tubes, "gaps", convective rolls), which were motivated in part by the observed filamentary structure of penumbrae or the necessity of explaining the substantial energy transport required to maintain the penumbral brightness. However, none of these models could self-consistently explain all aspects of penumbral structure (energy transport, filamentation, Evershed flow). In recent years, 3D radiative MHD simulations have been advanced dramatically to the point at which models of complete sunspots with sufficient resolution to capture sunspot fine structure are feasible. Here, overturning convection is the central element responsible for energy transport, filamentation leading to fine structure, and the driving of strong outflows. On the larger scale these models are also in the progress of addressing the subsurface structure of sunspots as well as sunspot formation. With this shift in modeling capabilities and the recent advances in high resolution observations, the future research will be guided by comparing observation and theory

Sunspot Modeling: From Simplified Models to Radiative MHD Simulations

We review our current understanding of sunspots from the scales of their fine structure to their large scale (global) structure including the processes of their formation and decay. Recently, sunspot models have undergone a dramatic change. In the past, several aspects of sunspot structure have been addressed by static MHD models with parametrized energy transport. Models of sunspot fine structure have been relying heavily on strong assumptions about flow and field geometry (e.g., flux-tubes, "gaps", convective rolls), which were motivated in part by the observed filamentary structure of penumbrae or the necessity of explaining the substantial energy transport required to maintain the penumbral brightness. However, none of these models could self-consistently explain all aspects of penumbral structure (energy transport, filamentation, Evershed flow). In recent years, 3D radiative MHD simulations have been advanced dramatically to the point at which models of complete sunspots with sufficient resolution to capture sunspot fine structure are feasible. Here overturning convection is the central element responsible for energy transport, filamentation leading to fine-structure and the driving of strong outflows. On the larger scale these models are also in the progress of addressing the subsurface structure of sunspots as well as sunspot formation. With this shift in modeling capabilities and the recent advances in high resolution observations, the future research will be guided by comparing observation and theory

Can overturning motions in penumbral filaments be detected?

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.Comment: 8 pages, 12 figures, accepted for publication in ApJ on 28th Ju

Transport of toroidal magnetic field by the meridional flow at the base of the solar convection zone

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.Comment: 8 pages, 6 figure

Subsurface magnetic field and flow structure of simulated sunspots

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.Comment: 20 pages, 14 figures, 4 animations, accepted for publication in Ap