Search

The length of the solar activity cycle fluctuates considerably. The temporal evolution of the corresponding cycle phase, that is, the deviation of the epochs of activity minima or maxima from strict periodicity, provides relevant information concerning the physical mechanism underlying the cyclic magnetic activity. An underlying strictly periodic process (akin to a perfect "clock"), with the observer seeing a superposition of the perfect clock and a small random phase perturbation, leads to long-term phase stability in the observations. Such behavior would be expected if cycles were synchronized by tides caused by orbiting planets or by a hypothetical torsional oscillation in the solar radiative interior. Alternatively, in the absence of such synchronization, phase fluctuations accumulate and a random walk of the phase ensues, which is a typical property of randomly perturbed dynamo models. Based on the sunspot record and the reconstruction of solar cycles from cosmogenic C14, we carried out rigorous statistical tests in order to decipher whether there exists phase synchronization or random walk. Synchronization is rejected at significance levels of between 95% (28 cycles from sunspot data) and beyond 99% (84 cycles reconstructed from C14, while the existence of random walk in the phases is consistent with all data sets. This result strongly supports randomly perturbed dynamo models with little inter-cycle memory., Comment: Astronomy & Astrophysics, in press

The polarity of the toroidal magnetic field in the solar convection zone periodically reverses in the course of the 11/22-year solar cycle. Among the various processes that contribute to the removal of `old-polarity' toroidal magnetic flux is the emergence of flux loops forming bipolar regions at the solar surface. We quantify the loss of subsurface net toroidal flux by this process. To this end, we determine the contribution of an individual emerging bipolar loop and show that it is unaffected by surface flux transport after emergence. Together with the linearity of the diffusion process this means that the total flux loss can be obtained by adding the contributions of all emerging bipolar magnetic regions. The resulting total loss rate of net toroidal flux amounts to 1.3 X 10^{15} Mx/s during activity maxima and 6.1 X 10^{14} Mx/s during activity minima, to which ephemeral regions contribute about 90% and 97%, respectively. This rate is consistent with the observationally inferred loss rate of toroidal flux into interplanetary space and corresponds to a decay time of the subsurface toroidal flux of about 12 years, also consistent with a simple estimate based on turbulent diffusivity. Consequently, toroidal flux loss by flux emergence is a relevant contribution to the budget of net toroidal flux in the solar convection zone. That the toroidal flux loss rate due to flux emergence is consistent with what is expected from turbulent diffusion, and that the corresponding decay time is similar to the length of the solar cycle are important constraints for understanding the solar cycle and the Sun's internal dynamics., Comment: A&A accepted

The frequency spectrum of the hemispheric asymmetry of solar activity shows enhanced power for the period ranges around 8.5 years and between 30 and 50 years. This can be understood as the sum and beat periods of the superposition of two dynamo modes: a dipolar mode with a (magnetic) period of about 22 years and aquadrupolar mode with a period between 13 and 15 years. An updated Babcock-Leighton-type dynamo model with weak driving as indicated by stellar observations shows an excited dipole mode and a damped quadrupole mode in the correct range of periods. Random excitation of the quadrupole by stochastic fluctuations of the source term for the poloidal field leads to a time evolution of activity and asymmetry that is consistent with the observational results., Comment: Astronomy & Astrophysics, accepted

Context. The solar dynamo consists of a process that converts poloidal field to toroidal field followed by a process which creates new poloidal field from the toroidal field. Aims. Our aim is to observe the poloidal and toroidal fields relevant to the global solar dynamo and see if their evolution is captured by a Babcock-Leighton dynamo. Methods. We use synoptic maps of the surface radial field from the KPNSO/VT and SOLIS observatories to construct the poloidal field as a function of time and latitude, and Wilcox Solar Observatory and SOHO/MDI full disk images to infer the longitudinally averaged surface azimuthal field. We show that the latter is consistent with an estimate of that due to flux emergence and therefore closely related to the subsurface toroidal field. Results. We present maps of the poloidal and toroidal magnetic field of the global solar dynamo. The longitude-averaged azimuthal field observed at the surface results from flux emergence. At high latitudes this component follows the radial component of the polar fields with a short time lag (1-3 years). The lag increases at lower latitudes. The observed evolution of the poloidal and toroidal magnetic fields is described by the (updated) Babcock-Leighton dynamo model., Comment: A&A

In 1969 Leighton developed a quasi-1D mathematical model of the solar dynamo, building upon the phenomenological scenario of Babcock(1961). Here we present a modification and extension of Leighton's model. Using the axisymmetric component of the magnetic field, we consider the radial field component at the solar surface and the radially integrated toroidal magnetic flux in the convection zone, both as functions of latitude. No assumptions are made with regard to the radial location of the toroidal flux. The model includes the effects of turbulent diffusion at the surface and in the convection zone, poleward meridional flow at the surface and an equatorward return flow affecting the toroidal flux, latitudinal differential rotation and the near-surface layer of radial rotational shear, downward convective pumping of magnetic flux in the shear layer, and flux emergence in the form of tilted bipolar magnetic regions. While the parameters relevant for the transport of the surface field are taken from observations, the model condenses the unknown properties of magnetic field and flow in the convection zone into a few free parameters (turbulent diffusivity, effective return flow, amplitude of the source term, and a parameter describing the effective radial shear). Comparison with the results of two-dimensional flux transport dynamo codes shows that the model captures the essential features of these simulations. We carry out a parameter study over the four-dimensional parameter space and identify the parameter ranges that provide solar-like solutions. Dipole parity is always preferred and solutions with periods around 22 years and a correct phase difference between flux emergence in low latitudes and the strength of the polar fields are found for a return flow speed around 2~m/s, turbulent diffusivity <80~km^2/s, and dynamo excitation not too far above the threshold (linear growth rate <0.1/yr)., Comment: A&A accepted

In order to match observed properties of the solar cycle, flux-transport dynamo models require the toroidal magnetic flux to be stored in a region of low magnetic diffusivity, typically located at or below the bottom of the convection zone. We infer the turbulent magnetic diffusivity affecting the toroidal field on the basis of empirical data. We consider the time evolution of mean latitude and width of the activity belts of solar cycles 12--23 and their dependence on cycle strength. We interpret the decline phase of the cycles as a diffusion process. The activity level of a given cycle begins to decline when the centers of its equatorward propagating activity belts come within their width (at half maximum) from the equator. This happens earlier for stronger cycles because their activity belts are wider. From that moment on, the activity and the belt width decrease in the same manner for all cycles, independent of their maximum activity level. In terms of diffusive cancellation of opposite-polarity toroidal flux across the equator, we infer the turbulent diffusivity experienced by the toroidal field, wherever it is located, to be in the range 150--450~km$^2$s$^{-1}$. Strong diffusive latitudinal spreading of the toroidal flux underneath the activity belts can be inhibited by an inflow towards the toroidal field bands in the convection zone with a magnitude of several meters per second. The inferred value of the turbulent magnetic diffusivity affecting the toroidal field agrees, to order of magnitude, with estimates based on mixing-length models for the solar convection zone. This is at variance with the requirement of flux-transport dynamo models. The inflows required to keep the toroidal field bands together before they approach the equator are similar to the inflows towards the activity belts observed with local helioseismology., Comment: A&A, accepted

A number of observational and theoretical aspects of solar magnetoconvection are considered in this review. We discuss recent developments in our understanding of the small-scale structure of the magnetic field on the solar surface and its interaction with convective flows, which is at the centre of current research. Topics range from plage areas in active regions over the magnetic network shaped by supergranulation to the ubiquituous `turbulent' internetwork fields. On the theoretical side, we focus upon magnetic field generation by small-scale dynamo action., Comment: Space Science Reviews, Springer, in press

We discuss recent observational, theoretical and modeling progress made in understanding the Sun's internal dynamics, including its rotation, meridional flow, convection and overshoot. Over the past few decades, substantial theoretical and observational effort has gone into appreciating these aspects of solar dynamics. A review of these observations, related helioseismic methodology and inference and computational results in relation to these problems is undertaken here., Comment: 31 pages, 10 figures, Space Science Reviews

Magnetic activity in Sun-like and low-mass stars causes X-ray coronal emission, which is stronger for more rapidly rotating stars. This relation is often interpreted in terms of the Rossby number, i.e., the ratio of rotation period to convective overturn time. We reconsider this interpretation on the basis of the observed X-ray emission and rotation periods of 821 stars with masses below 1.4 Msun. A generalized analysis of the relation between X-ray luminosity normalized by bolometric luminosity, L_X/L_bol, and combinations of rotational period, P, and stellar radius, R, shows that the Rossby formulation does not provide the solution with minimal scatter. Instead, we find that the relation L_X/L_bol ~ P^{-2}R^{-4} optimally describes the non-saturated fraction of the stars. This relation is equivalent to L_X ~ P^{-2}, indicating that the rotation period alone determines the total X-ray emission. Since L_X is directly related to the magnetic flux at the stellar surface, this means that the surface flux is determined solely by the star's rotation and is independent of other stellar parameters. While a formulation in terms of a Rossby number would be consistent with these results if the convective overturn time scales exactly as L_bol^{-1/2}, our generalized approach emphasizes the need to test a broader range of mechanisms for dynamo action in cool stars., Comment: 9 pages, accepted for publication in The Astrophysical Journal

The tilt angles of sunspot groups represent the poloidal field source in Babcock-Leighton-type models of the solar dynamo and are crucial for the build-up and reversals of the polar fields in Surface Flux Transport (SFT) simulations. The evolution of the polar field is a consequence of Hale's polarity rules, together with the tilt angle distribution which has a systematic component (Joy's law) and a random component (tilt-angle scatter). We determine the scatter using the observed tilt angle data and study the effects of this scatter on the evolution of the solar surface field using SFT simulations with flux input based upon the recorded sunspot groups. The tilt angle scatter is described in our simulations by a random component according to the observed distributions for different ranges of sunspot group size (total umbral area). By performing simulations with a number of different realizations of the scatter we study the effect of the tilt angle scatter on the global magnetic field, especially on the evolution of the axial dipole moment. The average axial dipole moment at the end of cycle 17 (a medium-amplitude cycle) from our simulations was 2.73G. The tilt angle scatter leads to an uncertainty of 0.78 G (standard deviation). We also considered cycle 14 (a weak cycle) and cycle 19 (a strong cycle) and show that the standard deviation of the axial dipole moment is similar for all three cycles. The uncertainty mainly results from the big sunspot groups which emerge near the equator. In the framework of Babcock-Leighton dynamo models, the tilt angle scatter therefore constitutes a significant random factor in the cycle-to-cycle amplitude variability, which strongly limits the predictability of solar activity., Comment: 23 pages, 7 figures, accepted for publication in ApJ

Bright points (BPs) in the solar photosphere are radiative signatures of magnetic elements described by slender flux tubes located in the darker intergranular lanes. They contribute to the ultraviolet (UV) flux variations over the solar cycle and hence may influence the Earth's climate. Here we combine high-resolution UV and spectro-polarimetric observations of BPs by the SUNRISE observatory with 3D radiation MHD simulations. Full spectral line syntheses are performed with the MHD data and a careful degradation is applied to take into account all relevant instrumental effects of the observations. It is demonstrated that the MHD simulations reproduce the measured distributions of intensity at multiple wavelengths, line-of-sight velocity, spectral line width, and polarization degree rather well. Furthermore, the properties of observed BPs are compared with synthetic ones. These match also relatively well, except that the observations display a tail of large and strongly polarized BPs not found in the simulations. The higher spatial resolution of the simulations has a significant effect, leading to smaller and more numerous BPs. The observation that most BPs are weakly polarized is explained mainly by the spatial degradation, the stray light contamination, and the temperature sensitivity of the Fe I line at 5250.2 \AA{}. The Stokes $V$ asymmetries of the BPs increase with the distance to their center in both observations and simulations, consistent with the classical picture of a production of the asymmetry in the canopy. This is the first time that this has been found also in the internetwork. Almost vertical kilo-Gauss fields are found for 98 % of the synthetic BPs. At the continuum formation height, the simulated BPs are on average 190 K hotter than the mean quiet Sun, their mean BP field strength is 1750 G, supporting the flux-tube paradigm to describe BPs., Comment: Accepted for publication in Astronomy & Astrophysics on May 30 2014

Within the Babcock-Leighton framework for the solar dynamo, the strength of a cycle is expected to depend on the strength of the dipole moment or net hemispheric flux during the preceding minimum, which depends on how much flux was present in each hemisphere at the start of the previous cycle and how much net magnetic flux was transported across the equator during the cycle. Some of this transport is associated with the random walk of magnetic flux tubes subject to granular and supergranular buffeting, some of it is due to the advection caused by systematic cross-equatorial flows such as those associated with the inflows into active regions, and some crosses the equator during the emergence process. We aim to determine how much of the cross-equatorial transport is due to small-scale disorganized motions (treated as diffusion) compared with other processes such as emergence flux across the equator. We measure the cross-equatorial flux transport using Kitt Peak synoptic magnetograms, estimating both the total and diffusive fluxes. Occasionally a large sunspot group, with a large tilt angle emerges crossing the equator, with flux from the two polarities in opposite hemispheres. The largest of these events carry a substantial amount of flux across the equator (compared to the magnetic flux near the poles). We call such events cross-equatorial flux plumes. There are very few such large events during a cycle, which introduces an uncertainty into the determination of the amount of magnetic flux transported across the equator in any particular cycle. As the amount of flux which crosses the equator determines the amount of net flux in each hemisphere, it follows that the cross-equatorial plumes introduce an uncertainty in the prediction of the net flux in each hemisphere. This leads to an uncertainty in predictions of the strength of the following cycle., Comment: A&A, accepted

It is proposed that the observed near-surface inflows towards the active regions and sunspot zones provide a nonlinear feedback mechanism that limits the amplitude of a Babcock-Leighton-type solar dynamo and determines the variation of the cycle strength. This hypothesis is tested with surface flux transport simulations including converging latitudinal flows that depend on the surface distribution of magnetic flux. The inflows modulate the build-up of polar fields (represented by the axial dipole) by reducing the tilt angles of bipolar magnetic regions and by affecting the cross-equator transport of leading-polarity magnetic flux. With flux input derived from the observed record of sunspot groups, the simulations cover the period between 1874 and 1980 (corresponding to solar cycles 11 to 20). The inclusion of the inflows leads to a strong correlation of the simulated axial dipole strength during activity minimum with the observed amplitude of the subsequent cycle. This in agreement with empirical correlations and in line with what is expected from a Babcock-Leighton-type dynamo. The results provide evidence that the latitudinal inflows are a key ingredient in determining the amplitude of solar cycles., Comment: accepted in A&A

Aims: We study the differences between non-magnetic and magnetic regions in the flow and thermal structure of the upper solar photosphere. Methods: Radiative MHD simulations representing a quiet region and a plage region, respectively, which extend into the layers around the temperature minimum, are analyzed. Results: The flow structure in the upper photospheric layers of the two simulations is considerably different: the non-magnetic simulation is dominated by a pattern of moving shock fronts while the magnetic simulation shows vertically extended vortices associated with magnetic flux concentrations. Both kinds of structures induce substantial local heating. The resulting average temperature profiles are characterized by a steep rise above the temperature minimum due to shock heating in the non-magnetic case and by a flat photospheric temperature gradient mainly caused by Ohmic dissipation in the magnetic run. Conclusions: Shocks in the quiet Sun and vortices in the strongly magnetized regions represent the dominant flow structures in the layers around the temperature minimum. They are closely connected with dissipation processes providing localized heating., Comment: Accepted for publicaton in A&A

Radiative hydrodynamic simulations of solar and stellar surface convection have become an important tool for exploring the structure and gas dynamics in the envelopes and atmospheres of late-type stars and for improving our understanding of the formation of stellar spectra. We quantitatively compare results from three-dimensional, radiative hydrodynamic simulations of convection near the solar surface generated with three numerical codes CO5BOLD, MURaM, and STAGGER and different simulation setups in order to investigate the level of similarity and to cross-validate the simulations. For all three simulations, we considered the average stratifications of various quantities (temperature, pressure, flow velocity, etc.) on surfaces of constant geometrical or optical depth, as well as their temporal and spatial fluctuations. We also compared observables, such as the spatially resolved patterns of the emerging intensity and of the vertical velocity at the solar optical surface as well as the center-to-limb variation of the continuum intensity at various wavelengths. The depth profiles of the thermodynamical quantities and of the convective velocities as well as their spatial fluctuations agree quite well. Slight deviations can be understood in terms of differences in box size, spatial resolution and in the treatment of non-gray radiative transfer between the simulations. The results give confidence in the reliability of the results from comprehensive radiative hydrodynamic simulations., Comment: Astronomy & Astrophysics, in press

We report on the occurrence of small-scale vortices in simulations of the convective solar surface. Using an eigenanalysis of the velocity gradient tensor, we find the subset of high vorticity regions in which the plasma is swirling. The swirling regions form an unsteady, tangled network of filaments in the turbulent downflow lanes. Near-surface vertical vortices are underdense and cause a local depression of the optical surface. They are potentially observable as bright points in the dark intergranular lanes. Vortex features typically exist for a few minutes, during which they are moved and twisted by the motion of the ambient plasma. The bigger vortices found in the simulations are possibly, but not necessarily, related to observations of granular-scale spiraling pathlines in "cork animations" or feature tracking., Comment: 11 pages, 13 figures, accepted for publication in A&A, complementary movies at http://www.mps.mpg.de/homes/moll/strudel/papermovies/

We quantify possible differences between turbulent dynamo action in the Sun and the dynamo action studied in idealized simulations. For this purpose we compare Fourier-space shell-to-shell energy transfer rates of three incrementally more complex dynamo simulations: an incompressible, periodic simulation driven by random flow, a simulation of Boussinesq convection, and a simulation of fully compressible convection that includes physics relevant to the near-surface layers of the Sun. For each of the simulations studied, we find that the dynamo mechanism is universal in the kinematic regime because energy is transferred from the turbulent flow to the magnetic field from wavenumbers in the inertial range of the energy spectrum. The addition of physical effects relevant to the solar near-surface layers, including stratification, compressibility, partial ionization, and radiative energy transport, does not appear to affect the nature of the dynamo mechanism. The role of inertial-range shear stresses in magnetic field amplification is independent from outer-scale circumstances, including forcing and stratification. Although the shell-to-shell energy transfer functions have similar properties to those seen in mean-flow driven dynamos in each simulation studied, the saturated states of these simulations are not universal because the flow at the driving wavenumbers is a significant source of energy for the magnetic field., Comment: 16 pages, 9 figures, accepted for publication in ApJ

We have used semi-synthetic records of emerging sunspot groups based on sunspot number data as input for a surface flux transport model to reconstruct the evolution of the large-scale solar magnetic field and the open heliospheric flux from the year 1700 onward. The statistical properties of the semi-synthetic sunspot group records reflect those of the observed the Royal Greenwich Observatory photoheliographic results. These include correlations between the sunspot numbers and sunspot group latitudes, longitudes, areas and tilt angles. The reconstruction results for the total surface flux, the polar field, and the heliospheric open flux (determined by a current sheet source surface extrapolation) agree well with the available observational or empirically derived data and reconstructions. We confirm a significant positive correlation between the polar field during activity minimum periods and the strength of the subsequent sunspot cycle, which has implications for flux transport dynamo models for the solar cycle. Just prior to the Dalton minimum, at the end of the 18th century, a long cycle was followed by a weak cycle. We find that introducing a possibly `lost' cycle between 1793 and 1800 leads to a shift of the minimum of the open flux by 15 years which is inconsistent with the cosmogenic isotope record., Comment: A&A, accepted

Changes in the solar surface area covered by small-scale magnetic elements are thought to cause long-term changes in the solar spectral irradiance, which are important for determining the impact on Earth's climate. To study the effect of small-scale magnetic elements on total and spectral irradiance, we derive their contrasts from 3-D MHD simulations of the solar atmosphere. Such calculations are necessary since measurements of small-scale flux tube contrasts are confined to a few wavelengths and suffer from scattered light and instrument defocus, even for space observations. To test the contrast calculations, we compare rms contrasts from simulations with those obtained with the broad-band filter imager mounted on the Solar Optical Telescope (SOT) onboard the Hinode satellite and also analyse centre-to-limb variations (CLV). The 3-D MHD simulations include the interaction between convection and magnetic flux tubes. They have been run with non-grey radiative transfer using the MURaM code. Simulations have an average vertical magnetic field of 0G, 50G, and 200G. Emergent intensities are calculated with the spectral synthesis code ATLAS9 and are convolved with a theoretical point-spread function to account for the properties of the observations' optical system. We find reasonable agreement for simulated and observed intensity distributions in the visible continuum bands. Agreement is poorer for the CN and G-Bands. The analysis of the simulations exhibits a potentially more realistic centre-to-limb behaviour than calculations based on 1-D model atmospheres. We conclude that starting from 3-D MHD simulations represents a powerful approach to obtaining intensity contrasts for a wide wavelength coverage and for different positions on the solar disk. This also paves the way for future calculations of facular and network contrast as a function of magnetic fluxes., Comment: 9 pages, 8 figures

Bright points (BPs) are manifestations of small magnetic elements in the solar photosphere. Their brightness contrast not only gives insight into the thermal state of the photosphere (and chromosphere) in magnetic elements, but also plays an important role in modulating the solar total and spectral irradiance. Here we report on simultaneous high-resolution imaging and spectropolarimetric observations of BPs using Sunrise balloon-borne observatory data of the quiet Sun at disk center. BP contrasts have been measured between 214 nm and 525 nm, including the first measurements at wavelengths below 388 nm. The histograms of the BP peak brightness show a clear trend toward broader contrast distributions and higher mean contrasts at shorter wavelengths. At 214 nm we observe a peak brightness of up to five times the mean quiet-Sun value, the highest BP contrast so far observed. All BPs are associated with a magnetic signal, although in a number of cases it is surprisingly weak. Most of the BPs show only weak downflows, the mean value being 240 m/s, but some display strong down- or upflows reaching a few km/s., Comment: Accepted for publication in The Astrophysical Journal Letters on September 08 2010

Until today, the small size of magnetic elements in quiet Sun areas has required the application of indirect methods, such as the line-ratio technique or multi-component inversions, to infer their physical properties. A consistent match to the observed Stokes profiles could only be obtained by introducing a magnetic filling factor that specifies the fraction of the observed pixel filled with magnetic field. Here, we investigate the properties of a small magnetic patch in the quiet Sun observed with the IMaX magnetograph on board the balloon-borne telescope Sunrise with unprecedented spatial resolution and low instrumental stray light. We apply an inversion technique based on the numerical solution of the radiative transfer equation to retrieve the temperature stratification and the field strength in the magnetic patch. The observations can be well reproduced with a one-component, fully magnetized atmosphere with a field strength exceeding 1 kG and a significantly enhanced temperature in the mid- to upper photosphere with respect to its surroundings, consistent with semi-empirical flux tube models for plage regions. We therefore conclude that, within the framework of a simple atmospheric model, the IMaX measurements resolve the observed quiet-Sun flux tube., Comment: Accepted for publication in The Astrophysical Journal Letters on Aug 11 2010

The Sunrise balloon-borne solar observatory consists of a 1m aperture Gregory telescope, a UV filter imager, an imaging vector polarimeter, an image stabilization system and further infrastructure. The first science flight of Sunrise yielded high-quality data that reveal the structure, dynamics and evolution of solar convection, oscillations and magnetic fields at a resolution of around 100 km in the quiet Sun. After a brief description of instruments and data, first qualitative results are presented. In contrast to earlier observations, we clearly see granulation at 214 nm. Images in Ca II H display narrow, short-lived dark intergranular lanes between the bright edges of granules. The very small-scale, mixed-polarity internetwork fields are found to be highly dynamic. A significant increase in detectable magnetic flux is found after phase-diversity-related reconstruction of polarization maps, indicating that the polarities are mixed right down to the spatial resolution limit, and probably beyond., Comment: accepted by ApJL

We present the first observations of the transverse component of photospheric magnetic field acquired by the imaging magnetograph Sunrise/IMaX. Using an automated detection method, we obtain statistical properties of 4536 features with significant linear polarization signal. Their rate of occurrence is 1-2 orders of magnitude larger than values reported by previous studies. We show that these features have no characteristic size or lifetime. They appear preferentially at granule boundaries with most of them being caught in downflow lanes at some point in their evolution. Only a small percentage are entirely and constantly embedded in upflows (16%) or downflows (8%)., Comment: Accepted for the Sunrise Special Issue of ApJL

The solar meridional flow is an important ingredient in Babcock-Leighton type models of the solar dynamo. Global variations of this flow have been suggested to explain the variations in the amplitudes and lengths of the activity cycles. Recently, cycle-related variations in the amplitude of the $P_2^1$ term in the Legendre decomposition of the observed meridional flow have been reported. The result is often interpreted in terms of an overall variation in the flow amplitude during the activity cycle. Using a semi-empirical model based upon the observed distribution of magnetic flux on the solar surface, we show that the reported variations of the $P_2^1$ term can be explained by the observed localized inflows into the active region belts. No variation of the overall meridional flow amplitude is required., Comment: Accepted in ApJ

We analyse mesogranular flow patterns in a three-dimensional hydrodynamical simulation of solar surface convection in order to determine its characteristics. We calculate divergence maps from horizontal velocities obtained with the Local Correlation Tracking (LCT) method. Mesogranules are identified as patches of positive velocity divergence. We track the mesogranules to obtain their size and lifetime distributions. We vary the analysis parameters to verify if the pattern has characteristic scales. The characteristics of the resulting flow patterns depend on the averaging time and length used in the analysis. We conclude that the mesogranular patterns do not exhibit intrinsic length and time scales.

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., Comment: Accepted for publication in ApJ

We model the surface magnetic field and open flux of the Sun from 1913 to 1986 using a surface flux transport model, which includes the observed cycle-to-cycle variation of sunspot group tilts. The model reproduces the empirically derived time evolution of the solar open magnetic flux, and the reversal times of the polar fields. We find that both the polar field and the axial dipole moment resulting from this model around cycle minimum correlate with the strength of the following cycle., Comment: Accepted for publication by ApJ

We studied the effect of the perturbation of the meridional flow in the activity belts detected by local helioseismology on the development and strength of the surface magnetic field at the polar caps. We carried out simulations of synthetic solar cycles with a flux transport model, which follows the cyclic evolution of the surface field determined by flux emergence and advective transport by near-surface flows. In each hemisphere, an axisymmetric band of latitudinal flows converging towards the central latitude of the activity belt was superposed onto the background poleward meridional flow. The overall effect of the flow perturbation is to reduce the latitude separation of the magnetic polarities of a bipolar magnetic region and thus diminish its contribution to the polar field. As a result, the polar field maximum reached around cycle activity minimum is weakened by the presence of the meridional flow perturbation. For a flow perturbation consistent with helioseismic observations, the polar field is reduced by about 18% compared to the case without inflows. If the amplitude of the flow perturbation depends on the cycle strength, its effect on the polar field provides a nonlinearity that could contribute to limiting the amplitude of a Babcock-Leighton type dynamo., Comment: 16 pages, 5 figures, accepted by ApJ

We obtain information about the magnetic flux present in the quiet Sun by comparing radiative MHD simulations with Hinode/SP observations, with particular emphasis on the role of surface dynamo action. Simulation runs with different magnetic Reynolds numbers (Rm) are used together with observations at different heliocentric angles with different levels of noise. The results show that simulations with an imposed mixed-polarity field and Rm below the threshold for dynamo action reproduce the observed vertical flux density, but do not display a sufficiently high horizontal flux density. Surface dynamo simulations at the highest Rm feasible at the moment yield a ratio of the horizontal and vertical flux density consistent with observational results, but the overall amplitudes are too low. Based on the properties of the local dynamo simulations, a tentative scaling of the magnetic field strength by a factor 2 - 3 reproduces the signal observed in the internetwork regions. We find an agreement with observations at different heliocentric angles. The mean field strength in internetwork, implied by our analysis, is roughly 170 G at the optical depth unity. Our study shows that surface dynamo could be responsible for most of the magnetic flux in the quiet Sun outside the network given that the extrapolation to higher Rm is valid., Comment: accepted in A&A

By coupling a solar surface flux transport model with an extrapolation of the heliospheric field, we simulate the evolution of the Sun's open magnetic flux and the heliospheric current sheet (HCS) based on observational data of sunspot groups since 1976. The results are consistent with measurements of the interplanetary magnetic field near Earth and with the tilt angle of the HCS as derived from extrapolation of the observed solar surface field. This opens the possibility for an improved reconstruction of the Sun's open flux and the HCS into the past on the basis of empirical sunspot data., Comment: 16 pages, 5 figures, Accepted for publication in ApJ

Realistic 3D radiative MHD simulations reveal the magneto-convective processes underlying the formation of the photospheric fine structure of sunspots, including penumbral filaments and umbral dots. Here we provide results from a statistical analysis of simulated umbral dots and compare them with reports from high-resolution observations. A multi-level segmentation and tracking algorithm has been used to isolate the bright structures in synthetic bolometric and continuum brightness images. Areas, brightness, and lifetimes of the resulting set of umbral dots are found to be correlated: larger umbral dots tend to be brighter and live longer. The magnetic field strength and velocity structure of umbral dots on surfaces of constant optical depth in the continuum at 630 nm indicate that the strong field reduction and high velocities in the upper parts of the upflow plumes underlying umbral dots are largely hidden from spectro-polarimetric observations. The properties of the simulated umbral dots are generally consistent with the results of recent high-resolution observations. However, the observed population of small, short-lived umbral dots is not reproduced by the simulations, possibly owing to insufficient spatial resolution., Comment: Accepted for publication in A&A

We study the observational signature of flux emergence in the photosphere using synthetic data from a 3D MHD simulation of the emergence of a twisted flux tube. Several stages in the emergence process are considered. At every stage we compute synthetic Stokes spectra of the two iron lines Fe I 6301.5 {\AA} and Fe I 6302.5 {\AA} and degrade the data to the spatial and spectral resolution of Hinode's SOT/SP. Then, following observational practice, we apply Milne-Eddington-type inversions to the synthetic spectra in order to retrieve various atmospheric parameters and compare the results with recent Hinode observations. During the emergence sequence, the spectral lines sample different parts of the rising flux tube, revealing its twisted structure. The horizontal component of the magnetic field retrieved from the simulations is close to the observed values. The flattening of the flux tube in the photosphere is caused by radiative cooling, which slows down the ascent of the tube to the upper solar atmosphere. Consistent with the observations, the rising magnetized plasma produces a blue shift of the spectral lines during a large part of the emergence sequence., Comment: A&A Letter, 3 figures

The structure and dynamics of small vertical photospheric magnetic flux concentrations has been often treated in the framework of an approximation based upon a low-order truncation of the Taylor expansions of all quantities in the horizontal direction, together with the assumption of instantaneous total pressure balance at the boundary to the non-magnetic external medium. Formally, such an approximation is justified if the diameter of the structure (a flux tube or a flux sheet) is small compared to all other relevant length scales (scale height, radius of curvature, wavelength, etc.). The advent of realistic 3D radiative MHD simulations opens the possibility of checking the consistency of the approximation with the properties of the flux concentrations that form in the course of a simulation. We carry out a comparative analysis between the thin flux tube/sheet models and flux concentrations formed in a 3D radiation-MHD simulation. We compare the distribution of the vertical and horizontal components of the magnetic field in a 3D MHD simulation with the field distribution in the case of the thin flux tube/sheet approximation. We also consider the total (gas plus magnetic) pressure in the MHD simulation box. Flux concentrations with super-equipartition fields are reasonably well reproduced by the second-order thin flux tube/sheet approximation. The differences between approximation and simulation are due to the asymmetry and the dynamics of the simulated structures.

Recent spectro-polarimetric observations have provided detailed measurements of magnetic field, velocity and intensity during events of magnetic field intensification in the solar photosphere. We consider the temporal evolution of the relevant physical quantities for three cases of magnetic field intensification in a numerical simulation. We determine the evolution of the intensity, magnetic flux density and zero-crossing velocity derived from the synthetic Stokes parameters by taking into account the spectral and spatial resolution of the spectropolarimeter (SP) on board Hinode. The three events considered show a similar evolution: advection of magnetic flux to a granular vertex, development of a strong downflow, evacuation of the magnetic feature, increase of the field strength and the appearance of the bright point. We find that synthetic and real observations are qualitatively consistent and, for one of the cases considered, agree very well also quantitatively. The effect of finite resolution (spatial smearing) is most pronounced in the case of small features, for which the synthetic Hinode/SP observations miss the bright point formation and also the high-velocity downflows during the formation of the smaller magnetic features., Comment: accepted in A&A

Long-term trends in the solar spectral irradiance are important to determine the impact on Earth's climate. These long-term changes are thought to be caused mainly by changes in the surface area covered by small-scale magnetic elements. The direct measurement of the contrast to determine the impact of these small-scale magnetic elements is, however, limited to a few wavelengths, and is, even for space instruments, affected by scattered light and instrument defocus. In this work we calculate emergent intensities from 3-D simulations of solar magneto-convection and validate the outcome by comparing with observations from Hinode/SOT. In this manner we aim to construct the contrast at wavelengths ranging from the NUV to the FIR., Comment: Proceedings paper, IAU XXVII, Symposium 264, 3 pages

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., Comment: 19 pages, 8 figures, main Science article + supporting online material combined into one file

We present results from numerical modeling of emerging flux regions on the solar surface. The modeling was carried out by means of 3D radiative MHD simulations of the rise of buoyant magnetic flux tubes through the convection zone and into the photosphere. Due to the strong stratification of the convection zone, the rise results in a lateral expansion of the tube into a magnetic sheet, which acts as a reservoir for small-scale flux emergence events at the scale of granulation. The interaction of the convective downflows and the rising magnetic flux undulates it to form serpentine field lines emerging into the photosphere. Observational characteristics including the pattern of emerging flux regions, the cancellation of surface flux and associated high speed downflows, the convective collapse of photospheric flux tubes, the appearance of anomalous darkenings, the formation of bright points and the possible existence of transient kilogauss horizontal fields are discussed in the context of new observations from the Hinode Solar Optical Telescope. Implications for the local helioseismology of emerging flux regions are also discussed.

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., Comment: 10 pages, 11 figures. Added references, corrected typos, minor revisions. To be published in ApJ

The contrast of granulation is an important quantity characterizing solar surface convection. We compare the intensity contrast at 630 nm, observed using the Spectro-Polarimeter (SP) aboard the Hinode satellite, with the 3D radiative MHD simulations of V{\"o}gler & Sch{\"u}ssler (2007). A synthetic image from the simulation is degraded using a theoretical point-spread function of the optical system, and by considering other important effects. The telescope aperture and the obscuration by the secondary mirror and its attachment spider, reduce the simulated contrast from 14.4 % to 8.5 %. A slight effective defocus of the instrument brings the simulated contrast down to 7.5 %, close to the observed value of 7.0 %. A proper consideration of the effects of the optical system and a slight defocus, lead to sufficient degradation of the synthetic image from the MHD simulation, such that the contrast reaches almost the observed value. The remaining small discrepancy can be ascribed to straylight and slight imperfections of the instrument, which are difficult to model. Hence, Hinode SP data are consistent with a granulation contrast which is predicted by 3D radiation MHD simulations., Comment: 5 pages, 4 figures, to be published in A&A

Observations with the Hinode spectro-polarimeter have revealed strong horizontal internetwork magnetic fields in the quiet solar photosphere. We aim at interpreting the observations by means of results from numerical simulations. Radiative MHD simulations of dynamo action by near-surface convection are analyzed with respect to the relation between vertical and horizontal magnetic field components. The dynamo-generated fields show a clear dominance of the horizontal field in the height range where the spectral lines used for the observations are formed. The ratio between the averaged horizontal and vertical field components is consistent with the values derived from the observations. This behavior results from the intermittent nature of the dynamo field with polarity mixing on small scales in the surface layers. Our results provide further evidence that local near-surface dynamo action contributes significantly to the solar internetwork fields., Comment: Astronomy & Astrophysics, in press

Flow-induced instabilities are relevant for the storage and dynamics of magnetic fields in stellar convection zones and possibly also in other astrophysical contexts. We continue the study started in the first paper of this series by considering the stability properties of longitudinal flows along magnetic flux tubes. A linear stability analysis was carried out to determine criteria for the onset of instability in the framework of the approximation of thin magnetic flux tubes. In the non-dissipative case, we find Kelvin-Helmholtz instability for flow velocities exceeding a critical speed that depends on the Alfv{\'e}n speed and on the ratio of the internal and external densities. Inclusion of a friction term proportional to the relative transverse velocity leads to a friction-driven instability connected with backward (or negative energy) waves. We discuss the physical nature of this instability. In the case of a stratified external medium, the Kelvin-Helmholtz instability and the friction-driven instability can set in for flow speeds significantly lower than the Alfv{\'e}n speed. Dissipative effects can excite flow-driven instability below the thresholds for the Kelvin-Helmholtz and the undulatory (Parker-type) instabilities. This may be important for magnetic flux storage in stellar convection zones and for the stability of astrophysical jets., Comment: accepted by Astronomy & Astrophysics

We present results of synthetic spectro-polarimetric diagnostics of radiative MHD simulations of solar surface convection with magnetic fields. Stokes profiles of Zeeman-sensitive lines of neutral iron in the visible and infrared spectral ranges emerging from the simulated atmosphere have been calculated in order to study their relation to the relevant physical quantities and compare with observational results. We have analyzed the dependence of the Stokes-I line strength and width as well as of the Stokes-V signal and asymmetries on the magnetic field strength. Furthermore, we have evaluated the correspondence between the actual velocities in the simulation with values determined from the Stokes-I (Doppler shift of the centre of gravity) and Stokes-V profiles (zero-crossing shift). We confirm that the line weakening in strong magnetic fields results from a higher temperature (at equal optical depth) in the magnetic flux concentrations. We also confirm that considerable Stokes-V asymmetries originate in the peripheral parts of strong magnetic flux concentrations, where the line of sight cuts through the magnetopause of the expanding flux concentration into the surrounding convective donwflow., Comment: Astronomy & Astrophysics, in press

Context: Observations indicate that the `quiet' solar photosphere outside active regions contains considerable amounts of magnetic energy and magnetic flux, with mixed polarity on small scales. The origin of this flux is unclear. Aims: We test whether local dynamo action of the near-surface convection (granulation) can generate a significant contribution to the observed magnetic flux. Methods: We have carried out MHD simulations of solar surface convection, including the effects of strong stratification, compressibility, partial ionization, radiative transfer, as well as an open lower boundary. Results: Exponential growth of a weak magnetic seed field (with vanishing net flux through the computational box) is found in a simulation run with a magnetic Reynolds number of about 2600. The magnetic energy approaches saturation at a level of a few percent of the total kinetic energy of the convective motions. Near the visible solar surface, the (unsigned) magnetic flux density reaches at least a value of about 25 G. Conclusions: A realistic flow topology of stratified, compressible, non-helical surface convection without enforced recirculation is capable of turbulent local dynamo action near the solar surface., Comment: accepted by Astronomy&Astrophysics (Letter)

We study the origin of the predictive skill of some methods to forecast the strength of solar activity cycles. A simple flux transport model for the azimuthally averaged radial magnetic field at the solar surface is used, which contains a source term describing the emergence of new flux based on observational sunspot data. We consider the magnetic flux diffusing over the equator as a predictor, since this quantity is directly related to the global dipole field from which a Babcock-Leighton dynamo generates the toroidal field for the next activity cycle. If the source is represented schematically by a narrow activity belt drifting with constant speed over a fixed range of latitudes between activity minima, our predictor shows considerable predictive skill with correlation coefficients up to 0.95 for past cycles. However, the predictive skill is completely lost when the actually observed emergence latitudes are used. This result originates from the fact that the precursor amplitude is determined by the sunspot activity a few years before solar minimum. Since stronger cycles tend to rise faster to their maximum activity (known as the Waldmeier effect), the temporal overlapping of cycles leads to a shift of the minimum epochs that depends on the strength of the following cycle. This information is picked up by precursor methods and also by our flux transport model with a schematic source. Therefore, their predictive skill does not require a memory, i.e., a physical connection between the surface manifestations of subsequent activity cycles., Comment: Astrophys. Journal, in press

We study the structure and reveal the physical nature of the reversed granulation pattern in the solar photosphere by means of 3-dimensional radiative hydrodynamics simulations. We used the MURaM code to obtain a realistic model of the near-surface layers of the convection zone and the photosphere. The pattern of horizontal temperature fluctuations at the base of the photosphere consists of relatively hot granular cells bounded by the cooler intergranular downflow network. With increasing height in the photosphere, the amplitude of the temperature fluctuations diminishes. At a height of z=130-140 km in the photosphere, the pattern of horizontal temperature fluctuations reverses so that granular regions become relatively cool compared to the intergranular network. Detailed analysis of the trajectories of fluid elements through the photosphere reveal that the motion of the fluid is non-adiabatic, owing to strong radiative cooling when approaching the surface of optical depth unity followed by reheating by the radiation field from below. The temperature structure of the photosphere results from the competition between expansion of rising fluid elements and radiative heating. The former acts to lower the temperature of the fluid whereas the latter acts to increase it towards the radiative equilibrium temperature with a net entropy gain. After the fluid overturns and descends towards the convection zone, radiative energy loss again decreases the entropy of the fluid. Radiative heating and cooling of fluid elements that penetrate into the photosphere and overturn do not occur in equal amounts. The imbalance in the cumulative heating and cooling of these fluid elements is responsible for the reversal of temperature fluctuations with respect to height in the photosphere.

Results from a realistic simulation of 3D radiative magneto-convection in a strong background magnetic field corresponding to the conditions in sunspot umbrae are shown. The convective energy transport is dominated by narrow upflow plumes with adjacent downflows, which become almost field-free near the surface layers. The strong external magnetic field forces the plumes to assume a cusp-like shape in their top parts, where the upflowing plasma loses its buoyancy. The resulting bright features in intensity images correspond well (in terms of brightness, size, and lifetime) to the observed umbral dots in the central parts of sunspot umbrae. Most of the simulated umbral dots have a horizontally elongated form with a central dark lane. Above the cusp, most plumes show narrow upflow jets, which are driven by the pressure of the piled-up plasma below. The large velocities and low field strengths in the plumes are effectively screened from spectroscopic observation because the surfaces of equal optical depth are locally elevated, so that spectral lines are largely formed above the cusp. Our simulations demonstrate that nearly field-free upflow plumes and umbral dots are a natural result of convection in a strong, initially monolithic magnetic field., Comment: Accepted by Astrophysical Journal Letters

Customary two-dimensional flux transport models for the evolution of the magnetic field at the solar surface do not account for the radial structure and the volume diffusion of the magnetic field. When considering the long-term evolution of magnetic flux, this omission can lead to an unrealistic long-term memory of the system and to the suppression of polar field reversals. In order to avoid such effects, we propose an extension of the flux transport model by a linear decay term derived consistently on the basis of the eigenmodes of the diffusion operator in a spherical shell. A decay rate for each eigenmode of the system is determined and applied to the corresponding surface part of the mode evolved in the flux transport model. The value of the volume diffusivity associated with this decay term can be estimated to be in the range 50--100 km^2/s by considering the reversals of the polar fields in comparison of flux transport simulations with observations. We show that the decay term prohibits a secular drift of the polar field in the case of cycles of varying strength, like those exhibited by the historical sunspot record., Comment: for further information visit: http://solweb.oma.be/users/baumann/

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., Comment: 11 pages, 5 figures, accepted by A&A

We address the question whether the properties of the observed latitude-time diagram of sunspot occurence (the butterfly diagram) provide evidence for the operation of a flux-transport dynamo, which explains the migration of the sunspot zones and the period of the solar cycle in terms of a deep equatorward meridional flow. We show that the properties of the butterfly diagram are equally well reproduced by a conventional dynamo model with migrating dynamo waves, but without transport of magnetic flux by a flow. These properties seem to be generic for an oscillatory and migratory field of dipole parity and thus do not permit an observational distinction between different dynamo approaches., Comment: 4 pages, 1 figure