Source author record

R. H. Cameron

R. H. Cameron appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

32works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

32 published item(s)

preprint2022arXiv

Faculae cancel out on the surfaces of active Suns

Surfaces of the Sun and other cool stars are filled with magnetic fields, which are either seen as dark compact spots or more diffuse bright structures like faculae. Both hamper detection and characterisation of exoplanets, affecting stellar brightness and spectra, as well as transmission spectra. However, the expected facular and spot signals in stellar data are quite different, for instance they have distinct temporal and spectral profiles. Consequently, corrections of stellar data for magnetic activity can greatly benefit from the insight on whether the stellar signal is dominated by spots or faculae. Here, we utilise a surface flux transport model (SFTM) to show that more effective cancellation of diffuse magnetic flux associated with faculae leads to spot area coverages increasing faster with stellar magnetic activity than that by faculae. Our calculations explain the observed dependence between solar spot and facular area coverages and allow its extension to stars more active than the Sun. This extension enables anticipating the properties of stellar signal and its more reliable mitigation, leading to a more accurate characterisation of exoplanets and their atmospheres.

preprint2022arXiv

Impact of spatially correlated fluctuations in sunspots on metrics related to magnetic twist

The twist of the magnetic field above a sunspot is an important quantity in solar physics. For example, magnetic twist plays a role in the initiation of flares and coronal mass ejections (CMEs). Various proxies for the twist above the photosphere have been found using models of uniformly twisted flux tubes, and are routinely computed from single photospheric vector magnetograms. One class of proxies is based on $α_z$, the ratio of the vertical current to the vertical magnetic field. Another class of proxies is based on the so-called twist density, $q$, which depends on the ratio of the azimuthal field to the vertical field. However, the sensitivity of these proxies to temporal fluctuations of the magnetic field has not yet been well characterized. We aim to determine the sensitivity of twist proxies to temporal fluctuations in the magnetic field as estimated from time-series of SDO/HMI vector magnetic field maps. To this end, we introduce a model of a sunspot with a peak vertical field of 2370 Gauss at the photosphere and a uniform twist density $q= -0.024$ Mm$^{-1}$. We add realizations of the temporal fluctuations of the magnetic field that are consistent with SDO/HMI observations, including the spatial correlations. Using a Monte-Carlo approach, we determine the robustness of the different proxies to the temporal fluctuations. The temporal fluctuations of the three components of the magnetic field are correlated for spatial separations up to 1.4 Mm (more than expected from the point spread function alone). The Monte-Carlo approach enables us to demonstrate that several proxies for the twist of the magnetic field are not biased in each of the individual magnetograms. The associated random errors on the proxies have standard deviations in the range between $0.002$ and $0.006$ Mm$^{-1}$, which is smaller by approximately one order of magnitude than the mean value of $q$.

preprint2022arXiv

The crucial role of surface magnetic fields for stellar dynamos: Epsilon Eridani, 61 Cygni A, and the Sun

Cool main-sequence stars, such as the Sun, have magnetic fields which are generated by an internal dynamo mechanism. In the Sun, the dynamo mechanism produces a balance between the amounts of magnetic flux generated and lost over the Sun's 11-year activity cycle and it is visible in the Sun's different atmospheric layers using multi-wavelength observations. We used the same observational diagnostics, spanning several decades, to probe the emergence of magnetic flux on the two close by, active- and low-mass K dwarfs: 61 Cygni A and Epsilon Eridani. Our results show that 61 Cygni A follows the Solar dynamo with a regular cycle at all wavelengths, while Epsilon Eridani represents a more extreme level of the Solar dynamo, while also showing strong Solar-like characteristics. For the first time we show magnetic butterfly diagrams for stars other than the Sun. For the two K stars and the Sun, the rate at which the toroidal field is generated from surface poloidal field is similar to the rate at which toroidal flux is lost through flux emergence. This suggests that the surface field plays a crucial role in the dynamos of all three stars. Finally, for Epsilon Eridani, we show that the two chromospheric cycle periods, of ~3 and ~13 years, correspond to two superimposed magnetic cycles.

preprint2020arXiv

Loss of toroidal magnetic flux by emergence of bipolar magnetic regions

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.

preprint2020arXiv

Power spectra of solar brightness variations at different inclinations

Magnetic features on the surfaces of cool stars cause variations of their brightness. Such variations have been extensively studied for the Sun. Recent planet-hunting space telescopes allowed measuring brightness variations in hundred thousands of other stars. The new data posed the question of how typical is the Sun as a variable star. Putting solar variability into the stellar context suffers, however, from the bias of solar observations being made from its near-equatorial plane, whereas stars are observed at all possible inclinations. We model solar brightness variations at timescales from days to years as they would be observed at different inclinations. In particular, we consider the effect of the inclination on the power spectrum of solar brightness variations. The variations are calculated in several passbands routinely used for stellar measurements. We employ the Surface Flux Transport Model (SFTM) to simulate the time-dependent spatial distribution of magnetic features on both near- and far-sides of the Sun. This distribution is then used to calculate solar brightness variations following the SATIRE (Spectral And Total Irradiance REconstruction) approach. We have quantified the effect of the inclination on solar brightness variability at timescales down to a day. Thus, our results allow making solar brightness records directly comparable to those obtained by the planet-hunting space telescopes. Furthermore, we decompose solar brightness variations into the components originating from the solar rotation and from the evolution of magnetic features.

preprint2020arXiv

Power spectrum of turbulent convection in the solar photosphere

The solar photosphere provides us with a laboratory for understanding turbulence in a layer where the fundamental processes of transport vary rapidly and a strongly superadiabatic region lies very closely to a subadiabatic layer. Our tools for probing the turbulence are high-resolution spectropolarimetric observations such as have recently been obtained with the two sunrise missions, and numerical simulations. Our aim is to study photospheric turbulence with the help of Fourier power spectra that we compute from observations and simulations. We also attempt to explain some properties of the photospheric overshooting flow with the help of its governing equations and simulations. We find that quiet-Sun observations and smeared simulations exhibit a power-law behavior in the subgranular range of their Doppler velocity power spectra with an index of$~\approx -2$. The unsmeared simulations exhibit a power-law index of$~\approx -2.25$. The smearing considerably reduces the extent of the power-law-like portion of the spectra. Therefore, the limited spatial resolution in some observations might eventually result in larger uncertainties in the estimation of the power-law indices. The simulated vertical velocity power spectra as a function of height show a rapid change in the power-law index from the solar surface to $300$~km above it. A scale-dependent transport of the vertical momentum occurs. At smaller scales, the vertical momentum is more efficiently transported sideways than at larger scales. This results in less vertical velocity power transported upward at small scales than at larger scales and produces a progressively steeper vertical velocity power law below $180$ km. Above this height, the gravity work progressively gains importance at all scales, making the atmosphere progressively more hydrostatic and resulting in a gradually less steep power law.

preprint2020arXiv

Reply to the comment of T. Metcalfe and J. van Saders on the Science report "The Sun is less active than other solar-like stars"

This is our reply to the comment of T. Metcalfe and J. van Saders on the Science report "The Sun is less active than other solar-like stars" by T. Reinhold, A. I. Shapiro, S. K. Solanki, B. T. Montet, N. A. Krivova, R. H. Cameron, E. M. Amazo-Gomez. We hope that both the comment and our reply lead to fruitful discussions which of the two presented scenarios is more likely.

preprint2020arXiv

Rossby modes in slowly rotating stars: depth dependence in distorted polytropes with uniform rotation

Large-scale Rossby waves have recently been discovered from measurements of horizontal surface and near-surface solar flows (Löptien at al. 2018). We are interested in understanding why only the sectoral modes are seen in the observations and also in modelling the radial structure of the observed modes. To do so, we characterise here the radial eigenfunctions of r modes for slowly-rotating polytropes in uniform rotation. We find that for free-surface boundary conditions on a spheroid of non-vanishing surface density, r modes can only exist for $\ell=m$ spherical harmonics in the inviscid case, and we compute their depth dependence and frequencies to leading order. For quasi-adiabatic stratification the sectoral modes with no radial nodes are the only modes which are almost toroidal and the depth dependence of the corresponding horizontal motion scales as $r^m$. For all r modes except the zero radial order sectoral ones, non-adiabatic stratification plays a crucial role in the radial force balance. The lack of quasi-toroidal solutions when stratification is close to neutral, except for the sectoral modes without nodes in radius, follows from the statement that the system needs to be in both horizontal and radial force balance. In the absence of super- or subadiabatic stratification and viscosity, both the horizontal and radial force balances independently determine the pressure perturbation. The only quasi-toroidal cases in which the two determinations of the pressure perturbation are consistent are the special cases where $\ell=m$, and the horizontal displacement scales with $r^m$.

preprint2020arXiv

Simulations Show that Vortex Flows could Heat the Chromosphere in Solar Plage

The relationship between vortex flows at different spatial scales and their contribution to the energy balance in the chromosphere is not yet fully understood. We perform three-dimensional (3D) radiation-magnetohydrodynamic (MHD) simulations of a unipolar solar plage region at a spatial resolution of 10 km using the MURaM code. We use the swirling-strength criterion that mainly detects the smallest vortices present in the simulation data. We additionally degrade our simulation data to smooth-out the smaller vortices, so that also the vortices at larger spatial scales can be detected. Vortex flows at various spatial scales are found in our simulation data for different effective spatial resolutions. We conclude that the observed large vortices are likely clusters of much smaller ones that are not yet resolved by observations. We show that the vertical Poynting flux decreases rapidly with reduced effective spatial resolutions and is predominantly carried by the horizontal plasma motions rather than vertical flows. Since the small-scale horizontal motions or the smaller vortices carry most of the energy, the energy transported by vortices deduced from low resolution data is grossly underestimated. In full resolution simulation data, the Poynting flux contribution due to vortices is more than adequate to compensate for the radiative losses in plage, indicating their importance for chromospheric heating.

preprint2016arXiv

MESA meets MURaM: Surface effects in main-sequence solar-like oscillators computed using three-dimensional radiation hydrodynamics simulations

... [C]urrent stellar models predict oscillation frequencies that are systematically affected by simplified modelling of the near-surface layers. We use three-dimensional radiation hydrodynamics simulations to better model the near-surface equilibrium structure of dwarfs with spectral types F3, G2, K0 and K5, and examine the differences between oscillation mode frequencies. ... We precisely match stellar models to the simulations' gravities and effective temperatures at the surface, and to the temporally- and horizontally-averaged densities and pressures at their deepest points. We then replace the near-surface structure with that of the averaged simulation and compute the change in the oscillation mode frequencies. We also fit the differences using several parametric models currently available in the literature. The surface effect in the stars of solar-type and later is qualitatively similar and changes steadily with decreasing effective temperature. In particular, the point of greatest frequency difference decreases slightly as a fraction of the acoustic cut-off frequency and the overall scale of the surface effect decreases. The surface effect in the hot, F3-type star follows the same trend in scale (i.e. it is larger in magnitude) but shows a different overall variation with mode frequency. We find that the two-term fit by Ball & Gizon (2014) is best able to reproduce the surface terms across all four spectral types, although the scaled solar term and a modified Lorentzian function also match the three cooler simulations reasonably well. ... Our simplified results suggest that the surface effect is generally larger in hotter stars (and correspondingly smaller in cooler stars) and of similar shape in stars of solar type and cooler. However, we cannot presently predict whether this will remain so when other components of the surface effect are included.

preprint2016arXiv

The turbulent diffusion of toroidal magnetic flux as inferred from properties of the sunspot butterfly diagram

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.

preprint2015arXiv

Limitations of force-free magnetic field extrapolations: revisiting basic assumptions

Force-free extrapolations are widely used to study the magnetic field in the solar corona based on surface measurements. The extrapolations assume that the ratio of internal energy of the plasma to magnetic energy, the plasma-beta is negligible. Despite the widespread use of this assumption observations, models, and theoretical considerations show that beta is of the order of a few percent to more than 10%, and thus not small. We investigate what consequences this has for the reliability of extrapolation results. We use basic concepts starting with the force and the energy balance to infer relations between plasma-beta and free magnetic energy, to study the direction of currents in the corona with respect to the magnetic field, and to estimate the errors in the free magnetic energy by neglecting effects of the plasma (beta<<1). A comparison with a 3D MHD model supports our basic considerations. If plasma-beta is of the order of the relative free energy (the ratio of the free magnetic energy to the total magnetic energy) then the pressure gradient can balance the Lorentz force. This is the case in the solar corona, and therefore the currents are not properly described. In particular the error in terms of magnetic energy by neglecting the plasma is of the order of the free magnetic energy, so that the latter can not be reliably determined by an extrapolation. While a force-free extrapolation might capture the magnetic structure and connectivity of the coronal magnetic field, the derived currents and free magnetic energy are not reliable. Thus quantitative results of extrapolations on the location and amount of heating in the corona (through current dissipation) and on the energy storage of the magnetic field (e.g. for eruptive events) are limited.

preprint2014arXiv

Effects of the scatter in sunspot group tilt angles on the large-scale magnetic field at the solar surface

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.

preprint2014arXiv

Magnetic Flux Transport at the Solar Surface

After emerging to the solar surface, the Sun's magnetic field displays a complex and intricate evolution. The evolution of the surface field is important for several reasons. One is that the surface field, and its dynamics, sets the boundary condition for the coronal and heliospheric magnetic fields. Another is that the surface evolution gives us insight into the dynamo process. In particular, it plays an essential role in the Babcock-Leighton model of the solar dynamo. Describing this evolution is the aim of the surface flux transport model. The model starts from the emergence of magnetic bipoles. Thereafter, the model is based on the induction equation and the fact that after emergence the magnetic field is observed to evolve as if it were purely radial. The induction equation then describes how the surface flows -- differential rotation, meridional circulation, granular, supergranular flows, and active region inflows -- determine the evolution of the field (now taken to be purely radial). In this paper, we review the modeling of the various processes that determine the evolution of the surface field. We restrict our attention to their role in the surface flux transport model. We also discuss the success of the model and some of the results that have been obtained using this model.

preprint2014arXiv

Migration of Ca II H bright points in the internetwork

The migration of magnetic bright point-like features (MBP) in the lower solar atmosphere reflects the dispersal of magnetic flux as well as the horizontal flows of the atmospheric layer they are embedded in. We analyse trajectories of the proper motion of intrinsically magnetic, isolated internetwork Ca II H MBPs (mean lifetime 461 +- 9 s) to obtain their diffusivity behaviour. We use seeing-free high spatial and temporal resolution image sequences of quiet-Sun, disc-centre observations obtained in the Ca II H 3968 Å passband of the Sunrise Filter Imager (SuFI) onboard the Sunrise balloon-borne solar observatory. Small MBPs in the internetwork are automatically tracked. The trajectory of each MBP is then calculated and described by a diffusion index (γ) and a diffusion coefficient (D). We further explore the distribution of the diffusion indices with the help of a Monte Carlo simulation. We find γ = 1.69 +- 0.08 and D = 257 +- 32 km^2/s averaged over all MBPs. Trajectories of most MBPs are classified as super-diffusive, i.e., γ > 1, with the determined γ being to our knowledge the largest obtained so far. A direct correlation between D and time-scale (τ) determined from trajectories of all MBPs is also obtained. We discuss a simple scenario to explain the diffusivity of the observed, relatively short-lived MBPs while they migrate within a small area in a supergranule (i.e., an internetwork area). We show that the scatter in the γ values obtained for individual MBPs is due to their limited lifetimes. The super-diffusive MBPs can be well-described as random walkers (due to granular evolution and intergranular turbu- lence) superposed on a large systematic (background) velocity, caused by granular, mesogranular and supergranular flows.

preprint2014arXiv

Simulated Magnetic Flows in the Solar Photosphere

Recent Sunrise/IMaX observations have revealed the existence of supersonic magnetic flows. Our aim is to determine the origin of such flows by using realistic MHD simulations. We simulate cancellation and emergence of magnetic flux through the solar photosphere. Our first numerical experiment starts with magnetic field of both polarities. To simulate emergence into a region with pre-existing field, we introduce a large-scale horizontally uniform sheet of horizontal field. We follow the subsequent evolution, creating synthetic polarimetric observations, including known instrumental effects of the Sunrise/IMaX and Hinode/SP instruments. We compare the simulated and observed spectropolarimetric signals. Strongly blue- and redshifted Stokes V signals are produced in locations where strong line-of-sight velocities coincide with the strong line-of-sight component of magnetic field. The size and strength of simulated events is smaller than observed and they are mostly associated with downflows, contrary to observations. In a few cases where they appear above a granule, single blue lobed Stokes V are produced due to strong gradients in magnetic field and velocity. No change of magnetic field sign is detected along the line of sight in these instances. More high-speed magnetized flows occur in the case where emergence is simulated than in the case where no horizontal field was added. The simulations indicate that the observed events result from magnetic flux emergence, where reconnection may take place but does not seem to be necessary.

preprint2014arXiv

The Role of Subsurface Flows in Solar Surface Convection: Modeling the Spectrum of Supergranular and Larger Scale Flows

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.

preprint2013arXiv

Limits to solar cycle predictability: Cross-equatorial flux plumes

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.

preprint2013arXiv

Modeling Solar Cycles 15 to 21 Using a Flux Transport Dynamo

Context: The Sun's polar fields and open flux around the time of activity minima have been considered to be strongly correlated with the strength of the subsequent maximum of solar activity. Aims: We aim to investigate the behavior of a Babcock-Leighton dynamo with a source poloidal term that is based on the observed sunspot areas and tilts. In particular, we investigate whether the toroidal fields at the base of convection zone from the model are correlated with the observed solar cycle activity maxima. Methods: We used a flux transport dynamo model that includes convective pumping and a poloidal source term based on the historical record of sunspot group areas, locations, and tilt angles to simulate solar cycles 15 to 21. Results: We find that the polar fields near minima and the toroidal flux at the base of the convection zone are both highly correlated with the subsequent maxima of solar activity levels (r = 0.85 and r = 0.93, respectively). Conclusions: The Babcock-Leighton dynamo is consistent with the observationally inferred correlations.

preprint2012arXiv

Are the strengths of solar cycles determined by converging flows towards the activity belts?

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.

preprint2012arXiv

Break up of returning plasma after the 7 June 2011 filament eruption by Rayleigh-Taylor instabilities

A prominence eruption on 7 June 2011 produced spectacular curtains of plasma falling through the lower corona. At the solar surface they created an incredible display of extreme ultraviolet brightenings. The aim is to identify and analyze some of the local instabilities which produce structure in the falling plasma. The structures were investigated using SDO/AIA 171A and 193A images in which the falling plasma appeared dark against the bright coronal emission. Several instances of the Rayleigh-Taylor instability were investigated. In two cases the Alfven velocity associated with the dense plasma could be estimated from the separation of the Rayleigh-Taylor fingers. A second type of feature, which has the appearance of self-similar branching horns, is also discussed.

preprint2012arXiv

Simulations of the solar near-surface layers with the CO5BOLD, MURaM, and Stagger codes

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.

preprint2012arXiv

Surface flux evolution constraints for flux transport dynamos

The surface flux transport (SFT) model of solar magnetic fields involves empirically well-constrained velocity and magnetic fields. The basic evolution of the Sun's large-scale surface magnetic field is well described by this model. The azimuthally averaged evolution of the SFT model can be compared to the surface evolution of the flux transport dynamo (FTD), and the evolution of the SFT model can be used to constrain several near-surface properties of the FTD model. We compared the results of the FTD model with different upper boundary conditions and diffusivity profiles against the results of the SFT model. Among the ingredients of the FTD model, downward pumping of magnetic flux, related to a positive diffusivity gradient, has a significant effect in slowing down the diffusive radial transport of magnetic flux through the solar surface. Provided the pumping was strong enough to give rise to a downflow of a magnetic Reynolds number of 5 in the near-surface boundary layer, the FTD using a vertical boundary condition matches the SFT model based on the average velocities above the boundary layer. The FTD model with a potential field were unable to match the SFT results.

preprint2012arXiv

Vortices, shocks, and heating in the solar photosphere: effect of a magnetic field

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.

preprint2012arXiv

Waves as the source of apparent twisting motions in sunspot penumbrae

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.

preprint2011arXiv

EUV jets, type III radio bursts and sunspot waves investigated using SDO/AIA observations

Images from the Solar Dynamics Observatory (SDO) at 211A are used to identify the solar source of the type III radio bursts seen in WIND/WAVES dynamic spectra. We analyse a 2.5 hour period during which six strong bursts are seen. The radio bursts correlate very well with the EUV jets coming from the western side of a sunspot in AR11092. The EUV jet emission also correlates well with brightening at what looks like their footpoint at the edge of the umbra. For 10-15 min after strong EUV jets are ejected, the footpoint brightens at roughly 3 min intervals. In both the EUV images and the extracted light curves, it looks as though the brightening is related to the 3-min sunspot oscillations, although the correlation coefficient is rather low. The only open field near the jets is rooted in the sunspot. We conclude that active region EUV/X-ray jets and interplanetary electron streams originate on the edge of the sunspot umbra. They form along a current sheet between the sunspot open field and closed field connecting to underlying satellite flux. Sunspot running penumbral waves cause roughly 3-min jet footpoint brightening. The relationship between the waves and jets is less clear.

preprint2011arXiv

Quenching of the alpha effect in the Sun -- what observations are telling us

The Babcock-Leighton type of dynamo has received recent support in terms of the discovery in the observational records of systematic cycle-to-cycle variations in the tilt angle of sunspot groups. It has been proposed that these variations might be the consequence of the observed inflow into the active region belt. Furthermore simulations have shown that such inflows restrict the creation of net poloidal flux, in effect acting to quench the alpha effect associated with the Coriolis force acting on rising flux tubes. In this paper we expand on these ideas.

preprint2011arXiv

Transport of magnetic flux from the canopy to the internetwork

Recent observations have revealed that 8% of linear polarization patches in the internetwork quiet Sun are fully embedded in downflows. These are not easily explained with the typical scenarios for the source of internetwork fields which rely on flux emergence from below. We explore using radiative MHD simulations a scenario where magnetic flux is transported from the magnetic canopy overlying the internetwork into the photosphere by means of downward plumes associated with convective overshoot. We find that if a canopy-like magnetic field is present in the simulation, the transport of flux from the canopy is an important process for seeding the photospheric layers of the internetwork with magnetic field. We propose that this mechanism is relevant for the Sun as well, and it could naturally explain the observed internetwork linear polarization patches entirely embedded in downflows.

preprint2011arXiv

Universality of the Small-Scale Dynamo Mechanism

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.

preprint2011arXiv

Vortices in simulations of solar surface convection

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.

preprint2010arXiv

Surface flux transport modeling for solar cycles 15--21: effects of cycle-dependent tilt angles of sunspot groups

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.

preprint2010arXiv

The effect of activity-related meridional flow modulation on the strength of the solar polar magnetic field

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.