Showing 1–50 of 78 results for author: Cameron, R H

The mean solar butterfly diagram and poloidal field generation rate at the surface of the Sun

Authors: Simon Cloutier, Robert H. Cameron, Laurent Gizon

Abstract: The difference between individual solar cycles in the magnetic butterfly diagram can mostly be ascribed to the stochasticity of the emergence process. We aim to obtain the expectation value of the butterfly diagram from observations of four cycles. This allows us to further determine the generation rate of the surface radial magnetic field. We use data from Wilcox Solar Observatory to generate tim… ▽ More The difference between individual solar cycles in the magnetic butterfly diagram can mostly be ascribed to the stochasticity of the emergence process. We aim to obtain the expectation value of the butterfly diagram from observations of four cycles. This allows us to further determine the generation rate of the surface radial magnetic field. We use data from Wilcox Solar Observatory to generate time-latitude diagrams spanning cycles 21 to 24 of the surface radial and toroidal magnetic fields, symmetrize them across the equator and cycle-average them. From the mean butterfly diagram and surface toroidal field we then infer the mean poloidal field generation rate at the surface of the Sun. The averaging procedure removes realization noise from individual cycles. The amount of emerging flux required to account for the evolution of the surface radial field is found to match that provided by the observed surface toroidal field and Joy's law. Cycle-averaging butterfly diagrams removes realization noise and artefacts due to imperfect scale separation, and corresponds to an ensemble average that can be interpreted in the mean-field framework. The result can then be directly compared to $αΩ$-type dynamo models. The Babcock-Leighton $α$-effect is consistent with observations, a result that can be appreciated only if the observational data is averaged in some way. △ Less

Submitted 27 May, 2024; originally announced May 2024.

Comments: 5 pages, 4 figures, submitted to A&A

The Sun's differential rotation is controlled by high-latitude baroclinically unstable inertial modes

Authors: Yuto Bekki, Robert H. Cameron, Laurent Gizon

Abstract: Rapidly rotating fluids have a rotation profile which depends only on the distance from the rotation axis, in accordance with the Taylor-Proudman theorem. Although the Sun was expected to be such a body, helioseismology showed that the rotation rate in the convection zone is closer to constant on radii. It has been postulated that this deviation is due to the poles being warmer than the equator by… ▽ More Rapidly rotating fluids have a rotation profile which depends only on the distance from the rotation axis, in accordance with the Taylor-Proudman theorem. Although the Sun was expected to be such a body, helioseismology showed that the rotation rate in the convection zone is closer to constant on radii. It has been postulated that this deviation is due to the poles being warmer than the equator by a few degrees. Using numerical simulations, we show that the pole-to-equator temperature difference cannot exceed 7 Kelvin as a result of the back-reaction of the high-latitude baroclinically unstable inertial modes. The observed amplitudes of the modes further indicate that this maximum temperature difference is reached in the Sun. We conclude that the Sun's latitudinal differential rotation reaches its maximum allowed value. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Comments: 52 pages, 21 figures, 2 tables. Published in Science Advances on 27 Mar 2024. The authors' version

Journal ref: Science Advances, 10 (13), eadk5643 (2024)

Stellar Surface Magnetic Fields Impact Limb Darkening

Authors: Nadiia M. Kostogryz, Alexander I. Shapiro, Veronika Witzke, Robert H. Cameron, Laurent Gizon, Natalie A. Krivova, Hans-G. Ludwig, Pierre F. L. Maxted, Sara Seager, Sami K. Solanki, Jeff Valenti

Abstract: Stars appear darker at their limbs than at their disk centers because at the limb we are viewing the higher and cooler layers of stellar photospheres. Limb darkening derived from state-of-the-art stellar atmosphere models systematically fails to reproduce recent transiting exoplanet light curves from the Kepler, TESS, and JWST telescopes -- stellar brightness obtained from measurements drops less… ▽ More Stars appear darker at their limbs than at their disk centers because at the limb we are viewing the higher and cooler layers of stellar photospheres. Limb darkening derived from state-of-the-art stellar atmosphere models systematically fails to reproduce recent transiting exoplanet light curves from the Kepler, TESS, and JWST telescopes -- stellar brightness obtained from measurements drops less steeply towards the limb than predicted by models. All previous models assumed atmosphere devoid of magnetic fields. Here we use our new stellar atmosphere models computed with the 3D radiative magneto-hydrodynamic code MURaM to show that small-scale concentration of magnetic fields on the stellar surface affect limb darkening at a level that allows us to explain the observations. Our findings provide a way forward to improve the determination of exoplanet radii and especially the transmission spectroscopy analysis for transiting planets, which relies on a very accurate description of stellar limb darkening from the visible through the infrared. Furthermore, our findings imply that limb darkening allows measuring the small-scale magnetic field on stars with transiting planets. △ Less

Submitted 29 February, 2024; originally announced March 2024.

Comments: 19 pages, 6 figures This is the submission version of manuscript to Nature Astronomy

Small-scale Dynamo in Cool Stars III. Changes in the photospheres of F3V to M0V stars

Authors: Tanayveer S. Bhatia, Robert H. Cameron, Hardi Peter, Sami K. Solanki

Abstract: Some of the quiet solar magnetic flux could be attributed to a small-scale dynamo (SSD) operating in the convection zone. An SSD operating in cool main-sequence stars is expected to affect the atmospheric structure, in particular the convection, and should have observational signatures. We aim to investigate the distribution of these fields as well as their effect on intensity characteristics, vel… ▽ More Some of the quiet solar magnetic flux could be attributed to a small-scale dynamo (SSD) operating in the convection zone. An SSD operating in cool main-sequence stars is expected to affect the atmospheric structure, in particular the convection, and should have observational signatures. We aim to investigate the distribution of these fields as well as their effect on intensity characteristics, velocities and spatial distribution of kinetic (KE) and magnetic energy (ME) in the lower photosphere of spectral types F3V, G2V, K0V and M0V using 3D radiative-MHD simulations. PDFs of field strength at the $τ=1$ surface are quite similar for all cases. The M0V star displays the strongest fields, but relative to the gas pressure, the fields on the F3V star reach the largest values. All stars display an excess of horizontal field relative to vertical field in the middle photosphere, with this excess becoming increasingly prominent towards later spectral types. These fields result in a decrease in upflow velocities, slightly smaller granules as well as the formation of bright points in intergranular lanes. The spatial distribution of KE and ME is also similar for all cases, implying a simple pressure scale height proportionality of important scales. SSD fields have rather similar effects on the photospheres of cool main-sequence stars, namely, a significant reduction in convective velocities as well as a slight reduction in granule size, and concentration of field to kG levels in intergranular lanes associated with the formation of bright points. The distribution of field strengths and energies is also rather similar. △ Less

Submitted 3 November, 2023; originally announced November 2023.

Comments: 11 pages, 14 figures, accepted for publication in A&A

A Babcock-Leighton dynamo model of the Sun incorporating toroidal flux loss and the helioseismically-inferred meridional flow

Authors: S. Cloutier, R. H. Cameron, L. Gizon

Abstract: We investigate whether the Babcock-Leighton flux-transport dynamo model remains in agreement with observations if the meridional flow profile is taken from helioseismic inversions. Additionally, we investigate the effect of the loss of toroidal flux through the solar surface. We employ the 2D flux-transport BL dynamo framework. We use the helioseismically-inferred meridional flow profile, and incl… ▽ More We investigate whether the Babcock-Leighton flux-transport dynamo model remains in agreement with observations if the meridional flow profile is taken from helioseismic inversions. Additionally, we investigate the effect of the loss of toroidal flux through the solar surface. We employ the 2D flux-transport BL dynamo framework. We use the helioseismically-inferred meridional flow profile, and include toroidal flux loss in a way that is consistent with the amount of poloidal flux generated by Joy's law. Our model does not impose a preference for emergences at low latitudes, but we do require that the model produces such a preference. We can find solutions in general agreement with observations, including the equatorward drift of the butterfly wings and the cycle's 11 year period. The most important free parameters in the model are the depth to which the radial turbulent pumping extends and the turbulent diffusivity in the lower half of the convection zone. We find that the pumping needs to extend to depths of about $0.80R_{\odot}$ and the bulk turbulent diffusivity needs to be around 10 km$^2$/s or less. We find that the emergences are restricted to low latitudes without the need to impose such a preference. The flux-transport BL model, incorporating the helioseismically inferred meridional flow and toroidal field loss term, is compatible with the properties of the observed butterfly diagram and with the observed toroidal loss rate. Reasonably tight constraints are placed on the remaining free parameters. The pumping needs to be to just below the depth corresponding to the location where the meridional flow changes direction. Our linear model does not however reproduce the observed "rush to the poles" of the diffuse surface radial field resulting from the decay of sunspots -- reproducing this might require the imposition of a preference for flux to emerge near the equator. △ Less

Submitted 29 September, 2023; originally announced September 2023.

Comments: 10 pages, 8 figures, accepted by A&A

Journal ref: A&A 680, A42 (2023)

No evidence for synchronization of the solar cycle by a "clock"

Authors: E. Weisshaar, R. H. Cameron, M. Schüssler

Abstract: 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 obser… ▽ More 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. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: Astronomy & Astrophysics, in press

Journal ref: A&A 671, A87 (2023)

Exploring the Solar Poles: The Last Great Frontier of the Sun

Authors: Dibyendu Nandy, Dipankar Banerjee, Prantika Bhowmik, Allan Sacha Brun, Robert H. Cameron, S. E. Gibson, Shravan Hanasoge, Louise Harra, Donald M. Hassler, Rekha Jain, Jie Jiang, Laurène Jouve, Duncan H. Mackay, Sushant S. Mahajan, Cristina H. Mandrini, Mathew Owens, Shaonwita Pal, Rui F. Pinto, Chitradeep Saha, Xudong Sun, Durgesh Tripathi, Ilya G. Usoskin

Abstract: Despite investments in multiple space and ground-based solar observatories by the global community, the Sun's polar regions remain unchartered territory - the last great frontier for solar observations. Breaching this frontier is fundamental to understanding the solar cycle - the ultimate driver of short-to-long term solar activity that encompasses space weather and space climate. Magnetohydrodyna… ▽ More Despite investments in multiple space and ground-based solar observatories by the global community, the Sun's polar regions remain unchartered territory - the last great frontier for solar observations. Breaching this frontier is fundamental to understanding the solar cycle - the ultimate driver of short-to-long term solar activity that encompasses space weather and space climate. Magnetohydrodynamic dynamo models and empirically observed relationships have established that the polar field is the primary determinant of the future solar cycle amplitude. Models of solar surface evolution of tilted active regions indicate that the mid to high latitude surges of magnetic flux govern dynamics leading to the reversal and build-up of polar fields. Our theoretical understanding and numerical models of this high latitude magnetic field dynamics and plasma flows - that are a critical component of the sunspot cycle - lack precise observational constraints. This limitation compromises our ability to observe the enigmatic kilo Gauss polar flux patches and constrain the polar field distribution at high latitudes. The lack of these observations handicap our understanding of how high latitude magnetic fields power polar jets, plumes, and the fast solar wind that extend to the boundaries of the heliosphere and modulate solar open flux and cosmic ray flux within the solar system. Accurate observation of the Sun's polar regions, therefore, is the single most outstanding challenge that confronts Heliophysics. This paper argues the scientific case for novel out of ecliptic observations of the Sun's polar regions, in conjunction with existing, or future multi-vantage point heliospheric observatories. Such a mission concept can revolutionize the field of Heliophysics like no other mission concept has - with relevance that transcends spatial regimes from the solar interior to the heliosphere. △ Less

Submitted 30 December, 2022; originally announced January 2023.

Comments: This White Paper was submitted in 2022 to the United States National Academies Solar and Space Physics (Heliophysics) Decadal Survey

Three-dimensional non-kinematic simulation of post-emergence evolution of bipolar magnetic regions and Babcock-Leighton dynamo of the Sun

Authors: Yuto Bekki, Robert H. Cameron

Abstract: The Babcock-Leighton (BL) flux-transport model is a widely-accepted dynamo model of the Sun. This dynamo model has been extensively studied in a two-dimensional (2D) mean-field framework in both kinematic and non-kinematic regimes. Recent three-dimensional (3D) models have been restricted to the kinematic regime. In these models, the surface poloidal flux is produced by the emergence of bipolar ma… ▽ More The Babcock-Leighton (BL) flux-transport model is a widely-accepted dynamo model of the Sun. This dynamo model has been extensively studied in a two-dimensional (2D) mean-field framework in both kinematic and non-kinematic regimes. Recent three-dimensional (3D) models have been restricted to the kinematic regime. In these models, the surface poloidal flux is produced by the emergence of bipolar magnetic regions (BMRs) that are tilted according to Joy's law. We investigate the prescription for emergence of a BMR in 3D non-kinematic simulations. We also report initial results of cyclic BL dynamo simulation. We extend a conventional 2D mean-field model of the BL flux-transport dynamo into 3D non-kinematic regime. The large-scale mean flows are driven by the parameterized $Λ$-effect in this model. For the induction equation, we use a BL source term by which the surface BMRs are produced in response to the dynamo-generated toroidal field inside the convection zone. We find that, in the 3D non-kinematic regime, the tilt angle of a newly-emerged BMR is very sensitive to the prescription for the subsurface structure of the BMR. Anti-Joy tilt angles are found unless the BMR is deeply embedded in the convection zone. We also find that the leading spot tends to become stronger than the following spot. The anti-Joy's law trend and the morphological asymmetry of the BMRs can be explained by the Coriolis force acting on the Lorentz-force-driven flows. Furthermore, we demonstrate that the solar-like magnetic cycles can be successfully obtained if the Joy's law is explicitly given in the BL $α$-effect. In these cyclic dynamo simulation, a strong Lorentz force feedback leads to cycle modulations in the differential rotation and meridional circulation. The non-axisymmetric components of the flows are found to exist as inertial modes such as the equatorial Rossby modes. △ Less

Submitted 7 December, 2022; v1 submitted 16 September, 2022; originally announced September 2022.

Comments: 18 pages, 19 figures, accepted for publication in A&A

Journal ref: A&A 670, A101 (2023)

Theory of solar oscillations in the inertial frequency range: Amplitudes of equatorial modes from a nonlinear rotating convection simulation

Authors: Yuto Bekki, Robert H. Cameron, Laurent Gizon

Abstract: Several types of inertial modes have been detected on the Sun. Properties of these inertial modes have been studied in the linear regime but have not been studied in nonlinear simulations of solar rotating convection. Comparing the nonlinear simulations, the linear theory, and the solar observations is important to better understand the differences between the models and the real Sun. We wish to d… ▽ More Several types of inertial modes have been detected on the Sun. Properties of these inertial modes have been studied in the linear regime but have not been studied in nonlinear simulations of solar rotating convection. Comparing the nonlinear simulations, the linear theory, and the solar observations is important to better understand the differences between the models and the real Sun. We wish to detect and characterize the modes present in a nonlinear numerical simulation of solar convection, in particular to understand the amplitudes and lifetimes of the modes. We developed a code with a Yin-Yang grid to carry out fully-nonlinear numerical simulations of rotating convection in a spherical shell. The stratification is solar-like up to 0.96R. The simulations cover a duration of about 15 solar years. Various large-scale modes at low frequencies are extracted from the simulation. Their characteristics are compared to those from the linear model and to the observations. Among other modes, both the equatorial Rossby modes and the columnar convective modes are seen in the simulation. The columnar convective modes contain most of the large-scale velocity power outside the tangential cylinder and substantially contribute to the heat and angular momentum transport. Equatorial Rossby modes with no radial node (n=0) are also found: They have the same spatial structures as the linear eigenfunctions. They are stochastically excited by convection and have the amplitudes of a few m/s and mode linewidths of about 20-30 nHz, which are comparable to those observed on the Sun. We also confirm the existence of the mixed modes between the equatorial Rossby modes and the columnar convective modes in our nonlinear simulation, as predicted by the linear eigenmode analysis. We also see the high-latitude mode with m=1 in our nonlinear simulation but its amplitude is much weaker than that observed on the Sun. △ Less

Submitted 23 August, 2022; originally announced August 2022.

Comments: 18 pages, 21 figures, 1 table. Astronomy & Astrophysics

Journal ref: A&A 666, A135 (2022)

Faculae cancel out on the surfaces of active Suns

Authors: N. -E. Nèmec, A. I. Shapiro, E. Işık, K. Sowmya, S. K. Solanki, N. A. Krivova, R. H. Cameron, L. Gizon

Abstract: 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… ▽ More 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. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: 10 pages, 6 Figures, accepted for publication in ApJL

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

Authors: C. Baumgartner, A. C. Birch, H. Schunker, R. H. Cameron, L. Gizon

Abstract: 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 proxi… ▽ More 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$. △ Less

Submitted 5 July, 2022; originally announced July 2022.

Comments: accepted in A&A

Journal ref: A&A 664, A183 (2022)

Small-scale Dynamo in Cool Stars: I. Changes in stratification and near-surface convection for main-sequence spectral types

Authors: Tanayveer S. Bhatia, Robert H. Cameron, Sami K. Solanki, Hardi Peter, Damien Przybylski, Veronika Witzke, Alexander Shapiro

Abstract: Some of the small-scale solar magnetic flux can be attributed to a small-scale dynamo (SSD) operating in the near-surface convection. The SSD fields have consequences for solar granular convection, basal flux, as well as chromospheric heating. A similar SSD mechanism is expected to be active in the near-surface convection of other cool main-sequence stars, but this has never been investigated. We… ▽ More Some of the small-scale solar magnetic flux can be attributed to a small-scale dynamo (SSD) operating in the near-surface convection. The SSD fields have consequences for solar granular convection, basal flux, as well as chromospheric heating. A similar SSD mechanism is expected to be active in the near-surface convection of other cool main-sequence stars, but this has never been investigated. We aim to investigate changes in stratification and convection due to inclusion of SSD fields for F3V, G2V, K0V and M0V spectral types in the near-surface convection. 3D magnetohydrodynamic (MHD) models of the four stellar boxes, covering the subsurface convection zone up to the lower photosphere in a small cartesian box, are studied using the \textit{MURaM} radiative-MHD simulation code. The SSD runs are compared against reference hydrodynamic runs. An SSD is found to efficiently produce magnetic field with energies ranging between 5\% to 80\% of the plasma kinetic energy at different depths. This ratio tends to be larger for larger $T_{\mathrm{eff}}$. The relative change in density and gas pressure stratification for the deeper convective layers due to SSD magnetic fields is negligible, except for the F-star. For the F-star, there is a substantial reduction in convective velocities due to Lorentz force feedback from magnetic fields, which, in turn, reduces the turbulent pressure. SSD in near-surface convection for cool main-sequence stars introduces small but significant changes in thermodynamic stratification (especially for the F-star) due to reduction in convective velocities. △ Less

Submitted 20 July, 2023; v1 submitted 31 May, 2022; originally announced June 2022.

Comments: Update (20.07.2023): an errata (appendix D) has been included to correct a small error in the derivation of turbulent pressure in appendix B. In addition, the text has been updated with the published version. 10 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 663, A166 (2022)

Theory of solar oscillations in the inertial frequency range: Linear modes of the convection zone

Authors: Yuto Bekki, Robert H. Cameron, Laurent Gizon

Abstract: On the one hand, several types of global-scale inertial modes of oscillation have been observed on the Sun. They include the equatorial Rossby modes, critical-latitude modes, and high-latitude modes. On the other hand, the columnar convective modes (predicted by simulations; also known as banana cells or thermal Rossby waves) remain elusive. We aim to investigate the influence of turbulent diffusi… ▽ More On the one hand, several types of global-scale inertial modes of oscillation have been observed on the Sun. They include the equatorial Rossby modes, critical-latitude modes, and high-latitude modes. On the other hand, the columnar convective modes (predicted by simulations; also known as banana cells or thermal Rossby waves) remain elusive. We aim to investigate the influence of turbulent diffusivities, non-adiabatic stratification, differential rotation, and a latitudinal entropy gradient on the linear global modes of the rotating solar convection zone. We solve numerically for the eigenmodes of a rotating compressible fluid inside a spherical shell. We identify modes in the inertial frequency range including the columnar convective modes, as well as modes of mixed character. The corresponding mode dispersion relations and eigenfunctions are computed for azimuthal orders $m \leq 16$. The three main results are as follows. Firstly, we find that, for $m \gtrsim 5$, the radial dependence of the equatorial Rossby modes with no radial node ($n=0$) is radically changed from the traditional expectation ($r^m$) for turbulent diffusivities $\gtrsim 10^{12}$ cm$^2$ s$^{-1}$. Secondly, we find mixed modes, i.e. modes that share properties of the equatorial Rossby modes with one radial node ($n=1$) and the columnar convective modes. Thirdly, we show that the $m=1$ high-latitude mode in the model is consistent with the solar observations when the latitudinal entropy gradient corresponding to a thermal wind balance is included (baroclinally unstable mode). To our knowledge, this work is the first realistic eigenvalue calculation of the global modes of the rotating solar convection zone. This calculation reveals a rich spectrum of modes in the inertial frequency range, which can be directly compared to the observations. In turn, the observed modes can inform us about the solar convection zone. △ Less

Submitted 8 March, 2022; originally announced March 2022.

Comments: 24 pages, 29 figures, submitted to Astronomy & Astrophysics on 20 Jan 2022

Journal ref: A&A 662, A16 (2022)

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

Authors: S. V. Jeffers, R. H. Cameron, S. C. Marsden, S. Boro Saikia, C. P. Folsom, M. M. Jardine, J. Morin, P. Petit, V. See, A. A. Vidotto, U. Wolter, M. Mittag

Abstract: 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 diagnosti… ▽ More 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. △ Less

Submitted 19 January, 2022; originally announced January 2022.

Comments: 8 pages, 5 figures: Accepted by A&A

Solar inertial modes: Observations, identification, and diagnostic promise

Authors: Laurent Gizon, Robert H. Cameron, Yuto Bekki, Aaron C. Birch, Richard S. Bogart, Allan Sacha Brun, Cilia Damiani, Damien Fournier, Laura Hyest, Kiran Jain, B. Lekshmi, Zhi-Chao Liang, Bastian Proxauf

Abstract: The oscillations of a slowly rotating star have long been classified into spheroidal and toroidal modes. The spheroidal modes include the well-known 5-min acoustic modes used in helioseismology. Here we report observations of the Sun's toroidal modes, for which the restoring force is the Coriolis force and whose periods are on the order of the solar rotation period. By comparing the observations w… ▽ More The oscillations of a slowly rotating star have long been classified into spheroidal and toroidal modes. The spheroidal modes include the well-known 5-min acoustic modes used in helioseismology. Here we report observations of the Sun's toroidal modes, for which the restoring force is the Coriolis force and whose periods are on the order of the solar rotation period. By comparing the observations with the normal modes of a differentially rotating spherical shell, we are able to identify many of the observed modes. These are the high-latitude inertial modes, the critical-latitude inertial modes, and the equatorial Rossby modes. In the model, the high-latitude and critical-latitude modes have maximum kinetic energy density at the base of the convection zone, and the high-latitude modes are baroclinically unstable due to the latitudinal entropy gradient. As a first application of inertial-mode helioseismology, we constrain the superadiabaticity and the turbulent viscosity in the deep convection zone. △ Less

Submitted 20 July, 2021; originally announced July 2021.

Comments: Forthcoming article, Astronomy and Astrophysics

Slow magneto-acoustic waves in simulations of a solar plage region carry enough energy to heat the chromosphere

Authors: Nitin Yadav, Robert H. Cameron, Sami K. Solanki

Abstract: We study the properties of slow magneto-acoustic waves that are naturally excited due to turbulent convection and investigate their role in the energy balance of a plage region using three dimensional (3D) radiation-MHD simulations. We calculate the horizontally averaged (over the whole domain) frequency power spectra for both longitudinal and vertical (i.e. the component perpendicular to the surf… ▽ More We study the properties of slow magneto-acoustic waves that are naturally excited due to turbulent convection and investigate their role in the energy balance of a plage region using three dimensional (3D) radiation-MHD simulations. We calculate the horizontally averaged (over the whole domain) frequency power spectra for both longitudinal and vertical (i.e. the component perpendicular to the surface) components of velocity. To compare our results with the observations we degrade the simulation data with Gaussian kernels having FWHM of 100 km and 200 km, and calculate horizontally averaged power spectra for the vertical component of velocity. The power spectra of the longitudinal component of velocity, averaged over field lines in the core of a kG magnetic flux concentration, reveal that the dominant period of oscillations shifts from around 6.5 minutes in the photosphere to around 4 minutes in the chromosphere. At the same time, the velocity power spectra, averaged horizontally over the whole domain, show that low frequency waves (approximately 6.5 minute period) may reach well into the chromosphere. Importantly, waves with frequencies above 5 mHz propagating along different field lines are found to be out of phase with each other even within a single magnetic concentration. The horizontally averaged power spectra of the vertical component of velocity at various effective resolutions show that the observed acoustic wave energy fluxes are underestimated, by a factor of three even if determined from observations carried out at a high spatial resolution of 200 km. Our results show that longitudinal waves carry (just) sufficient energy to heat the chromosphere in solar plage. We conjecture that current observations (with spatial resolution around 200 km) underestimate the energy flux by roughly a factor of three, or more if the observations have lower spatial resolution. △ Less

Submitted 6 May, 2021; originally announced May 2021.

Comments: Accepted for publication in A & A

Journal ref: A&A 652, A43 (2021)

Non-equilibrium Equation of State in stellar atmospheres

Authors: L. S. Anusha, M. van Noort, R. H. Cameron

Abstract: In the stellar chromospheres, radiative energy transport is dominated by only the strongest spectral lines. For these lines, the approximation of local thermodynamic equilibrium (LTE) is known to be very inaccurate, and a state of equilibrium cannot be assumed in general. To calculate the radiative energy transport under these conditions, the population evolution equation must be evaluated explici… ▽ More In the stellar chromospheres, radiative energy transport is dominated by only the strongest spectral lines. For these lines, the approximation of local thermodynamic equilibrium (LTE) is known to be very inaccurate, and a state of equilibrium cannot be assumed in general. To calculate the radiative energy transport under these conditions, the population evolution equation must be evaluated explicitly, including all time-dependent terms. We develop a numerical method to solve the evolution equation for the atomic-level populations in a time-implicit way, keeping all time-dependent terms to first order. We show that the linear approximation of the time dependence of the populations can handle very large time steps without losing the accuracy. We reproduce the benchmark solutions from earlier, well-established works in terms of non-LTE kinetic equilibrium solution and typical ionization/recombination time-scales in the solar chromosphere. △ Less

Submitted 28 April, 2021; originally announced April 2021.

Comments: 42 pages, 8 figures

Journal ref: ApJ, 911, 71 (2021)

A journey of exploration to the polar regions of a star: probing the solar poles and the heliosphere from high helio-latitude

Authors: Louise Harra, Vincenzo Andretta, Thierry Appourchaux, Frédéric Baudin, Luis Bellot-Rubio, Aaron C. Birch, Patrick Boumier, Robert H. Cameron, Matts Carlsson, Thierry Corbard, Jackie Davies, Andrew Fazakerley, Silvano Fineschi, Wolfgang Finsterle, Laurent Gizon, Richard Harrison, Donald M. Hassler, John Leibacher, Paulett Liewer, Malcolm MacDonald, Milan Maksimovic, Neil Murphy, Giampiero Naletto, Giuseppina Nigro, Christopher Owen , et al. (7 additional authors not shown)

Abstract: A mission to view the solar poles from high helio-latitudes (above 60$^\circ$) will build on the experience of Solar Orbiter as well as a long heritage of successful solar missions and instrumentation (e.g. SOHO \cite{SOHO}, STEREO \cite{stereo}, Hinode \cite{Hinode}, SDO \cite{SDO}), but will focus for the first time on the solar poles, enabling scientific investigations that cannot be done by an… ▽ More A mission to view the solar poles from high helio-latitudes (above 60$^\circ$) will build on the experience of Solar Orbiter as well as a long heritage of successful solar missions and instrumentation (e.g. SOHO \cite{SOHO}, STEREO \cite{stereo}, Hinode \cite{Hinode}, SDO \cite{SDO}), but will focus for the first time on the solar poles, enabling scientific investigations that cannot be done by any other mission. One of the major mysteries of the Sun is the solar cycle. The activity cycle of the Sun drives the structure and behaviour of the heliosphere and is, of course, the driver of space weather. In addition, solar activity and variability provides fluctuating input into the Earth climate models, and these same physical processes are applicable to stellar systems hosting exoplanets. One of the main obstructions to understanding the solar cycle, and hence all solar activity, is our current lack of understanding of the polar regions. In this White Paper, submitted to the European Space Agency in response to the Voyage 2050 call, we describe a mission concept that aims to address this fundamental issue. In parallel, we recognise that viewing the Sun from above the polar regions enables further scientific advantages, beyond those related to the solar cycle, such as unique and powerful studies of coronal mass ejection processes, from a global perspective, and studies of coronal structure and activity in polar regions. Not only will these provide important scientific advances for fundamental stellar physics research, they will feed into our understanding of impacts on the Earth and other planets' space environment. △ Less

Submitted 22 April, 2021; originally announced April 2021.

Comments: 29 pages, 9 figures, White paper for ESA Voyage 2050

Vortex Flow Properties in Simulations of Solar Plage Region: Evidence for their role in chromospheric heating

Authors: Nitin Yadav, Robert H. Cameron, Sami K. Solanki

Abstract: Vortex-flows exist across a broad range of spatial and temporal scales in the solar atmosphere. Small-scale vortices have been proposed to play an important role in energy transport in the solar atmosphere. However, their physical properties remain poorly understood due to the limited spatial resolution of the observations. We aim to explore and analyze the physical properties of small-scale vorti… ▽ More Vortex-flows exist across a broad range of spatial and temporal scales in the solar atmosphere. Small-scale vortices have been proposed to play an important role in energy transport in the solar atmosphere. However, their physical properties remain poorly understood due to the limited spatial resolution of the observations. We aim to explore and analyze the physical properties of small-scale vortices inside magnetic flux tubes using numerical simulations, and to investigate whether they contribute to heating the chromosphere in a plage region. Using the three-dimensional (3D) radiative magnetohydrodynamic (MHD) simulation code 'MURaM', we perform numerical simulations of a unipolar solar plage region. To detect and isolate vortices, we use the Swirling Strength criterion and select the locations where the fluid is rotating with an angular velocity greater than a certain threshold. We concentrate on small-scale as they are the strongest and carry most of the energy. We explore the spatial profiles of physical quantities viz. density, horizontal velocity, etc. inside these vortices. Moreover, to apprehend their general characteristics, a statistical investigation is performed. Magnetic flux tubes have a complex filamentary substructure harbouring an abundance of small-scale vortices. At the interfaces between vortices strong current sheets are formed that may dissipate and heat the solar chromosphere. Statistically, vortices have higher densities and higher temperatures than the average values at the same geometrical height in the chromosphere. We conclude that small-scale vortices are ubiquitous in solar plage regions, and they are denser and hotter structures that contribute to chromospheric heating, possibly by dissipation of the current sheets formed at their interfaces. △ Less

Submitted 28 October, 2020; originally announced October 2020.

Comments: accepted for publication in A&A

Journal ref: A&A 645, A3 (2021)

Power spectrum of turbulent convection in the solar photosphere

Authors: L. Yelles Chaouche, R. H. Cameron, S. K. Solanki, T. L. Riethmüller, L. S. Anusha, V. Witzke, A. I. Shapiro, P. Barthol, A. Gandorfer, L. Gizon, J. Hirzberger, M. van Noort, J. Blanco Rodríguez, J. C. Del Toro Iniesta, D. Orozco Suárez, W. Schmidt, V. Martínez Pillet, M. Knölker

Abstract: 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… ▽ More 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. △ Less

Submitted 18 October, 2020; originally announced October 2020.

Comments: 10 pages, 7 figures, Accepted in A and A

Journal ref: A&A 644, A44 (2020)

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"

Authors: T. Reinhold, A. I. Shapiro, S. K. Solanki, B. T. Montet, N. A. Krivova, R. H. Cameron, E. M. Amazo-Gómez

Abstract: 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. 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. △ Less

Submitted 9 July, 2020; originally announced July 2020.

Average motion of emerging solar active region polarities II: Joy's law

Authors: Hannah Schunker, Christian Baumgartner, Aaron C. Birch, Robert H. Cameron, Douglas C. Braun, Laurent Gizon

Abstract: The tilt of solar active regions described by Joy's law is essential for converting a toroidal field to a poloidal field in Babcock-Leighton dynamo models. In thin flux tube models the Coriolis force causes Joy's law, acting on east-west flows as they rise towards the surface. Our goal is to measure the evolution of the average tilt angle of hundreds of active regions as they emerge, so that we ca… ▽ More The tilt of solar active regions described by Joy's law is essential for converting a toroidal field to a poloidal field in Babcock-Leighton dynamo models. In thin flux tube models the Coriolis force causes Joy's law, acting on east-west flows as they rise towards the surface. Our goal is to measure the evolution of the average tilt angle of hundreds of active regions as they emerge, so that we can constrain the origins of Joy's law. We measured the tilt angle of the primary bipoles in 153 emerging active regions in the Solar Dynamics Observatory Helioseismic Emerging Active Region survey. We used line-of-sight magnetic field measurements averaged over 6 hours to define the polarities and measure the tilt angle up to four days after emergence. We find that at the time of emergence the polarities are on average aligned east-west, and that neither the separation nor the tilt depends on latitude. We do find, however, that ARs at higher latitudes have a faster north-south separation speed than those closer to the equator at the emergence time. After emergence, the tilt angle increases and Joy's law is evident about two days later. The scatter in the tilt angle is independent of flux until about one day after emergence, when higher-flux regions have a smaller scatter in tilt angle than lower-flux regions. Our finding that active regions emerge with an east-west alignment is surprising since thin flux tube models predict that tilt angles of rising flux tubes are generated below the surface. Previously reported tilt angle relaxation of deeply anchored flux tubes can be largely explained by the change in east-west separation. We conclude that Joy's law is caused by an inherent north-south separation speed present when the flux first reaches the surface, and that the scatter in the tilt angle is consistent with buffeting of the polarities by supergranulation. △ Less

Submitted 9 June, 2020; originally announced June 2020.

Comments: accepted in Astronomy and Astrophysics

Journal ref: A&A 640, A116 (2020)

The Sun is less active than other solar-like stars

Authors: Timo Reinhold, Alexander I. Shapiro, Sami K. Solanki, Benjamin T. Montet, Nathalie A. Krivova, Robert H. Cameron, Eliana M. Amazo-Gomez

Abstract: Magnetic activity of the Sun and other stars causes their brightness to vary. We investigate how typical the Sun's variability is compared to other solar-like stars, i.e. those with near-solar effective temperatures and rotation periods. By combining four years of photometric observations from the Kepler space telescope with astrometric data from the Gaia spacecraft, we measure photometric variabi… ▽ More Magnetic activity of the Sun and other stars causes their brightness to vary. We investigate how typical the Sun's variability is compared to other solar-like stars, i.e. those with near-solar effective temperatures and rotation periods. By combining four years of photometric observations from the Kepler space telescope with astrometric data from the Gaia spacecraft, we measure photometric variabilities of 369 solar-like stars. Most of the solar-like stars with well-determined rotation periods show higher variability than the Sun and are therefore considerably more active. These stars appear nearly identical to the Sun, except for their higher variability. Their existence raises the question of whether the Sun can also experience epochs of such high variability. △ Less

Submitted 4 May, 2020; originally announced May 2020.

Comments: Accepted for publication in Science. 3 (main) + 10 (supplementary) figures

Journal ref: Science 01 May 2020: Vol. 368, Issue 6490, pp. 518-521

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

Authors: Nitin Yadav, R. H. Cameron, S. K. Solanki

Abstract: 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 smalles… ▽ More 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. △ Less

Submitted 29 April, 2020; originally announced April 2020.

Comments: 8 pages, 5 figures, accepted in ApJL

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

Authors: C. Damiani, R. H. Cameron, A. C. Birch, L. Gizon

Abstract: 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 polytro… ▽ More 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$. △ Less

Submitted 11 March, 2020; originally announced March 2020.

Comments: 10 pages, 8 figures, accepted in A&A

Journal ref: A&A 637, A65 (2020)

Power spectra of solar brightness variations at different inclinations

Authors: N. -E. Nèmec, A. I. Shapiro, N. A. Krivova, S. K. Solanki, R. V. Tagirov, R. H. Cameron, S. Dreizler

Abstract: 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, h… ▽ More 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. △ Less

Submitted 25 February, 2020; originally announced February 2020.

Comments: 11 pages, 9 figures plus 3 additional figures in the appendix, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 636, A43 (2020)

Loss of toroidal magnetic flux by emergence of bipolar magnetic regions

Authors: R. H. Cameron, M. Schuessler

Abstract: 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,… ▽ More 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. △ Less

Submitted 13 February, 2020; originally announced February 2020.

Comments: A&A accepted

Journal ref: A&A 636, A7 (2020)

Average motion of emerging solar active region polarities I: Two phases of emergence

Authors: Hannah Schunker, Aaron C. Birch, Robert H. Cameron, Douglas C. Braun, Laurent Gizon, Raymond B. Burston

Abstract: Our goal is to constrain models of active region formation by tracking the average motion of active region polarity pairs as they emerge onto the surface. We measured the motion of the two main opposite polarities in 153 emerging active regions (EARs) using line-of-sight magnetic field observations from the Solar Dynamics Observatory Helioseismic Emerging Active Region (SDO/HEAR) survey (Schunker… ▽ More Our goal is to constrain models of active region formation by tracking the average motion of active region polarity pairs as they emerge onto the surface. We measured the motion of the two main opposite polarities in 153 emerging active regions (EARs) using line-of-sight magnetic field observations from the Solar Dynamics Observatory Helioseismic Emerging Active Region (SDO/HEAR) survey (Schunker et al. 2016). We first measured the position of each of the polarities eight hours after emergence and tracked their location forwards and backwards in time. We find that, on average, the polarities emerge with an east-west orientation and the separation speed between the polarities increases. At about 0.1 days after emergence, the average separation speed reaches a peak value of 229 +/- 11 m/s, and then starts to decrease, and about 2.5 days after emergence the polarities stop separating. We also find that the separation and the separation speed in the east-west direction are systematically larger for active regions with higher flux. Our results reveal two phases of the emergence process defined by the rate of change of the separation speed as the polarities move apart. Phase 1 begins when the opposite polarity pairs first appear at the surface, with an east-west alignment and an increasing separation speed. We define Phase 2 to begin when the separation speed starts to decrease, and ends when the polarities have stopped separating. This is consistent with the picture of Chen, Rempel, & Fan (2017): the peak of a flux tube breaks through the surface during Phase 1. During Phase 2 the magnetic field lines are straightened by magnetic tension, so that the polarities continue to move apart, until they eventually lie directly above their anchored subsurface footpoints. △ Less

Submitted 30 April, 2019; v1 submitted 28 March, 2019; originally announced March 2019.

Comments: accepted A&A

Journal ref: A&A 625, A53 (2019)

Starspot rotation rates vs. activity cycle phase: Butterfly diagrams of Kepler stars are unlike the Sun's

Authors: M. B. Nielsen, L. Gizon, R. H. Cameron, M. Miesch

Abstract: During the solar magnetic activity cycle the emergence latitudes of sunspots change, leading to the well-known butterfly diagram. This phenomenon is poorly understood for other stars as starspot latitudes are generally unknown. The related changes in starspot rotation rates caused by latitudinal differential rotation can however be measured. Using the set of 3093 Kepler stars with activity cycles… ▽ More During the solar magnetic activity cycle the emergence latitudes of sunspots change, leading to the well-known butterfly diagram. This phenomenon is poorly understood for other stars as starspot latitudes are generally unknown. The related changes in starspot rotation rates caused by latitudinal differential rotation can however be measured. Using the set of 3093 Kepler stars with activity cycles identified by Reinhold et al. (2017), we aim to study the temporal change in starspot rotation rates over magnetic activity cycles, and how this relates to the activity level, mean rotation rate, and effective temperature of the star. We measure the photometric variability as a proxy for the magnetic activity and the spot rotation rate in each quarter over the duration of the Kepler mission. We phase-fold these measurements with the cycle period. We perform averages over stars with comparable mean rotation rates and effective temperature at fixed activity-cycle phases. We detect a clear correlation between the variation of activity level and the variation of the starspot rotation rate. The sign and amplitude of this correlation depends on the mean stellar rotation and, to a lesser extent, on the effective temperature. For slowly rotating stars (with periods between 15-28 days) the starspot rotation rates are clearly anti-correlated with the level of activity during the activity cycles. A transition is seen at periods of 10-15 days, where stars with effective temperature above 4200K instead show positive correlation. Our results can be interpreted in terms of a stellar "butterfly diagram", but these appear different from the Sun's as the starspot rotation rates are either in phase or anti-phase with the activity level. Alternatively, the activity cycles seen by Kepler are short (around 2.5 years) and may therefore be secondary cycles, perhaps analogous to the solar quasi-biennial oscillations. △ Less

Submitted 16 December, 2018; originally announced December 2018.

Comments: Accepted for publication in Astronomy & Astrophysics. 11 pages, 8 figures

Journal ref: A&A 622, A85 (2019)

Origin of the hemispheric asymmetry of solar activity

Authors: M. Schüssler, R. H. Cameron

Abstract: 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-t… ▽ More 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. △ Less

Submitted 26 July, 2018; originally announced July 2018.

Comments: Astronomy & Astrophysics, accepted

Journal ref: A&A 618, A89 (2018)

The nature of solar brightness variations

Authors: A. I. Shapiro, S. K. Solanki, N. A. Krivova, R. H. Cameron, K. L. Yeo, W. K. Schmutz

Abstract: The solar brightness varies on timescales from minutes to decades. Determining the sources of such variations, often referred to as solar noise, is of importance for multiple reasons: a) it is the background that limits the detection of solar oscillations, b) variability in solar brightness is one of the drivers of the Earth's climate system, c) it is a prototype of stellar variability which is an… ▽ More The solar brightness varies on timescales from minutes to decades. Determining the sources of such variations, often referred to as solar noise, is of importance for multiple reasons: a) it is the background that limits the detection of solar oscillations, b) variability in solar brightness is one of the drivers of the Earth's climate system, c) it is a prototype of stellar variability which is an important limiting factor for the detection of extra-solar planets. Here we show that recent progress in simulations and observations of the Sun makes it finally possible to pinpoint the source of the solar noise. We utilise high-cadence observations from the Solar Dynamic Observatory and the SATIRE model to calculate the magnetically-driven variations of solar brightness. The brightness variations caused by the constantly evolving cellular granulation pattern on the solar surface are computed with the MURAM code. We find that surface magnetic field and granulation can together precisely explain solar noise on timescales from minutes to decades, i.e. ranging over more than six orders of magnitude in the period. This accounts for all timescales that have so far been resolved or covered by irradiance measurements. We demonstrate that no other sources of variability are required to explain the data. Recent measurements of Sun-like stars by CoRoT and Kepler uncovered brightness variations similar to that of the Sun but with much wider variety of patterns. Our finding that solar brightness variations can be replicated in detail with just two well-known sources will greatly simplify future modelling of existing CoRoT and Kepler as well as anticipated TESS and PLATO data. △ Less

Submitted 26 November, 2017; v1 submitted 11 November, 2017; originally announced November 2017.

Comments: This is the submitted version of the paper published in Nature Astronomy

Observing and modeling the poloidal and toroidal fields of the solar dynamo

Authors: R. H. Cameron, T. L Duvall Jr., M. Schüssler, H. Schunker

Abstract: 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 fie… ▽ More 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. △ Less

Submitted 19 October, 2017; originally announced October 2017.

Comments: A&A

Journal ref: A&A 609, A56 (2018)

Evidence for photometric activity cycles in 3203 Kepler stars

Authors: Timo Reinhold, Robert H. Cameron, Laurent Gizon

Abstract: In recent years it has been claimed that the length of stellar activity cycles is determined by the stellar rotation rate. It is observed that the cycle period increases with rotation period along the so-called active and inactive sequences. In this picture the Sun occupies a solitary position in between the two sequences. Our goal is to measure cyclic variations of the stellar light curve amplitu… ▽ More In recent years it has been claimed that the length of stellar activity cycles is determined by the stellar rotation rate. It is observed that the cycle period increases with rotation period along the so-called active and inactive sequences. In this picture the Sun occupies a solitary position in between the two sequences. Our goal is to measure cyclic variations of the stellar light curve amplitude and the rotation period using four years of Kepler data. Periodic changes of the light curve amplitude or the stellar rotation period are associated with an underlying activity cycle. Using the McQuillan et al. 2014 sample we compute the rotation period and the variability amplitude for each Kepler quarter and search for periodic variations of both time series. To test for periodicity in each stellar time series we consider Lomb-Scargle periodograms and use a selection based on a False Alarm Probability (FAP). We detect amplitude periodicities in 3203 stars between 0.5-6 years covering rotation periods between 1-40 days. Given our sample size of 23,601 stars and our selection criteria that the FAP is less than 5%, this number is almost three times higher than that expected from pure noise. We do not detect periodicities in the rotation period beyond those expected from noise. Our measurements reveal that the cycle period shows a weak dependence on rotation rate, slightly increasing for longer rotation period. We further show that the shape of the variability deviates from a pure sine curve, consistent with observations of the solar cycle. Our measurements do not support the existence of distinct sequences in the P_rot-P_cyc plane, although there is some evidence for the inactive sequence for rotation periods between 5-25 days. Unfortunately, the total observing time is too short to draw sound conclusions on activity cycles with similar length as the solar cycle. △ Less

Submitted 9 May, 2017; originally announced May 2017.

Comments: 8 pages, 9 figures. We will publish a table containing the measured cycle periods and the corresponding KIC numbers

Journal ref: A&A 603, A52 (2017)

Evolution of the Sun's non-axisymmetric toroidal field

Authors: David Martin-Belda, Robert H. Cameron

Abstract: Aims. We aim to infer the sub-surface distribution of the Sun's non-axisymmetric azimuthal magnetic flux from observable quantities, such as the surface magnetic field and the large scale plasma flows. Methods. We have built a kinematic flux transport model of the solar dynamo based on the Babcock-Leighton framework. We constructed the source term for the poloidal field using SOLIS magnetograms… ▽ More Aims. We aim to infer the sub-surface distribution of the Sun's non-axisymmetric azimuthal magnetic flux from observable quantities, such as the surface magnetic field and the large scale plasma flows. Methods. We have built a kinematic flux transport model of the solar dynamo based on the Babcock-Leighton framework. We constructed the source term for the poloidal field using SOLIS magnetograms spanning three solar cycles. Based on this source we calculated the azimuthal flux below the surface. The flux transport model has two free parameters which we constrain using sunspot observations from cycle 22. We compared the model results with observations from cycle 23. Results. The structure of the azimuthal field is mainly axisymmetric. The departures from axisymmetry represent, on average, ~3% of the total azimuthal flux. Owing to its relative weakness, the non-axisymmetric structure of the azimuthal field does not have a significant impact on the location in which the emergences appear or on the amount of flux contained in them. We find that the probability of emergence is a function of the ratio between the flux content of an active region and the underlying azimuthal flux. △ Less

Submitted 29 March, 2017; originally announced March 2017.

Journal ref: A&A 603, A53 (2017)

High-frequency Oscillations in Small Magnetic Elements Observed with Sunrise/SuFI

Authors: Shahin Jafarzadeh, S. K. Solanki, M. Stangalini, O. Steiner, R. H. Cameron, S. Danilovic

Abstract: We characterize waves in small magnetic elements and investigate their propagation in the lower solar atmosphere from observations at high spatial and temporal resolution. We use the wavelet transform to analyze oscillations of both horizontal displacement and intensity in magnetic bright points found in the 300 nm and the Ca II H 396.8 nm passbands of the filter imager on board the Sunrise balloo… ▽ More We characterize waves in small magnetic elements and investigate their propagation in the lower solar atmosphere from observations at high spatial and temporal resolution. We use the wavelet transform to analyze oscillations of both horizontal displacement and intensity in magnetic bright points found in the 300 nm and the Ca II H 396.8 nm passbands of the filter imager on board the Sunrise balloon-borne solar observatory. Phase differences between the oscillations at the two atmospheric layers corresponding to the two passbands reveal upward propagating waves at high frequencies (up to 30 mHz). Weak signatures of standing as well as downward propagating waves are also obtained. Both compressible and incompressible (kink) waves are found in the small-scale magnetic features. The two types of waves have different, though overlapping, period distributions. Two independent estimates give a height difference of approximately 450+-100 km between the two atmospheric layers sampled by the employed spectral bands. This value, together with the determined short travel times of the transverse and longitudinal waves provide us with phase speeds of 29+-2 km/s and 31+-2 km/s, respectively. We speculate that these phase speeds may not reflect the true propagation speeds of the waves. Thus, effects such as the refraction of fast longitudinal waves may contribute to an overestimate of the phase speed. △ Less

Submitted 30 March, 2017; v1 submitted 28 November, 2016; originally announced November 2016.

Comments: 14 pages, 7 figures

Journal ref: The Astrophysical Journal Supplement Series, 229, 10 (2017)

An update of Leighton's solar dynamo model

Authors: R. H. Cameron, M. Schuessler

Abstract: 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, bot… ▽ More 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). △ Less

Submitted 28 November, 2016; originally announced November 2016.

Comments: A&A accepted

Journal ref: A&A 599, A52 (2017)

Kinematics of Magnetic Bright Features in the Solar Photosphere

Authors: Shahin Jafarzadeh, S. K. Solanki, R. H. Cameron, P. Barthol, J. Blanco Rodriguez, J. C. del Toro Iniesta, A. Gandorfer, L. Gizon, J. Hirzberger, M. Knoelker, V. Martinez Pillet, D. Orozco Suarez, T. L. Riethmueller, W. Schmidt, M. van Noort

Abstract: Convective flows are known as the prime means of transporting magnetic fields on the solar surface. Thus, small magnetic structures are good tracers of the turbulent flows. We study the migration and dispersal of magnetic bright features (MBFs) in intergranular areas observed at high spatial resolution with Sunrise/IMaX. We describe the flux dispersal of individual MBFs as a diffusion process whos… ▽ More Convective flows are known as the prime means of transporting magnetic fields on the solar surface. Thus, small magnetic structures are good tracers of the turbulent flows. We study the migration and dispersal of magnetic bright features (MBFs) in intergranular areas observed at high spatial resolution with Sunrise/IMaX. We describe the flux dispersal of individual MBFs as a diffusion process whose parameters are computed for various areas in the quiet Sun and the vicinity of active regions from seeing-free data. We find that magnetic concentrations are best described as random walkers close to network areas (diffusion index, gamma=1.0), travelers with constant speeds over a supergranule (gamma=1.9-2.0), and decelerating movers in the vicinity of flux emergence and/or within active regions (gamma=1.4-1.5). The three types of regions host MBFs with mean diffusion coefficients of 130 km^2/s, 80-90 km^2/s, and 25-70 km^2/s, respectively. The MBFs in these three types of regions are found to display a distinct kinematic behavior at a confidence level in excess of 95%. △ Less

Submitted 29 March, 2017; v1 submitted 24 October, 2016; originally announced October 2016.

Comments: 8 pages, 4 figures

Journal ref: The Astrophysical Journal Supplement Series, 229, 8 (2017)

Inflows towards active regions and the modulation of the solar cycle: a parameter study

Authors: David Martin-Belda, Robert H. Cameron

Abstract: Aims: We aim to investigate how converging flows towards active regions affect the surface transport of magnetic flux, as well as their impact on the generation of the Sun's poloidal field. The inflows constitute a potential non-linear mechanism for the saturation of the global dynamo and may contribute to the modulation of the solar cycle in the Babcock-Leighton framework. Methods: We build a s… ▽ More Aims: We aim to investigate how converging flows towards active regions affect the surface transport of magnetic flux, as well as their impact on the generation of the Sun's poloidal field. The inflows constitute a potential non-linear mechanism for the saturation of the global dynamo and may contribute to the modulation of the solar cycle in the Babcock-Leighton framework. Methods: We build a surface flux transport code incorporating a parametrized model of the inflows and run simulations spanning several cycles. We carry out a parameter study to assess how the strength and extension of the inflows affect the build-up of the global dipole field. We also perform simulations with different levels of activity to investigate the potential role of the inflows in the saturation of the global dynamo. Results: We find that the interaction of neighbouring active regions can lead to the occasional formation of single-polarity magnetic flux clumps inconsistent with observations. We propose the darkening caused by pores in areas of high magnetic field strength as a plausible mechanism preventing this flux-clumping. We find that inflows decrease the amplitude of the axial dipole moment by a $\sim30\,\%$, relative to a no-inflows scenario. Stronger (weaker) inflows lead to larger (smaller) reductions of the axial dipole moment. The relative amplitude of the generated axial dipole is about $9\%$ larger after very weak cycles than after very strong cycles. This supports the inflows as a non-linear mechanism capable of saturating the global dynamo and contributing to the modulation of the solar cycle within the Babcock-Leighton framework. △ Less

Submitted 5 September, 2016; originally announced September 2016.

Journal ref: A&A 597, A21 (2017)

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

Authors: W. H. Ball, B. Beeck, R. H. Cameron, L. Gizon

Abstract: ... [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… ▽ More ... [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. △ Less

Submitted 8 June, 2016; originally announced June 2016.

Comments: 8 pages, 7 figures. Accepted for publication in Astronomy & Astrophysics

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

Authors: R. H. Cameron, M. Schuessler

Abstract: 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 an… ▽ More 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. △ Less

Submitted 25 April, 2016; originally announced April 2016.

Comments: A&A, accepted

Solar cycle 25: another moderate cycle?

Authors: Robert H. Cameron, Jie Jiang, Manfred Schuessler

Abstract: Surface flux transport simulations for the descending phase of cycle 24 using random sources (emerging bipolar magnetic regions) with empirically determined scatter of their properties provide a prediction of the axial dipole moment during the upcoming activity minimum together with a realistic uncertainty range. The expectation value for the dipole moment around 2020 $(2.5\pm1.1\,$G) is comparabl… ▽ More Surface flux transport simulations for the descending phase of cycle 24 using random sources (emerging bipolar magnetic regions) with empirically determined scatter of their properties provide a prediction of the axial dipole moment during the upcoming activity minimum together with a realistic uncertainty range. The expectation value for the dipole moment around 2020 $(2.5\pm1.1\,$G) is comparable to that observed at the end of cycle 23 (about $2\,$G). The empirical correlation between the dipole moment during solar minimum and the strength of the subsequent cycle thus suggests that cycle 25 will be of moderate amplitude, not much higher than that of the current cycle. However, the intrinsic uncertainty of such predictions resulting from the random scatter of the source properties is considerable and fundamentally limits the reliability with which such predictions can be made before activity minimum is reached. △ Less

Submitted 6 May, 2016; v1 submitted 18 April, 2016; originally announced April 2016.

Comments: 13 papges, 4 figures,Accepted for publication in ApJL

The Global Solar Dynamo

Authors: R. H. Cameron, M. Dikpati, A. Brandenburg

Abstract: A brief summary of the various observations and constraints that underlie solar dynamo research are presented. The arguments that indicate that the solar dynamo is an alpha-omega dynamo of the Babcock-Leighton type are then shortly reviewed. The main open questions that remain are concerned with the subsurface dynamics, including why sunspots emerge at preferred latitudes as seen in the familiar b… ▽ More A brief summary of the various observations and constraints that underlie solar dynamo research are presented. The arguments that indicate that the solar dynamo is an alpha-omega dynamo of the Babcock-Leighton type are then shortly reviewed. The main open questions that remain are concerned with the subsurface dynamics, including why sunspots emerge at preferred latitudes as seen in the familiar butterfly wings, why the cycle is about 11 years long, and why the sunspot groups emerge tilted with respect to the equator (Joy's law). Next, we turn to magnetic helicity, whose conservation property has been identified with the decline of large-scale magnetic fields found in direct numerical simulations at large magnetic Reynolds numbers. However, magnetic helicity fluxes through the solar surface can alleviate this problem and connect theory with observations, as will be discussed. △ Less

Submitted 4 February, 2016; originally announced February 2016.

Comments: Accepted in Space Science Reviews NORDITA-2015-97

Report number: NORDITA-2015-97

Journal ref: Spa. Sci. Rev. 210, 367-395 (2017)

Surface flux transport simulations: Effect of inflows toward active regions and random velocities on the evolution of the Sun's large-scale magnetic field

Authors: David Martin-Belda, Robert H. Cameron

Abstract: Aims: We aim to determine the effect of converging flows on the evolution of a bipolar magnetic region (BMR), and to investigate the role of these inflows in the generation of poloidal flux. We also discuss whether the flux dispersal due to turbulent flows can be described as a diffusion process. Methods: We developed a simple surface flux transport model based on point-like magnetic concentrati… ▽ More Aims: We aim to determine the effect of converging flows on the evolution of a bipolar magnetic region (BMR), and to investigate the role of these inflows in the generation of poloidal flux. We also discuss whether the flux dispersal due to turbulent flows can be described as a diffusion process. Methods: We developed a simple surface flux transport model based on point-like magnetic concentrations. We tracked the tilt angle, the magnetic flux and the axial dipole moment of a BMR in simulations with and without inflows and compared the results. To test the diffusion approximation, simulations of random walk dispersal of magnetic features were compared against the predictions of the diffusion treatment. Results: We confirm the validity of the diffusion approximation to describe flux dispersal on large scales. We find that the inflows enhance flux cancellation, but at the same time affect the latitudinal separation of the polarities of the bipolar region. In most cases the latitudinal separation is limited by the inflows, resulting in a reduction of the axial dipole moment of the BMR. However, when the initial tilt angle of the BMR is small, the inflows produce an increase in latitudinal separation that leads to an increase in the axial dipole moment in spite of the enhanced flux destruction. This can give rise to a tilt of the BMR even when the BMR was originally aligned parallel to the equator. △ Less

Submitted 14 December, 2015; v1 submitted 8 December, 2015; originally announced December 2015.

Journal ref: A&A 586, A73 (2016)

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

Authors: H. Peter, J. Warnecke, L. P. Chitta, R. H. Cameron

Abstract: 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… ▽ More 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. △ Less

Submitted 15 October, 2015; originally announced October 2015.

Comments: Accepted for publication in A&A (Oct.14,2015), 8 pages, 4 figures

Journal ref: A&A 584, A68 (2015)

The cause of the weak solar cycle 24

Authors: Jie Jiang, Robert H. Cameron, Manfred Schuessler

Abstract: The ongoing 11-year cycle of solar activity is considerably less vigorous than the three cycles before. It was preceded by a very deep activity minimum with a low polar magnetic flux, the source of the toroidal field responsible for solar magnetic activity in the subsequent cycle. Simulation of the evolution of the solar surface field shows that the weak polar fields and thus the weakness of the p… ▽ More The ongoing 11-year cycle of solar activity is considerably less vigorous than the three cycles before. It was preceded by a very deep activity minimum with a low polar magnetic flux, the source of the toroidal field responsible for solar magnetic activity in the subsequent cycle. Simulation of the evolution of the solar surface field shows that the weak polar fields and thus the weakness of the present cycle 24 are mainly caused by a number of bigger bipolar regions emerging at low latitudes with a `wrong' (i.e., opposite to the majority for this cycle) orientation of their magnetic polarities in the North-South direction, which impaired the growth of the polar field. These regions had a particularly strong effect since they emerged within $\pm10^\circ$ latitude from the solar equator. △ Less

Submitted 7 July, 2015; originally announced July 2015.

Comments: 15 pages, 5 figures, accepted for publication in ApJL

Three-dimensional simulations of near-surface convection in main-sequence stars. IV. Effect of small-scale magnetic flux concentrations on centre-to-limb variation and spectral lines

Authors: Benjamin Beeck, Manfred Schüssler, Robert H. Cameron, Ansgar Reiners

Abstract: Magnetic fields affect the local structure of the photosphere of stars. They can considerably influence the radiative properties near the optical surface, flow velocities, and the temperature and pressure profiles. We aim at understanding qualitatively the influence of small magnetic flux concentrations in unipolar plage regions on the centre-to-limb variation of the intensity and its contrast and… ▽ More Magnetic fields affect the local structure of the photosphere of stars. They can considerably influence the radiative properties near the optical surface, flow velocities, and the temperature and pressure profiles. We aim at understanding qualitatively the influence of small magnetic flux concentrations in unipolar plage regions on the centre-to-limb variation of the intensity and its contrast and on the shape of spectral line profiles in cool main-sequence stars. We analyse the bolometric and continuum intensity and its angular dependence of 24 radiative magnetohydrodynamic simulations of the near-surface layers of main-sequence stars with six different sets of stellar parameters (spectral types F to early M) and four different average magnetic field strengths (including the non-magnetic case). We also calculated disc-integrated profiles of three spectral lines. The small magnetic flux concentrations formed in the magnetic runs of simulations have a considerable impact on the intensity and its centre-to-limb variation. Spectral lines are not only broadened owing to the Zeeman effect, but are also strongly affected by the modified thermodynamical structure and flow patterns. This indirect magnetic impact on the line profiles is often bigger than that of the Zeeman effect. The effects of the magnetic field on the radiation leaving the star can be considerable and is not restricted to spectral line broadening and polarisation by the Zeeman effect. The inhomogeneous structure of the magnetic field on small length scales and its impact on (and spatial correlation with) the local thermodynamical structure and the flow field near the surface influence the measurement of the global field properties and stellar parameters. These effects need to be taken into account in the interpretation of observations. △ Less

Submitted 18 May, 2015; originally announced May 2015.

Comments: 16 pages, 13+3 figures, 1 appendix, accepted for publication in A&A

Journal ref: A&A 581, A43 (2015)

Three-dimensional simulations of near-surface convection in main-sequence stars. III. The structure of small-scale magnetic flux concentrations

Authors: Benjamin Beeck, Manfred Schüssler, Robert H. Cameron, Ansgar Reiners

Abstract: The convective envelopes of cool main-sequence stars harbour magnetic fields with a complex global and local structure. These fields affect the near-surface convection and the outer stellar atmospheres in many ways and are responsible for the observable magnetic activity of stars. Our aim is to understand the local structure in unipolar regions with moderate average magnetic flux density. These co… ▽ More The convective envelopes of cool main-sequence stars harbour magnetic fields with a complex global and local structure. These fields affect the near-surface convection and the outer stellar atmospheres in many ways and are responsible for the observable magnetic activity of stars. Our aim is to understand the local structure in unipolar regions with moderate average magnetic flux density. These correspond to plage regions covering a substantial fraction of the surface of the Sun (and likely also the surface of other Sun-like stars) during periods of high magnetic activity. We analyse the results of 18 local-box magnetohydrodynamics simulations covering the upper layers of the convection zones and the photospheres of cool main-sequence stars of spectral types F to early M. The average vertical field in these simulations ranges from 20 to 500G. We find a substantial variation of the properties of the surface magnetoconvection between main-sequence stars of different spectral types. As a consequence of a reduced efficiency of the convective collapse of flux tubes, M dwarfs lack bright magnetic structures in unipolar regions of moderate field strength. The spatial correlation between velocity and the magnetic field as well as the lifetime of magnetic structures and their sizes relative to the granules vary significantly along the model sequence of stellar types. △ Less

Submitted 18 May, 2015; originally announced May 2015.

Comments: 15 pages, 12 figures, accepted for publication in A&A

Journal ref: A&A 581, A42 (2015)

Simulated Magnetic Flows in the Solar Photosphere

Authors: S. Danilovic, R. H. Cameron, S. K. Solanki

Abstract: 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… ▽ More 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. △ Less

Submitted 26 August, 2014; originally announced August 2014.

Comments: accepted in A&A movies can be found at http://www2.mps.mpg.de/data/outgoing/danilovic/ss/

Magnetic Flux Transport at the Solar Surface

Authors: J. Jiang, D. H. Hathaway, R. H. Cameron, S. K. Solanki, L. Gizon, L. Upton

Abstract: 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… ▽ More 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. △ Less

Submitted 13 August, 2014; originally announced August 2014.

Comments: 39 pages, 15 figures, accepted for publication in Space Sci. Rev

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

Authors: J. W. Lord, R. H. Cameron, M. P. Rast, M. Rempel, T. Roudier

Abstract: 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 turb… ▽ More 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. △ Less

Submitted 11 July, 2014; v1 submitted 8 July, 2014; originally announced July 2014.

Comments: 36 pages, 6 figures